Image restoration principleIf everything was perfect (perfect optics - perfect seeing - perfect telescope tracking) the image from a star would be a single pixel on the CCD, like this:
In fact, because of all imperfections, the observed star image is spread over several pixels. This is known as the Point Spread Function (PSF). A star image, together with the star luminosity function along horizontal axis typically looks like this:
The transform between the expected image (the single pixel) and the observed image (the PSF) being known, one could expect that applying the invert transform would result in a perfectly restored image. This is not true because of the noise in the image, which will be strongly emphasized by the invert transform. Taking an analogy in the audio field, suppose you make a tape recording with treble control at minimum. When replaying the tape, if you attempt to "restore" it by setting the treble control at maximum, you will get hiss from the tape. This is why iterative restoration algorithms have been developped. At each pass, they tend to ameliorate the PSF towards a single pixel. When the best compromise between image detail enhancement and noise has been reached, the iterations should be stopped. Two image restoration algorithms are widely used in astronomy: Lucy-Richardson and Maximum Entropy. Lucy-Richardson is a linear algorithm, in the sense that it equally restores high and low luminosity portions of the image. Maximum Entropy first restores the high luminosity portions of
the image then, on successive passes, the low luminosity portions. On deep
sky images, this means that bright stars are first restored then the core
details of the object. The following experiments use non-astronomical images. A reference image is used with an added single pixel "perfect" star. This image is "deteriorated" in several ways, which transforms the perfect star into a PSF. Then, using the PSF, attempts are done to restore the image.
ConclusionImage restoration can be a powerfull tool to restore image reality if used properly. The big enemy to a good restoration is random image noise. With astronomy images, noise can be reduced by longer exposures, dark frame subtraction, lower CCD temperature and averaging of multiple exposures. Low noise images which exhibit some blurr due to optic limitations and/or poor seeing can be efficiently restored using Maximum Entropy. Low noise images with poor guiding can be efficiently restored by using Lucy-Richardson deconvolution.
|