基于異軌DSM提取技術(shù)的滑坡體積估計(jì)

庾 露,黃 艷 霞

(1.南寧師范大學(xué) 地理科學(xué)與規(guī)劃學(xué)院,廣西 南寧 530100; 2.廣西國(guó)土資源規(guī)劃設(shè)計(jì)集團(tuán)有限公司,廣西 南寧 530029)

0 引 言

滑坡是常見(jiàn)的自然災(zāi)害之一,極易誘發(fā)次生災(zāi)害,對(duì)人類的生存環(huán)境造成極大的威脅,一直是災(zāi)害防治工作中的重中之重[1]。據(jù)自然資源部統(tǒng)計(jì),2020年全國(guó)共發(fā)生地質(zhì)災(zāi)害7 840起,其中滑坡4 810起,占總數(shù)的61%,滑坡災(zāi)害呈常態(tài)化趨勢(shì)[2]。獲取滑坡的體積是研究滑坡災(zāi)害的一個(gè)重要內(nèi)容,滑坡體積對(duì)滑坡的破壞力、展布范圍及其可能引發(fā)的二次災(zāi)害如堰塞湖和洪水具有指示意義[3-5],精確的滑坡測(cè)量能降低滑坡相關(guān)的危險(xiǎn)度和風(fēng)險(xiǎn)性分析中的不確定因素并且有助于進(jìn)一步的防災(zāi)減災(zāi)[6]。

數(shù)字表面模型(Digital Surface Model,DSM)差值計(jì)算是測(cè)量滑坡體積方法之一,通常使用遙感測(cè)繪技術(shù)獲取滑坡前后的高程分布,利用高程差和像元面積的乘積獲取滑坡體積,具有覆蓋范圍較大、時(shí)間和人力成本較低的優(yōu)點(diǎn)。新一代遙感衛(wèi)星和無(wú)人機(jī)平臺(tái)的成熟并大規(guī)模投入使用,以及LiDAR、InSAR等測(cè)量技術(shù)的發(fā)展,為地質(zhì)災(zāi)害監(jiān)測(cè)提供了多樣的數(shù)據(jù)來(lái)源和分析手段[7]。如MA[8],彭大雷[9],CHANG[10]等使用低空無(wú)人機(jī)測(cè)量計(jì)算出滑坡體積,此方法基本不受云層遮擋,但需要攜帶設(shè)備到達(dá)現(xiàn)場(chǎng)開(kāi)展工作。同時(shí),受飛行平臺(tái)穩(wěn)定性和載荷量限制,無(wú)人機(jī)遙感系統(tǒng)所獲取的圖像往往像幅較小、數(shù)量較多、畸變較大,圖像處理難度較高[11]。衛(wèi)星遙感具有覆蓋范圍大、成本低、精度高的特點(diǎn),也可作為滑坡體積估算的手段[12],但由于滑坡發(fā)生后,局部通常伴隨長(zhǎng)時(shí)間降雨,會(huì)對(duì)衛(wèi)星觀測(cè)造成諸多不利影響。傳統(tǒng)的光學(xué)立體測(cè)繪衛(wèi)星受固有重訪周期限制、拍攝調(diào)度安排和降雨天氣云層的遮擋,往往在單顆衛(wèi)星過(guò)境的短暫間隙內(nèi)無(wú)法獲取有效影像。星載重復(fù)軌道InSAR測(cè)量在成像過(guò)程中,大氣水汽含量的時(shí)空變化會(huì)引起雷達(dá)信號(hào)傳播延遲,給干涉相位帶來(lái)較大的附加相位變化,導(dǎo)致干涉測(cè)量結(jié)果出現(xiàn)誤差[13]。因此,有必要根據(jù)研究區(qū)所收集的影像數(shù)據(jù)和地形環(huán)境情況,進(jìn)一步拓展滑坡地形測(cè)量的手段。

