Stacking InSAR和IPTA技術(shù)應(yīng)用于探測寧夏隆德縣滑坡隱患

王 川，涂 寬，諶 華，耿 丹，王文龍，李樵民

（1.二十一世紀(jì)空間技術(shù)應(yīng)用股份有限公司，北京 100096；2.寧夏回族自治區(qū)遙感調(diào)查院，寧夏銀川 750021；3.高分辨率對地觀測系統(tǒng)寧夏數(shù)據(jù)與應(yīng)用中心，寧夏銀川 750021）

引言

滑坡是在重力、降雨、地震等因素的作用下，坡體作剪切運動的自然現(xiàn)象，通?？蓜澐譃橥埔剖交?、土流型滑坡、土質(zhì)滑坡和古崩滑坡堆積體等［1-2］。其產(chǎn)生條件包含內(nèi)因和外因兩類，內(nèi)因為地質(zhì)構(gòu)造和地貌特征，外因則為降雨、地震和人類工程活動等。近年來由于全球氣候波動，異常天氣條件增多，導(dǎo)致更多滑坡災(zāi)害在各地發(fā)生［2-4］。

由于滑坡災(zāi)害的發(fā)生嚴(yán)重威脅了人類的生命和財產(chǎn)安全，對人類社會環(huán)境、經(jīng)濟(jì)環(huán)境、生態(tài)環(huán)境造成了巨大破壞，因此需要對滑坡隱患進(jìn)行及時的探測和監(jiān)測。通常，對滑坡進(jìn)行探測和監(jiān)測的手段有：常規(guī)三角網(wǎng)測量、伸縮計深部位移監(jiān)測、應(yīng)力應(yīng)變測量、GPS遙測監(jiān)測網(wǎng)、地震監(jiān)測、近景攝影測量以及多種技術(shù)的聯(lián)合等［2，5］。上述探測和監(jiān)測手段的準(zhǔn)確性高，測量結(jié)果可靠，但測量主要集中在小區(qū)域，不利于開展大規(guī)?；码[患探測。

近年來，隨著科技的不斷進(jìn)步，出現(xiàn)了以測量可見光光譜、熱輻射、微波輻射等為目標(biāo)的多種遙感技術(shù)［6-7］。特別是合成孔徑雷達(dá)（Synthetic Aperture Radar，SAR）技術(shù)的出現(xiàn)。微波遙感干涉合成孔徑雷達(dá)（Interferometry Synthetic Aperture Radar，InSAR）測量地表形變技術(shù)測量結(jié)果準(zhǔn)確度高，覆蓋面大，在滑坡隱患探測中得到了普遍應(yīng)用［8］。最先發(fā)展起來的應(yīng)用電磁波相位值對地表起伏進(jìn)行量測的是差分干涉技術(shù)（Differential Interferometry Synthetic Aperture Radar，DInSAR），但是單獨使用兩景影像進(jìn)行干涉易產(chǎn)生各種去相干（時間、空間、體）等，且易受大氣、系統(tǒng)熱噪聲、處理誤差等影響，因此目前主流技術(shù)是多個干涉對共同測 量 的 時 序InSAR技 術(shù)（Multi-Temporal Interferometry Synthetic Aperture Radar，MT-InSAR）［2，9］。在MT-InSAR中，典型的是：（1）永久散射體技術(shù)（Persistent Scatter InSAR，PS-InSAR），該方法只選擇回波信號穩(wěn)定的點（永久散射體）進(jìn)行分析，可以擺脫時空基線的限制，相干性保持好；（2）小基線集技術(shù)（Small Baseline Subsets InSAR，SBAS-InSAR），該技術(shù)選擇分布式散射體（如沙地、草地等）進(jìn)行分析，更適用于自然環(huán)境［10-12］。以此兩種思想為核心，產(chǎn)生了相干點目標(biāo)分析（Interferometry Point IPTA）、SqueeSAR等技術(shù)，同時也出現(xiàn)了以計算形變速率為目標(biāo)的干涉疊加技術(shù)（Stacking）和不依賴時序相位特征測量地表形變的多孔徑雷達(dá)技術(shù)（Multiple Aperture Interferometry，MAI）［13-19］。InSAR技術(shù)不可避免的會受到植被、地形起伏和大氣等的影響，使用單一技術(shù)往往只能克服部分缺陷而導(dǎo)致不能更全面的識別滑坡隱患，同時基于光學(xué)影像對滑坡隱患進(jìn)行專家判讀缺少坡體滑移的數(shù)據(jù)支持，部分老（古）滑坡已趨于穩(wěn)定。因此需要聯(lián)合多種InSAR方法并結(jié)合光學(xué)影像進(jìn)行滑坡隱患探測。

IPTA技術(shù)針對點目標(biāo)進(jìn)行回歸分析以分離相位（大氣、地形、非線性形變等），計算效率和精度顯著提高，而Stacking技術(shù)利用加權(quán)平均以抑制噪聲使其能夠快速獲取大范圍的形變信息，在現(xiàn)有各類InSAR技術(shù)中，此2種技術(shù)可操作性最強，同時結(jié)果可互相補充。隆德縣境內(nèi)以黃土丘陵區(qū)為主，較少出現(xiàn)大型滑坡隱患，小型滑坡隱患分布廣泛，利用Stacking技術(shù)能夠?qū)崿F(xiàn)對大部分滑坡隱患的探測，基于IPTA技術(shù)能夠去除各種誤差，探測滑坡隱患的可靠性高。

