Opened 11 years ago

Closed 9 years ago

#549 closed defect (wontfix)

r.mfilter should _respect_ i.fft outputs for use with i.ifft

Reported by: nikos Owned by: grass-dev@…
Priority: major Milestone: 6.4.0
Component: Raster Version: svn-develbranch6
Keywords: i.fft, r.mfilter, i.ifft Cc: nikos.alexandris@…
CPU: Unspecified Platform: Unspecified

Description

# fft of some raster image
i.fft in=a real=a_real imaginary=a_imaginary

# do something with the real_fft_image, e.g. r.mfilter
r.mfilter in=a_real filter=some_filter out=a_real_filtered

# perform i.ifft
i.ifft real=a_real_filtered imaginary=a_imaginary output=a_via_fft_filtered

ERROR: Unable to open real-image in the 'cell_misc' directory. Raster map probably wasn't created by i.fft

# questions

Is it fair enough to ask for this support?
Or is there a reason why _r.mfilter_ alters the i.fft-products in such a way that _i.ifft_ can't handle?

Nikos

Change History (4)

comment:1 Changed 11 years ago by nikos

Well, I bombed (again) in an issue that I had seen in the past but did not remember of it anymore.

Read Question about MODIS band 5 destriping?

comment:2 in reply to:  description Changed 11 years ago by glynn

Replying to nikos:

Is it fair enough to ask for this support?

It's already present in 7.0, but the interface had to be changed (separate maps for the real and imaginary components), so it's unlikely to be back-ported.

Or is there a reason why _r.mfilter_ alters the i.fft-products in such a way that _i.ifft_ can't handle?

i.fft stores the actual frequency-domain data in auxiliary data files. The output map only exists for viewing with d.rast. Apart from anything else, a raster map can only store real values, and the FFT results are complex values.

i.ifft operates upon the auxiliary data files, not the map itself, so it only works with maps which were generated by i.fft. The only processing which can be performed on the frequency-domain data is to apply a mask (i.ifft explicitly applies the MASK raster if it exists).

comment:3 Changed 9 years ago by nikos

Should this ticket be marked as will-not-be-fixed for grass64?

Nikos

comment:4 Changed 9 years ago by neteler

Resolution: wontfix
Status: newclosed
Note: See TracTickets for help on using tickets.