基于光學(xué)衛(wèi)星的異軌立體測(cè)量技術(shù),是對(duì)傳統(tǒng)光學(xué)測(cè)繪衛(wèi)星同軌測(cè)量方法的有效補(bǔ)充。該技術(shù)可將未設(shè)計(jì)用于立體測(cè)繪的衛(wèi)星數(shù)據(jù)充分利用,借助單顆或多顆衛(wèi)星重復(fù)過(guò)境縮短重訪周期的優(yōu)勢(shì),多期覆蓋研究區(qū),將可能存在的有效影像篩選出來(lái)并組成異軌立體像對(duì),生成DSM產(chǎn)品。已有學(xué)者使用該技術(shù)提取了地物高程,如胡芬等[14]對(duì)比分析了異軌和同軌立體影像生成的DSM,發(fā)現(xiàn)二者精度趨近。倪文儉等[15]使用高分2號(hào)衛(wèi)星(GF-2)數(shù)據(jù)形成異軌立體像對(duì),提取森林高度,與LiDAR結(jié)果的均方根誤差為3.6 m。

上述研究成果驗(yàn)證了該技術(shù)的可行性,但目前的研究仍主要集中于單一成像時(shí)間點(diǎn)上的DSM提取,尚缺乏考慮滑坡事件發(fā)生前后高程劇烈變化產(chǎn)生的數(shù)據(jù)一致性問(wèn)題,以及在地形起伏大、無(wú)實(shí)地控制點(diǎn)約束下的精度控制情況。因此,有必要就如何將異軌立體測(cè)量技術(shù)應(yīng)用于滑坡體積提取場(chǎng)景做更深入的研究。本文以2020年7月11日湖北省恩施市馬者村沙子壩滑坡為例,重點(diǎn)分析使用高分2號(hào)衛(wèi)星(GF-2)通過(guò)異軌立體成像立體測(cè)量技術(shù)在無(wú)實(shí)地控制點(diǎn)的場(chǎng)景下,提取滑坡后DSM的可行性和精度。同時(shí),結(jié)合資源3號(hào)02星(ZY-302)影像提取的滑坡前DSM,估算滑坡體積,并對(duì)滑坡的構(gòu)成進(jìn)行分析。

1 研究區(qū)概況

受梅雨期連續(xù)降雨影響,2020年7月20～21 日湖北省恩施市清江上游馬者村沙子壩(東經(jīng)109°18′,北緯30°22′)滑坡垮塌,發(fā)生縱向長(zhǎng)1 200～1 500 m,橫向?qū)?20～580 m的特大型滑坡,導(dǎo)致清江干流河道形成堰塞湖,水位的急速上升導(dǎo)致存在潰壩風(fēng)險(xiǎn)[16]。研究區(qū)地處巫山余脈和武陵山北上的余支交會(huì)部,地勢(shì)整體為西低東高,最大落差達(dá)1 380 m,地形崎嶇復(fù)雜,多為海拔 1 000 m 以上的高山。在氣候上,研究區(qū)屬亞熱帶大陸性季風(fēng)氣候,雨量充沛,易發(fā)生暴雨洪澇,造成滑坡。具體的影像范圍與研究區(qū)范圍如圖1所示。

圖1 研究區(qū)域位置和影像范圍Fig.1 Location and image coverage of the study area

2 數(shù)據(jù)獲取與研究方法

2.1 數(shù)據(jù)來(lái)源

滑坡前的DSM數(shù)據(jù)由2020年3月19日的資源三號(hào)02(ZY-302)衛(wèi)星影像生成,這是距離發(fā)生滑坡時(shí)間較近、成像質(zhì)量較高的衛(wèi)星影像數(shù)據(jù)。ZY-3和后續(xù)發(fā)射的ZY-302星均是中國(guó)民用高分辨率立體測(cè)繪衛(wèi)星,通過(guò)星載三線陣光學(xué)相機(jī),以同軌方式實(shí)現(xiàn)地形圖繪制?；潞蟾采w研究區(qū)域的可用影像較少,以ZY-3和ZY-302為例,在滑坡發(fā)生后180 d內(nèi)僅有3景影像覆蓋,且滑坡區(qū)受到云層遮擋。而在此期間,高分二號(hào)(GF-2)衛(wèi)星以不同角度拍攝了多景覆蓋滑坡區(qū)的影像,雖然該衛(wèi)星并不是為立體測(cè)繪任務(wù)而設(shè)計(jì)的,但具有最高0.8 m的高空間分辨率和最大35°的側(cè)擺能力,在不同觀測(cè)角度和較短重復(fù)觀測(cè)周期的前提下可構(gòu)成異軌立體觀測(cè)像對(duì),從而實(shí)現(xiàn)DSM提取[15]。經(jīng)過(guò)反復(fù)試驗(yàn),本文選取了距離滑坡發(fā)生后時(shí)間較接近的2020年11月08日和2020年12月22日的兩幅GF-2影像構(gòu)成立體像對(duì),用以提取滑坡后DSM。