因此，本研究利用Sentinel-1衛(wèi)星影像，采用IPTA和Stacking這2種InSAR技術(shù)對隆德縣全境滑坡隱患進(jìn)行探測，分析和評價2種InSAR技術(shù)的探測結(jié)果。該項工作對隆德縣滑坡隱患普查及監(jiān)測具有十分重要的指導(dǎo)意義。

1 研究區(qū)概況及研究數(shù)據(jù)

研究區(qū)隆德縣位于寧夏回族自治區(qū)的南部，六盤山西側(cè)，全縣總面積992.39 km2。在地形地貌特征上，隆德縣處于昆侖秦嶺地槽褶皺區(qū)東段，曾受到燕山期、喜馬拉雅期構(gòu)造運動的影響，其主要地貌類型有：紅層巖丘陵地貌、中低山地貌、河谷平原地貌、黃土丘陵地貌。在氣候特征上，隆德縣屬于溫帶大陸性季風(fēng)氣候，受東高西低的地勢差異影響，致使東部濕潤寒冷而西部干燥溫暖，冬季干燥而夏季多雨。每年7～9月降雨最為集中，約占全年降水量的60%。降雨和降雪是隆德縣地質(zhì)災(zāi)害發(fā)生的重要誘因。由于人類工程活動的影響，如毀林開荒、削坡修路、興修水利和過度放牧等，破壞了生態(tài)系統(tǒng)平衡，加劇了地質(zhì)環(huán)境惡化，更進(jìn)一步的促使了滑坡災(zāi)害的發(fā)生。

隆德縣滑坡包含黃土滑坡和巖質(zhì)滑坡2種，其中黃土滑坡包含黃土層內(nèi)滑坡和黃土—基巖滑坡，巖質(zhì)滑坡主要是新近系泥巖殘坡積物滑坡?；乱?guī)模通常為小型滑坡，中、大型滑坡發(fā)育較少。滑坡和滑坡隱患多分布于渝河、南川河的沖蝕岸邊及其支流或支溝兩側(cè)，呈條帶狀或樹枝狀分布，在地貌類型上多集中在丘陵地貌區(qū)［20］?；乱詼\層和中層為主，厚層滑坡較少。

Sentinel-1衛(wèi)星影像幅寬廣（IW成像模式可達(dá)250 km），重訪周期相對較短且固定（單星重訪周期12 d），因此文中研究選擇2019年1月至2020年5月在重訪周期過境的全部衛(wèi)星影像用于測量地表形變，表1詳細(xì)列出了所使用的衛(wèi)星影像參數(shù)信息。圖1為研究區(qū)概況圖和Sentinel-1衛(wèi)星影像覆蓋范圍。此外，為了糾正地形影響和提取地形信息，文中研究還同時獲取了SRTM 30 m的DEM作為輔助數(shù)據(jù)。

