Inverse FFT of Stack

Previous Topic Next Topic
 
classic Classic list List threaded Threaded
2 messages Options
Reply | Threaded
Open this post in threaded view
|

Inverse FFT of Stack

Roger_Almeida
Hello to everyone.

I'm just trying to automatize the Inverse FFT process. Anyone know some Macro/Plugin to apply the Inverse FFT process in a Stack of images?

I'm using this Macro to get the FFt result:

a=getTitle(); // your stack
for(i=1;i<=nSlices;i++){
 setSlice(i);
 run("FFT");
 selectImage(a);
}
run("Images to Stack", "name=Stack title=[FFT of] use");


So I made this Macro to get the Inverse result:

a=getTitle(); // your stack
for(i=1;i<=nSlices;i++){
 setSlice(i);
 run("Inverse FFT");
 selectImage(a);
}
run("Images to Stack", "name=Stack title=[FFT of] use");

So ImageJ give this: "Frequency domain image reguired", I got this box even using the basic command (Inverse FFT).
Thanks to all.


Reply | Threaded
Open this post in threaded view
|

Re: Inverse FFT of Stack

Michael Schmid
Hi Roger,

the FFT image that you see is just an 8-bit version of the power spectrum (with logarithmic grayscale), the actual Fourier transform (needed for back transform) is internally stored and invisible.
So you can't have a stack of 'true' FFTs, only a stack of the 8-bit power spectra.

Technical details: The FFT (actually an FHT) is an invisible image property of the ImagePlus.
An ImagePlus can have only one FHT, so you can't have a stack (which is contained in a single ImagePlus) of multiple FFTs that can be transformed back.

Michael
________________________________________________________________
On Aug 5, 2015, at 19:33, Roger_Almeida wrote:

> Hello to everyone.
>
> I'm just trying to automatize the Inverse FFT process. Anyone know some
> Macro/Plugin to apply the Inverse FFT process in a Stack of images?
>
> I'm using this Macro to get the FFt result:
>
> a=getTitle(); // your stack
> for(i=1;i<=nSlices;i++){
> setSlice(i);
> run("FFT");
> selectImage(a);
> }
> run("Images to Stack", "name=Stack title=[FFT of] use");
>
>
> So I made this Macro to get the Inverse result:
>
> a=getTitle(); // your stack
> for(i=1;i<=nSlices;i++){
> setSlice(i);
> run("Inverse FFT");
> selectImage(a);
> }
> run("Images to Stack", "name=Stack title=[FFT of] use");
>
> So ImageJ give this: "Frequency domain image reguired", I got this box even
> using the basic command (Inverse FFT).
> Thanks to all.
>
>
>
>
>
>
> --
> View this message in context: http://imagej.1557.x6.nabble.com/Inverse-FFT-of-Stack-tp5013910.html
> Sent from the ImageJ mailing list archive at Nabble.com.
>
> --
> ImageJ mailing list: http://imagej.nih.gov/ij/list.html

--
ImageJ mailing list: http://imagej.nih.gov/ij/list.html