Next Article in Journal
Traffic Light Recognition Based on Binary Semantic Segmentation Network
Previous Article in Journal
Characterization of Gas Absorption Modules Based on Flexible Mid-Infrared Hollow Waveguides
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Laser Vegetation Detecting Sensor: A Full Waveform, Large-Footprint, Airborne Laser Altimeter for Monitoring Forest Resources

1
College of Forestry, Beijing Forestry University, Beijing 100083, China
2
Academy of Forest Inventory and Planning, National Forestry and Grassland Administration, Beijing 100714, China
3
USDA Forest Service, Northern Research Station, Newtown Square, PA 19073, USA
4
Geography Information and Tourism College, Chuzhou University, Chuzhou 239000, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(7), 1699; https://doi.org/10.3390/s19071699
Submission received: 26 January 2019 / Revised: 4 April 2019 / Accepted: 8 April 2019 / Published: 10 April 2019
(This article belongs to the Section Remote Sensors)

Abstract

:
The use of satellite-borne large-footprint LiDAR (light detection and ranging) systems allows for the acquisition of forest monitoring data. This paper mainly describes the design, use, operating principles, installation and data properties of the new Laser Vegetation Detecting Sensor (LVDS), a LiDAR system designed and developed at the Academy of Forest Inventory and Planning (AFIP) and the Beijing Institute of Telemetry (BIT). Data from LVDS were used to calculate the mean height of forest trees on sample plots using data collected in the Hunan province of China. The results show that the full waveform data obtained by LVDS has the ability to accurately characterize forest height. The mean absolute percentage error of mean forest height per plot in flat areas was 6.8%, with a mean absolute deviation of 0.78 m. The airborne LVDS system provides prototype data sets and a platform for instrument proof-of-concept studies for China’s Terrestrial Ecosystem Carbon Monitoring (TECM) mission, which is an Earth remote sensing satellite due for launch in 2020. The information produced by LVDS allows for forest structure studies with high accuracy and coverage of large areas.

1. Introduction

Forests play an important role in the ecosystem. With the development of China’s economy, forest resource monitoring has become more and more important as pressures on forest resources have increased. Mean canopy height is an important parameter that reflects the vertical structure of forests [1,2], and is of great importance for the quantitative estimation of forest habitat extent [3], biomass quantity, and the area affected by landscape changes [4,5,6]. In traditional field-based forest surveys, ground measurements on inventory plots have been used to estimate mean forest height. Not only is this method time- and effort-consuming, but also tree height collection is often unreliable, especially in areas where forests are dense or otherwise difficult to traverse. Therefore, traditional forest survey instruments cannot meet the need for low-cost, reliable forest canopy height information.
Optical remote sensing images from passive sensors have been used extensively to estimate the horizontal distribution of forest. For example, spaceborne sensor systems like the Landsat series [7,8,9], MODIS [10,11,12,13], GF series [14,15,16], ZY-3 [17,18], and Sentinel series [16,19,20] provide spectral information that is used as the foundation for image classification globally. However, due to the complex spatial structure of forests and limitations of optical data, forest cover maps from these sources often suffer from inaccuracy [21], particularly for the classification of vertical forest structure [22]. Microwave remote sensing, a technology that uses an active sensor system, has often been used for tree height estimation, but its signal tends to saturate in forests with dense canopies [23]. Compared to other technologies, light detection and ranging (LiDAR) have been shown to be very useful for detecting vertical structural attributes like canopy height and tree crown density in forested areas [24]. The advanced technology, coupled with the GNSS (Global Navigation Satellite System) and IMU (inertial measurement unit), can improve the horizontal and vertical positional accuracy of airborne remote sensing data to the decimeter level [25].
The first Scanning LiDAR Imager of Canopies by Echo Recovery (SLICER), which had a 15 m footprint, was designed by NASA in the early 1990s. A series of studies about the equipment and data products followed [26,27,28,29]. A later advance in laser remote sensing was the launch of the Laser Vegetation Imaging Sensor (LVIS) [30], designed as the airborne simulator for the Vegetation Canopy LiDAR (VCL) mission. Performance tests and analyses of the properties of these data were conducted, leading to a growth in interest in the use of laser imaging technology [31]. In January 2003, the Geoscience Laser Altimeter System (GLAS), carried on the ICESat instrument, was launched. This sensor provided full waveform LiDAR data for a set of well-distributed, 70-m diameter footprints for much of the globe [32,33]. GLAS data were used in various ways to calibrate models of vegetation structure parameters [34,35,36,37,38,39,40,41]. The Global Ecosystem Dynamics Investigation (GEDI), carried by the SpaceX Commercial Resupply Mission 16, launched successfully on December 5th, 2018 from Cape Canaveral in Florida. GEDI will map ecosystem structures including carbon and nutrient cycling, habitat quality and biodiversity, forest health and productivity, hydrologic cycling, and effects of natural and human caused disturbances. At the same time, some researchers have proved the GEDI’s ability to estimate ground elevation and forest parameters (maximum forest height and forest aboveground biomass) through simulated GEDI data [42,43,44,45].
The successful application of spaceborne LiDAR in forest monitoring provides an impetus for researchers in China to improve the country’s forest resource monitoring systems. Motivated by this interest, the Terrestrial Ecosystem Carbon Monitoring Satellite Project, which consists of a large-footprint LiDAR system, was formally approved. The mission will launch at the end of 2020, with the goal of providing continuous forest resource monitoring of China. In order to test the sensor system prior to implementation of the spaceborne instrument, a ¼ scale, large-footprint airborne LiDAR system based on the spaceborne system was designed by the Academy of Forest Inventory and Planning, State Forestry Administration of China.
To the best of our knowledge, no study has introduced this new airborne large-footprint LiDAR system and tested the application of the system to estimate characteristics of different landscape features such as forest height. The aim of this study, therefore, was to deploy and test a newly designed large-footprint LiDAR system that was developed to test technology that will eventually be deployed in space. In forested areas, we examined the relationship between LiDAR height estimates and mean forest height based on ground plot data. Our goal was to better understand the capabilities of the sensor and to develop a validation procedure for this type of data.

2. The Laser Vegetation Detecting Sensor (LVDS)

2.1. System Overview