本文使用ESRI公司提供的全球衛(wèi)星影像網(wǎng)絡(luò)服務(wù)作為平面位置參考,使用GLO30 DSM數(shù)據(jù)作為高程參考。前者在研究區(qū)范圍內(nèi)提供了分辨率為1 m,經(jīng)過(guò)地形校正和拼接處理的衛(wèi)星影像底圖,其分辨率與GF-2全色波段影像的分辨率接近,并優(yōu)于ZY-302影像的分辨率,可用于生成控制點(diǎn)和驗(yàn)證點(diǎn)的平面坐標(biāo);后者為歐空局哥白尼地球觀測(cè)計(jì)劃的數(shù)字表面高程產(chǎn)品,數(shù)據(jù)來(lái)源于TanDEM-X任務(wù)期間獲得的SAR衛(wèi)星資料,GLO30的像元尺寸為30 m,絕對(duì)垂直精度<4 m,絕對(duì)水平精度<6 m[17],研究區(qū)最新數(shù)據(jù)更新至2014年可用于生成控制點(diǎn)和驗(yàn)證點(diǎn)的高程值。

2.2 研究方法

2.2.1技術(shù)流程

對(duì)于滑坡區(qū)DSM提取問(wèn)題,受其所處地理位置、地形地貌和自身不穩(wěn)定地質(zhì)環(huán)境等多方面因素制約,可能不利于開(kāi)展實(shí)測(cè)工作。為此,本文提出完全基于衛(wèi)星遙感數(shù)據(jù)源進(jìn)行異軌DSM提取和精度校正的新思路,并在后文的論述中討論該方法的可行性。本次研究的技術(shù)流程如圖2所示,使用1景ZY-3 02衛(wèi)星影像提取滑坡前DSM,使用2景異軌的GF-2衛(wèi)星影像提取滑坡后DSM,使用ESRI衛(wèi)星全球衛(wèi)星影像作為平面位置參考,GLO30作為高程參考對(duì)提取的DSM進(jìn)行位置和高程校正和精度評(píng)價(jià),最后計(jì)算滑坡體積。

圖2 技術(shù)路線Fig.2 Technical route

2.2.2異軌立體測(cè)量基本原理

異軌立體測(cè)量是基于光學(xué)衛(wèi)星所具備的側(cè)擺功能,由單臺(tái)相機(jī)在不同軌道觀測(cè)獲取立體影像的方法。如圖3所示,在兩個(gè)不同軌道上分別對(duì)目標(biāo)成像,產(chǎn)生類似于雙目立體視覺(jué),通過(guò)獲取兩個(gè)軌道的基線、焦距以及視差等信息可計(jì)算出地物的三維點(diǎn)云分布。異軌立體測(cè)量的精度包括平面精度和高程精度。其中平面精度和同軌立體測(cè)量方式接近[18-20],可通過(guò)衛(wèi)星數(shù)據(jù)供應(yīng)商提供的RPC(Rational Polynomial Cofficient)參數(shù),運(yùn)用成像幾何模型,并配合一定數(shù)量的地面控制點(diǎn)建立圖像坐標(biāo)和地理坐標(biāo)的幾何關(guān)系。而高程精度則主要受飛行平臺(tái)的軌道精度、影像空間分辨率、成像紋理清晰度、基高比等影響。對(duì)于異軌立體的測(cè)量,通常前3個(gè)因素基本相同或接近,而基高比可根據(jù)所收集的影像進(jìn)行較靈活的調(diào)整。理論上基高比越大,交會(huì)條件越好,高程測(cè)量的精度越高。根據(jù)公式(1)基高比B/H由兩軌影像的側(cè)視角α和β決定:

圖3 基高比示意Fig.3 Baseline-to-height ratio schematic

(1)

式中:B為軌道間基線長(zhǎng)度,H為攝站高度,α和β分別為側(cè)視角。經(jīng)計(jì)算,GF-2異軌像對(duì)的基高比約為0.50,ZY-302衛(wèi)星立體像對(duì)的基高比約為0.83。

2.2.3利用GF-2提取DSM

由于本次滑坡體積提取實(shí)驗(yàn)未到達(dá)現(xiàn)場(chǎng)采集地面控制點(diǎn),因此用于高精度幾何定位的控制點(diǎn)和用于DSM結(jié)果精度評(píng)估的驗(yàn)證點(diǎn)均通過(guò)已公開(kāi)的免費(fèi)遙感數(shù)據(jù)產(chǎn)品采集獲得。本文使用ArcGIS Desktop載入全球衛(wèi)星影像網(wǎng)絡(luò)服務(wù)和GLO30 DSM數(shù)據(jù),采集地形起伏較小、地表紋理清晰穩(wěn)定、易于辨認(rèn)的人造地物(如道路交叉口),作為控制點(diǎn)和驗(yàn)證點(diǎn)。使用Trimble Inpho軟件,導(dǎo)入控制點(diǎn),對(duì)GF-2衛(wèi)星異軌立體影像采用區(qū)域網(wǎng)平差的方式進(jìn)行高精度的幾何定位,其流程如圖4所示。 首先,通過(guò)特征匹配方法自動(dòng)提取大量均勻分布的連接點(diǎn),實(shí)現(xiàn)立體影像密集匹配,提高匹配的效率和可靠性,降低誤匹配率;通過(guò)自由網(wǎng)平差定向,對(duì)觀測(cè)像點(diǎn)的殘差進(jìn)行統(tǒng)計(jì)分析,進(jìn)而采用選權(quán)迭代法實(shí)現(xiàn)誤匹配連接點(diǎn)的自動(dòng)識(shí)別與剔除。然后,進(jìn)行控制點(diǎn)半自動(dòng)量測(cè),通過(guò)區(qū)域網(wǎng)平差定向和定向殘差分析,迭代計(jì)算得到區(qū)域網(wǎng)平差定向結(jié)果。最后,在自由網(wǎng)平差或區(qū)域網(wǎng)平差的基礎(chǔ)上,利用檢查點(diǎn)對(duì)立體影像無(wú)控制或有控制條件下的幾何定位誤差進(jìn)行統(tǒng)計(jì)分析。

圖4 DSM提取流程Fig.4 DSM extraction process

通過(guò)上述方法,在GF-2異軌像對(duì)重疊區(qū)共選取了13個(gè)控制點(diǎn)、78個(gè)連接點(diǎn),生成分辨率為4.88 m 的滑坡后DSM。將相同方法應(yīng)用在ZY-302的同軌像對(duì)上,共選取15個(gè)控制點(diǎn)、495個(gè)連接點(diǎn),生成分辨率為14.19 m的滑坡前DSM?？紤]到后續(xù)體積估算疊加分析,將滑坡前后的DSM統(tǒng)一重采樣為15 m,結(jié)果如圖5所示?；潞驞SM影像中東北角和清江南岸存在一系列噪聲點(diǎn),這是由于2020年12月22日的GF-2影像存在云層產(chǎn)生的測(cè)量誤差,不在滑坡區(qū)域內(nèi),對(duì)后續(xù)實(shí)驗(yàn)并無(wú)影響。

圖5 滑坡前后DSMFig.5 DSM before and after the landslide

2.2.4DSM校正與精度評(píng)價(jià)

