牛宗碩

（昆明理工大學(xué)國(guó)土資源工程學(xué)院，云南昆明 650093）

?

衛(wèi)星雷達(dá)數(shù)據(jù)與處理軟件在德欽縣地質(zhì)災(zāi)害監(jiān)測(cè)中的最優(yōu)選取研究

牛宗碩

（昆明理工大學(xué)國(guó)土資源工程學(xué)院，云南昆明 650093）

摘 要：德欽縣地質(zhì)災(zāi)害監(jiān)測(cè)項(xiàng)目主要運(yùn)用的技術(shù)手段為合成孔徑雷達(dá)干涉技術(shù)，簡(jiǎn)稱InSAR，該技術(shù)具有監(jiān)測(cè)范圍廣、不受時(shí)間空間限制、監(jiān)測(cè)精度高等特點(diǎn)。當(dāng)前，雷達(dá)衛(wèi)星和處理軟件多種多樣，根據(jù)德欽縣地質(zhì)環(huán)境特征，選取合適的雷達(dá)衛(wèi)星和處理軟件，是保證地質(zhì)項(xiàng)目得以順利實(shí)施，取得良好監(jiān)測(cè)結(jié)果的重要環(huán)節(jié)。本文通過(guò)對(duì)比覆蓋德欽縣的雷達(dá)衛(wèi)星數(shù)據(jù)，探討數(shù)據(jù)處理軟件的優(yōu)缺點(diǎn)，采用最優(yōu)原則，分析適合德欽縣地質(zhì)災(zāi)害監(jiān)測(cè)的雷達(dá)衛(wèi)星和處理軟件。

關(guān)鍵詞：地質(zhì)災(zāi)害；雷達(dá)衛(wèi)星；合成孔徑雷達(dá)干涉技術(shù)；差分雷達(dá)干涉技術(shù)

文獻(xiàn)著錄格式:牛宗碩.衛(wèi)星雷達(dá)數(shù)據(jù)與處理軟件在德欽縣地質(zhì)災(zāi)害監(jiān)測(cè)中的最優(yōu)選取研究［J］.浙江農(nóng)業(yè)科學(xué)，2016，57（6）: 926-930.

德欽縣位于青藏高原南延的橫斷山脈縱谷部位，境內(nèi)高山峽谷，高差較大；冰峰林立，雪山縱橫；構(gòu)造復(fù)雜，地形陡峻。特殊的地質(zhì)環(huán)境特征給該區(qū)帶來(lái)了嚴(yán)重的地質(zhì)災(zāi)害，主要有泥石流、滑坡、崩塌等。地質(zhì)災(zāi)害給德欽縣人民造成了極大的威脅，嚴(yán)重制約了城市建設(shè)及國(guó)民經(jīng)濟(jì)的發(fā)展。德欽縣地質(zhì)災(zāi)害監(jiān)測(cè)項(xiàng)目主要運(yùn)用合成孔徑雷達(dá)干涉技術(shù)（InSAR），為保證項(xiàng)目實(shí)施效果，需遵從“最優(yōu)原則”，對(duì)覆蓋德欽縣的雷達(dá)衛(wèi)星和數(shù)據(jù)處理軟件進(jìn)行選擇。通過(guò)對(duì)比多種雷達(dá)衛(wèi)星，分析雷達(dá)衛(wèi)星數(shù)據(jù)的可用性和最優(yōu)性，選擇適合德欽縣的雷達(dá)衛(wèi)星；通過(guò)對(duì)比多款數(shù)據(jù)處理軟件，分析軟件的可行性和實(shí)用性，選擇適合覆蓋德欽縣的雷達(dá)衛(wèi)星數(shù)據(jù)的處理軟件。

1 在軌民用雷達(dá)衛(wèi)星簡(jiǎn)介

雷達(dá)衛(wèi)星是遙感衛(wèi)星的一種，屬于微波遙感衛(wèi)星，技術(shù)原理為雷達(dá)干涉原理。雷達(dá)衛(wèi)星的發(fā)展始于20世紀(jì)70年代，由美國(guó)率先進(jìn)行發(fā)射和研究，主要用來(lái)進(jìn)行軍事偵察和外星探索。隨著該技術(shù)的不斷成熟和普及，越來(lái)越多的國(guó)家研究并發(fā)射了雷達(dá)衛(wèi)星，雷達(dá)衛(wèi)星的應(yīng)用也從單一的軍事機(jī)密領(lǐng)域發(fā)展到現(xiàn)在的多個(gè)領(lǐng)域。

