Creating Your Own Deblurring Functions

Creating Your Own Deblurring Functions

All the toolbox deblurring functions perform deconvolution in the frequency domain, where the process becomes a simple matrix multiplication. To work in the frequency domain, the deblurring functions must convert the PSF you provide into an optical transfer function (OTF), using the psf2otf function. The toolbox also provides a function to convert an OTF into a PSF, otf2psf. The toolbox makes these functions available in case you want to create your own deblurring functions. (In addition, to aid this conversion between PSFs and OTFs, the toolbox also makes the padding function padarray available.)

13-23

13 Image Deblurring

Avoiding Ringing in Deblurred Images

The discrete Fourier transform (DFT), used by the deblurring functions, assumes that the frequency pattern of an image is periodic. This assumption creates a high-frequency drop-off at the edges of images. In the figure, the shaded area represents the actual extent of the image; the unshaded area represents the assumed periodicity.

This high-frequency drop-off can create an effect called boundary related ringing in deblurred images. In this figure, note the horizontal and vertical patterns in the image.

To avoid ringing, use the edgetaper function to preprocess your images before passing them to the deblurring functions. The edgetaper function removes the high-frequency drop-off at the edge of an image by blurring the entire image and then replacing the center pixels of the blurred image with the original image. In this way, the edges of the image taper off to a lower frequency.

13-24