目前的研究成果通常只涉及單一時(shí)間點(diǎn)的DSM提取,且研究區(qū)多以低緩丘陵為主,適合開(kāi)展實(shí)地精度控制工作。如胡芬等[14]利用異軌立體像對(duì)提取DSM時(shí),使用了地面控制點(diǎn)半自動(dòng)量測(cè),并通過(guò)增加控制點(diǎn)數(shù)量提高精度。相似的,倪文儉等[15]在使用高分2號(hào)提取DSM的過(guò)程中,利用無(wú)人機(jī)搭載高精度激光雷達(dá)提取了研究區(qū)的DSM作為驗(yàn)證數(shù)據(jù),該區(qū)域80%為坡度15°的緩坡。而本文需要對(duì)滑坡前后兩個(gè)時(shí)間點(diǎn)提取的DSM進(jìn)行差值計(jì)算,由于不同的影像數(shù)據(jù)源所生成的DSM結(jié)果存在高程基準(zhǔn)差異,有必要將其統(tǒng)一校正到參考DSM的基準(zhǔn)上。為此,本文提出在區(qū)域網(wǎng)平差的基礎(chǔ)上對(duì)滑坡前后的DSM均進(jìn)行二次高程校正,從而進(jìn)一步降低在復(fù)雜山地地形測(cè)量上出現(xiàn)的誤差。首先,采集一定數(shù)量的校正點(diǎn),以GLO30數(shù)據(jù)作為統(tǒng)一的參考DSM,使用最小二乘擬合二次多項(xiàng)式分別對(duì)滑坡前后的DSM進(jìn)行校正;然后再另采集一定數(shù)量的驗(yàn)證點(diǎn)對(duì)校正后的DSM做精度評(píng)價(jià)。所有校正點(diǎn)和驗(yàn)證點(diǎn)均以ESRI全球衛(wèi)星影像網(wǎng)絡(luò)服務(wù)為底圖參考,通過(guò)目視解譯,在地物變化較少、無(wú)植被覆蓋且遠(yuǎn)離滑坡崩塌區(qū)的硬質(zhì)水泥路面上選取,以減少因各種因素造成的地物變形而產(chǎn)生的高程偏差。

如圖6所示,本文共選擇了44個(gè)校正點(diǎn)和24個(gè)驗(yàn)證點(diǎn),主要均勻分布在研究區(qū)四周的道路,高程較為穩(wěn)定的區(qū)域。式(2)和式(3)分別為滑坡前后DSM與參考GLO30 數(shù)據(jù)的擬合方程:

圖6 校正點(diǎn)和驗(yàn)證點(diǎn)分布Fig.6 Layout of correction points and verification points

Y=0.9837X1+6.7636

(2)

Y=0.9816X2+11.2100

(3)

式中:Y為GLO30高程值,X1和X2分別為滑坡前和滑坡后DSM原始值,上述兩式中的一次項(xiàng)系數(shù)均十分接近于1,說(shuō)明滑坡前后的DSM與GLO30的高程差異基本為一個(gè)常數(shù)平面。

經(jīng)過(guò)高程校正后,滑坡前后的DSM高程中誤差分別為4.64 m和5.59 m,相較于校正前中誤差分別降低了5.64 m和2.45 m?？紤]到GLO30 數(shù)據(jù)在非極地地區(qū)的平均誤差為1.92 m[17],滑坡前后的區(qū)域DSM誤差基本<10 m,滿足《基礎(chǔ)地理信息數(shù)字成果1∶5 000、1∶10 000、1∶25 000、1∶50 000、1∶100 000數(shù)字高程模型》[21]規(guī)定的1∶10 000 DEM高山地的測(cè)圖要求,可用于滑坡體積的估算。

2.2.5滑坡物質(zhì)體積估算

滑坡發(fā)生后會(huì)形成滑坡凹地和滑坡堆積區(qū)。通過(guò)對(duì)滑坡后與滑坡前的高程差分,不但可以清晰地標(biāo)識(shí)出二者的位置和范圍,還能通過(guò)式(4)估算出滑坡物質(zhì)由凹地輸出或在堆積區(qū)匯集的體積。

(4)