目前，主要的在軌民用雷達(dá)衛(wèi)星技術(shù)參數(shù)詳見(jiàn)表1。除表中列出的以外，擁有雷達(dá)衛(wèi)星的國(guó)家還包括美國(guó)、俄羅斯、中國(guó)、印度、以色列等，由于這些國(guó)家的雷達(dá)衛(wèi)星數(shù)據(jù)只用于軍事或者只限本國(guó)使用，所以相關(guān)衛(wèi)星參數(shù)不對(duì)外公布，因此，不將其作為德欽縣InSAR數(shù)據(jù)選擇的備選衛(wèi)星。

表1　在軌民用雷達(dá)衛(wèi)星技術(shù)參數(shù)

2 德欽縣雷達(dá)衛(wèi)星數(shù)據(jù)獲取

2.1 德欽縣雷達(dá)數(shù)據(jù)采集要求

根據(jù)德欽縣城周邊地質(zhì)災(zāi)害可行性研究報(bào)告的要求，InSAR數(shù)據(jù)采集應(yīng)有以下幾點(diǎn)要求:監(jiān)測(cè)范圍為德欽縣城周邊約60 km2；監(jiān)測(cè)時(shí)間為2012—2016年，約5期數(shù)據(jù)，其中，在2013年8月28日的地震前后分別獲取2期歷史存檔數(shù)據(jù)；空間分辨率要求≤10 m。

2.2 衛(wèi)星參數(shù)對(duì)地質(zhì)災(zāi)害監(jiān)測(cè)項(xiàng)目的影響

在衛(wèi)星選取過(guò)程中，除了考慮德欽縣地質(zhì)災(zāi)害監(jiān)測(cè)的要求外，還需要對(duì)衛(wèi)星自身的參數(shù)進(jìn)行考慮。

空間分辨要求≤10 m，這是可研報(bào)告的要求。一般發(fā)射時(shí)間較早的衛(wèi)星精度較低，如歐洲的Envisat衛(wèi)星，其最大的空間分辨率為10 m，加上其他因素的影響，定會(huì)造成精度的降低，達(dá)不到技術(shù)要求，而且該衛(wèi)星發(fā)射時(shí)間早，技術(shù)相對(duì)落后，不是最佳的衛(wèi)星選擇。

波段波長(zhǎng)是影響監(jiān)測(cè)結(jié)果的一個(gè)重要因素，衛(wèi)星發(fā)射的波段越長(zhǎng)，其穿透能力越好。德欽縣是一個(gè)低緯度、高海拔的地區(qū)，夏季植被發(fā)育茂密、多云多霧多雨，而且地形多高山峽谷的特征均會(huì)導(dǎo)致微波的干涉效果欠佳，所以要盡量選擇波長(zhǎng)較長(zhǎng)的衛(wèi)星，如日本的ALOS衛(wèi)星，采用L波段，波長(zhǎng)為一般衛(wèi)星的5～10倍。

中央入射角、重復(fù)周期、極化方式是影響編程數(shù)據(jù)的重要因素。中心入射角越大、越接近正射，干涉效果就越好。重復(fù)周期越短，就可以獲得越多的景數(shù)據(jù)，利于差分處理。進(jìn)行干涉差分處理時(shí)，極化方式通常要選擇HH方式或者VV方式。

歷史存檔數(shù)據(jù)是制約項(xiàng)目數(shù)據(jù)選取的重要因素，由于德欽縣地質(zhì)災(zāi)害監(jiān)測(cè)項(xiàng)目要求在2013年8月28日前后要有2期數(shù)據(jù)，這就要求衛(wèi)星在該時(shí)間點(diǎn)前后均要有歷史存檔數(shù)據(jù)，而且差分處理至少需要兩景存檔數(shù)據(jù)。據(jù)調(diào)查，記錄了德欽縣歷史存檔數(shù)據(jù)的雷達(dá)衛(wèi)星主要有意大利的Cosmo-Skymed衛(wèi)星、德國(guó)的TerraSAR-X衛(wèi)星、加拿大的Radarsat-2衛(wèi)星。

