收稿日期：2023-04-08；接受日期：2023-08-29

基金項目：國家自然科學(xué)基金項目（61966040）

作者簡介：韓俊，男，碩士研究生，研究方向為圖像處理和深度學(xué)習(xí)。E-mail：455804886@qq.com

通信作者：

王保云，男，副教授，博士，研究方向為機器學(xué)習(xí)及圖像處理。E-mail：wspbmly@163.com

Editorial Office of Yangtze River. This is an open access article under the CC BY-NC-ND 4.0 license.

文章編號：1001-4179（2024） 03-0123-11

引用本文：韓俊，王保云，徐繁樹.基于自校正原型網(wǎng)絡(luò)的泥石流災(zāi)害易發(fā)性評價：以怒江州為例［J］.人民長江，2024，55（3）：123-133.

摘要：

為解決泥石流易發(fā)性評價中因子選擇不一致造成的評價差異問題，以及目前神經(jīng)網(wǎng)絡(luò)不能有效提取泥石流特征以提升易發(fā)性評價正確率問題，提出了基于自校正原型網(wǎng)絡(luò)的泥石流災(zāi)害易發(fā)性評價方法。以溝谷為評價單元，提取溝谷的DEM、高分一號和Google Earth遙感影像作為訓(xùn)練數(shù)據(jù)，引入注意力機制和空洞空間卷積池化金字塔結(jié)構(gòu)構(gòu)建原型網(wǎng)絡(luò)的特征提取器，并使用自校正的方法優(yōu)化原型網(wǎng)絡(luò)的計算，將未發(fā)生泥石流的溝谷圖像輸入改進后的模型，計算其泥石流災(zāi)害易發(fā)性指數(shù)從而得出泥石流評價等級。運用該模型對怒江州的溝谷進行評價，并與歷史災(zāi)害數(shù)據(jù)進行對比。結(jié)果表明：模型分類正確率達到86.32%，評價結(jié)果中的易發(fā)區(qū)和高易發(fā)區(qū)均與歷史泥石流溝谷的空間分布較為吻合；相比于傳統(tǒng)評價方法，該方法能夠較好地自動學(xué)習(xí)遙感影像中泥石流特征，并實現(xiàn)災(zāi)害區(qū)域的快速識別與評價。研究成果可為泥石流災(zāi)害的研究提供新的思路。

關(guān)鍵詞：泥石流； 易發(fā)性評價； 原型網(wǎng)絡(luò)； 小樣本學(xué)習(xí)； 遙感影像； 怒江州

中圖法分類號： P954

文獻標(biāo)志碼： A" " " " " " " " " "DOI：10.16232j.cnki.1001-4179.2024.03.017

0引 言

泥石流是一種自然地質(zhì)災(zāi)害，發(fā)生時會對當(dāng)?shù)鼐用竦娜松砗拓敭a(chǎn)安全造成重大的危害。中國泥石流研究防治工作始于20世紀(jì)60年代，唐邦興等［1-2］通過對中國山地泥石流的考察研究，繪制了早期1∶6 000 000泥石流災(zāi)害分布與危險區(qū)劃圖。陳寧生［3-5］、崔鵬［6-8］等對四川、云南、新疆、西藏等特定地區(qū)的泥石流進行了考察，并提出了有效的泥石流判別方法和防治策略。隨著遙感技術(shù)的發(fā)展，從遙感產(chǎn)品中提取泥石流相關(guān)因子并使用動力學(xué)模型FLO-2D［9］、Debris2D［10］、DAN3D［11］等模擬特定溝谷的泥石流動力學(xué)過程，能夠較準(zhǔn)確地研究溝谷泥石流的形成機理并提出預(yù)警和防治方法。

