Scale factor correction for Gaussian beam truncation in second moment beam radius measurements

Lucas R. Hofer; Rocco V. Dragone; Andrew D. MacGregor

doi:10.1117/1.OE.56.4.043110

25 April 2017 Scale factor correction for Gaussian beam truncation in second moment beam radius measurements

Lucas R. Hofer, Rocco V. Dragone, Andrew D. MacGregor

Author Affiliations +

Optical Engineering, Vol. 56, Issue 4, 043110 (April 2017). https://doi.org/10.1117/1.OE.56.4.043110

Abstract

Charged-couple devices (CCD) and complementary metal oxide semiconductor (CMOS) image sensors, in conjunction with the second moment radius analysis method, are effective tools for determining the radius of a laser beam. However, the second moment method heavily weights sensor noise, which must be dealt with using a thresholding algorithm and a software aperture. While these noise reduction methods lower the random error due to noise, they simultaneously generate systematic error by truncating the Gaussian beam’s edges. A scale factor that is invariant to beam ellipticity and corrects for the truncation of the Gaussian beam due to thresholding and the software aperture has been derived. In particular, simulations showed an order of magnitude reduction in measured beam radius error when using the scale factor—irrespective of beam ellipticity—and further testing with real beam data demonstrated that radii corrected by the scale factor are independent of the noise reduction parameters. Thus, through use of the scale factor, the accuracy of beam radius measurements made with a CCD or CMOS sensor and the second moment are significantly improved.

1. Introduction

Charged-couple devices (CCD) and complementary metal oxide semiconductor (CMOS) image sensors have become standard beam radius measurement tools due to their ability to return a two-dimensional (2-D) beam intensity profile. The first-order moment (arithmetic centroid) has been used with CCD sensors to track the position of lasers,¹ Shack–Hartman wavefronts,² and stars.³ From values such as second moment beam diameter, beam ellipticity, and beam orientation, more complex propagation related beam quality metrics such as $M^{2}$ may be determined.⁴^–⁶ We note that for a pure ${TEM}_{00}$ Gaussian beam, the second moment radius equals the $1 / e^{2}$ (13.53%) radius determined by percent-of-peak intensity. Although previous literature has dealt primarily with the CCD sensor, CMOS sensors are quickly replacing CCD sensors as they are more economical to produce and require less power.⁷

The second moment beam radius measurement technique, in conjunction with the CCD or CMOS sensor, improves on simpler beam profiling devices, which return only a one-dimensional (1-D) beam intensity profile through use of a scanning slit or knife edge. The beam radius on the scanning beam profilers is calculated by finding the radius at which the beam’s 1-D intensity profile intersects a predetermined percentage (e.g., 13.53%) of the beam’s peak intensity.⁸ Although a similar method can be used with CCD or CMOS sensors, the second moment utilizes the full 2-D beam profile, providing a more accurate beam measurement technique. The radius of the monochromatic beam $w_{σ} (z)$ is given by the second moment as⁵

Eq. (1)

w_{σ} (z) = \sqrt{2} σ_{r} (z) .

σ_{r} {(z)}^{2}

is defined as

Eq. (2)

σ_{r}^{2} (z) = \frac{\iint r^{2} I (r, z) r d r d θ}{\iint I (r, z) r d r d θ},

where

r

is the radius extending from the center of the beam,

z

is the beam’s axis of propagation, and

I (r, z)

is the beam’s time-averaged intensity distribution. This paper addresses only monochromatic beams, thus forgoing the use of the second moment polychromatic beam radius definition.⁹

In the past, image sensors were avoided for beam measurement since high noise levels from background light photon noise, readout noise, background level, and sampling error¹⁰ were problematic when coupled with the second moment measurement method, which is highly susceptible to error from noise. However, the last two decades have seen improvements to image sensor hardware and the development of noise reduction techniques, which enable image sensors to provide accurate beam radius measurements. The noise reduction techniques can be understood by considering the second moment integrals [see Eq. (2)], which heavily weight noise at high values of $r$ .¹¹ If—for a theoretical sensor—the limits of $σ_{r}^{2} (z)$ extended toward infinity, then the second moment integral would be completely dominated by noisy pixels.¹² Even with finite dimensions, noise at high $r$ values and outside the beam area can adversely affect the accuracy of small beam measurements. Two techniques are generally used to reduce noise. In the first method, pixels below a certain threshold value are removed from the second moment calculation (thresholding) in an attempt to retain only valid beam data.¹³ However, noisy pixels falling above the threshold value can remain. Therefore, the second technique limits the region about the centroid over which the summation occurs (software aperture) to further reduce the second moment error due to noise.¹⁴

The thresholding algorithm is applied to the image sensor data before using a software aperture or calculating the second moment.¹³ Thresholding entails first sampling an unilluminated portion of the sensor (e.g., the four corners) to determine the mean $μ$ and standard deviation $σ_{n}$ of the background noise. The mean and standard deviation are then used to calculate a threshold value against which each pixel’s intensity value is compared. Pixels falling below the threshold value are considered unilluminated, whereas pixels above the threshold have the mean subtracted from their intensity value and are then used in the second moment calculation. A pixel $I_{p} (x, y)$ is considered illuminated when