2.3 雷達(dá)衛(wèi)星的數(shù)據(jù)對(duì)比

根據(jù)以上分析，在考慮德欽縣地震時(shí)間點(diǎn)的情況下，意大利Cosmo-Skymed衛(wèi)星、德國(guó)的TerraSAR-X衛(wèi)星、加拿大的Radarsat-2衛(wèi)星被定為備選衛(wèi)星。下面，根據(jù)衛(wèi)星的詳細(xì)信息進(jìn)行備選分析。

德國(guó)的TerraSAR-X衛(wèi)星。該衛(wèi)星的歷史存檔數(shù)據(jù)為2011年，分辨率為16 m，不符合本項(xiàng)目精度的要求，所以德國(guó)的TerraSAR-X衛(wèi)星不適合本項(xiàng)目的選擇。

加拿大的Radarsat-2衛(wèi)星。如圖1所示，加拿大的Radarsat-2衛(wèi)星有2種不同分辨率的歷史存檔數(shù)據(jù)，左側(cè)為3 m分辨率的3景數(shù)據(jù)，右側(cè)為8 m分辨率的5景數(shù)據(jù)。要想得到效果較好的差分干涉圖，要求衛(wèi)星的采集模式、中央入射角度、分辨率、軌道、極化方式等參數(shù)保持一致。但從圖1可以看出，無(wú)論是3 m還是8 m分辨率的數(shù)據(jù)，它們的采集情況都不相同（相同的話則兩圖框應(yīng)重合），即在2009年這一期不能進(jìn)行差分干涉處理，只能選取其中一景與編程數(shù)據(jù)進(jìn)行處理。另外，3 m分辨率的數(shù)據(jù)不能全部覆蓋德欽縣城；8 m分辨率的數(shù)據(jù)雖然能夠覆蓋，但精度較低。另外，存檔數(shù)據(jù)為一景的話，往往不能與編程數(shù)據(jù)進(jìn)行差分干涉處理，原因是數(shù)據(jù)太少，時(shí)間過(guò)長(zhǎng)。綜上，加拿大的Radarsat-2衛(wèi)星不適用于德欽縣城的地質(zhì)災(zāi)害監(jiān)測(cè)。

圖1　加拿大的Radarsat-2衛(wèi)星存檔數(shù)據(jù)

意大利的Cosmo-Skymed衛(wèi)星。意大利的Cosmo-Skymed衛(wèi)星在德欽縣的歷史存檔數(shù)據(jù)共有11景，空間分辨率為3 m。2009年降軌數(shù)據(jù)為5景數(shù)據(jù)，其中有2組數(shù)據(jù)為同一參數(shù)的存檔數(shù)據(jù)，可用來(lái)InSAR處理；2009年升軌數(shù)據(jù)為6景數(shù)據(jù)，其中有1組數(shù)據(jù)為同一參數(shù)的存檔數(shù)據(jù)，可用來(lái)D-InSAR處理。

綜上，意大利衛(wèi)星存檔數(shù)據(jù)可形成3個(gè)干涉組:干涉組合1，升軌，日期分別為2009-09-03和2009-06-07，中心入射角51.12°；干涉組合2，降軌，日期分別為2009-09-14和2009-06-02，中心入射角: 48.08°；干涉組合3，降軌，日期分別為2009-05-31和2009-04-13，中心入射角29.36°。中心入射角過(guò)小，會(huì)降低數(shù)據(jù)信號(hào)的相干性，不利于InSAR分析，因此建議選擇干涉組合3進(jìn)行研究。

2.4 雷達(dá)衛(wèi)星數(shù)據(jù)評(píng)估與總結(jié)

在考慮德欽縣地震時(shí)間點(diǎn)的情況下，通過(guò)對(duì)比各衛(wèi)星，初步確定意大利Cosmo-Skymed衛(wèi)星較適應(yīng)于德欽縣城地質(zhì)災(zāi)害監(jiān)測(cè)項(xiàng)目。除了數(shù)據(jù)要求合格以外，數(shù)據(jù)價(jià)格也在預(yù)算之內(nèi)。

