MathGroup Archive 2010

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: 2-D Fourier Transform?

  • To: mathgroup at
  • Subject: [mg113784] Re: 2-D Fourier Transform?
  • From: "Nasser M. Abbasi" <nma at>
  • Date: Fri, 12 Nov 2010 05:28:23 -0500 (EST)
  • References: <ibgivj$ea9$>
  • Reply-to: nma at

On 11/11/2010 3:10 AM, Matthias Odisio wrote:

> CenteredFourier[img_] := Module[{data = ImageData[img], dim},
>     dim = Dimensions[data];
>     Fourier[data*(-1)^Table[i + j, {i, First[dim]}, {j, Last[dim]}]]
>     ];

Hello Matthias;

I liked your way to center the image more than the way
I did it. Even though it is the same effect, your's is
little bit more clear, and clear code is important for

So, if you do not mind, I borrowed your line above to
center the image, and use it instead of mine.

I also now use Image to display the spectrum, not ArrayPlot.

Here is the updated complete script to the 2D spectrum of an image
(Mathematica 7)

----------Spectrum of 2D image ---------


data=ImageData[img];  (*get data*)

d=data[[All,All,2]]; (*get channel 2 to FFT but center it first*)

(*make FFT, and look at spectrum and phase*)

fw= Fourier[d,FourierParameters->{1,1}];
fudgeFactor=100; (*adjust for better viewing as needed*)
abs=fudgeFactor * Log[1+Abs@fw];
Labeled[Image[abs/Max[abs],ImageSize->500],Style["Magnitude spectrum", 18]]

Labeled[Image[arg/Max[arg],ImageSize->500],Style["Phase spectrum", 18]]



  • Prev by Date: Re: Balance point of a solid
  • Next by Date: Re: command to save as .m file
  • Previous by thread: Re: 2-D Fourier Transform?
  • Next by thread: Address to the Mathematica-users Wiki Homepage