為了實現(xiàn)大面積區(qū)域和更加精確的評價，李益敏［12］、孫濱［13］、趙巖［14］等將研究區(qū)域通過匯水面積或其他閾值劃分為小流域單元，統(tǒng)計各單元的泥石流相關(guān)因子，通過模型學(xué)習(xí)泥石流流域因子蘊含的特征來完成其他流域的危險性分類預(yù)測。然而，各學(xué)者劃分流域和選擇因子不同導(dǎo)致評價結(jié)果存在差異性，為了規(guī)避這些“主觀”選擇因子的過程，劉坤香等［15］改進殘差神經(jīng)網(wǎng)絡(luò)自動學(xué)習(xí)泥石流的DEM圖像特征進行溝谷的分類預(yù)測。楊小兵等［16］使用神經(jīng)網(wǎng)絡(luò)訓(xùn)練和識別泥石流堆積扇遙感圖像，以此尋找更多的泥石流隱患點。區(qū)域尺度上的災(zāi)害數(shù)據(jù)屬于小樣本問題，小樣本學(xué)習(xí)在遙感影像較少的情況下能夠取得更好的分類效果［17-18］，部分學(xué)者通過改進小樣本學(xué)習(xí)中經(jīng)典的原型網(wǎng)絡(luò)來進一步提高模型在處理遙感數(shù)據(jù)時的分類性能［19-21］。目前將小樣本方法用于泥石流災(zāi)害遙感的研究中還較為少見。

為了解決動力學(xué)模擬研究中無法快速高效進行大面積評價和泥石流因子評價中選擇因子不一導(dǎo)致的評價結(jié)果差異問題，本文提出了一種基于自校正原型網(wǎng)絡(luò)的泥石流災(zāi)害易發(fā)性評價方法。該方法以溝谷為評價單元劃分流域，利用溝谷的DEM、高分遙感和Google Earth遙感構(gòu)建溝谷泥石流災(zāi)害遙感數(shù)據(jù)集，并采用小樣本中經(jīng)典的原型網(wǎng)絡(luò)作為訓(xùn)練模型，自動學(xué)習(xí)溝谷遙感圖像中的泥石流特征并進行溝谷的泥石流災(zāi)害易發(fā)性評價。為了提升模型性能，引入注意力機制構(gòu)建一個高分辨率特征提取網(wǎng)絡(luò)作為原型網(wǎng)絡(luò)的特征提取器，同時，采用自校正的方法改進原型計算，并改進余弦退火學(xué)習(xí)率，使用漸進式焦點損失函數(shù)對模型進行優(yōu)化。最后，利用ArcGIS繪制溝谷易發(fā)性分區(qū)圖，并對結(jié)果進行了詳細分析。

1研究區(qū)域概況及泥石流遙感數(shù)據(jù)集構(gòu)建

云南省怒江傈僳族自治州（以下簡稱怒江州）是全國發(fā)生泥石流較嚴(yán)重的州（市）之一。怒江州地處云南省西北部橫斷山脈縱向嶺谷區(qū)，轄瀘水市、福貢縣、貢山縣和蘭坪縣，州內(nèi)高黎貢山、怒山山脈、云嶺山脈近南北向延伸，瀾滄江、怒江、獨龍江等深切峽谷平行南下。全州氣候具有年溫差小、日溫差大、干濕季分明、四季之分不明顯等低緯度高原季風(fēng)氣候特點。同時，因受地形地貌和緯度的影響，又具有北部冷、中部溫暖、南部熱、高山寒冷、半山溫暖、江邊炎熱等獨特立體氣候特征。高山峽谷的特殊地形、每年較高的降水量及其他因素導(dǎo)致怒江州泥石流災(zāi)害頻發(fā)。怒江州的地理位置和地形分布情況如圖1所示。

以溝谷流域為單元提取遙感影像構(gòu)建泥石流災(zāi)害遙感數(shù)據(jù)集，提取的溝谷圖像數(shù)據(jù)包含DEM、高分一號遙感和Google Earth遙感，詳細信息如表1所列。

數(shù)據(jù)的提取和處理如圖2所示。使用DEM提取河網(wǎng)，在河網(wǎng)中按照溝谷流域的溝口位置確定傾瀉點，使用分水嶺功能按照傾瀉點計算溝谷掩膜，對掩膜進行修正后按掩膜批量提取校正后的高分遙感影像、無偏移的Google Earth遙感影像和DEM圖像。將提取到的高分圖像和Google Earth圖像按照DEM圖像大小進行重采樣，并按DEM中最大溝谷圖像尺寸進行填0處理，將溝谷圖像統(tǒng)一填充為1 164×1 164大小。最后將同一溝谷的3種圖像進行疊加，共8個圖層通道。

數(shù)據(jù)集包含正樣本、負樣本和待評價樣本。正樣本從《云南減災(zāi)年鑒》［22］、相關(guān)文獻和網(wǎng)絡(luò)中獲得發(fā)生泥石流的山谷和村莊的精確地點信息，從而確定發(fā)生泥石流的溝谷。負樣本確定依據(jù)為在歷史記錄中不能檢索到溝谷發(fā)生泥石流信息、影像中沒有明顯的沖積痕跡和堆積扇且周圍村莊密集的溝谷。最終提取到正樣本82條、負樣本108條、其余溝谷樣本482條。

