Institut für Astronomie und AstrophysikAbteilung AstronomieWaldhäuser Str. 64, D-72076 Tübingen, Germany |
MAX_ENTROPY
Deconvolution of data by Maximum Entropy analysis, given the PSF
Deconvolution of data by Maximum Entropy analysis, given the instrument point spread response function (spatially invariant psf). Data can be an observed image or spectrum, result is always positive. Default is convolutions using FFT (faster when image size = power of 2).
for i=1,Niter do begin Max_Entropy, image_data, psf, image_deconv, multipliers, FT_PSF=psf_ft
data = observed image or spectrum, should be mostly positive, with mean sky (background) near zero. psf = Point Spread Function of instrument (response to point source, must sum to unity). deconv = result of previous call to Max_Entropy, multipliers = the Lagrange multipliers of max.entropy theory (on first call, set = 0, giving flat first result).
deconv = deconvolution result of one more iteration by Max_Entropy. multipliers = the Lagrange multipliers saved for next iteration.
FT_PSF = passes (out/in) the Fourier transform of the PSF, so that it can be reused for the next time procedure is called, /NO_FT overrides the use of FFT, using the IDL function convol() instead. /LINEAR switches to Linear convergence mode, much slower than the default Logarithmic convergence mode. LOGMIN = minimum value constraint for taking Logarithms (default=1.e-9).
function convolve( image, psf ) for convolutions using FFT or otherwise.
Iteration with PSF to maximize entropy of solution image with constraint that the solution convolved with PSF fits data image. Based on paper by Hollis, Dorband, Yusef-Zadeh, Ap.J. Feb.1992, which refers to Agmon, Alhassid, Levine, J.Comp.Phys. 1979. A more elaborate image deconvolution program using maximum entropy is available at http://sohowww.nascom.nasa.gov/solarsoft/gen/idl/image/image_deconvolve.pro
written by Frank Varosi at NASA/GSFC, 1992. Converted to IDL V5.0 W. Landsman September 1997
[Home Page] [Software, Documentation] [IDL Documentation] [Quick Reference] [Feedback]