意大利的Cosmo-Skymed系統(tǒng)是由4顆雷達(dá)衛(wèi)星組成的衛(wèi)星星座，每顆衛(wèi)星擁有一個(gè)多模式高分辨率合成孔徑雷達(dá)。在德欽縣城地質(zhì)災(zāi)害監(jiān)測(cè)項(xiàng)目中建議采用的束波模式為條帶模式（STRIPMAP）中的Himage成像模式，分辨率3 m，幅寬分別為40×40 km2。采用的產(chǎn)品級(jí)別為L(zhǎng)evel1A產(chǎn)品，單視復(fù)斜距SLC（SCS），數(shù)據(jù)中包含了強(qiáng)度信息和相位信息，用于生成DEM或做差分干涉測(cè)量。

但該衛(wèi)星存在以下缺點(diǎn)。出售該衛(wèi)星數(shù)據(jù)的商家比較少，目前，只了解到1家，商家過(guò)少不利于對(duì)比選優(yōu)。該衛(wèi)星采用的是X波段，相對(duì)較短，對(duì)植被、云霧等條件的穿透能力較差，干涉誤差較大，選用時(shí)需避開(kāi)雨季進(jìn)行數(shù)據(jù)采集。該衛(wèi)星的歷史存檔數(shù)據(jù)為2009年，按照項(xiàng)目計(jì)劃，2015年底進(jìn)行數(shù)據(jù)編程，相隔6年之久可能會(huì)造成干涉效果不佳。

如果不考慮地震時(shí)間點(diǎn)，即選用2014年以后的數(shù)據(jù)，可以采用日本的ALOS衛(wèi)星，該衛(wèi)星費(fèi)用較高，超過(guò)預(yù)算要求，但其采用長(zhǎng)波段L，穿透能力較好，且該衛(wèi)星無(wú)須用戶要求，時(shí)時(shí)刻刻采集數(shù)據(jù)，較多的景數(shù)據(jù)易于進(jìn)行差分干涉處理，相對(duì)精度也更高，變形監(jiān)測(cè)效果更好。

3 雷達(dá)衛(wèi)星數(shù)據(jù)處理軟件比選

3.1 德欽縣雷達(dá)數(shù)據(jù)處理要求

D-InSAR是在InSAR的基礎(chǔ)上發(fā)展起來(lái)的一項(xiàng)技術(shù)，以InSAR復(fù)數(shù)據(jù)提供的相位信息（干涉圖）為信息源，采用多幅雷達(dá)影像的相位差，通過(guò)重復(fù)軌道干涉測(cè)量的原理消除地形信息從而得到地表形變。根據(jù)去除地形信息方法的不同，D-InSAR技術(shù)可以分為二軌法、三軌法和四軌法等，德欽縣城地質(zhì)災(zāi)害監(jiān)測(cè)項(xiàng)目采用常規(guī)的二軌法和三軌法進(jìn)行數(shù)據(jù)處理。

二軌法。首先利用同一地區(qū)前后2幅影像生成干涉紋圖，然后再利用外部獲取DEM數(shù)據(jù)模擬干涉紋圖，最后從兩干涉紋圖中去除地形信息得到地表變化信息。

三軌法。首先利用同一地區(qū)三景影像生成2幅干涉條紋圖，一幅反映地形信息，一幅反映地表形變信息，然后進(jìn)行去地平和相位解纏，最后利用差分干涉測(cè)量原理計(jì)算得到地表信息。

選擇合適有效的干涉影像數(shù)據(jù)集（SLC數(shù)據(jù)）后，從外部獲取DEM數(shù)據(jù)，以其中一幅影像作為主影像，另外一幅或幾幅影像作為從影像，便可以開(kāi)始干涉測(cè)量的數(shù)據(jù)差分處理。數(shù)據(jù)處理的過(guò)程包括2部分，第一部分是InSAR數(shù)據(jù)處理，主要包括:影像的配準(zhǔn)、影像重采樣、干涉圖的生成、參考相位（平地相位）的去除、地形相位的去除（D-InSAR處理用到）、干涉圖的濾波、相位解纏（三軌法）和地理編碼（坐標(biāo)系統(tǒng)的投影轉(zhuǎn)換）等，得到DEM和干涉圖；第2部分是差分處理，經(jīng)過(guò)差分干涉處理、相位解纏（二軌法）、相高轉(zhuǎn)化等步驟，生成形變圖。