在正負樣本中，形態(tài)各異的溝谷潛在危險性存在特征差異，直接將溝谷分為正負樣本兩類訓(xùn)練，在忽略特征差異的同時網(wǎng)絡(luò)型也不穩(wěn)定。李益敏等［12］在泥石流小流域研究中，按Melton［23］指數(shù)對研究區(qū)域泥石流點的分布進行統(tǒng)計，0.141～0.441分段內(nèi)大約有60個泥石流點，0.441～0.741分段內(nèi)大約250個，0.741～1.341分段內(nèi)大約50個，其余兩個分段只有不到10個。趙巖［14］、張書豪［24］等的研究中泥石流點也有著同樣的分布情況。同時研究發(fā)現(xiàn)Melton因子對泥石流影響最為重要［12，14］，由此可見，Melton指數(shù)是泥石流溝谷的一個重要特征。使用K-means聚類的方法按Melton指數(shù)將正負樣本溝谷各預(yù)分3類，0，1，2類屬于正樣本，3，4，5類屬于負樣本，為保持?jǐn)?shù)據(jù)平衡，將個別類別樣本進行旋轉(zhuǎn)增強數(shù)據(jù)并按聚類標(biāo)簽設(shè)置監(jiān)督訓(xùn)練標(biāo)簽，結(jié)果如表2所列。

2基于自校正原型網(wǎng)絡(luò)的泥石流災(zāi)害易發(fā)性評價

2.1基于原型網(wǎng)絡(luò)的評價流程

基于原型網(wǎng)絡(luò)的溝谷泥石流災(zāi)害易發(fā)性評價流程如圖3所示，流程包括數(shù)據(jù)收集、溝谷提取、流域劃分、數(shù)據(jù)集劃分、數(shù)據(jù)集預(yù)分類、模型訓(xùn)練和預(yù)測、易發(fā)性評價和易發(fā)性分區(qū)圖等。

2.2原型網(wǎng)絡(luò)分類模型

如圖4所示，原型網(wǎng)絡(luò)通過特征提取器（feature extractor）提取支持集（Support Set）中每一類溝谷的特征計算該類的原型中心（Ck），將查詢集（Query Set）中待評價溝谷x^的特征（Qx^）與原型中心進行距離度量來完成待評價溝谷的分類和預(yù)測，并通過Softmax函數(shù)計算待評價溝谷所屬類別概率。

在訓(xùn)練和測試時，原型中心計算如式（1）所示：.

式中：Ck為每一類溝谷的原型中心，Sk為訓(xùn)練批次中每一類溝谷的樣本集合，Sk為該類的樣本數(shù)量；xi，yi分別為某一樣本和對應(yīng)的標(biāo)簽；fφ為特征提取器，φ為提取器的參數(shù)。

原型網(wǎng)絡(luò)使用余弦或歐氏距離度量方式作為分類器，通過計算溝谷樣本x^的特征與每類原型中心的距離大小來完成分類，并通過式（2）計算溝谷樣本x^屬于某一類溝谷的概率：

式中：pφ表示某一待評價溝谷樣本x^屬于真實類別k的概率，d為距離度量函數(shù)，Qx^為樣本x^通過fφ提取的特征，k′表示同一批訓(xùn)練的每一類別。在訓(xùn)練優(yōu)化階段，通過將式（3）中的目標(biāo)函數(shù)J（φ）使用反向傳播和隨機

梯度下降來對網(wǎng)絡(luò)模型參數(shù)φ進行優(yōu)化。.

2.3引入注意力機制的高分辨率特征提取器

在經(jīng)典的網(wǎng)絡(luò)資源、效率和性能無法滿足的情況下，結(jié)合GoogleNet中Inception多分支結(jié)構(gòu)和ResNet中殘差結(jié)構(gòu)的優(yōu)點，構(gòu)造了雙路卷積（two path convolution，TPC）結(jié)構(gòu)，并加入坐標(biāo)注意力機制（Coordinate Attention，CA）［25］和空洞空間卷積池化金字塔（Atrous Spatial Pyramid Pooling，ASPP）［26］結(jié)構(gòu)，提出基于注意力機制的高分辨率特征提取器（Attention-based High-Resolution Feature Extractor，AHRFE），結(jié)構(gòu)如圖5所示。

