0
$\begingroup$
  1. I am imaging a sample using CMOS camera. The pixel ratio is 6um x 6um and a total resolution {480, 752}. I understand that each pixel on the camera sensor is 6um x 6um in size (Have I understood this correctly ?).
  2. If I have to find the signal-noise ratio of the 2D image to determine the image quality. In doing so I have to convert it into Fourier domain using FT. How can I know the sampling interval to calculate the FFT of the image data ?. Does the sampling interval depends on the pixel size ?
  3. Suppose this is the Power spectral density of the image. How can I know the noise level and signal level from the figure 2?. I have assumed the sampling interval to be DelT =[(2*Pi)/PixelSize]

Image : Image PSD: PSD

$\endgroup$
1

1 Answer 1

1
$\begingroup$

Q1. You are right.

Q2. The sampling interval in both directions is $6{\rm \mu\,m}$. This yields the peak frequency of the Fourier transform in each case: namely $2\pi/(6\times 10^6)\,{\rm rad\,m^{-1}}$. Note that the Nyquist frequency is half this: the signal is uniquely determined by the samples as long as it is bandlimited to below the Nyquist frequency: the frequency "folds" around the Nyquist frequency so that the FT points above the Nyquist frequency are actually the superposition weights for negative frequency Fourier components. For real valued signals like yours, these negative frequency weights are complex conjugates of the corresponding positive frequency ones, so that they are redundant information. The resolution of the Fourier transform gives frequency spacing. You will almost certainly be doing a radix 2 fast Fourier transform, so you'll need to pad your image out with noughts to make its resolution $2^{N_x}\times 2^{N_y}$ i.e. the resolution components are whole number powers of two. If you have software that doesn't need you to do this, it will almost certainly be doing this automatically for you. So let's assume you embed the image centrally in either a $512\times 1024$ or $1024\times 1024$ complex dataset. In the first case, your frequency samples will be spaced at $2\pi/(512\times6\times 10^6)\,{\rm rad\,m^{-1}}$ in the 512 samples direction and $2\pi/(1024\times6\times 10^6)\,{\rm rad\,m^{-1}}$ in the other. In the $1024\times1024$ case, the frequency spacing is $2\pi/(1024\times6\times 10^6)\,{\rm rad\,m^{-1}}$ in both directions.

Q3. This question has no definite answer. I would use an "experimental" approach. Set up a signal processing sequence that takes the FFT, set all the components at spatial frequencies higher that a cutoff value to nought and transform back with an inverse FFT. You'll see a smoothed image. Try lowering the cutoff frequency and see how far you can go and still keep the image "understandable", i.e. sharply enough resolved to make out the "meaningful" details. (Ask yourself, "Does the image still mean the same thing for practical purposes". From your frequency spread graphs, it looks as though this bandwidth will be somewhere between about 0.05 and 0.1 on your horizontal axis: this is the bandwidth above which the Fourier weights are between 2 and 3 orders of magnitude smaller than the lowest frequency components.

$\endgroup$
2
  • $\begingroup$ awesome explaination. But why is the sampling interval 6*10^6 ? Should it not be 6*10^-6 ? $\endgroup$ Commented Nov 3, 2013 at 18:07
  • $\begingroup$ @abhilashsukumari yes, that's a typo! $\endgroup$ Commented Nov 3, 2013 at 18:30

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.