LVDS is a pulsed laser altimeter, which measures the distance between the instrument and the target surface. The entire time history of the outgoing and return laser pulses is digitized using a single detector, digital converter, and timing clock to clearly describe the vertical distribution of surfaces within each laser footprint. The main system is composed of a laser sensor unit, a signal processing and control unit, a two-dimensional stabilized platform, a GNSS/INS integrated navigation unit, and a charge coupled device (CCD) camera unit (Figure 1).
The signal processing and control unit controls the emissions from the single-wavelength solid state laser, as well as conducts proportional beam splitting of the transmitted beam through the optical beam splitting device. A low energy beam is coupled to the relay optical device through a fiber optic cable. The laser beam is then coupled to the detector by the coupling lens of the relay optical device. A high energy beam is then expanded and sent to the ground by the transmitting optical beam splitting device. The laser light signal reflected by the ground targets is then received by the telescope. This received signal is filtered and focused by the relay optical device, and the detector transforms the laser light signals into electrical signals.
The signal processing and control unit collects, processes and uploads these electrical signals, collects data on the working state of system, and sends these to the real-time processing and display software of the upper computer. The upper computer’s real-time processing and display software integrates and processes the received signals from the sensors and control unit, allowing for real-time data display and storage.
The GNSS/INS combined navigation unit and laser sensor unit are affixed to the frame of the unit and used to obtain the current position and attitude data of the laser sensor unit in real-time. The laser sensor unit and aerial camera are mounted on the two-dimensional stabilized platform. The two-dimensional stabilized platform holds the gyroscope sensor and laser sensor units. Based on the real-time attitude data of the gyroscope sensor, the two-dimensional stabilized platform attitude is adjusted dynamically, so as to keep the appropriate pointing angle of the laser and to absorb shock experienced by the laser sensor unit. The aerial camera points in the same direction as the laser sensor unit to complete the image acquisition of ground targets.

2.1.1. Laser Sensor Unit

The laser sensor unit (Figure 2) consists of a single-wavelength solid-state laser, a transmitting optical beam splitting device, a telescope, a relay optical device, and a detector. By applying narrow pulse laser emission technology, the single-pulse energy is reduced and the echo signal-to-noise ratio will be high. At the same time, the overall application safety and reliability of the system can be improved, and the volume of the system (1000 × 870 × 550 mm) can be effectively reduced. The output wavelength of the laser is 1064 nm, the single-pulse energy is 2 mJ, the pulse repetition frequency is 40 Hz, and the pulse width is 1.5 ns. The beam splitting ratio of the transmitting optical beam splitting device is 1:100. The divergence angle of the emitting laser beam after beam expansion is 5 mrad, and the footprint size is 15 m at altitudes up to 3 km above ground level.
The receiving telescope adopts the transmitted configuration, with the advantage of no center shielding in the field of view, high optical efficiency, compact structure, and reducing the difficulty of machining and setting. The telescope consists of four BK7 lenses, with an effective diameter of 100 mm, a total length of 298.8 mm, a receiving field of view angle of 6 mrad, an effective focal length of 148.8 mm, and an image surface diameter of 3 mm.
The relay optical unit is mainly composed of a speculum, a collimating lens, a focusing lens and a detector triggering light path. The system length is shortened by placing a speculum in front of the focal length of the telescope. The optical spot diameter of the telescope mirror surface is about 3 mm, and the diameter of the photosensitive surface of the detector is about 0.8 mm. Collimating and focusing the echo signal received by the telescope can reduce energy loss and improve the detection efficiency of the system. The detector in the relay optical unit triggers the light path, processes the energy introduced by the fiber optic, and triggers the detector. The field of view of the collimating lens was designed to be ±20 mrad, and the diameter of the image surface (3 mm) is consistent with that of the telescope mirror surface. The collimating lens is composed of three BK7 lenses, which has a total length of about 100 mm and a posterior intercept of 62 mm. The imaging quality is close to the diffraction limit. The function of the focus lens is to focus all echo signals into the photosensitive surface of the detector. The focus lens is made up of three BK7 lenses, with a total length of about 25 mm and a diameter of 0.3 mm for the light spot on the image surface. The trigger optical channel can adjust the trigger energy, simplify the structure design and shield stray light interference by placing a customized prism at the rear end of the collimating lens.
The detector is used to convert the laser echoes received by the optical receiving system into electrical signals for the measurement of signal amplitude, waveform and flight time. It uses a silicon avalanche diode photoelectric receiving module as the detection unit. A semiconductor cooling unit is integrated into the sealed package to improve the heat dissipation performance and stability of the detector. The diameter of the photosensitive surface of the detector is 0.8 mm, and its responsiveness to the wavelength of 1064 nm is 200 KV/W.

2.1.2. Signal Processing and Control Unit

The signal processing and control unit (Figure 3) is composed of a high-speed waveform acquisition collector, data processing and interface controller, state acquisition controller, controller, upper-computer real-time processing display software, and a laser controller. All units are integrated into a single structural platform, which effectively reduces the number and complexity of devices and improves their reliability, convenience, and usability.
The laser controller controls the laser’s emission time, frequency, and energy. The controller is responsible for the signal processing and control unit, command sending and receiving, data storage, and output of synchronous signals. The data storage capacity can reach 1 TB. The high-speed waveform collector completes high-speed sampling of the electrical signal output by the detector, triggered by the system synchronization signal output by the controller. The sampling rate is 1 GHz, and the effective sampling bits can reach 10 bits. The data processing and interface control complete the processing and packaging of sampled data, and upload measurement data in real-time through an ethernet interface and controller. The state acquisition and controller are used to monitor the system state in real time and complete the setting of working parameters and uploading of working state detection. The upper-computer real-time processing display software is installed in the controller system for system parameter configuration, system switching machine control, real-time display of state monitoring parameters, and real-time display of sampled waveforms.

2.1.3. Two-Dimensional Stability Platform

The two-dimensional stability platform is designed with a two-dimensional orthogonal rotating shaft structure. The two-way rotating shaft is driven by an independent servo motor. The system integrates the two-axis gyro inertial sensor (INS) and attitude automatic controller, which can measure the attitude change of the laser radar installation platform in real time and calculate the motor rotation angle to be adjusted in real time. The two-dimensional stability platform can respond quickly to the attitude change of the system platform to ensure the stability of the system during the measurement process. At the same time, the shock absorber is connected at the bottom of the high stability platform to filter high-frequency vibration due to aircraft flight motion and reduce its impact on measurements. The stability range (azimuth and pitching axis) is greater than or equal to ±8 deg, the maximum angular velocity (azimuth and pitching axis) is greater than or equal to 50 deg/s, the maximum angular acceleration (azimuth and pitching axis) is greater than or equal to 200 deg/s2, the dynamic error (azimuth and pitching axis) is less than 0.08 deg (3 deg, 3 Hz), the vertical deviation is less than 0.01 deg in automatic leveling mode and the effective load is 50 kg.

2.1.4. GNSS/INS Integrated Navigation Unit

The INS module provides platform roll, pitch, and bearing information at a 200 Hz sampling rate with 0.005 deg resolution. The GNSS module provides real-time geographic position of the instrument which is fed into a pilot assistance system that allows precision flying (to within 0.008 deg) along predetermined flight tracks. Data from the on-board GNSS are processed post-flight with data from ground-based, GNSS receivers to provide an airplane trajectory with vertical accuracy of ~30 cm.

2.1.5. CCD Camera Unit