AHRFE首先通過卷積來對高分辨率遙感溝谷圖像進行下采樣，使用坐標(biāo)注意力機制在主干網(wǎng)絡(luò)中保留位置信息特征，接著使用空洞卷積池化金字塔結(jié)構(gòu)和雙路卷積結(jié)構(gòu)的多尺度特點來進行細節(jié)特征的提取，最后將特征進行降維和正則化，得到溝谷的深層特征。

2.4自校正原型計算改進

在樣本量極少的情況下，原型網(wǎng)絡(luò)比一般網(wǎng)絡(luò)有著更好的分類性能。但也存在著缺點：① 計算訓(xùn)練樣本特征的平均值作為原型，不能準(zhǔn)確地代表訓(xùn)練樣本情況；② 提取器的泛化能力較弱，對不同的數(shù)據(jù)魯棒性差。Cheng等［19］計算原型中心、訓(xùn)練樣本、測試樣本兩兩之間的損失來減少類別內(nèi)差異，擴大類別間差異，以提高分類性能。Liu等［20］將預(yù)測置信度較高的測試樣本加入到訓(xùn)練集中以減少類內(nèi)特征帶來的偏差從而提高網(wǎng)絡(luò)性能。Li等［27］通過優(yōu)化提取器的分類損失和度量損失來增加類間差異，最大化原型的代表能力。本文通過優(yōu)化原型與訓(xùn)練樣本及測試樣本之間的損失來對原型計算進行校正，改進損失計算和學(xué)習(xí)率策略來優(yōu)化網(wǎng)絡(luò)。其次，構(gòu)建基于注意力機制網(wǎng)絡(luò)結(jié)構(gòu)作為原型網(wǎng)絡(luò)的特征提取器，改進后提出自校正原型網(wǎng)絡(luò)（Self-Calibration Prototypes Net，SCPNet）。

2.4.1原型校正

如圖6所示，計算訓(xùn)練集（Support Set）的每類原型Ck，求得測試集（Query Set）與Ck之間的損失Lqc。同時計算訓(xùn)練集每類中樣本與該類Ck的損失Lsc，即自校正（Self-Calibration）過程。通過最小化Lqc和Lsc來達到校正Ck的目的，提高分類性能。

Lqc為測試集樣本與Ck之間的損失，計算如式（4）所示：.

式中：N為類別數(shù)，KQ為測試集樣本個數(shù)，pφi為測試樣本x^i屬于真實類別k的概率。

Lsc為訓(xùn)練集樣本與Ck之間的損失，即自校正損失，計算如式（5）所示：.

式中：KS為訓(xùn)練集樣本個數(shù)。

則總的損失為.

式中：λ為自校正的損失權(quán)重，λ∈R+。

將式（6）使用隨機梯度下降對模型參數(shù)φ進行優(yōu)化，使用優(yōu)化后的參數(shù)φ′按式（7）重新計算原型Ck′：.

按照式（1）～（7），在訓(xùn)練中循環(huán)計算損失及原型，最終不斷地優(yōu)化參數(shù)φ′，對原型的計算進行校正。

2.4.2漸進式焦點損失優(yōu)化

漸進式焦點損失函數(shù)（Progressive Focal Loss，PFL）［28］是在焦點損失函數(shù)［29］基礎(chǔ)上改進得到的，其對樣本不平衡問題和困難樣本問題有著更好的優(yōu)化能力。PFL的計算如式（8）～（13）所示。.

式中：αad表示不平衡樣本的調(diào)整權(quán)重，γad表示困難樣本的調(diào)整權(quán)重，n表示每個batchsize樣本數(shù)，pt表示樣本xi的onehot概率。

2.4.3改進余弦退火優(yōu)化

使用梯度下降算法優(yōu)化目標(biāo)函數(shù)時，按余弦方式對學(xué)習(xí)進行衰減被證明有著良好的優(yōu)化效果［30］。圖7中實線為學(xué)習(xí)率按余弦方式周期性變化，周期變化能夠防止損失函數(shù)陷入局部最小值。點線表示1/4個余弦周期重啟學(xué)習(xí)率衰減，在防止損失函數(shù)陷入局部最小值的同時減少收斂輪數(shù)。點劃線表示取半個余弦周期進行學(xué)習(xí)率衰減，有助于使網(wǎng)絡(luò)趨于平穩(wěn)，經(jīng)過實驗，如圖中改進衰減策略虛線所示，在1/4個周期上衰減到一定學(xué)習(xí)率后，穩(wěn)定學(xué)習(xí)率能夠帶來更好的學(xué)習(xí)效果。

