基于深度学习和高分遥感影像的青藏高原东部融冻泥流阶地提取研究
1.
2.
Extraction of solifluction terraces in the eastern Qinghai-Tibet Plateau based on deep learning and high-resolution remote sensing images
1.
2.
通讯作者:
收稿日期: 2022-05-01 修回日期: 2022-06-01
基金资助: |
|
Received: 2022-05-01 Revised: 2022-06-01
作者简介 About authors
王习敏,硕士生,主要从事冰缘地貌遥感智能解译研究.E-mail:
关键词:
Keywords:
本文引用格式
王习敏, 黄荣刚, 徐志达, 焦志平, 江利明.
WANG Ximin, HUANG Ronggang, XU Zhida, JIAO Zhiping, JIANG Liming.
0 引言
目前,针对青藏高原融冻泥流阶地观测较少,大范围制图研究基本空白。在先前的研究中,受到卫星遥感影像分辨率过低的限制,实地调绘和航空遥感是融冻泥流阶地制图的两种主要方式[7-11]。例如,1984—1990年,郭东信等[11]对青藏公路风火山垭口盆地融冻泥流阶地先后进行了两次地面调绘,现场测量还在很多泥流阶地站点开展过实验,如格陵兰东北部[12]、斯瓦尔巴群岛[13]、瑞典阿比斯库[14]、科罗拉多前沿山脉[15]等。自20世纪90年代开始,航空遥感与摄影测量技术逐渐用于泥流阶地提取和制图。Walsh等[10]使用ADAR-5500航空影像绘制了美国冰川国家公园的泥流阶地;Ridefelt等[7-8]利用航空摄影测量手段在瑞典阿比斯库山进行了多期融冻泥流阶地识别与分析。结合实地勘测调查和航空影像解译,陈游东[16]获取了中天山典型地区融冻泥流阶地的分布信息,并分析了其空间分布规律及对公路工程的灾害性影响。但是,以上手段只针对地块和小范围的泥流阶地,难以满足大范围制图的需求。
近年来,随着遥感技术的发展,高分辨率卫星遥感技术在融冻泥流阶地大范围识别中具有良好的应用潜力。然而,受融冻泥流阶地几何形态不一、表面覆盖多样(如裸地、草甸和稀疏植被)等因素影响,融冻泥流阶地纹理复杂,基于遥感影像的自动提取依然面临很大的挑战。快速发展的深度学习方法,如U-Net[17],SegNet[18]和DeepLab[19],在遥感影像地物分类和目标识别等方面取得了突破性的进展,在地球科学和遥感领域得到了广泛的应用[20-23]。其中,DeepLab V3+[24]深度学习模型,能获取多尺度上下文信息和清晰的边界,已用于青藏高原热融滑塌制图[25-26],有望为复杂环境下基于遥感影像的大范围泥流阶地自动提取提供一种新途径。
因此,本文提出了一种基于DeepLab V3+和高分辨率遥感影像的融冻泥流阶地提取方法,并选择甘孜州新都桥周边为研究区,在青藏高原开展融冻泥流阶地的遥感自动识别研究。首先,利用研究区Google Earth高分辨率遥感影像,通过人工勾勒方式绘制了泥流阶地几何边界,建立深度学习的样本数据集;然后,构建DeepLab V3+泥流阶地深度学习模型,并将训练好的模型用于高分辨率遥感影像进行泥流阶地自动提取;最后,对提取结果进行精度评定,并分析了泥流阶地空间分布特征及其与地形因子的关系。本研究为融冻泥流阶地的大范围遥感制图提供了新思路。
1 研究区与数据
研究区位于四川省甘孜藏族自治州,如图1,总面积为900 km2。该区域平均海拔3 945 m;年降水量约553 mm,降水量主要集中于6—7月,日降雨量最高达到50 mm;温度范围为-15~28 ℃,最高温度分布于7月,最低温度主要在1月和12月。该区域地表覆盖主要是草地、少量灌木和农田等。区域内融冻泥流阶地分布较多,主要集中在318国道、215省道沿线。
图1
本文选用Google Earth遥感影像作为实验数据,空间分辨率为0.59 m,数据采集时间为2016—2019年。此外,选用2008年的12.5 m分辨率ALOS PALSAR DEM数据,用于泥流阶地的地形因子分析。2020年30 m分辨率的地表覆盖数据用于分析融冻泥流阶地与地表下垫面的关系,数据来源于国家基础地理信息中心全球地表覆盖数据产品服务网站。
2 研究方法
2.1 DeepLab V3+模型
DeepLab V3+网络结构如图2所示,采用下采样编码和上采样解码结构。在编码阶段,通过骨干网络,如Xception等[29],结合多尺度空洞卷积,实现多尺度深度特征感知。在解码阶段,通过融合骨干网络低层次特征和空洞卷积高层次特征,并采用简单的上采样操作,逐步提高语义分割的分辨率,实现与输入图像一致大小的分割结果。其流程主要包括:采用骨干网络(如Xception65[29])进行深度可分离卷积运算,然后通过ASPP模块获取多尺度特征(空洞率分别为1, 6, 12, 18),并将ASPP空洞卷积层特征与全局池化层进行特征的融合,以此得到多尺度的高级语义特征输出结果;将从主干网络中提取的低级语义特征和编码阶段提取的高层次多尺度特征图进行融合,并采用3×3的卷积运算和4倍上采样,提高其空间分辨率。
图2
2.2 基于DeepLab V3+的泥流阶地自动提取
基于DeepLab V3+的泥流阶地自动提取处理流程如图3所示。首先,训练样本准备,制作泥流阶地数据集。然后,设置参数,对深度学习网络DeepLab V3+进行训练。其次,利用训练好的模型对研究区域Google Earth高分辨率影像进行预测、后处理。最后,对提取的泥流阶地进行精度评估。
图3
图3
基于DeepLab V3+的泥流阶地提取流程图
Fig. 3
Flowchart of solifluction terraces extraction based on DeepLab V3+
2.2.1 多类型泥流阶地正负样本制作
综合泥流阶地明暗相间的纹理特征、长条状几何形态、集群分布等特点,结合Google Earth三维地形分析,采用ArcGIS矢量化工具勾勒泥流阶地前沿陡坎边缘线作为样本。在区域A影像上勾勒泥流阶地陡坎边界,共计4 752条,主要包括高陡坎强纹理的泥流阶地、低陡坎弱纹理的泥流阶地、植被覆盖的泥流阶地和裸地覆盖的泥流阶地4种类型,如图4(a)~4(d)所示。此外,本文还勾勒出易混淆的其他目标作为负样本,如山脊线、山谷线、田地间陡坎、冲沟等线状目标区域,如图4(e)~4(h)所示,共计110个多边形。最终勾勒的泥流阶地和非泥流阶地样本分布如图1所示。基于人工勾勒的正负样本矢量多边形,结合对应的Google Earth高分辨率影像,生成大小为513×513的正负样本切片数据集,所得正负样本切片数量分别为3 046和1 863。
图4
图4
泥流阶地人工勾勒正负样本示例
Fig. 4
Examples of manually delineated positive and negative samples: solifluction terraces with high risers and strong textures (a); solifluction terraces with low risers and weak textures (b); solifluction terraces covered by bare soil (c); solifluction terraces covered by vegetation (d); ridge lines (e); valley lines (f); raised paths through farmlands (g); gullies (h)
2.2.2 DeepLab V3+模型训练与泥流阶地提取
基于生成的泥流阶地数据集,对DeepLab V3+模型进行迭代训练,获取其模型参数。在训练过程中,参数设置如下:批处理量为12;迭代次数为6×104次;初始学习率为1×10-3;权重衰减系数设置为4×10-5;采用的优化器是Adam;骨干网络为Xception65;ASPP结构中空洞率分别为[1,6,12,18],ASPP模块可有效提取不同尺度下的语义特征,一定程度上有助于提高不同尺度泥流阶地的提取效果。
在模型训练好后,将其用于整个实验区的融冻泥流阶地预测提取。首先,将实验区所有Google Earth高分辨率影像切割成513×513大小的切片;然后,利用训练好的模型对各切片进行泥流阶地提取,将切片分割成泥流阶地陡坎和其他背景像素。最后,使用GDAL进行镶嵌处理,输出完整的预测结果。
为了剔除错误提取的泥流阶地,结合几何特征、集群存在特征、与等高线相对平行的走势关系等先验知识,对预测结果进行后处理。后处理的规则主要包括:(1)剔除面积小于40 m2、长度小于10 m、宽度小于4 m的提取目标;(2)删除条数少于10条的零散分布提取目标;(3)去除个别与等高线走势相悖的目标。
2.3 精度评估
为评估泥流阶地提取方法的性能,采用人工勾绘的结果作为参考,分别根据公式(
式中:
3 结果与分析
3.1 泥流阶地提取结果
图5为基于深度学习的融冻泥流阶地深度学习提取结果,共9 203条,总面积为2.12 km2。据提取结果发现,在本研究区内,泥流阶地沿山谷线两侧分布,主要位于国道等不同等级道路两侧,呈现了聚集分布。
图5
图5
基于DeepLab V3+模型的泥流阶地自动提取结果
Fig. 5
Results of solifluction terraces based on DeepLab V3+
对于不同覆盖类型具有明显特征的泥流阶地,本文方法在目标提取率和边界吻合程度等方面均取得很好的效果,泥流阶地的边界提取效果较好,但类型间也存在差别。如图6所示,植被覆盖的泥流阶地,条状形态、纹理不大明显的部分,提取结果有待提高;陡坎较高的泥流阶地,边界位置相对精准;陡坎纹理不清的泥流阶地,深度学习模型仅能提取到部分结果。此外,相比较于形态尺度较小的泥流阶地陡坎,形态尺度较大的泥流阶地陡坎提取效果更好。
图6
图6
泥流阶地提取结果与人工判读比对
Fig. 6
Comparison between automated extraction solifluction terraces and the interpretation results: solifluction terraces covered by vegetation (a); solifluction terraces with high risers (b); solifluction terraces covered by meadow (c); solifluction terraces covered by bare soil (d)
表1 融冻泥流阶地提取结果精度评估
Table 1
评估指标 | TP/km2 | FP/km2 | FN/km2 | Precision | Recall | F1 |
---|---|---|---|---|---|---|
区域A精度 | 1.350 | 0.557 | 0.401 | 0.708 | 0.771 | 0.738 |
区域B精度 | 0.774 | 0.280 | 0.446 | 0.733 | 0.634 | 0.680 |
综上,泥流阶地提取方法性能分析表明,DeepLab V3+深度学习模型能够有效提取研究区泥流阶地陡坎的边界,且拥有良好的泛化能力,在形态、纹理相似的泥流阶地区域(如G318国道沿线)具有推广应用潜力。
3.2 实地考察验证
图7
图7
提取结果的实地验证
Fig. 7
Field validation of extraction results: on-site photos (a); results comparison (b)
4 讨论
4.1 地形因子统计分析
为分析泥流阶地的空间分布特征,本文讨论了研究区内泥流阶地与地形因子的统计关系。图8为区域A、区域B的地形因子(坡度、坡向、高程)统计图。结果显示:
图8
图8
泥流阶地的地形因子统计分析
Fig. 8
Statistical analysis of topographic factors of solifluction terraces: aspect (a), elevation (b), slope (c) of Area A and B
(2)在研究区B内,泥流阶地坡向主要集中在正西、西北和东北等阴坡方向,与研究区A内的泥流阶地具有类似的坡向分布特征[图8(a2)]。在高程分布方面,泥流阶地主要集中在3 650~3 750 m分布范围内,占研究区B内泥流阶地的43%;区域内泥流阶地平均高程为3 674.694 m,标准差为86.240 m。在坡度上,集中分布在10°~30°,并且在20°~25°区间的泥流阶地数量最多。
总体而言,本研究区域内泥流阶地以西北、北等阴坡方向为主;在坡度上,集中分布在10°~30°之间;在高程上,3 650~3 750 m区间为主要分布范围。
4.2 几何因子统计分析
此外,本文还对泥流阶地陡坎的几何因子(面积、周长、长度、宽度)进行了统计分析和讨论,如图9所示。统计结果显示:
图9
图9
泥流阶地几何因子统计
Fig. 9
Statistics of geometric factors of solifluction terraces: area (a); perimeter (b); length (c); width (d)
(1)在研究区A中,泥流阶地面积平均值为369 m2,标准差为330.79 m2。绝大部分泥流阶地面积小于1 000 m2,占总数的95%。其中,小于500 m2的泥流阶地占79%。泥流阶地周长平均为128.46 m,标准差为79.59 m,最大值和最小值分别为890.65 m和20.50 m。长度、宽度大部分小于500 m,均值分别为54 m、26 m。
(2)在研究区B中,泥流阶地面积均值463 m2,标准差为412.46 m2,小于1 000 m2的占91%;周长平均为146.84 m,标准差为86.36 m,100~200 m区间最为集中;长度、宽度均值分别为62 m、29 m。
总体而言,研究区A和研究区B内的泥流阶地在几何形态上具有较高的相似性,但也存在一定差异。例如,研究区B的平均长度和宽度均大于研究区A。
4.3 地表覆盖
图10
图10
泥流阶地地表覆盖空间分布图
Fig. 10
Spatial distribution map of surface cover on solifluction terraces
图11
图11
泥流阶地地表覆盖类型统计
Fig. 11
Statistics on surface coverage types of solifluction terraces
5 结论
基于DeepLab V3+深度学习网络和Google Earth高分辨率卫星影像,本文提出了一种大范围融冻泥流阶地遥感自动提取的方法,并在四川省甘孜州新都桥地区开展了实验研究,分析了该区泥流阶地分布特征及其与地形因子、几何因子、地表覆盖等因素的关系。主要结论如下:
(1)与人工解译结果对比,本方法提取结果的综合精度达到0.68以上,并经野外调查验证了其有效性。陡坎较高的泥流阶地提取精度最高,对于弱纹理的泥流阶地精度较低但能实现其定位。
(2)共识别了9 203条融冻泥流阶地,主要分布在新都桥镇附近的山谷两侧,海拔高程主要分布在3 650~3 750 m区间,且主要位于20°~25°的阴坡上。
(3)泥流阶地面积和周长分别在1 000 m2、150 m以内,地表覆盖类型主要是草地,占比高达95%。
参考文献
Solifluction rates, processes and landforms: a global review
[J].
Evolution of permafrost and environmental changes of cold regions in eastern and interior Qinghai-Tibetan Plateau since the Holocene
[J].
青藏高原中、东部全新世以来多年冻土演化及寒区环境变化
[J].
Evolution of permafrost in China during the last 20 ka
[J].
两万年来的中国多年冻土形成演化
[J].
Principal characteristics of permafrost in China
[J].
我国多年冻土的主要特征
[J].
The impact of latitude and altitude on the extent of permafrost during the Last Permafrost Maximum (LPM) in North China
[J].
Arctic soil patterns analogous to fluid instabilities
[J].
Observations on regional variation in solifluction landform morphology and environment in the Abisko region, northern Sweden
[J].
Measurement of solifluction rates using multi-temporal aerial photography
[J].
Spatial analysis of solifluction landforms and process rates in the Abisko Mountains, northern Sweden
[J].
Solifluction steps and risers, Lee Ridge, Glacier National Park, Montana, USA: a scale and pattern analysis
[J].
A preliminary of soliflution terraces in Fenghuoshan Pass Basin on Qinghai-Xizang Plateau
[J].
青藏公路风火山垭口盆地融冻泥流阶地初步研究
[J].
Instrumental observations of mass wasting in the Mesters Vig District, Northeast Greenland
[J].
Relations between slow slope processes and active-layer thickness 1972—2002, Kapp Linné, Svalbard
[J].
56 years of solifluction measurements in the abisko mountains, northern Sweden: analysis of temporal and spatial variations of slow soil surface movement
[J].
Downslope soil movement in a Colorado alpine region: rates, processes, and climatic significance
[J].
Study on periglaciation and the countermeasures of highway disaster in middle Mount Tianshan
[D].
中天山地区冰缘作用与公路灾害防治对策研究
[D].
U-net: convolutional networks for biomedical image segmentation
[C]//
Segnet: a deep convolutional encoder-decoder architecture for image segmentation
[J].
Semantic image segmentation with deep convolutional nets and fully connected CRFs
[J].
Machine learning for the geosciences: challenges and opportunities
[J].
Deep learning for remote sensing data: A technical tutorial on the state of the art
[J].
Deep learning in remote sensing: a comprehensive review and list of resources
[J].
Automatic mapping of thermokarst landforms from remote sensing images using deep learning: A case study in the Northeastern Tibetan Plateau
[J].
Encoder-decoder with atrous separable convolution for semantic image segmentation
[C]//
Using deep learning to map retrogressive thaw slumps in the Beiluhe region (Tibetan Plateau) from CubeSat images
[J].
The characteristics and patterns of retrogressive thaw slumps developed in permafrost region of the Qinghai-Tibet Plateau
[J].
青藏高原多年冻土区热融滑塌发育特征及规律
[J].
Image recognition method for transmission line based on the DeepLab V3+ deep convolutional network
[J].
基于DeepLab V3+深度卷积网络的输电导线图像识别方法
[J].
Road extraction of high resolution remote sensing imagery based on DeepLab V3
[J].
利用Deeplab V3提取高分辨率遥感影像道路
[J].
Xception: deep learning with depthwise separable convolutions
[C]//
A bidirectional analysis method for extracting glacier crevasses from airborne LiDAR point clouds
[J].
/
〈 |
|
〉 |