Eq. (3)

I_{p} (x, y) > μ + n σ_{n},

where

2 < n < 4

according to the ISO 11146 standard.¹³ If

n

is set too low, noisy pixels are considered illuminated, but if

n

is set too high then valid beam data are lost (see Fig. 1). Therefore, there exists an optimum value of

n

for which noise is excluded from the beam measurement but at which the beam truncation is minimal.¹⁰ Regardless of the value

n

, some portion of a noisy beam is truncated by the thresholding procedure.

Fig. 1

(a) 1-D Gaussian intensity distribution (solid line). Pixels below the threshold level [dotted line, Eq. (3)] are considered unilluminated. The corresponding threshold radius [dash-dot line, Eq. (8)] denotes the truncation of the Gaussian beam by thresholding, so the beam intensity data to the right of the threshold radius is not used in subsequent calculations—including the calculation of the software aperture. The software aperture radius [dashed line, Eq. (12)] shows the final truncation of the Gaussian beam after both thresholding and the software aperture are applied. Only beam intensity data to the left of the software aperture radius is used in the second moment calculations. The thresholding level, threshold radius and software aperture radius are set to arbitrary values to demonstrate truncation. In actuality, the threshold level moves up-down (denoted by double-headed arrows) depending on the present noise levels, which give a corresponding left-right shift to the threshold radius. Furthermore, the software aperture radius moves left-right depending on the software aperture conditions set in the beam analysis software.

After thresholding, the software aperture can be applied to exclude extraneous noise from the second moment calculation yet retain valid beam data. The software aperture is centered on the beam and extends a finite distance from the centroid [see Fig. 2(b)]. While the software aperture removes noise, it also truncates the tails of the Gaussian beam (see Fig. 1), and a balance must be struck between eliminating extraneous noise and removing valid data.¹⁵ Theoretical work and simulations have shown that the optimal software aperture size varies with both the value $n$ used in thresholding¹⁴ and the radius of the beam.¹¹ Although $n$ is predetermined, the beam radius is generally not known before measurement, and an iterative algorithm must determine the optimal software aperture size for the measured beam.⁵

Fig. 2