表1 Sentinel-1衛(wèi)星影像參數(shù)表Table 1 Information of Sentinel-1

注：在Sentinel-1假彩色合成圖中，綠色區(qū)域反映出地物主要由體散射比較強烈的散射體如植被構(gòu)成，深粉色區(qū)域反映出地物主要由粗糙的表面如裸地構(gòu)成，淺粉色區(qū)域表示地物由去極化能力非常弱且后向散射強烈的建筑物組成，水體由于鏡面散射而呈現(xiàn)黑色。

本次實驗所獲取的SAR數(shù)據(jù)能夠完整覆蓋2019年全年，時間和空間基線滿足形變測量需求，配準(zhǔn)精度達(dá)到InSAR干涉測量要求（1/8像元精度），DEM數(shù)據(jù)空間分辨率與SAR數(shù)據(jù)接近且已對空值進(jìn)行插值，能夠利用全部衛(wèi)星影像和DEM數(shù)據(jù)開展后續(xù)研究工作。

2 基于Stacking和IPTA技術(shù)的滑坡體探測方法

2.1 InSAR成像原理

根據(jù)圖2幾何關(guān)系，P點的基線分解為：

同樣，參考平面的點P0的基線分解為：

假設(shè)2次觀測時沒有地表起伏變化，同時忽略地物散射和輻射傳輸中產(chǎn)生的隨機噪聲，則P點的干涉相位可以視作2次傳輸路徑之差所對應(yīng)的相位：

對路徑差ΔR求導(dǎo)，可知相位隨路徑的變化關(guān)系為：

對于星載SAR衛(wèi)星，當(dāng)2次觀測視角的差異較小時，基于微波輻射傳輸?shù)倪h(yuǎn)場平行近似［21］，可以將ΔR等效替換為水平基線B∥，當(dāng)觀測距離不變且視角由θ0開始變化產(chǎn)生的干涉相位變化為：

若定義參考面點P0的干涉相位為參考相位，其與P點到傳感器S1的距離相同，因此可以視其為P點的平地相位（排除地形起伏的影響，僅由觀測視角θ0造成）：

根據(jù)圖像的幾何關(guān)系：

則視角由θ0開始變化產(chǎn)生的高程變化：

結(jié)合式（7）和式（10）可得：

圖2 InSAR原理圖Fig.2 InSAR schematic diagram

此處??if的變化僅與高程相關(guān)（P0為參考面起始觀察點），即由P和P0點的干涉相位差分得到，所以稱其為高程相位?h（忽略大氣、形變等噪聲相位）。

若P點在觀測期間產(chǎn)生了位移（地表形變），則干涉相位可表示為：

式中：?d表示地表形變產(chǎn)生的相位，與之對應(yīng)的雷達(dá)視線向的形變量為ΔRd，由地表形變相位與視線向形變量的關(guān)系可以看出，對于C波段SAR衛(wèi)星，其地表形變的測量精度通常能達(dá)到毫米級。

2.2 研究方法

本研究選擇Stacking和IPTA這2種較為成熟的InSAR方法用于測量隆德縣境內(nèi)的地表形變值。Stacking技術(shù)采用對干涉相位加權(quán)平均的方式計算相位的形變速率繼而得到地表延視線向的形變速率，其處理方法相對簡單且能夠在一定程度上對大氣相位等噪聲進(jìn)行抑制，因此能夠快速的對大區(qū)域在長時間序列上得到地表形變速率的測量結(jié)果［22］。IPTA方法是一種基于點目標(biāo)的干涉測量方法，通過迭代回歸分析將大氣相位、高程相位殘余等進(jìn)行分離，能夠同時獲取地表形變速率和地表非線性形變，與SBAS和PS等方法相比具有計算效率高、節(jié)省存儲空間的特點［23］。

