1. Introduction
One of the most readily apparent features of the Moon’s surface is the presence of impact craters [
1]. The study of impact craters is of significant importance for the understanding of the lunar morphology, the rock structure of the Moon and the origin of the Moon [
2,
3,
4]. The study of impact craters can be traced back to many years ago. In recent years, the scientific study of impact craters has been greatly facilitated by the ongoing lunar and deep space exploration missions [
5]. A significant corpus of high-resolution remote sensing images of the lunar surface as well as other classes of scientific data have been and continue to be acquired, a large number of analyses of the morphological characteristics of impact craters have been conducted, and a database of impact craters of different sizes on the surfaces of Earth-like planets has been established [
6]. Furthermore, researchers conducted an analysis of the impact crater populations, examining the degradation patterns of impact crater populations based on their degraded morphological features. The research is particularly concerned with the current production and degradation processes on the lunar surface, taking into account the characteristics of small young lunar impact craters [
7]. This is significant for speculating the age of the geological units on the lunar surface as well as the evolution of the structure of the surface materials [
8].
Previous studies of impact crater degradation were primarily based on images obtained from the lunar exploration satellites of the Chang’e project. The successful launch of Chang’e-4 has enabled the landing camera of the lander and the panoramic camera of the Yutu-2 rover to provide a substantial number of high-resolution images for in situ exploration and roving surveys of the lunar surface [
9]. Consequently, the image data from disparate platforms enables the examination of the deterioration patterns of the impact crater populations from a multitude of vantages and in a comprehensive manner.
A substantial body of research has been conducted on the degradation patterns of impact craters. Currently, the method of impact craters’ morphology–age relationship analysis is the principal research method employed by the majority of researchers. For instance, Craddock and Haruyama et al. investigated the process of impact crater degradation and proposed a method for determining the relative age of geological units based on the concept of a degradation parameter [
10,
11]; Wood, Ivanov et al. classified the degradation through the impact craters’ morphology [
12,
13]; and Craddock established a degradation model and estimated the degradation and erosion rate of impact craters by using Clementine data [
10]. The saturated equilibrium state of impact crater density has been studied; it is the state in which a crater population smaller than a given diameter is considered to be in equilibrium when the production rate of craters is balanced by the removal rate. For example, Hartman and James studied cratering equilibrium using computer modelling methods [
14,
15]. Cratering equilibrium is a crater population with a diameter smaller than a given diameter when the production rate of craters is balanced by the removal rate. Statistical studies of the morphology of degradation classes and other relevant features are crucial for elucidating the correlation between the morphological evolution of impact craters and the age of the corresponding lunar surface in the area. The classification of impact crater degradation has been predominantly based on visual interpretation. In 1963, Arthur proposed a classification of four degradation classes for the Moon. Class C1 comprises fresh impact craters, Class C2 initially degraded impact craters with smooth rims, Class C3 impact craters with severely eroded morphology, and Class C4 completely flattened impact craters. However, distinguishing the four degradation classes through visual interpretation may be much affected by subjective judgement and by crater size [
16]. In 1976, five degradation classes were proposed for the Moon based on the study of McGill and Wise, Basilevsky et al. These comprised three main classes (A, B, and C) and two transitional classes (AB and BC) [
17]. A refers to the fresh impact craters, Class B refers to the initially degraded impact craters (with smooth edges), and Class C refers to the completely flattened impact craters. The distinction between the five degradation classes may be significantly influenced by subjective judgement.
The current research on the degradation patterns of impact crater populations on the lunar surface is primarily conducted using the scientific data provided by the CCD stereo camera, which is equipped on the lunar exploration satellite; the scientific data provided by the landing camera, which is equipped on the lander; and the scientific data provided by the panoramic camera, which is equipped on the rover. In previous studies, satellite data were primarily employed to investigate kilometre- and sub-kilometre-sized impact craters on the lunar surface. Some researchers have demonstrated that the depth-to-diameter ratio of fresh simple craters on the Moon is approximately 0.2. Nevertheless, for fresh simple impact craters with a diameter of less than 1 km, the depth-to-diameter ratio is only 0.13, which may be related to the strength of the shallow lunar surface layer [
18,
19,
20,
21]. The database of Losiak et al. indicates that the depth-to-diameter ratio of impact craters with diameters between 1 m and 15 km is approximately 0.17 [
22]. The observed low value of the ratio is likely the result of the degradation of some impact craters. The lunar frontal images captured by the Apollo 11 rover demonstrate the presence of a considerable number of fresh impact craters and severely degraded impact craters in the vicinity of the Apollo 11 landing site. Some researchers also compared the degradation patterns of impact craters on the Apollo 11 landing site with those in the Chang’e-4 landing area. Their findings indicated that the degradation of impact craters initially occurred at a rapid pace, but the degradation rate subsequently decreased significantly as the degree of degradation increased [
23].
In addition, the researchers who created the images from the rover also invested a great deal of effort into the project. Yang et al. employed the panoramic camera of the Yutu-2 rover to analyse the statistical pattern of secondary impact craters in the landing zone. The study demonstrated that the majority of impact craters exhibited a depth-to-diameter ratio within the range of 0.10 to 0.16. An analysis of the size–frequency distribution (SFD) of impact craters reveals that impact craters with diameters of tens of metres are typically in a saturated state, exhibiting a range of degradation patterns [
24]. Hu et al. validated the degradation classes of impact craters by utilising the ground-based topographic camera and the panoramic camera on the Yutu-2, with a particular focus on analysing centimetre to metre impact craters, and assessed the equilibrium state of the centimetre to metre impact crater populations in the Chang’e-4 landing area. The conclusion was that more than 90% of the impact craters have moderately or severely degraded, and as with the degradation pattern of large impact crater populations, the smaller the size of the crater, the severer the degradation is. Moreover, Hu et al. demonstrated that it is more challenging to achieve equilibrium for small impact craters in 2022, with cumulative SFD slopes of −2.696 and −2.392 observed in these two areas. This may be attributed to the presence of secondary craters and/or resurfacing events resulting from ejecta from adjacent craters [
25].
Existing studies generally agree that for impact craters of the same size, the earlier they were formed, i.e., the older they are, the more obvious their desegregation is, and the smaller the impact craters are, the likelier they are to disappear faster as a result of degradation. However, few studies have integrated satellite image data, landing camera image data, and rover image data for analysis. Consequently, in this study, we integrated three distinct datasets from disparate platforms with the objective of analysing the statistical degradation patterns of the impact crater populations from multiple platforms. This approach enabled us to achieve a comprehensive evaluation of the impact craters in the area, a feat that had not been accomplished in previous studies.
Section 2 mainly presents the experimental data and data processing methods employed in this paper.
Section 3 presents the manual extraction results and classification results of impact craters from the data respectively from the satellite, landing camera, and rover. This includes a separate and comprehensive analysis of the outcomes from the three platforms.
Section 4 compares the data from multiple perspectives and discusses the overall degradation and saturation status.
2. Data and Methodology
2.1. Methodology
This paper explores the distribution and degradation patterns of multi-platform impact craters by using downloaded image data from the satellite and the landing camera, as well as the processed image data from the panoramic camera. The image data include the digital elevation model (DEM) and the digital orthophoto map (DOM). Based on the image data, all of the recognizable impact craters in the images were manually extracted, and the degradation classification was carried out by taking into account the spatial morphologies of the impact craters. The specific process is shown in
Figure 1. The detailed description of the method is as below.
2.2. Data and Processing
The image data are from Chang’e-4 satellite, Chang’e-4 landing camera, and Yutu-2 panoramic camera. Among them, the image data from Chang’e-4 satellite were downloaded from the SELENE (KAGUYA) dataset [
26]. This DOM image’s (M1298916428) resolution is 0.8 m. The data from Chang’e-4 landing camera were produced on the basis of the Chang’e-4 landing camera dataset. The ground resolution of the images was heterogeneous, the resolution of the landing camera increased from 0.09 m/pixel to 0.05 m/pixel when the camera descended along the main optical axis, and the georeferenced coordinates (GCS_Moon_2000) were provided. The image data from the panoramic camera were produced on the basis of the Yutu-2 panoramic camera dataset, and the georeferenced coordinates (GCS_Moon_2000) were provided [
27].
The panoramic camera has two modes—colour imaging and panchromatic imaging, which can be switched to each other during the shooting process. The image resolution of the colour mode is 2352 × 1728 pixels, while the image resolution of the panchromatic imaging mode is 1176 × 864 pixels, and the camera field-of-view angle is 19.6° × 14.5° [
9]. The panoramic camera relies on the up and down, left and right rotation of the mast to capture 360° images around the station, and 56 images can be captured in each turn, including 28 pairs of images in different directions at two altitude angles, thus realizing the three-dimensional stereo imaging of the topography and landscape of the landing area and the roving area [
28]. By 20 September 2023, Chang’e-4 had completed its 62nd lunar day of work, and the data used in this paper are the images taken by the panoramic camera and produced from 1st to 40th lunar day.
The processing of images from the panoramic camera was divided into three steps. Firstly, the images were classified according to the time of production, and were divided into 66 stations, with each station including 122 image data files. The data from one station refers to the image data obtained from the Yutu-2 rover after it ceased forward motion, remained stationary, and was photographed by rotating the panoramic camera in the four cardinal directions. Images from each station were imported into MetashapePro1.7.4, where camera calibration and zooming were performed, and a dense point cloud was created. The focal length of the panchromatic camera is 50.29 mm, with the pixel being 0.0148 mm; the focal length of the colour camera is 50.274 mm, with the pixel being 0.0074 mm, and the zooming was based on the track width of Yutu-2—15 cm. Finally, through the dense point cloud, the grid image, the DEM image, and the DOM image were produced and the 3D model and corresponding texture were obtained at the same time. In particular, as the PCAM of Yutu-2 can achieve clear imaging at an infinite distance from the detector foot to the field of view, the DEM and DOM resolutions are variable, ranging from 0.0011–0.004 m/pixel.
2.3. Identification Methods
Impact craters are depressed and circular, so they can be identified through the contrast between light and dark under lighting conditions. In this paper, CraterTools [
29] was used to identify and mark the impact craters in the generated DOM and DEM images.
In the identification process, we adopted a multi-source data fusion method, as shown in
Figure 2. We manually identified the impact craters in the DOM and DEM images by identifying the depressed circles and the contrast of illuminated and shadowed areas. In the manual identification of impact craters, the approximate location of the impact crater is extracted, and it is not required to ensure that the identification results fit the centre and edges of the impact crater. The exact location of the impact crater and the extraction of morphological parameters were obtained based on an extraction algorithm, which results in improved extraction speed and accuracy. In addition, it is evident from
Figure 2b,d,f that we created full-map contour lines on the basis of the DEM images to assist the identification of impact craters in the DOM images, where the interval between contour lines is 0.01 m. The purpose of creating contour lines is to better identify the complete morphology of the craters, and to more clearly find the impact craters that do not have obvious changes of depth in the DEM and DOM images. According to
Figure 2a,b, image data from the satellite are mainly based on the impact craters identified through the DOM images and contour information established using DEM. As shown in
Figure 2c,d, the data of impact craters from the landing camera (LCAM) were mainly identified by means of the 9 cm DOM images and the contour lines generated from the DEM images, similar to impact crater identification in satellite images. The identification of impact craters’ images in the panoramic camera (PCAM) was more complicated. Due to the limitation of the imaging angle of the panoramic camera, the bottom of some impact craters failed to be photographed, so these incomplete images had to be abandoned. The metre-sized impact craters are small and shallow, without obvious contrast between illuminated and shadowed areas; as in
Figure 2e,f, we used the contour based on the DEM images and DOM images generated to identify the impact craters, in which the elevation interval between contour lines is 0.01 m. As the completeness of the fit to the centre of the craters could not be clearly observed in the DOM images generated by the panoramic camera of the rover, when identifying these impact craters, we compared the DOM images with the original images taken by the rover to ensure the morphological completeness of such impact craters.
Figure 2g illustrates the original images from the rover which correspond to
Figure 2e.
2.4. Morphological Parameter Extraction
In this paper, an automatic extraction algorithm [
30] involved the use of a technique designed to extract multi-morphological parameters from the impact craters, and the algorithm includes three steps, as shown in
Figure 3. The algorithm firstly corrected the original data using the method of circle fitting, including the correction of the centroid offset and the position of the impact crater edge. The identification result of the offset is closer to the real position of the impact crater, so that the more accurate position and parameter information of the impact crater can be obtained and improve the data quality.
Subsequently, the least squares algorithm was employed to fit high-order polynomials to the profile lines of the impact crater. This resulted in the acquisition of high-order polynomial functions with a superior degree of fit, which were then used to obtain the fitted profile lines and the polynomial model. Finally, the morphology of the impact crater profiles was described from a mathematical perspective. Finally, the automatic extraction algorithm, combined with the morphological method to analyse the impact crater profile, designed eight morphological descriptors, including single morphological parameters—depth, internal and external slope, pit lip width, pit lip height, etc.—and composite morphological parameters—roundness, depth-to-diameter ratio, and body shape ratio. The automatic extraction algorithm was designed to implement a range of extraction algorithms, each tailored to extract a specific morphological parameter. This paper primarily examines the morphological parameters, including impact crater diameter, depth, and depth-to-diameter ratio. These parameters offer insights into the shape, size, and evolution of impact craters.
2.5. Degradation Classification
According to the study of Basilevsky et al. [
17], impact craters were classified into five classes—A, AB, B, BC, and C—according to their different degrees of degradation. Among them, Class A refers to the freshest impact craters, and Class C refers to the mostly degraded impact craters.
Table 1 illustrates examples of different types of impact craters on different platforms: Class A has obvious light or clear edges, deep crater, and sharp contrast between light and dark; Class AB has slight degradation with obvious contrast between light and dark; Class B has clear and intact edges, degradation traces, and obvious shadows; Class BC has incomplete edges, obvious degradation traces, and shallow craters with unobvious contrast between illuminated and shadowed areas; and Class C has blurred and incomplete edges, severe degradation, shallow crater, no shadows, and a circular outline that can be distinguished when the image is stretched. In addition, satellite image, LCAM image, and PCAM image were all taken on a lunar day with better illumination, which resulted in fewer impact crater shadows because of the lighting conditions and less impact on degradation classification.
In this paper, the degradation classes of impact craters were identified manually and qualitatively. Prior to the degradation classification, our researchers performed an overall identification of the impact craters in the different resolution images, for example, by judging the range of brightness and grey scale values of the images. The results of the overall identification are taken into account for impact crater classification, which reduces the impact of image resolution and solar elevation on the classification results. Manual and qualitative classification is usually affected by subjective factors, so judgement mistakes are inevitable. For instance, Class BC in
Table 1 is analogous to Class C, which may result in misclassification. Consequently, a multiple voting procedure was applied to all impact craters. In the initial phase of the voting process, all impact craters were classified as degraded by three researchers. Subsequently, the degradation classes were collated. For any given impact crater, the degradation class would be confirmed only if the number of votes was 60% or more. In the event that the impact crater was not classified with a specific degradation class, it would be considered to have an uncertain degradation class and would require reclassification. Finally, all impact crater degradation classes received 60% or more support, thus allowing for the classification of the impact crater degradation to be completed.
3. Results
Based on the above research methods, we manually identified and extracted the impact craters’ morphological parameters in the Chang’e-4 satellite image, Chang’e-4 LCAM image, and Yutu-2 PCAM images, and conducted statistics and analyses on their degradation degree. The completeness of the impact crater population statistics was mainly subject to the image resolution and illumination conditions, and the impact craters with a diameter longer than 10× pixels had a higher confidence coefficient [
31]. In order to provide a robust statistical analysis of crater degradation, it was deemed appropriate to consider 95% of the craters in the centre of each area to have a high confidence coefficient. All impact craters within the 95% confidence interval were subsequently presented in the results section. A total of 12,089 impact craters were manually extracted and located within the confidence interval. The minimum diameter of the impact crater identified in the panoramic camera image was 0.067m, the minimum diameter of the impact crater identified in the landing camera image was 1.55m, and the minimum diameter of the impact crater identified in the satellite imagery was 56.04 m, which are all greater than 10× pixels.
3.1. Extraction
In order to extract the impact craters from the images captured by the Yutu-2 rover panoramic camera, we reconstructed the scenes of 66 stations and generated DEM images and DOM images, as illustrated in
Figure 4. The results demonstrated the high image resolution and high completeness of the impact craters, which provided a robust foundation for the study of impact craters along the trajectory of the Yutu-2 rover.
The digital orthophotos of each station were employed as the fundamental study area, which was a circular area with an average outer diameter of 11 m and an inner diameter of 3 m, situated at the centre of the rover’s position, with an average area of 297 m
2.
Figure 5 illustrates the distribution of impact craters in the DOM image. It can be observed that there is a radial distribution around the rover, with the density of impact craters decreasing with increasing distance from the rover. A small number of complete impact craters can be seen to be scattered outside the DOM image. In general, the impact craters are numerous, relatively small in size, and distributed in close proximity to the centre of the rover. The distribution pattern of all impact craters extracted from all stations is essentially consistent with that shown in
Figure 5.
The extraction of the impact craters from the images captured by the landing camera of Chang’e-4 is illustrated in
Figure 6. This figure also presents the image range and the extracted results. The total imaging area is 25,868 m
2, and the imaging range includes some images captured by the panoramic camera.
Figure 7 illustrates the image range and extraction results for the extraction of impact craters from Chang’e-4 satellite images. The total imaging area is 217.4965 km
2, with a strip-like shape. The imaging range encompasses the images captured by the landing camera and the panoramic camera.
The images were analysed to determine the number of impact craters visible in them. This was achieved by utilizing the panoramic camera of the rover. A total of 11,070 impact craters were counted, and the histogram of the size frequency distribution is presented in
Figure 8a. Among these impact craters, the smallest diameter is 0.1005 m, the largest diameter is 6.376 m, and 70 craters have a diameter of more than 2.0 m. The number of impact craters with diameters of 0.2 m to 0.3 m was counted. The impact craters with diameters between 0.2 m and 0.3 m are the greatest in quantity at 2736, followed by those with diameters between 0.3 m and 0.4 m, with a quantity of 2434. A review of the available data indicates that more than 80% of the impact craters have a diameter between 0.2 and 0.9 m. Among these, 61.4% have a diameter between 0.2 and 0.5 m. The impact craters extracted from the images captured with the panoramic camera are predominantly sub-metre sized, with a high degree of uniformity in size range.
The images captured by the Chang’e-4 landing camera were analysed to determine the number of impact craters present. A total of 96 impact craters were identified, and the histogram of the size frequency distribution is presented in
Figure 8b. For these craters, the smallest diameter is 1.554 m, and the largest diameter is 50.870 m. A total of 64 impact craters have diameters between 2.5 m and 7.5 m, and more than 80% of the impact craters have diameters between 2.5 m and 15.0 m. It can be seen that the impact craters extracted from the images photographed with the landing camera are mainly metre-sized impact craters. Due to the limitation of image resolution, the extracted impact craters are small in number but similar in size range.
We counted 923 impact craters extracted from the satellite images, and the histogram of the size frequency distribution is shown in
Figure 8c. For these craters, the smallest diameter is 56.04 m, the largest diameter is 791.975 m, and more than 80% of the craters have diameters between 80.0 m and 320.0 m. The impact craters with a diameter between 120 and 160 m have the highest frequency (218), followed by those with a diameter between 80 and 120 m (152). From the data presented in the figure, it can be observed that in contrast to the smaller impact craters, those with a diameter greater than 100 m do not exhibit a similar size range. This may be influenced by the image range and resolution.
3.2. Extraction Results of Morphological Parameters
A frequency analysis was conducted on the depth-to-diameter ratio (h/D) of impact craters extracted from images captured by the satellite, landing camera, and the panoramic camera of the rover, as illustrated in
Figure 9. It can be observed that the depth-to-diameter ratios of the impact craters with varying diameters from different platforms exhibit a normal distribution. The depth-to-diameter ratios of the impact craters extracted from the images captured by the panoramic camera ranged from 0.059 to 0.261. Those extracted from the images captured by the landing camera exhibited a similar range, from 0.054 to 0.249. The depth-to-diameter ratios of the impact craters extracted from the satellite images spanned a narrower range, from 0.022 to 0.19. The mean values are 0.122, 0.115, and 0.079, respectively. In over 80% of the total impact craters photographed through platforms, the aforementioned values are 0.07–0.19, 0.04–0.16, and 0.01–0.13. Both the information in the figure and the mean values of the depth-to-diameter ratios can conclude that the average depth-to-diameter ratio of the impact craters extracted from satellite images is the smallest at 0.079, that of the landing camera is the second largest, and that of the panoramic camera is the largest at 0.122. We can conclude, based on the information in
Figure 8, that 80% of the impact craters extracted from satellite images are concentrated in the range between 80.0 m and 320.0 m in diameter, 80% of those from the landing camera are concentrated in the range between 2.5 m and 15.0 m in diameter, and 80% of those from the panoramic camera are concentrated in the range between 0.2 m and 0.9 m in diameter. From this, the depth-to-diameter ratio decreases as the diameter increases, and craters of similar sizes have similar depth-to-diameter ratios, basically showing a normal distribution.
The relationship between the diameter of an impact crater and the depth-to-diameter ratio is illustrated in
Figure 10. We can see that for impact craters with different diameters, the values of depth-to-diameter ratios are dispersed, and as a whole, the smallest, largest and average values of the depth-to-diameter ratios show a decreasing trend as the diameters increases, with the decreasing smallest values being mainly distributed in the vicinity of 0.05, followed by the average value, and the largest values showing the most obvious decreasing trend. For the craters with diameters between 3 m and 10 m, the depth-to-diameter ratios do not show a decreasing trend and fluctuate in the plot. This could be due to chance, because of the small amount of data in this diameter range. The incorporation of additional data within this diameter range may facilitate a more precise delineation of the evolving trend in depth-to-diameter ratios as diameters fluctuate.
3.3. Degradation Classification Results
We classified the degradation of the impact craters by voting. The PCAM images do not contain any classified A and AB impact craters. Instead, they contain 279, 2905, and 7886 craters classified as Class B, BC, and C, respectively. These craters represent 2.54%, 26.20%, and 71.26% of all craters extracted from the PCAM images. Among the impact craters extracted from the LCAM images, there are 3, 4, 5, 8, and 76 impact craters classified as Class A, AB, B, BC, and C, respectively. These represent 1.03%, 4.12%, 5.15%, 8.25%, and 81.44% of the impact craters extracted from the LCAM images. Among the impact craters extracted from the satellite image, there are 1, 46, 105, 211, and 382 craters classified as Class A, AB, B, BC, and C, respectively. These represent 0.64%, 6.75%, 14.24%, 27.09%, and 51.28% of the impact craters extracted from the satellite image. When combined with the impact craters identified by multiple platforms and their classification results, it can be concluded that the majority of impact craters are classified as Class BC and C, indicating that the majority of impact craters are severely degraded. The majority of impact craters observed through the three imaging platforms were classified as Class C, representing over 50% of the total. As the diameter of the impact craters increases, the probability of fresh impact craters shows an increasing trend. However, this trend is not significant due to the small number of impact craters in the metre-sized and 100-m-sized range that have been recorded by LCAM images and satellite images. In general, among impact craters with diameters from 0.1 m to 500 m, the smaller the size of the crater, the higher the degree of their degradation.