**259**±

**7**is reflected off of a silver-coated diffuse surface. Speckle data are gathered for 16 surface slope angles corresponding to speckle contrast between about 0.55 and 1. Taking the measured data as truth, both equations show error mean and standard deviation of less than 3%. Thus, the theory is validated over the range of this experiment.

*μ*m## 1.

## Introduction

In recent years, the high energy laser (HEL) community has largely shifted focus to short to medium range (tactical) mission goals such as counter-UAV (unmanned aerial vehicle) and counter-RAM (rocket, artillery, and mortar). Concurrently, there have been numerous advances in tactically relevant HEL modeling, including advances related to wave optics simulation, nonlinear thermal blooming, combined turbulence and thermal blooming effects, and lower atmosphere conditions.^{1}2.3.4.5.6.^{–}^{7} However, even at shorter tactical ranges, HEL tracking requirements are very demanding.^{8} Yet little published research exists regarding the potential tracking advantages of active (laser) illumination and the resulting detrimental speckle effects at tactical ranges.^{9}^{,}^{10} Potential advantages of active illumination include eliminating dependence upon target emission and reflection of natural light, increasing signal-to-noise ratio (SNR), and increasing resolution over longer wavelength thermal imaging. However, the coherence of laser illumination also introduces detrimental speckle phenomena. Speckle causes the image of the target to break up into bright and dark patches, which degrades the ability of the tracking algorithm to estimate target position, regardless of whether that algorithm is centroid, correlation, or edge/feature based.^{10}^{,}^{11} Thus, speckle is one source of noise in active tracking, and neglecting it can cause large inaccuracies in performance estimates. This means that speckle effects need to be well understood and characterized so that the potential advantages of active tracking can be weighed against the disadvantages of speckle.

This work investigates the statistics of tactical active-tracking speckle involving partially temporally coherent (polychromatic) illumination. The focus is on speckle contrast, $C$, a measure of the strength of speckle effects defined as $C={\sigma}_{I}/\overline{I}$, where ${\sigma}_{I}$ is irradiance standard deviation and $\overline{I}$ is mean irradiance. Speckle contrast is 1 for fully coherent light and 0 for completely incoherent illumination. It is reduced from unity in tactical engagements due to source bandwidth, surface slope and roughness, and the polarization properties of the reflected light. A fast and accurate numerical model for speckle contrast is developed by improving a prior engineering derivation by Huntley.^{12} Improvements involve the treatment of surface roughness and polarization and the inclusion of a previously derived geometrical correction factor. The resulting modified Huntley equation matches the more exact Hu theory^{13} very well over a wide range of tactically relevant conditions. An experiment is also conducted which validates both Hu and modified Huntley theory over a range of speckle contrast of about 0.55 to 1.

## 2.

## Modeling Approach

Investigation of polychromatic speckle phenomena and reduced contrast speckle began in earnest almost as early as the development of the first commercial lasers in the 1960s.^{14} Polarization diversity was investigated by early researchers, often in combination with other effects.^{15} The impact of surface roughness for polychromatic illumination was investigated both experimentally and theoretically in the 1970s.^{16}17.^{–}^{18} Also, during that timeframe, the effects of geometry on polychromatic speckle for nonimaging systems were investigated by Pederson and Parry.^{19}^{,}^{20} McKechnie then extended the theory to imaging systems.^{21} Much later, Hu derived an expression for polychromatic speckle contrast, which allowed for any imaging and illumination geometry except those involving near glancing angles with the surface.^{13} The recent work involved experimentation, extension, and engineering simplification of this theory.^{12}^{,}^{22}23.24.^{–}^{25}