A Phase One iXU-R180 was selected as the aerial camera. Its main performance parameters are shown in Table 1. Further details on the camera specifications can be found at https://industrial.phaseone.com/.

3. Calibration of Prototype

3.1. Laboratory Test of Laser Emission Pulse

A PIN (positive intrinsic-negative) photodetector was used to measure the laser pulse waveform and the laser pulse duration in the laboratory (Figure 4). The laser pulse size and energy distribution characteristics were measured with the array CCD detector.
The results of the high-speed PIN photodetector test show that the full width at half maximum of the laser emission pulse is 1 ns. By using array CCD detector, the power/energy and the space of the laser emission pulse captured by the CCD are showed in Table 2, which reflects the quality of the laser emission spot is good. The laser beam on the X-axis and the Y-axis are Quasi-Gaussian pulses. Divergence angle of the laser beams is calculated by CCD camera technique measuring spot size, and the result is 5 mrad.

3.2. Laboratory Test of the Laser Echo Signal

A second laboratory test was performed to assess the distance measurement capabilities of the laser-sensor system. By pointing the laser at a fixed hard target with a distance of 944.1 m, the time delay between laser pulse emission and detection was measured (Figure 5). After more than 10 min of cumulative measurement, 450 pulses were randomly selected for analysis using MATLAB R2014a software manufactured by MathWorks Inc, Natick, MA, USA.
The analyzed result based on MATLAB R2014a showed that the noise (jitter) of echo delay time was less than 1 ns (Figure 6), and the ranging accuracy was thus less than 0.15 m (Figure 7), which indicated that the large-footprint LiDAR equipment had good performance.

4. Experimentation Outdoors

4.1. Field Testing of the Instrument

As part of the calibration and validation activities for the TECM mission, LVDS was used in two operational tests, the first of which was conducted in Shijiazhuang, Hebei province in April 2017, and the second in the Northwestern Hunan province in December 2017. This paper only addresses results from the second test.
The aerial campaign for the second LVDS forest resources detection test began at the Lotus airport near Zhangjiajie city, which is in Hunan Province, in December 2017. Two flights were carried out. The first flight was between 10:52 and 15:35 on 24 December 2017. During the 283 min flight, 19 files containing a total of 28 GB of large-footprint LiDAR data were collected. The second flight occurred between 11:11 and 16:22 on 25 December 2017. During this 191 min flight, 43 files contained 33 GB of large-footprint LiDAR data were collected. The flight routes were chosen in order to cover landscape types and species assemblages of interest in Taoyuan, Huayuan, Baojing, Yongshun, Yuanling, and Guzhang counties, and Jishou City in Northwestern Hunan province, and were more than 980 km in length. The planes flew at elevations of 3000 m over the ground.
In the field sample plot data collection, forest parameters including species type, canopy height, and diameter at breast height (DBH) were measured. A total of 306 circular sample plots with a diameter of 15 m (with an area of 177 m2) were measured in mixed forests between December 2017 and January 2018. The DBHs of all trees greater than 5.0 cm were measured in each plot. In order to improve the accuracy of the investigation [46], plot centers were geographically located using a Galaxy 6 RTK GNSS receiver manufactured by SOUTH Inc, Surveying and Mapping, Guangzhou, China. These units have a horizontal accuracy better than 0.25 m. For each of the sample plots, GNSS points were collected for more than 30 min using a one second logging rate. Forest canopy heights were measured using a VERTEX LASER VL5 manufactured by Haglof, Sweden [47]. These devices have a vertical accuracy of 0.1 m. All geographic data were collected using the WGS84 datum and the UTM zone 49N coordinate system (Figure 8).

4.1.1. Estimation of Mean Forest Height Based on LVDS Data

The first step in the data analysis workflow for the forest height estimation tests was to pre-process the raw LVDS data. This included reducing the data volume through an approximately 1/100 subsample. First we removed acquisition sites that were not in forest or were not on flat ground, and then we subsampled the dataset to arrive at a total of 36 sites that met criteria in Taoyuan County; these were the points used in subsequent steps, including LVDS height estimate validation.

Large-Footprint LiDAR Data Processing

To calculate the mean forest height for each LVDS site, the waveforms were pre-processed, which included noise estimation, threshold calculation, identification of wave initiation, and bare ground identification. In order improve on processing speed, a 500 points subset around the maximum return value (200 points before and 300 points after) was chosen, as this range consistently contained the relevant information needed for canopy height estimation. Hereafter, all references to the waveform will be to this subset dataset.
(1) Noise estimation
Before the device detects the signal from the ground or the crown reflection, the background noise of the signal is received. After analyzing 200 waveforms, we developed the following heuristics with which to identify the location in the waveform at which the signal associated with reflectance from an object on the ground begins and ends. First, the mean and variance of the return energy associated with noise were estimated before the relevant part of the return signal by using the first 50 observations comprising the waveform. That is:
N ¯ b = i = 1 50 w i 50
and:
S N b = i = 1 50 ( w i N ¯ b ) 2 49
where wi is the voltage value of the ith return value of the waveform. Next, the mean and variance of the noise after the relevant portion of the waveform signal were calculated by using the last 50 data points of the waveform signal:
N ¯ e = i = 451 500 w i 50
and:
S N e = i = 451 500 ( w i N ¯ e ) 2 49
(2) Signal threshold calculation
The noise detection threshold T used for LiDAR signal analysis is commonly chosen by multiplying N b ¯ and N e ¯ by an appropriate number n of standard deviations [48]. The appropriate value of n varies in different studies (n = 2 [49], n = 3 [50], n = 3.5 [51], n = 4 [52], n = 4.5 [53]). In this study, n = 2 was used, i.e.,
T b = N b ¯ + 2 S N b
and:
T e = N e ¯ + 2 S N e
Tb was used to extract the initial position of signal, while Te was used to extract the ground position.
(3) Estimating initial signal position
The initial position of the waveform signal (Sb) represents the location of the highest tree in the study area. It defines the location of the signal’s starting position, and was defined as where the values of three consecutive frames (from the front to the back) are all larger than the begin-threshold value. That is, for each LiDAR return wj (j = 1…500):
Sb = wj if (wj-2 < Tbwj-1 < Tbwj < Tb) is true
(4) Estimating ground position
Due to a high proportion of canopy cover and low ground solar energy reflection in the study area, the return energy diminishes when the pulse is reflected by the ground. Using the same logic as above, it follows that the signal associated with the ground (Sg) is where the values of three consecutive frames (from the back to the front) are all larger than the end-threshold value. That is, for each LiDAR return wj (j = 1...500):
Sg = wj if (wj-2 > Tewj-1 > Tewj > Te) is true.
(5) Estimating mean forest height
Mean forest height associated with a laser pulse (HLVDS) is calculated by multiplying half of the distance light travels in one nanosecond (30 cm) by the difference between Sg and Sc. The initial signal position is the emission energy associated with reflection from the top of the canopy. As the canopy area increases, more and more energy is returned, and reaches the peak at the mean forest height location. Therefore, the first peak position (Sc) after the initial signal position is defined as the mean forest height position. From the mean forest location and the ground location, mean forest height can be estimated as (Figure 9):
HLVDS = Rv(SgSc)
HLVDS is the mean forest height and Sg is the location of the recorded time of the ground position, Sc is the recorded time of the mean forest location, and Rv is a resolution of 15 cm (the signal can, thus, be written as the time interval of 1 ns, which gives a resolution of 15 cm.).