3.2 常用雷達(dá)數(shù)據(jù)處理軟件

目前，InSAR數(shù)據(jù)處理軟件有許多種，有商業(yè)的、開(kāi)源的、試用版的，等等。商業(yè)的專門(mén)用來(lái)處理InSAR數(shù)據(jù)，有特定的模塊，處理方式比較簡(jiǎn)單，但費(fèi)用較高；開(kāi)源的軟件需要進(jìn)行程序編碼才能使用，需耗費(fèi)時(shí)間與精力；試用版的軟件僅適合個(gè)人學(xué)習(xí)和使用，不適合單位集體的研究。具體分述如下。

Gamma。強(qiáng)大的雷達(dá)成像、干涉處理、差分處理軟件，擁有DInSAR，PSInSAR及地理編碼等專業(yè)處理功能，屬于商業(yè)軟件。

SARscape。強(qiáng)大的雷達(dá)圖像數(shù)據(jù)處理與分析軟件，擁有極化干涉處理模塊，屬于商業(yè)軟件。但市面上亦存在盜版。

ERDAS。具有雷達(dá)圖像數(shù)據(jù)導(dǎo)入功能，有IFSAR，StereoSAR，OrthoRadar模塊，但不是專業(yè)的SAR處理軟件，屬于盜版軟件。

ROI-PACK。干涉SAR領(lǐng)域開(kāi)源中最為強(qiáng)大的軟件，屬于開(kāi)源軟件。

Doris。面向全球的開(kāi)源軟件，有人專門(mén)開(kāi)發(fā)了InSAR/DInSAR等功能，屬于開(kāi)源軟件。

3.3 數(shù)據(jù)處理軟件對(duì)比

常用的InSAR數(shù)據(jù)處理中，最常用、且最好用的軟件為Gamma和SARscape，下面對(duì)其各自特點(diǎn)進(jìn)行簡(jiǎn)單介紹。

Gamma。該軟件能夠從軌道數(shù)據(jù)中進(jìn)行基線估計(jì)、干涉圖像的配準(zhǔn)、干涉紋圖的生成、相干系數(shù)的估計(jì)、去平效應(yīng)、自適應(yīng)濾波等步驟。如果從外部獲得了DEM，可以使用二軌法，利用這個(gè)模塊模擬干涉紋圖，利用干涉紋圖進(jìn)行差分處理就可以消除地形相位的影響。此外，該軟件還具有永久散射體分析模塊，可針對(duì)特定目標(biāo)點(diǎn)的相位進(jìn)行分析，達(dá)到高精度誤差分析。

SARscape。ENVISARscape軟件是在ENVI遙感圖像處理軟件的基礎(chǔ)之上研發(fā)的，功能主要包括:雷達(dá)干涉測(cè)量、雷達(dá)強(qiáng)度圖像處理、極化雷達(dá)干涉測(cè)量。支持多種雷達(dá)產(chǎn)品數(shù)據(jù)和原始數(shù)據(jù)，提供SAR數(shù)據(jù)的數(shù)據(jù)導(dǎo)入、基線估計(jì)、多視處理、幾何校正、輻射校正、去噪影像、地表變形分析等基本功能；同時(shí)，支持極化SAR和極化干涉SAR數(shù)據(jù)的處理等功能，處理結(jié)果可以達(dá)到mm級(jí)的位移。

3.4 數(shù)據(jù)處理軟件評(píng)估與總結(jié)

