[OTF,PSF,HO,HA] = AtmosphericOTF(f,Type,Model,r0,di,D, ... lambda,alpha,Geom,Atm)
Computes the optical transfer function (OTF) and point spread function (PSF) for ensembleaveraged exposure through Kolmogorov turbulence.
If not specified, the following default values are used in modeling: lambda = 1 micron, D = 1 m, r0 = D/4, di = 1 m, alpha =1.
Parameters 
Description 
f [matrix] 
Spatial frequencies in image (cycles/m) 
Type [scalar] 
Switch [1/0] to calculate long/short exposure OTF. 1  long exposure, 0  short exposure 
Model [string] 
Modeling method for OTF/PSF. Supported models are: 'FRIED'  Fried's theoretical propagation (default), 'GAUSSIAN'  Gaussian approximation to theory 
r0 [scalar] 
(Optional) Spherical wave coherence diameter (m) 
di [scalar] 
(Optional) Image distance from system pupil (m) 
D [scalar] 
(Optional) Aperture diameter (m) 
lambda [scalar] 
(Optional) Mean wavelength of light in image (m) 
alpha [scalar] 
(Optional) Turbulence multiplier 
Geom [struct] 
(Optional) Geometry structure from GeomStruct specifying platform height, target height, downrange. If Geom and Atm are input, r0 is calculated, and di is assumed to be slant range given by the geometry. 
Atm [struct] 
(Optional) Atmospheric structure from AtmStruct. alpha, Geom, and Atm are specified in order to calculate an appropriate r0 value. 
Return Values 
Description 
OTF [matrix] 
Optical transfer function for imaging system in spatial frequency domain 
PSF [matrix] 
Point spread function for imaging system in spatial domain 
HO [matrix] 
OTF component due to optics only 
HA [matrix] 
OTF component due to atmosphere only 
>> f = linspace(0,1,100);
>> [OTF,PSF] = AtmosphericOTF(f/550e9,1,'FRIED',0.18,1,1,550e9)
>> [OTF,PSF] = AtmosphericOTF(f/550e9,0,'FRIED',1,1,550e9,Atm)
>> [OTF,PSF] = AtmosphericOTF(f/550e9/Atm.L,1,'GAUSSIAN',1,G,Atm)
Copyright (c) 2009. All rights reserved.

What do you think about this topic? Send feedback!