Height Modeling

The height modeling workflow is describe in overview in Figure 10. Although the center locations of the footprints for the large-footprint test site acquisitions were georeferenced, it was impossible to collocate our ground plots with the large footprint LiDAR sites (e.g., Figure 8) in order to calibrate height models or validate estimates. Therefore, we used height estimates modeled from small-footprint LiDAR data collected from another study as inputs to the LVDS validation process. Small-footprint LiDAR data offers the unique opportunity to bridge the gap between large-footprint LiDAR data and sample plot data. The existing small-footprint LiDAR data with which we had to work consisted of a 7 km × 90 km study site that covered the LVDS acquisition area. These data were collected on 20 December 2017 with a point density of 5.3 returns/m2. Georeferenced point clouds containing information on height above ground level were created using standard processing methods.
GIS files of the ground sample plot footprints (where tree heights were measured) were used to create subsets of the small-footprint LiDAR data contained within them. For each of several candidate small footprint LiDAR within-plot point cloud height value percentiles, simple linear regression models of the relationship between sample plot-based mean forest height (dependent variable) and small footprint LiDAR height value percentiles (independent variable) were constructed using 306 ground plots. A total of 214 of the 306 samples (96 samples whose dominant species were coniferous and 118 samples whose dominant species were deciduous) were used as the training samples and another 92 samples (41 coniferous and 51 deciduous) were treated as test samples for accuracy evaluation. After assessing model fit statistics for the candidate point cloud height percentiles, the 70th percentile point cloud height value was selected as a surrogate for measured mean forest height on a plot, and the associated regression models estimating forest height from the small footprint LiDAR point cloud data were used in subsequent analysis.
Mean forest height was estimated using the linear regression model from the previous step within the boundary for each of the 36 LVDS sites in our test. To assess accuracy of the height model, Mean Absolute Deviation (MAD) and mean absolute percentage error (MAPE) were calculated as follows:
MAD = 1 n i = 1 n | x i x ¯ |
and:
MAPE = 1 n i = 1 n | x i x ¯ x i | · 100
where x i is the deviation (ground plot canopy height—70th percentile height) associated with plot i and x ¯ is the average of x i across all 36 observations.

Comparison of LVDS and Height Modeled from Small Footprint LiDAR Regression Equations

The chosen linear model of the regression of mean forest height (HSFL) on the 70th percentile height value of the small-footprint LiDAR data ( X 70 th ) was of the form
H S F L = 0.8822 · X 70 th + 2.3147
The model’s R2 was 0.92 and the root mean square error (RMSE) was 1.58 m (Figure 11).
Data from 92 sample plots not used to develop the height model were used to evaluate the estimation accuracy of Equation (12). Using deviations calculated from the difference between the LiDAR-estimated and actual (ground) heights at these 92 locations, the MAD was found to be 1.32 m and the MAPE was found to be 9.4%. It is difficult to say, however, if the deviations occurred due to inaccuracies in the ground height measurements or the small-footprint LiDAR data aggregation approach. The relationship in Figure 11 does appear to be linear, however, with a nearly 1:1 relationship and no evidence of heteroscedasticity, suggesting that the model performance is acceptable.
In order to validate height estimates derived from the LVDS instrument (Equation (9), Figure 9), Equation (12) was applied to the 70th percentile of the small footprint LiDAR height values associated with the polygonal LVDS footprints (a 15 m diameter circle) for the 36 validation plots. In other words, the 70th percentile of the small footprint LiDAR height values within the LVDS validation plot footprint was calculated and compared to the LVDS-derived height estimate. Results are shown in Figure 12. The MAD was 0.78 m, and the MAPE was 6.8%. The relationship does follow the 1:1 line, and there does not appear to be any systematic bias or heteroscedasticity. These results strongly suggest that the LVDS system performs as expected when measuring forest canopy height, and can be a valuable tool for monitoring both the status of and trends in forest cover in China and elsewhere.

5. Discussion

5.1. The Advantage of the Equipment

As with other studies conducted with the same type of equipment [30], an advantage of this approach is the integration of a large laser spot emission and a large field of view. A laser beam is used to couple the laser pulse to a beam splitter through an optical fiber. The part of the beam with less energy is coupled to the relay optical device through the optical fiber. The remainder of the beam, which has more energy, is sent to the ground target after it is expanded. The placement of the refractive mirror in front of the telescope shortens the system length, collimating and focusing the reflected signal and maintaining its strength, improving the detection efficiency of the system. The light path through the detector uses a custom prism, allowing for the appropriate attenuation of the signal to correspond with the limits of the detector. In our system design, the combination of this prism and laser output attenuator thus adjusts the trigger energy.
In addition, the device adopts a high-precision time difference measurement mode governed by the transmitting and receiving optical paths. The laser emitting pulse interacts with the beam splitter through the optical fiber. The beam splitter lens diverts a small energy emitting trigger laser beam to a relay optical device through an optical fiber. The relay optical device sends the input emitting laser trigger pulse along the detector trigger light path. By using a custom prism, the transmitting pulse and the echo signal received by the telescope are coupled to detect the transmitted light pulses and echo signals directly, which avoids errors in the measurement of time difference caused by some factors, like laser dithering and system timing error, thus improving the accuracy of system time and distance measurement.

5.2. The Positional Accuracy of the Sample Plots

Some studies use handheld GPS units to locate sample plot centers [54]. It is generally known that with many commercially-available GPS units, the horizontal error is about 5 m where the signal is good [55,56]. In dense forest, the error may reach 30 m or more [57]. Therefore, if using handheld GPS to locate the sample plots with the aim of collocating them with the footprints of aerial laser pulses, large errors can occur. Real-time positioning that relies on GNSS units can be unreliable as well; when used outside the signal range, GNSS units’ positioning accuracy is not much different from that of a normal handheld GPS.
In our study, the total station and GNSS are used together to determine the position of sample plot center and all trees in the sample plot. First, GNSS1 was set up in a place with good mobile phone signal. After arriving at the sample plot site, the total station was set up at the center of the selected surveyed sample plot. The origin of coordinates of the total station was set as (0,0,0). The direction was selected and the virtual coordinate system was established. At the same time, GNSS2 and GNSS3 were set up in the relatively empty position of the sample plot, and the GNSS2/3 could clearly be seen by the total station. All the three GNSS base station units were selected static mode. All GNSS base stations worked at least 30 min at a time, and no more than 30 km apart. In addition, GNSS1 can obtain highly accurate location coordinates by selecting CORS (continuously operating reference stations) service. Taking the high accuracy coordinate as the control point, the differential data was processed. GNSS2 and GNSS3 within the sample plot have both the virtual coordinates collected by the total station and the high accuracy coordinates after the differential data processing. Through the coordinate transformation relationship, the accurate coordinates of the sample plot center and all trees can be obtained.
In other words, through setting up three GNSS base station units (one of which is located using the CORS service), the sample plot positional accuracy was greatly enhanced.