Three key speckle reduction factors relevant to tactical active tracking are investigated in this work, namely polarization diversity, surface roughness, and surface slope. Other factors include spatial and temporal integration. Spatial integration reduces image speckle contrast appreciably if the speckle size is comparable to or less than the pixel pitch in the detector array, as multiple speckles will influence each pixel’s measurement.^{14}^{,}^{26}27.^{–}^{28} Additionally, target motion can cause temporal decorrelation of the speckle. For example, a linearly translating target will cause time variation in the angle between the surface normal and the line of sight (LOS) of the tracker, which, in turn, causes translation of the speckle field across the aperture, along with other lesser effects. When the speckle field across the aperture decorrelates, the image speckle will also decorrelate. If the speckle changes appreciably over the integration time of the detector, speckle contrast is reduced.^{15} These two factors, spatial and temporal integration, are independent of the speckle reduction factors included in the model developed here, so they can easily be added to the model. Such inclusion will not be discussed further.

The first factor considered in this work, polarization diversity, can reduce speckle contrast by up to a factor of $\sqrt{2}$.^{15} Under polarized illumination, the rough surface may have the effect of separating the reflected light into two linear polarization states, each of which can produce an independent speckle field. Allowing for partial reflection degree of polarization, ${P}_{\mathrm{r}}$, the speckle contrast for polarization diversity is^{14}

The second speckle reduction factor, surface roughness, also has a fairly minor impact. It reduces contrast by up to about a factor of two under realistic conditions. Under the assumptions of a Gaussian source spectrum, Gaussian surface height distribution, and no depolarization, speckle contrast for the surface reflection case is given by

## (2)

$$C={\{1+2{\pi}^{2}{\left(\frac{\mathrm{\Delta}\lambda}{{\lambda}_{0}}\right)}^{2}{\left(\frac{{\sigma}_{\mathrm{h}}}{{\lambda}_{0}}\right)}^{2}{[\mathrm{cos}({\theta}_{\mathrm{ill}})+\mathrm{cos}({\theta}_{\mathrm{im}})]}^{2}\}}^{(-1/4)},$$^{14}An extreme, but realistic, example might involve a wavelength of $1\text{\hspace{0.17em}\hspace{0.17em}}\mu \mathrm{m}$, coherence length of 1 mm, a Gaussian spectrum, collinear illumination and imaging rays, and surface height standard deviation of $400\text{\hspace{0.17em}\hspace{0.17em}}\mu \mathrm{m}$.

^{29}For this case, the minimum speckle contrast is 0.52 at zero slope angle, or zero angle between the imaging ray and the surface normal.

While polarization diversity and surface roughness usually reduce speckle contrast by less than a factor of two, surface slope often has a more significant impact. Hu developed a physical optics equation for speckle contrast reduction due to arbitrary illumination and imaging directions with respect to the target surface normal as

## (3)

$${C}^{2}({\theta}_{\mathrm{im}},{\theta}_{\mathrm{ill}})={[1/{\int}_{-\infty}^{\infty}S(\lambda )\mathrm{d}\lambda ]}^{2}{\iint}_{-\infty}^{\infty}\mathrm{d}{\lambda}_{1}\text{\hspace{0.17em}}\mathrm{d}{\lambda}_{2}S({\lambda}_{1})S({\lambda}_{2})\mathrm{exp}\{-4{\pi}^{2}{\sigma}^{2}{[\mathrm{cos}({\theta}_{\mathrm{im}})+\mathrm{cos}({\theta}_{\mathrm{ill}})]}^{2}{(1/{\lambda}_{2}-1/{\lambda}_{1})}^{2}\}\times \frac{{\left|{\iint}_{-\infty}^{\infty}H\right\{{\lambda}_{1}R[x-\frac{\mathrm{sin}({\theta}_{\mathrm{ill}})}{{\lambda}_{1}\text{\hspace{0.17em}}\mathrm{cos}({\theta}_{\mathrm{im}})}],{\lambda}_{1}Ry\left\}{H}^{*}\right\{{\lambda}_{2}R[x-\frac{\mathrm{sin}({\theta}_{\mathrm{ill}})}{{\lambda}_{2}\text{\hspace{0.17em}}\mathrm{cos}({\theta}_{\mathrm{im}})}-\frac{\mathrm{tan}({\theta}_{\mathrm{im}})}{{\lambda}_{1}}+\frac{\mathrm{tan}({\theta}_{\mathrm{im}})}{{\lambda}_{2}}],{\lambda}_{2}Ry\left\}\mathrm{d}x\text{\hspace{0.17em}}\mathrm{d}y\right|}^{2}}{{\iint}_{-\infty}^{\infty}{|H({\lambda}_{1}Rx,{\lambda}_{1}Ry)|}^{2}\mathrm{d}x\text{\hspace{0.17em}}\mathrm{d}y{\iint}_{-\infty}^{\infty}{|H({\lambda}_{2}Rx,{\lambda}_{2}Ry)|}^{2}\mathrm{d}x\text{\hspace{0.17em}}\mathrm{d}y},$$^{13}Here, Gaussian distributed surface height is assumed. This equation uses few approximations, but it is not directly suitable for most real-time applications or numerical simulations of active imaging, as it requires considerable computation time due to its integral form. However, it will be taken as truth in some comparisons later on.

