The Waveform Model of Laser Altimeter System with Flattened Gaussian Laser
 Author: Ma Yue, Wang Mingwei, Yang Fanlin, Li Song
 Publish: Journal of the Optical Society of Korea Volume 19, Issue4, p363~370, 25 Aug 2015

ABSTRACT
The current waveform model of a laser altimeter is based on a Gaussian laser beam of fundamental mode, while the flattened Gaussian beam has many advantages such as nearly constant energy distribution on the center of the crosssection. Following the theory of the flattened Gaussian beam and the waveform theory of the laser altimeter, some of the primary parameters of the received waveform were derived, and a laser altimetry waveform simulator and waveform processing software were programmed and improved under the circumstance of a flattened Gaussian beam. The result showed that the bias between theoretical and simulated waveforms was less than 3% for every order mode, the waveform width and range error would increase as target slope or order number rose. Under higher order mode, the shapes of the received waveforms were no longer Gaussian, and could be fitted more precisely as a generalized Gaussian function with power bigger than 2. The flattened beam got much better performance for a multisurface target, especially when the small surface is far from the center of the laser footprint. This article provides the waveform theoretical basis for the use of a flattened Gaussian beam in a laser altimeter.

KEYWORD
Remote sensing , Laser altimeter , Flattened Gaussian beam , Waveform model , Generalized Gaussian