在進(jìn)行InSAR處理前，需要首先對Sentinel-1衛(wèi)星影像進(jìn)行預(yù)處理，主要包括：讀取原始數(shù)據(jù)并基于外部DEM數(shù)據(jù)對SLC圖像進(jìn)行配準(zhǔn)，將配準(zhǔn)后的圖像裁剪至研究區(qū)范圍，對DEM進(jìn)行重投影至SAR坐標(biāo)系。兩種InSAR方法的主要處理流程如下：

2.2.1 Stacking InSAR技術(shù)流程

（1）計算各個像對間的垂直基線，設(shè)置閾值為-200 m～200 m，計算時間基線，設(shè)置閾值50 d，生成基線連接圖和干涉對列表（圖3）。

（2）根據(jù)干涉對列表對SLC圖像進(jìn)行差分干涉，對差分干涉圖進(jìn)行濾波以提高相干性，對干涉圖進(jìn)行解纏，對解纏圖進(jìn)行一定的空間濾波部分消除大氣、軌道誤差、高程誤差等噪聲相位。解纏方法：最小費用流法（MCF）［24］。

（3）基于加權(quán)平均的方法計算相位形變速率：

式中：φrate為相位形變速率；Δtj表示第j個干涉相對的時間間隔；N表示干涉像對數(shù)，然后按照λ/4π將相位形變速率換算為延視線向的地表形變速率（單位：mm/年）。由于單軌干涉測量存在著多種去相干（如時間、幾何、體等），在解纏過程中需要對相干性較低的像元進(jìn)行掩膜（文中研究相干性閾值為0.2），使得在進(jìn)行Stacking時，各像素點所使用的樣本數(shù)（干涉像對數(shù)）均不同，為了使結(jié)果更具可靠性，文中研究設(shè)定最小可接受樣本數(shù)為108，即當(dāng)像素點干涉像對數(shù)低于108則對結(jié)果進(jìn)行掩膜。圖4為Stacking InSAR總體處理技術(shù)流程圖。

圖3 Sentinel-1衛(wèi)星影像基線連接圖Fig.3 Sentinel-1 perpendicular base connection diagram

圖4 Stacking InSAR處理技術(shù)流程圖Fig.4 The flow chart of Stacking InSAR

2.2.2 IPTA技術(shù)流程

（1）分別基于像素點在時序上的穩(wěn)定性和相干性挑選像元并生成點目標(biāo)文件。計算兩兩像對間的垂直基線，設(shè)置閾值為-200～200 m，計算時間基線，設(shè)置閾值50 d，生成基線連接圖和干涉對列表（圖3）。根據(jù)干涉對列表對點目標(biāo)進(jìn)行差分干涉。

（2）在第j幅差分干涉圖，點目標(biāo)i其相位可以表示為：

圖5 IPTA處理技術(shù)流程圖Fig.5 The flow chart of IPTA

式中：Δ?表示差分干涉相位；vi表示第i個點目標(biāo)的地表形變速率；Δtj表示第j個干涉對的時間間隔；δi,j表示第i個點目標(biāo)在Δtj時間段內(nèi)的非線性地表形變值；B⊥j表示第j個干涉像對的垂直基線；Δzi表示第i個點目標(biāo)的高程誤差；λ表示波長；ρi表示微波從雷達(dá)到達(dá)第i個點所對應(yīng)的地物再返回的路徑（斜距）；θi表示入射角；?atm表示大氣相位；n表示噪聲相位。通過非線性回歸分析分別求解出其中的差分大氣相位和高程誤差，對地表高程進(jìn)行補償并對差分大氣相位進(jìn)行去除后再次進(jìn)行回歸分析，進(jìn)行多次迭代后可將差分大氣相位逐步分離，并對地表高程進(jìn)行充分補償。

（3）通過回歸分析獲得解纏的差分干涉相位后，使用奇異值分解（SVD）的方法對多主影像的相位進(jìn)行求解計算出延時間序列的相位形變［25］，對其在空間上和時間上進(jìn)行濾波，并換算成地表沉降值（單位：mm）。在得到時序相位形變后，再次進(jìn)行回歸分析計算出相位的形變速率，同樣換算為地表形變速率（單位：mm/年）。圖5為IPTA總體處理技術(shù)流程圖。

