Deconvolution tips

Revision as of 09:07, 23 May 2011 by White (talk | contribs) (add Elizabeth's content and format it a little)

A page for deconvolution tips and how-tos

Compiled by E. Crowell and D. J. White

What the hell is the point spread function (PSF) and why should I care?

The PSF is the image of a point source of light as imaged through a lens (or set of lenses) It is the function (shape) describing the diffraction of light in an optical system. Foir a micrscope systems, it can be derived from processed images of sub-resolution fluorescent beads, or an approximation can be calculated from optical theory equations.

Three different types of PSF images are useful:

  1. The raw measured PSF obtained by imaging sub-resolution size beads using the optical system.
  2. The measured PSF obtained by acquiring and processing multiple z-stacks of sub-resolution size beads.
  3. The theoretical PSF generated by a computer program, after input of values describing the optical system.

It is necessary to measure the PSF of the optical system to verify that no abnormalities are present. Imaging sub-resolution beads is a sensitive method to detect any problems in the optical system. If the Airy rings in the raw PSF are highly asymmetrical, it indicates a problem, such as an objective lens that might be fractured, scratched, misaligned, or simply dirty. Or the coverglass could be at an an angle other than exactly flat.

The effective numerical aperture (NA) of the objective can also be estimated by measuring the angle of the cone of light in the raw PSF. If the raw PSF has the expected characteristics of a good quality optical system, then a theoretical PSF can be used for deconvolution and the time-consuming processing steps to generate a useful measured PSF (see below) might be avoided. However, since no optical system is perfect, a high-quality measured PSF will typically give better results than a theoretical PSF.

A normal measured PSF image consists of a three-dimensional Airy pattern of rings of increasing diameter which radiate out from the spherical bead as one increases or decreases the z-position. The Airy rings should be radially symmetrical circles if the lenses are in good condition. When viewed orthogonally, two cones of light are visible around the fluorescent bead (point source of light). The cones of light will not necessarily be perfectly symmetrical in the z-axis, depending on the quality of the optical system. This asymmetry is one characteristic of the measured PSF which makes it superior to the theoretical PSF for deconvolution applications. A very asymmetric in z PSF is indicative of the presence of spherical aberration, bue to wrong coverglass thickness and/or refractive index mismatch between the sample mounting and lens immersion media.

In order to successfully measure the PSF, the fluorescent object being imaged should behave as a point source of light. For this to be possible, the fluorescent beads should be smaller than the optical resolution of the system (e.g. 100 - 170 nm diameter). The optical resolution = 0.61 . lambda / NA , where lambda is the emission wavelength of the fluorophore and NA is the numerical aperture of the objective. All things being perfect, a 1.4 NA objective has an optical resolution of ~ 250 nm. In optimal conditions, the pixel spacing in the digitized image should be 2.3-3 times smaller, in the range of 80 - 100 nm or a bit less, to fulfill the Nyquist - Shannon sampling criterion and to avoid missing information, oversampling and under sampling.

Fluorescent beads larger than the optical resolution do not behave as point sources and deform the light passing through them due to their refractive index being higher than that of the mounting medium.

Before acquiring z-stacks to generate a measured PSF, one should consider the acquisition conditions that will be routinely used to generate images for deconvolution. The measured PSF image should be sufficiently sampled so that no information is lost, yet must be of dimensions that are equal to (or smaller than) the dimensions of the images to deconvolve. The PSF image must "fit" into the image to deconvolve. The PSF image should be large enough in xy to accomodate the full Airy rings, and large enough in z to make the cones of light visible and complete.

The PSF shape depends on three criteria:

  1. The emission wavelength of the fluorophore
  2. The objective used (magnification and NA)
  3. The refractive index of the objective lens immersion medium

For proper deconvolution, the z-step interval is also an important factor, since the PSF image and the z-stack to be processed should ideally have the same z-step interval. PSF images with a smaller z-step interval than that of the z-stack to be processed can also be used.

Therefore, a distinct PSF must be acquired/created for every combination of these criteria routinely used when acquiring images.

Post-processing of acquired PSF images to generate a measured PSF

A high-quality PSF image is difficult to generate by acquiring images of sub-resolution fluorescent beads. Such small beads seldom contain enough fluorescent dye to yield a high signal-to-noise (SNR) ratio. Therefore, z-stacks should be acquired using a relatively high exposure time (> 200 ms). Caution should be used to avoid increasing the exposure time or lamp power to a level which will cause rapid photobleaching of the bead. If photobleaching occurs during acquisition of the PSF z-stack, it will result in a false PSF due to loss of signal caused by photobleaching.

The beads should be sufficiently diluted so that on average only one bead is visible per field of view. (Neighboring beads will contribute their own light diffraction patterns and deform the measured PSF.) Appropriate beads for measuring PSFs should have a low tendency to aggregate. Nevertheless, with subresolution size beads, one can never know if one is imaging a single bead, or a pair or small cluster of beads. Therefore, it is necessary to image multiple fluorescent beads and to compare the obtained PSFs. The ones that show nicely symmetrical Airy rings represent the useful PSFs. If none of the PSFs obtained shows the expected pattern, there is likely a problem in the optics (see above).

A PSF should be obtained for every combination of z-step intervals, wavelengths, objectives, and immersion medium. However, there is no sense in measuring the PSF of low magnification objectives which give undersampled images, or stainings like DAPI which do not need to be deconvolved.

The last step to generating a useful measured PSF for deconvolution is to subtract the background signal from the stacks so that some pixel values will be equal to zero where there is no real signal. To accomplish this, simply measure the mean intensity of a black region in the stack, outside of the Airy rings. This mean value represents the background signal in the image, which usually originates from the positive offset/bias of a CCD camera chip. Subtract the mean value from all the pixels in the stack (use ImageJ/Process/Math/Subtract...). If this step is omitted, square deformation artifacts may be visible in the deconvolved image.

Obtaining images optimized for deconvolution

It is best to start the z-stack well above the object of interest, and to continue well below, in order to cover the PSF of every interesting point in the sample. If only the "in focus" part of the sample is imaged, the deconvolution result will be ugly since information is missing.

The appropriate z-step interval can be calculated from the Abbe/Rayleigh diffraction limit equation Dz = (lambda * RI) / (NA * NA) where Dz = the ideal z-step interval in nanometers, lambda = the emission wavelength, RI = the refractive index of the mounting medium, and NA = the numerical aperture of the objective. For example, a z-step interval of ~ 350 nm is appropriate for imaging FITC fluorescence using a 1.4 NA objective and a mounting medium with an RI of 1.5.

Always use the same conditions that were used for acquisition of the measured PSF, or generate a theoretical PSF that corresponds to your exact conditions.

How to do the deconvolution

Several plugins for imageJ exist. The most sophisticated are the DeconvolutionLab (BIG - EPFL) and the parallel iterative and spectral deconvolution plugins.


For further very nicely written and informative descriptions of this topic, see http://www.svi.nl/HomePage.