[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (16)

Search Parameters:
Keywords = DTM quality estimation

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
23 pages, 3263 KiB  
Article
Application of Diverse Testing to Improve Integrated Circuit Test Yield and Quality
by Chung-Huang Yeh, Shou-Rong Chen and Kan-Hsiang Liao
Eng 2024, 5(4), 3517-3539; https://doi.org/10.3390/eng5040183 - 20 Dec 2024
Viewed by 428
Abstract
This paper utilizes the digital integrated circuit testing model to compute the test yield curve of future wafers and explore the influence of test guardband (TGB) on quality and yield. With the passage of three years since the COVID-19 pandemic disrupted semiconductor production [...] Read more.
This paper utilizes the digital integrated circuit testing model to compute the test yield curve of future wafers and explore the influence of test guardband (TGB) on quality and yield. With the passage of three years since the COVID-19 pandemic disrupted semiconductor production lines, the semiconductor manufacturing industry still faces chip shortages. Although initiatives such as the CHIPS and Science Act in the United States have helped stabilize chip supply chains, manufacturers still face inventory shortages and delayed deliveries. Moreover, the backwardness and inaccuracy of semiconductor test equipment have led to a decline in both test yield and wafer quality, resulting in reduced shipments. Therefore, to mitigate yield losses and enhance the test yield and shipment volume of semiconductor products, this paper proposes a diverse test method (DTM) to improve test outcomes through the alteration of the testing strategy and TGB adjustment. Furthermore, according to the wafer estimation table published in the IEEE International Roadmap for Devices and Systems (2023), the proposed DTM can effectively enhance the test yield of wafers and improve the testing capabilities of ATE testers (automatic test equipment). Consequently, suppliers can stabilize the chip supply chain and enhance their companies’ profits and reputation by improving chip test yield. Full article
(This article belongs to the Section Electrical and Electronic Engineering)
Show Figures

Figure 1

Figure 1
<p>Errors in semiconductor manufacturing and testing processes.</p>
Full article ">Figure 2
<p>Distribution of chip delay times after wafer fabrication.</p>
Full article ">Figure 3
<p>Threshold test module.</p>
Full article ">Figure 4
<p>Impact of test guardband selection on test yield and test quality.</p>
Full article ">Figure 5
<p>Influence of test guardband on wafer yield and test quality during testing.</p>
Full article ">Figure 6
<p>Normal distribution of chip delay time after manufacturing and normal distribution of integrated circuit (IC) tester capability (OTA).</p>
Full article ">Figure 7
<p>Influence of automated test equipment (ATE) tester accuracy parameter (OTA) on test yield and quality.</p>
Full article ">Figure 8
<p>Test decision chart and decision-making methodology for diverse testing methods.</p>
Full article ">Figure 9
<p>Comparison of test results between the diverse test method (DTM) and traditional testing method (TTM).</p>
Full article ">Figure 10
<p>Improving test outcomes (300 ppm) with the diverse test method according to the IRDS 2023 table.</p>
Full article ">Figure 11
<p>Enhancement of the yield of high-quality product testing (10 ppm) through the diverse testing method.</p>
Full article ">
19 pages, 5704 KiB  
Article
Error Analysis and Accuracy Improvement in Forest Canopy Height Estimation Based on GEDI L2A Product: A Case Study in the United States
by Yi Li, Shijuan Gao, Haiqiang Fu, Jianjun Zhu, Qing Hu, Dong Zeng and Yonghui Wei
Forests 2024, 15(9), 1536; https://doi.org/10.3390/f15091536 - 31 Aug 2024
Cited by 3 | Viewed by 1241
Abstract
Various error factors influence the inversion of forest canopy height using GEDI full-waveform LiDAR data, and the interaction of these factors impacts the accuracy of forest canopy height estimation. From an error perspective, there is still a lack of methods to fully correct [...] Read more.
Various error factors influence the inversion of forest canopy height using GEDI full-waveform LiDAR data, and the interaction of these factors impacts the accuracy of forest canopy height estimation. From an error perspective, there is still a lack of methods to fully correct the impact of various error factors on the retrieval of forest canopy height from GEDI. From the modeling perspective, establishing clear coupling models between various environments, collection parameters, and GEDI forest canopy height errors is challenging. Understanding the comprehensive impact of various environments and collection parameters on the accuracy of GEDI data is crucial for extracting high-quality and precise forest canopy heights. First, we quantitatively assessed the accuracy of GEDI L2A data in forest canopy height inversion and conducted an error analysis. A GEDI forest canopy height error correction model has been developed, taking into account both forest density and terrain effects. This study elucidated the influence of forest density and terrain on the error in forest canopy height estimation, ultimately leading to an improvement in the accuracy of forest canopy height inversion. In light of the identified error patterns, quality control criteria for GEDI footprints are formulated, and a correction model for GEDI forest canopy height is established to achieve high-precision inversion. We selected 19 forest areas located in the United States with high-accuracy Digital Terrain Models (DTMs) and Canopy Height Models (CHMs) to analyze the error factors of GEDI forest canopy heights and assess the proposed accuracy improvement for GEDI forest canopy heights. The findings reveal a decrease in the corrected RMSE value of forest canopy height from 5.60 m to 4.19 m, indicating a 25.18% improvement in accuracy. Full article
Show Figures

Figure 1

Figure 1
<p>Overview of the study area and distribution of each site.</p>
Full article ">Figure 2
<p>The flowchart of the pre-processing.</p>
Full article ">Figure 3
<p>Correlation coefficients between GEDI forest canopy height and the reference value.</p>
Full article ">Figure 4
<p>Scatter plot and difference statistics of forest canopy height estimation of GEDI. (<b>a</b>) Scatter plot between GEDI forest heights and NEON forest heights. (<b>b</b>) The frequency of difference statistics of forest canopy height estimation of GEDI.</p>
Full article ">Figure 5
<p>Performance statistics of forest canopy height for each site. (<b>a</b>) Numbers. (<b>b</b>) Bais. (<b>c</b>) RMSE. (<b>d</b>) R2. (<b>e</b>) %RMSE. (<b>f</b>) MAE.</p>
Full article ">Figure 6
<p>Effects of beam type and acquisition time on the accuracy of GEDI forest canopy height.</p>
Full article ">Figure 7
<p>The impact of error factors on GEDI forest canopy height estimation. (<b>a</b>) number of peaks; (<b>b</b>) sensitivity; (<b>c</b>) canopy cover; (<b>d</b>) slope.</p>
Full article ">Figure 8
<p>Importance of error factors.</p>
Full article ">Figure 9
<p>Comparison of scatter plots before and after filtering. (<b>a</b>) before; (<b>b</b>) after.</p>
Full article ">Figure 10
<p>Flowchart of forest canopy height correction.</p>
Full article ">Figure 11
<p>Variation in model error with the number of decision trees.</p>
Full article ">Figure 12
<p>Comparison of validation results before and after correction. (<b>a</b>) before; (<b>b</b>) after.</p>
Full article ">Figure 13
<p>Accuracy comparison before and after correction for each NEON site. (<b>a</b>) RMSE; (<b>b</b>) R<sup>2</sup>.</p>
Full article ">
21 pages, 20260 KiB  
Article
Assessment of Leica CityMapper-2 LiDAR Data within Milan’s Digital Twin Project
by Marica Franzini, Vittorio Marco Casella and Bruno Monti
Remote Sens. 2023, 15(21), 5263; https://doi.org/10.3390/rs15215263 - 6 Nov 2023
Cited by 1 | Viewed by 2122
Abstract
The digital twin is one of the most promising technologies for realizing smart cities in terms of planning and management. For this purpose, Milan, Italy, has started a project to acquire aerial nadir and oblique images and LiDAR and terrestrial mobile mapping data. [...] Read more.
The digital twin is one of the most promising technologies for realizing smart cities in terms of planning and management. For this purpose, Milan, Italy, has started a project to acquire aerial nadir and oblique images and LiDAR and terrestrial mobile mapping data. The Leica CityMapper-2 hybrid sensor has been used for aerial surveys as it can capture precise and high-resolution multiple data (imagery and LiDAR). The surveying activities are completed, and quality checks are in progress. This paper concerns assessing aerial LiDAR data of a significant part of the metropolitan area, particularly evaluating the accuracy, precision, and congruency between strips and the point density estimation. The analysis has been conducted by exploiting a ground control network of GNSS and terrestrial LiDAR measurements created explicitly for this purpose. The vertical component has an accuracy root mean square error (RMSE) of around 5 cm, and a horizontal component of around 12 cm. Meanwhile, the precision RMSE ranges from 2 to 8 cm. These values are suitable for generating products such as DSM/DTM. Full article
(This article belongs to the Special Issue Lidar Sensing for 3D Digital Twins)
Show Figures

Figure 1

Figure 1
<p>The metropolitan area of Milan covers an area of 1776 square kilometers; Zone1 and Zone2 are colored in green and blue, respectively. In the accompanying frames, the location of the site within the Lombardy region and Italy is displayed.</p>
Full article ">Figure 2
<p>Oblique LiDAR scanning pattern (Bacher, 2022 [<a href="#B19-remotesensing-15-05263" class="html-bibr">19</a>]).</p>
Full article ">Figure 3
<p>The RGB representations of the acquired point clouds.</p>
Full article ">Figure 4
<p>The CIR representations of the acquired point clouds.</p>
Full article ">Figure 5
<p>Location of the 200 ARCOs; red dots represent the Special-ARCOs (polygons color have the same meaning as <a href="#remotesensing-15-05263-f001" class="html-fig">Figure 1</a>). The small frame shows an example of the relative position of the two ARCO benchmarks (the example reports the same ARCO shown in <a href="#remotesensing-15-05263-f006" class="html-fig">Figure 6</a>).</p>
Full article ">Figure 6
<p>An example of ARCOs: (<b>a</b>) TypeA benchmark stuck in a large concrete curb; (<b>b</b>) TypeB photogrammetric marker positioned in a flat area, without obstacles or slope variations in a neighborhood of 2 m.</p>
Full article ">Figure 7
<p>Examples of Special-ARCOs: (<b>a</b>) ARCO #7; (<b>b</b>) ARCO #49.</p>
Full article ">Figure 8
<p>Overlapping between adjacent strips: (<b>a</b>) histogram of overlap; (<b>b</b>) example of overlapping between two point cloud tiles belonging to adjacent strips.</p>
Full article ">Figure 9
<p>Location of ARCO_50B within the LiDAR strip; the small frame shows the uniform acquisition pattern near the strip’s central line.</p>
Full article ">Figure 10
<p>Location of ARCO_93B within the LiDAR strip; the small frame shows the effect of the circular acquisition pattern near the strip’s border.</p>
Full article ">Figure 11
<p>Density trend according to distance from strip’s central line; blue asterisks depict the mean density in each considered ARCO.</p>
Full article ">Figure 12
<p>Heuristic explanation of vertical error estimation.</p>
Full article ">Figure 13
<p>Histogram of vertical distances between ARCOs-typeB vertices and LiDAR clouds.</p>
Full article ">Figure 14
<p>Examples of the extracted flat surfaces for vertical error analysis: (<b>a</b>) a portion of a road close to the scan location; (<b>b</b>) a portion of a rooftop. ALS and TLS data are reported in green and magenta, respectively.</p>
Full article ">Figure 15
<p>Histogram of vertical distances between Special-ARCOs’ selected horizontal surfaces and LiDAR clouds.</p>
Full article ">Figure 16
<p>Examples of the extracted flat surfaces for horizontal error analysis. ALS and TLS data are reported in green and magenta, respectively. (<b>a</b>,<b>b</b>) Two examples of building facades used for planimetric accuracy evaluation.</p>
Full article ">Figure 17
<p>Histogram of horizontal distances between Special-ARCOs’ selected vertical surfaces and LiDAR clouds.</p>
Full article ">Figure 18
<p>Examples of road marking visual comparisons for the assessment of horizontal accuracy. TLS datasets are colored in accordance with their native RGB information, while the ALS-selected high-reflectivity points are depicted with orange dots. (<b>a</b>–<b>c</b>) Some examples of zebra crossing used for planimetric accuracy visual evaluation.</p>
Full article ">
16 pages, 21054 KiB  
Article
Resolving Ambiguities in SHARAD Data Analysis Using High-Resolution Digital Terrain Models
by Léopold Desage, Alain Herique, Sylvain Douté, Sonia Zine and Wlodek Kofman
Remote Sens. 2023, 15(3), 764; https://doi.org/10.3390/rs15030764 - 28 Jan 2023
Cited by 1 | Viewed by 2195
Abstract
The SHAllow RADar (SHARAD) onboard Mars Reconnaissance Orbiter (MRO) is a 20 MHz Synthetic Aperture Radar (SAR) that probes the first hundreds of meters of the Martian subsurface. In order to interpret the detection of subsurface interfaces with ground penetrating radars, simulations using [...] Read more.
The SHAllow RADar (SHARAD) onboard Mars Reconnaissance Orbiter (MRO) is a 20 MHz Synthetic Aperture Radar (SAR) that probes the first hundreds of meters of the Martian subsurface. In order to interpret the detection of subsurface interfaces with ground penetrating radars, simulations using Digital Terrain Models (DTM) are necessary. This methodology paper focuses on the analysis of the first tens of meters of the Martian subsurface with SHARAD, comparing the use of different high-resolution DTMs for radar simulation, namely, from the High-Resolution Stereo Camera (HRSC) onboard the Mars Express and from the Context Camera (CTX) onboard MRO. The region of Terra Cimmeria was chosen as a demonstration area. It is a highly cratered southern midlatitude region, where, as will be discussed, the higher resolution of the aforementioned terrain models is mandatory to describe the surface at an acceptable level of detail for shallow subsurface radar interpretation. With a DTM corrected by photoclinometry using CTX imagery, we show that a reflector that was visible on SHARAD data but not on the simulation made with an HRSC DTM is, in fact, a surface echo that was not reproduced by the HRSC surface model. We also show that, unlike laser altimetry DTMs, optical DTMs are prone to artifacts that can make radar analysis more complicated for some scenarios. Reciprocally, we show that the comparison between radar and its corresponding simulated data is a way of assessing a DTM’s quality, which is especially useful in missions where ground control points are lacking, unlike Martian observations. Full article
(This article belongs to the Special Issue Radar for Planetary Exploration)
Show Figures

Figure 1

Figure 1
<p>Single echo of SHARAD data compared to an echo of the simulation. The largest peaks of SHARAD data are well matched to the simulation.</p>
Full article ">Figure 2
<p>(<b>a</b>) Output of the simulation of a 150 km portion of the SHARAD dataset n°5128501 in Terra Cimmeria using a MOLA DTM. (<b>b</b>) Same profile after SAR synthesis. Parabolas in the simulation are the result of the trajectory of the spacecraft relative to each reflector.</p>
Full article ">Figure 3
<p>(<b>a</b>) A 40 km long portion of the SHARAD profile n°0806502 sounding Deuteronilus Mensae; (<b>b</b>) corresponding simulation using a MOLA DTM; (<b>c</b>) HRSC image of the terrain at nadir in the area where the radargram was taken. The two reflectors circled in white are present in the SHARAD data, but not in the simulation, meaning that they probably originate from the subsurface. Note that the two top echoes (left and right of the radargram) are off-nadir echoes originating from plateaus, thus arriving before the stronger nadir echo.</p>
Full article ">Figure 4
<p>(<b>a</b>) A 40 km portion of the SHARAD radargram n°01708401 sounding Hellas Planitia; (<b>b</b>) corresponding simulation using a MOLA DTM. (<b>c</b>) CTX image of the terrain at the nadir on the area where the radargram was acquired. In areas where the roughness is higher, we can see that the identification of subsurface echoes is much harder with MOLA, as the faintest reflectors are not reproduced by the simulation.</p>
Full article ">Figure 5
<p>(<b>a</b>) A 65 km portion of the SHARAD profile n°5128501; (<b>b</b>) corresponding simulation using a MOLA DTM. The reflector circled in white is present in the SHARAD data but not in the MOLA simulation.</p>
Full article ">Figure 6
<p>(<b>a</b>) MOLA DTM on Tarq crater at a 463 m per pixel resolution. (<b>b</b>) The corresponding HRSC DTM at 50 m per pixel.</p>
Full article ">Figure 7
<p>(<b>a</b>) A 12 km portion of the SHARAD dataset n°5128501 on Terra Cimmeria. (<b>b</b>) Simulation using a HRSC DTM. (<b>c</b>) Simulation using a MOLA DTM. Finer details are reproduced by HRSC, but a higher level of artifacts is present. The reflector circled in white is not visible in either of the simulations presented on the two rightmost images.</p>
Full article ">Figure 8
<p>Visual comparison of the reconstruction of the edge of Tarq crater (southern midlatitudes) with a HRSC DTM at 50 m per pixel (<b>left</b>) and a CTX DTM at 12 m per pixel (<b>right</b>).</p>
Full article ">Figure 9
<p>Comparison between a 40 km portion of the SHARAD dataset n°5128501 and the simulation of it using different DTMs. The results using the downsampled CTX DTM (<b>d</b>) at 100 m per pixel represents the echoes more accurately than the simulation using the HRSC DTM at 50 m per pixel (finer details visible and lower level of artifacts). Parabolas on the right images highlight the echoes that are most significatively improved by the CTX simulation, as the left images are there to visualize the echoes without figures on top.</p>
Full article ">Figure 10
<p>Mapping the identified echo on a HRSC DTM with a CTX image mapped onto it. The area of origin of the echo is materialized by the white dotted area in the center. If the echo came from the surface, it would come from the edge of a plateau.</p>
Full article ">Figure 11
<p>Comparison between vertical exaggeration of topographies (exaggeration factor of 20) for the HRSC DTM with a CTX image mapped on it (<b>a</b>) and the HCPC DTM (<b>b</b>). The edge of the circled plateau was straightened by photoclinometry.</p>
Full article ">Figure 12
<p>Data-processing flowchart for photoclinometry using the HRSC DTM and the CTX image.</p>
Full article ">Figure 13
<p>(<b>a</b>,<b>b</b>) Comparison of the angles between the facets and the spacecraft for the original HRSC DTM (<b>a</b>) and for the HCPC DTM (<b>b</b>). (<b>c</b>) Outlines of the two areas used for statistical comparison on the northern (cyan) and southern (yellow) plateaus. The southern plateau area geometry was chosen to match the calculations of the potential origin of the reflector (<a href="#remotesensing-15-00764-f010" class="html-fig">Figure 10</a>). Both areas are about 1.5 km<sup>2</sup>. (<b>d</b>) CTX image for context. The angles are shallower in the refined DTM in the area where the reflector is thought to come from (circled in green).</p>
Full article ">Figure 14
<p>Comparison between a 18 km portion of the SHARAD dataset n°5128501 (<b>a</b>) and simulations using the corrected DTM with photoclinometry. The simulation in (<b>b</b>) was performed with the original DTM refined by photoclinometry at 6m. We do not see the surface on the simulated radargrams due to the DTM not covering the nadir part of the trajectory. (<b>c</b>) Wavelet-transformed DTM filtered at 320 m. (<b>d</b>) RGB composition with the SHARAD radargram on the red channel, and the simulation using the filtered DTM refined by photoclinometry on the cyan channel. The reflector that we are looking for is circled in white.</p>
Full article ">Figure 15
<p>(<b>a</b>–<b>c</b>) Mapping of the across-track round-trip distance between the points of the DTMs and the spacecraft. The slant range was mapped modulo lambda using a cyclic color scale that can be directly linked to the evolution of the signal’s phase across track. The plateau that we want to correct by photoclinometry is located in the white circle. The areas where the phase is stationary are the areas where the radar signal reflects on, and the larger the area, the brighter the echo. (<b>a</b>) Range map in the MOLA DTM. (<b>b</b>) Range map in the 6m HCPC DTM. (<b>c</b>) Range map of the 320 m scale wavelet-transform of the HCPC DTM. (<b>d</b>) HCPC-shaded DTM to visualize the terrain details.</p>
Full article ">
21 pages, 7985 KiB  
Article
Assessing the Quality of the Leica BLK2GO Mobile Laser Scanner versus the Focus 3D S120 Static Terrestrial Laser Scanner for a Preliminary Study of Garden Digital Surveying
by Graziella Del Duca and Carol Machado
Heritage 2023, 6(2), 1007-1027; https://doi.org/10.3390/heritage6020057 - 25 Jan 2023
Cited by 15 | Viewed by 3860
Abstract
Gardens play a key role in the definition of the cultural landscape since they reflect the culture, identity, and history of a people. They also contribute to the ecological balance of the city. Despite the fact that gardens have an historic and social [...] Read more.
Gardens play a key role in the definition of the cultural landscape since they reflect the culture, identity, and history of a people. They also contribute to the ecological balance of the city. Despite the fact that gardens have an historic and social value, they are not protected as much as the rest of the existing heritage, such as architecture and archaeological sites. While methods of built-heritage mapping and monitoring are increasing and constantly improving to reduce built-heritage loss and the severe impact of natural disasters, the documentation and survey techniques for gardens are often antiquated. In addition, inventories are typically made by non-updated/updateable reports, and they are rarely in digital format or in 3D. This paper presents the results of a comprehensive study on the latest technology for laser scanning in gardens. We compared static terrestrial laser scanning and mobile laser scanning point clouds generated by the Focus 3D S120 and the Leica BLK2GO, respectively, to evaluate their quality for documentation, estimate tree attributes, and terrain morphology. The evaluation is based on visual observation, C2C comparisons, and terrain information extraction capabilities, i.e., M3C2 comparisons for topography, DTM generation, and contour lines. Both methods produced useful outcomes for the scope of the research within their limitations. Terrestrial laser scanning is still the method that offers accurate point clouds with a higher point density and less noise. However, the more recent mobile laser scanning is able to survey in less time, significantly reducing the costs for site activities, data post-production, and registration. Both methods have their own restrictions that are amplified by site features, mainly the lack of plans for the geometric alignment of scans and the simultaneous location and mapping (SLAM) process. We offer a critical description of the issues related to the functionality of the two sensors, such as the operative range limit, light dependency, scanning time, point cloud completeness and size, and noise level. Full article
Show Figures

Figure 1

Figure 1
<p>Site localization in red. The studied area is located in Monsanto Park, near the Faculty of Architecture in Lisbon.</p>
Full article ">Figure 2
<p>Top view of the surveyed areas with TLS (<b>a</b>) and MLS (<b>b</b>).</p>
Full article ">Figure 3
<p>Point clouds of the scope area obtained with the TLS (<b>a</b>) and MLS (<b>b</b>). Difference in density, detail accuracy, and level of noise is already noticeable from this top view.</p>
Full article ">Figure 4
<p>Comparison between the two point clouds. (<b>a</b>) The two point clouds are overlayed in a horizontal slice having a height of 10 cm: in red, that of Faro Focus, and in green, that of BLK2GO. (<b>b</b>) Histogram of C2C comparison. (<b>c</b>) C2C analysis.</p>
Full article ">Figure 5
<p>Point clouds of the scope area obtained with the TLS (<b>a</b>) and MLS (<b>b</b>), respectively.</p>
Full article ">Figure 6
<p>Maximum crown height of the trees. Tree crown profile is complete in both MLS (<b>a</b>) and TLS (<b>b</b>) scans. Sections S1, S2, and S3 display point clouds as captured by the two scanners and they are visualized by z-scalar field. Section box extensions in (<b>c</b>). Point distribution diagrams in (<b>d</b>,<b>e</b>); the highest point for Faro dataset is at 13.73 m, and it is at 13.94 m for BLK2GO dataset.</p>
Full article ">Figure 7
<p>Terrain data comparison. (<b>a</b>) Cloud-to-mesh comparison between the reference data (Faro Focus point cloud) and the Digital Terrain Model generated from the segmented terrain cloud acquired with the Leica BLK2GO scanner. (<b>b</b>) Histogram describing the deviation distribution. (<b>c</b>) Cloud-to-cloud comparison; Faro Focus 3D point cloud is set as reference data. (<b>d</b>) M3C2 cloud-to-cloud Z-distance analysis.</p>
Full article ">Figure 8
<p>Contour lines every 50 cm are obtained from the Faro Focus’s point cloud (<b>a</b>) and BLK2GO’s (<b>b</b>).</p>
Full article ">Figure 9
<p>Horizontal sections of trunks obtained at a height of 2 m from the ground.</p>
Full article ">Figure 10
<p>Horizontal sections of the trees on the periphery area obtained by TLS and MLS. Although there is more noise, MLS is able to better define the sections in less scanning time.</p>
Full article ">Figure 11
<p>5cm tick point cloud slices, respectively, at 1 m (blue), 2 m (dark green), 3 m (light green), 4 m (yellow), and 5 m (red) height from the terrain. (<b>a</b>) 3D visualization of the sections. (<b>b</b>) Two-dimensional plan view of the trunk sections; those included in the perimeter area are the sections where the point cloud describes the entire trunk perimeter at each level, for both Faro and BLK2GO data.</p>
Full article ">Figure 12
<p>Perspective view of the total surveyed area: TLS (<b>a</b>) and MLS (<b>b</b>).</p>
Full article ">
40 pages, 74193 KiB  
Article
Revealing Active Mars with HiRISE Digital Terrain Models
by Sarah S. Sutton, Matthew Chojnacki, Alfred S. McEwen, Randolph L. Kirk, Colin M. Dundas, Ethan I. Schaefer, Susan J. Conway, Serina Diniega, Ganna Portyankina, Margaret E. Landis, Nicole F. Baugh, Rodney Heyd, Shane Byrne, Livio L. Tornabene, Lujendra Ojha and Christopher W. Hamilton
Remote Sens. 2022, 14(10), 2403; https://doi.org/10.3390/rs14102403 - 17 May 2022
Cited by 24 | Viewed by 5366
Abstract
Many discoveries of active surface processes on Mars have been made due to the availability of repeat high-resolution images from the High Resolution Imaging Science Experiment (HiRISE) onboard the Mars Reconnaissance Orbiter. HiRISE stereo images are used to make digital terrain models (DTMs) [...] Read more.
Many discoveries of active surface processes on Mars have been made due to the availability of repeat high-resolution images from the High Resolution Imaging Science Experiment (HiRISE) onboard the Mars Reconnaissance Orbiter. HiRISE stereo images are used to make digital terrain models (DTMs) and orthorectified images (orthoimages). HiRISE DTMs and orthoimage time series have been crucial for advancing the study of active processes such as recurring slope lineae, dune migration, gully activity, and polar processes. We describe the process of making HiRISE DTMs, orthoimage time series, DTM mosaics, and the difference of DTMs, specifically using the ISIS/SOCET Set workflow. HiRISE DTMs are produced at a 1 and 2 m ground sample distance, with a corresponding estimated vertical precision of tens of cm and ∼1 m, respectively. To date, more than 6000 stereo pairs have been acquired by HiRISE and, of these, more than 800 DTMs and 2700 orthoimages have been produced and made available to the public via the Planetary Data System. The intended audiences of this paper are producers, as well as users, of HiRISE DTMs and orthoimages. We discuss the factors that determine the effective resolution, as well as the quality, precision, and accuracy of HiRISE DTMs, and provide examples of their use in time series analyses of active surface processes on Mars. Full article
Show Figures

Figure 1

Figure 1
<p>HiRISE focal plane layout schematic projected onto the ground for a nadir imaging scenario (not to scale). Inset illustrates stereo imaging, where MRO rolls perpendicular to the along-track direction to point HiRISE at the target on different orbits. Emission angles are <math display="inline"><semantics> <msub> <mi>θ</mi> <mn>1</mn> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>θ</mi> <mn>2</mn> </msub> </semantics></math>, shown here as opposite side roll angles.</p>
Full article ">Figure 2
<p>Simplified flowchart of DTM and orthoimage workflow. Numbers correspond to subsections in <a href="#sec2dot1-remotesensing-14-02403" class="html-sec">Section 2.1</a>. There are circumstances when this workflow would be modified, or iterated on more than once or twice, but the decision path and steps are essentially those followed at the HiRISE Operations Center (HiROC).</p>
Full article ">Figure 3
<p>Plots of pixel offsets in the sample, or cross-track (<b>top</b>), and line, or along–track (<b>bottom</b>) directions, measured in the overlapping region of the RED4 and RED5 CCDs over the length of observation ESP_049009_1520. Before jitter correction (purple), pixel offsets show an amplitude of up to 4 px. After correction (light blue), offset amplitudes are reduced to ∼0.5 px in both the sample and line directions.</p>
Full article ">Figure 4
<p>Examples of typical artifacts found in HiRISE DTMs, shown in terrain shaded relief maps, which emphasize the surface texture of the model. In all panels illumination is from upper left, with north oriented to the top of the page. Arrows indicate artifacts. (<b>a</b>) Tiling “boxes” typical of NGATE, and linear artifacts along CCD seams (DTEPC_009689_2645_010084_2645). (<b>b</b>) Isolated noisy boxes possibly due to jitter (DTEEC_057321_1220_057453_1220_A01). (<b>c</b>) Triangular facets in an area of deep shadows in the stereo images (DTEPC_041121_0985_041029_0985_A01). (<b>d</b>) Areas of failed stereo correlation over dark dunes with ripples that were active between the acquisition of the two stereo images (DTEEC_050438_1890_051071_1890_A01). (<b>e</b>) Large triangular facets due to deeply shadowed areas with differing shadow boundaries, as well as an obscured surface in one image due to viewing geometry (DTEEC_069071_2020_063847_2020_A01). (<b>f</b>) Ripple pattern characteristic of jitter in the cross-track direction (DTEEC_039216_1835_039849_1835_A01). (<b>g</b>) Linear features along CCD seam boundaries due to along-track jitter (DTEED_063209_1800_063420_1800_A01). (<b>h</b>) Artifacts typical of those that occur where dark shadow boundaries shifted between stereo image acquisitions (DTEEC_002486_1860_001985_1860_U01).</p>
Full article ">Figure 5
<p>Error maps output from <span class="html-italic">autoTriangulation</span>. Dots represent the difference in elevation in meters between individual MOLA PEDRs and the corresponding location in the HiRISE DTM (not rendered). Dots are color-coded in 5-m increments, with yellow indicating the desired error range of 0 to ±5 m. Horizontal and vertical axes are longitude and latitude, respectively. The left panels (“Actual Error”) show the measured differences in elevation between MOLA and the HiRISE model. The right panels (“Predicted Error”) show the elevation differences after translation and rotation of the HiRISE model within <span class="html-italic">autoTriangulation</span> to achieve a minimum RMS error: (<b>a</b>) Actual error after the initial (relative +Z-control) solution from SOCET Set, and the corresponding predicted error from <span class="html-italic">autoTriangulation</span>. (<b>b</b>) Error maps after applying the transformation from <span class="html-italic">autoTriangulation</span> to the control points in SOCET Set and generating a new DTM. The actual error and the predicted error now agree closely, and the actual error is within the expected range. At this point, the solution is considered acceptable and the MST step is completed. Error maps shown for DTEEC_005161_1720_016237_1720_A01.</p>
Full article ">Figure 6
<p>Example of a Figure of Merit (FOM) map, legend, and orthoimage. North is up and illumination is from the upper left in all images. (<b>a</b>) FOM map showing the classified map including deeply shadowed areas, which resulted in no reliable elevation data. (<b>b</b>) One of the corresponding orthoimages from the stereo pair. (<b>a.1</b>) Detail of shadowed crater in the FOM map. (<b>b.1</b>) Detail of orthoimage illustrating deeply shadowed portion of the small crater that resulted in the “Shadowed” classification in the FOM. (<b>a.2</b>) Detail of FOM map showing areas of Low Correlation, Suspicious, and Manually Edited areas. (<b>b.2</b>) Corresponding area in the image to illustrate that the low correlation FOM values likely arose from the featureless, bland, low-contrast surface.</p>
Full article ">Figure 7
<p>PDS “Extras”, or browse product examples for DTEEC_052299_2150_052510_2150. All images shown at same scale with north up, illumination from the upper left. (<b>a</b>) Colorized altimetry draped over shaded relief with legend from annotated version. (<b>b</b>) Shaded relief map. (<b>c</b>) Grayscale elevation. (<b>d</b>) FOM map, the only Extra made at full scale in JPEG2000 format with a detached legend. (<b>e</b>) Orthoimage of one half of the stereo pair, ESP_052299_2150. (<b>f</b>) Corresponding enhanced color (near–infrared, red and blue–green, IRB) orthoimage.</p>
Full article ">Figure 8
<p>Global map of HiRISE DTM locations available in the PDS (and at <a href="https://uahirise.org/dtm" target="_blank">https://uahirise.org/dtm</a>, accessed on 10 March 2022) as of December 2021. The symbol size corresponds to the number of unique observations orthorectified to each DTM, not the footprint of the DTMs, which are too small to be distinct at this scale. The smallest circles represent a DTM (or a typical stereo pair; i.e., two orthoimages). Larger circles represent monitoring sites, with labels indicating the number of orthoimages for sites with the most images. The base map is the global MOLA shaded relief topography, in equirectangular projection, positive east longitude 0<math display="inline"><semantics> <msup> <mrow/> <mo>°</mo> </msup> </semantics></math>–360<math display="inline"><semantics> <msup> <mrow/> <mo>°</mo> </msup> </semantics></math>.</p>
Full article ">Figure 9
<p>HiRISE DTM mosaic made from three adjacent, overlapping HiRISE stereo pairs (ESP_037162_1880 and ESP_036384_1880, ESP_036740_1880 and ESP_037030_1880, ESP_035962_1880 and ESP_036529_1880). (<b>a</b>) Image footprints overlaid with symbols indicating locations of tie and control points used in SOCET Set MST. (<b>b</b>) Shaded relief image with autoTriangulation error map showing comparison to MOLA. Error “e” is elevation difference from MOLA in meters. The mean elevation difference is −0.71 m, with a standard deviation 4.34 m. (<b>c</b>) Colorized shaded relief image illustrating the continuity of elevation values across seams. Tie points in overlapping regions allowed images to be controlled relatively to within the vertical precision of each DTM, resulting in differences at DTM seams of &lt;1 m.</p>
Full article ">Figure 10
<p>Perspective view of a detail of the orthoimage (<b>a</b>) draped over a 3D rendering of the DTM (DTEPC_009689_2645_010084_2645_A01). The shadows in the shaded relief map (<b>b</b>), created with simulated illumination angles matching those of the source image closely resemble the shadows in the image, providing a qualitative measure of the DTM resolution. X and Y values are distance (m). Z values are elevation (m).</p>
Full article ">Figure 11
<p>Lyot crater central peak (50.4<math display="inline"><semantics> <msup> <mrow/> <mo>°</mo> </msup> </semantics></math> N, 29.3<math display="inline"><semantics> <msup> <mrow/> <mo>°</mo> </msup> </semantics></math> E) difference of DTMs DTEEC_ 008823_2310_009245_2310_A01 and DTEEC_027376_2310_027297_2310_A01 displayed over the shaded relief map, illuminated from upper left. The difference map and profile show long-wavelength differences, ranging within ±2 m amplitude, due to low-frequency jitter in one or both models. The shaded relief map is shown on the right without the difference map overlay, for clarity.</p>
Full article ">Figure 12
<p>Example slope map derived from a DTM of Tivat crater (45.9<math display="inline"><semantics> <msup> <mrow/> <mo>°</mo> </msup> </semantics></math> S, 9.54<math display="inline"><semantics> <msup> <mrow/> <mo>°</mo> </msup> </semantics></math> N, DTEEC_012991_1335_013624_1335_A02) overlaid on the orthoimage ESP_013624_1335, demonstrating the effects of multi-step smoothing and downsampling to minimize intrinsic noise in HiRISE DTMs. Recurring slope lineae (RSL) in Tivat crater were shown by [<a href="#B58-remotesensing-14-02403" class="html-bibr">58</a>] to originate in areas of steep (&gt;35<math display="inline"><semantics> <msup> <mrow/> <mo>°</mo> </msup> </semantics></math>) rocky slopes (white oval) and to terminate on slopes within the range of the angle of repose for unconsolidated granular material (black oval).</p>
Full article ">Figure 13
<p>Dune migration in Nili Patera (8.70<math display="inline"><semantics> <msup> <mrow/> <mo>°</mo> </msup> </semantics></math> N, 67.35<math display="inline"><semantics> <msup> <mrow/> <mo>°</mo> </msup> </semantics></math> E) over three Mars years. The upper (solid lines) and lower (dashed lines) edges of stoss slopes traced on orthoimages ESP_017762_1890 (orange lines) and ESP_043779_1980 (red lines). For clarity, only the later image is shown (ESP_043779_1980). Small, white arrows indicate direction of bedform movement. See also <a href="#app1-remotesensing-14-02403" class="html-app">Figure S2</a>.</p>
Full article ">Figure 14
<p>(<b>a</b>) Co-registered DTMs at the site dubbed Inca City (81.4<math display="inline"><semantics> <msup> <mrow/> <mo>°</mo> </msup> </semantics></math> S, 295.8<math display="inline"><semantics> <msup> <mrow/> <mo>°</mo> </msup> </semantics></math> E), with overlapping area in color showing the difference map. Overall, the elevation differences are within ±1 m, indicating a satisfactory alignment of the two stereo models. Black arrows indicate interpolation artifacts due to dark shadows where the stereo correlator failed (e.g., <a href="#remotesensing-14-02403-f004" class="html-fig">Figure 4</a>c). (<b>b</b>) Detail of orthophoto ESP_041121_0985, located at the small black square in (<b>a</b>), showing an araneiform with no apparent changes between MY 30 and MY 32. Profile A–A’ is shown from both DTMs. Locally, the differences between the two models are &lt;1 m, with a mean difference of −0.39 m. Gray shaded regions on the difference plot illustrate ±1 and 2 × RSS EP of both DTMs about the local mean difference.</p>
Full article ">Figure 15
<p>Study of ice accumulation changes at a 60 m diameter crater on the north polar residual cap, 85.6<math display="inline"><semantics> <msup> <mrow/> <mo>°</mo> </msup> </semantics></math> N, 58.4<math display="inline"><semantics> <msup> <mrow/> <mo>°</mo> </msup> </semantics></math> E. Each image was taken during northern summer when seasonal frost was minimal, allowing for the precise measurement of changes in the persistent ice. Images are orthorectified to DTEPC_044872_2655_044728_2655_A01. Illumination is from the lower right in all frames. A relative contrast stretch has been applied to each image.</p>
Full article ">Figure 16
<p>Gully activity in Gasa crater (35.74<math display="inline"><semantics> <msup> <mrow/> <mo>°</mo> </msup> </semantics></math> S, 129.42<math display="inline"><semantics> <msup> <mrow/> <mo>°</mo> </msup> </semantics></math> E) observed in HiRISE IRB images separated by nearly one Mars year. The color makes the changes more obvious. However, the mass movement is not resolvable in topographic data even though small topographic changes are visible. Both images are orthorectified to DTEEC_021584_1440_022217_1440_A01. North is to the top of the page, illumination from upper left.</p>
Full article ">Figure 17
<p>DTM (DTEED_016907_1330_016973_1330_U01) and orthoimage (ESP_016907_1330_ RED_D_01_ORTHO, illumination from upper left) featuring a megabarchan dune in Kaiser crater (46.74<math display="inline"><semantics> <msup> <mrow/> <mo>°</mo> </msup> </semantics></math> S, 20.15<math display="inline"><semantics> <msup> <mrow/> <mo>°</mo> </msup> </semantics></math> E). The west-facing slipface of the dune has many annually active gullies.</p>
Full article ">Figure 18
<p>(<b>a</b>) Estimated number of refereed articles per year that use HiRISE DTMs based on a full text query of the Astrophysics Data System (<a href="https://ui.adsabs.harvard.edu" target="_blank">https://ui.adsabs.harvard.edu</a>, accessed on 10 March 2022). (<b>b</b>) Cumulative number of HiRISE DTMs in the PDS per year.</p>
Full article ">
28 pages, 7805 KiB  
Article
Measurement of Soil Tillage Using UAV High-Resolution 3D Data
by Carla Rebelo and João Nascimento
Remote Sens. 2021, 13(21), 4336; https://doi.org/10.3390/rs13214336 - 28 Oct 2021
Cited by 4 | Viewed by 2812
Abstract
Remote sensing methodologies could contribute to a more sustainable agriculture, such as monitoring soil preparation for cultivation, which should be done properly, according to the topographic characteristics and the crop’s nature. The objectives of this work are to (1) demonstrate the potential of [...] Read more.
Remote sensing methodologies could contribute to a more sustainable agriculture, such as monitoring soil preparation for cultivation, which should be done properly, according to the topographic characteristics and the crop’s nature. The objectives of this work are to (1) demonstrate the potential of unmanned aerial vehicle (UAV) technology in the acquisition of 3D data before and after soil tillage, for the quantification of mobilised soil volume; (2) propose a methodology that enables the co-registration of multi-temporal DTMs that were obtained from UAV surveys; and (3) show the relevance of quality control and positional accuracy assessment in processing and results. An unchanged-area-matching method based on multiple linear regression analysis was implemented to reduce the deviation between the Digital Terrain Models (DTMs) to calculate a more reliable mobilised soil volume. The production of DTMs followed the usual photogrammetric-based Structure from Motion (SfM) workflow; the extraction of fill and cut areas was made through raster spatial modelling and statistical tools to support the analysis. Results highlight that the quality of the differential DTM should be ensured for a reliable estimation of areas and mobilised soil volume. This study is a contribution to the use of multi-temporal DTMs produced from different UAV surveys. Furthermore, it demonstrates the potential of UAV data in the understanding of soil variability within precision agriculture. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Methodology for the quantification of mobilised soil volume.</p>
Full article ">Figure 2
<p>Location of the study area in Portugal, Alentejo: aerial survey and area of interest (<b>a</b>); onsite images captured in June and August (<b>b</b>).</p>
Full article ">Figure 3
<p>UAV system: DJI Phantom 4 Pro (<b>a</b>) and sensor (<b>b</b>).</p>
Full article ">Figure 4
<p>UAV flight for June (<b>a</b>) and August (<b>b</b>). The characterisation of each UAV flight with the flight time and the number of aerial images acquired. The red areas in the flight indicate “small disturbances” on image coverage. On the right side, the pictures show image overlaps, where the colour coding represents the number of images where the camera location appears.</p>
Full article ">Figure 5
<p>Artificial targets on the ground, spatial distribution of artificial GCPs on the ground for June and August, and location of the base station for each survey.</p>
Full article ">Figure 6
<p>Photogrammetric workflow.</p>
Full article ">Figure 7
<p>Volume change workflow.</p>
Full article ">Figure 8
<p>Differential DTM profile: DTM before the tillage operation (dashed line); DTM of soil tillage (filled line); and fill and cut volumes.</p>
Full article ">Figure 9
<p>GCPs (red triangles) used for SfM processing and ICPs (yellow circles) used for accuracy assessment. Disposition of GCPs and ICPs in the June block (<b>a</b>) and August block (<b>b</b>).</p>
Full article ">Figure 10
<p>3D models obtained from the June image block processing: dense point cloud, sparse point cloud, DTM as a triangular mesh and shaded relief, and orthoimage over DTM.</p>
Full article ">Figure 11
<p>Location of unchanged areas. Terrain profile analysis for each DTM crossing two unchanged areas (A1 and A3) and the soil tillage area.</p>
Full article ">Figure 12
<p>3D scatterplots: vertical errors in metres for each unchanged area elevation profile (<b>a</b>) and a fitting plane performed through a planar regression of these errors on a 30-metre grid (<b>b</b>). In these plots, only the elevation points with a distance of 3 m from each other were represented, resulting in 534 points out of a total of 20,124.</p>
Full article ">Figure 13
<p>Residual plots for each unchanged area (point distance of profiles were filtered for a distance of 3 m): north vertical error profiles (<b>a</b>), east vertical error profiles (<b>b</b>), south vertical error profiles (<b>c</b>), and west vertical error profiles (<b>d</b>). Black dots represent the vertical errors between the corrected August DTM (estimated) and the June DTM, and red dots represent the vertical errors between the August DTM (observed) and the June DTM.</p>
Full article ">Figure 14
<p>Contour lines (1-metre) derived from each DTM. June contour lines (JCL) (<b>a</b>); JCL and August contour lines (black line) overlaid (<b>b</b>); JCL and August contour lines derived from DTM corrected using model of Equation (6) (brown line) overlaid (<b>c</b>).</p>
Full article ">Figure 15
<p>DTM and digital slope model of the area of interest in June and August. For the representation of slopes, a quantile classification was used.</p>
Full article ">Figure 16
<p>Differential DTMs: differential DTM without correction (<b>a</b>), and differential DTM obtained from the corrected August DTM (<b>b</b>).</p>
Full article ">Figure 17
<p>Fill and cut areas estimation, without and with DTM co-registration; (<b>a</b>) negative elevation difference values or cut areas; and (<b>b</b>) positive elevation difference values, the soil has been “filled in”.</p>
Full article ">Figure 18
<p>Volume change (m<sup>3</sup> per pixel) obtained from the multiplication of the pixel area (0.49 cm<sup>2</sup>) by elevation difference value.</p>
Full article ">
21 pages, 12221 KiB  
Article
Individual Tree Detection and Qualitative Inventory of a Eucalyptus sp. Stand Using UAV Photogrammetry Data
by André Almeida, Fabio Gonçalves, Gilson Silva, Adriano Mendonça, Maria Gonzaga, Jeferson Silva, Rodolfo Souza, Igor Leite, Karina Neves, Marcus Boeno and Braulio Sousa
Remote Sens. 2021, 13(18), 3655; https://doi.org/10.3390/rs13183655 - 13 Sep 2021
Cited by 15 | Viewed by 4493
Abstract
Digital aerial photogrammetry (DAP) data acquired by unmanned aerial vehicles (UAV) have been increasingly used for forest inventory and monitoring. In this study, we evaluated the potential of UAV photogrammetry data to detect individual trees, estimate their heights (ht), and monitor [...] Read more.
Digital aerial photogrammetry (DAP) data acquired by unmanned aerial vehicles (UAV) have been increasingly used for forest inventory and monitoring. In this study, we evaluated the potential of UAV photogrammetry data to detect individual trees, estimate their heights (ht), and monitor the initial silvicultural quality of a 1.5-year-old Eucalyptus sp. stand in northeastern Brazil. DAP estimates were compared with accurate tree locations obtained with real time kinematic (RTK) positioning and direct height measurements obtained in the field. In addition, we assessed the quality of a DAP-UAV digital terrain model (DTM) derived using an alternative ground classification approach and investigated its performance in the retrieval of individual tree attributes. The DTM built for the stand presented an RMSE of 0.099 m relative to the RTK measurements, showing no bias. The normalized 3D point cloud enabled the identification of over 95% of the stand trees and the estimation of their heights with an RMSE of 0.36 m (11%). However, ht was systematically underestimated, with a bias of 0.22 m (6.7%). A linear regression model, was fitted to estimate tree height from a maximum height metric derived from the point cloud reduced the RMSE by 20%. An assessment of uniformity indices calculated from both field and DAP heights showed no statistical difference. The results suggest that products derived from DAP-UAV may be used to generate accurate DTMs in young Eucalyptus sp. stands, detect individual trees, estimate ht, and determine stand uniformity with the same level of accuracy obtained in traditional forest inventories. Full article
(This article belongs to the Special Issue Applications of Individual Tree Detection (ITD))
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Location of the study area, highlighting tree identification and Ground Control Points (GCPs).</p>
Full article ">Figure 2
<p>(<b>a</b>) Orthomosaic of the study area; (<b>b</b>) DTM generated from RTK measurements; (<b>c</b>) DTM generated from unsupervised classification of the point cloud; (<b>d</b>) DTM generated from supervised classification of the point cloud; and (<b>e</b>–<b>g</b>) Spatial distribution of differences between the DTMs.</p>
Full article ">Figure 3
<p>(<b>a</b>) RGB orthomosaic of a young <span class="html-italic">Eucalyptus</span> stand indicating tree crowns mapped through visual interpretation (black outlined polygons), tree bases obtained in the field with RTK (black dots), and three tops detected with three different normalized point clouds: (<b>b</b>) NPC_RTK (yellow dots), (<b>c</b>) NPC_UAV<math display="inline"><semantics> <msub> <mrow/> <mn>1</mn> </msub> </semantics></math> (white circles), and (<b>d</b>) NPC_UAV<math display="inline"><semantics> <msub> <mrow/> <mn>2</mn> </msub> </semantics></math> (orange circles).</p>
Full article ">Figure 4
<p>Scatterplots (<b>a</b>–<b>c</b>) and boxplots (<b>d</b>) of observed (<math display="inline"><semantics> <mrow> <mi>h</mi> <msub> <mi>t</mi> <mrow> <mi>F</mi> <mi>i</mi> <mi>e</mi> <mi>l</mi> <mi>d</mi> </mrow> </msub> </mrow> </semantics></math>) versus estimated (<math display="inline"><semantics> <msub> <mi>H</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> </semantics></math>) tree heights from UAV photogrammetry, considering three normalized point clouds (NPC). The dashed line in (<b>d</b>) represents the mean value of <math display="inline"><semantics> <mrow> <mi>h</mi> <msub> <mi>t</mi> <mrow> <mi>F</mi> <mi>i</mi> <mi>e</mi> <mi>l</mi> <mi>d</mi> </mrow> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 5
<p>Scatterplot (<b>a</b>) and boxplot (<b>b</b>) of observed (<math display="inline"><semantics> <mrow> <mi>h</mi> <msub> <mi>t</mi> <mrow> <mi>F</mi> <mi>i</mi> <mi>e</mi> <mi>l</mi> <mi>d</mi> </mrow> </msub> </mrow> </semantics></math>) versus estimated (<math display="inline"><semantics> <mrow> <mi>h</mi> <msub> <mi>t</mi> <mrow> <mi>U</mi> <mi>A</mi> <mi>V</mi> </mrow> </msub> </mrow> </semantics></math>) tree heights from UAV photogrammetry in a young <span class="html-italic">Eucalyptus</span> stand. The dashed line in (<b>b</b>) represents the mean value of <math display="inline"><semantics> <mrow> <mi>h</mi> <msub> <mi>t</mi> <mrow> <mi>F</mi> <mi>i</mi> <mi>e</mi> <mi>l</mi> <mi>d</mi> </mrow> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 6
<p>Lorenz curve and uniformity metrics of a young <span class="html-italic">Eucalyptus</span> stand. The metrics are based on tree heights (<span class="html-italic">ht</span>) estimated from (<b>a</b>) traditional field inventory, and (<b>b</b>) UAV photogrammetry data (<math display="inline"><semantics> <msub> <mi>H</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> </semantics></math> in black and <math display="inline"><semantics> <mrow> <mi>h</mi> <msub> <mi>t</mi> <mrow> <mi>U</mi> <mi>A</mi> <mi>V</mi> </mrow> </msub> </mrow> </semantics></math> in red).</p>
Full article ">Figure 7
<p>(<b>a</b>) RGB point cloud of a young <span class="html-italic">Eucalyptus</span> stand; (<b>b</b>) the associated vertical profile; and maps detailing the spatial distribution of tree height as (<b>c</b>) measured in the field (<math display="inline"><semantics> <mrow> <mi>h</mi> <msub> <mi>t</mi> <mrow> <mi>F</mi> <mi>i</mi> <mi>e</mi> <mi>l</mi> <mi>d</mi> </mrow> </msub> </mrow> </semantics></math>) and (<b>d</b>,<b>e</b>) estimated from UAV photogrammetry (<math display="inline"><semantics> <msub> <mi>H</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>h</mi> <msub> <mi>t</mi> <mrow> <mi>U</mi> <mi>A</mi> <mi>V</mi> </mrow> </msub> </mrow> </semantics></math>, respectively).</p>
Full article ">
20 pages, 19601 KiB  
Article
The Effect of LiDAR Sampling Density on DTM Accuracy for Areas with Heavy Forest Cover
by Mihnea Cățeanu and Arcadie Ciubotaru
Forests 2021, 12(3), 265; https://doi.org/10.3390/f12030265 - 25 Feb 2021
Cited by 20 | Viewed by 3105
Abstract
Laser scanning via LiDAR is a powerful technique for collecting data necessary for Digital Terrain Model (DTM) generation, even in densely forested areas. LiDAR observations located at the ground level can be separated from the initial point cloud and used as input for [...] Read more.
Laser scanning via LiDAR is a powerful technique for collecting data necessary for Digital Terrain Model (DTM) generation, even in densely forested areas. LiDAR observations located at the ground level can be separated from the initial point cloud and used as input for the generation of a Digital Terrain Model (DTM) via interpolation. This paper proposes a quantitative analysis of the accuracy of DTMs (and derived slope maps) obtained from LiDAR data and is focused on conditions common to most forestry activities (rough, steep terrain with forest cover). Three interpolation algorithms were tested: Inverse Distance Weighted (IDW), Natural Neighbour (NN) and Thin-Plate Spline (TPS). Research was mainly focused on the issue of point data density. To analyze its impact on the quality of ground surface modelling, the density of the filtered data set was artificially lowered (from 0.89 to 0.09 points/m2) by randomly removing point observations in 10% increments. This provides a comprehensive method of evaluating the impact of LiDAR ground point density on DTM accuracy. While the reduction of point density leads to a less accurate DTM in all cases (as expected), the exact pattern varies by algorithm. The accuracy of the LiDAR-derived DTMs is relatively good even when LiDAR sampling density is reduced to 0.40–0.50 points/m2 (50–60 % of the initial point density), as long as a suitable interpolation algorithm is used (as IDW proved to be less resilient to density reductions below approximately 0.60 points/m2). In the case of slope estimation, the pattern is relatively similar, except the difference in accuracy between IDW and the other two algorithms is even more pronounced than in the case of DTM accuracy. Based on this research, we conclude that LiDAR is an adequate method for collecting morphological data necessary for modelling the ground surface, even when the sampling density is significantly reduced. Full article
Show Figures

Figure 1

Figure 1
<p>Location of test area.</p>
Full article ">Figure 2
<p>Spatial and frequency distribution of the initial LiDAR point cloud. The point density grid was generated at a 10 m resolution, to reduce the amount of noise.</p>
Full article ">Figure 3
<p>Spatial and frequency distribution of the filtered ground point cloud. The point density grid was generated at a 10 m resolution, to reduce the amount of noise.</p>
Full article ">Figure 4
<p>General research workflow.</p>
Full article ">Figure 5
<p>Mean unsigned elevation error as a measure of DTM accuracy, across different point densities. Whiskers represent the 95% confidence interval of mean values. Points are horizontally offset for clarity.</p>
Full article ">Figure 6
<p>Percentage of validation points with an elevation error over 0.25 m.</p>
Full article ">Figure 7
<p>Percentage of validation points with an elevation error over 0.50 m.</p>
Full article ">Figure 8
<p>Mean unsigned elevation error as a measure of slope map accuracy, across different point densities. Whiskers represent the 95% confidence interval of mean values. Points are horizontally offset for clarity.</p>
Full article ">Figure 9
<p>Number of validation points with a slope error over the 5-degree threshold.</p>
Full article ">Figure 10
<p>Number of validation points with a slope error over the 10-degree threshold.</p>
Full article ">
18 pages, 1994 KiB  
Article
Impact of Calibrating Filtering Algorithms on the Quality of LiDAR-Derived DTM and on Forest Attribute Estimation through Area-Based Approach
by Diogo N. Cosenza, Luísa Gomes Pereira, Juan Guerra-Hernández, Adrián Pascual, Paula Soares and Margarida Tomé
Remote Sens. 2020, 12(6), 918; https://doi.org/10.3390/rs12060918 - 12 Mar 2020
Cited by 11 | Viewed by 3230
Abstract
Ground point filtering of the airborne laser scanning (ALS) returns is crucial to derive digital terrain models (DTMs) and to perform ALS-based forest inventories. However, the filtering calibration requires considerable knowledge from users, who normally perform it by trial and error without knowing [...] Read more.
Ground point filtering of the airborne laser scanning (ALS) returns is crucial to derive digital terrain models (DTMs) and to perform ALS-based forest inventories. However, the filtering calibration requires considerable knowledge from users, who normally perform it by trial and error without knowing the impacts of the calibration on the produced DTM and the forest attribute estimation. Therefore, this work aims at calibrating four popular filtering algorithms and assessing their impact on the quality of the DTM and the estimation of forest attributes through the area-based approach. The analyzed filters were the progressive triangulated irregular network (PTIN), weighted linear least-squares interpolation (WLS) multiscale curvature classification (MCC), and the progressive morphological filter (PMF). The calibration was established by the vertical DTM accuracy, the root mean squared error (RMSE) using 3240 high-accuracy ground control points. The calibrated parameter sets were compared to the default ones regarding the quality of the estimation of the plot growing stock volume and the dominant height through multiple linear regression. The calibrated parameters allowed for producing DTM with RMSE varying from 0.25 to 0.26 m, against a variation from 0.26 to 0.30 m for the default parameters. The PTIN was the least affected by the calibration, while the WLS was the most affected. Compared to the default parameter sets, the calibrated sets resulted in dominant height equations with comparable accuracies for the PTIN, while WLS, MCC, and PFM reduced the models’ RMSE by 6.5% to 10.6%. The calibration of PTIN and MCC did not affect the volume estimation accuracy, whereas calibrated WLS and PMF reduced the RMSE by 3.4% to 7.9%. The filter calibration improved the DTM quality for all filters and, excepting PTIN, the filters increased the quality of forest attribute estimation, especially in the case of dominant height. Full article
(This article belongs to the Special Issue LiDAR Remote Sensing of Forest Resources and Wildland Fires)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>(<b>a</b>) Map of plots distributed over the area. (<b>b</b>) Histogram of mean terrain slope within plots used on the filtering calibration assessment (41 plots) and forest modeling assessments (25 plots). (<b>c</b>) Exemplification of a digital terrain model (DTM) with their respective control points within a plot. (<b>d</b>) A showcase example of an eroded DTM for the plot (<b>c</b>).</p>
Full article ">Figure 2
<p>Root mean square error (RMSE) values of the digital terrain models (DTMs) derived from the calibration of the parameters of the progressive triangulated irregular network (PTIN). Grey tiles represent settings that produced eroded DTMs.</p>
Full article ">Figure 3
<p>Root mean square error (RMSE) values of the digital terrain models (DTMs) derived from the calibration of the parameters of the weighted linear least-squares interpolation (WLS). Grey tiles represent settings that produced eroded DTMs.</p>
Full article ">Figure 4
<p>Root mean square error (RMSE) values of the digital terrain models (DTMs) derived from the calibration of the parameters of the multiscale curvature classification (MCC).</p>
Full article ">Figure 5
<p>Root mean square error (RMSE) values of the digital terrain models (DTMs) derived from the calibration of the parameters of the progressive morphological filter (PMF). Grey tiles represent settings that produced eroded DTMs.</p>
Full article ">
16 pages, 3370 KiB  
Article
Evaluating the Temperature Difference Parameter in the SSEBop Model with Satellite-Observed Land Surface Temperature Data
by Lei Ji, Gabriel B. Senay, Naga M. Velpuri and Stefanie Kagone
Remote Sens. 2019, 11(16), 1947; https://doi.org/10.3390/rs11161947 - 20 Aug 2019
Cited by 9 | Viewed by 4547
Abstract
The Operational Simplified Surface Energy Balance (SSEBop) model uses the principle of satellite psychrometry to produce spatially explicit actual evapotranspiration (ETa) with remotely sensed and weather data. The temperature difference (dT) in the model is a predefined parameter quantifying the difference [...] Read more.
The Operational Simplified Surface Energy Balance (SSEBop) model uses the principle of satellite psychrometry to produce spatially explicit actual evapotranspiration (ETa) with remotely sensed and weather data. The temperature difference (dT) in the model is a predefined parameter quantifying the difference between surface temperature at bare soil and air temperature at canopy level. Because dT is derived from the average-sky net radiation based primarily on climate data, validation of the dT estimation is critical for assuring a high-quality ETa product. We used the Moderate Resolution Imaging Spectroradiometer (MODIS) data to evaluate the SSEBop dT estimation for the conterminous United States. MODIS data (2008–2017) were processed to compute the 10-year average land surface temperature (LST) and normalized difference vegetation index (NDVI) at 1 km resolution and 8-day interval. The observed dT (dTo) was computed from the LST difference between hot (NDVI < 0.25) and cold (NDVI > 0.7) pixels within each 2° × 2° sampling block. There were enough hot and cold pixels within each block to create dTo timeseries in the West Coast and South-Central regions. The comparison of dTo and modeled dT (dTm) showed high agreement, with a bias of 0.8 K and a correlation coefficient of 0.88 on average. This study concludes that the dTm estimation from the SSEBop model is reliable, which further assures the accuracy of the ETa estimation. Full article
(This article belongs to the Special Issue Remote Sensing: 10th Anniversary)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Map of the study area showing the 2° × 2° (approximately 200 km) sampling blocks and examples of the 10-year (2008–2017) average land surface temperature (LST) and normalized difference vegetation index (NDVI) images for 12–19 July. (<b>a</b>) LST. (<b>b</b>) NDVI. Data source: Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 6 MYD11A2 and MYD09A1 (<a href="#remotesensing-11-01947-t001" class="html-table">Table 1</a>).</p>
Full article ">Figure 2
<p>Maps of hot and cold pixels created based on the 10-year average (2008–2017) LST and NDVI images. The label of each map indicates the dates of the beginning and the end of the 8-day interval.</p>
Full article ">Figure 3
<p>Sampling blocks used for retrieving LST from hot and cold pixels to calculate observed <span class="html-italic">dT</span> (<span class="html-italic">dT<sub>o</sub></span>). The date of the hot and cold pixels in the map is 12–19 July. (<b>a</b>) West Coast region. (<b>b</b>) South-Central region.</p>
Full article ">Figure 4
<p>Timeseries plots of <span class="html-italic">dT<sub>o</sub></span> and <span class="html-italic">dT</span> (<span class="html-italic">dT<sub>m</sub></span>) for the sample blocks in the West Coast region. The label in each plot indicates the name of the sampling block shown in <a href="#remotesensing-11-01947-f003" class="html-fig">Figure 3</a>.</p>
Full article ">Figure 5
<p>Timeseries plots of <span class="html-italic">dT<sub>o</sub></span> and <span class="html-italic">dT<sub>m</sub></span> for the sample blocks in the South-Central region. The labels in each plot indicates the paired sampling blocks shown in <a href="#remotesensing-11-01947-f003" class="html-fig">Figure 3</a>.</p>
Full article ">
23 pages, 6439 KiB  
Article
Filtering Airborne LiDAR Data Through Complementary Cloth Simulation and Progressive TIN Densification Filters
by Shangshu Cai, Wuming Zhang, Xinlian Liang, Peng Wan, Jianbo Qi, Sisi Yu, Guangjian Yan and Jie Shao
Remote Sens. 2019, 11(9), 1037; https://doi.org/10.3390/rs11091037 - 1 May 2019
Cited by 65 | Viewed by 9478
Abstract
Separating point clouds into ground and non-ground points is a preliminary and essential step in various applications of airborne light detection and ranging (LiDAR) data, and many filtering algorithms have been proposed to automatically filter ground points. Among them, the progressive triangulated irregular [...] Read more.
Separating point clouds into ground and non-ground points is a preliminary and essential step in various applications of airborne light detection and ranging (LiDAR) data, and many filtering algorithms have been proposed to automatically filter ground points. Among them, the progressive triangulated irregular network (TIN) densification filtering (PTDF) algorithm is widely employed due to its robustness and effectiveness. However, the performance of this algorithm usually depends on the detailed initial terrain and the cautious tuning of parameters to cope with various terrains. Consequently, many approaches have been proposed to provide as much detailed initial terrain as possible. However, most of them require many user-defined parameters. Moreover, these parameters are difficult to determine for users. Recently, the cloth simulation filtering (CSF) algorithm has gradually drawn attention because its parameters are few and easy-to-set. CSF can obtain a fine initial terrain, which simultaneously provides a good foundation for parameter threshold estimation of progressive TIN densification (PTD). However, it easily causes misclassification when further refining the initial terrain. To achieve the complementary advantages of CSF and PTDF, a novel filtering algorithm that combines cloth simulation (CS) and PTD is proposed in this study. In the proposed algorithm, a high-quality initial provisional digital terrain model (DTM) is obtained by CS, and the parameter thresholds of PTD are estimated from the initial provisional DTM based on statistical analysis theory. Finally, PTD with adaptive parameter thresholds is used to refine the initial provisional DTM. These contributions of the implementation details achieve accuracy enhancement and resilience to parameter tuning. The experimental results indicate that the proposed algorithm improves performance over their direct predecessors. Furthermore, compared with the publicized improved PTDF algorithms, our algorithm is not only superior in accuracy but also practicality. The fact that the proposed algorithm is of high accuracy and easy-to-use is desirable for users. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Most ground points on discontinuities and hilltops are not preserved (they are misclassified as “off-terrain”) by progressive TIN densification filtering (PTDF) because of the coarse initial terrain. (<b>a</b>) Reference digital terrain model (DTM), (<b>b</b>) DTM generated from PTDF.</p>
Full article ">Figure 2
<p>Dense point cloud. (<b>a</b>) Original point cloud, (<b>b</b>) DSM.</p>
Full article ">Figure 3
<p>Flowchart of the proposed algorithm.</p>
Full article ">Figure 4
<p>Procedure of cloth simulation (CS): (<b>a</b>) Initial state. A cloth is placed above the inverted measurements. (<b>b</b>) Each particle drops to the inverted measurements under the external force, and their displacement is computed based on the second law of Newton. If the cloth particles collide with the measurements, these particles will stop falling and be labeled unmovable particles. (<b>c</b>) Except for unmovable particles, each particle is moved according to the internal force produced by neighboring particles. Steps (<b>b</b>) and (<b>c</b>) are repeated until the maximum height variation of all particles is sufficiently small.</p>
Full article ">Figure 5
<p>Parameterization of rigidness.</p>
Full article ">Figure 6
<p>Extraction of ground seed points by collision detection.</p>
Full article ">Figure 7
<p>Selection result of ground seed points. (<b>a</b>) Original point cloud, (<b>b</b>) ground seed points.</p>
Full article ">Figure 8
<p>Construction result of initial provisional DTM. (<b>a</b>) Ground seed points, (<b>b</b>) initial provisional DTM.</p>
Full article ">Figure 9
<p>Key parameters in PTD. (<b>a</b>) Initial provisional DTM, (<b>b</b>) key parameters.</p>
Full article ">Figure 10
<p>The discrete histogram of terrain slope.</p>
Full article ">Figure 11
<p>Simple demonstration of the maximum distance threshold setting.</p>
Full article ">Figure 12
<p>Mirrored operation.</p>
Full article ">Figure 13
<p>Final DTM.</p>
Full article ">Figure 14
<p>Total error average on 15 samples for 32 filters, in which CR value in the proposed algorithm is set to 0.5 m. The proposed algorithm is marked in red.</p>
Full article ">Figure 15
<p>Total error average on 15 samples for 32 filters, in which CR value in the proposed algorithm is set to 1 m. The proposed algorithm is marked in red.</p>
Full article ">Figure 16
<p>Accuracy comparison of our results and optimally tuned results under different CR values.</p>
Full article ">Figure 17
<p>Results of each group (S11 and S52 are selected as representatives). First column: DSMs, second column: DTMs generated from the references, third column: DTMs produced by the proposed algorithm under CR = 0.5 m, fourth column: DTMs produced by the proposed algorithm under CR = 1 m, fifth column: spatial distributions of the type I and type II errors produced by the proposed algorithm under CR = 0.5 m, sixth column: spatial distributions of the type I and type II errors produced by the proposed algorithm under CR = 1 m.</p>
Full article ">Figure 18
<p>Three types of errors on data with high point density for PTDF, CSF, and our algorithm.</p>
Full article ">Figure 19
<p>Comparison of DTM generated by PTDF, CSF and our algorithm. (<b>a</b>) Reference DTM, (<b>b</b>) DTM generated from PTDF, (<b>c</b>) DTM generated from CSF, (<b>d</b>) DTM generated from our algorithm. Filtering results for the cross-section for (<b>e</b>) reference data, (<b>f</b>) PTDF, (<b>g</b>) CSF, and (<b>h</b>) our algorithm.</p>
Full article ">Figure 20
<p>Accuracy comparison of the proposed algorithm under different point density.</p>
Full article ">Figure 21
<p>Results of representative sample: (<b>a</b>) Overlay result of reference DTM and ground seed points (red points) produced by PTDF, (<b>b</b>) overlay result of reference DTM and ground seed points (red points) produced by the proposed algorithm, (<b>c</b>) initial provisional DTM generated from ground seed points of (<b>a</b>), (<b>d</b>) initial provisional DTM generated from ground seed points of (<b>b</b>).</p>
Full article ">Figure 22
<p>Cloth simulation time and total error at different <math display="inline"><semantics> <mrow> <mi>CR</mi> </mrow> </semantics></math> values.</p>
Full article ">Figure 23
<p>The graph of the relationship between the total error and parameter configuration index.</p>
Full article ">
21 pages, 3360 KiB  
Article
UAV-Based Digital Terrain Model Generation under Leaf-Off Conditions to Support Teak Plantations Inventories in Tropical Dry Forests. A Case of the Coastal Region of Ecuador
by Fernando J. Aguilar, José R. Rivas, Abderrahim Nemmaoui, Alberto Peñalver and Manuel A. Aguilar
Sensors 2019, 19(8), 1934; https://doi.org/10.3390/s19081934 - 25 Apr 2019
Cited by 29 | Viewed by 5942
Abstract
Remote sensing is revolutionizing the way in which forests studies are conducted, and recent technological advances, such as Structure from Motion (SfM) photogrammetry from Unmanned Aerial Vehicle (UAV), are providing more efficient methods to assist in REDD (Reducing Emissions from Deforestation and forest [...] Read more.
Remote sensing is revolutionizing the way in which forests studies are conducted, and recent technological advances, such as Structure from Motion (SfM) photogrammetry from Unmanned Aerial Vehicle (UAV), are providing more efficient methods to assist in REDD (Reducing Emissions from Deforestation and forest Degradation) monitoring and forest sustainable management. The aim of this work was to develop and test a methodology based on SfM from UAV to generate high quality Digital Terrain Models (DTMs) on teak plantations (Tectona grandis Linn. F.) situated in the Coastal Region of Ecuador (dry tropical forest). UAV overlapping images were collected using a DJI Phantom 4 Advanced© quadcopter during the dry season (leaf-off phenological stage) over 58 teak square plots of 36 m side belonging to three different plantations located in the province of Guayas (Ecuador). A workflow consisting of SfM absolute image alignment based on field surveyed ground control points, very dense point cloud generation, ground points filtering and outlier removal, and DTM interpolation from labeled ground points, was accomplished. A very accurate Terrestrial Laser Scanning (TLS) derived ground points were employed as ground reference to estimate the UAV-SfM DTM vertical error in each reference plot. The plot-level obtained DTMs presented low vertical bias and random error (−3.1 cm and 11.9 cm on average, respectively), showing statistically significant greater error in those reference plots with basal area and estimated vegetation coverage above 15 m2/ha and 60%, respectively. To the best of the authors’ knowledge, this is the first study aimed at monitoring of teak plantations located in dry tropical forests from UAV images. It provides valuable information that recommends carrying out the UAV image capture during the leaf-off season to obtain UAV-SfM derived DTMs suitable to serve as ground reference in supporting teak plantations inventories. Full article
(This article belongs to the Special Issue UAV-based 3D Mapping)
Show Figures

Figure 1

Figure 1
<p>Situation map of the three teak plantations located in the province of Guayas (Ecuador): Morondava, El Tecal and Allteak.</p>
Full article ">Figure 2
<p>Flowchart representing the applied UAV-based photogrammetric workflow for 3D point cloud and grid DTM generation.</p>
Full article ">Figure 3
<p>Reference plots geometry and TLS-derived point cloud. (<b>a</b>) Square and circular reference plots for UAV-DAP derived DTM and TLS field work, respectively. The four scans pattern is shown; (<b>b</b>) Semi-automatically segmented TLS point cloud showing Ground (brown) and Vegetation (green) classified points. Case study: reference plot number 1 (El Tecal plantation).</p>
Full article ">Figure 4
<p>Example of an RGB original image taken at a flying height of 50 m (AGL) over the reference plot number 2 (Morondava). Estimated vegetation coverage = 0%. Note in the zoom-in area “a” the completely leaf-off stage of the teak trees (yellow ellipses showing trees inclined due to the perspective view), with a large amount of fallen leaves on the ground, and a detail of a rectangular panel (orange ellipse) used to mark the GCPs. Also note the small size of the trees compared to the two people who are setting one of the TLS scanning positions in the zoom-in area “b”.</p>
Full article ">Figure 5
<p>Orthoimages (5 cm ground pixel size) computed for the reference plot number 13 (Allteak). Estimated vegetation coverage = 83.96%. (<b>a</b>) True color composite (R-G-B) orthoimage. (<b>b</b>) Nir false color composite (Nir-R-G) orthoimage composed of images taken from a Parrot Sequoia© (Parrot SA, Paris, France) multispectral sensor on-board the UAV DJI Phantom 4 Advance© used in this work.</p>
Full article ">Figure 6
<p>Exploded pie diagrams of the distribution of mean values (<b>a</b>) and standard deviations (<b>b</b>) computed from the DTM vertical errors in the 58 reference plots.</p>
Full article ">Figure 7
<p>Z-differences distribution at two different reference plots. The corresponding normal distribution is overlaid in red. (<b>a</b>) Reference plot number 2 (Morondava plantation). (<b>b</b>) Reference plot number 13 (Allteak plantation). Units of horizontal axis in meters.</p>
Full article ">Figure 8
<p>Relationship between estimated vegetation coverage and vertical random error of UAV-DAP DTM. (<b>a</b>) Standard deviation of DTM z-differences. (<b>b</b>) 90th (LE90) percentile error of DTM z-differences.</p>
Full article ">Figure 9
<p>Perspective block diagram of the original photogrammetrically derived ground points (red) and the corresponding DTM interpolation (mesh in black) for the reference plot number 13 (Allteak plantation). Estimated vegetation coverage = 83.96%.</p>
Full article ">
22 pages, 17755 KiB  
Article
Impact of the Acquisition Geometry of Very High-Resolution Pléiades Imagery on the Accuracy of Canopy Height Models over Forested Alpine Regions
by Livia Piermattei, Mauro Marty, Wilfried Karel, Camillo Ressl, Markus Hollaus, Christian Ginzler and Norbert Pfeifer
Remote Sens. 2018, 10(10), 1542; https://doi.org/10.3390/rs10101542 - 25 Sep 2018
Cited by 29 | Viewed by 6003
Abstract
This work focuses on the accuracy estimation of canopy height models (CHMs) derived from image matching of Pléiades stereo imagery over forested mountain areas. To determine the height above ground and hence canopy height in forest areas, we use normalised digital surface models [...] Read more.
This work focuses on the accuracy estimation of canopy height models (CHMs) derived from image matching of Pléiades stereo imagery over forested mountain areas. To determine the height above ground and hence canopy height in forest areas, we use normalised digital surface models (nDSMs), computed as the differences between external high-resolution digital terrain models (DTMs) and digital surface models (DSMs) from Pléiades image matching. With the overall goal of testing the operational feasibility of Pléiades images for forest monitoring over mountain areas, two questions guide this work whose answers can help in identifying the optimal acquisition planning to derive CHMs. Specifically, we want to assess (1) the benefit of using tri-stereo images instead of stereo pairs, and (2) the impact of different viewing angles and topography. To answer the first question, we acquired new Pléiades data over a study site in Canton Ticino (Switzerland), and we compare the accuracies of CHMs from Pléiades tri-stereo and from each stereo pair combination. We perform the investigation on different viewing angles over a study area near Ljubljana (Slovenia), where three stereo pairs were acquired at one-day offsets. We focus the analyses on open stable and on tree covered areas. To evaluate the accuracy of Pléiades CHMs, we use CHMs from aerial image matching and airborne laser scanning as reference for the Ticino and Ljubljana study areas, respectively. For the two study areas, the statistics of the nDSMs in stable areas show median values close to the expected value of zero. The smallest standard deviation based on the median of absolute differences (σMAD) was 0.80 m for the forward-backward image pair in Ticino and 0.29 m in Ljubljana for the stereo images with the smallest absolute across-track angle (−5.3°). The differences between the highest accuracy Pléiades CHMs and their reference CHMs show a median of 0.02 m in Ticino with a σMAD of 1.90 m and in Ljubljana a median of 0.32 m with a σMAD of 3.79 m. The discrepancies between these results are most likely attributed to differences in forest structure, particularly tree height, density, and forest gaps. Furthermore, it should be taken into account that temporal vegetational changes between the Pléiades and reference data acquisitions introduce additional, spurious CHM differences. Overall, for narrow forward–backward angle of convergence (12°) and based on the used software and workflow to generate the nDSMs from Pléiades images, the results show that the differences between tri-stereo and stereo matching are rather small in terms of accuracy and completeness of the CHM/nDSMs. Therefore, a small angle of convergence does not constitute a major limiting factor. More relevant is the impact of a large across-track angle (19°), which considerably reduces the quality of Pléiades CHMs/nDSMs. Full article
Show Figures

Figure 1

Figure 1
<p>Study areas and imaging geometries of the Pléiades data sets (Google Earth preview of the footprints and the satellite’s position).</p>
Full article ">Figure 2
<p>General workflow to generate the nDSMs from Pléiades images.</p>
Full article ">Figure 3
<p>For (<b>a</b>) the Ticino and (<b>b</b>) the Ljubljana study areas, left: the 3D point cloud. Centre: the orthophoto generated from the 4-bands Pléiades images, visualised as true colour RGB, and overlaid with the GCPs (red circles) and CPs (orange circles). The yellow rectangles represent the common regions of interest for all scene combinations. Right: the colour coded, reconstructed DSMs. The blue squares are the 500 × 500 m selected areas.</p>
Full article ">Figure 4
<p>nDSM statistics for stable areas, with results before (<b>top</b>) and after (<b>bottom</b>) LSM. (<b>Left</b>) colour coding of the Pléiades (FNB) and aerial nDSMs (non-stable areas shown in white). (<b>Right</b>) distribution of nDSM heights for the aerial images, and for the Pléiades tri-stereo and all stereo pairs.</p>
Full article ">Figure 5
<p>Height differences (ΔH) between Pléiades and aerial CHMs before (<b>top</b>) and after (<b>bottom</b>) LSM. (<b>Left</b>) spatial distribution for FNB as colour coding (areas outside forest mask shown in white). Centre: distributions of ΔH for tri-stereo and all stereo pairs as histograms. (<b>Right</b>) the same distributions shown as box plots.</p>
Full article ">Figure 6
<p>Box plots of the distributions of ∆H after LSM for different classes of (<b>a</b>) reference canopy height, and (<b>b</b>) terrain slope. Respective cell counts are plotted in grey and refer to the right axes.</p>
Full article ">Figure 7
<p>Profiles of 1 m width of the tri-stereo and stereo Pléiades point clouds (<b>top</b>) and of the Pléiades nDSMs for each stereo pair in comparison to the reference aerial nDSM (<b>bottom</b>). The red line in the orthophoto indicates the position of the profiles.</p>
Full article ">Figure 8
<p>(<b>a</b>) Histograms of the absolute height errors (Abs∆H) for each Pléiades CHM (MinAbs∆H i.e., optimal selection from stereo DSMs, NB, FN, FB, and FB-NB-FN), and (<b>b</b>) the spatial distribution of the stereo pairs with locally minimum absolute differences (MinAbs∆H).</p>
Full article ">Figure 9
<p>The spatial distribution of the nDSMs of ALS (<b>left</b>, <b>top</b>) and each Pléiades stereo scene before (<b>right</b>, <b>top</b>) and after (<b>righ</b>, <b>bottom</b>) LSM. At the (<b>left</b>, <b>bottom</b>), the orthophoto derived from the forward Pléiades image of 27 July is shown.</p>
Full article ">Figure 10
<p>Distribution of nDSM heights in stable areas (<b>left</b>) and CHM errors (∆H) before (<b>right</b>, <b>top</b>) and after (<b>right</b>, <b>bottom</b>) LSM for each stereo pair.</p>
Full article ">Figure 11
<p>(<b>a</b>) The spatial distribution of the three reference CHM tree height classes, disregarding areas with large height variation. Distributions of the CHM error (∆H) for each stereo pair after LSM grouped by (<b>b</b>) tree height class and (<b>c</b>) terrain slope.</p>
Full article ">Figure 12
<p>For the selected (<b>a</b>) area 1 and (<b>b</b>) area 2, the reference CHM (left), the Pléiades CHMs for each stereo pair (centre), and the orthophoto generated from the four bands Pléiades image (right). The dashed black line in ((<b>a</b>), left) indicates the position of the profiles of 1 m width of reference and Pléiades CHMs shown on the left of (<b>c</b>). On the right of (<b>c</b>), the aerial orthophoto for area 2 that was acquired simultaneously to the reference data.</p>
Full article ">Figure 13
<p>2D histograms of Pléiades CHM errors and reference CHM for the selected (<b>a</b>) area 1 and (<b>b</b>) area 2, and each stereo pair (left, centre, right).</p>
Full article ">Figure A1
<p>For Ticino study area (<b>a</b>) the image residual vectors (scale 5000) for GCPs (red circle), CPs (yellow circle) and tie points (blue circle) at their respective image positions after the RPC correction. (<b>b</b>–<b>d</b>) The 2D scatter plot of the image residuals of the tie points of each image in pixel units.</p>
Full article ">Figure A2
<p>For Ljubljana study area (<b>a</b>) Image residual vectors (scale 5000) for GCPs (red circle), CPs (yellow circle) and tie points (blue circle) at their respective image positions after the RPC correction. (<b>b</b>) The 2D distribution of image residuals of tie points of each image.</p>
Full article ">
6203 KiB  
Article
Impacts of Airborne Lidar Pulse Density on Estimating Biomass Stocks and Changes in a Selectively Logged Tropical Forest
by Carlos Alberto Silva, Andrew Thomas Hudak, Lee Alexander Vierling, Carine Klauberg, Mariano Garcia, António Ferraz, Michael Keller, Jan Eitel and Sassan Saatchi
Remote Sens. 2017, 9(10), 1068; https://doi.org/10.3390/rs9101068 - 23 Oct 2017
Cited by 55 | Viewed by 9566
Abstract
Airborne lidar is a technology well-suited for mapping many forest attributes, including aboveground biomass (AGB) stocks and changes in selective logging in tropical forests. However, trade-offs still exist between lidar pulse density and accuracy of AGB estimates. We assessed the impacts of lidar [...] Read more.
Airborne lidar is a technology well-suited for mapping many forest attributes, including aboveground biomass (AGB) stocks and changes in selective logging in tropical forests. However, trade-offs still exist between lidar pulse density and accuracy of AGB estimates. We assessed the impacts of lidar pulse density on the estimation of AGB stocks and changes using airborne lidar and field plot data in a selectively logged tropical forest located near Paragominas, Pará, Brazil. Field-derived AGB was computed at 85 square 50 × 50 m plots in 2014. Lidar data were acquired in 2012 and 2014, and for each dataset the pulse density was subsampled from its original density of 13.8 and 37.5 pulses·m−2 to lower densities of 12, 10, 8, 6, 4, 2, 0.8, 0.6, 0.4 and 0.2 pulses·m−2. For each pulse density dataset, a power-law model was developed to estimate AGB stocks from lidar-derived mean height and corresponding changes between the years 2012 and 2014. We found that AGB change estimates at the plot level were only slightly affected by pulse density. However, at the landscape level we observed differences in estimated AGB change of >20 Mg·ha−1 when pulse density decreased from 12 to 0.2 pulses·m−2. The effects of pulse density were more pronounced in areas of steep slope, especially when the digital terrain models (DTMs) used in the lidar derived forest height were created from reduced pulse density data. In particular, when the DTM from high pulse density in 2014 was used to derive the forest height from both years, the effects on forest height and the estimated AGB stock and changes did not exceed 20 Mg·ha−1. The results suggest that AGB change can be monitored in selective logging in tropical forests with reasonable accuracy and low cost with low pulse density lidar surveys if a baseline high-quality DTM is available from at least one lidar survey. We recommend the results of this study to be considered in developing projects and national level MRV systems for REDD+ emission reduction programs for tropical forests. Full article
(This article belongs to the Special Issue Biomass Remote Sensing in Forest Landscapes)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Location of the study area. (<b>A</b>) South America and Brazil; (<b>B</b>) States of Pará and Paragominas city; (<b>C</b>) Paragominas city; (<b>D</b>) Airborne lidar coverage; (<b>E</b>) Field plots on the lidar-derived CHM. Reduced-impact logging (RIL).</p>
Full article ">Figure 2
<p>A 3D illustration of airborne lidar pulse density reduction at the plot level (0.25 ha) in 2014.</p>
Full article ">Figure 3
<p>Flowchart of the lidar data processing for AGB stocks and AGB change estimation in tropical forest. The green panel to the left shows the lidar data processing (<b>a</b>) and the gray panel to the right shows the AGB stocks and change estimation steps (<b>b</b>).</p>
Full article ">Figure 4
<p>Lidar-derived HMEAN (m) (<b>a1-2</b>); Standard deviation of HMEAN (m) for the sample plots (30 repetitions) (<b>b1-2</b>); Reliability Ratio for HMEAN (<b>c1-2</b>); 2012 (<b>a1</b>–<b>c1</b>) and 2014 (<b>a2</b>–<b>c2</b>); (<span class="html-italic">n</span> = 84). DS1 (orange): DTM scenario 1; DS2 (green): DTM scenario 2.</p>
Full article ">Figure 5
<p>Boxplot of R<sup>2</sup> (<b>a</b>), relative (<b>b1</b>–<b>c1</b>) and absolute (<b>b2</b>–<b>c2</b>) RMSE and bias for the AGB leave-one-out cross validation – LOOCV models. DS1 (orange): DTM Scenario 1; DS2 (green): DTM Scenario 2.</p>
Full article ">Figure 6
<p>Boxplot of the AGB estimates for 2012 and 2014 (<b>a1</b>,<b>b1</b>), and AGB change (<b>c1</b>). Standard deviation of AGB stock in 2012 (<b>a2</b>), 2014 (<b>b2</b>) and AGB change (<b>c2</b>) (30 repetitions) (<span class="html-italic">n</span> = 84).</p>
Full article ">Figure 7
<p>Map of AGB stock in 2012 (<b>a</b>), 2014 (<b>b</b>) and AGB change (<b>c</b>) at 12 pulse·m<sup>−2</sup> in DS2. Zoom in the AGB change maps derived at 0.2 and 12 pulse·m<sup>−2</sup> in DS1 (<b>d</b>,<b>e</b>) and DS2 (<b>f</b>,<b>g</b>) in an unlogged and logged unit. The maps were calculated as the mean of the 30 replicates.</p>
Full article ">Figure 8
<p>Digital terrain model (<b>a</b>) and Slope (%) (<b>b</b>) maps of the study area at 12 pulses·m<sup>−2</sup> in 2012; Standard deviation of AGB change at 0.2 pulses·m<sup>−2</sup> for DS1 (<b>c</b>–<b>c1</b>) and DS2 (<b>d</b>–<b>d1</b>).</p>
Full article ">Figure 9
<p>Boxplot of the differences in predicted AGB change at stand level of 12 m<sup>−2</sup> degraded to 0.2, 0.4, 0.6, 0.8, 2, 4, 6, 8 and 10 pulses·m<sup>−2</sup> in areas with slopes ranging from 0 to 12% (<b>a</b>), 12–24% (<b>b</b>) and 24–36% (<b>c</b>), under DS1 (orange) and DS2 (dark green).</p>
Full article ">
Back to TopTop