[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
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (225)

Search Parameters:
Keywords = digital orthophoto

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
29 pages, 96249 KiB  
Article
SAR-MINF: A Novel SAR Image Descriptor and Matching Method for Large-Scale Multidegree Overlapping Tie Point Automatic Extraction
by Shuo Li, Xiongwen Yang, Xiaolei Lv and Jian Li
Remote Sens. 2024, 16(24), 4696; https://doi.org/10.3390/rs16244696 - 16 Dec 2024
Viewed by 324
Abstract
The automatic extraction of large-scale tie points (TPs) for Synthetic Aperture Radar (SAR) images is crucial for generating SAR Digital Orthophoto Maps (DOMs). This task involves matching SAR images under various conditions, such as different resolutions, incidence angles, and orbital directions, which is [...] Read more.
The automatic extraction of large-scale tie points (TPs) for Synthetic Aperture Radar (SAR) images is crucial for generating SAR Digital Orthophoto Maps (DOMs). This task involves matching SAR images under various conditions, such as different resolutions, incidence angles, and orbital directions, which is highly challenging. To address the feature extraction challenges of different SAR images, we propose a Gamma Modulated Phase Congruency (GMPC) model. This improved phase congruency model is defined by a Gamma Modulation Filter (GMF) and an adaptive noise model. Additionally, to reduce layover interference in SAR images, we introduce a GMPC-Harris feature point extraction method with layover perception. We also propose a matching method based on the SAR Modality Independent Neighborhood Fusion (SAR-MINF) descriptor, which fuses feature information from different neighborhoods. Furthermore, we present a graph-based overlap extraction algorithm and establish an automated workflow for large-scale TP extraction. Experiments show that the proposed SAR-MINF matching method increases the Correct Match Rate (CMR) by an average of 31.2% and the matching accuracy by an average of 57.8% compared with other prevalent SAR image matching algorithms. The proposed TP extraction algorithm can extract full-degree TPs with an accuracy of less than 0.5 pixels for more than 98% of 2-degree TPs and over 95% of multidegree TPs, meeting the requirements of DOM production. Full article
Show Figures

Figure 1

Figure 1
<p>The Fourier series expansions for (<b>a</b>) square waves and (<b>b</b>) triangular waves are depicted, where the black solid line represents the original function, the blue solid line represents the sum of the first four terms of the Fourier series expansion, and the dashed lines represent the individual Fourier series expansion terms.</p>
Full article ">Figure 2
<p>The two-dimensional Gamma Modulation Filter consists of (<b>a</b>) the Gamma kernel function part, (<b>b</b>) the odd part <math display="inline"><semantics> <mrow> <mi>G</mi> <mi>M</mi> <msub> <mi>F</mi> <mi>o</mi> </msub> </mrow> </semantics></math>, and (<b>c</b>) the even part <math display="inline"><semantics> <mrow> <mi>G</mi> <mi>M</mi> <msub> <mi>F</mi> <mi>e</mi> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>(<b>a</b>) Simulated SAR images with multiplicative noise and the definitions of (<b>b</b>) <math display="inline"><semantics> <msub> <mi>o</mi> <mrow> <mi>s</mi> <mi>a</mi> <mi>r</mi> </mrow> </msub> </semantics></math>, (<b>c</b>) <math display="inline"><semantics> <msub> <mi>e</mi> <mrow> <mi>s</mi> <mi>a</mi> <mi>r</mi> </mrow> </msub> </semantics></math>, and (<b>d</b>) SAR local energy <math display="inline"><semantics> <msub> <mi>E</mi> <mrow> <mi>s</mi> <mi>a</mi> <mi>r</mi> </mrow> </msub> </semantics></math> with two-dimensional Gamma Modulation Filter at <math display="inline"><semantics> <mrow> <mi>θ</mi> <mo>=</mo> <mn>0</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>Comparison of (<b>b</b>) GMPC, (<b>c</b>) PC, and (<b>d</b>) SAR-PC on (<b>a</b>) real SAR images.</p>
Full article ">Figure 5
<p>SAR-MIND Matching Process.</p>
Full article ">Figure 6
<p>Comparison of different keypoint extraction algorithms: (<b>a</b>) GMPC-Harris, (<b>b</b>) SAR-Harris, (<b>c</b>) Harris, (<b>d</b>) SURF, (<b>e</b>) FAST.</p>
Full article ">Figure 7
<p>(<b>a</b>) The layover area in SAR image and (<b>b</b>) its geometric relationships.</p>
Full article ">Figure 8
<p>Comparison of <math display="inline"><semantics> <msub> <mi>H</mi> <mrow> <mi>G</mi> <mi>M</mi> <mi>P</mi> <mi>C</mi> </mrow> </msub> </semantics></math> between layover and normal area. (<b>b</b>,<b>d</b>,<b>f</b>,<b>h</b>) is the GMPC orientation histogram of the yellow region in (<b>a</b>,<b>c</b>,<b>e</b>,<b>g</b>), respectively.</p>
Full article ">Figure 9
<p>GMPC-Harris keypoints based on Layover perception. (<b>a</b>) Mountain area. (<b>b</b>) Building area.</p>
Full article ">Figure 10
<p>Sampling space for different feature descriptors: (<b>a</b>) DAISY, (<b>b</b>) GLOH, (<b>c</b>) SAR-MINF.</p>
Full article ">Figure 11
<p>The construction process of the SAR-MINF descriptor.</p>
Full article ">Figure 12
<p>Flowchart of large-scale tie point automatic extraction.</p>
Full article ">Figure 13
<p>Experimental images for SAR-MINF matching. (<b>a</b>–<b>h</b>) Pair A–H.</p>
Full article ">Figure 14
<p>(<b>a</b>) CMR and (<b>b</b>) RMSE of different algorithms.</p>
Full article ">Figure 15
<p>Data from five experiments of RIFT with pair G.</p>
Full article ">Figure 16
<p>Matching results of Pair F: (<b>a</b>) SAR-MINF, (<b>b</b>) NCC, (<b>c</b>) SAR-PC, (<b>d</b>) SAR-SIFT, (<b>e</b>) KAZE-SAR, (<b>f</b>) RIFT.</p>
Full article ">Figure 17
<p>Enlarged checkboard mosaic sub-images of pair F: (<b>a</b>) SAR-MINF, (<b>b</b>) NCC, (<b>c</b>) SAR-PC, (<b>d</b>) SAR-SIFT, (<b>e</b>) KAZE-SAR, (<b>f</b>) RIFT.</p>
Full article ">Figure 18
<p>Checkboard mosaic images and enlarged sub-images of all images under the SAR-MINF algorithm. (<b>a</b>–<b>h</b>) Pair A–H.</p>
Full article ">Figure 19
<p>(<b>a</b>) CMR and (<b>b</b>) RMSE of the varying search radius.</p>
Full article ">Figure 20
<p>Ablation experiment for neighborhood fusion. (<b>a</b>) CMR. (<b>b</b>) RMSE.</p>
Full article ">Figure 21
<p>Geographical distribution of (<b>a</b>) the first and (<b>b</b>) the second set of TPs test data.</p>
Full article ">Figure 22
<p>The multidegree overlapping graph of the first set of images.</p>
Full article ">Figure 23
<p>Histogram of RMSE distribution for TPs. (<b>a</b>) First group of 2-degree TPs. (<b>b</b>) First group of multi-degree TPs. (<b>c</b>) Second group of 2-degree TPs. (<b>d</b>) Second group of multi-degree TPs.</p>
Full article ">Figure 24
<p>Partial multidegree overlapping TPs’ slices. (<b>a</b>–<b>c</b>) 3-degree overlapping TPs, (<b>d</b>–<b>f</b>) 4-degree overlapping TPs, (<b>g</b>–<b>i</b>) 5-degree overlapping TPs.</p>
Full article ">
22 pages, 24817 KiB  
Article
Construction of Mining Subsidence Basin and Inversion of Predicted Subsidence Parameters Based on UAV Photogrammetry Products Considering Horizontal Displacement
by Jinqi Zhao, Yufen Niu, Zhengpei Zhou, Zhong Lu, Zhimou Wang, Zhaojiang Zhang, Yiyao Li and Ziheng Ju
Remote Sens. 2024, 16(22), 4283; https://doi.org/10.3390/rs16224283 - 17 Nov 2024
Viewed by 526
Abstract
Constructing high-precision subsidence basins is of paramount importance for mining subsidence monitoring. Traditional unmanned aerial vehicle (UAV) photogrammetry techniques typically construct subsidence basins by directly differencing digital elevation models (DEMs) from different monitoring periods. However, this method often neglects the influence of horizontal [...] Read more.
Constructing high-precision subsidence basins is of paramount importance for mining subsidence monitoring. Traditional unmanned aerial vehicle (UAV) photogrammetry techniques typically construct subsidence basins by directly differencing digital elevation models (DEMs) from different monitoring periods. However, this method often neglects the influence of horizontal displacement on the accuracy of the subsidence basin. Taking a mining area in Ordos, Inner Mongolia, as an example, this study employed the normalized cross-correlation (NCC) matching algorithm to extract horizontal displacement information between two epochs of a digital orthophoto map (DOM) and subsequently corrected the horizontal position of the second-epoch DEM. This ensured that the planar positions of ground feature points remained consistent in the DEM before and after subsidence. Based on this, the vertical displacement in the subsidence area (the subsidence basin) was obtained via DEM differencing, and the parameters of the post-correction subsidence basin were inverted using the probability integral method (PIM). The experimental results indicate that (1) the horizontal displacement was influenced by the gully topography, causing the displacement within the working face to be segmented on both sides of the gully; (2) the influence of the terrain on the subsidence basin was significantly reduced after correction; (3) the post-correction surface subsidence curve was smoother than the pre-correction curve, with abrupt error effects markedly diminished; (4) the accuracy of the post-correction subsidence basin increased by 43.12% compared with the total station data; and (5) comparing the measured horizontal displacement curve with that derived using the probability integral method revealed that the horizontal displacement on the side of an old goaf adjacent to the newly excavated working face shifted toward the advancing direction of the new working face as mining progressed. This study provides a novel approach and insights for using low-cost UAVs to construct high-precision subsidence basins. Full article
Show Figures

Figure 1

Figure 1
<p>Schematic diagram of the study area location. (<b>a</b>) Map of China; (<b>b</b>) DEM of Ordos; (<b>c</b>) study area.</p>
Full article ">Figure 2
<p>Technical flow chart of this research.</p>
Full article ">Figure 3
<p>Schematic diagram of the DEM correction process.</p>
Full article ">Figure 4
<p>(<b>a</b>) East–west displacement; (<b>b</b>) north–south displacement.</p>
Full article ">Figure 5
<p>Illustration of the relationship between horizontal displacement and topography. (<b>a</b>,<b>b</b>) are cross-sectional views of profile A-A′; (<b>c</b>,<b>d</b>) are cross-sectional views of profile B-B′.</p>
Full article ">Figure 6
<p>Horizontal displacement in gully topography. (<b>a</b>) A-A′ cross-section; (<b>b</b>) local displacement field.</p>
Full article ">Figure 7
<p>Subsidence basin. (<b>a</b>) Pre-correction subsidence basin; (<b>b</b>) post-correction subsidence basin.</p>
Full article ">Figure 8
<p>Local maps of areas I and II. (<b>a</b>) Magnified view of area I pre-correction; (<b>b</b>) magnified view of area I post-correction; (<b>c</b>) magnified view of area II pre-correction; (<b>d</b>) magnified view of area II post-correction; (<b>e</b>) 1-1′ cross-section; (<b>f</b>) 2-2′ cross-section.</p>
Full article ">Figure 9
<p>Subsidence curves of pre-correction and post-correction. (<b>a</b>) A-A′ cross-section; (<b>b</b>) C-C′ cross-section.</p>
Full article ">Figure 10
<p>Inverted subsidence basin.</p>
Full article ">Figure 11
<p>Measured subsidence basin.</p>
Full article ">Figure 12
<p>(<b>a</b>) Strike main profile; (<b>b</b>) dip main profile.</p>
Full article ">Figure 13
<p>Horizontal displacement of strike main profile. (<b>a</b>) Strike main profile; (<b>b</b>) partial enlarged detail.</p>
Full article ">Figure 14
<p>Horizontal displacement of dip main profile. (<b>a</b>) Dip main profile; (<b>b</b>) partial enlarged detail.</p>
Full article ">Figure 15
<p>Horizontal displacement error.</p>
Full article ">Figure 16
<p>Statistical chart of residuals for subsidence basin.</p>
Full article ">Figure 17
<p>Statistical chart of strike residuals.</p>
Full article ">Figure 18
<p>Statistical chart of dip residuals.</p>
Full article ">Figure 19
<p>Statistical analysis of errors in subsidence basin.</p>
Full article ">
18 pages, 4434 KiB  
Article
Monitoring of Heracleum sosnowskyi Manden Using UAV Multisensors: Case Study in Moscow Region, Russia
by Rashid K. Kurbanov, Arkady N. Dalevich, Alexey S. Dorokhov, Natalia I. Zakharova, Nazih Y. Rebouh, Dmitry E. Kucher, Maxim A. Litvinov and Abdelraouf M. Ali
Agronomy 2024, 14(10), 2451; https://doi.org/10.3390/agronomy14102451 - 21 Oct 2024
Viewed by 864
Abstract
Detection and mapping of Sosnowsky’s hogweed (HS) using remote sensing data have proven effective, yet challenges remain in identifying, localizing, and eliminating HS in urban districts and regions. Reliable data on HS growth areas are essential for monitoring, eradication, and control measures. Satellite [...] Read more.
Detection and mapping of Sosnowsky’s hogweed (HS) using remote sensing data have proven effective, yet challenges remain in identifying, localizing, and eliminating HS in urban districts and regions. Reliable data on HS growth areas are essential for monitoring, eradication, and control measures. Satellite data alone are insufficient for mapping the dynamics of HS distribution. Unmanned aerial vehicles (UAVs) with high-resolution spatial data offer a promising solution for HS detection and mapping. This study aimed to develop a method for detecting and mapping HS growth areas using a proposed algorithm for thematic processing of multispectral aerial imagery data. Multispectral data were collected using a DJI Matrice 200 v2 UAV (Dajiang Innovation Technology Co., Shenzhen, China) and a MicaSense Altum multispectral camera (MicaSense Inc., Seattle, WA, USA). Between 2020 and 2022, 146 sites in the Moscow region of the Russian Federation, covering 304,631 hectares, were monitored. Digital maps of all sites were created, including 19 digital maps (orthophoto, 5 spectral maps, and 13 vegetation indices) for four experimental sites. The collected samples included 1080 points categorized into HS, grass cover, and trees. Student’s t-test showed significant differences in vegetation indices between HS, grass, and trees. A method was developed to determine and map HS-growing areas using the selected vegetation indices NDVI > 0.3, MCARI > 0.76, user index BS1 > 0.10, and spectral channel green > 0.14. This algorithm detected HS in an area of 146.664 hectares. This method can be used to monitor and map the dynamics of HS distribution in the central region of the Russian Federation and to plan the required volume of pesticides for its eradication. Full article
Show Figures

Figure 1

Figure 1
<p>Study area. Moscow region is highlighted in red.</p>
Full article ">Figure 2
<p>UAV DJI Matrice 200 v2 with a GNSS receiver Topodrone PPK.</p>
Full article ">Figure 3
<p>Additional equipment: (1) GNSS Emlid Reach RS2; (2) ground control point; (3) calibration panel for multispectral camera; and (4) MicaSense Altum multispectral camera.</p>
Full article ">Figure 4
<p>File structure of data storage folders.</p>
Full article ">Figure 5
<p>The studied samples: (<b>a</b>) Sosnowsky’s hogweed; (<b>b</b>) grass; (<b>c</b>) trees.</p>
Full article ">Figure 6
<p>Algorithm for identification of HS plants on digital maps with UAVs.</p>
Full article ">Figure 7
<p>Display of HS thickets: (<b>a</b>) a custom BS1 index with values above 0.11; (<b>b</b>) an orthophoto, with HS marked in red; (<b>c</b>) the highlighting of single HS plants by the MCARI index; and (<b>d</b>) the highlighting of single HS plants by the spectral channel green.</p>
Full article ">Figure 8
<p>Spectral characteristics of Sosnowsky’s hogweed, trees, and grass.</p>
Full article ">Figure 9
<p>Comparison of vegetation indices values.</p>
Full article ">
15 pages, 12314 KiB  
Article
Low-Cost Photogrammetry for Detailed Documentation and Condition Assessment of Earthen Architectural Heritage: The Ex-Hotel Oasis Rouge in Timimoun as a Case Study
by Haroune Ben Charif, Ornella Zerlenga and Rosina Iaderosa
Buildings 2024, 14(10), 3292; https://doi.org/10.3390/buildings14103292 - 17 Oct 2024
Viewed by 784
Abstract
Earthen architecture holds deep historical, cultural, and ecological value, forming an essential component of our global cultural heritage. However, these structures face numerous threats, including climate change, socio-economic shifts, and, in many cases, neglection, which accelerate their deterioration. This study introduces a photogrammetry-based [...] Read more.
Earthen architecture holds deep historical, cultural, and ecological value, forming an essential component of our global cultural heritage. However, these structures face numerous threats, including climate change, socio-economic shifts, and, in many cases, neglection, which accelerate their deterioration. This study introduces a photogrammetry-based methodology adapted for the digital documentation and preservation of earthen architecture within the context of developing countries. We focus on the Ex-Hotel Oasis Rouge in Timimoun, an iconic earthen building in southwestern Algeria and the current headquarters of CAPTERRE (Algerian Centre for Earthen Built Cultural Heritage). This paper details our approach to using photogrammetry to capture both the interior and exterior of the building, produce detailed orthophotos for archiving the unique earthen bas-reliefs, and conduct a condition assessment. Our findings demonstrate the effectiveness of photogrammetry as a cost-effective tool for heritage documentation, highlighting its potential to assist in the ongoing preservation and informed restoration of earthen architecture. Full article
Show Figures

Figure 1

Figure 1
<p>Satellite images showing the rapid transformation of the Ksar of Timimoun from 2006 (<b>left</b>) to 2023 (<b>right</b>), demonstrating the rapid changes in the urban landscape due to shifts from the use of local materials to industrial building materials. Source: Google Earth.</p>
Full article ">Figure 2
<p>Methodology overview.</p>
Full article ">Figure 3
<p>Position of the city of Timimoun within its Algerian and north African context. Source: D-maps.com (accessed on 20 August 2024), elaborated by authors.</p>
Full article ">Figure 4
<p>View from the interior of the Ex-Hotel Oasis Rouge, into the room that once hosted Grand Duchess Charlotte of Luxembourg. The door features a photograph from 1926, showing the Grand Duchess with Prince Sixtus of Bourbon-Parma on the terrace of the Hotel Transatlantique.</p>
Full article ">Figure 5
<p>(<b>Left</b>): Historical photograph of Hotel SATT, around 1950. This anonymous photograph was cited in [<a href="#B23-buildings-14-03292" class="html-bibr">23</a>]. (<b>Right</b>): Recent photograph of the Ex-Hotel Oasis Rouge, now the headquarters of CAPTERRE.</p>
Full article ">Figure 6
<p>Interior wall finishes of the Ex-Hotel Oasis Rouge, highlighting the earthen hand-crafted bas-relief geometric motifs.</p>
Full article ">Figure 7
<p>Views of the collapsed arches during the 2016 maintenance work. The arches have been supported by temporary formwork since the collapse, which remains in place to this day.</p>
Full article ">Figure 8
<p>Point clouds and 3D model of the Ex-Hotel Oasis Rouge, showcasing details from the interior and exterior.</p>
Full article ">Figure 9
<p>Three-dimensional view of the point cloud of the Ex-Hotel Oasis Rouge.</p>
Full article ">Figure 10
<p>One of the resulting orthophotos of the main corridor of the Ex-Hotel Oasis Rouge.</p>
Full article ">Figure 11
<p>Plan view of the Ex-Hotel Oasis Rouge (<b>bottom</b>) showing locations of orthophotos within corresponding zones. (<b>Above</b>), orthophoto B2-13 is displayed with ×2 and ×4 enlarged details, highlighting the resolution and detail captured.</p>
Full article ">Figure 12
<p>Example of a detailed condition assessment for an earthen bas-relief based on an orthophoto (B2-12).</p>
Full article ">
19 pages, 12162 KiB  
Article
DRA-UNet for Coal Mining Ground Surface Crack Delineation with UAV High-Resolution Images
by Wei Wang, Weibing Du, Xiangyang Song, Sushe Chen, Haifeng Zhou, Hebing Zhang, Youfeng Zou, Junlin Zhu and Chaoying Cheng
Sensors 2024, 24(17), 5760; https://doi.org/10.3390/s24175760 - 4 Sep 2024
Viewed by 791
Abstract
Coal mining in the Loess Plateau can very easily generate ground cracks, and these cracks can immediately result in ventilation trouble under the mine shaft, runoff disturbance, and vegetation destruction. Advanced UAV (Unmanned Aerial Vehicle) high-resolution mapping and DL (Deep Learning) are introduced [...] Read more.
Coal mining in the Loess Plateau can very easily generate ground cracks, and these cracks can immediately result in ventilation trouble under the mine shaft, runoff disturbance, and vegetation destruction. Advanced UAV (Unmanned Aerial Vehicle) high-resolution mapping and DL (Deep Learning) are introduced as the key methods to quickly delineate coal mining ground surface cracks for disaster prevention. Firstly, the dataset named the Ground Cracks of Coal Mining Area Unmanned Aerial Vehicle (GCCMA-UAV) is built, with a ground resolution of 3 cm, which is suitable to make a 1:500 thematic map of the ground crack. This GCCMA-UAV dataset includes 6280 images of ground cracks, and the size of the imagery is 256 × 256 pixels. Secondly, the DRA-UNet model is built effectively for coal mining ground surface crack delineation. This DRA-UNet model is an improved UNet DL model, which mainly includes the DAM (Dual Dttention Dechanism) module, the RN (residual network) module, and the ASPP (Atrous Spatial Pyramid Pooling) module. The DRA-UNet model shows the highest recall rate of 77.29% when the DRA-UNet was compared with other similar DL models, such as DeepLabV3+, SegNet, PSPNet, and so on. DRA-UNet also has other relatively reliable indicators; the precision rate is 84.92% and the F1 score is 78.87%. Finally, DRA-UNet is applied to delineate cracks on a DOM (Digital Orthophoto Map) of 3 km2 in the mining workface area, with a ground resolution of 3 cm. There were 4903 cracks that were delineated from the DOM in the Huojitu Coal Mine Shaft. This DRA-UNet model effectively improves the efficiency of crack delineation. Full article
(This article belongs to the Special Issue Smart Image Recognition and Detection Sensors)
Show Figures

Figure 1

Figure 1
<p>Overview of the study area: (<b>a</b>) location of the study area, (<b>b</b>) UAV mosaiced imagery, (<b>c</b>) imagery of some of the ground surface cracks captured by the UAV, (<b>d</b>) damage to the building from the ground surface crack, and (<b>e</b>) imagery of a ground crack.</p>
Full article ">Figure 2
<p>Construction of the ground crack dataset in the coal mining area.</p>
Full article ">Figure 3
<p>A schematic diagram of the DRA-UNet model.</p>
Full article ">Figure 4
<p>Structure of the residual network module.</p>
Full article ">Figure 5
<p>Structure of the DAM module.</p>
Full article ">Figure 6
<p>Structure of the ASPP module.</p>
Full article ">Figure 7
<p>Crack delineation by different models using the GCCMA-UAV dataset.</p>
Full article ">Figure 8
<p>Crack delineation by different models using the Crack500 dataset.</p>
Full article ">Figure 9
<p>Crack delineation of the coal mining ground surface by UAV images: (<b>a</b>) cracks in the coal mining ground surface, (<b>b</b>) cracks in the mining area, (<b>c</b>) cracks in concrete pavement, (<b>d</b>) cracks in the woods.</p>
Full article ">
18 pages, 7033 KiB  
Article
Pseudo-Spectral Spatial Feature Extraction and Enhanced Fusion Image for Efficient Meter-Sized Lunar Impact Crater Automatic Detection in Digital Orthophoto Map
by Huiwen Liu, Ying-Bo Lu, Li Zhang, Fangchao Liu, You Tian, Hailong Du, Junsheng Yao, Zi Yu, Duyi Li and Xuemai Lin
Sensors 2024, 24(16), 5206; https://doi.org/10.3390/s24165206 - 11 Aug 2024
Viewed by 1487
Abstract
Impact craters are crucial for our understanding of planetary resources, geological ages, and the history of evolution. We designed a novel pseudo-spectral spatial feature extraction and enhanced fusion (PSEF) method with the YOLO network to address the problems encountered during the detection of [...] Read more.
Impact craters are crucial for our understanding of planetary resources, geological ages, and the history of evolution. We designed a novel pseudo-spectral spatial feature extraction and enhanced fusion (PSEF) method with the YOLO network to address the problems encountered during the detection of the numerous and densely distributed meter-sized impact craters on the lunar surface. The illumination incidence edge features, isotropic edge features, and eigen frequency features are extracted by Sobel filtering, LoG filtering, and frequency domain bandpass filtering, respectively. Then, the PSEF images are created by pseudo-spectral spatial techniques to preserve additional details from the original DOM data. Moreover, we conducted experiments using the DES method to optimize the post-processing parameters of the models, thereby determining the parameter ranges for practical deployment. Compared with the Basal model, the PSEF model exhibited superior performance, as indicated by multiple measurement metrics, including the precision, recall, F1-score, mAP, and robustness, etc. Additionally, a statistical analysis of the error metrics of the predicted bounding boxes shows that the PSEF model performance is excellent in predicting the size, shape, and location of impact craters. These advancements offer a more accurate and consistent method to detect the meter-sized craters on planetary surfaces, providing crucial support for the exploration and study of celestial bodies in our solar system. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) An illustration of an impact crater slice. (<b>b</b>) The number–diameter distribution of 34,876 impact craters with diameters less than 14 m in our cropped 59 slices from the CE-4 landing site.</p>
Full article ">Figure 2
<p>Schematic diagram of the PSEF method for detecting meter-sized impact craters. The parameter preceding the ‘@’ represents the number of channels in the image. The parameters <span class="html-italic">W</span> and <span class="html-italic">H</span> following ‘@’ denote the width and height of the image, respectively. In this figure, <span class="html-italic">W</span> and <span class="html-italic">H</span> are both set to 320 pixels, e.g., ‘1@4<span class="html-italic">W</span>×4<span class="html-italic">H</span>’ under the PHS image indicates that the PHS image is single-channel with a size of 1280 × 1280 pixels. The grey arrows indicate that the images connected by these arrows represent the same image at different stages of the process.</p>
Full article ">Figure 3
<p>Illustration of images in the pseudo-spectral space. (<b>a</b>) The original PHS image. (<b>b</b>) Incidence features of crater rims obtained by the Sobel operator along the direction of incident light on the PHS image. (<b>c</b>) Isotropic features of meter-sized craters obtained by the LoG operator on the PHS image. (<b>d</b>) Eigen frequency features of meter-sized craters obtained by the bandpass filter in the frequency domain on the PHS image. (<b>e</b>) The PMS image obtained by amalgamating the images (<b>b</b>–<b>d</b>). (<b>f</b>) The PSEF image obtained by panchromatic sharpening on the PMS and PHS images.</p>
Full article ">Figure 4
<p>Frequency domain maps of PHS images with various mask regions. (<b>a</b>) Only the central 1–11 pixels annular mask region in the frequency domain of the PHS image is allowed to pass through. (<b>b</b>) The pseudo-spectrum obtained after the inverse DFT exhibiting eigen frequency features of impact craters with diameters in the hundreds of meters. (<b>c</b>) The annular mask region of 25–50 pixels in the frequency domain. (<b>d</b>) Eigen frequency features of impact craters with diameters over 20 m. (<b>e</b>) The concentric dual annular masks of 215–270 pixels and 350–400 pixels in the frequency domain. (<b>f</b>) Eigen frequency features for meter-sized impact craters.</p>
Full article ">Figure 5
<p>Illustrations for the optimization of the confidence threshold (<span class="html-italic">σ</span>) and IoU threshold (<span class="html-italic">τ</span>) using the discrete exhaustive search (DES) method. The <span class="html-italic">F</span><sub>1</sub>-score in the best post-processing parameters for the PSEF model and the Basal model are highlighted with blue circles.</p>
Full article ">Figure 6
<p>Visualization of partial test dataset slices and their corresponding automatic identification results. (<b>a</b>–<b>c</b>) are images from the test dataset. (<b>d</b>–<b>f</b>) are the identification results from the PSEF model. (<b>g</b>–<b>i</b>) display the identification results of the Basal model. Among them, the green, red, yellow, and blue circles refer to true positives, ground truth, false negatives, and false positives, respectively.</p>
Full article ">Figure 7
<p>Statistical marginal distribution of the diameter relative error (<span class="html-italic">δ<sub>D</sub></span>), the eccentricity error (Δ<span class="html-italic">e</span>), and the location error (Δ<span class="html-italic">L</span>) between true positives and the corresponding ground truths in the test dataset under the optimal post-processing parameters for the Basal model (<b>a</b>–<b>c</b>) and the PSEF model (<b>d</b>–<b>f</b>). The scatter points represent the distribution of error metrics—diameter for each true positive. The top and bottom of each box plot represent 75% and 25% of the data, respectively. Within the box, the solid line indicates the median, while the dot denotes the mean value. The whiskers extending from the boxes represent the interval within one standard deviation of the mean. Outliers that exceed the whiskers are marked individually.</p>
Full article ">
17 pages, 12277 KiB  
Article
Is Your Training Data Really Ground Truth? A Quality Assessment of Manual Annotation for Individual Tree Crown Delineation
by Janik Steier, Mona Goebel and Dorota Iwaszczuk
Remote Sens. 2024, 16(15), 2786; https://doi.org/10.3390/rs16152786 - 30 Jul 2024
Cited by 2 | Viewed by 941
Abstract
For the accurate and automatic mapping of forest stands based on very-high-resolution satellite imagery and digital orthophotos, precise object detection at the individual tree level is necessary. Currently, supervised deep learning models are primarily applied for this task. To train a reliable model, [...] Read more.
For the accurate and automatic mapping of forest stands based on very-high-resolution satellite imagery and digital orthophotos, precise object detection at the individual tree level is necessary. Currently, supervised deep learning models are primarily applied for this task. To train a reliable model, it is crucial to have an accurate tree crown annotation dataset. The current method of generating these training datasets still relies on manual annotation and labeling. Because of the intricate contours of tree crowns, vegetation density in natural forests and the insufficient ground sampling distance of the imagery, manually generated annotations are error-prone. It is unlikely that the manually delineated tree crowns represent the true conditions on the ground. If these error-prone annotations are used as training data for deep learning models, this may lead to inaccurate mapping results for the models. This study critically validates manual tree crown annotations on two study sites: a forest-like plantation on a cemetery and a natural city forest. The validation is based on tree reference data in the form of an official tree register and tree segments extracted from UAV laser scanning (ULS) data for the quality assessment of a training dataset. The validation results reveal that the manual annotations detect only 37% of the tree crowns in the forest-like plantation area and 10% of the tree crowns in the natural forest correctly. Furthermore, it is frequent for multiple trees to be interpreted in the annotation as a single tree at both study sites. Full article
(This article belongs to the Section AI Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Mapping informal settlements in Nairobi, Kenya with manual annotations. Each colored line indicates a different annotator’s delineation of the same area [<a href="#B16-remotesensing-16-02786" class="html-bibr">16</a>]: (<b>a</b>) boundary deviation due to generalization of informal settlements and (<b>b</b>) deviation resulting from inclusion or exclusion of fringe [<a href="#B26-remotesensing-16-02786" class="html-bibr">26</a>] (adapted from Elemes et al. [<a href="#B16-remotesensing-16-02786" class="html-bibr">16</a>] with permission from Kohli et al. [<a href="#B26-remotesensing-16-02786" class="html-bibr">26</a>]).</p>
Full article ">Figure 2
<p>The four validation areas (red outlines) of study site 1.</p>
Full article ">Figure 3
<p>Nadir 3D point cloud in RGB color scheme (<b>a</b>) and derived 2D segments (<b>b</b>), which represent the single tree reference data for the validation process of study site 2.</p>
Full article ">Figure 4
<p>Example annotation images with 512 × 512 pixel resolution based on the digital orthophoto (<b>a</b>) and the satellite image from WorldView-3 (<b>b</b>).</p>
Full article ">
22 pages, 6799 KiB  
Article
Detection of Cliff Top Erosion Drivers through Machine Learning Algorithms between Portonovo and Trave Cliffs (Ancona, Italy)
by Nicola Fullin, Michele Fraccaroli, Mirko Francioni, Stefano Fabbri, Angelo Ballaera, Paolo Ciavola and Monica Ghirotti
Remote Sens. 2024, 16(14), 2604; https://doi.org/10.3390/rs16142604 - 16 Jul 2024
Viewed by 1231
Abstract
Rocky coastlines are characterised by steep cliffs, which frequently experience a variety of natural processes that often exhibit intricate interdependencies, such as rainfall, ice and water run-off, and marine actions. The advent of high temporal and spatial resolution data, that can be acquired [...] Read more.
Rocky coastlines are characterised by steep cliffs, which frequently experience a variety of natural processes that often exhibit intricate interdependencies, such as rainfall, ice and water run-off, and marine actions. The advent of high temporal and spatial resolution data, that can be acquired through remote sensing and geomatics techniques, has facilitated the safe exploration of otherwise inaccessible areas. The datasets that can be gathered from these techniques, typically combined with data from fieldwork, can subsequently undergo analyses employing/applying machine learning algorithms and/or numerical modeling, in order to identify/discern the predominant influencing factors affecting cliff top erosion. This study focuses on a specific case situated at the Conero promontory of the Adriatic Sea in the Marche region. The research methodology entails several steps. Initially, the morphological, geological and geomechanical characteristics of the areas were determined through unmanned aerial vehicle (UAV) and conventional geological/geomechanical surveys. Subsequently, cliff top retreat was determined within a GIS environment by comparing orthophotos taken in 1978 and 2022 using the DSAS tool (Digital Shoreline Analysis System), highlighting cliff top retreat up to 50 m in some sectors. Further analysis was conducted via the use of two Machine Learning (ML) algorithms, namely Random Forest (RF) and eXtreme Gradient Boosting (XGB). The Mean Decrease in Impurity (MDI) methodology was employed to assess the significance of each factor. Both algorithms yielded congruent results, emphasising that cliff top erosion rates are primarily influenced by slope height. Finally, a validation of the ML algorithm results was conducted using 2D Limit Equilibrium Method (LEM) codes. Ten sections extracted from the sector experiencing the most substantial cliff top retreat, as identified by DSAS, were utilised for 2D LEM analysis. Factor of Safety (FS) values were identified and compared with the cliff height of each section. The results from the 2D LEM analyses corroborated the outputs of the ML algorithms, showing a strong correlation between the slope instability and slope height (R2 of 0.84), with FS decreasing with slope height. Full article
Show Figures

Figure 1

Figure 1
<p>Map displaying the position of the study area: (<b>a</b>) Satellite image showing the study area and the three sectors called Portonovo, Mezzavalle and Trave. (image taken from GeoEye satellite database, 2020). (<b>b</b>) Location of the study area along the Italian Adriatic coast.</p>
Full article ">Figure 2
<p>Workflow sketch: (1) fieldwork; (2) data analyses and surveys. (3) Parameters’ extraction. (4) Machine learning analysis. (5) Slope stability analysis.</p>
Full article ">Figure 3
<p>Representation of Random Forest.</p>
Full article ">Figure 4
<p>Representation of eXtreme Gradient Boosting.</p>
Full article ">Figure 5
<p>Extracted sections are highlighted in violet and numbered. Sections used for the LEM analysis are represented in violet from 1–10.</p>
Full article ">Figure 6
<p>Geological/geomorphological setting of the study area and “Sectors”: (1) Portonovo; (2) Mezzavalle; (3) Trave. Bedrock legend: SCH (Schlier Fm., Lower Miocene-Upper Miocene), GNOa (Sapigno Fm. Upper Miocene), FCO (Colombacci Fm., Upper Miocene), Tv (Trave horizon, Lower Pliocene), FAA (Argille Azzurre Fm., Lower Pliocene-Lower Pleistocene).</p>
Full article ">Figure 7
<p>Pictures portraying the different configurations of the cliff base. (<b>a</b>) Portonovo sector; (<b>b</b>) Mezzavalle sector; (<b>c</b>) Trave sector.</p>
Full article ">Figure 8
<p>Horizontal georeferencing uncertainty, measured between the 1978 and 2022 orthophotos.</p>
Full article ">Figure 9
<p>NSM values calculated along transects, referring to the period 1978–2022: Portonovo and Trave sector showed the highest values of retreat. Computed transects (<b>a</b>) at Portonovo sector; (<b>b</b>) at Mezzavalle sector; (<b>c</b>) at Trave sector.</p>
Full article ">Figure 10
<p>ML results and analysis performed using RF and XGB algorithms. (1) Sketch displaying the parameters used in the analysis. (2) Illustration of the ML algorithm used. (3) Sketch illustrating ML results. (3a) Graph showing feature importance resulting from RF algorithm. (3b) Confusion matrix of RF elaboration. (3c) Graph illustrating the feature importance extracted using XGB algorithm. (3d) Confusion matrix of XGB analysis.</p>
Full article ">Figure 11
<p>Graph comparing FS and cliff height values in the ten extracted sections. The number below the dot refers to the section number in <a href="#remotesensing-16-02604-t005" class="html-table">Table 5</a>.</p>
Full article ">
25 pages, 20831 KiB  
Article
Digital Tools for the Preventive Conservation of Built Heritage: The Church of Santa Ana in Seville
by Estefanía Chaves, Jaime Aguilar, Alberto Barontini, Nuno Mendes and Víctor Compán
Heritage 2024, 7(7), 3470-3494; https://doi.org/10.3390/heritage7070164 - 30 Jun 2024
Cited by 3 | Viewed by 1490
Abstract
Historic Building Information Modelling (HBIM) plays a pivotal role in heritage conservation endeavours, offering a robust framework for digitally documenting existing structures and supporting conservation practices. However, HBIM’s efficacy hinges upon the implementation of case-specific approaches to address the requirements and resources of [...] Read more.
Historic Building Information Modelling (HBIM) plays a pivotal role in heritage conservation endeavours, offering a robust framework for digitally documenting existing structures and supporting conservation practices. However, HBIM’s efficacy hinges upon the implementation of case-specific approaches to address the requirements and resources of each individual asset and context. This paper defines a flexible and generalisable workflow that encompasses various aspects (i.e., documentation, surveying, vulnerability assessment) to support risk-informed decision making in heritage management tailored to the peculiar conservation needs of the structure. This methodology includes an initial investigation covering historical data collection, metric and condition surveys and non-destructive testing. The second stage includes Finite Element Method (FEM) modelling and structural analysis. All data generated and processed are managed in a multi-purpose HBIM model. The methodology is tested on a relevant case study, namely, the church of Santa Ana in Seville, chosen for its historical significance, intricacy and susceptibility to seismic action. The defined level of detail of the HBIM model is sufficient to inform the structural analysis, being balanced by a more accurate representation of the alterations, through linked orthophotos and a comprehensive list of alphanumerical parameters. This ensures an adequate level of information, optimising the trade-off between model complexity, investigation time requirements, computational burden and reliability in the decision-making process. Field testing and FEM analysis provide valuable insight into the main sources of vulnerability in the building, including the connection between the tower and nave and the slenderness of the columns. Full article
(This article belongs to the Special Issue Architectural Heritage Management in Earthquake-Prone Areas)
Show Figures

Figure 1

Figure 1
<p>General view of Seville. Author Ambrosius Brambilla 1585. Santa Ana church is marked with a circle. National library of Spain.</p>
Full article ">Figure 2
<p>Aerial view of the church. Main and Gospel nave façade.</p>
Full article ">Figure 3
<p>Interior of the church, main altar.</p>
Full article ">Figure 4
<p>Evolution of the plans from the original design until the mid-19th century: (<b>a</b>) original plan, (<b>b</b>) mid-16th century, (<b>c</b>) mid-19th century. Grey: original plant, dark red: new modifications, light red: previous modifications. Adapted from [<a href="#B33-heritage-07-00164" class="html-bibr">33</a>].</p>
Full article ">Figure 5
<p>Seismic hazard map of Spain, Seville marked with a blue dot [<a href="#B35-heritage-07-00164" class="html-bibr">35</a>].</p>
Full article ">Figure 6
<p>Photogrammetric model of the exterior of the church.</p>
Full article ">Figure 7
<p>Exterior and interior photogrammetric models’ horizontal cross section.</p>
Full article ">Figure 8
<p>Photogrammetric model of a transversal arch.</p>
Full article ">Figure 9
<p>Interior orthographic photo of the vaults.</p>
Full article ">Figure 10
<p>A 360° view of the church. Second bay of the central nave.</p>
Full article ">Figure 11
<p>Types of masonry in the church: (<b>a</b>) vertical elements; (<b>b</b>) nave vaults; (<b>c</b>) altar vaults.</p>
Full article ">Figure 12
<p>Cracks in the internal face of the perimetral walls of the church: (<b>a</b>) Epistle nave wall; (<b>b</b>) main façade; (<b>c</b>) Gospel nave wall.</p>
Full article ">Figure 13
<p>Damage map of the Epistle nave wall.</p>
Full article ">Figure 14
<p>Damage map of the intrados of the vaults.</p>
Full article ">Figure 15
<p>Transversal section of the photogrammetry model of the arch and the vault. Measurements in meters.</p>
Full article ">Figure 16
<p>Plan view of the accelerometer locations and setups.</p>
Full article ">Figure 17
<p>Six experimental vibration modes of the global model.</p>
Full article ">Figure 18
<p>Five experimental modes of the portico model.</p>
Full article ">Figure 19
<p>BIM parametric objects adopted in the project: examples of native ones in green; in orange, new ones created ad hoc.</p>
Full article ">Figure 20
<p>Comparison of the idealised geometry and the point cloud.</p>
Full article ">Figure 21
<p>Section of the 3D model showing the parametric objects that includes no geometric information of the model.</p>
Full article ">Figure 22
<p>FEM model of Santa Ana church: (<b>a</b>) high-fidelity model; (<b>b</b>) simplified model.</p>
Full article ">Figure 23
<p>Elements of the model involved in the calibration.</p>
Full article ">Figure 24
<p>Comparison of dynamic properties of the experimental (OMA) and numerical model (FEM). Chapels of the numerical model hidden for a clearer visual comparison.</p>
Full article ">Figure 25
<p>Capacity curves obtained from the pushover analysis and location of control points.</p>
Full article ">Figure 26
<p>Tensile principal strains for the maximum displacements plotted in <a href="#heritage-07-00164-f025" class="html-fig">Figure 25</a> for each direction.</p>
Full article ">
18 pages, 4924 KiB  
Article
LOD2-Level+ Low-Rise Building Model Extraction Method for Oblique Photography Data Using U-NET and a Multi-Decision RANSAC Segmentation Algorithm
by Yufeng He, Xiaobian Wu, Weibin Pan, Hui Chen, Songshan Zhou, Shaohua Lei, Xiaoran Gong, Hanzeyu Xu and Yehua Sheng
Remote Sens. 2024, 16(13), 2404; https://doi.org/10.3390/rs16132404 - 30 Jun 2024
Viewed by 1172
Abstract
Oblique photography is a regional digital surface model generation technique that can be widely used for building 3D model construction. However, due to the lack of geometric and semantic information about the building, these models make it difficult to differentiate more detailed components [...] Read more.
Oblique photography is a regional digital surface model generation technique that can be widely used for building 3D model construction. However, due to the lack of geometric and semantic information about the building, these models make it difficult to differentiate more detailed components in the building, such as roofs and balconies. This paper proposes a deep learning-based method (U-NET) for constructing 3D models of low-rise buildings that address the issues. The method ensures complete geometric and semantic information and conforms to the LOD2 level. First, digital orthophotos are used to perform building extraction based on U-NET, and then a contour optimization method based on the main direction of the building and the center of gravity of the contour is used to obtain the regular building contour. Second, the pure building point cloud model representing a single building is extracted from the whole point cloud scene based on the acquired building contour. Finally, the multi-decision RANSAC algorithm is used to segment the building detail point cloud and construct a triangular mesh of building components, followed by a triangular mesh fusion and splicing method to achieve monolithic building components. The paper presents experimental evidence that the building contour extraction algorithm can achieve a 90.3% success rate and that the resulting single building 3D model contains LOD2 building components, which contain detailed geometric and semantic information. Full article
Show Figures

Figure 1

Figure 1
<p>Diagram of the survey area.</p>
Full article ">Figure 2
<p>MVS point cloud and Digital Orthophoto Map (DOM) image in the scene.</p>
Full article ">Figure 3
<p>DOM labeled blocks.</p>
Full article ">Figure 4
<p>Building monolith point cloud extraction flowchart.</p>
Full article ">Figure 5
<p>Statistical weighting method to identify the main direction of the building.</p>
Full article ">Figure 6
<p>Contour line optimization method.</p>
Full article ">Figure 7
<p>Fine building modeling RANSAC segmentation flowchart.</p>
Full article ">Figure 8
<p>Extraction results of this paper’s method.</p>
Full article ">Figure 9
<p>Segmentation result of the single building point cloud.</p>
Full article ">Figure 10
<p>Roof, pat, eaves point cloud segmentation result.</p>
Full article ">Figure 11
<p>Balcony point cloud segmentation results.</p>
Full article ">Figure 12
<p>Façade point cloud segmentation results.</p>
Full article ">Figure 13
<p>Building 3D model construction in the paper.</p>
Full article ">
24 pages, 24274 KiB  
Article
Multi-Platform Integrated Analysis of the Degradation Patterns of Impact Crater Populations on the Lunar Surface
by Meixi Chen, Xinyu Ma, Teng Hu, Zhizhong Kang and Meng Xiao
Remote Sens. 2024, 16(13), 2359; https://doi.org/10.3390/rs16132359 - 27 Jun 2024
Viewed by 869
Abstract
Following the processing of the Chang’e-4 satellite images, Chang’e-4 landing camera images, and Yutu-2 panoramic camera images, data were obtained in a variety of resolutions, including digital elevation model (DEM) and digital orthophoto map (DOM). By determining the morphological parameters, including the depths [...] Read more.
Following the processing of the Chang’e-4 satellite images, Chang’e-4 landing camera images, and Yutu-2 panoramic camera images, data were obtained in a variety of resolutions, including digital elevation model (DEM) and digital orthophoto map (DOM). By determining the morphological parameters, including the depths and diameters of impact craters in the study area, as well as their degradation classes based on surface texture features, we conducted a comprehensive analysis of the morphological parameters and population degradation patterns of impact craters from multiple platforms. The data from three platforms were employed to identify 12,089 impact craters with diameters ranging from 0.1 m to 800.0 m, which were then classified into five degradation classes based on their morphology in the images. This study indicates that as the size of impact craters increases, the population within them experiences a greater degree of degradation. However, the severe degradation of impact craters with diameters of less than 1 m or even 2 m is influenced by the rapid rate of degradation of the crater and the low solidity of the crater lips. The results of the equilibrium state of impact craters indicate that for sub-metre-sized impact craters (with diameters below 2.0 m), it is challenging to reach equilibrium. Furthermore, the smaller the impact crater, the more difficult it is to achieve equilibrium, which is probably the result of simpler generation conditions and the faster degradation of small impact craters. Full article
(This article belongs to the Special Issue Remote Sensing and Photogrammetry Applied to Deep Space Exploration)
Show Figures

Figure 1

Figure 1
<p>Flow chart of impact craters’ data processing.</p>
Full article ">Figure 2
<p>Impact crater identification method. The red outlined areas correspond to impact craters identified.</p>
Full article ">Figure 3
<p>Crater morphological parameter extraction algorithm.</p>
Full article ">Figure 4
<p>PCAM images. (<b>a</b>) DOM image. (<b>b</b>) DEM image.</p>
Full article ">Figure 5
<p>Identified impact craters in DOM images.</p>
Full article ">Figure 6
<p>LCAM images and identified impact craters.</p>
Full article ">Figure 7
<p>Satellite images (M1298916428) and identified impact craters.</p>
Full article ">Figure 8
<p>Histogram of frequencies versus diameters of craters. (<b>a</b>) Results of PCAM image data. (<b>b</b>) Results of LCAM image data. (<b>c</b>) Results of satellite image data.</p>
Full article ">Figure 8 Cont.
<p>Histogram of frequencies versus diameters of craters. (<b>a</b>) Results of PCAM image data. (<b>b</b>) Results of LCAM image data. (<b>c</b>) Results of satellite image data.</p>
Full article ">Figure 9
<p>Histogram of frequency versus depth-to-diameter ratio of craters. (<b>a</b>) Results of PCAM image data. (<b>b</b>) Results of LCAM image data. (<b>c</b>) Results of satellite image data. The red lines indicate the frequency curve of the data.</p>
Full article ">Figure 9 Cont.
<p>Histogram of frequency versus depth-to-diameter ratio of craters. (<b>a</b>) Results of PCAM image data. (<b>b</b>) Results of LCAM image data. (<b>c</b>) Results of satellite image data. The red lines indicate the frequency curve of the data.</p>
Full article ">Figure 10
<p>Trend plot of the diameter versus the depth-to-diameter ratio of craters.</p>
Full article ">Figure 11
<p>Histogram of frequencies versus diameters of all craters.</p>
Full article ">Figure 12
<p>Depth versus diameter values of craters.</p>
Full article ">Figure 13
<p>Histogram of frequency versus depth-to-diameter ratio of craters.</p>
Full article ">Figure 14
<p>Crater degradation class.</p>
Full article ">Figure 15
<p>Histogram of percentage versus diameter and degradation class of craters.</p>
Full article ">Figure 15 Cont.
<p>Histogram of percentage versus diameter and degradation class of craters.</p>
Full article ">Figure 16
<p>Examples of impact craters, including results of morphology and depth-to-diameter ratio in DOM images. The yellow outlined areas correspond to impact craters identified.</p>
Full article ">Figure 17
<p>Craterstats2-program-calculated equilibrium lines. (<b>a</b>) The three areas as a whole. (<b>b</b>) The satellite imagery area. (<b>c</b>) The landing area. (<b>d</b>) The area along the route of the rover. EF refer to the equilibrium function of Hartmann [<a href="#B36-remotesensing-16-02359" class="html-bibr">36</a>]. PF and CF refer to the production function and the chronology function of Neukum et al. [<a href="#B4-remotesensing-16-02359" class="html-bibr">4</a>], respectively.</p>
Full article ">Figure 18
<p>Craterstats2-program-calculated equilibrium lines of the area along the rover. EF refer to the equilibrium function of Hartmann [<a href="#B36-remotesensing-16-02359" class="html-bibr">36</a>]. PF and CF refer to the production function and the chronology function of Neukum et al. [<a href="#B4-remotesensing-16-02359" class="html-bibr">4</a>], respectively.</p>
Full article ">
21 pages, 7548 KiB  
Article
Photogrammetric Measurement of Grassland Fire Spread: Techniques and Challenges with Low-Cost Unmanned Aerial Vehicles
by Marián Marčiš, Marek Fraštia, Tibor Lieskovský, Martin Ambroz and Karol Mikula
Drones 2024, 8(7), 282; https://doi.org/10.3390/drones8070282 - 22 Jun 2024
Viewed by 1204
Abstract
The spread of natural fires is a complex issue, as its mathematical modeling needs to consider many parameters. Therefore, the results of such modeling always need to be validated by comparison with experimental measurements under real-world conditions. Remote sensing with the support of [...] Read more.
The spread of natural fires is a complex issue, as its mathematical modeling needs to consider many parameters. Therefore, the results of such modeling always need to be validated by comparison with experimental measurements under real-world conditions. Remote sensing with the support of satellite or aerial sensors has long been used for this purpose. In this article, we focused on data collection with an unmanned aerial vehicle (UAV), which was used both for creating a digital surface model and for dynamic monitoring of the spread of controlled grassland fires in the visible spectrum. We subsequently tested the impact of various processing settings on the accuracy of the digital elevation model (DEM) and orthophotos, which are commonly used as a basis for analyzing fire spread. For the DEM generated from images taken during the final flight after the fire, deviations did not exceed 0.1 m compared to the reference model from LiDAR. Scale errors in the model with only approximal WGS84 exterior orientation parameters did not exceed a relative accuracy of 1:500, and possible deformations of the DEM up to 0.5 m in height had a minimal impact on determining the rate of fire spread, even with oblique images taken at an angle of 45°. The results of the experiments highlight the advantages of using low-cost SfM photogrammetry and provide an overview of potential issues encountered in measuring and performing photogrammetric processing of fire spread. Full article
(This article belongs to the Special Issue Unconventional Drone-Based Surveying 2nd Edition)
Show Figures

Figure 1

Figure 1
<p>Location of the testing site in central Slovakia (<b>left</b>) and details of the meadow used for the experimental controlled fire in highlighted yellow region with red point corresponding to the displayed ETRS89 coordinates (<b>right</b>) (source: <a href="http://google.com/maps" target="_blank">google.com/maps</a>, accessed on 25 February 2024).</p>
Full article ">Figure 2
<p>Burned control line from the eastern side of the specific location—wider view including numbers of GCPs and a yellow box (<b>left</b>) to which the detailed image from UAV (<b>right</b>) pertains.</p>
Full article ">Figure 3
<p>Flight at the beginning of the fire ignition—camera network configuration (<b>left</b>) and an example of an oblique image (<b>right</b>).</p>
Full article ">Figure 4
<p>Camera network from dynamic monitoring of fire development at a 2 s interval (red images)—the images were aligned towards the original image block from the first flight (blue images).</p>
Full article ">Figure 5
<p>Image block created immediately after the completion of dynamic monitoring (<b>left</b>) and an example of an oblique image from the UAV (<b>right</b>).</p>
Full article ">Figure 6
<p>Graphic representation of the workflow of the basic experiment.</p>
Full article ">Figure 7
<p>Digital surface model from the 3rd aerial survey after the burning of the entire meadow, containing approximately 0.5 million triangles.</p>
Full article ">Figure 8
<p>Overview of static subsets 1–12 within the image block from the 2nd flight (monitoring—red color).</p>
Full article ">Figure 9
<p>Example of orthorectified images at intervals of 60 s, showing the change in camera perspective during dynamic fire monitoring.</p>
Full article ">Figure 10
<p>Result of manual vectorization of fire spread at an 8 s interval. In the background is an orthophoto mosaic of the burned meadow from the last flight.</p>
Full article ">Figure 11
<p>Example of the distribution of tie points used for orienting the selected image. Blue points represent those for which a corresponding pair was found in subsequent images, while gray points represent those for which no match was found.</p>
Full article ">Figure 12
<p>Vertical deviations of the photogrammetric 3D model of the burned meadow from ALS point cloud—GCPs active (<b>left</b>) and inactive (<b>right</b>) during bundle adjustment (BA). Gray color represents values outside the range of ±0.5 m.</p>
Full article ">Figure 13
<p>Confidence of points in the point cloud used for generating the 3D model determined based on image overlap or the number of depth maps used for point reconstruction. Blue represents the highest reliability, and red represents the lowest reliability.</p>
Full article ">Figure 14
<p>Vertical deviations of the photogrammetric 3D model created solely from the side monitoring flight compared to the ALS point cloud. The gray color represents values outside the range of ±1.0 m.</p>
Full article ">Figure 15
<p>Illustration of length errors (yellow values) on the orthophotomosaic georeferenced using approximate onboard coordinates of projection centers in the WGS84 coordinate system. All lengths (green values) between GCPs (white numbers) were shorter than they should be.</p>
Full article ">Figure 16
<p>Visualization of the impact of orthoimages generated at 2 s intervals from vertically displaced 3D models on the accuracy of determining relative changes, from which the rate of fire spread is calculated.</p>
Full article ">Figure 17
<p>Selected UAV image (<b>left</b>) and its orthorectified version (<b>right</b>). The green rectangle in the right image indicates the area near point 204, which is detailed in <a href="#drones-08-00282-f018" class="html-fig">Figure 18</a>.</p>
Full article ">Figure 18
<p>Detail of the orthoimage rectified based on the 3D model from the post-fire flight (<b>left</b>) and the deformed model generated from the side monitoring flight (<b>right</b>). Differences in the position of the green markers against the background texture are visible in the images. Deviations reached a maximum of 0.3 m, which corresponds to the vertical deviation of the model in the corresponding area in <a href="#drones-08-00282-f014" class="html-fig">Figure 14</a>.</p>
Full article ">
18 pages, 10669 KiB  
Article
Accuracy of Determination of Corresponding Points from Available Providers of Spatial Data—A Case Study from Slovakia
by Slavomir Labant, Patrik Petovsky, Pavel Sustek and Lubomir Leicher
Land 2024, 13(6), 875; https://doi.org/10.3390/land13060875 - 18 Jun 2024
Viewed by 1130
Abstract
Mapping the terrain and the Earth’s surface can be performed through non-contact methoYes, that is correct.ds such as laser scanning. This is one of the most dynamic and effective data collection methods. This case study aims to analyze the usability of spatial data [...] Read more.
Mapping the terrain and the Earth’s surface can be performed through non-contact methoYes, that is correct.ds such as laser scanning. This is one of the most dynamic and effective data collection methods. This case study aims to analyze the usability of spatial data from available sources and to choose the appropriate solutions and procedures for processing the point cloud of the area of interest obtained from available web applications. The processing of the point cloud obtained by airborne laser scanning results in digital terrain models created in selected software. The study also included modeling of different types of residential development, and the results were evaluated. Different data sources may have compatibility issues, which means that the position of the same object from different spatial data databases may not be identical. To address this, deviations of the corresponding points were determined from various data sources such as Real Estate Cadaster, ZBGIS Buildings, LiDAR point cloud, orthophoto mosaic, and geodetic measurements. These deviations were analyzed according to their size and orientation, with the average deviations ranging from 0.22 to 0.34 m and standard deviations from 0.11 to 0.20 m. The Real Estate Cadaster was used as the correct basis for comparison. The area of the building was also compared, with the slightest difference being present between the Real Estate Cadaster and geodetic measurement. The difference was zero after rounding the area to whole numbers. The maximum area difference was +5 m2 for ZBGIS Buildings. Full article
(This article belongs to the Special Issue Geospatial Technology for Landscape Design)
Show Figures

Figure 1

Figure 1
<p>Graphical visualization of different measurement methods in the relationship between (<b>a</b>) object size and number of points on object and (<b>b</b>) object size and measurement accuracy.</p>
Full article ">Figure 2
<p>Localization of the selected area of study from Slovakia based on the regions from which data are provided from the ALS: (<b>a</b>) cadastral unit of Horné Srnie village; (<b>b</b>) orthophoto mosaic of the study area; and (<b>c</b>) classified point cloud of the study area.</p>
Full article ">Figure 3
<p>Methodology of the case study of the positional deviations of the corresponding points.</p>
Full article ">Figure 4
<p>Levels of detail (LOD) in modeling 3D objects.</p>
Full article ">Figure 5
<p>Point cloud visualization of the study area: (<b>a</b>) display of colored point cloud according to classification in <a href="#land-13-00875-t003" class="html-table">Table 3</a>; (<b>b</b>) display of point cloud according to reflection intensity; and (<b>c</b>) DTM displayed with Smooth Modeling visualization style.</p>
Full article ">Figure 6
<p>Visualization of the subject area: (<b>a</b>) DTM display with orthophoto mosaic connection and (<b>b</b>) view of DTM with orthophoto mosaic overlay.</p>
Full article ">Figure 7
<p>Display of cropped point cloud of the open pit mine surface: (<b>a</b>) display of point cloud Ground (open pit mine surface—brown color), Unassigned (embankments—white color), and hight noise (loaders—black color); (<b>b</b>) view of the merged point cloud (Ground and Unassigned); and (<b>c</b>) DTM overlaid with orthophoto mosaic in 3D iso view.</p>
Full article ">Figure 8
<p>Comparison of models of the warehouse building with a flat roof: (<b>a</b>) view of the classified point cloud around the building and the wireframe; (<b>b</b>) filling the wireframe without the point cloud; (<b>c</b>) generated building model from the point cloud as DSM; and (<b>d</b>) 3D building model visualized in ZBGIS.</p>
Full article ">Figure 9
<p>Comparison of models of the hipped-roof apartment building: (<b>a</b>) view of the classified point cloud around the building; (<b>b</b>) filling the wireframe without the point cloud; (<b>c</b>) generated building model from the point cloud; and (<b>d</b>) 3D building model visualized in ZBGIS.</p>
Full article ">Figure 10
<p>Positional comparison of buildings on orthophoto mosaic and Real Estate Cadaster: (<b>a</b>) warehouse building with a flat roof and two height levels and (<b>b</b>) apartment building with a hipped roof.</p>
Full article ">Figure 11
<p>Compatibility of corresponding points—corners of the selected building together with a detailed display of their positional deviations: (<b>a</b>) left—points 1 and 6 and (<b>b</b>) right—points 2 and 5.</p>
Full article ">
31 pages, 34074 KiB  
Article
The Generation of High-Resolution Mapping Products for the Lunar South Pole Using Photogrammetry and Photoclinometry
by Pengying Liu, Xun Geng, Tao Li, Jiujiang Zhang, Yuying Wang, Zhen Peng, Yinhui Wang, Xin Ma and Qiudong Wang
Remote Sens. 2024, 16(12), 2097; https://doi.org/10.3390/rs16122097 - 10 Jun 2024
Viewed by 1212
Abstract
High-resolution and high-accuracy mapping products of the Lunar South Pole (LSP) will play a vital role in future lunar exploration missions. Existing lunar global mapping products cannot meet the needs of engineering tasks, such as landing site selection and rover trajectory planning, at [...] Read more.
High-resolution and high-accuracy mapping products of the Lunar South Pole (LSP) will play a vital role in future lunar exploration missions. Existing lunar global mapping products cannot meet the needs of engineering tasks, such as landing site selection and rover trajectory planning, at the LSP. The Lunar Reconnaissance Orbiter (LRO)’s narrow-angle camera (NAC) can acquire submeter images and has returned a large amount of data covering the LSP. In this study, we combine stereo-photogrammetry and photoclinometry to generate high-resolution digital orthophoto maps (DOMs) and digital elevation models (DEMs) using LRO NAC images for a candidate landing site at the LSP. The special illumination and landscape characteristics of the LSP make the derivation of high-accuracy mapping products from orbiter images extremely difficult. We proposed an easy-to-implement shadow recognition and contrast stretching method based on the histograms of the LRO NAC images, which is beneficial for photogrammetric and photoclinometry processing. In order to automatically generate tie points, we designed an image matching method considering LRO NAC images’ features of long strips and large data volumes. The terrain and smoothness constraints were introduced into the cost function of photoclinometry adjustment, excluding pixels in shadow areas. We used 61 LRO NAC images to generate mapping products covering an area of 400 km2. The spatial resolution of the generated DOMs was 1 m/pixel, and the grid spacing of the derived DEMs was 1 m (close to the spatial resolution of the original images). The generated DOMs achieved a relative accuracy of better than 1 pixel. The geometric accuracy of the DEM derived from photoclinometry was consistent with the lunar orbiter laser altimeter (LOLA) DEM with a root mean square error of 0.97 m and an average error of 0.17 m. Full article
(This article belongs to the Special Issue Remote Sensing and Photogrammetry Applied to Deep Space Exploration)
Show Figures

Figure 1

Figure 1
<p>Flowchart of the geometric processing for the LRO NAC images.</p>
Full article ">Figure 2
<p>Contrast stretching based on the statistics of the image histogram for a typical LRO NAC image.</p>
Full article ">Figure 3
<p>Shadow extraction and contrast stretching results for M174930921RE. (<b>a</b>) Contrast stretching based on the original minimum and maximum gray values. (<b>b</b>) Contrast stretching based on the proposed method. (<b>c</b>) Shadow recognition results based on the proposed method. Red areas in the right figure indicate the shadows.</p>
Full article ">Figure 4
<p>Illustration of the proposed image matching method for automatic tie point extraction.</p>
Full article ">Figure 5
<p>Comparison results of image matching between raw images and approximate orthophotos for a stereo pair composed of M174984904LE and M174971411RE. n means the number of features used for matching.</p>
Full article ">Figure 6
<p>Flowchart of the generation of tie points for LRO NAC images.</p>
Full article ">Figure 7
<p>Flowchart of SPC reconstruction based on the photogrammetric processing results.</p>
Full article ">Figure 8
<p>Coverage of LRO NAC images. The red rectangle indicates the candidate landing area.</p>
Full article ">Figure 9
<p>Observation and illumination conditions of raw LRO NAC images.</p>
Full article ">Figure 10
<p>Distribution of the tie points on the images of local areas. The green crosses represent matched tie points.</p>
Full article ">Figure 11
<p>Overall distribution of tie points in the images. The purple crosses represent the absolute elevation control points acquired from the LOLA DEM, and the green crosses represent matched tie points. Overlapping areas of images without tie points are shadowed areas.</p>
Full article ">Figure 12
<p>Results of the evaluation of geometric accuracy of the bundle adjustment.</p>
Full article ">Figure 13
<p>Residual vector maps in the image space of the tie points. (<b>a</b>,<b>c</b>) are residual vectors, and (<b>b</b>,<b>d</b>) are the corresponding original image tie points.</p>
Full article ">Figure 14
<p>Residual vector maps of elevation control points. The explanation of the symbols in the upper left corner is a quantitative description of the magnitude corresponding to the length of the arrow.</p>
Full article ">Figure 15
<p>DOM mosaic effect. The green rectangle indicates the candidate landing area.</p>
Full article ">Figure 16
<p>Initial DEM chips for the SPC terrain reconstruction. The spatial resolution of the DEM was 1 m/pixel. The size of each DEM chip was 2 km by 2 km, with an overlap of 200 m between adjacent chips. The blue rectangle indicates the candidate landing area.</p>
Full article ">Figure 17
<p>Color-shaded map of the generated DEM mosaic (1 m/pixel). The blue rectangle indicates the candidate landing area.</p>
Full article ">Figure 18
<p>Analysis of the illumination conditions of the available images. (<b>a</b>) Average solar illumination of SPC-generated DEM area. (<b>b</b>) Valid pixels in the candidate landing site. Pixels displayed in black areas have not been illuminated. The blue rectangle indicates the candidate landing area.</p>
Full article ">Figure 19
<p>Qualitative evaluation of the generated DEM. (<b>a</b>) Comparison of the SPC-generated DEM (left side) and the DOM of an LRO NAC image (right side). The hill-shaded DEM was drawn on top of the DOM of the LRO NAC image (M105925266LE). The gray rectangle on the right side indicates the same area corresponding to the generated DEM on the left side. Red rectangles indicate distorted terrain areas. (<b>b</b>) Comparison between the SPC-generated DEM (1 m/pixel) and resampled LOLA DEM (1 m/pixel). The left side shows the hill-shaded map of the SPC-generated DEM, and the right side shows the hill-shaded map of the LOLA DEM.</p>
Full article ">Figure 20
<p>Quantitative evaluation results of the height differences in the overlapping areas of adjacent DEM chips. (<b>a</b>) Positions of the sample points in the overlapping area of two adjacent DEM chips. The red line indicates the sample points for calculating the height differences. (<b>b</b>) Elevation profiles of the adjacent DEM chips. (<b>c</b>) Illustration of the height differences between the adjacent DEM chips.</p>
Full article ">Figure 21
<p>Results of a comparison between the SPC-generated DEM and LOLA DEM. The red line in the top-right subfigure indicates the location of the sample points.</p>
Full article ">Figure 22
<p>Pixel-by-pixel height difference map between SPC-generated DEM and LOLA DEM.</p>
Full article ">
14 pages, 21676 KiB  
Technical Note
A Catalogue of Impact Craters and Surface Age Analysis in the Chang’e-6 Landing Area
by Yexin Wang, Jing Nan, Chenxu Zhao, Bin Xie, Sheng Gou, Zongyu Yue, Kaichang Di, Hong Zhang, Xiangjin Deng and Shujuan Sun
Remote Sens. 2024, 16(11), 2014; https://doi.org/10.3390/rs16112014 - 4 Jun 2024
Cited by 2 | Viewed by 1518
Abstract
Chang’e-6 (CE-6) is the first sample-return mission from the lunar farside and will be launched in May of 2024. The landing area is in the south of the Apollo basin inside the South Pole Aitken basin. Statistics and analyses of impact craters in [...] Read more.
Chang’e-6 (CE-6) is the first sample-return mission from the lunar farside and will be launched in May of 2024. The landing area is in the south of the Apollo basin inside the South Pole Aitken basin. Statistics and analyses of impact craters in the landing area are essential to support safe landing and geologic studies. In particular, the crater size–frequency distribution information of the landing area is critical to understanding the provenance of the CE-6 lunar samples to be returned and can be used to verify and refine the lunar chronology model by combining with the radioisotope ages of the relevant samples. In this research, a digital orthophoto map (DOM) mosaic with resolution of 3 m/pixel of the CE-6 landing area was generated from the 743 Narrow Angle Camera of the Lunar Reconnaissance Orbiter Camera. Based on the DOM, craters were extracted by an automated method and checked manually. A total of 770,731 craters were extracted in the whole area of 246 km × 135 km, 511,484 craters of which were within the mare area. Systematic analyses of the crater distribution, completeness, spatial density, and depth-to-diameter ratio were conducted. Geologic model age estimation was carried out in the mare area that was divided into three geologic units according to the TiO2 abundance. The result showed that the east part of the mare had the oldest model age of μ3.270.045+0.036 Ga, and the middle part of the mare had the youngest model age of μ2.490.073+0.072 Ga. The crater catalogue and the surface model age analysis results were used to support topographic and geologic analyses of the pre-selected landing area of the CE-6 mission before the launch and will contribute to further scientific researches after the lunar samples are returned to Earth. Full article
(This article belongs to the Special Issue Planetary Geologic Mapping and Remote Sensing (Second Edition))
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The base map of the CE-6 landing area generated from LROC NAC images with a pixel size of 3 m. The Lambert conformal conic projection was adopted.</p>
Full article ">Figure 2
<p>Comparison between the diameters automatically measured from extracted craters and diameters of the manually confirmed craters.</p>
Full article ">Figure 3
<p>Crater dating areas (i.e., A1, A2, and A3) of the basalt inside the Apollo basin. The basemap is the TiO<sub>2</sub> abundance [<a href="#B30-remotesensing-16-02014" class="html-bibr">30</a>] overlying the LROC Wide Angle Camera global mosaic [<a href="#B32-remotesensing-16-02014" class="html-bibr">32</a>]. The black polygon is the mare extent defined by [<a href="#B31-remotesensing-16-02014" class="html-bibr">31</a>]. The red, green, and magenta polygons are the crater counting areas of A1, A2, and A3. The black patches inside A2 and A3 are regions excluded for counting due to contamination of secondary craters.</p>
Full article ">Figure 4
<p>The (<b>a</b>) incremental and (<b>b</b>) cumulative size–frequency distributions of craters in the crater catalogue with the diameter internal of <math display="inline"><semantics> <mrow> <msqrt> <mn>2</mn> </msqrt> </mrow> </semantics></math>D in a log-log plot. Note that the diameter axis is the crater median in each bin. (<b>c</b>) The incremental size–frequency distribution established by robust kernel density estimation [<a href="#B33-remotesensing-16-02014" class="html-bibr">33</a>].</p>
Full article ">Figure 5
<p>The mapped craters annotated in red with diameters larger than 200 m in the CE-6 landing area overlaying on the LROC NAC DOM mosaic.</p>
Full article ">Figure 6
<p>Spatial densities of craters in the CE-6 landing area: (<b>a</b>) spatial density of the crater diameters within 30 m to 1 km excluding possible secondary craters and (<b>b</b>) spatial density of craters with D ≥ 1 km excluding possible secondary craters. The Lambert conformal conic projection is adopted. Note that the colors represent different density values in different plots.</p>
Full article ">Figure 7
<p>Crater depth (D &gt; 480 m) distribution of CE-6 landing area in a log-log plot. There are 1547 craters catalogued with a depth interval of 50 m.</p>
Full article ">Figure 8
<p>The relationships between the crater diameter and (<b>a</b>) the crater depth (in a log-log plot) and (<b>b</b>) the d-D ratio. The results for mare craters are shown in red, while the results for highland are shown in black.</p>
Full article ">Figure 9
<p>AMAs of the crater dating areas: (<b>a</b>) A1; (b) A2; (<b>c</b>) A3; (<b>d</b>) A1 + A2. The standard lunar crater equilibrium line is from [<a href="#B37-remotesensing-16-02014" class="html-bibr">37</a>], the PF is from [<a href="#B38-remotesensing-16-02014" class="html-bibr">38</a>], and the CF is from [<a href="#B8-remotesensing-16-02014" class="html-bibr">8</a>].</p>
Full article ">Figure 10
<p>AMAs of the crater dating areas: (<b>a</b>) A1; (<b>b</b>) A2; (<b>c</b>) A3; (<b>d</b>) A1 + A2. The standard lunar crater equilibrium line is from [<a href="#B37-remotesensing-16-02014" class="html-bibr">37</a>], the PF is from [<a href="#B5-remotesensing-16-02014" class="html-bibr">5</a>], and the CF is from [<a href="#B5-remotesensing-16-02014" class="html-bibr">5</a>].</p>
Full article ">Figure A1
<p>Different depth-to-diameter ratios of craters at various degradation stages: (<b>a</b>) a fresh crater (location: 156.07°W, 41.59°S) with diameter of 1769.7 m and depth-to-diameter ratio of 0.177, (<b>b</b>) a degraded crater (location: 152.80°W, 40.58°S) with diameter of 453.2 m and depth-to-diameter ratio of 0.07, and (<b>c</b>) a heavily degraded crater (location: 158.02°W, 41.59°S) with diameter of 364 m and depth-to-diameter ratio of 0.023.</p>
Full article ">
Back to TopTop