Gamma，SARscape這2款軟件都可以進(jìn)行InSAR數(shù)據(jù)處理，相比起來(lái)，Gamma更加專業(yè)，功能更加強(qiáng)大，操作簡(jiǎn)單，但費(fèi)用較高，國(guó)內(nèi)擁有授權(quán)的公司很少。SARscape軟件正式版費(fèi)用較高，與Gamma軟件功能大體相同，但是該軟件有試用版，適合個(gè)人學(xué)習(xí)和研究。這2款軟件的主要區(qū)別是操作系統(tǒng)不同和處理效率不同: Gamma主要使用的操作系統(tǒng)為L(zhǎng)inux系統(tǒng)，處理方式為輸入命令進(jìn)行處理，可以進(jìn)行批量處理；SARscape采用的是Windows操作系統(tǒng)，數(shù)據(jù)處理采用可視化窗口模式，處理速度較慢，不能進(jìn)行批量處理。據(jù)此，建議綜合使用Gamma軟件和SARscape軟件來(lái)進(jìn)行數(shù)據(jù)處理，Gamma軟件負(fù)責(zé)批量處量，SARscape軟件負(fù)責(zé)個(gè)人學(xué)習(xí)和研究。

4 小結(jié)

德欽縣地質(zhì)災(zāi)害監(jiān)測(cè)工作若要順利開(kāi)展，勢(shì)必在設(shè)計(jì)階段對(duì)收集到的雷達(dá)衛(wèi)星資料進(jìn)行分析和評(píng)估，為監(jiān)測(cè)項(xiàng)目的順利實(shí)施提供理論參考和技術(shù)依據(jù)，為該項(xiàng)目的立項(xiàng)奠定可行性基礎(chǔ)。結(jié)合選取的雷達(dá)衛(wèi)星和處理軟件，采用D-InSAR技術(shù)能夠快速獲取德欽縣地表形變區(qū)域內(nèi)的地表形變信息，能夠滿足對(duì)地表形變引起的泥石流的監(jiān)測(cè)要求。通過(guò)比選，初步確定意大利Cosmo-Skymed衛(wèi)星為德欽縣地質(zhì)災(zāi)害監(jiān)測(cè)的初選衛(wèi)星數(shù)據(jù)，Gamma為初選軟件。

參考文獻(xiàn):

［1］ CARNECC，MASSONNETD，KINGC.Twoexamplesofthe useofSARinterferometryon displacementfieldsofsmallspatial extent［J］.GeophysicalResearch Letters，1996，23（24）: 3579-3582.

［2］ SINGHROY V，MOLCH K.Characterizingand monitoring rockslidesfrom SAR techniques［J］.Advancesin Space Research，2004，33（3）: 290-295.

［3］ 王桂杰，謝謨文，邱騁，等.D-INSAR技術(shù)在大范圍滑坡監(jiān)測(cè)中的應(yīng)用［J］.巖土力學(xué)，2010，31（4）: 1337-1344.

［4］ 梁軍.雷達(dá)干涉測(cè)量及其在青藏鐵路沿線地面形變監(jiān)測(cè)中的應(yīng)用研究［D］.成都:成都理工大學(xué)，2007.

［5］ 馮珂.InSAR復(fù)圖像配準(zhǔn)方法研究［D］.長(zhǎng)沙:長(zhǎng)沙理工大學(xué)，2013.

［6］ 匡薇.Radarsat-2數(shù)據(jù)處理方法研究及其在地面沉降中的應(yīng)用［D］.北京:中國(guó)地質(zhì)大學(xué)，2014.

［7］ 葛大慶，王艷，范景輝，等.地表形變D-InSAR監(jiān)測(cè)方法及關(guān)鍵問(wèn)題分析［J］.國(guó)土資源遙感，2007（4）: 14-22.

［8］ 張佳磊.基于相干目標(biāo)的D-InSAR方法研究及在北京地面沉降監(jiān)測(cè)中的應(yīng)用［D］.北京:中國(guó)地質(zhì)大學(xué)，2013.

（責(zé)任編輯：高 峻）

中圖分類號(hào)：P694

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：0528-9017（2016）06-0926-04

DOI：10.16178/j.issn.0528-9017.20160642

收稿日期:2016-03-04

作者簡(jiǎn)介:牛宗碩（1989—），男，河北張家口人，碩士研究生，從事衛(wèi)星遙感研究工作，E-mail: nizongshuo123@163.com。