2.2.3 有效探測滑坡隱患和未有效探測滑坡隱患的定義及有效探測率

在開展文中研究工作之前，已先期對隆德縣滑坡隱患進(jìn)行了調(diào)查，此次工作聯(lián)合多源遙感技術(shù)對隆德縣全境滑坡隱患進(jìn)行了普查，在此基礎(chǔ)上對所識別出的滑坡隱患已完成了實地的調(diào)查和驗證工作。該項工作由“寧夏科技廳重點研發(fā)計劃項目”支持開展［26］。

根據(jù)InSAR地表形變測量結(jié)果識別滑坡隱患，若圈定的滑坡隱患與隆德縣滑坡隱患遙感監(jiān)測工作成果中的滑坡隱患相對應(yīng)，則稱其為該InSAR技術(shù)（Stacking和IPTA）的“有效探測滑坡隱患”，對于使用該InSAR測量方法未解譯出但在隆德縣滑坡隱患遙感監(jiān)測工作成果中出現(xiàn)的其他滑坡隱患，則稱其為“未有效探測滑坡隱患”。對于一種InSAR技術(shù)（IPTA或Stacking），其有效探測滑坡隱患占所有滑坡隱患的百分比文中研究稱其為“有效探測率”。

3 Stacking和IPTA滑坡隱患探測結(jié)果

通過與隆德縣滑坡隱患調(diào)查結(jié)果進(jìn)行核實，2種InSAR方法對47處滑坡隱患進(jìn)行了有效探測［26］（圖6，圖7（a）、（b）），在形變測量中正值表示滑坡隱患在延視線向靠近衛(wèi)星滑動（滑坡隱患和衛(wèi)星間的距離縮短），負(fù)值表示滑坡隱患在延視線向遠(yuǎn)離衛(wèi)星滑動（滑坡隱患和衛(wèi)星間的距離增大）。表2詳細(xì)列出了滑坡隱患的基本信息，其中滑坡隱患面積達(dá)到2.0×106m2的共6處，最大滑坡隱患面積6.2×106m2。

Stacking技術(shù)獲得的形變結(jié)果中有效探測滑坡隱患44處，有效探測率93.6%，IPTA有效探測滑坡隱患共計37處，有效探測率78.7%，Stacking技術(shù)略高。2種InSAR技術(shù)共有效探測滑坡隱患34處，在相同數(shù)據(jù)源情況下監(jiān)測重合度高，對最大形變速率進(jìn)行相關(guān)性分析（圖8），線性相關(guān)性強（R20.56），測量結(jié)果一致性強（例如圖7（c）和（d）所示滑坡隱患），可相互印證2種InSAR技術(shù)測量結(jié)果可靠。

圖6 隆德縣境內(nèi)滑坡隱患分布圖Fig.6 Distribution map of potential landslide in Longde County

圖7 隆德縣Stacking和IPTA測量形變結(jié)果圖及典型形變滑坡隱患Fig.7 The displacement map measured by Stacking and IPTA and the typical potential landslide in Longde County

表2 隆德縣解譯滑坡隱患基本信息表Table 2 Information of interpreting potential landslide in Longde County

續(xù)表

圖8 不同InSAR方法測量滑坡隱患最大形變量線性相關(guān)性分析Fig.8 Correlation analysis of the maximum displacement for different InSAR methods

圖9 有效探測滑坡隱患最大形變速率分布特征Fig.9 Characteristics of the maximum displacement of effectively detected potential landslides

3.1 Stacking和IPTA有效探測滑坡隱患結(jié)果分析

3.1.1 有效探測滑坡隱患的最大形變速率分布特征