To address needs for a faster model of sloped surface speckle contrast, Huntley developed a simplified model for engineering applications.^{12} Huntley’s equation defines the number of independent coherence regions seen by each pixel based on geometrical approximations as

Further, Huntley’s values for $\alpha $ and $\beta $ were based on a curve fit of his equation to Hu’s over a wide range of possible conditions.^{12}

Huntley’s equation is exceptionally computationally efficient, but some additional accuracy can be gained at the expense of some of that efficiency. First, speckle reduction due to roughness can be handled separately with more accuracy using Eq. (2), setting $h$ to 0 in Eq. (4). With $h=0$, the Huntley equation immediately reduces to

where a factor of 1.22 has been lumped into the fitting parameter $\alpha $ and the leading unity term has been removed, as it is just an approximation to the geometrical corrections which will be discussed next.Second, previously derived geometrical correction factors can be applied.^{15} While these factors were derived for the one-dimensional (1-D) case of integration of speckle irradiance over time, they are still relevant to the Huntley equation. In Huntley’s geometrical approximation, the $N$ coherence areas are established as bands across the sloped target surface due to the coherence length of the illuminator separating each pixel’s viewing area on target into coherence regions based on apparent surface depth (Fig. 1). When slope causes apparent depth to change by one coherence length, a new coherence region forms. Because the coherence areas are bands with long side perpendicular to the slope gradient, the problem is still 1-D, although the weight of each band is a function of the circularly approximated diffraction-limited pixel field of view. Thus, when the correction factors are applied to Huntley’s $N$, they do not correct it perfectly. The form of the correction factor depends on the source distribution. For a Gaussian source spectrum, it is^{15}

## (7)

$${N}_{\mathrm{eff}}={\left\{\frac{1}{N}\mathrm{erf}\right(\sqrt{\pi}N)-\frac{1}{\pi}{N}^{-2}[1-\mathrm{exp}(-\pi {N}^{2})]\}}^{-1}.$$For a Lorentzian distribution, it is^{15}

These correction factors relate the speckle degrees of freedom from the geometrical model, $N$, to a corrected number of degrees of freedom, ${N}_{\mathrm{eff}}$, which produces more accurate contrast values.

The three speckle reduction effects are then combined. This goal is accomplished through the straightforward multiplication of the numbers of degrees of freedom due to the independent reduction effects of polarization diversity and combined slope and roughness,^{14} while the dependent effects of surface slope and roughness are combined in a root of the sum of the squares (RSS) fashion according to

## (9)

