Eng-Tips is the largest engineering community on the Internet

Intelligent Work Forums for Engineering Professionals

Low-pass filtering via FFTW: wrong scaling 1

Status
Not open for further replies.

norpoth

Electrical
May 5, 2006
2
0
0
DE
Dear members,
in order to convolve a greyscale image with a gaussian I used FFTW code (release 3.1.1) to work in fourier domain. In detail my initial data is lying on a 1000x10000 point grid, correspondingly I put the gaussian bell on a 1000x1000 grid too. The gaussian distribution is normalized (i.e. integration over complete 2D equals 1) and centered on the grid around (500.5, 500.5). After transforming both via complex-to-complex 2D dft to fourier domain, complex multiplication and subsequent backward trafo are performed. Since the input is real (data and gaussian), the output is real too. The output shows the favoured, low-pass filtered features like unsharpening, but one problem occurs:
the absolute values of the output data is several orders of magnitude larger than in the input (i.e. the average intensity of the image is enormously increased), even though I scaled the output with 1/N^2 (where N=1000 is the number of grid points per axis). I feel like this point has sth. to do with the variance of the chosen gaussian.

Thanks in advance for any help on this.
Greetings
J.Norpoth
 
Replies continue below

Recommended for you

Please don't double post... original thread here:
thread236-154028

Dan - Owner
Footwell%20Animation%20Small.gif
 
Status
Not open for further replies.
Back
Top