2.5易發(fā)性指數(shù)

式中：Si表示樣本xi屬于第i類的預(yù)測概率也即相似概率，vi表示樣本特征向量，j表示分類數(shù)。

訓(xùn)練時，將正負樣本按Melton指數(shù)分布各分為3類進行訓(xùn)練學(xué)習(xí)，評價時，待評價樣本與正樣本中0，1，2類的其中一類的相似概率越高，說明待評價溝谷泥石流災(zāi)害易發(fā)性越高。使用式（15）將相似概率轉(zhuǎn)化為易發(fā)性指數(shù)，將正樣本每類概率之和與負樣本每類之和作差，差值大于等于0，表示樣本與正樣本相似程度大，取正樣本中的最大值作為易發(fā)性指數(shù)，值越大易發(fā)性越大。相反，差值小于0，取負樣本中最大值的相反數(shù)為易發(fā)性指數(shù)，值越小易發(fā)性越小。.

式中：Ii表示評價樣本xi的易發(fā)性指數(shù)，SP和SN分別表示xi樣本的正樣本概率之和和負樣本概率之和，m和u分別表示正樣本和負樣本的類別數(shù)，文中m=u=3。

3結(jié)果及分析

3.1分類性能

本次實驗實現(xiàn)的硬件為：CPU為Intel Xeon E5-2650 v3，內(nèi)存為64GB，GPU為NVIDIA GeForce RTX3090。軟件為Ubuntu 18.06、Pytorch 1.9.0、python 3.8、CUDA 11.1。實驗時，提取溝谷相關(guān)的泥石流因子，使用常見的機器學(xué)習(xí)算法學(xué)習(xí)因子特征并進行預(yù)測。同時選取常見深度學(xué)習(xí)網(wǎng)絡(luò)模型結(jié)構(gòu)，包括卷積結(jié)構(gòu)的網(wǎng)絡(luò)：Conv6、VGG、GoogleNet、ResNet，AHRFE加上分類器構(gòu)造的AHRNet分類模型進行學(xué)習(xí)和預(yù)測。然后，將這些網(wǎng)絡(luò)作為原型網(wǎng)絡(luò)的特征提取器進行實驗。進行5折交叉驗證后計算模型評價指標(biāo)，如表3所列。

可以看出，改進的特征提取器AHRFE構(gòu)建的網(wǎng)絡(luò)AHRNet比常見機器學(xué)習(xí)算法和一般的深度學(xué)習(xí)有著更好的分類預(yù)測性能。將AHRFE作為原型網(wǎng)絡(luò)的特征提取器，構(gòu)建的自校正原型網(wǎng)絡(luò)（SCPNet）總體正確率、查準(zhǔn)率、召回率、Kappa系數(shù)都優(yōu)于其他網(wǎng)絡(luò)。選取表3中3種評價方法中較好的分類模型，畫出正負二分類的ROC曲線如圖8所示。在ROC曲線中，SCPNet的AUC達到0.87，說明了SCPNet在識別泥石流發(fā)生情況上有著良好的分類性能。

3.2度量方式和損失計算的影響

文獻［19-20，27］使用余弦距離度量得到了更好的分類效果，文中對兩種度量進行對比，在本次實驗中，歐氏距離表現(xiàn)了更好的分類效果，如表4所列。

同時，文中使用了自校正原型優(yōu)化方式來校正原型的計算，按照式（6）對損失函數(shù)中的λ參數(shù)進行了多輪試驗發(fā)現(xiàn)，當(dāng)λ取6時，即總損失L=Lqc+6Lsc時，自校正效果最好，六分類和二分類正確率最高，如圖9所示。

3.3易發(fā)性評價

3.3.1易發(fā)性分區(qū)

