1. Introduction
Landslides are the most common geological hazards in mountainous areas, which are characterized by large numbers, wide distribution, and severe destruction. The factors that induce landslides include rainfall, earthquake, snow melt, human activities, etc. Among them, landslides that are induced by strong earthquakes are higher in number, distribution range, and scale than any other types of landslides, and they have the characteristics of a concentrated and continuous distribution. Earthquake-triggered landslides (ETLs) often cause great damages to roads, houses, farmlands, oil and gas pipelines, and water management facilities in the earthquake area, resulting in large casualties and seriously affecting the post-earthquake rescue and post-disaster reconstruction. For example, the 2008 Wenchuan Ms 8.0 earthquake in Sichuan triggered more than 15,000 geological hazards [
1], which caused more than 20,000 deaths. Landslides accounted for up to 40% of all potential hazard sites induced by the earthquake. At some point, the extent and impact of earthquake-induced geological hazards exceed the hazards that are directly caused by the earthquake [
2].
The investigation of ETLs is one of the most urgent tasks for post-earthquake emergency relief and subsequent post-disaster reconstruction. Visual detection is the mostly used method to catalog ETLs, which is time-consuming and laborious, and often cannot meet the needs of earthquake emergency relief. Therefore, there is an urgent demand for automatic interpretation for investigating ETLs.
For example, Wen et al. used a Landsat series remote sensing data and digital elevation data (DEM) to conduct an in-depth analysis of the feature characteristics of landslide areas such as Mao County in Sichuan and Nayong in Guizhou [
3]. Given the fact that the damage to topography and landforms by landslide disasters shows a similar pattern in remote sensing data, all the landslides had been detected by analyzing such a pattern [
4,
5,
6]. Ramdhoni et al. performed a landslide extraction by the Smorph method using a slope and slope shape to establish the corresponding transformation matrix, and the landslide extraction accuracy reached 79.54% [
7]. In the field of landslide detection automation, machine learning techniques such as support vector machines, artificial neural networks, and deep learning have been widely used. Sameen et al. proposed a method for landslide detection that used residual networks to enhance the same designed network through fusing feature information to obtain better performance [
8]. Long et al. used the Mianyuan River basin, the extreme seismic area of the Wenchuan earthquake, as the study area to automatically detect geological hazards’ information in the basin. They used the maximum likelihood method and random forest algorithm based on high-resolution multitemporal satellite images such as RapidEye to quantitatively analyze the spatial and temporal evolution trends of geological hazards after the earthquake [
9]. Ye et al. used hyperspectral remote sensing data to detect landslides based on a deep learning framework with constraints (DLWC) [
10]. Their method improved the detection results through logistic regression (LR) by combining the detected image features with additional constraints as input. In recent years, in the field of machine learning on ETL detection, the neural network approach is more common [
11,
12,
13,
14,
15,
16]. The advanced development of deep learning techniques such as the convolutional neural network (CNN) has been widely successful in extracting information from images and it surpasses the other traditional methods [
17]. Liu et al. used a U-Net neural network and combined it with a GEE large remote sensing platform to achieve the rapid detection of large-area co-seismic landslides with a detection accuracy of about 70% [
18]. Most of the articles used medium-resolution images, and since ETLs often cause severe surface damage and land cover changes, all the above studies have shown the feasibility of remote sensing automatic detection technology in ETL investigation. However, an accurate comparison of multiple types of image effects with machine learning models is lacking.
With the introduction of different discriminative models, despite the fact that landslide automatic detection technology has been developed, there are still some problems, especially that most of the current studies are based on a single kind of remote sensing images for detection research. For ETLs, due to their wide distribution, a large number of remote sensing images of different types and resolutions are often used in emergency investigations [
19,
20,
21]. Additionally, such differences in landslide detection from multisource remote sensing images are rarely addressed.
On 5 September 2022, at 1252 h, an MS6.8 magnitude earthquake occurred in Luding County, Ganzi Prefecture, Sichuan Province, China, with a depth of 16 km. The epicenter was 29.59° N, 102.08° E. The seismic source mechanism of the earthquake is resolved as a walking-slip type rupture. As of 17:00 on 11 September 2022, a total of 93 people were killed by the earthquake. The earthquake-affected area is located in the Hengduan Mountains on the southeastern edge of the Tibet Plateau, which is a typical alpine valley area characterized by short and steep slopes, high altitudes, and broken rocks. These geomorphological features caused a large number of landslides after the earthquake. In this paper, the Wandong village in the Luding earthquake area was selected as the study area. Remote sensing images with different resolutions, such as UAV images, GF-6 satellite images, and Landsat 8 satellite images were used to automatically detect landslides using the full convolutional neural network (FCN). The differences in the automatic ETL detection accuracy from different images were quantitatively compared, and the reasons for these accuracy differences were analyzed. This study provides support in terms of the applicability and effectiveness of deep learning methods for further improvement in the field of ETL detection. It also holds a reference value for assessing the integrity of acquired landslide catalogs based on remote sensing images with different resolutions.
5. Comparison of Individual Detection of Landslide
The auto-detection only captured the area where landslides occurred but it did not distinguish different individual landslides. The area where ETLs developed would have a large number of landslides concentrated in patches. In order to study the effect of resolution on automatic landslide detection, this study used a manual method to distinguish different individual landslides from the automatically interpreted ones, especially by the criteria of geomorphological features to differentiate the landslides in mixed patches. The differences in automatically detected landslides based on remote sensing images of different resolutions were further analyzed for cataloguing the segmented landslides.
5.1. Number and Area of Individual Landslides
A total of 434 landslides were visually interpreted in the study area. Landslide areas were mainly concentrated in the range of 100–1000 m2 and 1000–10,000 m2, with 199 and 179 landslides, respectively, accounting for 87% of the total number of landslides. A total of 15 landslides were smaller than 100 m2 and 41 landslides were larger than 10,000 m2. Among them, the smallest landslide was 29 m2 and the largest landslide area was 84,191 m2. The smallest landslide that could be deciphered by UAV images was 46.91 m2, 73.54 m2 by GF-6 images, and 105.14 m2 by Landsat 8 images.
In general, the higher the resolution of the image, the greater the number of landslides detected. However, for the 199 landslides of 100–1000 m
2, UAV images detected 197 landslides, accounting for 98.99% of the landslides, while for GF-6 images and Landsat 8 images, the detection rate was improved, with 77 and 83 landslides, accounting for 39.07% and 41.71%, respectively. The details can be found in
Table 3.
When detecting landslides, although landslides can be detected by remote sensing images of different resolutions, the detected areas of landslides can also vary. For example, the actual area of a landslide located at 102°9′6″ E, 29°31′52″ N was 21,680 m2, while the area detected by UAV images was 20,762 m2, 14,248 m2 by GF-6 images, and only 10,975 m2 by Landsat 8 images. The area proportion follows the area trend. For the landslides of 1000–10,000 m2, the total area of landslides is 529,482.48 m2, the UAV image can detect 80.16% of the landslides, and the GF-6 image and Landsat 8 image are still lower, 37.63% and 32.54%, respectively.
According to the above data, we found that the larger the landslide area was, the better the image detection would be in all three types of different-resolution images. The larger the landslide area is, the more obvious the change in the surface coverage is, so it is easier to be identified and interpreted by FCN. In order to analyze the effectiveness of landslide detection, the average detection rate was used, which is the average of all individual landslide detection rates within the given range.
In our study, Landsat 8 images could not detect landslides below 100 m
2. Only landslides between 100–1000 m
2 could be identified and interpreted, but only in small amounts. For landslides of 1000–10,000 m
2, the average detection rate of UAV images was 35.58% higher than that of GF-6 images and 36.86% higher than Landsat 8 images. For landslides larger than 10,000 m
2, the average detection rate of UAV images was 16.93% higher than GF-6 images and 21.81% higher than Landsat 8 images. (
Figure 6)UAV images were more effective than the other two images in landslide detection regardless of area scales, and the smaller the area scale, the better the results would be. At the same time, the three types of images conformed to the understanding that the higher the image resolution, the better the detection. Only in detecting the area of 100–1000 m
2, the overall detection rate of GF-6 images was 17.11%, slightly lower than that of the Landsat 8 images, while the average detection rate is 9.08% lower due to some large-scale landslides in the Landsat 8 images that were not fully detected.
5.2. Frequency–Area Curves
Landslide area–frequency curves can represent the proportional characteristics of landslides at different sizes and also reflect the overall magnitude of landslide events triggered by earthquakes [
25,
26]. The catalogued landslide data detected from UAV, GF-6, and Landsat 8 images were plotted for producing area–frequency curves to compare with the reference landslide data (
Figure 7). The results showed that when a landslide area was larger than 2000 m
2, the area–frequency curves of the landslide detected by UAV, GF-6, and Landsat 8 images were more similar to those of the reference landslides with small differences.
When the landslide area was less than 2000 m2, the area–frequency curves of landslides detected by different images showed larger differences. The inflection point of the area–frequency curve of the reference landslide data appeared at the position of 250 m2, while for UAV images, the inflection point was located at the position of 100 m2. Additionally, GF-6 and Landsat 8 images showed two inflection points. The landslide frequency from UAV image interpretation was higher than the reference curve in the interval between 250 m2 and 2000 m2 regarding the landslide area. The landslide frequency from GF-6 and Landsat 8 images was lower than that of the reference, it was mostly because some landslides in the area could not be detected, which thus reduced the frequency. Additionally, since some landslides could not be detected as complete, the broken and incomplete landslides were therefore detected and counted in a smaller range.
As for the landslides with an area of less than 100 m
2, there was a higher frequency of landslides detected by GF-6 and Landsat 8 images. According to the conclusion above, it was caused by incomplete detection of landslides in larger areas.
Table 4 shows that GF-6 detected only one landslide within 100 m
2, while Landsat 8 images did not detect any landslides. Therefore, it suggests that the frequency of landslides in this part was composed of landslides that were incompletely detected in the range of 250 m
2 to 2000 m
2.
6. Error Analysis
6.1. Types of Automatic Remote Sensing Detection Errors
The types of errors that existed in landslide detection influenced by different resolutions were analyzed and summarized. Seven types of errors that existed in this study were found.(
Figure 8) Considering the fact that landslide detection areas are hardly identical, the error in the interval of area 95%–105% is relatively small, and was therefore, not included in this paper for error statistics. Landslide detection areas that were smaller than the actual area had more errors in the three images. Especially for GF-6 and Landsat 8 images, due to the relatively low resolution of the images and the large image element size, the mixed image elements in the landslide boundary area had a greater influence on landslide detection. This error type existed in 77 places in UAV image-based landslide detection, 83 places in GF-6 images, and 97 places in Landsat 8 images, accounting for 35.00%, 23.31%, and 24.49% of the error quantity, respectively (
Table 5).
The appearance of the error of the landslide detection area larger than the actual area on UAV images was mainly a case of misjudging the gully at the foot of the landslide as a landslide, thus making the landslide detection area larger than the actual area. This type of error existed in 17 places in UAV image-based landslide detection, reaching 7.73% of its error quantity. Additionally, this type of error was found in GF-6 and Landsat 8 images, mainly because of the misjudgment of similar feature attributes, especially in areas where the landslide area was more similar to the surrounding environmental feature attributes, such as roads, open spaces, etc. This type of error was 19 and 35 on GF-6 and Landsat 8 images, reaching 5.34% and 8.84% of their error numbers, respectively. The number of landslide errors identified was highest for GF-6 and Landsat 8 images, with 170 and 169 errors, reaching 47.75% and 42.68% of their errors, while this error was less frequent in UAV images.
6.2. Analysis of Spatial Distribution Errors
The relative position of landslide identification errors was analyzed by comparing the unidentified landslide areas, the misidentified areas, and the mass center of the landslide areas.
The distance between the mass center of the undetected landslide area and that of the corresponding landslide was calculated (
Figure 9). The origin position was the mass center of the landslide. The results showed that the area of the undetected landslides in UAV images was relatively small and mainly located in the area around the landslide far from the mass center of the landslide with an average distance of 17.61 m, and the farthest of 160.86 m. The area of the undetected areas in GF-6 and Landsat 8 images was relatively large and closer to the mass center of the landslide. When the undetected landslide area was located at the origin position, it means that the landslide was not detected as a whole. The average distance of the undetected area from the mass center of the corresponding landslide in GF-6 images was 11.76 m, and the farthest distance was 180.28 m. For Landsat 8 images, it was 13.65 m and 211.63 m, respectively.
The average distance between the misidentified area center and the corresponding landslide center (
Figure 10) in UAV images was 35.07 m, with the nearest distance being 2 m and the farthest distance being 486.04 m; 36.45 m, 1.41 m, and 186.10 m for GF-6 images; and 103.14 m, 21.63 m, and 194.65 m for Landsat 8 images. UAV image-based data of each type were the largest among the three types of images, and the distribution was the most dispersed, which suggested that the misidentified parts in UAV images were those farthest from the boundary part of the landslide area.
6.3. Landslide Detection and Land Cover Type
The accuracy of automatic ETL detection is influenced by land cover types of the landslide occurrence area. The five main land cover types included in this study are forestland, high-coverage grassland, low-coverage grassland, built land, and water area (
Table 6).
In this paper, the number of landslides on each land cover type was counted. Some of the landslide areas were under various land cover types, and the area counts were performed for the areas within multiple different land cover types in this paper.
In the four feature types, the correct interpretation rate decreased as the resolution decreased. The effect of resolution on the interpretation results was particularly evident on the built land, for example, industrial and mining, and urban and rural residential areas. The accuracy of interpretation in UAV images in this feature type reached 89.38%, which was the highest among the four feature types; while for GF-6 images, it was only 42.35%, the lowest among the four feature types; and for Landsat 8 images, the accuracy rate decreased to 7.37%, which was much lower than the accuracy of interpretation in the other three feature types.
As the resolution decreased, the percentage of landslide interpretation increased under the feature types of high-coverage grassland. This is because when landslides occur in this feature type, the surface damage is more obvious and the landslide area changes significantly from the surrounding non-landslide areas, so it is easier to be interpreted. Under the feature type of high-coverage grassland, the interpretation accuracy in UAV images was 17.33% higher than that in GF-6 images and 46.14% higher than that in Landsat 8 images. The interpretation accuracy in GF-6 images was 28.81% higher than that in Landsat 8 images.
As the resolution decreased, the percentage of landslides interpreted under the forest feature type decreased continuously. The percentage of landslides interpreted under the forest feature type in Landsat 8 images was 9.35% lower than that in UAV images, and 1.95% lower than that in GF-6 images. (
Figure 11) This is due to the fact that when landslides occur under this feature type, the surface damage is more obvious and the landslide area changes significantly from the surrounding area where no landslides occur, so it is easier to be detected. In contrast, for landslides under two feature types: low-coverage grassland and built land, Landsat 8 images detected 4.22% more landslides under the low-coverage grassland feature type than UAV images, and 3.34% more than GF-6 images; under the build land, Landsat 8 image-based detection results were higher than those of UAV images. Under the type of built land, Landsat 8 image-based detection result was 2.81% higher than that of UAV images and 1.58% higher than that of GF-6 images. The results suggested that as the resolution decreased, the proportion of landslide interpretation area with low vegetation coverage increased.
7. Conclusions and Discussions
We quantitatively compared and analyzed the ability and error types of remote sensing images with different resolutions: UAV images (resolution of 0.5 m), GF-6 images (resolution of 2 m), and Landsat 8 images (resolution of 15 m) in detecting the area and number of ETLs. The conclusions are as follows:
1. Differences in detecting ETL areas.
The resolution of the image plays a major role in automatic ETL detection. The higher the resolution of the image is, the higher the accuracy of landslide detection. In this study, the overall accuracy rates of UAV, GF-6, and Landsat 8 images were 94.45%, 88.78%, and 85.22%, respectively, and the precision rates were 86.01%, 70.54%, and 54.96%, respectively. The ETL area that was correctly detected in UAV images was 82.00%, against 52.67% and 48.32% when compared with GF-6 and Landsat 8 images, respectively.
2. Differences in detecting individual ETLs.
The ability of remote sensing images with different resolutions to detect individual ETLs varies greatly. The minimum landslide areas that were automatically detected in UAV, GF-6, and Landsat 8 images were 46.91 m2, 73.54 m2, and 105.14 m2, respectively, in this study. The quantity of ETLs detected in UAV images was 99.07%, against 60.83% and 61.06%, compared with GF-6 and Landsat 8 images, respectively.
The detection rate of remote sensing images for landslides of different scales is quite different as well. In general, the larger the landslide area is, the higher the detection accuracy. The study showed that when the area of landslide was smaller than 100 m2, 66.39% of the landslide was detected by UAV images, while for GF-6, it was only 0.70%. No landslide was detected by the Landsat 8 images. For landslide areas of 100–1000 m2, 74.29% landslides were detected in UAV images, while for GF-6 and Landsat 8 images the rate was only 17.11% and 21.06%, respectively. For 1000–10,000 m2, 80.16% of landslides were detected in UAV images, and only 37.63% and 32.54% for GF-6 and Landsat 8 images, respectively. For areas larger than 10,000 m2, 84.21% landslides were detected by UAV images, and only 64.78% and 61.40% for GF-6 and Landsat 8 images.
3. Errors of detecting ETLs.
The main errors in automatic ETLs detection by remote sensing include six types: the detected area of a single landslide was smaller than the actual area, larger than the actual area, partially overlaps with the actual area, a single landslide was detected as multiple landslides, multiple landslides were detected as one single landslide, and multiple landslides corresponded to multiple landslides. These types accounted for different proportions in remote sensing images of different resolutions. For example, the main error in UAV images was that the detected area of a single landslide was smaller than the actual area, while for GF-6 and Landsat 8 images, the errors were mainly about undetected landslides. In addition, the positions of the undetected and incorrectly detected landslide areas relative to the actual landslide are also different. The higher the image resolution, the closer the undetected landslide area is to the actual landslide, while the incorrectly detected landslide area is to the actual landslide.
The land cover has a great influence on landslide detection. ETLs are mainly distributed in low-coverage grassland, high-coverage grassland, and forest areas. In this study, low-coverage grassland, high-coverage grassland, and forest areas accounted for 38.10%, 30.37%, and 27.12% of the total landslide area, respectively, while landslides in built land occurred less. The GF-6 image had a higher detection accuracy of landslides in forestland, followed by low-coverage grassland and high-coverage grassland, and a lower detection accuracy of landslides in built areas; Landsat 8 images had a relatively high detection rate of landslides in low-grass covered areas, followed by grassland and forestland, but the detection rate was relatively low. The most prominent is that the detection of landslides on built land was particularly low, which also shows that Landsat 8 images cannot distinguish well between built land and landslides. Compared with GF-6 and Landsat 8 images, the image of land cover types detected in UAV images was less affected by the land cover type, including landslides on built land, which also had a high detection rate.
When carrying out ETLs surveys, the number and area of actual seismic landslides can be inferred using the results of this study based on the number of landslides interpreted from low- and medium-resolution remote sensing images for some areas where aerial photography is not available. In addition, this method is mainly based on the detection of land cover, so it is also applicable to rainfall-type landslides after rapid sliding, but for slowly deforming landslides, the detection method in this paper is not applicable because of their insignificant surface changes.
Obtaining a complete ETLs inventory by automatic detection is more difficult than visual detection. At present, many factors affect the accuracy of landslide detection [
27,
28]. For example, the FCN detection method takes pixels as the basic unit, thus producing a large number of debris polygons. It is still time-consuming and laborious to integrate and eliminate them manually. The detection accuracy of ETLs from medium-resolution images still needs to be improved. Further improvements to the algorithm model are needed. This study illustrated the ability and accuracy of ETL detection by remote sensing images of different resolutions on the one hand, and statistically analyzed the types of errors in ETL detection by remote sensing images with different resolutions. On the other hand, the model can be further optimized in order to provide the detection accuracy of ETLs by remote sensing. Additionally, the algorithm for debris polygons generated by landslide edges can be improved by adding some edge detection algorithms to make landslide identification more accurate.