Compute the human optical transfer function [otf, williamsFactor] = humanCore(wave,nWave,sampleSF,p,D0) The computation of the human OTF is performed here. It is placed in this routine to permit compilation. The OTF is computed at the wavelengths in wave and at sampleSF spatial frequencies in units of cycles/deg The returned williamsFactor is a general loss of contrast, independent of wavelength, that we apply to bring the data into alignment with human optics. That factor comes from empirical measurements out of the Williams' lab. p: Pupil radius in meters D0: Dioptric power (accomodation), usually around 60 sampleSF: Spatial frequencies in cycles/deg wave: Wavelength in nanometers Example: wave = 400:10:700; sampleSF = [0:1:30]; p = 0.0015; D0 = 60; otf = humanCore(wave,sampleSF,p,D0); mesh(sampleSF,wave,otf) Copyright ImagEval Consultants, LLC, 2005.