$${N}_{\text{total}}={N}_{\mathrm{P}}[1+\sqrt{{({N}_{\mathrm{eff}}-1)}^{2}+{({N}_{\mathrm{R}}-1)}^{2}}],$$The fitting coefficient, $\alpha $, for the modified Huntley equation was found through least squares percent error fitting to Hu’s theory over a wide range of tactically relevant conditions. Those conditions are shown in Table 1. They include aperture diameters of 0.3 to 1.0 m, ranges of 500 to 20,000 m, wavelengths of 1.0 or $1.5\text{\hspace{0.17em}\hspace{0.17em}}\mu \mathrm{m}$, coherence lengths of 0.1 to 2.0 mm, and both Gaussian and Lorentzian source spectrums for a total of 24 fitting cases. For each case, the modified Huntley equation was fitted to Hu’s theory for slopes of 0 deg to $60+$ deg as shown in Fig. 2. Slope angle is the angle between the surface normal and the illumination and observation rays, which are assumed to be collinear ($\psi ={\theta}_{\mathrm{ill}}=-{\theta}_{\mathrm{im}}$). This assumption is typically well justified for tactical tracking when the illumination and imaging aperture separation is much less than the range to target. The overall least squares error coefficient of $\alpha =1.749$ yielded 0.37% mean error with 2.02% standard deviation over the range of conditions tested.

## Table 1

Range of conditions for modified Huntley equation fitting.

Aperture Dia. (m) | 0.3 | 0.3 | 0.3 | 0.6 | 0.6 | 1.0 | 1.0 | 1.0 | 1.0 | 1.0 | 1.0 | 2.72 mm |

Object Dist. (m) | 500 | 500 | 5k | 10k | 20k | 500 | 500 | 500 | 500 | 5k | 10k | 2.934 |

Wavelength ($\mu \mathrm{m}$) | 1.000 | 1.500 | 1.000 | 1.000 | 1.000 | 1.000 | 1.500 | 1.000 | 1.500 | 1.500 | 1.500 | 0.671 |

Coh. Lgth. (mm) | 0.100 | 1.000 | 1.000 | 1.000 | 2.000 | 0.100 | 0.100 | 1.000 | 1.000 | 0.500 | 1.000 | 0.259 |

Figure 2 shows a comparison of the Hu, Huntley, and modified Huntley equations under conditions for which slope dominates speckle contrast reduction. Here, the surface height standard deviation is $5\text{\hspace{0.17em}\hspace{0.17em}}\mu \mathrm{m}$, a reasonable value, but one for which speckle reduction due to surface roughness is insignificant. All three models make the assumption that surface features are delta correlated spatially in the plane transverse to the surface normal. Combined with surface height standard deviation much greater than the wavelength, this surface model can be considered optically rough.^{14} Hu and modified Huntley results agree well, while the original Huntley results show noticeable disagreement for the larger slope angles. The original Huntley equation was fit to Hu’s theory to minimize error over a wide range of conditions, but the fit compromised large slope accuracy with tactical conditions in exchange for excellent agreement with a range of other conditions. In the modification, the inclusion of the geometrical correction factor allows small error over all slope angles.

In the previous analysis, surface roughness was not a significant factor, as is often the case for tactical active illumination. However, because Hu’s theory models both slope and roughness, it is possible to compare it against both original and modified Huntley theory for conditions in which roughness is significant. Such a comparison is shown in Fig. 3. Here, both coherence length and surface height standard deviation are 1.0 mm. The illumination and imaging angles are again assumed to be collinear. The modified Huntley equation agrees well with Hu’s theory, while the original Huntley equation shows moderate error. Since both slope and roughness are significant, the geometrical correction factor still explains some of the improvement of the modified Huntley equation over the original, while the more exact treatment of surface roughness explains the rest. For 18 cases involving 30 cm aperture diameter, 5 km range, slopes of 0 deg to 60 deg, coherence lengths of 1 to 16 mm, and surface height standard deviations of 0.25 to 4 mm, the mean difference between modified Huntley and Hu was 1.98% with standard deviation of 3.33% for a Gaussian spectrum. For a Lorentzian spectrum, the mean and standard deviation of the differences were 1.36% and 2.87%, respectively. Thus, the modified Huntley equation developed here is highly accurate both when slope dominates the speckle contrast reduction (Fig. 2) and when slope and roughness both contribute to the reduction (Fig. 3).

## 3.

## Experimental Validation