式中:V為區(qū)域內(nèi)滑坡物質(zhì)體積總量,A,N和ΔDi分別為區(qū)域內(nèi)單個(gè)像元面積、像元總數(shù)和第i個(gè)像元的高程差。經(jīng)計(jì)算,沙子壩滑坡的凹地體積為237.01萬(wàn)m3,堆積區(qū)體積為137.24萬(wàn)m3,滑坡物質(zhì)體積變化分布情況如圖7所示,其中負(fù)值代表滑坡凹地,正值代表滑坡堆積區(qū)。

圖7 滑坡體積變化Fig.7 Landslide volume change

理論上,滑坡凹地輸出的物質(zhì)總量應(yīng)等于滑坡堆積區(qū)匯集的物質(zhì)總量,這種等量關(guān)系反映了“物質(zhì)平衡原理”[22],而從本文計(jì)算結(jié)果可看出,堆積區(qū)體積遠(yuǎn)小于凹地區(qū)體積,分析其原因:①該滑坡是典型的河岸滑坡,滑坡堆積區(qū)有很大一部分物質(zhì)已經(jīng)傾瀉入清江河道內(nèi),并逐漸累積形成滑坡壩,通過(guò)光學(xué)立體測(cè)量手段只能獲取滑坡壩露出水面的部分體積,因此出現(xiàn)堆積區(qū)體積偏低的情況;②GF-2影像拍攝時(shí)間距滑坡發(fā)生已間隔5個(gè)月以上,期間堆積物的清理工作已持續(xù)展開(kāi),大量堆積物已被人工清理,或被清江河水自然沖蝕搬運(yùn),致使測(cè)量獲得的堆積體積進(jìn)一步減少。

綜合上述分析,要獲得較準(zhǔn)確的滑坡物質(zhì)體積,應(yīng)從地面變化較小且較容易測(cè)量的滑坡凹地體積入手。在原凹地體積的基礎(chǔ)上,還應(yīng)考慮夾帶和刮擦效應(yīng),即滑動(dòng)物質(zhì)在路徑運(yùn)動(dòng)過(guò)程中會(huì)夾帶和液化飽和土壤,增大實(shí)際的滑坡物質(zhì)體積[8]。參考前人的研究成果[23-25],假設(shè)這一增量為25%,最終估算得到的滑坡物質(zhì)體積為296.28萬(wàn)m3。

3 結(jié)果分析與討論

根據(jù)圖7的體積變化分布結(jié)果,對(duì)沙子壩滑坡構(gòu)造特征進(jìn)行了劃分,其中:A為滑源區(qū),是滑坡的啟動(dòng)區(qū),也是主要的物質(zhì)輸出區(qū),投影面積為7.29萬(wàn)m2,平均高程差為-21.78 m,輸出物質(zhì)體積占總量的65.05%。B為滑動(dòng)帶,是滑坡物質(zhì)移動(dòng)的主要通道,后緣連接滑源區(qū),前緣位于清江北岸的高陡坡,前后緣高程差為300 m,長(zhǎng)約1 200 m,投影面積為24.88萬(wàn)m2,整體體積變化輕微,呈中部減少兩側(cè)增加的分布。其原因是:滑坡體在通過(guò)該區(qū)域的過(guò)程中發(fā)生了鏟刮和攜帶效應(yīng),裹挾位于滑動(dòng)帶中部的物質(zhì),造成中部體積損失;而部分位于滑坡體邊緣的物質(zhì)由于速度和能量的衰減,逐漸停滯并堆積,這些堆積物沿滑坡東西側(cè)壁分布,造成了該處體積增加。C為滑坡壩,也是本次滑坡的主要堆積區(qū),沿清江河道展布,長(zhǎng)約1 000 m。受衛(wèi)星攝影測(cè)量方法的局限性,圖7結(jié)果只反映了滑坡壩露出水面的體積,約為137.24萬(wàn)m3。D為陡坡,位于清江北岸,連接滑動(dòng)帶和滑坡壩。在已公開(kāi)的對(duì)本次滑坡研究的文獻(xiàn)中,對(duì)滑坡體積的估算主要采用以下兩種方法:①平均厚度乘以面積法,估算得到的滑坡體積約為1 000萬(wàn)m3[26];②現(xiàn)場(chǎng)估算法得到滑坡體積約為280萬(wàn)m3[27]。方法①實(shí)施最簡(jiǎn)單,但通常存在極大誤差,方法②則更具可信度,因此由方法②估算獲得的滑坡體積也被后續(xù)研究所采納。本文估算獲得的結(jié)果與方法②接近,偏差僅為5%,表明本文提出的滑坡體積估算方法不但具有較高的可行性和可靠的精度,且能直觀地反映滑坡體積變化的空間分布情況。

