I have an interferometric image which shows fringes at a certain frequency, basically a dominantly periodic signal overlaying the information that I want to extracthere are the fringes that I want to remove. The size of the image is 1024x4067 pixels. So using matlab I performed fft2(im,n) with n = size of image, giving me this fft-image, with maximum values given at (9,24) and it's symetric (1017,4045). I am totally confused now as of how to interprete this result. Does that mean I have a frequency of 9 in y and 24 in x direction making it some sort of Acos(2pi/9y+2pi/24x)+Bisin(2pi/9y+2pi/24x) wave!? I'm guessing this it totally wrong and stupid, but I'm really confused.. :-/ thanks for your help!
- Remove the mean from the image you linked before analysing it with fft2.
- Use fftshift to get a more readable image.
- Use a logarithmic scale on the 2d spectrum if you don't see anything.
- Please label the x and y axis with normalized frequency.
- Using the insight on the noise spectrum you got from the fft2, design the filters you need to remove the periodic noise.
- You noise isn't stationary withing the image, so your filters won't work awesome.
You could also try the following : Consider your columns as realization of a random vector of size approx 1000 (number of rows). Take your first 100 columns, apply PCA with covariance matric proportional to V*V' with V of size [100*1000], the vertical periodic pattern will be your most powerfull eigen vector. Project your data on the 99 remaining vector. Repeat this process on others 100 columns batches, be carefull with energy normalization.