圖9描述了滑坡隱患最大形變速率分布特征。箱形圖中間箱體上邊沿表示上四分位數(shù)（3/4分位數(shù)），下邊沿表示下四分位數(shù)（1/4分位數(shù)），上下四分位數(shù)之差為四分位距（IQR），箱體中間橫線表示中位數(shù)，箱體中圓圈表示均值，由箱體延伸出上下邊緣線（箱體外橫線），上邊緣線為上四分位數(shù)加1.5倍四分位距，下邊緣線為下四分位數(shù)減1.5倍四分位距，超過上下邊緣線的視為異常值。2種InSAR方法對不同形變速率滑坡隱患探測能力基本一致。其中一半的滑坡隱患最大形變速率集中在-20～20 mm/年之間（箱的上邊沿和下邊沿分別表示3/4和1/4分位數(shù)）。Stacking所識別滑坡隱患形變速率中位數(shù)（箱內(nèi)橫線）高于IPTA。

3.1.2 有效探測滑坡隱患的地形分布特征

研究主要針對地形中的坡度和坡向進(jìn)行分析，其中坡度劃分為平坡、緩坡、斜坡和陡坡4個等級（表3）［27］。坡向劃分為北、東北、東、東南、南、西南、西和西北8個方向（表4）［27］。IPTA和Stacking有效探測滑坡隱患均集中在緩坡和斜坡上，對于坡度較大的平坡（1.7°～5°）和坡度較小的陡坡（25°～31.6°）亦有少量分布（圖10（a）），這與InSAR的識別原理相符合，當(dāng)坡度較小時滑移速率較慢無法滿足InSAR測量精度，坡度過大則滑移速率過快易產(chǎn)生失相干導(dǎo)致無法測量。北坡方向未出現(xiàn)有效探測滑坡隱患，與北坡所受光照輻射少其微氣候環(huán)境不易出現(xiàn)滑坡相關(guān)，在其他方向上未出現(xiàn)明顯集中分布特征（圖10（b）），因此綜合2種InSAR方法能夠較全面的對各個方向的滑坡隱患進(jìn)行有效探測。

表4 地形-坡向統(tǒng)計結(jié)果Table 4 Statistics of topographic aspects

表3 地形-坡度統(tǒng)計結(jié)果Table 3 Statistics of topographic slopes

坡度級劃分來源于“基于3PGS-MTCLIM模型模擬根河林區(qū)火后植被凈初級生產(chǎn)力恢復(fù)及其影響因子”。

圖10 滑坡隱患地形分布特征Fig.10 Characters of topography of potential landslides

3.2 滑坡隱患滑動特征分析

圖11 IPTA有效探測滑坡隱患的時序形變量Fig.11 The displacement of potential landslides detected by IPTA method

IPTA技術(shù)不僅能夠獲取滑坡隱患的年平均形變速率，同時能夠獲取時序累計形變量，因此可以通過IPTA技術(shù)了解滑坡隱患的滑動狀態(tài)。文中研究對IPTA獲取的37處滑坡隱患的滑動過程做一基本描述（圖11）。

LDhp1、LDhp3、LDhp4、LDhp5、LDhp9、LDhp13、LDhp26、LDhp28、LDhp37、LDhp44均以一定的滑動速率發(fā)生線性形變，滑動過程在整個監(jiān)測期內(nèi)較穩(wěn)定，沒有表現(xiàn)出受其他因素干擾而滑動過程受到擾動的情況。

LDhp2、LDhp6、LDhp8、LDhp11、LDhp12、LDhp18、LDhp22、LDhp30、LDhp32、LDhp43在前期滑動不明顯或者滑動速率保持穩(wěn)定，在經(jīng)過一段時間后滑動速度突然加劇并且會維持至監(jiān)測結(jié)束，這一加速過程發(fā)生在2019年6月～9月時間段內(nèi)，這一時間段是一年中降水集中時期，因此推測是由于降雨沖刷坡體導(dǎo)致滑動過程加劇。

LDhp7、LDhp10、LDhp14、LDhp16、LDhp24、LDhp25、LDhp27、LDhp29、LDhp39同樣在2019年6月～9月出先了較明顯的加速過程，但隨后在2019年10月～12月份加速過程突然減緩進(jìn)入新的相對慢速滑動的過程，這一特征應(yīng)該與該地區(qū)溫度降低導(dǎo)致土層的穩(wěn)定性發(fā)生改變相關(guān)，其滑動過程是否隨季候變化呈現(xiàn)周期性還需要在更長時間周期內(nèi)監(jiān)測進(jìn)行驗證。