(a) A ${TEM}_{00}$ elliptical Gaussian intensity distribution in region $S$ surrounded by the software aperture (solid line). (b) The elliptical Gaussian intensity distribution from (a) transformed into a radially symmetric Gaussian distribution in region $S^{'}$ with the inclusion region (solid line) transformed as well. A radially symmetric intensity distribution and software aperture are needed in the scale factor derivation to return analytic solutions from the second moment integrals.

An elliptically invariant scale factor that corrects for the truncation of a Gaussian beam by both thresholding and the software aperture is derived in this paper. Since a scale factor must be based on the expected beam,¹⁵ we use the ubiquitous ${TEM}_{00}$ Gaussian, which has a time-averaged intensity profile $I (r, z)$ represented by

Eq. (4)

I (r, z) = I_{0} e^{- γ r^{2}},

where

I_{0}

is the peak intensity of the beam,

γ

is defined as

Eq. (5)

γ = \frac{2}{w^{2} (z)},

and

w (z)

is the radius of the beam when the intensity is

1 / e^{2}

.

2. Scale Factor Derivation

The scale factor derived herein corrects for the truncation of the Gaussian beam by both thresholding and the software aperture. To derive the scale factor, three main steps are needed. First, the radius (threshold radius) at which the Gaussian beam is truncated by thresholding alone must be determined. This value is then used to calculate the software aperture radius—where the beam is truncated by both noise reduction techniques. Finally, the software aperture radius is applied to the second moment integrals, which yields an analytic solution and allows the scale factor to be determined through algebraic manipulation.

2.1.

Thresholding

Thresholding is performed on the sensor data before the software aperture is applied or the second moment is calculated. For ease of calculation, the standard deviation of the noise is rewritten as the peak beam intensity multiplied by a scalar $k$

Eq. (6)

σ_{n} = k I_{0} .

After the threshold comparison has been applied, the mean is subtracted from the illuminated pixels. The radius

r_{t}

at which the Gaussian intensity distribution is truncated by thresholding can then be found by setting

I (r_{t}, z)

equal to the threshold value

n k I_{0}

Eq. (7)

I_{0} e^{- γ r_{t}^{2}} = n k I_{0} .

Solving Eq. (7) for

r_{t}

yields

Eq. (8)

r_{t} = \sqrt{\frac{1}{γ} \ln (\frac{1}{n k})} .

2.2.

Software Aperture

The software aperture is applied after thresholding to include valid beam data and exclude extraneous noise from the second moment calculations. Although previous work has focused on a square or rectangular software aperture,¹²^,¹³ a radially symmetric software aperture (for a radially symmetric beam) has the advantage of excluding extraneous noise found in the corners of a rectangle or square. We define the software aperture as a circular region centered on the beam, which includes a portion $β$ of the sensors total power $P_{T}$ after thresholding. $β$ is given as

Eq. (9)

β = \frac{P_{S}}{P_{T}},

where

P_{S}

is the power within the software aperture. To calculate

P_{S}

and

P_{T}

, the equation for power is used

Eq. (10)

P = \int_{0}^{2 π} \int_{0}^{R} I_{0} e^{- γ r^{2}} r d r d θ = \frac{π I_{0}}{γ} (1 - e^{- γ R^{2}}) .

Nominally, the radius limit

R

would extend to infinity for

P_{T}

; however, thresholding truncates the beam before the software aperture’s application, so the beam’s remaining power is contained within the radius

r_{t}

. Thus, for the total power,

R = r_{t}

. For the power in the software aperture,

R = r_{s}

, where

r_{s}

is the radius of the software aperture. Now Eq. (9) can be rewritten

Eq. (11)

β = \frac{P_{S}}{P_{T}} = \frac{\frac{π I_{0}}{γ} (1 - e^{- γ r_{s}^{2}})}{\frac{π I_{0}}{γ} (1 - e^{- γ r_{t}^{2}})} = \frac{1 - e^{- γ r_{s}^{2}}}{1 - n k},

with Eq. (8) replacing

r_{t}

. Solving Eq. (11) for

r_{s}

yields

Eq. (12)

r_{s} = \sqrt{\frac{1}{γ} \ln (\frac{1}{1 - ν})},

where

ν

is defined as

Eq. (13)

ν = β (1 - n k) .

Equation (12) gives the radius at which the Gaussian beam is truncated by both thresholding and the software aperture (see Fig. 1).

2.3.

Scale Factor

To find the radius of the beam, Eq. (2) is first rewritten with integration limits defined by the software aperture

Eq. (14)

σ_{r}^{2} (z) = \frac{\int_{0}^{2 π} \int_{0}^{r_{s}} r^{2} I (r, z) r d r d θ}{\int_{0}^{2 π} \int_{0}^{r_{s}} I (r, z) r d r d θ} .

Evaluating the integral in the numerator gives

Eq. (15)

\int_{0}^{2 π} \int_{0}^{r_{s}} r^{3} I_{0} e^{- γ r^{2}} d r d θ = \frac{π I_{0}}{γ} (\frac{1 - e^{- γ r_{s}^{2}} (γ r_{s}^{2} + 1)}{γ}),

whereas the integral in the denominator yields

Eq. (16)

\int_{0}^{2 π} \int_{0}^{r_{s}} I (r) r d r d θ = \frac{π I_{0}}{γ} (1 - e^{- γ r_{s}^{2}}) .

Substituting Eqs. (15) and (16) into Eq. (14) and simplifying gives

Eq. (17)

σ_{r}^{2} (z) = \frac{1}{γ} - \frac{r_{s}^{2}}{e^{γ r_{s}^{2}} - 1} .

Next, we solve Eq. (17) for

σ_{r} (z)

in terms of

ν

by substituting Eq. (12) for

r_{s}

Eq. (18)

σ_{r} (z) = \frac{1}{γ} {[1 + \frac{(1 - ν) \ln (1 - ν)}{ν}]}^{\frac{1}{2}} .

Now, the second moment radius

w_{σ} (z)

is found by substituting Eq. (18) into Eq. (1) and replacing

γ

with Eq. (5)

Eq. (19)

w_{σ} (z) = w (z) {[1 + \frac{(1 - ν) \ln (1 - ν)}{ν}]}^{\frac{1}{2}} .

Solving Eq. (19) for the actual beam radius

w (z)

in terms of

w_{σ} (z)

gives

Eq. (20)

w (z) = w_{σ} (z) ψ (ν),

where

ψ (ν)

is the scale factor defined by

Eq. (21)

ψ (ν) = {[1 + \frac{(1 - ν) \ln (1 - ν)}{ν}]}^{- \frac{1}{2}} .

Thus, the second moment radius—subject to systematic error from truncation—can be multiplied by the scale factor (see Fig. 3) to give the beam radius devoid of truncation error. The variables in

ν

are either predetermined or found using beam analysis software, thereby allowing the scale factor to be calculated for a measured beam.

Fig. 3

(a) Scale factor as a function of $ν$ [see Eq. (21)], which corrects for truncation of the Gaussian beam by both thresholding and the software aperture. Values of $ν < 0.5$ represent a significant uncorrected error that should be addressed by either tuning the noise reduction settings or reducing the actual noise on the sensor. Note that, as expected, the scale factor goes to one as $ν$ goes to one [see Eq. (24)].

Although the scale factor accounts for both thresholding and the software aperture, it simplifies when only one noise reduction technique is used (either thresholding or the software aperture). If only thresholding is used, $β = 1$ and $ν = 1 - n k$ , so the scale factor reduces to

Eq. (22)

ψ (n, k) = {[1 + \frac{n k \ln (n k)}{(1 - n k)}]}^{- \frac{1}{2}} .

Conversely, when only the software aperture is applied,

n = 0

and

ν = β

. Thus, the scale factor becomes

Eq. (23)

ψ (β) = {[1 + \frac{(1 - β) \ln (1 - β)}{β}]}^{- \frac{1}{2}} .

If neither thresholding nor the software aperture is used, the limit as

β

goes to one and

n

goes to zero yields

Eq. (24)

\lim_{(β, n) \to (1,0)} ψ (β, n) = 1 .

Equation (20) can then be rewritten with

ψ (ν)

replaced by Eq. (24) as

Eq. (25)

w (z) = w_{σ} (z),

which is the expected result since the Gaussian beam is truncated by neither thresholding nor the software aperture.

It can be shown (see Appendix A) that the scale factor is invariant to the ellipticity of the beam. However, to find analytical solutions to the second moment integrals, the scale factor derivation must be performed with a radially symmetric Gaussian intensity distribution and the polar second moment equations. Therefore, an elliptical intensity distribution in region $S$ is rewritten as a radially symmetric distribution in the new region $S^{'}$ (see Fig. 2) through use of a linear transform. The scale factor derivation demonstrated above is then performed in region $S^{'}$ with the polar second moment equations. After transforming from region $S^{'}$ back to $S$ , the scale factor has the same form as Eq. (21) and corrects the second moment calculated radii along both the $x$ and $y$ axes.

3. Testing

Both simulations and real beam data were used to test the scale factor. Whereas simulations provided a more in-depth analysis of the scale factor’s effectiveness—by comparing the measured beam radii to the actual beam radii, determined a priori—the analysis of real beam data was needed to demonstrate the efficacy of the scale factor in real-world applications. However, since the radius of a real beam cannot be predetermined, the noise reduction parameters ( $n$ and $β$ ) were varied during measurement to show that radii corrected by the scale factor are invariant to $n$ and $β$ as derived.

3.1.

Simulations

Simulated laser beams were generated to test the scale factor’s effectiveness and invariance to ellipticity by comparing the generated radii $w (z)$ to the measured radii $w_{σ} (z)$ and determining the error in the measured radii. The beams were simulated on a 16-bit sensor of $512 \times 512 pixels$ , with each pixel’s length and width equal to $p_{w} = 5.5 μ m$ . Although most image sensors are 12 bit, error from analog-to-digital quantization can be ignored when using a 16-bit sensor.¹⁴ The beams were generated with the beam’s peak amplitude scaled to $0.9 a_{m}$ , where $a_{m}$ is the total analog-to-digital converter levels for a sensor (i.e., $a_{m} = 2^{16}$ for a 16-bit sensor).

Beam analysis software calculated the second moment radii of the generated beams and allowed thresholding and the software aperture to be turned on or off independently. Thresholding was applied with $n$ (the scalar of the standard deviation of the noise) set to 4 (highest allowable value per the ISO 11146 standard). Generally, $n$ is optimized to sufficiently threshold noise yet minimize beam truncation; however, the scale factor corrects error due to truncation, allowing high values of $n$ —which remove more noise—to be used. The software aperture radius (or major and minor axes for an ellipse) was iterated by the software until

Eq. (26)

0.990 P_{T} \leq P_{S} \leq 0.992 P_{T},

which corresponds to an aperture radius roughly 1.5 times the beam radius. When both thresholding and the software aperture were turned on, the software aperture radius was calculated after thresholding had been applied.

Due to the layout of the image sensor’s pixels as a rectangular grid, the second moment equations in Cartesian coordinate form, rather than polar coordinate form, were used. Since, for the beams generated, the angle between the measurement axes and the beam axes was zero, the second moment equations for the beam radii simplified to¹⁶

Eq. (27)

w_{σ_{x}} (z) = 2 σ_{x} (z),

Eq. (28)

w_{σ_{y}} (z) = 2 σ_{y} (z),

where

σ_{x}^{2} (z)

and

σ_{y}^{2} (z)

for a discretized intensity distribution

I_{p} (x, y)

are given by

Eq. (29)

σ_{x}^{2} (z) = \frac{\sum_{}^{} {(x - \bar{x})}^{2} I_{p} (x, y) Δ x Δ y}{\sum_{}^{} I_{p} (x, y) Δ x Δ y},

Eq. (30)

σ_{y}^{2} (z) = \frac{\sum_{}^{} {(y - \bar{y})}^{2} I_{p} (x, y) Δ x Δ y}{\sum_{}^{} I_{p} (x, y) Δ x Δ y} .

The

\bar{x}

and

\bar{y}

centroids are defined as

Eq. (31)

\bar{x} = \frac{\sum_{}^{} x I_{p} (x, y) Δ x Δ y}{\sum_{}^{} I_{p} (x, y) Δ x Δ y},

Eq. (32)

\bar{y} = \frac{\sum_{}^{} y I_{p} (x, y) Δ x Δ y}{\sum_{}^{} I_{p} (x, y) Δ x Δ y},

with

Δ x

and

Δ y

each equal to

p_{w}

. Equations (27) and (28) were each multiplied by the scale factor (calculated via the beam analysis software for each beam) to give the corrected beam radius. The analysis software’s calculated radii (with and without the scale factor) were compared to the generated beam’s actual radii to determine the percent error. We postulated that the calculated radii corrected by the scale factor would return a smaller percent error than the uncorrected radii and that the scale factor would correct the radii values irrespective of the beam’s ellipticity.

Three sets of radially symmetric beams were generated and analyzed with different noise reduction settings to characterize the effectiveness of the scale factor. Whereas the first set of beams used only thresholding and the second set of beams applied only the software aperture, the third set of beams employed both thresholding and the software aperture together. Gaussian noise with a mean of $μ = 0.0174 a_{m}$ and a standard deviation of $σ_{n} = 0.0036 a_{m}$ was added to the first and third set of beams to simulate noise on the image sensor. The values for $μ$ and $σ_{n}$ were chosen to reflect the noise levels seen when measuring real beam data. Since, the second set of beams was analyzed without thresholding, the scale factor reduced to Eq. (23) and had no dependence on noise; therefore, for these beams, no noise was added.

The error of the corrected radii was over an order of magnitude lower than the error of the uncorrected radii across all three sets of calculated beam radii. The error for the first set of calculated beam radii exhibited exponentially decaying behavior and became asymptotic at a radius of $30 p_{w}$ [see Fig. 4(a)] with an average error of 3.46% and 0.07% for the uncorrected and corrected beams, respectively. The smaller beams had fewer illuminated pixels, and the noise left after the thresholding process dominated the second moment calculation, which led to high errors. For the second set of beams, the error was approximately constant with an average of 2.41% for the uncorrected radii and 0.20% for the corrected radii [see Fig. 4(b)]. The third set of beams showed exponential decay, which reached an asymptote around $30 p_{w}$ with an average uncorrected error of 5.24% and a corrected error of 0.05% [see Fig. 4(c)]. Although the large error for the smaller beams was similar to the first set of beams (thresholding only), the software aperture removed outlying noisy pixels, which significantly reduced the error compared to thresholding alone.

Fig. 4

(a) The radius percent error as a function of the beam radius with only thresholding applied. Both the uncorrected error (triangle) and the corrected error (diamond) showed decaying exponential behavior, but quickly became asymptotic around $30 p_{w}$ with average errors of 3.46% and 0.07%, respectively. (b) The radius percent error as a function of the beam radius using only the software aperture. The uncorrected error (triangle, average error 2.41%) and the corrected error (diamond, average error 0.20%) were approximately constant for the various beam radii. (c) The radius percent error as a function of the beam radius with both thresholding and the software aperture applied. The error decayed exponentially, reaching the asymptote around $30 p_{w}$ , but with an error significantly lower than in (a). Although the uncorrected error (triangle, average error 5.24%) for the asymptotic region was higher than in (a) and (b), the corrected error (diamond, average error 0.05%) was of similar magnitude. For (a)–(c), the error from beam radius measured along the $x$ axis and the error from the beam radius measured along the $y$ axis were averaged to give a single error value. In (a) and (c), the beams were generated and measured 10 times at each point to give the average and standard deviation (error bars). The beams in (b) were only measured once since they contained no noise.

Whereas the first and second tests demonstrated that the scale factor independently corrects for beam truncation due to either thresholding or the software aperture, the third test showed it can correct for them simultaneously. Indeed, thresholding and the software aperture were needed in conjunction since thresholding alone was unable to return accurate beam measurements for the small beams. Although the uncorrected error was significantly higher when both thresholding and the software aperture were used (due to the combination of systematic errors), the corrected error was of similar magnitude as the corrected error when only thresholding or the software aperture were employed.

Finally, a set of elliptical beams was generated with Gaussian noise ( $μ = 0.0174 a_{m}$ and $σ_{n} = 0.0036 a_{m}$ ) to determine whether the scale factor is, as derived, invariant to the ellipticity of the ${TEM}_{00}$ Gaussian beam (see Fig. 5). The generated beams had ellipticities ( $ε$ ) ranging from 0.08 to 1.00, and their second moment radii were calculated after both thresholding and the software aperture had been applied. The error for the corrected radii was approximately constant as the beam ellipticities changed (see Fig. 5), and the scale factor reduced the error from an average of 5.23% (uncorrected) to 0.06% (corrected). Although slight variations were seen in the corrected error, these were due to the random Gaussian noise in the beams rather than changes in ellipticity. The significant reduction in beam radius error coupled with the consistency of the corrected error, with respect to the changes in ellipticity, demonstrated that the scale factor corrects for the truncation of a Gaussian beam irrespective of ellipticity.

Fig. 5

The radius percent error as a function of the beam ellipticity with both thresholding and the software aperture applied. Both the uncorrected error (triangle, average percent error 5.23%) and the corrected error (diamond, average percent error 0.06%) were approximately constant with varying beam ellipticity. The beams were generated and measured ten times at each point to give the average and standard deviation (error bars). Finally, the error from beam radius measured along the $x$ axis and the error from the beam radius measured along the $y$ axis were averaged to give a single error value.

3.2.

Real Beam Data

Since the radius of a real laser beam is determined only through measurement and not known a priori, the measured radius cannot be compared to the exact radius as with the simulated beams. However, noting that the scale factor corrects for the second moment measurement’s dependence on the noise reduction parameters $n$ and $β$ [see Eq. (20)], we postulated that systematically varying the noise parameters would result in a varied uncorrected radius and a constant corrected radius. Furthermore, the invariance of the corrected radius to $n$ and $β$ reflected the physical description of a real laser beam—the radius of which does not change based on software parameters—and allowed the scale factor to be tested on real beam data with a predicted result.

Two tests were implemented using a single mode ${TEM}_{00}$ Gaussian beam with a wavelength of 675 nm and a beam radius of $\sim 145 μ m$ measured on a 12-bit image sensor of $512 \times 512 pixels$ (pixel height and width $5.5 μ m$ ) with the beam analysis software described in the simulation section. Whereas the first test used only thresholding with $n$ varied from 3.5 to 4 (for $n > 3.5$ the effect of the scale factor is unobscured by noise), the second test used both thresholding and the software aperture with $n$ held constant at $n = 4$ and $β$ evenly incremented from 0.985 to 0.995. The measured radii (both uncorrected and corrected) were divided by their initial radius values to show deviations from the initial values (Fig. 6).

Fig. 6

(a) Radius as a function of $n$ divided by the initial radius $r_{0}$ for both the uncorrected (triangles) and corrected radii (diamonds) with only thresholding applied. As $n$ was incremented, the uncorrected beam radius decreased and the corrected beam radius remained approximately constant. (b) Radius as a function of $β$ divided by the initial radius $r_{0}$ for both the uncorrected (triangles) and corrected radii (diamonds). Thresholding and the software aperture were both used, but $n$ was held constant at $n = 4$ . The uncorrected beam radii increased as the value of $β$ increased (larger software aperture), whereas the corrected radii remained approximately constant. For both (a) and (b), the error from beam radius measured along the $x$ axis and the error from the beam radius measured along the $y$ axis were averaged to give a single error value. Additionally, the beam was measured 10 times at each point to give the average and standard deviation (error bars).

In both tests, the uncorrected radii decreased as the truncation of the beam increased (larger $n$ and smaller $β$ ); however, the corrected radii remained approximately constant, demonstrating that the radii corrected by the scale factor are invariant to changes in $n$ and $β$ as predicted. Furthermore, the tests—which were performed on real beam data—showed that the scale factor effectively corrects the second moment radius measurements in real-world application.

4. Conclusion

An image sensor (CCD or CMOS) in conjunction with the second moment radius analysis method can accurately measure the radius of a laser beam; however, the second moment method heavily weights noise and noise reduction techniques such as thresholding and a software aperture must be used before calculating the second moment beam radius. Even as these techniques reduce the random error due to noise, they simultaneously generate systematic error by truncating the tails of the Gaussian beam. A scale factor, invariant to the ellipticity of the beam, was derived to correct the second moment beam radius in response to this systematic error.

Both simulations and real beam data were used to test the scale factor. Simulated laser beams tested the effectiveness of the scale factor and its elliptical invariance by comparing the measured beam radii—corrected and uncorrected—to the generated beams’ actual radii. The measured beam radii were calculated with three different noise reduction settings: only thresholding, only the software aperture, and both thresholding and the software aperture. With all three methods, the scale factor reduced the error of the measured radii by over an order of magnitude; thus, as predicted, the scale factor significantly improved the accuracy of the second moment beam radius measurements. Next, a set of elliptical beams was generated and the radii were measured with both thresholding and the software aperture applied. The corrected error of the beam radii was approximately constant as the ellipticity changed, proving that the scale factor is, as derived, invariant to beam ellipticity. Finally, the scale factor was tested with real beam data by systematically varying the values used in the noise reduction techniques. The observed behavior—the corrected radii remained constant, while the uncorrected radii changed—matched both the derived and physical predictions, thus demonstrating the scale factor’s effectiveness in real-world applications.

Further simulations could be done to test the scale factor when either the noise levels are varied or the size of the simulated sensor is changed (e.g., 256, 1024, or 2048 pixels across). Additional theoretical work could be done to derive a scale factor that is valid for higher-order modes (i.e., Hermite–Gaussian beams and Laguerre–Gaussian beams).

Appendices

Appendix A:

Scale Factor Invariance for Elliptical Gaussian Beams

The scale factor is invariant to the ellipticity of the beam as can be shown by transforming an elliptical beam in region $S$ to a radially symmetric beam in region $S^{'}$ and performing the scale factor derivation there. The time-averaged intensity $I (x, y, z)$ of an elliptical Gaussian beam is given in a Cartesian coordinate system $S$ as

Eq. (33)

I (x, y, z) = I_{0} e^{- 2 (\frac{x^{2}}{w_{x}^{2} (z)} + \frac{y^{2}}{w_{y}^{2} (z)})},

where

I_{0}

is the peak intensity,

w_{x} (z)

is the radius of the beam along the

x

axis, and

w_{y} (z)

is the radius of the beam along the

y

axis. A new region

S^{'}

with coordinates

x^{'}

and

y^{'}

is formed from

S

with a transform

Eq. (34)

| \begin{matrix} \frac{w_{y} (z)}{w_{x} (z)} & 0 \\ 0 & 1 \end{matrix} | | \begin{matrix} x \\ y \end{matrix} | = | \begin{matrix} x^{'} \\ y^{'} \end{matrix} |,

which yields

Eq. (35)

x^{'} = \frac{w_{y} (z)}{w_{x} (z)} x

and

Eq. (36)

y^{'} = y .

Substituting Eqs. (35) and (36) into Eq. (33) gives

Eq. (37)

I (x, y, z) = I (x^{'}, y^{'}, z) = I_{0} e^{\frac{- 2}{w_{y}^{2} (z)} (x^{' 2} + y^{' 2})} .

In polar coordinates, Eq. (37) can be written

Eq. (38)

I (r^{'}, z) = I_{0} e^{- τ r^{' 2}},

where

τ

is defined as

Eq. (39)

τ = \frac{2}{w_{y}^{2} (z)} .

Thus, the radially symmetric Gaussian intensity distribution in Eq. (4) is recovered with

γ

replaced by

τ

(see Fig. 2).

Integrals are transformed when moving from $S$ to $S^{'}$ with the following equation:¹⁷

Eq. (40)

\iint_{S} f (x, y) d x d y = \iint_{S^{'}} f (x^{'}, y^{'}) | \frac{\partial (x, y)}{\partial (x^{'}, y^{'})} | d x^{'} d y^{'},

where

| \frac{\partial (x, y)}{\partial (x^{'}, y^{'})} |

is the Jacobian determinant given by

Eq. (41)

| \frac{\partial (x, y)}{\partial (x^{'}, y^{'})} | = | \begin{matrix} \frac{\partial x}{\partial x^{'}} & \frac{\partial y}{\partial x^{'}} \\ \frac{\partial x}{\partial y^{'}} & \frac{\partial y}{\partial y^{'}} \end{matrix} | = | \begin{matrix} α & 0 \\ 0 & 1 \end{matrix} | = α,

and

α

is defined as

Eq. (42)

α = \frac{w_{x} (z)}{w_{y} (z)} .

Equation (40) can be rewritten by changing to polar coordinates and substituting Eq. (41) for the Jacobian determinant

Eq. (43)

α \iint_{S^{'}} f (r^{'}, θ^{'}) r d r^{'} d θ^{'} .

Therefore, the transformation from region

S

to region

S^{'}

requires the integral to be multiplied by the scalar

α

. However, the integrals in the scale factor derivation are always found as a ratio [Eqs. (2) and (9)], so the factor

α

cancels out when the scale factor derivation is performed in region

S^{'}

with the transformed elliptical Gaussian. Using the results of Eqs. (38) and (43), Eq. (18) becomes

Eq. (44)

σ_{r}^{'} (z) = \frac{1}{\sqrt{τ ψ^{2} (ν)}},

and using Eq. (1)

Eq. (45)

w_{σ}^{'} (z) = \sqrt{2} σ_{r}^{'} = \frac{1}{ψ (ν)} \sqrt{\frac{2}{τ}} .

Substituting Eq. (39) for

τ

and solving for

w_{y} (z)

gives

Eq. (46)

w_{y} (z) = w_{σ^{'}} (z) ψ (ν) .

Using Eqs. (42), (46) becomes

Eq. (47)

w_{x} (z) = α w_{σ^{'}} (z) ψ (ν) .

Transforming

w_{σ^{'}} (z)

from region

S^{'}

back to region

S

along the

x

and

y

axes yields

Eq. (48)

w_{σ^{'}} (z) = \frac{1}{α} w_{σ_{x}} (z),

Eq. (49)

w_{σ^{'}} (z) = w_{σ_{y}} (z),

where

w_{σ_{x}}

and

w_{σ_{y}}

are the measured second moment radii along the

x

and

y

axes, respectively. Substituting Eqs. (48) into (47) and Eqs. (49) into (46) gives

Eq. (50)

w_{x} (z) = w_{σ_{x}} (z) ψ (ν),

Eq. (51)

w_{y} (z) = w_{σ_{y}} (z) ψ (ν) .

Thus, the scale factor is invariant to the ellipticity of a

{TEM}_{00}

Gaussian beam.

Acknowledgments

The authors would like to thank Steve Garvey for his contributions to the numeric scale factor preceding the analytic scale factor and John Cadigan for his help with Python simulations. This work was supported by DataRay Inc.

References

1.

R. P. Loce and R. E. Jodoin, “Sampling theorem for geometric moment determination and its application to a laser beam position detector,” Appl. Opt., 29 (26), 3835 –3843 (1990). http://dx.doi.org/10.1364/AO.29.003835 APOPAI 0003-6935 Google Scholar

2.

M. Xia et al., “Shack-Hartmann wavefront sensor with large dynamic range,” J. Biomed. Opt., 15 (2), 026009 (2010). http://dx.doi.org/10.1117/1.3369810 JBOPFO 1083-3668 Google Scholar

3.

E. Dennison and R. Stanton, “Ultra-precise star tracking using charge coupled devices (CCDs),” in 24th Annual Technical Symp., 54 –63 (1980). Google Scholar

4.

A. Siegman, “How to (maybe) measure laser beam quality,” DPSS Lasers: Applications and Issues, 184 –199 Optical Society of America, Washington, D.C. (1998). Google Scholar

5.

“Lasers and laser-related equipment-test methods for laser beam widths, divergence angles and beam propagation ratios,” (2005). Google Scholar

6.

A. E. Siegman, “New developments in laser resonators,” Proc. SPIE, 1224 2 –14 (1990). http://dx.doi.org/10.1117/12.18425 PSISDG 0277-786X Google Scholar

7.

N. Waltham, “CCD and CMOS sensors,” Observing Photons in Space, 423 –442 Springer, New York (2013). Google Scholar

8.

A. E. Siegman, M. Sasnett and T. Johnston, “Choice of clip levels for beam width measurements using knife-edge techniques,” IEEE J. Quantum Electron., 27 (4), 1098 –1104 (1991). http://dx.doi.org/10.1109/3.83346 IEJQA7 0018-9197 Google Scholar

9.

L. Martí-López and O. Mendoza-Yero, “Effect of chromatic aberration on gaussian beams: non-dispersive laser resonators,” Opt. Laser Technol., 31 (3), 239 –245 (1999). http://doi.org/10.1016/S0030-3992(99)00047-X Google Scholar

10.

X. Ma, C. Rao and H. Zheng, “Error analysis of CCD-based point source centroid computation under the background light,” Opt. Express, 17 (10), 8525 –8541 (2009). http://dx.doi.org/10.1364/OE.17.008525 OPEXFF 1094-4087 Google Scholar

11.

T. S. Ross and W. P. Latham, “Appropriate measures and consistent standard for high energy laser beam quality (postprint),” Journal of Directed Energy, 2 22 –58 (2006). Google Scholar

12.

T. S. Ross, Laser Beam Quality Metrics, SPIE Press, Bellingham, Washington (2013). Google Scholar

13.

“Lasers and laser-related equipment-test methods for laser beam widths, divergence angles and beam propagation ratios,” (2004). Google Scholar

14.

H. Yuanxing and L. Xinyanga, “Error analysis of Gaussian spot width measured with CCD sensor,” Proc. SPIE, 8417 841720 (2012). http://dx.doi.org/10.1117/12.970553 PSISDG 0277-786X Google Scholar

15.

G. Mana, E. Massa and A. Rovera, “Accuracy of laser beam center and width calculations,” Appl. Opt., 40 (9), 1378 –1385 (2001). http://dx.doi.org/10.1364/AO.40.001378 APOPAI 0003-6935 Google Scholar

16.

“Lasers and laser related equipment-Test methods for laser beam parameters-Beam widths, divergence angles and beam propagation ratios,” (1999). Google Scholar

17.

M. L. Boas, Mathematical Methods in the Physical Sciences, Wiley, New York (1966). Google Scholar

Biography

Lucas R. Hofer is a graduate student at the Universität Stuttgart, holding a fellowship from the International Max Planck Research School for Condensed Matter and was formerly a technical scientist at DataRay Inc. He received his BS degree in physics from the University of California, Los Angeles, USA, in 2015. He is the author of one journal paper detailing an innovative Mach–Zehnder interferometer control system.

Rocco V. Dragone is the CTO and VP of engineering at DataRay Inc. He received his BS and MS degrees in electrical engineering from Drexel University in 2008. His current area of interest is laser beam profiling.

Andrew D. MacGregor has a BSc degree in physics and a PhD in instrumentation for IR Astronomy from Imperial College, University of London. After work on IR imaging systems at Marconi Avionics and in fiber optic sensors, he joined RCA New Products in Quebec. Avalanche photodiode teamwork led to record sensitivity for intersatellite digital optical receivers and the first TTL output single photon counting module. Since 1998, he has proudly supported DataRay’s laser beam profiling instrumentation. He is a member of SPIE, OSA, and IEEE LEOS.

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation Download Citation

Lucas R. Hofer, Rocco V. Dragone, and Andrew D. MacGregor "Scale factor correction for Gaussian beam truncation in second moment beam radius measurements," Optical Engineering 56(4), 043110 (25 April 2017). https://doi.org/10.1117/1.OE.56.4.043110

Received: 16 November 2016; Accepted: 10 April 2017; Published: 25 April 2017

Access the abstract

JOURNAL ARTICLE
9 PAGES

DOWNLOAD PAPER SAVE TO MY LIBRARY