5.3. Accuracy of Large Laser Spot Data Validation

Due to the influence of the sensor itself, atmospheric scattering, and other factors, a certain amount of noise often exists in the received echo waveform. The accuracy of noise estimation and waveform data processing algorithms affects the estimation of mean forest height. Laboratory tests showed that the measurement jitter of echo delay time was less than 1 ns, and the range precision was better than 0.15 m for the measurement distance of a fixed target (Section 3.2). Based on these laboratory tests, customized noise detection thresholds were determined. Future research will assess the universality of these thresholds under different operating conditions.
Previous studies on extracting mean forest height included areas with sloping topography [39,58,59]. However, based on the need for a timely performance verification of the LVDS system, we devised a new method of calculating mean forest height over level terrain. It is worth noting that the method proposed in this paper is thus mainly suitable for level areas. Terrain slope will have an impact on the pulse broadening [52], thus directly affecting height estimation. Therefore, future research will account for this and include assessing the performance of the newly designed system in areas with undulating topography.

6. Conclusions

Our study of the LVDS system provides a test platform for the large-footprint LiDAR system to be used on the TECM satellite-borne LiDAR mission to be launched in 2020. We demonstrated through both laboratory testing of the individual components and field testing of the integrated system that the height measurements provided by the system are acceptable. The deviance between the instrument’s height measurements and those obtained in the field are small, and the impacts of this deviance on tree biomass estimation is small. We, therefore, conclude that data from the sensor system has much potential to improve forest monitoring and other land cover assessment projects in China and elsewhere.

Author Contributions

Y.H., F.W., and Z.S. conceived the idea of the study. Y.H. and F.W. conducted the field work. Y.H., Z.S., and A.L. processed the data and analyzed the results. All authors contributed to the writing and revision of the manuscript. Y.H. and F.W. are co-first authors of the article.

Funding