I. INTRODUCTION
High resolution and accurate elevation of Earth or other planets can be acquired using spaceborne laser altimeters. The range between the reference center on a satellite and the target surface is calculated by measuring the roundtrip propagation time of short laser pulses, and other information such as surface roughness, slope, and reflectivity can also be measured through processing the received waveform. The Geoscience Laser Altimeter System (GLAS) on the Ice, Cloud and land Elevation Satellite (ICESat) with small footprint and along track spacing provided more accurate surface elevation and higher horizontal resolution on Earth observation, and was widely used to monitor ice sheet, sea ice, vegetation, and other changing environments [13].
A Nd:YAG laser of fundamental mode was used in the GLAS system, and the cross section of the transmitted laser pulse satisfies a 2dimensional Gaussian distribution, which means the majority of total energy focuses on the small central region, and decreases rapidly with the rising radius. Under this kind of energy distribution, the reflection of surface covering which locates on the center of the laser footprint will play a much bigger role in the received waveform. The energy distribution on the crosssection of a flattened Gaussian laser beam is more uniform than that of the fundamental mode, and the energy is nearly a constant in more than half of the total crosssection under the higher order beam. It is more suitable to detect the complex coverings of the Earth’s surface.
Currently, the waveform model of fundamental mode laser used in GLAS was established by Gardner, Tsai, etc [46]. In this paper, the expressions of the main parameters of the laser altimeter waveform under flattened Gaussian laser were derived, and these results were used to compute total received energy, pulse width, variances of the propagation time delay and ranging error. A waveform simulator and parameter extraction program were used to verify the analytical expressions, which were universally used for both the flattened Gaussian beam and conventional fundamental mode ones. The effects of laser speckle, shot noise, and surface profile of the target were considered, and the influence of varying laser mode was quantitatively analyzed.
II. FLATTENED GAUSSIAN BEAM
For the fundamental mode laser, the laser crosssection and waveform of the transmitted laser pulse are both assumed to be Gaussian in shape:
Where
ρ is the radius on the target surface measured from the center of the laser footprint,z is the altitude of the altimeter measured from target surface,Q is the total energy of the transmitted laser pulse,θ_{T} is the divergence angle of the laser beam, andσ_{f} is the pulse width of the transmitted laser pulse. For the fundamental mode lasers, the laser energy decreases rapidly as the radius increases, while for the highorder mode flattened Gaussian, the energy is nearly identical on the central crosssection of the laser footprint. According to the description of the reference [7, 8], the crosssection of a highorder mode flattened Gaussian beam can be represented as:According to Eq. (1), Eq. (3) is normalized processed, and is shown in Eq. (4).
III. ALTIMETER RECEIVED MODEL
The transmitted laser pulse undergoes Fresnel diffraction during propagation, and then it is reflected by the Earth’s surface. After another Fresnel diffraction, the incident laser light on the telescope is focused onto the photoelectric detector, and the received signal power is given by: [9]
η is the efficiency of the receiver optics and detector,
A_{R} is the area of the receiver aperture,T_{a} is oneway intensity transmittance of the atmosphere,hv is the energy of one signal photon,β(ρ) is the intensity reflection coefficient of the ground target, atρ , andψ is the time delay accounted for surface fluctuation and footprint reflection, where:ξ(ρ) is the surface elevation, atρ , andc is the light speed in a vacuum.The essential parameters of the received waveform consist of the expected number of detected signal photons denoted by
N , the mean square width of the received pulse denoted byσ_{s} , the received signal time delay denoted byT_{s} , and the variance ofT_{s} denoted by Var(T_{s} ), which affect the total received energy, waveform shape, received pulse centroid, and the uncertainty of laser ranging, respectively.If we assume the reflection coefficient
β(ρ) is a constant within the laser footprint, it can be denoted simply byβ . If the target surface with fluctuant profile is quasiflat with slope, and by generalizing the theory of Filin [10], the surface slope in the local level coordinate system should be converted into ‘modified surface slope’φ in altimeter coordinate system, which is determined as a function of the satellite attitude, laser pointing and the surface slope itself. Then,N can be acquired by integrating signal powerS(t) overt , and the integration result is:Eq. (7) indicates that the received energy is irrelevant to transmitted laser mode, but related to transmitted laser energy, which means if the transmitted laser energy, the target surface, and the parameters of atmosphere and detector are the same, altering lasers from fundamental mode to flattened Gaussian does not actually change
N .σ_{s} is expressed as the secondorder normalized central moment, defined by: [9]T_{s} can be calculated usingT_{s} = 2z (1+tan^{2}θ_{T} )/c . We substitute Eqs. (1), (4) and (5) into Eq. (8) to derive the expressions of mean square width of fundamental mode and flattened Gaussian laser, and the results are given by Eq. (9) and Eq. (10), respectively.In both Eq. (9) and Eq. (10), the received pulse widths are composed of 4 items. The transmitted pulse width is described as the first term, and the broadened effects produced by surface roughness of the ground target, modified slope, and the beam curvature are shown as the second, third and fourth terms, respectively. For a beam divergence of tens of μrad on a common laser altimeter system, the broadening effect of the last term is relatively weak and can be ignored. Actually, when the order in Eq.(4) is
k = 0, which means flattened Gaussian laser changes to fundamental mode, Eq. (9) and Eq. (10) are identical, so Eq. (10) is common to both fundamental mode and flattened Gaussian laser.According to the reference [9], the variance of
T_{s} can be expressed in Eq. (11), and the first term caused by shot noise decreases as the signal energy increases, while the second term due to speckle can only be reduced by increasingK_{s} .K_{s} is the ratio of the receiver area to the speckle correlation, and is usually called the speckle signaltonoise ratio, given by Eq. (12). Whereλ is the laser wavelength andF is the excess noise factor of the detector.Under weak signal circumstances such as a spaceborne laser altimeter system, shot noise is the dominant noise source. In Eq. (11),
b_{n} (p ,z ) is the normalized effective cross section of the laser beam.By substituting Eqs.(1), (4) and (13) into Eq. (11) and integrating over ρ, the variance of
T_{s} can be expressed as Eq. (14) for fundamental mode laser and Eq. (15) for flattened Gaussian laser. In addition, if the order isk = 0, the two expressions will be the same.As for the GLAS system, the satellite altitude is
z = 600 km, the telescope diameter isd = 1 m, the laser wavelength isλ = 1064 nm, the divergence angle isθ_{T} = 110 μ rad, the transmitted energy of the laser isQ = 75 mJ, and the transmitted pulse width isσ_{f} = 2.37 ns. Then,K_{s} is nearly 10^{6}, which is considerably larger than N which is approximately 10^{4}, so the effect due to speckle can be ignored.IV. RESULTS
According to the document related to the laser altimeter simulator given by NASA [11], the waveform simulator software was programmed in MATLAB codes, and was improved to run under a flattened Gaussian laser. The received waveform could be simulated by substituting the system parameters and assumed or actual surface profiles. Firstly, according to Eq. (1), the continuous laser waveform is discretized into a time series, and the sampling interval is 0.1ns typically. Then, as shown in Eq. (2) or (4), the cross section intensity of the transmitted beam is Gaussian or flattened Gaussian in the farfield, and on the target surface, the total energy of the incident beam is divided into a finite number of grids, 0.2 m × 0.2 m for a typical GLAS laser footprint. The reflected intensity and range delay are calculated independently for each grid, and the height and diffuse reflectivity can be specified independently for each grid on the surface profile. We use the GLAS system parameters as shown above, and set grid bin of surface profile to 0.2 m × 0.2 m and time bin to 0.1 ns, then the simulated waveforms are shown in Fig. 1.
In Eq. (7), we proved that if the transmitted laser energy and measuring environment were kept the same, the laser mode would not affect the received laser energy. Consequently, in Fig. 1 the received laser energy of different laser modes were normalized. When modified slope is
φ = 0 in Fig. 1 (a), which means normal incidence and flat surface with zero slope, the simulated waveforms are nearly the same in different modes of lasers, which corresponds well to the models in Eqs. (9) and (10). The reason is that the first two terms are the same, the third term influenced by slope is input by 0 slope, and the last term is relatively weak in simulation. When the slope is φ = 0.05, which is approximately 2.86°, the waveforms of flattened Gaussian with different orders are illustrated in Fig. 1 (b), and zero order is the same as for the fundamental mode laser. In the comparison of Fig. 1, increasing slope will decrease the amplitude of received waveform and broaden the pulse width severely, especially for higherorder flattened beams.0.05 is a typical slope in the icesheet region, because the slope of icesheet is slight in inland Antarctica, approximately ~1:1000, and increases from inland to coastal area. The surface slope is less than 1:300 in more than half of the total region, and is less than 3:200 in 90% of the total. [12] As shown in Fig. 1, the pulse widths broaden as the order of flattened Gaussian laser increases, whose trend agrees with Eq. (10), while the peak values of waveforms are inversely proportional to the order.
In order to evaluate whether or not the simulated waveform precisely corresponds to the analytic expression shown as Eq. (10), we first should ensure that the pulse widths extracted from simulated waveforms are correct. So, one measured waveform by GLAS in Antarctica in March 2003 was chosen, illustrated in Fig. 2, to compare the main parameters given by NASA official GLA05 records and extracted by our waveform processing program based on the theory and method described in reference [6, 13]. In GLAS Algorithm Theoretical Basis Document, the transmitted and received waveform were treated as single Gaussian function or a sum of several, and the Gaussian parameters were extracted by Gaussian fit based on least squares algorithm.
The transmitted width with
σ_{f} = 2.72 ns (GLA05 i_parmTr) and received width withσ_{s} = 4.40 ns (GLA05 i_parm) were the same as the results extracted by our program. This waveform corresponded to the reflection from a very flat surface, which was typical of the Antarctica inland region. Then, we used the processing program to extract the pulse widths of the simulated waveforms which were illustrated in Fig. 1, and calculated the theoretical pulse widths according to Eq. (10). The comparisons between theoretical and simulated waveform widths under different order beams are listed in Table 1, and two groups of waveforms are drawn in Fig. 3.As shown in Table 1, the simulated and theoretical pulse widths are nearly the same with all biases less than 1ns or 3%, which verifies the expression in Eq. (10). The comparison between Fig. 3 (a) and (b) shows that under the same laser of fundamental mode, the pulse width of waveform reflected by a flat surface is only about onetenth than that of a surface with slope φ = 0.05, while the peak value is approximately ten times as high. Both shapes of waveform are Gaussian functions.
The domains of coordinate axes in Fig. 3 (b)(f) are the same. We can see clearly that the peak value decreases, the width increases as the order
k ascends. From Table 1, the width under the condition of slopeφ = 0.05 and orderk = 4 is nearly twice that of slopeφ = 0.05 and orderk = 0.Besides, as the order
k increases, the shape of waveform deforms, and cannot be supposed to be a strict Gaussian function in Fig. 3 (b)(f). We use the generalized Gaussian function, shown in Eq. (16), to fit these waveforms, and the new fitted shapes illustrated in Fig. 3 (a)(f) by red dash lines.Where
A is the amplitude,B is the center of Gaussian function,n is the power of generalized Gaussian function and σ is the generalized Gaussian width. As we can see in Fig. 3, compared to Gaussian function fit, the generalized Gaussian model significantly improves the fit results. Generalized Gaussian fittings influence the amplitudeA in Eq. (16) slightly, however, as shown in Table 2, owing to the rising ofn in Eq. (16), which serves as the power parameter, the widths increase drastically as the beam order rises. And, the derivation process in this paper is irrelevant to the generalized Gaussian function, so the physical interpretation needs to be researched in the future.When laser ranging is based on the time of flight, the ranging error or ranging precision arising from the laser itself and the target will be the product of uncertainty of timing and half of the velocity of light, and it can be expressed in Eq. (17), where ‘std’ means standard deviation.
We substitute
N , derived in Eq. (7), into Eqs. (15) and (17). As analyzed above, the shot noise serves as the primary effect in the spaceborne laser altimeter system, so after simplification we get the expression for ranging error.If we assume atmospheric transmittance is
T_{a} = 0.7, detector efficiency is η = 0.5, excess noise factor isF = 5, reflection coefficient isβ = 0.3, the expected photons will beN = 20478, which is approximately 3.83 fJ and agreed with the result given in reference [14]. Under any orders of flattened Gaussian lasers, when slope isφ = 0, the ranging error will be less than 1 cm and can be nearly ignored. However, under 04th beam orders, when slope isφ = 0.05, it will become 5.2 cm, 6. 8cm, 8.2 cm, 9.3 cm and 10.3 cm, respectively. Consequently, the ranging error also increases with the rising of mode orderk and slopeφ , and if the slope isφ = 0.05, the error of the 4^{th} order flattened Gaussian laser will be twice than that of fundamental mode laser.In the above content, all the targets are assumed to be a single quasiflat surface, however, the actual target would be much more complicated. As a result, if we assume that there are two quasiflat surfaces on the target, which is shown in Fig. 4 (a) and can be considered as a building with 5m height on the ground. In Fig. 5 (a), the simulator shows the corresponding waveforms of fundamental mode and 4^{th} order flattened beam. Then, we assume that there are two buildings with different heights and a pit on the ground, shown in Fig. 4 (b), and finally the ground is changed to be an inclined plane with 0.01 slope, shown in Fig. 4 (c). As for GLAS parameters, the divergence
θ_{T} is equal to 110 μrad, which means that more than 99% of laser energy locates within a 200 m radius circle.The transmitted laser energies are the same, and we can clearly find that the flattened beam get much better performance for a multisurface target. In Fig. 5 (a), both lasers indicate that there are two surfaces on the target, but the first waveform peak of fundamental mode with approximately 0.01 amplitude, which corresponds to the 5 m height building, is much less than that of flattened beam with over 0.03 amplitude. And the time interval between two peaks of waveform could invert the distance of two surfaces of target. In Fig. 5 (b) and 5 (c), when the areas of new surfaces become smaller and further from the laser center, we can hardly recognize the corresponding peaks in waveform, and the last peak of the fundamental mode has completely disappeared. The reason is that the energy of the fundamental mode decreases rapidly with the radius of laser footprint rising, while flattened Gaussian laser beam is more uniform than that of fundamental mode and can reflect the effects of target surface far from the laser center.
Although the ranging error of flattened Gaussian beams is bigger than that of the fundamental mode, typically for 10.3 cm under 4^{th} order compared to less than 1 cm with 0.05 slope, the flattened beams could detect more details of the target. And if the surface is flat without slope, the ranging error is the same under flattened Gaussian beams and the fundamental mode. Besides, the derived ranging error expression in this article is arising from the laser itself and the target, which is only part of the total ranging error, and another main influencing factor is pointing errors, 1.5″ of which makes 21.8 cm ranging error with 0.05 slope for typical GLAS parameters [15, 16]. According to law of error propagation, the total ranging error is equal to the root of the sum of every error squared.
V. CONCLUSION
The main parameters of received pulse under flattened Gaussian laser were derived and verified through waveform simulator and processing, and the bias between simulated and theoretical pulse widths were less than 3%. These new expressions of flattened Gaussian laser were universal compared to the former fundamental mode, and could be used not only under flattened Gaussian lasers but also fundamental mode when the order was set to be zero.
Both the pulse width and ranging error increase as the modified slope and flattened Gaussian order grow; therefore, laser altimeter will get more accurate ranging performance under small nadir angles, meaning nearly normal incidence, and small surface slope. As for received waveform of higher order flattened Gaussian, the shape deforms and is no longer a strict Gaussian function, because the incident energy of laser footprint does not satisfy 2dimensional Gaussian distribution. Generalized Gaussian model was used to express these curves and got very good fit results, and the physical interpretation of generalized Gaussian model needs to be researched in the future.
The flattened beam got much better performance under the circumstance of multisurface target, especially when the small surface is far from the center of the laser footprint, and it is more suitable to detect the complex coverings of the Earth’s surface. In addition, the established expressions and waveform simulator are meaningful to the system design of the flattened Gaussian laser altimeter.

