humanLSF

 Calculate the human linespread function

   [lsf,xDim,wave] = humanLSF([pupilRadius=0.0015],[dioptricPower=59.9404],[unit='degree']);

 The human line spread function includes optical defocus and chromatic
 aberration. The spatial extent (and spatial frequency) range are
 determined by the spatial extent and sampling density of the original
 scene.

 The returned units can be degrees (default) or mm or um (unit = 'mm',
 unit = 'um').
  
 See also:  humanOTF and the discussion therein.

 Reference:
   Marimont & Wandell (1994 --  J. Opt. Soc. Amer. A,  v. 11,
   p. 3113-3122 -- see also Foundations of Vision by Wandell, 1995.

 Examples:
  [lsf,xDim,wave]  = humanLSF; colormap(jet); mesh(xDim,wave,lsf);
  xlabel('wave'); ylabel('mm')
  [lsf,xDim,wave]  = humanLSF([],[],[],'mm'); colormap(jet); mesh(xDim,wave,lsf)

  radius = 0.003/2;  % In meters
  dioptricPower = 60;% In diopters (1/m)
  unit = 'mm';
  [lsf,xDim,wave]  = humanLSF([],radius,dioptricPower,'mm'); colormap(jet); mesh(xDim,wave,lsf)

 Copyright ImagEval Consultants, LLC, 2005.