The authors gratefully acknowledge the financial support provided by the National Key Research and Development Program of China under grant no. 2016YFD0600205 and the Land ecological system carbon monitoring satellite ground comprehensive experiment under grant no. 2017-21-4 ** and the Key Project of Natural Science Research of the Anhui Education Department under grant no. KJ2017A413. We also thank the anonymous reviewers as well as the editors for their valuable comments improving the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ni, X.; Zhou, Y.; Cao, C.; Wang, X.; Shi, Y.; Park, T.; Choi, S.; Myneni, R.B. Height over Continental China Using Multi-Source Remote Sensing Data. Remote Sens. 2015, 7, 8436–8452. [Google Scholar] [CrossRef]
  2. Næsset, E. Determination of mean tree height of forest stands using airborne laser scanner data. ISPRS J. Photogramm. Remote Sens. 1997, 52, 49–56. [Google Scholar] [CrossRef]
  3. Ashcroft, M.B.; Gollan, J.R.; Ramp, D. Creating vegetation density profiles for a diverse range of ecological habitats using terrestrial laser scanning. Methods Ecol. Evol. 2014, 5, 263–272. [Google Scholar] [CrossRef] [Green Version]
  4. Nelson, R.; Hill, R.; Suárez, R. Model effects on GLAS-based regional estimates of forest biomass and carbon. Int. J. Remote Sens. 2010, 31, 1359–1372. [Google Scholar] [CrossRef] [Green Version]
  5. Feldpausch, T.R.; Lloyd, J.; Lewis, S.L.; Brienen, R.J.W.; Gloor, M.; Mendoza, A.M.; Lopezgonzalez, G.; Banin, L.; Salim, K.A.; Affumbaffoe, K. Tree height integrated into pantropical forest biomass estimates. Biogeosciences 2012, 9, 2567–2622. [Google Scholar] [CrossRef]
  6. Calders, K.; Newnham, G.; Burt, A.; Murphy, S.; Raumonen, P.; Herold, M.; Culvenor, D.; Avitabile, V.; Disney, M.; Armston, J. Nondestructive estimates of above-ground biomass using terrestrial laser scanning. Methods Ecol. Evol. 2015, 6, 198–208. [Google Scholar] [CrossRef]
  7. Janssen, T.A.J.; Ametsitsi, G.K.D.; Collins, M.; Adu-Bredu, S.; Oliveras, I.; Mitchard, E.T.A.; Veenendaal, E.M. Extending the baseline of tropical dry forest loss in Ghana (1984–2015) reveals drivers of major deforestation inside a protected area. Biol. Conserv. 2018, 218, 163–172. [Google Scholar] [CrossRef]
  8. Rahman, M.M.; Islam, M.S.; Pramanik, M. Monitoring of changes in woodlots outside forests by multi-temporal Landsat imagery. iForest Biogeosci. For. 2018, 11, 162–170. [Google Scholar] [CrossRef]
  9. Melnikova, I.; Awaya, Y.; Saitoh, T.; Muraoka, H.; Sasai, T. Estimation of Leaf Area Index in a Mountain Forest of Central Japan with a 30-m Spatial Resolution Based on Landsat Operational Land Imager Imagery: An Application of a Simple Model for Seasonal Monitoring. Remote Sens. 2018, 10, 179. [Google Scholar] [CrossRef]
  10. Zhang, Y.; Li, X.; Ling, F.; Atkinson, P.M.; Ge, Y.; Shi, L.; Du, Y. Updating Landsat-based forest cover maps with MODIS images using multiscale spectral-spatial-temporal superresolution mapping. Int. J. Appl. Earth Obs. 2017, 63, 129–142. [Google Scholar] [CrossRef]
  11. Fang, X.; Zhu, Q.; Ren, L.; Chen, H.; Wang, K.; Peng, C. Large-scale detection of vegetation dynamics and their potential drivers using MODIS images and BFAST: A case study in Quebec, Canada. Remote Sens. Environ. 2018, 206, 391–402. [Google Scholar] [CrossRef]
  12. Kahiu, M.N.; Hanan, N.P. Estimation of Woody and Herbaceous Leaf Area Index in Sub-Saharan Africa Using MODIS Data. J. Geophys. Res. Biogeosci. 2018, 123, 3–17. [Google Scholar] [CrossRef] [Green Version]
  13. Testa, S.; Soudani, K.; Boschetti, L.; Borgogno Mondino, E. MODIS-derived EVI, NDVI and WDRVI time series to estimate phenological metrics in French deciduous forests. Int. J. Appl. Earth Obs. 2018, 64, 132–144. [Google Scholar] [CrossRef]
  14. Wang, H.; Wang, C.; Wu, H. Using GF-2 Imagery and the Conditional Random Field Model for Urban Forest Cover Mapping. Remote Sens. Lett. 2016, 7, 378–387. [Google Scholar] [CrossRef]
  15. Guifen, S.; Xianlin, Q.; Lingyu, Y.; Shuchao, L.; Zeng-yuan, L.; Xiaozhong, C.; Xiangqing, Z. Changes Analysis of Post-Fire Vegetation Spectrum and Index Based on Time Series GF-1 WFV Images. Spectrosc. Spectr. Anal. 2018, 38, 511–517. [Google Scholar] [CrossRef]
  16. Zhou, T.; Zhao, M.; Sun, C.; Pan, J. Exploring the Impact of Seasonality on Urban Land-Cover Mapping Using Multi-Season Sentinel-1A and GF-1 WFV Images in a Subtropical Monsoon-Climate Region. ISPRS Int. J. Geo-Inf. 2018, 7, 3. [Google Scholar] [CrossRef]
  17. Zhao, J.; Li, J.; Liu, Q.; Wang, H.; Chen, C.; Xu, B.; Wu, S. Comparative Analysis of Chinese HJ-1 CCD, GF-1 WFV and ZY-3 MUX Sensor Data for Leaf Area Index Estimations for Maize. Remote Sens. 2018, 10, 68. [Google Scholar] [CrossRef]
  18. Xu, K.; Jiang, Y.; Zhang, G.; Zhang, Q.; Wang, X. Geometric Potential Assessment for ZY3-02 Triple Linear Array Imagery. Remote Sens. 2017, 9, 658. [Google Scholar] [CrossRef]
  19. Laurin, G.V.; Balling, J.; Corona, P.; Mattioli, W.; Papale, D.; Puletti, N.; Rizzo, M.; Truckenbrodt, J.; Urban, M. Above-ground biomass prediction by Sentinel-1 multitemporal data in central Italy with integration of ALOS2 and Sentinel-2 data. J. Appl. Remote Sens. 2018, 12. [Google Scholar] [CrossRef]
  20. Mura, M.; Bottalico, F.; Giannetti, F.; Bertani, R.; Giannini, R.; Mancini, M.; Orlandini, S.; Travaglini, D.; Chirici, G. Exploiting the capabilities of the Sentinel-2 multi spectral instrument for predicting growing stock volume in forest ecosystems. Int. J. Appl. Earth Obs. 2018, 66, 126–134. [Google Scholar] [CrossRef]
  21. Bye, I.J.; North, P.R.J.; Los, S.O.; Kljun, N.; Rosette, J.A.B.; Hopkinson, C.; Chasmer, L.; Mahoney, C. Estimating forest canopy parameters from satellite waveform LiDAR by inversion of the FLIGHT three-dimensional radiative transfer model. Remote Sens. Environ. 2017, 188, 177–189. [Google Scholar] [CrossRef]
  22. Sun, G.Q.; Ranson, K.J.; Zhang, Z.J. Forest Vertical Parameters from Lidar and Multi-Angle Imaging Spectrometer Data. J. Remote Sens. 2006, 10, 523–530. [Google Scholar]
  23. Sandberg, G.; Ulander, L.M.H.; Fransson, J.E.S.; Holmgren, J.; Toan, T.L. L- and P-band backscatter intensity for biomass retrieval in hemiboreal forest. Remote Sens. Environ. 2011, 115, 2874–2886. [Google Scholar] [CrossRef]
  24. Pang, Y.; Li, Z.; Ju, H.; Lu, H.; Jia, W.; Si, L.; Guo, Y.; Liu, Q.; Li, S.; Liu, L.; et al. LiCHy: The CAF’s LiDAR, CCD and Hyperspectral Integrated Airborne Observation System. Remote Sens. 2016, 8, 398. [Google Scholar] [CrossRef]
  25. Ip, A.; El-Sheimy, N.; Mostafa, M. Performance Analysis of Integrated Sensor Orientation. Photogramm. Eng. Remote Sens. 2007, 73, 89–97. [Google Scholar] [CrossRef] [Green Version]
  26. Kotchenova, S.Y. Modeling lidar waveforms with time-dependent stochastic radiative transfer theory for remote estimations of forest structure. J. Geophys. Res. 2003, 108. [Google Scholar] [CrossRef] [Green Version]
  27. Harding, D.J.; Lefsky, M.A.; Parker, G.G.; Blair, J.B. Laser altimeter canopy height profiles Methods and validation for closed-canopy, broadleaf forests. Remote Sens. Environ. 2001, 76, 283–297. [Google Scholar] [CrossRef]
  28. Lefsky, M.A.; Cohen, W.B.; Acker, S.A.; Parker, G.G.; Spies, T.A.; Harding, D. Lidar Remote Sensing of the Canopy Structure and Biophysical Properties of Douglas-Fir Western Hemlock Forests. Remote Sens. Environ. 1999, 70, 339–361. [Google Scholar] [CrossRef]
  29. Lefsky, M.A.; Harding, D.; Cohen, W.B.; Parker, G.; Shugart, H.H. Surface Lidar Remote Sensing of Basal Area and Biomass in Deciduous Forests of Eastern Maryland, USA. Remote Sens. Environ. 1999, 67, 83–98. [Google Scholar] [CrossRef] [Green Version]
  30. Blair, J.B.; Rabine, D.L.; Hofton, M.A. The Laser Vegetation Imaging Sensor: A medium-altitude, digitisation-only, airborne laser altimeter for mapping vegetation and topography. ISPRS J. Photogramm. Remote Sens. 1999, 54, 115–122. [Google Scholar] [CrossRef]
  31. Drake, J.B.; Dubayah, R.O.; Clark, D.B.; Knox, R.G.; Blair, J.B.; Hofton, M.A.; Chazdon, R.L.; Weishampel, J.F.; Prince, S. Estimation of tropical forest structural characteristics using large-footprint lidar. Remote Sens. Environ. 2002, 79, 305–319. [Google Scholar] [CrossRef]
  32. Harding, D.J. ICESat waveform measurements of within-footprint topographic relief and vegetation vertical structure. Geophys. Res. Lett. 2005, 32. [Google Scholar] [CrossRef] [Green Version]
  33. Zwally, H.J.; Schutz, B.; Abdalati, W.; Abshire, J.; Bentley, C.; Brenner, A.; Bufton, J.; Dezio, J.; Hancock, D.; Harding, D.; et al. ICESat’s laser measurements of polar ice, atmosphere, ocean, and land. J. Geodyn. 2002, 34, 405–445. [Google Scholar] [CrossRef]
  34. Duncanson, L.I.; Niemann, K.O.; Wulder, M.A. Estimating forest canopy height and terrain relief from GLAS waveform metrics. Remote Sens. Environ. 2010, 114, 138–154. [Google Scholar] [CrossRef]
  35. Enßle, F.; Heinzel, J.; Koch, B. Accuracy of vegetation height and terrain elevation derived from ICESat/GLAS in forested areas. Int. J. Appl. Earth Obs. 2014, 31, 37–44. [Google Scholar] [CrossRef]
  36. Luo, S.; Wang, C.; Li, G.; Xi, X. Retrieving leaf area index using ICESat/GLAS full-waveform data. Remote Sens. Lett. 2013, 4, 745–753. [Google Scholar] [CrossRef]
  37. Mahoney, C.; Hopkinson, C.; Kljun, N.; van Gorsel, E. Estimating Canopy Gap Fraction Using ICESat GLAS within Australian Forest Ecosystems. Remote Sens. 2017, 9, 59. [Google Scholar] [CrossRef]
  38. Tang, H.; Dubayah, R.; Brolly, M.; Ganguly, S.; Zhang, G. Large-scale retrieval of leaf area index and vertical foliage profile from the spaceborne waveform lidar (GLAS/ICESat). Remote Sens. Environ. 2014, 154, 8–18. [Google Scholar] [CrossRef]
  39. Wang, X.; Huang, H.; Gong, P.; Liu, C.; Li, C.; Li, W. Forest Canopy Height Extraction in Rugged Areas with ICESat/GLAS Data. IEEE Trans. Geosci. Sens. 2014, 52, 1650–4657. [Google Scholar] [CrossRef]
  40. Tian, J.; Wang, L.; Li, X.; Shi, C.; Gong, H. Differentiating Tree and Shrub LAI in a Mixed Forest with ICESat/GLAS Spaceborne LiDAR. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2017, 10, 87–94. [Google Scholar] [CrossRef]
  41. Pourrahmati, M.R.; Baghdadi, N.N.; Darvishsefat, A.A.; Namiranian, M.; Fayad, I.; Bailly, J.; Gond, V. Capability of GLAS/ICESat Data to Estimate Forest Canopy Height and Volume in Mountainous Forests of Iran. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 5246–5261. [Google Scholar] [CrossRef]
  42. Hancock, S.; Armston, J.; Hofton, M.; Sun, X.; Tang, H.; Duncanson, L.I.; Kellner, J.R.; Dubayah, R. The GEDI simulator: A large-footprint waveform lidar simulator for calibration and validation of spaceborne missions. Earth Space Sci. 2019, 6, 294–310. [Google Scholar] [CrossRef]
  43. Qi, W.; Lee, S.; Hancock, S.; Luthcke, S.; Tang, H.; Armston, J.; Dubayah, R. Improved forest height estimation by fusion of simulated GEDI Lidar data and TanDEM-X InSAR data. Remote Sens. Environ. 2019, 221, 621–634. [Google Scholar] [CrossRef]
  44. Silva, C.A.; Saatchi, S.; Garcia, M.; Labriere, N.; Klauberg, C.; Ferraz, A.; Meyer, V.; Jeffery, K.J.; Abernethy, K.; White, L.; et al. Comparison of Small- and Large-Footprint Lidar Characterization of Tropical Forest Aboveground Structure and Biomass: A Case Study from Central Gabon. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 3512–3526. [Google Scholar] [CrossRef]
  45. Saarela, S.; Holm, S.; Healey, S.; Andersen, H.; Petersson, H.; Prentius, W.; Patterson, P.; Næsset, E.; Gregoire, T.; Ståhl, G. Generalized Hierarchical Model-Based Estimation for Aboveground Biomass Assessment Using GEDI and Landsat Data. Remote Sens. 2018, 10, 1832. [Google Scholar] [CrossRef]
  46. Milenković, M.; Schnell, S.; Holmgren, J.; Ressl, C.; Lindberg, E.; Hollaus, M.; Pfeifer, N.; Olsson, H. Influence of footprint size and geolocation error on the precision of forest biomass estimates from space-borne waveform LiDAR. Remote Sens. Environ. 2017, 200, 74–88. [Google Scholar] [CrossRef]
  47. Popescu, S.C.; Wynne, R.H.; Nelson, R.F. Estimating plot-level tree heights with lidar: Local filtering with a canopy-height based variable window size. Comput. Electron. Agric. 2002, 37, 71–95. [Google Scholar] [CrossRef]
  48. Wang, M.; Sun, R.; Xiao, Z. Estimation of Forest Canopy Height and Aboveground Biomass from Spaceborne LiDAR and Landsat Imageries in Maryland. Remote Sens. 2018, 10, 344. [Google Scholar] [CrossRef]
  49. Yong, P.; Xinfang, Y.; Zengyuan, L.; Guoqing, S.; Erxue, C.; Bingxiang, T. Waveform Length Extraction from ICEsat GLAS Data and Forest Application Analysis. Sci. Silvae Sin. 2006, 42, 137–140. [Google Scholar]
  50. Sun, G.; Ranson, K.J.; Kimes, D.S.; Blair, J.B.; Kovacs, K. Forest vertical structure from GLAS: An evaluation using LVIS and SRTM data. Remote Sens. Environ. 2008, 112, 107–117. [Google Scholar] [CrossRef]
  51. Hilbert, C.; Schmullius, C. Influence of Surface Topography on ICESat/GLAS Forest Height Estimation and Waveform Shape. Remote Sens. 2012, 4, 2210–2235. [Google Scholar] [CrossRef] [Green Version]
  52. Lefsky, M.A.; Harding, D.J.; Keller, M.; Cohen, W.B.; Carabajal, C.C.; Del Bom Espirito-Santo, F.; Hunter, M.O.; de Oliveira, R. Estimates of forest canopy height and aboveground biomass using ICESat. Geophys. Res. Lett. 2005, 32. [Google Scholar] [CrossRef] [Green Version]
  53. Fayad, I.; Baghdadi, N.; Gond, V.; Bailly, J.S.; Barbier, N.; El Hajj, M.; Fabre, F. Coupling potential of ICESat/GLAS and SRTM for the discrimination of forest landscape types in French Guiana. Int. J. Appl. Earth Obs. 2014, 33, 21–31. [Google Scholar] [CrossRef] [Green Version]
  54. Park, T.; Kennedy, R.; Choi, S.; Wu, J.; Lefsky, M.; Bi, J.; Mantooth, J.; Myneni, R.; Knyazikhin, Y. Application of Physically-Based Slope Correction for Maximum Forest Canopy Height Estimation Using Waveform Lidar across Different Footprint Sizes and Locations: Tests on LVIS and GLAS. Remote Sens. 2014, 6, 6566–6586. [Google Scholar] [CrossRef] [Green Version]
  55. Wing, M.G.; Eklund, A. Performance Comparison of a Low-Cost Mapping Grade Global Positioning Systems (GPS) Receiver and Consumer Grade GPS Receiver under Dense Forest Canopy. J. For. 2007, 105, 9–14. [Google Scholar]
  56. Hernández-Stefanoni, J.; Reyes-Palomeque, G.; Castillo-Santiago, M.; George-Chacón, S.; Huechacona-Ruiz, A.; Tun-Dzul, F.; Rondon-Rivera, D.; Dupuy, J. Effects of Sample Plot Size and GPS Location Errors on Aboveground Biomass Estimates from LiDAR in Tropical Dry Forests. Remote Sens. 2018, 10, 1586. [Google Scholar] [CrossRef]
  57. Abdia, E.; Mariva, H.S.; Deljoueia, A.; Sohrabib, H. Accuracy and precision of consumer-grade GPS positioning in an urban green space environment. Forest Sci. Technol. 2014, 10, 141–147. [Google Scholar] [CrossRef] [Green Version]
  58. Rosette, J.A.B.; North, P.R.J.; Suárez, J.C.; Los, S.O. Uncertainty within satellite LiDAR estimations of vegetation and topography. Int. J. Remote Sens. 2010, 31, 1325–1342. [Google Scholar] [CrossRef]
  59. Keller, M. Revised method for forest canopy height estimation from Geoscience Laser Altimeter System waveforms. J. Appl. Remote Sens. 2007, 1, 13537–13555. [Google Scholar] [CrossRef]