6. Brenner A. C., Zwally H. J., Bentley C. R., Csatho B. M., Harding D. J., Hofton M. A., Minster J.B., Roberts L., Saba J. L. 2011 Derivation of Range and Range Distributions from Laser Pulse Waveform Analysis for Surface Elevations, Roughness, Slope, and Vegetation Heights, AS Algorithm Theoretical Basis Document Version 5.0

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[FIG. 1.] Simulated waveforms of flattened Gaussian with different orders. (a) Surface with 0 slope; (b) Surface with 0.05 slope.

[FIG. 2.] The actual waveform and processed result of GLAS.

[TABLE 1.] Contrast between theoretical and simulated waveform widths under different order beams

[FIG. 3.] The contrast of original, fitted Gaussian, and fitted generalized Gaussian waveforms with different slopes and beam orders. (a) Modified slope with φ=0 and beam order with k=0; (b)(f) Modified slope with φ=0.05 and beam orders with k=0, 1, 2, 3, 4.

[]

[TABLE 2.] Parameters fitted by generalized Gaussian model

[]

[]

[FIG. 4.] Laser footprint illuminating different target surfaces. (a) Target with two flat planes; (b) Target with four flat planes; (c) Target with three flat surfaces and an inclined planes.

[FIG. 5.] Received waveform of different target surfaces (a) Waveforms corresponding to two flat planes; (b) Waveforms corresponding to four flat planes; (c) Waveforms corresponding to three flat surfaces and an inclined planes.