Statistics of A Geometric Representation

Download as pdf or txt
Download as pdf or txt
You are on page 1of 9

VOLUME 55, NUMBER 11 NOVEMBER 1965

JOURNAL OF THE OPTICAL SOCIETY OF AMERICA

Statistics of a Geometric Representation of Wavefront Distortion


D. L. FRIED
Electro-Optical Laboratory, North American Aviation, Inc., Space and Information Systems Division,
Torrance, California
(Received 26 February 1965)

A precise statistical definition is established for the geometric "shape" of a randomly distorted wavefront.
Relationships between the phase-structure function and the statistics governing the shape are derived. The
most significant portion of wavefront distortion caused by atmospheric turbulence is a random tilting of the
plane-wave front. A procedure is outlined for calculating the influence of wavefront distortion on optical
systems. Estimates are formed of the effect of wave-front distortion on photographic resolution and optical
heterodyne efficiency.

I. INTRODUCTION ful glimpse into the problem, providing quantitative


results, can be obtained if we have data on the statistical
OPTICAL wave propagation through a randomly "shape" of wavefront deformation. Similarly, the per-
inhomogeneous medium produces wavefront de- formance of an optical heterodyne receiver detecting a
formation. The spatial (as distinguished from the signal with a deformed wavefront can be estimated
temporal) statistics of this deformation are generally from wavefront "shape." 4
described by a quantity called the phase-structure func- The key to the whole concept is in formulating a
tion.' If, at points x and x' ,the phase variations associ- precise statistical definition of the deformed wavefront's
acted with the deformation are denoted by q(x) and "shape." In this paper, this has been accomplished by
¢(x'), respectively, then the phase-structure function is examining the wavefront over the circular aperture of
defined as an unspecified optical system. The deformation in that
ED(r)= ([+(x)-k(X')] 2 ), (1.1) region is represented by a series of orthonormal poly-
where nomials. These polynomials are chosen so that each
r=jx-x'j, (1.2) term or group of terms in the series can be associated
with a specific geometric shape. The shape of the wave-
and the brackets ( ) denote an ensemble average. front can then be related to the mean square value of
(Isotropy and homogeneity of the deformation statistics
the coefficients of the polynomials in the series. The
are implicit in this definition.) Theoretical and experi- mean square error for a truncated series provides a
mental procedures for the evaluation of D(r) have been
developed. 2
measure of the shape of the deformed wavefront.
Inasmuch as D(r) represents all the statistics of phase The polynomials used are the first six terms of an
fluctuation,' it should be possible to obtain information infinite sequence. Their sum can be cast in the form of
concerning the "shape" of the deformed wavefront a general quadratic; i.e.,
from D(r). Knowledge of the shape would provide
insight into the nature of the effects of wavefront de- a+bx+cy+dx2 +exy+fy2 .
formation on an optical system. For example, if wave
deformation can be approximated by a random tilting This can represent any wavefront distortion which can
of a plane wavefront, then shorter exposures will be built up from an average phase change, a tilt, a
clearly result in improved photographic resolution. The spherical deformation, and a hyperbolic deformation.
resolution can be computed directly from the phase- The extension of the results obtained in this paper to
structure function, 4 but a quick and physically meaning- include higher order types of distortion is a straight-
forward process.

1V. I. Tatarski, Wave Propagation in a Turbulent Medium II. NOTATION


(McGraw-Hill Book Company, Inc., New York, 1961).
2 See Ref. 1. R. E. Hufnagel and N. R. Stanley, J. Opt. Soc. Am. Consider a plane perpendicular to the nominal direc-
54, 52 (1964). D. L. Fried and J. D. Cloud, J. Opt. Soc. Am. 54,
574A (1964). tion of propagation of a distorted plane wave. Two
3It can be shown theoretically that phase fluctuation, because position vectors, x and x' are defined in terms of a co-
of the central limit theorem, has a Gaussian distribution. Since a
Gaussian distribution is completely described (except for a mean ordinate system in this plane. The vectors x and x' have
value) by its second moment, and since OD (r) is the second moment components (xy) and (x',y'), and magnitudes x and x'.
for differential-phase fluctuation (which difference has zero mean),
we conclude that D)(r) completely specifies the statistics of phase As indicated previously, the phase at x and x' is denoted
fluctuation. by k(x) or +(x'), for a particular realization of the
ID. L. Fried, "The Effect of Wave-Front Distortion on the
Performance of an Ideal Optical Heterodyne Receiver and an wave-front distortion. Brackets ( ) are used to denote
Ideal Camera," presented at the Conference on Atmospheric an ensemble average. The function W(x,D) defines a
Limitations to Optical Propagation at the U. S. National Bureau circular region of diameter D, radius R. It equals unity
of Standards CRPL, 18-19 March 1965.
1427
1428 D. L. FRIED Vol. 55

inside this circle and zero outside; i.e., +(x); i.e.,


A/\ (1/7rR2) dxW(xD)[0(x)- c(x) ]2 (3.2)
IV(xD)-={' if 1x1<R (2.1)
0, if !x>R.
To minimize A, a,, is chosen so that
The six orthonormal polynomials worked with in this
paper are defined in Eqs. (2.2a- f). The first polynomial (a/a,=)/A = 0. (3.3)

F,(x) - (7rR)- ~ ~ 2.2a) It is shown in Sec. IV that extrapolation to two-


dimensional form from the usual one-dimensional result
can be used to represent a change in the average phase for the coefficients of an orthonormal series is applicable;
over the circular region. The second pair of polynomials i.e., that
F 2(X)-(7rR4/4)-'(x) (2.2b) a,,= dxW(xD),0(x)F,(x). (4.3)
Fa3(X)- (7rR414)-1(y) (2.2c)
New coefficients ac, aL, as, and aQ, are defined by
can be used to represent an average tilt. The next
(aC)2 - (a 1)2, (3.4a)
F 4(x)= (7rR1/12)1(x 2+y 2 -R2/2) (2.2d)
(aL)2= (a2)2+ (a3)2, (3.4b)
can be used to represent a spherical deformation. Finally,
2
(as)2= (a4) (3.4c)
F5(X)-(7rR'/6)-1(X2-y2) (2.2e) 2
(aQ) = (a4)2+ (a5) 2+ (a6)2. (3.4d)
Fo(x) (7rRG/24)1(Xy), (2.2f)
Examining the polynomials Fl(x) to F 6(x) and visu-
can be used to represent a hyperbolic deformation. alizing the surface each represent, we see that the
These polynomials are closely related to Zernike's quantities ((ac) 2), ((aL)2), ((aS) 2), and ((aQ)2) can be
polynomials. It is easy to generate higher-order terms; considered measures of the mean square amount of
i.e., F 7(x) and up (from Table XXI of Ref. 5). It is easy average phase fluctuation over the circular region of
to verify the fact that these polynomials satisfy the interest, of the average tilt of the wavefront, of the
orthonormality condition extent of spherical deformation of the wavefront, and
of the amount of quadratic (i.e., spherical plus astig-
JdxW(x,D)F,(x)F,(x) = b, (2.3)
matic) deformation, respectively.
In addition to evaluating these four quantities, it
would be desirable to know how good an approximation
where 6,,, is the Kronecker delta to the deformed wavefront is provided by a finite
number of terms in the series cb(x) of Eq. (3.1). For this
t1 if bu=P purpose, define the finite series Ic(x), 4L(X), Cbs(x)
lo { if U :. (2.4) and cDQ(x), as

The integration in Eq. (2.3) and all other integrations =1j (3.5)
throughout this paper are understood to cover the
infinite plane, unless explicitly indicated otherwise. where j takes the values C, L, S, and Q and nj has the
Finite bounds for the integraion are provided, in effect, values 1, 3, 4, and 6, respectively. The geometric
by the aperture function. significance of the four Ij's defined above is obvious.
Now, let us define Ai as
III. WAVEFRONT APPROXIMATION
The distorted wavefront's phase 4 (x) is approximated Aj-=(l/rR2) dxW(xD)[0(x)-cPj(X)]2. (3.6)
by the infinite series

(D(x)= E a,,F.(x), (3.1) If the ensemble average (Aj) can be computed, we obtain
A1=1 a measure of the goodness of various possible finite
where the coefficients a. are chosen to optimize the series approximations to the distorted wave front.
approximation. The error in the approximation is A, the
aperture-averaged-square difference between >(x) and IV. DERIVATION OF UNAVERAGED
RELATIONSHIPS

5 M. Born and E. Wolf, Principlesof Optics (Pergamon Press, Substituting for /j(x) in Eq. (3.3) by using Eq. (3.1),
Inc., New York, 1959). expanding the square, and interchanging the order of
November 1965 STATISTICS OF WAVEFRONT DISTORTION 1429

summation and integration, we find From Eqs. (3.4b, c, and d) and Eq. (4.4') we can see that
2
(aL) == AC-AL, (4.6a)
dxW(x,D)q2 (x)
A=-
(as) 2 = AL- AS, (4.6b)
2 RJ (aQ) 2 = AL-AQ. (4.6c)
- aa,, dxW(x,D)O(x)F,(x)
=rR2 J The problem at this point is reduced to computation of
the ensemble average of Aj as given by Eq. (4.5') and
+ A aY a- b dxW(x,D)F,(x)F^(x). (4.1) relating the result to the phase-structure function.
$C=1 7rR 2
V. DERIVATION OF AVERAGED
Invoking the orthonormality condition, Eq. (2.3), we RELATIONSHIPS
can perform the last integration and then one of the FE~amining Eqs. (2.2a-f), we note that FI(x) is inde-
two sum~mations in the double sum, we obtain
pendent of x so that Jfdx'W(x',D)F 1 (x)F1 (x') equals
11- unity. Also for g not equal to one, JfdxW(x',D)F,(x)
'dxW(x,D)02(x) XF,(x') vanishes. Consequently,
7rR 2 J
2 .0 1 1
- 2
A a,; dxTW(x,D)Ok(x)F,,(x) - fdxW(x,D)k2(x)= ffdxdx'W(x,D)
,7rR2 '_1 7rR2 J rR 2 J
1 0
+ 2
E_ (a,,E2~. (4.2) XTV(x',D)Y F,,.(x)F,,(X'),02(X). (5.1)
7rR g=1 i1-
Substituting Eq. (4.2) into Eq. (3.2), differentiating, Symmetrizing between x and x' in the right-hand side
and solving for a,, we get of Eq. (5.1) gives

a, = dxW(x,D)k(x)F,(x), (4.3) fdxW(x,D),02(x)=


7rR2 J27rR2
f J
dxdx'W(x,D)

as indicated previously. Substituting Eq. (4.3) into


(4.2) to eliminate the second integration, we get XW(x',D), F(x)F,,(XI) [t,2(X) +02(X') ].
F, (5.2)
Y=1
Ir
1
A= |dxW(x,D),02(X)__YE (a,)2. (4.4 Substituting Eq. (5.2) into (4.5'), taking ensemble aver-
2
R2 J7rR 2=1 ages of both sides, converting the variables of integra-
tion from x, x' to r, r', where r, r' are defined as
If we substitute Eq. (4.3) into (4.2) to eliminate a,,,
make double integrals out of the product of two rf= 2 (X+X'), (5.3a)
integrals, and interchange the order of integration and r= x-x', (5.3b)
summation, we obtain
and recognizing the presence of the phase-structure
2
-dxW(x,D)0 (x)-- Ildxdx'W(x,D) function in the integrand, we find
irR 2 J r12 J
(Aj)=- ffdrdr'W( r'+ rj, D)W( r'-2rj, D)
X W(x',D) YF,(x)F.(x')0f(x)0f(x'). (4.5)
n?

Working with the finite series approximation, we can XE F,l(r'+ r)F,,,(r'- 1r) O(r) . (5.4)
8t=1
easily show the results corresponding to Eq. (4.4) and
(4.5) to be The r' integration in Eq. (5.4) can be performed ex-
plicitly. For this purpose, define the function 5#(r,D) as
1 1- 1 nlj
A= IdxW(x,D)0 2 (x)- Y(a.)2 (4.4')
7rR2 J 7rR 2
g=1
5#(r,D)= fdr'W( Ir'+±rf,D)W( r'-1rj,D)
A=-
1 1.1C
dxW(x,D)02(x)-- J dxdx'W(x,D) nX
XA F,(r'Jr r)F,,(r'- Ir). (5.5)
pull

X W(x',D)E F,,(x)F,,(x')0(x)0(x'). (4.5')


i=l This function is evaluated in Appendix A for j= C, L,
1430 D. L. FRIED Vol. 55

S, and Q. The results are ignore this factor. (A discussion of this matter is given
1 in Appendix D.)
5Fc(r,D)=(1/7r){2 cos- (r/D) In this paper, particular values of -A are not of con-
-2(r/D)[1-(r/D) 2 ]i} W(r,2D), (5.6a) cern. (The necessary data for computing A are provided
9L(r,D) = (1/7r){6 cos-'(r/D)- [14(r/D)-8(r/D)3 ] in Appendix C.) In fact, it is convenient to replace -A
with a new quantity r0 which has the dimensions of
X E1- (r/D)2]} W(r,2D), (5.6b) length and is defined by
gs(r,D)= (1/7r) (8 cos-1 (r1D) rOc-(6.88/eA)3 5. (6.2)
-E24(r/D)-(80/3)(r/D)3 +(32/3)(r/LD) 5 ]
X[1-(r/D)2 ]1 }2W(r,2D), (5.6c) Correspondingly, the phase structure function, ex-
pressed in terms of r0, is
9Q(r,D)= (l)1 2 cos-1 (r1D)
D(r) = 6.88(r/ro)5 13 . (6.3)
-E44(r/D)-64(r/D)+±32(r/D) 5]
X 1-(r/D)2 ]} W(r,2D). (5.6d) The apparently arbitrary constant 6.88 was chosen on
the basis of the analysis of the performance of an optical
The evaluation of 5:j(r,D) given in Eqs. (5.6a-d), ob- heterodyne detection system.' In Ref. 4, ro is shown to
tained in the Appendix through a rather extensive and be that diameter of a heterodyne collector for which
quite uninteresting calculation, is fundamental to the distortion effects begin to seriously limit performance.
program of this paper. With this in hand, the rest of the It is seen later that it is also that diameter for which
work is straightforward. Substituting from Eq. (5.5) (Ac) is essentially unity. Typical values of ro for visible
into Eq. (5.4) and noting that the integrand is isotropic and near infrared wavelengths and for approximately
in r, so that the angular integration can be performed, vertical propagation paths down through the atmos-
we find phere are of the order of several centimeters.
(Aj) =
1
f
rD
rdrgj(r,D)D(r), (5.7) VII. EVALUATIONS
Substituting Eq. (6.6) into (5.7), replacing the
1 D
variable r by it, where
((aL )2) = - rdr[9c(r,D)- 5L(r,D)]5t(2 ), (5.8a)
u= r/D, (7.1)
rD
and noting that
(
((as)2) =- rdr[5L(r,D)-5is(r,D)]°)(r), (5-8b) 3j(r,D) = 5j(ul), (7.2)
we get
1 rD (Aj)=27.5(D/ro)513 udu~j(u,1)u5 /3 . (7.3)
((aQ)2) = R J rdr[IEL(r,D)- 5:Q(r,D)]D(r). (5.8c)

Utilizing the known relations 7


VI. THE PHASE-STRUCTURE FUNCTION
Several theoretical studies of the phase-structure JZ2,,+l(l-Z2)Odz= 'B(a+1, 0+1), (7.4a)
function have been performed utilizing various ap-
proaches to the propagation problem. 2 Based on the
IKolomogoroff similarity theory of turbulence, 6 which
predicts a spatial correlation of turbulence which de-
creases proportionally to the two-thirds power of the
fZe2(re
Cos-1(z)dz= a
1)
B(- 2)
2 2
(7.4b)

spatial separation, it is possible to show that the phase where B(a,O) is the well-known beta function, and
structure function may be written as defining Ij as
5)(r) =,_*r5/3 (6.1) _Tj = it"15j(ul)du, (7.5)
whereXA is a constant determined by the path of propa-
gation, the wavelength, and the particular environ- we find that
mental conditions. Equation (6.1) is exact only in the Ic!3.68X 10-2, (7.6a)
near field; i.e., for short propagation paths. For longer IL-4. 73
X 10-3, (7.6b)
distances of propagation, there is an additional r-de-
pendent factor which varies from one half to one, which IS-3.96X 10-3, (7.6c)
is suppressed in (6.1). For the body of this paper, we IQc2.29X 10-. (7.6d)
61A. Kolomogoroff in Turbulence, Classic Papers on Statistical
Thzeory, edited by S. K. Friedlander and L. Topper (Interscience 7 W. Grobner and N. Hofreiter, Integraltafel (Springer-Verlag,
Publishers, Inc., New York, 1961), p. 151. Berlin/Vienna, 1961), Vol II, Eqs. (121.1) and (341.5a).
November 1965 STATISTICS OF WAVEFRONT DISTORTION 1431

Thus, Since the coefficient of the linear term is so much larger


(Ac)>1.013(D/ro) 5 13 , (7.7a) than that of the spherical or other quadratic terms, we
513 (7.7b) can conclude that a large part of the deformation con-
(AL)c0.1301 (D/ro) ,
*b sists of wavefront tilting. We might have expected this
5 3
(As)-0.1090(D/ro) /, (7.7c) result from the fact that, as proven in Appendix B,
6 3
(AQ)Ž0.0630(D/r,) I . (7.7d) 1 D R2/2 if j=C (8.1)
- Ydr5~j#,D)r6ll=
From these results we can see that R2J 0 if jC5 (8.1)
((aL) 2)c0.883(D/r0)5 13 , (7.8a)
Considering this in conjunction with Eqs. (5.7) and
((aS)2)-o.02l11(D/ro) 131, (7.8b) (5.8a-c), we see that a phase structure function which
((aQ) 2
)-o.0671(D/ro) 5 8
/ . (7.8c) has a six-thirds dependence on r [instead of the five-
thirds power of Eq. (6.1)], corresponds exactly to a
It is principally through ((aj)2 ) and (Z\j) that a precise randomly tilted, but otherwise undistorted, plane wave.
definition of the deformed wave front's "shape" is Noting that the spherical coefficient ((as) 2 ) corre-
developed. A quantitative physical insight into the sponds to a single degree of freedom in the series ap-
nature of the deformation can be obtained, however, by proximation, that the quadratic coefficient ((aQ)2)
computing the quantity Dj*. The quantity Dj* is the corresponds to three degrees of freedom, and noting
diameter of the circular region for which (As) takes on further that ((as)2)/((aQ)2) as determined from Eqs.
some critical value A*. Thus, it is the aperture diameter (7.8b c) is almost exactly one-third, we can conclude
over which the average deformation may be considered that spherical deformation of the wavefront is no more
to consist of no significant amount of deformations of and no less significant than the other two quadratic
higher order than j. What constitutes a significant forms of deformation.
amount is defined by A\*. The aperture diameter Dj* is Considering Eqs. (7.12a-d), we see that the area per
defined by Eq. (7.9). degree of freedom over which a series approximation
can provide a given quality match (as specified by A*)
4 f dxW(xDj*)E[(x) - 4j(x)]2> (7.9) to the wavefront jumps abruptly between the one-term
j= C approximation, and the three-term j= L approxi-
mation. If this jump had not occurred, the conclusion
The corresponding aperture area Aj* is simply would have to be drawn that the wavefront is distorted
A j*_ (7 r/4)(Dj*) 2 . (7.10) in a way which is not subject to a geometric interpreta-
tion. In fact, quite the contrary is the case.
The quantity Aj*lnj is of interest in that it provides a To get some insight into the photographic resolution
measure of the "average utility" of each degree of to be expected when the wavefront being collected is
freedom represented by one of the nj adjustable co- atmospherically deformed, note that the Strehl defini-
efficients in 4Žj. Utilizing Eqs. (7.7a-d), we can evaluate tion of a diffraction-limited optical system collecting a
D,* and Aj*/nj. deformed wave front is set by the mean square phase
Dc*c_0.992r0(A*) 81 5 , (7.11a)
variation over the aperture. 8 For a mean square phase
3 5 2
DL*_3.40rO(A*) / , (7.1lb) variation A*= 1 rad , the Strehl definition is about 30%,
2
i.e., exp[- (A*) ]. For long exposure photography, the
Ds*_3.79ro(A*)15, (7.11c) applicable mean square phase deviation is (Ac) and
DQ*'5.26ro(A*) l5. 3
(7.11d) should be computed from Eq. (7.11a). A 30% Strehl
definition is thus achieved with a lens diameter about
A c*/ c*Ž0.77ro2 (A*) 18 , (7.12a) equal to ro and an angular resolution of the order of
4
A L*/nLL3.03ro2 (A*)6/ 5 , (7.12b) X/ro. It can be shown that increasing the lens diameter
2
beyond ro increases the phase deviation so rapidly that
A s*/ns-2.82ro (A*) 6 /5 , (7.12c) the achievable angular resolution is not improved be-
A Q*/nQc%.3.63r02(A*)6 15. (7.12d) yond X/ro. For a very short exposure, wavefront tilt, as
distinguished from higher-order-type distortion, does
These are all the quantities that were to be computed. not reduce the resolution of the system. Wavefront tilt
In the next section, their physical significance is displaces the image of a point but does not blur the
discussed. image. Uncorrelated displacement of different points in
an image does, however, distort the picture. Con-
VIII. APPLICATION AND DISCUSSION sequently, for a very short-exposure photograph, the
OF RESULTS
8E. L. O'Neill, Introduction to Statistical Optics (Addison-
To get a quick picture of the shape of the deformed Wesley Publishing Company, Inc., Reading, Massachusetts, 1963),
wave front, we consider the coefficients in Eqs. (7.8a-c). pp. 87, 106.
1432 D. L. FRIED Vol. 55

mean square phase deviation used to compute the computed from the phase-structure function. Formulas
Strehl definition should be taken as (AL), not (Ac). have been provided for these computations for the first
From Eq. (7.11b) it is seen that 30% Strehl definition four types of deformation. It has been shown that
for very short exposures is achieved when the lens according to current theory of optical propagation in the
diameter is about 3. 40ro. The corresponding angular atmosphere a large part of the wavefront deformation
resolution is of the order of X/3.4ro, which is 3.4 times can be understood as wavefront tilting. Sample tech-
as much as the long-exposure resolution. (In a more niques for interpreting the deformation results in terms
precise calculation considering the effect of wavefront of performance of certain types of optical systems have
distortion on the integrated modulation-transfer func- been provided, the implication being that quantitative
tion, we have found that the resolution for high- estimate can be generated of the performance of any
speed exposure is increased by a factor of 2, not 3.4, type of optical system whose behavior is limited by
over the best that can be achieved with a long exposure. wavefront distortion.
The peak resolution is achieved at a lens diameter of
approximately 3.8ro.) ACKNOWLEDGMENTS
For an optical heterodyne receiver, as indicated in
I would like to take this opportunity to thank J. D.
Sec. VI, the performance starts to deteriorate signifi-
Cloud and M. F. Sternberg for assistance in tracking
cantly when the collector diameter becomes larger than down
errors in the evaluation of 5s(r,D) and 5Q(r,D),
ro. This can be understood from Eq. (7.7a) which indi- and for computation
of the 11's.
cates that the rms phase deviation for this size aperture
is about one radian and from consideration of the fact APPENDIX A. EVALUATION OF 5:j(r)
that if the rms difference is this large then there are
many pairs of points in the aperture where the phase In this appendix we outline the procedure used in
difference is 7r. The signal generated by these pairs of evaluating the function 5:(r,D) for j= C, L, S and Q.
points is very small, contributing very little to the total As intermediate steps in the process, we have to evaluate
performance of the detector. the fuictions K#(r,D) for i=0, 1, 2, and 3; J(m,n; r,D)
For a heterodyne receiver which can track at a suffi- for (m,n) = (0,0), (2,0), (0,2), (4,0), (2,2), and (0,4); and
ciently high rate any tilting of the signal wavefront, the L(m,n; r,D) for (m,n)= (4,0), (3,0), (2,0), (1,0), (0,0)
mean square phase deviation as far as the system per- (2,2), (1,2), (0,2), and (0,4). These functions are
formance is concerned, should be computed from Eq. defined as
(7.7b) for (AL) rather than from (Ac). The rras phase
deviation, according to Eq. (7.1 lb) does not become one
Ko(r,D) = fdr'W(J r'+Ir, D)W(I r'- r, D), (Ala)
radian until the collector diameter is 3.40ro. Thus a
tracking system could utilize a collector diameter 3.4
times as large as a nontracking system. This would 2
yield a shot-noise-limited signal-to-noise-ratio improve- Ko(r,D) = fdr'W(1r'+'rI,D)W(|r'--rJ, D)r' , (Alb)
ment over the nontracking system of 20 log3.4~Ž10 dB.
(A more precise calculation would probably yield a
somewhat smaller improvement.) K 2(r,D)= dr'W(|r'+4rj,D)W(Jr'- r D)r'4, (Atc) ,
By use of a rapidly adjustable variable focal-length
lens system, spherical deformation of the collected
wavefront as well as tilt could be tracked out. However, K 3(r,D) = fdr'W( Ir+ 2r |, D)
comparison of Eqs. (7.1lb, c) shows that the resulting J XW(|r'-4rJ, D)(r.r')2 , (Ald)
gain over a system which only tracked tilt would be
rR-!r r+[R2-(R+,',)2]4
trivial.
J(nnl; r,D) = 2] dp dqpn7ql, (A2)
Similar semiquantitative analyses of the effect of
J0 J-R2-(R+A'r)'] !
wavefront deformation on other types of optical sys-
tems can be computed in the same manner as the cases
treated in this section. L(n,,n; rD)=f dzvv(l- V2 )(n+l)I 2 . (A3)
In the discussion of the imaging and the heterodyne SbD

detection systems, it has been assumed that intensity Starting from Eqs. (5.3a, b), we can show that
variations across the aperture were not present, or could
be ignored. This consideration is discussed in AppendixD. 1ic 1

E Fi(x)Fi(x') = -, 2 (A4a)
IX. CONCLUSIONS i=1 7rR

It has been shown that a precise definition of the "7L I


shape of a deformed wavefront can be generated and E Fi(x)Fi(x')= (4r'2 -,2+R 2
), (A4b)
that the statistics of the deformation shape can be i=1 rR4
November 1965 STATISTICS OF WAVEFRONT DISTORTION 1433

ns 1 see that
2
Z Fj(x)Fj(x')= [l2r4+±(6r2- 8R2)r
E
i=7rR6 J(0,0; r,D)=4R2[L(O,O; r,D)], (AlOa)
4
J(2,0; r,D) = 4R [L(2,0; r,D)-2(r/D)L(lO;r,,D)
+ (3r4 -4R 2r2+4R 4)-12(r-r') 2 ], (A4c)
+(r/D)2 L(O,O; rD)], (AlOb)
flQ 1
J(0,2; rD)= !RIEL(0,2; rD)], (A10c)
Z Fj(x)Fj(x')=-(18r'
6
4 +(-3r2- 8R2)r'2
i=l1 7rR
J(4,0; r,D) = 4R6 [L(4,0; r,D)- 4(r/D)L(3,0; r,D)
+ [(9/8)r4-4R r +4R 2 2 4 ]-6(r. r') 2 }). (A4d) +6(r/D) L(2,0; r,D)-4(r/D)3 L(lO; r,D)
2

+(r/D)4L(O,O; r,D)], (A10d)


From this we see that
J(2,2; r,D)= !R'E(2,2; r,D)-2(r/D)L(1,2; r,.D)
gc(r,D)= (1/7rR 2)Ko(rD), (A5a) +(r/D)2L(0,2; r,D)], (AlOe)
4 2
5L(r,D) = (1/7rR )[4Kj(r,D)+ (R - r2)KO(r,D)], (A5b) J(0,4; rD)= '-R L(0,4; r,D)]. (AlOf)
i5s(r,D)= (l/7rR 6)[12K 2(r,D) + (6r2 -8R 2)Kl(r,D) The L(m,n; r,D) functions can be evaluated using
+ (3r-4R2r 2+4R 4)Ko(rD) integration formulas from Dwight's "Table of Inte-
-12K 3(r,D)], (A5c) grals". The results are
5Q(r,D) = (l/7R6) { 18K2 (r,D)+ (- 3r 2 - 8R 2)Kj(r,D) L(4,0; r,D)= '- cos-'(r/D)+[1-(r/D)2]-T- 6 (r/D)5
+ [(9/8)r 4-4R 2r 2+4R4]Ko(r,D) + (1/24)(r/D)3 + 1 (r/D)], (Al la)
-6K 3 (r,D)}. (AMd)
L(3,0; r,D) = [1-(r/D)2 ]E[- (r/D)4
We define a rectangular coordinate system with axes + (1/15)(r/D) 2+ (2/15)],
(Alb)
parallel and perpendicular to the vector r. Let p and q
L(2,0; r,D) = cos-'(r/D)+E(l-(r/D)2]1/2
be the components of r' in this system, with p denoting X[-'-(r/D) 3+8(r/D) ], (Al1c)
the coefficient parallel to r. Considering the range of p
and q in which W(jr'+±rj,D)W(|r'- rf,D) is non- L(l,O; r,D)=E[1-(r1D 'E2-3(rlD)2 2 (Al1d)
zero, i.e., the region of overlap of two circles of diameter
Dwhose centers are at (p ==r,q= O)and (p=-r, q=O) L(O,O; r,D) =2cos-1 (r1D)
we see that +El1-(rlD 'E2-2' (r/D)]J, (A 1le)
Ko(r,D) =J(O,O; r,D)W(r,2D), (A6a) L(2,2; r,D) =16 cos-1(r/D)+ [1- (r1D)2]1 (rlD)
-(7/24)(r/D) 3+ I16(r/D)], (Allf)
Kl(r,D) = [J(2,0; r,D)+J(0,2;r,D)]W(r,2D), (A6b)
K 2(r,D)==[J(4,0; r,D)±+2J(2,2;r,D)
L(1,2; r,D) El
[-(rlD)2]'.
X [R (r/D)4- 2 (r/D)2+ 5 ], (A1 lg)
+J(0,4; r,D)]W(r,2D), (A6c)
1
L(0,2; r,D) = I cos- (r/D)+ [1-(r/D)21-
K 3(r,D) =r2J(2,0; r,D)W(r,2D), (A6d) XE[1(r1D)3-A(r1D)J, (Allh)
where we have made use of the fact that
L(0,4; r,D) = 15 Cos-1(r1D)+E1-(r1D IE2-I(rlD)
r/ 2 =p2+q2, (A7a) +(13/24)(r/D) 3 -} 6 (r/D)]. (Alli)
4
r-4= p 4 +2p 2
q
2
+ q , (A7b) Substituting these equations into Eqs. (AlOa-f) and the
(r- r') =2
=r 2 2
p . (A7c) results into Eqs. (A6a-d), and that set of results into
Eqs. (A5a-d), we obtain Eqs. (5.6a-d).
Examining Eq. (A2), performing the q integration and
making the substitution APPENDIX B. A THEOREM
v= (2p+r)/D, (A8) In this Appendix we prove the theorem that
for p, we get 1 R 2 /2
rD if j=C
- I rdr5j(r,D)r2 = (B1)
J(mn,n; r,D) R
2
O if j AC.
M
4 r
= -R+n+2| dv( v- ) (1-v2)(n+l)12. (A9) We start by noting that
n+l Ir/D D 2 =- yf . ( _n /

Expanding [v-(r/D)]m and considering Eq. (A3) we = X2+y2+ X12+Y12 - 2xx'- 2yy'. (B1)
1434 D. L. FRIED Vol. 55

This can be expressed in terms of the functions F,(x) sequently depends on A because altitude may vary along
and F,(x') giving the path of propagation. CN2 may be written in terms
of the atmospheric refractive-index variance A and the
r 2= [ 2 R8 /12][1[l(x)F.4(X')+Fl(x')F
4 (X)] outer scale of turbulence Lo defined here as the distance
+ 7rR 4I,'1(x)F,(x') - (rR 4 /2) in which turbulence correlation falls to one-half of the
X [F2(x)F2 (x')+F3 (x)F 3 (x')]. (B2) maximum value. This expression is
We note that because of the orthonormality conditions CN2 =AL&-'. (C2)
A reasonable fit to available data for A and L0 is ob-
1
27rR 2 JJ
- dxdx'W(x,D)W(x',D)r'
n
E F,.(x)F,,(x') tained with expressions
'l = 6.7 X 10-'4 exp(-Iz/3200), (C3a)
2
R /2 if j=C L,= 2/0, (C3b)
0 if j C, where h is altitude, and all units of measure are under-
stood to be in mks units. Substituting (C3a, b) into
since the F1F4 term cannot contribute, the FF, term Eq. (C2), and the result into Eq. (Cl), we obtain a
always contributes R 2/2 to the result, and the F2F2 and closed-form result for -A. This result for the most general
F3,F terms each contribute - R2 /4 whenever nj is three type of propagation path (ignoring curvature of the
or greater; i.e., whenever j#C. In this latter case, the earth) is expressible in terms of the incomplete gamma
net contribution is zero. However, transforming to r and function of order two-thirds.
r' coordinates and considering the definition of 5F#(r,D)
as given in Eq. (5.5), we see that APPENDIX D. SCINTILLATION EFFECTS

1 f f This entire article thus far has considered wavefront


R-| J dxdx'W(x,D) IV(x',D)r2 E Fi(x)Fi(x') distortion with the implicit assumption that intensity is
uniform over the area of interest. For some problems
this is indeed the case and the results derived and con-
= - drr25#(r,D) clusions drawn need no qualifications. In this Appendix,
27R 2 J we sketch out a proof of the fact that the results also
1 rD apply when scintillation is present, but with a slightly
=- 2 I rdrr2 5#(r,D). (B4) different interpretation.
R J w Instead of working with only the phase deviation
O(x), let us also consider the log amplitude l(x), defined
From Eqs. (B3) and (B4) we see that as the natural logarithm of the ratio of the instantaneous
1 rD R2/2 if j=C amplitude at x, to what the amplitude would be if there
-
2 I rdr~j(r,D)r2 = (Bs) were no scintillation. Thus, 0(x)-il(x) is the complex
R J10 if phase which completely specifies the instantaneous
wave fluctuation at x. Now rather than consider the
thus proving the theorem. simple phase-structure function as defined in (1.1), we
Similar relationships, such as that consider what we may call the wave-structure function
for which we also use the symbol D(r).
i jDrdrSj(r,D)=.j (for all j) (B6) 5)(r) == IE0(x)-i1(x)1- E0(x')-i1(x')1 12)
= (Eo W-O(X')]2)+ OW _J(X1)]2)' (DI)
can be proven just as easily. where the vertical bars denote taking an absolute value.
We now go back through the paper asking that the
APPENDIX C. EVALUATION OF XA infinite series D(x) and the finite polynomials 1j(x),
define in (3.1) and (3.5), provide the best fit in the sense
It may be shown 2 that -A can be written as
of Eqs. (3.2) and (3.3), not simply to O(x) but to
k(x)-il(x). It is to be understood in (3.2) and appro-
A= 2.91l() ath of dACN2 priate places thereafter that absolute values are to be
propagation taken in squaring. Of course, the coefficients a,; are to
be allowed to take on complex values and Eq. (3.3)
where X is the optical wavelength, A is a variable must be understood as two equations for each value of
defining length along the path of propagation. The gu, one for the real and one for the imaginary part of
quantity CN is called the atmospheric refractive-index a,. The coefficients ac 2, aL2 , as 2 and aQ2, Eqs. (3.4a-d),
structure constant. It is a function of altitude and con- are to be understood as the sum of magnitudes.
November 1965 STATISTICS OF WAVEFRONT DISTORTION 1435

Equation (4.1) now takes the form Offhand, while pleased with the increased precision
of the structure function for propagation in the atmos-
phere, obtained by considering scintillation and working
fA=- 2 2
dxW(x,D)[0 (x)+1 (x)] with the complex phase we have to ask if we can now
assign any meaning to the complex polynomial used to
match the complex phase. For a problem in which there
- JdxW(x,D) is no scintillation, or when the propagation path is
-7rR2-=I
short, the results derived in the main body of the paper
X {f[(x) - il(x)]d,,+ [E(x)+il(x)]a,,}F,,(x) are correct with 5O(r) interpreted as the phase-structure
function and we can assign exact physical meaning to
1 /0 the quantities we have been evaluating. For long-path
+- E a,a^ dxW(x,D)F,(x)F,(x), (D2) propagation, where scintillation may be significant, we
7rR 2v~
make no more than a factor-of-two error by still inter-
where the overbar denotes a complex conjugate. Equa- preting aD(r) as the phase-structure function, and we
tion (4.2) takes a corresponding form. Now, applying retain some of the physical insight into the quantities
the modified form of (3.2) to determine the real and evaluated.
imaginary parts of a,, separately, we find that Eq. To the extent that only phase fluctuation is of im-
(4.3) takes the form portance in a process, we overestimate the extent of the
distortion by a factor of two. In those processes in which
f
a,,= dxW(x,D)[4(x)-il (x)]F,,(x),
intensity fluctuations also are significant we may not
be overestimating the effect of distortion. In the case
of the image-forming process of a lens, for instance,
as might have been expected. All the rest of the equa- intensity variations across the lens result in additional
tions of Sec. 4 follow immediately with k(x) replaced diffraction which limits resolution. If the intensity
by [O(x)-il(x)], and squares understood to involve variations are random, the log amplitude appears in
absolute values. image statistics in a form analogous to the phase so that
Likewise, all the equations of Sec. 5 follow with the the wave-structure function is the quantity which
phase 0 (x) replaced by the complex phase [4(x) - il(x)], appears. Thus, use of Dc*, as computed from (6.1) and
and with D(r), previously the phase-structure function (Cl), to estimate photographic resolution should be
as defined in (1.1), now understood to be the wave- more precise than the derivation, based on phase
structure function as defined in (Dl). Thus Eqs. (5.7) fluctuation alone, would justify. On the other hand, the
and (5.8a-c) are seen to be correct as stated when scintilla- use of DL*, when intensity fluctuations are significant,
tion is present except that 5)(r) has a different interpreta- to estimate high-speed resolution may be more opti-
tion than given in Eq. (1.1). mistic than warranted since the high-speed exposure,
The balance of the paper is better with this inter- freezing the distortion, can compensate only for wave-
pretation of 5(r) than it was initially. Equation (6.1) front tilt, but not for a linear variation of log amplitude.
and (Cl) are precise when D(r) denotes the wave- When intensity fluctuations are not significant, i.e.,
2
structure function. If D(r) is simply the phase-structure when ([l(x)-l(x')]2)«([k(x)-4(x')] ), the use of DL*
function then these two equations are accurate only to estimate high-speed photographic resolution should
when the propagation path is short. For a long path, the be quite accurate.
coefficient a contains an r-dependent factor which varies How to use the geometric statistics of wavefront
slowly and monotonically from a value of about one-half distortion, when to ignore, and how to consider the
at r = 0 to unity for large r. This factor was arbitrarily effect of intensity variations requires careful judgment
set equal to unity in (6.1) and (Cl). by the user.

You might also like