Figure 1. (a) Two-dimensional stabilized platform; (b) laser, CCD, IMU, and APD are installed in (a); (c) signal processing and control unit; (a,c) are mounted on the airplane adjacent to one another (d).
Figure 1. (a) Two-dimensional stabilized platform; (b) laser, CCD, IMU, and APD are installed in (a); (c) signal processing and control unit; (a,c) are mounted on the airplane adjacent to one another (d).
Sensors 19 01699 g001
Figure 2. Laser sensor unit design.
Figure 2. Laser sensor unit design.
Sensors 19 01699 g002
Figure 3. Signal processing and control unit structure.
Figure 3. Signal processing and control unit structure.
Sensors 19 01699 g003
Figure 4. Laboratory test of laser emission pulse.
Figure 4. Laboratory test of laser emission pulse.
Sensors 19 01699 g004
Figure 5. Laboratory test of laser echo signal.
Figure 5. Laboratory test of laser echo signal.
Sensors 19 01699 g005
Figure 6. Time jitter result.
Figure 6. Time jitter result.
Sensors 19 01699 g006
Figure 7. Corresponding distance measurement error.
Figure 7. Corresponding distance measurement error.
Sensors 19 01699 g007
Figure 8. A segment of the flightline on which large-footprint LiDAR data (red points), small-footprint LiDAR (large rectangular area) and ground plots (yellow points) were collected.
Figure 8. A segment of the flightline on which large-footprint LiDAR data (red points), small-footprint LiDAR (large rectangular area) and ground plots (yellow points) were collected.
Sensors 19 01699 g008
Figure 9. Estimation of mean forest height. The black line connects each of the 500 ns-level reflectance observations, and the red line is a filtered (fast Fourier transformation) representation of the same information.
Figure 9. Estimation of mean forest height. The black line connects each of the 500 ns-level reflectance observations, and the red line is a filtered (fast Fourier transformation) representation of the same information.
Sensors 19 01699 g009
Figure 10. Overview of height modeling and validation process for the estimation of mean forest height.
Figure 10. Overview of height modeling and validation process for the estimation of mean forest height.
Sensors 19 01699 g010
Figure 11. Relationship between 70th percentile small-footprint LiDAR height values and heights measured on sample plots (n = 214). The dark line is the linear regression line from Equation (12).
Figure 11. Relationship between 70th percentile small-footprint LiDAR height values and heights measured on sample plots (n = 214). The dark line is the linear regression line from Equation (12).
Sensors 19 01699 g011
Figure 12. Mean forest height comparison between estimated by Equation (12) and estimated by large-footprint LiDAR data.
Figure 12. Mean forest height comparison between estimated by Equation (12) and estimated by large-footprint LiDAR data.
Sensors 19 01699 g012
Table 1. iXU-R180 technical index.
Table 1. iXU-R180 technical index.
Index NameDesign Value
Resolution10,328 × 7760 (80 MP)
Dynamic range>72 db
Pixel size5.2 µm
CCD size effective53.7 mm × 40.4 mm
Aspect ratio4:3
Light sensitivity (ISO)35–800
Shutter speed1/1600 s
Camera lens focus50 mm f/4.0
FOV56.5° × 44°
InterfacesUSB 3.0
Minimum photo interval1.8 s
Data storage1 TB SSD storage (optional iX Controller)
CompactFlash card Type I/II including UDMA 6 and 7
Table 2. Power and space of laser emission pulse.
Table 2. Power and space of laser emission pulse.
NameValueUnit
Power/energy
Total energy445,221,505.16cnts
Peak value8777.03cnts
Minimum value−205.77cnts
Space
The coordinate of centroid on the X-axis1.484627 × 104μm
The coordinate of centroid on the Y-axis1.230543 × 104μm
D4 σ X2.744 × 104μm
D4 σ Y2.019 × 104μm