將482條待評價溝谷及正負樣本輸入到訓(xùn)練后的自校正原型網(wǎng)絡(luò)中，計算出易發(fā)性指數(shù)并使用自然斷點法將指數(shù)分為5個自然段：不易發(fā)區(qū)［-0.961 6，-0.754 2］、低易發(fā)區(qū)（-0.754 2，-0.543 4］、較低易發(fā)區(qū)（-0.543 4，-0.231 2］、易發(fā)區(qū)（-0.231 2，0.599 4］、高易發(fā)區(qū)（0.599 4，0.999 7］。在ArcGIS中繪制易發(fā)性分區(qū)圖，如圖10所示。從圖10中可以看出，評價為易發(fā)及高易發(fā)的溝谷分布在瀾滄江以西，易發(fā)和高易發(fā)的溝谷分布在面積大、溝長較長的溝谷中。

表5顯示按照易發(fā)等級統(tǒng)計各評價數(shù)量情況。從表5可以看出，正樣本數(shù)目分布隨著易發(fā)性的增加而增加，負樣本的數(shù)據(jù)分布隨著易發(fā)性的增加而減少，這說明了自校正原型網(wǎng)絡(luò)預(yù)測的正確性。與其他學(xué)者［12，24］ 以5 km2為閾值提取的河網(wǎng)和河網(wǎng)節(jié)點劃分小流域的評測不同，很多大溝谷被劃分成多個小流域，這使得小流域數(shù)量較多，出現(xiàn)按照易發(fā)性增加而所占面積減少的情況，但是文中以溝谷單元劃分流域，更加注重溝谷型泥石流的易發(fā)性評價，多數(shù)歷史泥石流均分布在較大溝谷中，表中呈現(xiàn)隨易發(fā)性增加面積增加的情況，說明以溝谷流域評價，結(jié)果是合理的。

3.3.2易發(fā)性分析

為了進一步地討論評價結(jié)果的情況，從DEM、高分遙感和其他模型產(chǎn)品中提取得到與泥石流相關(guān)的16個因子，包含一級指標(biāo)因子8個，二級指標(biāo)因子8個，二級指標(biāo)使用一級指標(biāo)間接計算得到，這些指標(biāo)描述了溝谷的各項形態(tài)和特征，如表6所列。

分析時，根據(jù)Max等［36］篩選因子的方法，選出高程差、水系總長、平均坡度、Metlon、流域相對切割度、圓狀率、1 h最大降雨量、植被歸一化指數(shù)、標(biāo)準(zhǔn)化陸地水儲量指數(shù)共9個較不相關(guān)指數(shù)因子，按易發(fā)性分區(qū)統(tǒng)計它們的分布情況，如圖11所示，作為對比，正樣本分布如圖中虛線空心小提琴所示。從圖11可以看出，評價為易發(fā)和高易發(fā)的溝谷在分布上和正樣本非常相似，數(shù)據(jù)密度分布也非常接近，因此，改進后的原型網(wǎng)絡(luò)很好地學(xué)習(xí)了溝谷的泥石流特征，評價為易發(fā)和高易發(fā)的溝谷與歷史泥石流分布相似，評價結(jié)果合理。

3.3.3典型溝谷分析

在《云南減災(zāi)年鑒》中檢索到4條發(fā)生人員傷亡或多次發(fā)生泥石流的溝谷，地理位置如圖12（a）中A、B、C、D處所示，4條溝谷的溝口泥石流堆積和沖刷情況如圖12（b）～（e）所示。從圖12可以看出，東月谷和咪谷河暴發(fā)的泥石流規(guī)模較大，形成較大的堆積扇，而臘咱村和色仲村由于發(fā)生規(guī)模相對前兩條溝較小，發(fā)生后已經(jīng)著手治理，修建了排洪溝渠。

4條泥石流溝谷在模型中的相似概率和易發(fā)性指數(shù)如表7所列。

可以看出，模型對4條溝谷的預(yù)測全部正確，東月谷和咪谷河溝谷與正樣本2類有著較高的相似特征，并且它們的Melton指數(shù)也包含在Melton指數(shù)劃分的第2類區(qū)間，色仲村和臘咱村溝谷與正樣本0類有著較高的相似特征，它們的Melton指數(shù)包含在Melton指數(shù)劃分的第0類區(qū)間，這也說明了使用Melton對溝谷進行預(yù)劃分具有合理性。4條溝谷在李益敏等［12］評價中的結(jié)果為東月谷為低易發(fā)、咪谷河為高易發(fā)、臘咱村為高易發(fā)、色仲村為低易發(fā)。對于4條已經(jīng)發(fā)生泥石流的溝谷，文中的評價與歷史情況更為符合。

4結(jié) 論