在方法設(shè)計(jì)上,本文所提出的數(shù)據(jù)處理思路可在數(shù)據(jù)缺乏、無(wú)實(shí)測(cè)控制點(diǎn)的情況下,充分挖掘已有數(shù)據(jù)資料,通過(guò)構(gòu)建適合的異軌影像對(duì),不依賴地面實(shí)測(cè)提取滑坡體積,可用于高山地等復(fù)雜地形的滑坡體積提取,也可進(jìn)一步推廣應(yīng)用于常規(guī)山地的高程測(cè)量。與此同時(shí),我們?nèi)詰?yīng)注意基高比對(duì)于異軌高程測(cè)量精度的影響。在數(shù)據(jù)選擇上應(yīng)優(yōu)先選擇大基高比的影像對(duì)。但由于衛(wèi)星的運(yùn)行軌道和載荷參數(shù)并未針對(duì)異軌測(cè)量而設(shè)計(jì)等多方面限制,大基高比的像對(duì)所覆蓋的公共區(qū)域通常較為狹窄,可能無(wú)法提取出足量的高質(zhì)量地面控制點(diǎn)。因此,還可采取適當(dāng)放寬基高比的策略,進(jìn)而依靠更高分辨率的影像,并結(jié)合本文提出的二次高程校正方法彌補(bǔ)基高比不足帶來(lái)的精度損失。

4 結(jié) 論

本文以如何充分利用衛(wèi)星遙感數(shù)據(jù)資源,并準(zhǔn)確提取滑坡物質(zhì)的體積為出發(fā)點(diǎn),以2020年7月發(fā)生在湖北省恩施市馬者村沙子壩的大型滑坡為研究對(duì)象,將異軌DSM提取方法應(yīng)用于傳統(tǒng)高分辨率光學(xué)衛(wèi)星GF-2影像,并利用其生成的滑坡后DSM與ZY-302衛(wèi)星生成的滑坡前同軌DSM進(jìn)行差分計(jì)算獲取了滑坡體積,結(jié)果表明:

(1) 經(jīng)過(guò)外部高程數(shù)據(jù)的二次校正精度后,異軌與同軌立體測(cè)量技術(shù)所生成的DSM精度相當(dāng),均滿足1∶10 000 DEM高山地的測(cè)圖要求?；虑昂驞SM與GLO30參考DSM均存在顯著的線性相關(guān)性,均滿足滑坡體積量算的精度要求。值得注意的是,GF-2全色影像的高分辨率,彌補(bǔ)了基高比偏低的不足,其精度接近于ZY-302同軌測(cè)量。

(2) 提取的滑坡體積為296.28萬(wàn)m3,與已公開(kāi)的文獻(xiàn)報(bào)道相吻合,具有較高的準(zhǔn)確性。在生成的滑坡體積變化分布圖中,清晰地呈現(xiàn)此次滑坡的構(gòu)造特征,可識(shí)別出滑源區(qū)、滑動(dòng)帶、陡坡和滑坡壩4部分。

本文提出基于異軌立體測(cè)量技術(shù)估算滑坡物質(zhì)體積的方法,充分利用已有衛(wèi)星遙感資料,以較低成本、較便捷的途徑,實(shí)現(xiàn)對(duì)滑坡地形的高精度遙感測(cè)量,對(duì)常規(guī)山地地形測(cè)繪具有參考價(jià)值;通過(guò)更精細(xì)地描繪出滑坡構(gòu)造,對(duì)進(jìn)一步認(rèn)識(shí)滑坡機(jī)理具有科學(xué)意義。

致 謝

本次研究在遙感數(shù)據(jù)獲取和處理工作中得到了廣西壯族自治區(qū)自然資源遙感院的大力支持和協(xié)助,在此表示感謝。