Exhaustive Photon Restoration
Exhaustive Photon Restoration (EPR) is BIG's 3D Deconvolution routine. It uses a point spread function (PSF) to deconvolve images from a microscope. The PSF is a 3D image of a point source, usually a small fluorescent bead. The PSF should be acquired under similar conditions as the images that will be deconvolved. The fourier transform of the PSF represents the optical transfer function (OTF) of a microscope and is used to deconvolve images in fourier space.
There are three main programs involved in 3D Deconvolution: preppsf, prepdata, and epr_i2i. They all use our own i2i image format. We use a program called mmtoi2i to convert µManager images acquired on TESM to i2i files.
preppsf prepares the psf for epr_i2i.
prepdata prepares the images for epr_i2i.
epr (or super_epr or gpu_epr) performs the 3D deconvolution.
An account is needed with the Biomedical Imaging Group to access these programs.
3D Deconvolution Steps
These steps need to done on a BIG's workstation or server. I recommend Mizar.umassmed.edu because it has a quad core CPU, and 32GB of memory. For GPU deconvolutions, Oxygen.umassmed.edu is available with permission from Karl.
The following is a list of example steps.
- ssh (or putty on Windows) into mizar.umassmed.edu or oxygen.umassmed.edu (with permission)
- ssh -l <username> mizar.umassmed.edu
- Change to the parent directory of where your images are located.
- cd /storage/big1/<username>/tesm/
- List all the directories
- Convert µManager PSF and Image tiff files to i2i format.
- mmtoi2i <directory>
- Prepare PSF: #preppsf
- -NA number
- NA of the objective,
- -n number
- The refractive index of the immersion medium
- -size number
- pixel size at the object in micrometers/pixel
- -spaced number
- Distance spacing among Z planes in micrometers
- -norm parameter
- Normalize total intensity method (infocus or value)
- -wavelength number
- Emission wavelength of light in micrometers
- preppsf -NA 1.4 -n 1.515 -size 0.187 -spaced 0.25 -norm infocus -mask 1 psf.i2i psf_prep.i2i
- Prepare Images: #predata
- Run Deconvolution Software: #epr_i2i
- This is a default version of EPR that runs on a CPUs.
- This is an EPR script wrapper that enables super-resolution 3D deconvolution. It is not GPU accelerated.
- This is the newest version of EPR that enables GPU acceleration for faster 3D Deconvolution. It is currently being updated to run using newest CUDA libraries. This version will be available soon.
All versions of EPR need a 3D PSF image, the 3D image to be deconvolved, and where the 3D Deconvolved image should be saved. All in i2i format.
- epr -psf psf_prep.i2i -smoothness 0.0005 -iterations 250 mycell_prep.i2i mycell_prep_r.i2i
Values for Convergence and Iterations
Two options for EPR are convergence (-co <number>) and iterations (-it <number>).
The number of iterations is the maximum number of iterations that EPR will run if it hasn't reached convergence.
- -it 300
Convergence is used to determine when to stop the deconvolution. The range for -co can be between 1 and 0.0001. 0.0001 = true convergence (100%), while 0.001 achieves 90-95% convergence. The smaller the value you use for -co equates to longer times for the deconvolution.
- -co 0.001
Unlike epr and gpu_epr, super_epr only allows setting -sm. Super_epr uses "-co 0.001" and "-it 300" internally.
Values for Smoothness
The most significant option for EPR is smoothness (-sm <number>).
The value passed as smoothness is referred to as alpha. Alpha is usually between 0 and RNL (Residual Noise Limit). RNL is reported by prepdata if you know the gain and read specifications of the camera used to acquire the images (see -rnl option for prepdata). These values can be measured empirically. A smaller alpha correspond to less smoothing. (RNL)^2 is usually a good starting choice for alpha. You do not want to overconvolve your images by being too aggressive with alpha (smaller values).
Run the command for full option list: preppsf -h
preppsf - prepares a 3-D point spread function (psf) image data set, acquired using a Digital Imaging (light) Microscope, for use with image restoration/3-D reconstruction. Normally psf images are first processed using prepdata (see below) to apply basic corrections and any background corrections, but not temporal corrections.
preppsf [options] before-image after-image
preppsf performs further corrections on a psf image data set based on the fluorescence D.I.M image formation/acquisition model. The before- image argument is the name of a 3-D psf image set, acquired (DIM-1, DIM-2, CELLscan, UFM, or other) and processed using prepdata (see below) without the -norm option. The psf is extracted as a symmetric sub-region (square in XY) centered at the psf origin (-center), optionally normalized for constant total intensity (-norm), and optionally masked (-mask) to exclude any extramural data. The after-image argument is the file name for the processed psf image set. All options must precede input/output image files
- Numerical Aperture of the objective.
- Index of refraction of immersion medium.
- Effective Pixel Size at the object in micrometers/pixel
- DIstance Spacing among Z planes in micrometers.
preppsf -NA 1.4 -n 1.515 -size 0.187 -spaced 0.25 -norm infocus -mask 1 psf_in psf_out
Run the command for full option list: prepdata -h
prepdata - prepares a 2-D, 3-D or 4-D image data sets, acquired using a Digital Imaging (light) Microscope, for further processing (e.g. image restoration/3-D reconstruction)
prepdata [options] before-image after-image
prepdata performs all necessary corrections on a image data set based on the fluorescence D.I.M image formation/acquisition model, and can be used to perform basic imaging corrections to many forms of digitally acquired image data. The before-image argument is the name of a 3-D image set as acquired (DIM-1, DIM-2, CELLscan, UFM, or other). The after-image argument is the file name for the corrected image set, ready for further processing. The normal order of application of corrections (see -before) is basic[->background[->temporal]]. All calculations are performed using single-precision floating-point arithmetic. All options must precede input and output image file names.
Run the command for full option list: epr_i2i -h
epr_i2i - Exhaustive Photon Replacement (EPR) restores contrast by removing residual out-of-focus light and improves resolution while maintaining numerical accuracy of 3-D images of specimens obtained with serial optical sectioning from wide-field or confocal light microscopy.
epr_i2i [options] before-image[.i2i] after-image[.i2i]
epr performs regularized, iterative image restoration with a non- negativity constraint. The before-image is a three dimensional (3-D) TIFF image composed of rectangular, regularly spaced, optical sections. Large images are decomposed into smaller, overlapping (in x and y only) image segments for restoration, and the restored segments are recomposed. The after-image is the restored 3-D image. All options must appear before the image file names, but the order of options is not important.
epr -psf mypsf_ -smoothness 0.0005 -iterations 250 mycell_ mycell_r