使用改進后的原型網(wǎng)絡(luò)對溝谷泥石流遙感影像數(shù)據(jù)進行學(xué)習(xí)，并用于泥石流災(zāi)害易發(fā)性評價，通過實驗，得到以下結(jié)論：

（1） 文中構(gòu)建基于注意力機制的高分辨率網(wǎng)絡(luò)結(jié)構(gòu)作為原型網(wǎng)絡(luò)的特征提起器，采用自校正的方法計算原型中心，改進后的原型網(wǎng)絡(luò)比基于因子的機器學(xué)習(xí)分類模型和常規(guī)卷積網(wǎng)絡(luò)分類模型有著更好的分類預(yù)測性能，在正負樣本的分類精度達到86.32%，說明改進后的原型網(wǎng)絡(luò)模型有著較好的評價性能。

（2） 以溝谷流域為評測單元，將溝谷泥石流災(zāi)害易發(fā)性分區(qū)分為不易發(fā)區(qū)、低易發(fā)區(qū)、較低易發(fā)區(qū)、易發(fā)區(qū)和高易發(fā)區(qū)5個等級。高易發(fā)區(qū)溝谷主要集中在瀾滄江以西，分布在貢山縣、福貢縣、瀘水市中，蘭坪縣較少，與李益敏等［12］基于泥石流因子統(tǒng)計分析的方法得到的易發(fā)性分布有相似的結(jié)果，這與歷史泥石流溝谷在空間的分布上較為吻合。

（3） 提取溝谷的泥石流因子，使用相關(guān)性分析選出9個具有代表性的因子，分析其在5個易發(fā)性等級中的分布情況，結(jié)果顯示高易發(fā)區(qū)和較高易發(fā)區(qū)溝谷因子與正樣本有著相似的分布，說明改進原型網(wǎng)絡(luò)模型較好地學(xué)習(xí)了溝谷圖像的泥石流特征，且評價結(jié)果較為合理。

參考文獻：

［1］唐邦興，杜榕桓，康志成，等.我國泥石流研究［J］.地理學(xué)報，1980（3）：259-264.

［2］唐邦興，柳素清，劉世建.我國山地災(zāi)害的研究［J］.山地研究，1984（1）：1-7.

［3］陳寧生，楊成林，李戰(zhàn)魯，等.汶川地震次生泥石流形成發(fā)展趨勢與防治對策［C］∥汶川大地震工程震害調(diào)查分析與研究，成都，2009.

［4］陳寧生，謝萬銀，李戰(zhàn)魯.中國西南山區(qū)的泥石流分區(qū)與預(yù)測［J］.高原氣象，2004（增1）：134-140.

［5］陳寧生，胡桂勝.天山天池景區(qū)地質(zhì)災(zāi)害防治與生態(tài)環(huán)境保護關(guān)鍵技術(shù)研究及應(yīng)用［J］.科技促進發(fā)展，2015（3）：380-384.

［6］崔鵬，陳曉清，程尊蘭，等.西藏泥石流滑坡監(jiān)測與防治［J］.自然雜志，2010，32（1）：19-25，66.

［7］崔鵬，莊建琦，陳興長，等.汶川地震區(qū)震后泥石流活動特征與防治對策［J］.四川大學(xué)學(xué)報（工程科學(xué)版），2010，42（5）：10-19.

［8］崔鵬，楊坤，陳杰.前期降雨對泥石流形成的貢獻：以蔣家溝泥石流形成為例［J］.中國水土保持科學(xué)，2003（1）：11-15.

［9］林文，周偉，李婧，等.基于Flow-R和FLO-2D的溝谷型泥石流危險性評價［J］.人民長江，2022，53（5）：143-148.

［10］LIU K F，WU Y H.The assessement of debris flow hazard in korea using debris 2D［C］∥Interpraevent-international Symposium in Pacific Rim，2010.

［11］盛豪，胡卸文，熊沖沖.基于DAN3D的五花海泥石流運動過程模擬研究［J］.四川水力發(fā)電，2020，39（3）：11-16，93.

［12］李益敏，楊蕾，魏蘇杭.基于小流域單元的怒江州泥石流災(zāi)害易發(fā)性評價［J］.長江流域資源與環(huán)境，2019，28（10）：2419-2428.

［13］孫濱，祝傳兵，康曉波，等.基于信息量模型的云南東川泥石流災(zāi)害易發(fā)性評價［J］.中國地質(zhì)災(zāi)害與防治學(xué)報，2022，33（5）：119-127.