LDhp15、LDhp17、LDhp19、LDhp21、LDhp23、LDhp38、LDhp42、LDhp47其最主要的特征是在2019年6月～12月相對于其他時間段始終保持緩慢滑動的狀態(tài)，而在2019年1月～5月份或是2020年1月～5月份則出現(xiàn)了明顯加速滑動的過程，這種在非生長季滑動過程加快應(yīng)于坡體上植被的季候變化相關(guān)。

4 結(jié)論

文中研究基于Sentinel-1衛(wèi)星影像進(jìn)行Stacking和IPTA處理，對隆德縣境內(nèi)的滑坡隱患進(jìn)行探測，核實探測到的滑坡隱患共計47處，其中，Stacking技術(shù)有效探測滑坡隱患44處，IPTA有效探測滑坡隱患37處，2種InSAR技術(shù)共同識別34處?；?種InSAR技術(shù)結(jié)果：

（1）Stacking和IPTA技術(shù)所探測到的滑坡隱患最大形變速率線性相關(guān)性強，形變測量結(jié)果可相互印證，可以認(rèn)為基于InSAR技術(shù)進(jìn)行滑坡隱患的探測具有較高的可靠性。

（2）由于IPTA需要選擇高質(zhì)量的點目標(biāo)從而導(dǎo)致部分滑坡隱患無點目標(biāo)覆蓋，而使用Stacking方法能夠?qū)@些滑坡隱患進(jìn)行有效識別。點目標(biāo)的缺失與野外環(huán)境復(fù)雜有關(guān)，如植被等地表覆蓋物造成相干性降低，地形起伏在局部地區(qū)出現(xiàn)大量地形相關(guān)的大氣相位。后續(xù)應(yīng)聯(lián)合多種InSAR技術(shù)，如采用DS InSAR方法增加同質(zhì)點的選?。?8］，采用BM3D進(jìn)行更好的濾波降噪，提高IPTA分析的點目標(biāo)覆蓋范圍和測量精度，從而提高滑坡隱患探測的可靠性。

（3）在形變速率上，隆德縣以黃土丘陵區(qū)為主要地貌特征，較少存在劇烈的坡體滑動過程，因此Stacking和IPTA技術(shù)能夠?qū)υ摰貐^(qū)不同形變速率的滑坡隱患進(jìn)行有效探測。同樣由于IPTA的選點機制導(dǎo)致滑動較的坡體候選點減少，其所能識別的滑坡隱患形變速率中位數(shù)略低于Stacking技術(shù)。

（4）在地形特征上，Stacking和IPTA所探測到的滑坡隱患主要集中在緩坡和斜坡上，同時除北坡不易出現(xiàn)滑坡外，在各個坡向的滑坡隱患均可探測。

（5）基于IPTA獲取的時序形變量表明，隆德縣境內(nèi)的滑坡隱患存在4種典型滑動類型：（a）以穩(wěn)定的速率持續(xù)滑動；（b）在2019年6月～9月份出現(xiàn)明顯加速滑動過程并持續(xù)到監(jiān)測結(jié)束；（c）在2019年6月～9月份出現(xiàn)明顯加速滑動過程，10月～12月份滑動過程逐漸減緩；（d）在2019年或2020年1月～5月滑動過程明顯，其它時間坡體較穩(wěn)定。

總之，綜合Stacking和IPTA這2種InSAR技術(shù)能夠全面的識別隆德縣境內(nèi)的滑坡隱患，在后續(xù)滑坡隱患探測中應(yīng)聯(lián)合多種InSAR方法，以獲取更可靠的結(jié)果。下一步將首先基于GNSS等技術(shù)對Stacking InSAR和IPTA技術(shù)得到的平均形變速率和時序形變結(jié)果進(jìn)行精度評價，并聯(lián)合其他InSAR技術(shù)提高IPTA的探測能力。