I'm trying to determine the point-spread-function (PSF) of a certain microscope by a deconvolution process of a measurement of a test pattern (representing the blurred image) and a simulation of this measurement (representing the non-blurred image). The "images" them self are lines out of a microscopic image of the pattern (a height profile so to say), in other words, I'm working in 1D space. (I then want to use the results on 2D problems, just to mention that by the way.)
My question now: Is there any option to utilize the deconvolution functions already given by MATLAB (deconvlucy
, deconvwnr
, deconvreg
...) to get to the PSF (their input demands a blurred image and a PSF normally) or another approach to get the PSF out of these two data sets? I already tried to do a fft
of the "images", dividing them and then doing a ifft
of the division (PSF = ifft((fft(measured))/(fft(sim)))
) but the PSF of this ansatz does not yield the desired result when used in a later deconvolution.