［14］趙巖.基于機器學(xué)習(xí)的白龍江流域潛在低頻泥石流溝識別［D］.蘭州：蘭州大學(xué)，2020.

［15］劉坤香，王保云，徐繁樹，等.基于殘差注意力機制的泥石流溝谷識別［J］.中國地質(zhì)災(zāi)害與防治學(xué)報，2022，33（6）：134-141.

［16］楊小兵，王登貴.基于多光譜影像和DEM的泥石流堆積扇識別研究：以白龍江流域武都段為例［J］.測繪與空間地理信息，2016，39（4）：118-121.

［17］茹穎.基于標(biāo)簽傳播與小樣本學(xué)習(xí)的高光譜圖像分類方法［D］.西安：西安電子科技大學(xué)，2021.

［18］張萌月.小樣本光學(xué)遙感影像目標(biāo)檢測識別技術(shù)研究［D］.北京：中國電子科技集團公司電子科學(xué)研究院，2021.

［19］CHENG G，CAI L，LANG C，et al.SPNet：siamese-prototype network for few-shot remote sensing image scene classification［J］.IEEE Transactions on Geoscience and Remote Sensing，2021，60：1-11.

［20］LIU J，SONG L，QIN Y.Prototype rectification for few-shot learning［C］∥Cham：European Conference on Computer Vision，2020.

［21］SNELL J，SWERSKY K，ZEMEL R S.Prototypical networks for few-shot learning［J］.arXiv Preprint arXiv：1703.05175，2017：1-11

［22］陳勤.云南減災(zāi)年鑒［M］.昆明：云南科技出版社，2018.

［23］MELTON M A.The geomorphic and paleoclimatic significance of alluvial deposits in southern arizona［J］.Journal of Geology，1965，73（1）：1-38.

［24］張書豪，吳光，張喬，等.基于子流域特征的泥石流災(zāi)害易發(fā)性評價［J］.水文地質(zhì)工程地質(zhì)，2018，45（2）：142-149.

［25］HOU Q，ZHOU D，F(xiàn)ENG J.Coordinate attention for efficient mobile network design［C］∥Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition，2021：13713-13722.

［26］CHEN L C，PAPANDREOU G，KOKKINOS I，et al.DeepLab：semantic image segmentation with deep convolutional nets，atrous convolution，and fully connected CRFs［J］.IEEE Transactions on Pattern Analysis and Machine Intelligence，2018，40（4）：834-848.

［27］LI H，CUI Z，ZHU Z，et al.RS-MetaNet：deep meta metric learning for few-shot remote sensing scene classification［J］.arXiv preprint arXiv：2009.13364，2020：1-12.

［28］WU B，YANG Y，YANG D，et al.Progressive hard-case mining across pyramid levels in object detection［J］.arXiv Preprint arXiv 2019.07217，2021：1-17.

［29］LIN T Y，GOYAL P，GIRSHICK R，et al.Focal loss for dense object detection［C］∥Proceedings of the IEEE International Conference on Computer Vision，2017：2980-2988.

［30］LOSHCHILOV I，HUTTER F.Sgdr：stochastic gradient descent with warm restarts［J］.arXiv Preprint arXiv：1608.03983，2016：1-16.

［31］方教勇.基于GIS的清平地區(qū)滑坡分布分形特征與危險性評價研究［D］.成都：成都理工大學(xué)，2018.

［32］劉欣.基于地形指標(biāo)的黃土高原坡面侵蝕溝發(fā)生與空間分異研究［D］.西安：西北大學(xué)，2021.

［33］王毅，方志策，牛瑞卿，等.基于深度學(xué)習(xí)的滑坡災(zāi)害易發(fā)性分析［J］.地球信息科學(xué)學(xué)報，2021，23（12）：2244-2260.

［34］連綱，郭旭東，傅伯杰，等.黃土高原縣域土壤養(yǎng)分空間變異特征及預(yù)測：以陜西省橫山縣為例［J］.土壤學(xué)報，2008，45（4）：577-584.

［35］BEVEN K J，KIRKBY M J.A physically based，variable contributing area model of basin hydrology［J］.Hydrological Sciences Journal，1979，24（1）：43-69.

［36］MAX K，KJELL J.Applied predictive modeling［M］.London：Springer，2013：84-87.

（編輯：劉 媛）