Apply shift-variant point spread function to irradiance image in oi outIrrad = rtApplyPSF(oi,PSF_ARRAY_ALL) rtApplyPSF performed PSF rotation and interpolation to optical image resolution dynamically for each pixel's angular location. To speed up computation: For each each sample wavelength, we have a set of rotated and sampled PSFs are computed and passed in. These are interpolated to 1 deg by rtPrecomputePSF (can this be an arbitrary spacing?) A pixel at an angle between indexes of the look up table is linearly interpolated between the floor and ceiling indices of the PSF LUT. Radial interpolation is then performed between the angularly interpolated PSFs as in rtApplyPSF. Examples Copyright ImagEval, LLC, 2005