Hu’s theory and the modified Huntley theory developed here are also validated experimentally using the layout shown in Fig. 4. A short coherence length diode laser (Toshiba T0LD9200) provides light with coherence length of $259\pm 7\text{\hspace{0.17em}\hspace{0.17em}}\mu \mathrm{m}$ at 56.0 mA current and 17.00°C ($3.5\pm 0.5\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{mW}$ output power). An aperture and prism pair reduce the elliptical beam to a 1-mm diameter circular beam, which then propagates almost 4 m and through a diverging lens to ensure coherence diameter $>3$ times the pixel field of view. The light scatters off of a silvered diffuser on a rotation mount for control of surface slope. Here, slope angle is the angle between the surface normal and the bisector of the illumination and imaging rays. The rotation mount both controls the slope and allows rotation of the diffuser about its surface normal, thus decorrelating the speckle between measurements and providing multiple independent speckle realizations for each slope angle. The diffuser’s surface height standard deviation was measured at $3.27\text{\hspace{0.17em}\hspace{0.17em}}\mu \mathrm{m}$ using a profilometer, too small to impact speckle contrast. Also, the correlation length was measured at only $122\text{\hspace{0.17em}\hspace{0.17em}}\mu \mathrm{m}$. These two measurements justify the surface as optically rough. The scattered light is passed through a linear polarizer to remove any effects of polarization diversity, for which Hu’s theory does not allow. It is then imaged onto a focal plane for speckle measurement. The detector is a Marlin F131B complementary metal oxide semiconductor sensor with the cover glass removed to eliminate any interference effects that could otherwise arise due to reflections in the cover glass or other coatings. Alternately, the beam can be diverted to a Michelson interferometer, which allows coherence length measurement immediately prior to each speckle measurement.

The experimental dimensions represent a scaled tactical engagement. Fresnel number, ${N}_{\mathrm{F}}={r}^{2}/\lambda R$, where $r$ is the aperture radius, is a common measure of diffractive conditions. For the laboratory setup, ${N}_{\mathrm{F}}=0.94$. For tactical tracking, reasonable dimensions are 20 km range, 30 cm aperture diameter, and $1.0\text{\hspace{0.17em}\hspace{0.17em}}\mu \mathrm{m}$ wavelength, which result in ${N}_{\mathrm{F}}=1.13$, similar to that of the lab setup. Matching speckle statistics is also desirable. For the previous conditions, coherence length of 2.0 mm, and target surface slope of 30 deg, speckle contrast is reduced to about 0.13. The experiment source power was insufficient to reach such low speckle contrast levels but did allow model validation down to about 0.55 contrast.

Figure 5 shows the measured speckle contrast compared to both Hu and Huntley theory for slope angles of 0 deg to 15 deg. The measurements have bars, which represent both vertical and horizontal error bounds. The vertical error bounds are the standard error of the mean, computed from many measurements at each slope angle and sometimes referred to as 68% confidence intervals. The width of each bar represents the uncertainty in measured slope angle. The Hu equation results are also plotted as a high and low bound representing the uncertainty in image distance and imaging aperture diameter. The figure shows excellent agreement between the measurements, Hu’s theory, and the modified Huntley equation. To quantify this observation, taking Hu’s theory as truth, the measurements have $-1.1\%$ mean difference with 2.9% standard deviation, while the modified Huntley equation has 1.4% mean difference with 1.0% standard deviation. The original Huntley equation also agrees well over this range, but begins to break away high around 0.55 contrast, a trend which was seen to continue in Fig. 2 due to the lack of the geometrical correction factor in this equation. From these results, both the Hu theory and the modified Huntley equation are valid models for sloped surface speckle contrast, at least over the contrast range of 0.55 to 1 measured here. Unfortunately, due to lack of sufficient light, the experiment could not be extended to larger angles to more clearly separate the original Huntley equation from the others. Rather, the original Huntley equation is also found to be valid for contrast greater than about 0.7.

## 4.

## Conclusion