Share and Cite

MDPI and ACS Style

Hu, Y.; Wu, F.; Sun, Z.; Lister, A.; Gao, X.; Li, W.; Peng, D. The Laser Vegetation Detecting Sensor: A Full Waveform, Large-Footprint, Airborne Laser Altimeter for Monitoring Forest Resources. Sensors 2019, 19, 1699. https://doi.org/10.3390/s19071699

AMA Style

Hu Y, Wu F, Sun Z, Lister A, Gao X, Li W, Peng D. The Laser Vegetation Detecting Sensor: A Full Waveform, Large-Footprint, Airborne Laser Altimeter for Monitoring Forest Resources. Sensors. 2019; 19(7):1699. https://doi.org/10.3390/s19071699

Chicago/Turabian Style

Hu, Yang, Fayun Wu, Zhongqiu Sun, Andrew Lister, Xianlian Gao, Weitao Li, and Daoli Peng. 2019. "The Laser Vegetation Detecting Sensor: A Full Waveform, Large-Footprint, Airborne Laser Altimeter for Monitoring Forest Resources" Sensors 19, no. 7: 1699. https://doi.org/10.3390/s19071699

APA Style

Hu, Y., Wu, F., Sun, Z., Lister, A., Gao, X., Li, W., & Peng, D. (2019). The Laser Vegetation Detecting Sensor: A Full Waveform, Large-Footprint, Airborne Laser Altimeter for Monitoring Forest Resources. Sensors, 19(7), 1699. https://doi.org/10.3390/s19071699

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop