Reply by divy...@gmail.com●January 22, 20092009-01-22
On Jan 21, 11:35�pm, "adam3914" <a...@lle.rochester.edu> wrote:
> I figured out that if I use the complex input fftw method it works just
> fine. �So now I pass in a complex image with the complex parts all zero,
> and everything works fine.
Cool. Glad to know you solved it.
Reply by adam3914●January 21, 20092009-01-21
>On Jan 20, 7:41=A0pm, "adam3914" <a...@lle.rochester.edu> wrote:
>> I have what I have to assume is a stupid mistake on my part. =A0I am
able=
> to
>> take the fft of an image and get back the half of the transformed image
a=
>s
>> I should. =A0When I try to do a inverse fft on the data from the fft I
ge=
>t
>> really big numbers like 10^71 size numbers. =A0When I would expect to
get=
> my
>> original image times 640*480, the size of my image. =A0My first
question =
>is
>> do I have to create the other half of the fft to pass to the inverse
fft?
>> Or do I need to scale the transformed data before doing an inverse
I figured out that if I use the complex input fftw method it works just
fine. So now I pass in a complex image with the complex parts all zero,
and everything works fine.
Reply by divy...@gmail.com●January 21, 20092009-01-21
On Jan 20, 7:41�pm, "adam3914" <a...@lle.rochester.edu> wrote:
> I have what I have to assume is a stupid mistake on my part. �I am able to
> take the fft of an image and get back the half of the transformed image as
> I should. �When I try to do a inverse fft on the data from the fft I get
> really big numbers like 10^71 size numbers. �When I would expect to get my
> original image times 640*480, the size of my image. �My first question is
> do I have to create the other half of the fft to pass to the inverse fft?
> Or do I need to scale the transformed data before doing an inverse fft?
> Thanks For any help.
> Adam
>
> My code
>
> double * image=(double *) malloc(640*480*sizeof(double));
> fftw_complex * out2d=(fftw_complex *)
> malloc(640*480*sizeof(fftw_complex));
> fftw_complex * out2c=(fftw_complex *)
> malloc(640*480*sizeof(fftw_complex));
>
> p2d = fftw_plan_dft_r2c_2d(480, 640, image,out2d,FFTW_MEASURE );
> c2c = fftw_plan_dft_2d(480, 640, out2d,out2c, FFTW_BACKWARD,FFTW_MEASURE
> );
> fftw_execute(p2d);
> fftw_execute(c2c);
Well, I don't use FFTW, so I can't help you with your code. But what
caught my eye was:
>take the fft of an image and get back the half of the transformed image as
>I should.
??? A 2D FFT on a 640x480 real input array (with zero valued imaginary
parts) results in a 640x480 complex output array. And that output
array has certain symmetry properties to it that differ from the 1D
FFT case (i.e.: for 1D, you have upper/lower half conjugates - for 2D,
it's a little more complicated). Are you presuming that you only need
half of the outputs, as would be the case if you only want the
positive frequency outputs for a 1D FFT on real inputs?
As for the inverse, you asked if you have to "create the other half of
the FFT ...." So I'm kind of presuming that perhaps you're thinking
about this as if it were the same as the 1D case, and it isn't. For a
640x480 2D FFT, forward or inverse, you should have a 640x480 array
that can hold complex numbers. For a real image on the input side,
your imaginary inputs should be zero. Your transformed result will be
640x480 complex (with non-zero imaginary parts). To do the inverse,
you should be starting out with a 640x480 complex input.
And yes, if your FFT is unscaled, then you should definitely scale by
dividing each output by the NxM size of the array at some place in the
processing. I usually do it after the forward and before the
inverse. Overflow is just way to easy to get when doing 2 or higher
dimensional FFTs, especially for integer processing.
Kevin
Reply by adam3914●January 20, 20092009-01-20
I have what I have to assume is a stupid mistake on my part. I am able to
take the fft of an image and get back the half of the transformed image as
I should. When I try to do a inverse fft on the data from the fft I get
really big numbers like 10^71 size numbers. When I would expect to get my
original image times 640*480, the size of my image. My first question is
do I have to create the other half of the fft to pass to the inverse fft?
Or do I need to scale the transformed data before doing an inverse fft?
Thanks For any help.
Adam
My code
double * image=(double *) malloc(640*480*sizeof(double));
fftw_complex * out2d=(fftw_complex *)
malloc(640*480*sizeof(fftw_complex));
fftw_complex * out2c=(fftw_complex *)
malloc(640*480*sizeof(fftw_complex));
p2d = fftw_plan_dft_r2c_2d(480, 640, image,out2d,FFTW_MEASURE );
c2c = fftw_plan_dft_2d(480, 640, out2d,out2c, FFTW_BACKWARD,FFTW_MEASURE
);
fftw_execute(p2d);
fftw_execute(c2c);