In this work, the statistical properties of speckle relevant to tactical active illumination were investigated with a focus on rapid numerical modeling. This focus addresses a current need in the HEL community to simulate such systems for system design, target recognition, aimpoint identification, and tracking algorithm development. The theory developed and validated here will allow understanding of both the detrimental effects of speckle as well as the possible utility of speckle for providing three-dimensional (3-D) information for identification and tracking. The later application involves the fact that speckle contrast often varies across a target due to changes in surface slope. The fast modified Huntley theory will find use in large trade studies and real-time applications for which the integral form of the Hu equation is too slow. Even so, both models were experimentally validated.

The Huntley equation for polychromatic speckle contrast, which uses geometrical approximations to reduce computation time, was improved through the incorporation of a geometrical correction factor and by independent treatment of surface roughness and polarization diversity. This modified Huntley equation was found to match the more exact Hu theory very well over a wide range of tactically relevant conditions. For 24 cases in which surface slope was the dominant speckle reduction factor, the modified Huntley equation showed mean error of 0.37% with 2.02% standard deviation. For 18 cases in which both surface roughness and slope were significant, modified Huntley showed 1.67% error with 3.12% standard deviation. Further, both the modified Huntley equation and Hu’s theory were experimentally validated over a contrast range of about 0.55 to 1. For this experiment, the measurements differed from Hu by $-1.1\%$ mean and 2.9% standard deviation, while the modified Huntley disagreed by only 1.4% mean with 1.0% standard deviation. Thus, both the new modified Huntley equation and Hu’s equation are validated over the speckle contrast range of 0.55 to 1.

Areas for future research include extension of the validation range, exploitation of the 3-D information, and investigation of the impact of optical turbulence. First, the SNR issues which prevented measurements below about 0.55 speckle contrast could be overcome to extend the model validation range either by implementing further noise suppression or using a source capable of maintaining very short coherence length at higher output power. Further, because speckle contrast varies with surface slope, it carries 3-D target information. Exploitation of that information could prove useful for target identification, aimpoint identification, or aimpoint tracking. Finally, because tactical illumination involves light propagation through the turbulent atmosphere, the effects of optical turbulence on polychromatic illumination speckle contrast should be investigated.

## Acknowledgments

The authors recognize the critical support of the High Energy Laser Joint Technology Office in Albuquerque, New Mexico, which sponsored this work as the first author’s master’s thesis research. This research was also supported in part by an appointment to the Postgraduate Research Participation Program at the Air Force Institute of Technology (AFIT) administered by the Oak Ridge Institute for Science and Education through an interagency agreement between the U.S. Department of Energy and AFIT. Gratitude is also extended to two unnamed reviewers whose comments and suggestions greatly improved the paper. The views expressed in this paper are those of the authors and do not necessarily reflect the official policy of the U.S. Air Force, the Department of Defense, or U.S. Government.

## References

## Biography

**Noah R. Van Zandt** is a PhD student with the Center for Directed Energy at the Air Force Institute of Technology (AFIT) in Dayton, Ohio. He received his BS degree in electrical engineering from Cedarville University in 2010 and his MS degree in optical science and engineering from AFIT in 2015. His research interests include atmospheric propagation of high energy lasers, laser beam combination, and active and passive target tracking. He is a student member of SPIE.

**Jack E. McCrae** received his PhD in physics from the AFIT in 1997, an MS in physics (optics) AFIT in 1993, and a BS degree in physics from the MIT in 1984. He is a retired Air Force colonel and currently a research assistant professor in the Engineering Physics Department at AFIT. His research interests include optics, lasers, quantum and non-linear optics, laser radar, atmospheric propagation, and imaging.

**Steven T. Fiorino** received his BS degree in geography and meteorology from Ohio State (1987) and Florida State (1989) Universities, respectively. He additionally holds an MS degree in atmospheric dynamics from Ohio State (1993) and a PhD in physical meteorology from Florida State (2002). He is an AFIT research associate professor of atmospheric physics with research interests in microwave remote sensing, weather signal processing algorithms, and atmospheric effects on military weapon systems. He is a member of SPIE.