基于InSAR相干系數(shù)的河湖岸線(xiàn)水事活動(dòng)監(jiān)測(cè)

陳子琪 陳亮 李建松

摘 要：及時(shí)發(fā)現(xiàn)河湖岸線(xiàn)水事活動(dòng)對(duì)監(jiān)管河湖岸線(xiàn)、維護(hù)河湖健康具有重要意義。選取黃河小北干流4期Radarsat-2遙感影像形成2期干涉像對(duì)，對(duì)干涉數(shù)據(jù)進(jìn)行數(shù)據(jù)預(yù)處理生成相干系數(shù)圖。利用地物的干涉相干特性，采用差值閾值法獲取河湖岸線(xiàn)水事活動(dòng)疑似“熱點(diǎn)”圖斑，以同期高分辨率可見(jiàn)光遙感影像為驗(yàn)證數(shù)據(jù)對(duì)結(jié)果進(jìn)行驗(yàn)證和分析。結(jié)果表明，基于相干系數(shù)變化可以成功發(fā)現(xiàn)河道臨河房屋建設(shè)、堆放、占用等水事活動(dòng)，利用該方法進(jìn)行河道水事活動(dòng)監(jiān)測(cè)是可行的。

關(guān)鍵詞：相干系數(shù);水事活動(dòng);InSAR;發(fā)現(xiàn)監(jiān)測(cè);河湖岸線(xiàn)

中圖分類(lèi)號(hào)：TP79 ? 文獻(xiàn)標(biāo)志碼：A

doi：10.3969/j.issn.1000-1379.2020.07.035

Abstract：Water matter activity monitoring is of great significance for river and lake shoreline supervision and maintenance of river and lake health. Four radarsat-2 remote sensing images of Xiaobeiganliu reach of the Yellow River were chosen to form phase 2 interference image pairs. Coherence coefficient difference maps were generated through data preprocessing of the interference images. Based on the coherence characteristics of the ground object in InSAR images， the suspected “hot spot” pattern of water matter activity in the river and lake shoreline were obtained by the threshold method. The monitoring results were verified by the high-resolution visible remote sensing image of the same period. The results show that the methods based on the change of coherence coefficient is feasible and effective to find water matter activity of river and lake shoreline such as house construction， stacking and occupation.

Key words： coherent coefficient; water matter activity; InSAR; monitoring; river and lake shoreline

1 引 言

河湖水系是水資源的重要載體，在防洪、供水、灌溉、發(fā)電、航運(yùn)、生態(tài)等方面具有不可替代的作用。長(zhǎng)期以來(lái)，河湖岸線(xiàn)是新老水問(wèn)題最為集中的區(qū)域，違規(guī)建設(shè)開(kāi)發(fā)、圍墾河湖、侵占河道、違法養(yǎng)殖、非法采砂等水事活動(dòng)頻發(fā)，河道和湖泊面積萎縮、水域空間減小、水質(zhì)惡化、生物棲息地破壞等問(wèn)題十分突出，嚴(yán)重影響了河湖健康。2016年和2017年中共中央辦公廳、國(guó)務(wù)院辦公廳先后印發(fā)了《關(guān)于全面推行河長(zhǎng)制的意見(jiàn)》和《關(guān)于在湖泊實(shí)施湖長(zhǎng)制的指導(dǎo)意見(jiàn)》，在全國(guó)江河湖泊全面推行河湖長(zhǎng)制，為維護(hù)河湖健康生命、實(shí)現(xiàn)河湖功能永續(xù)利用提供了制度保障。在全國(guó)范圍內(nèi)開(kāi)展的河湖清“四亂”專(zhuān)項(xiàng)行動(dòng)是推動(dòng)河長(zhǎng)制從“有名”向“有實(shí)”轉(zhuǎn)變的第一抓手，是水利行業(yè)強(qiáng)監(jiān)管的重要內(nèi)容，是對(duì)各級(jí)河長(zhǎng)、湖長(zhǎng)和水行政主管部門(mén)履職盡責(zé)的底線(xiàn)要求。

2019年9月，習(xí)近平總書(shū)記在黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展座談會(huì)講話(huà)中指出，當(dāng)前黃河流域仍存在一些突出困難和問(wèn)題，流域生態(tài)環(huán)境脆弱，水資源保障形勢(shì)嚴(yán)峻，發(fā)展質(zhì)量有待提高。要堅(jiān)持生態(tài)優(yōu)先、綠色發(fā)展，以水而定、量水而行，因地制宜、分類(lèi)施策，上下游、干支流、左右岸統(tǒng)籌謀劃，共同抓好大保護(hù)，協(xié)同推進(jìn)大治理。利用科技手段加強(qiáng)黃河干流及其支流河道監(jiān)管，快速發(fā)現(xiàn)河道“四亂”變化情況，第一時(shí)間制止河湖“四亂”問(wèn)題，防止事態(tài)擴(kuò)大、減小損失，對(duì)保護(hù)黃河流域生態(tài)環(huán)境、維護(hù)河湖健康生命具有重大意義。

利用遙感技術(shù)進(jìn)行城市、礦區(qū)等土地利用開(kāi)發(fā)活動(dòng)監(jiān)測(cè)已經(jīng)取得廣泛應(yīng)用。張翰超等[1]利用高分辨率可見(jiàn)光遙感影像對(duì)2000—2015年中國(guó)省會(huì)城市擴(kuò)張情況進(jìn)行高精度監(jiān)測(cè);翟浩然等[2]利用Landsat遙感影像對(duì)1990—2015年中國(guó)超大城市建設(shè)用地時(shí)空變化進(jìn)行監(jiān)測(cè);占昕等[3]結(jié)合Landsat遙感影像和DMSP/OLS夜間燈光數(shù)據(jù)進(jìn)行廈門(mén)市城市擴(kuò)張測(cè)度研究。以上成功應(yīng)用為城市規(guī)劃、發(fā)展決策提供支撐。趙鵬[4]、王海慶等[5]利用高分辨率可見(jiàn)光遙感影像進(jìn)行礦區(qū)開(kāi)發(fā)監(jiān)測(cè)，為礦區(qū)監(jiān)管提供依據(jù);李春雷等[6]利用高分辨率可見(jiàn)光遙感影像進(jìn)行河湖清“四亂”應(yīng)用的探索研究。河道水事活動(dòng)發(fā)展速度快，對(duì)監(jiān)測(cè)時(shí)效性要求高，可見(jiàn)光遙感技術(shù)受云、霧等天氣影響無(wú)法開(kāi)展有效監(jiān)測(cè)。合成孔徑雷達(dá)（SAR）具有全天時(shí)、全天候的監(jiān)測(cè)優(yōu)勢(shì)，可以不受天氣因素影響及時(shí)發(fā)現(xiàn)地表變化。干涉雷達(dá)（InSAR）可以獲取地表相干特性，在地表變化監(jiān)測(cè)方面具有獨(dú)特優(yōu)勢(shì)。盧婷婷[7]利用InSAR技術(shù)對(duì)城市和礦區(qū)地表變化進(jìn)行監(jiān)測(cè);Drezet等[8]利用ERS干涉相干系數(shù)圖像進(jìn)行森林覆蓋率提取;張舞燕等[9]利用ENVISAT ASAR干涉相干系數(shù)進(jìn)行豫北地區(qū)城市邊界提取，取得較好效果;郭唯娜等[10]采用兩個(gè)時(shí)期InSAR干涉相干系數(shù)取閾值的方法獲取東帕米爾高原冰川邊界，提取精度達(dá)到90%。地表變化引起相干性發(fā)生變化，為采用InSAR相干系數(shù)進(jìn)行河道水事活動(dòng)監(jiān)測(cè)提供了新監(jiān)測(cè)手段。

采用不同時(shí)期高分辨率SAR干涉像對(duì)生成相干系數(shù)圖，基于水事活動(dòng)發(fā)生后河道地表在雷達(dá)干涉影像上相干性變化特性，采用差值閾值法進(jìn)行水事活動(dòng)“熱點(diǎn)”圖斑監(jiān)測(cè)，及時(shí)發(fā)現(xiàn)水事活動(dòng)疑似點(diǎn)，為河湖監(jiān)管提供有力支撐。

2 研究區(qū)概況與數(shù)據(jù)

2.1 研究區(qū)概況

選取黃河小北干流河段為研究區(qū)，該河段位于黃河干流禹門(mén)口和潼關(guān)之間，河段長(zhǎng)132.5 km，為晉陜兩省界河。左岸為山西省運(yùn)城市所屬河津、萬(wàn)榮、臨猗、永濟(jì)、芮城5縣（市），右岸為陜西省渭南市所屬韓城、合陽(yáng)、大荔、潼關(guān)4縣（市）。河道穿行于汾、渭地塹谷凹地區(qū)，沿程有汾河、涺水、涑水河、渭河、北洛河等支流匯入，兩岸為高出河床10～200 m的黃土臺(tái)塬。黃河小北干流河道寬4～19 km，河道總面積1 107 km2，其中灘區(qū)面積657 km2，屬淤積性游蕩型河道。寬廣的河道灘區(qū)居住了數(shù)十萬(wàn)居民，為了生產(chǎn)生活需要在河道內(nèi)開(kāi)展大量生產(chǎn)建設(shè)活動(dòng)，河道內(nèi)違規(guī)建設(shè)開(kāi)發(fā)、侵占河道、非法采砂等水事活動(dòng)時(shí)有發(fā)生，是河道監(jiān)管的重點(diǎn)關(guān)注河段。

2.2 研究數(shù)據(jù)

以Radarsat-2雷達(dá)衛(wèi)星遙感影像為研究數(shù)據(jù)進(jìn)行水事活動(dòng)監(jiān)測(cè)研究。Radarsat-2雷達(dá)衛(wèi)星由加拿大太空署與MDA公司合作研發(fā)，2007年12月14日發(fā)射升空。衛(wèi)星軌高798 km，重訪(fǎng)周期為24 d，成像波段為C波段，具有HH、VV、HV、VH 4種極化方式和11種波束模式，成像空間分辨率為1～100 m，單景幅寬為20～500 km。選用3 m空間分辨率、超精細(xì)波束模式、VV極化方式的Radarsat-2雷達(dá)衛(wèi)星遙感影像，共采集了4期影像組成兩對(duì)干涉像對(duì)，分別為2016年12月28日和2017年1月21日的干涉像對(duì)、2017年8月1日和2017年8月25日的干涉像對(duì)，Radarsat-2 SAR影像相關(guān)參數(shù)見(jiàn)表1。

選用SRTM 90 m空間分辨率DEM數(shù)據(jù)作為參考DEM進(jìn)行干涉雷達(dá)影像預(yù)處理。選用2 m空間分辨率的可見(jiàn)光遙感影像進(jìn)行結(jié)果驗(yàn)證與對(duì)比分析。

影響相干系數(shù)大小的因素主要有系統(tǒng)熱噪聲、空間基線(xiàn)、時(shí)間基線(xiàn)、方位向多普勒中心偏移、雷達(dá)波穿透引起的體散射、復(fù)數(shù)影像配準(zhǔn)誤差等[10，12]。本研究暫不考慮系統(tǒng)熱噪聲和雷達(dá)波穿透引起的體散射影響。為減小其他因素對(duì)相干系數(shù)的影響，本研究采用以下處理辦法：選取時(shí)間間隔短、空間基線(xiàn)小的干涉像對(duì)以降低時(shí)間、空間的相干影響;數(shù)據(jù)處理時(shí)對(duì)干涉像對(duì)兩期影像進(jìn)行精配準(zhǔn)和濾波，消除配準(zhǔn)誤差和方位向多普勒中心偏移引起的相干影響。

不同地物具有不同的干涉相干特性，居民地等人工地物相干系數(shù)一般高于植被、水體等，且隨時(shí)間變化不大，根據(jù)這一特性可通過(guò)干涉數(shù)據(jù)進(jìn)行地物分類(lèi)[7，13-14]。在河道內(nèi)開(kāi)展占、采、堆、建等水事活動(dòng)，會(huì)引起地表地物變化，如植被、裸灘變成建筑物，水面、裸灘變成砂堆、棄土棄渣等。水事活動(dòng)在不同時(shí)期相干系數(shù)圖上表現(xiàn)為相同位置相干系數(shù)發(fā)生顯著變化，以此可以利用多期相干系數(shù)差異進(jìn)行水事活動(dòng)監(jiān)測(cè)。

3.2 研究方法

選用4期Radarsat-2遙感影像形成2對(duì)干涉像對(duì)，利用SARscape軟件進(jìn)行干涉數(shù)據(jù)處理，在對(duì)干涉像對(duì)分別進(jìn)行數(shù)據(jù)預(yù)處理的基礎(chǔ)上生成相干系數(shù)圖，采用差值閾值法進(jìn)行水事活動(dòng)信息提取，并采用高分辨率可見(jiàn)光遙感影像對(duì)監(jiān)測(cè)結(jié)果進(jìn)行驗(yàn)證和分析，技術(shù)路線(xiàn)見(jiàn)圖1。

3.2.1 數(shù)據(jù)預(yù)處理

由于干涉像對(duì)兩期影像成像時(shí)傳感器的位置不一樣，影像同名點(diǎn)之間存在位置偏移，因此采用重復(fù)軌道數(shù)據(jù)進(jìn)行干涉測(cè)量時(shí)需要對(duì)影像進(jìn)行精配準(zhǔn)，將配準(zhǔn)后的影像生成干涉圖，并去除干涉條紋。InSAR的相干相位包括平地分量和高差分量，通過(guò)干涉去平、Goldstein自適應(yīng)濾波等數(shù)據(jù)處理方法去除平地相位差、保留高程相位差，在此基礎(chǔ)上進(jìn)行相干性計(jì)算，并進(jìn)行地理編碼，生成相干系數(shù)圖。

3.2.2 信息提取

不同時(shí)期相干系數(shù)圖之間存在幾何位置偏差，在進(jìn)行信息提取前需要對(duì)兩期相干系數(shù)圖進(jìn)行影像配準(zhǔn)，使同名點(diǎn)在兩個(gè)時(shí)期圖像上具有相同的坐標(biāo)位置。閾值法是遙感影像變化監(jiān)測(cè)的常用方法[10，12]，本文采用差值閾值法進(jìn)行水事活動(dòng)監(jiān)測(cè)。將后一期與前一期相干系數(shù)圖進(jìn)行差值計(jì)算得到相干系數(shù)差值圖，將差值大于一定閾值的像元判定為水事活動(dòng)疑似“熱點(diǎn)”，根據(jù)差值圖和野外實(shí)地查勘資料將閾值設(shè)定為0.15。采用Majority鄰域后處理方法對(duì)初步結(jié)果圖像進(jìn)行后處理，窗口大小為3×3。最后去除噪聲引起的過(guò)小圖斑，并將柵格成果轉(zhuǎn)成矢量，得到水事活動(dòng)疑似“熱點(diǎn)”圖斑。

4 結(jié)果與分析

根據(jù)前文提出的監(jiān)測(cè)方法進(jìn)行水事活動(dòng)監(jiān)測(cè)，獲得水事活動(dòng)疑似“熱點(diǎn)”圖斑。為了方便結(jié)果驗(yàn)證和分析，選取龍門(mén)鎮(zhèn)部分監(jiān)測(cè)結(jié)果，相干系數(shù)差值和監(jiān)測(cè)結(jié)果如圖2所示。

選取2016年6月27日和2017年6月20日的2 m高分辨率衛(wèi)星遙感影像對(duì)該區(qū)域監(jiān)測(cè)結(jié)果進(jìn)行驗(yàn)證分析，見(jiàn)圖3。A區(qū)域?yàn)橐惶幮略龅呐R河房屋，圖2（a）中該區(qū)域相干系數(shù)明顯增大，與2016年影像相比，2017年區(qū)域內(nèi)部分土地進(jìn)行了平整，部分房屋已經(jīng)建成;B區(qū)域?yàn)橐惶帡壴逊艌?chǎng)，2016年該棄渣堆為黑色棄渣且堆放形態(tài)較為規(guī)整，2017年棄渣顏色變?yōu)辄S色，堆放形態(tài)也發(fā)生變化，為一個(gè)個(gè)散亂分布的小堆，導(dǎo)致相干系數(shù)發(fā)生變化;C區(qū)域?yàn)橐惶幟娣e擴(kuò)大的棄渣堆放場(chǎng)，與2016年相比，2017年棄渣堆形態(tài)發(fā)生較大變化，棄渣堆向南北兩側(cè)都有擴(kuò)張。由此可見(jiàn)，基于InSAR相干系數(shù)可以發(fā)現(xiàn)河道開(kāi)發(fā)建設(shè)、堆放、占用等水事活動(dòng)變化。但是，由于雷達(dá)影像存在噪聲、建筑物存在多次散射且散射情況復(fù)雜，因此結(jié)果中會(huì)存在一些偽變化點(diǎn)，需要進(jìn)一步判別處理。

5 結(jié) 語(yǔ)

采用4期Radarsat-2遙感影像形成2對(duì)干涉像對(duì)，提出了基于InSAR相干系數(shù)的河道水事活動(dòng)疑似“熱點(diǎn)”圖斑的監(jiān)測(cè)方法，在黃河小北干流河段進(jìn)行試驗(yàn)，成功發(fā)現(xiàn)河道臨河房屋建設(shè)、堆放、占用等水事活動(dòng)，表明利用相干系數(shù)進(jìn)行河道水事活動(dòng)發(fā)現(xiàn)監(jiān)測(cè)是可行的。利用該方法進(jìn)行河道水事活動(dòng)監(jiān)測(cè)能夠快速發(fā)現(xiàn)疑似對(duì)象，在此基礎(chǔ)上進(jìn)行水事活動(dòng)解譯，可以縮小解譯范圍，大大降低解譯工作量，為河湖岸線(xiàn)監(jiān)管提供重要支撐。

本研究的閾值選取主要依據(jù)人工經(jīng)驗(yàn)，不同遙感影像、不同作業(yè)區(qū)域需要重新確定閾值。雷達(dá)影像存在噪聲、建筑物存在多次散射造成監(jiān)測(cè)結(jié)果存在一些誤差。將來(lái)需要在閾值自動(dòng)確定和利用多源數(shù)據(jù)抑制雷達(dá)噪聲等方面開(kāi)展研究，進(jìn)一步提高監(jiān)測(cè)的自動(dòng)化程度和準(zhǔn)確度。

參考文獻(xiàn)：

[1] 張翰超，寧曉剛，王浩，等.基于高分辨率遙感影像的2000—2015年中國(guó)省會(huì)城市高精度擴(kuò)張監(jiān)測(cè)與分析[J].地理學(xué)報(bào)，2018，73（12）：2345-2363.

[2] 翟浩然，唐新明，王光輝，等.中國(guó)超大城市建設(shè)用地遙感監(jiān)測(cè)與時(shí)空格局[J].遙感信息，2019，34（4）：37-47.

[3] 占昕，潘文斌，蔡芫鑌，等.基于多源遙感數(shù)據(jù)的城市擴(kuò)張測(cè)度研究：以廈門(mén)市為例[J].福州大學(xué)學(xué)報(bào)（自然科學(xué)版），2017，45（3）：355-361.

[4] 趙鵬.新街臺(tái)格廟礦區(qū)開(kāi)發(fā)初期土地利用動(dòng)態(tài)遙感監(jiān)測(cè)[J].國(guó)土資源遙感，2015，27（4）：144-149.

[5] 王海慶，武明德，劉瓊，等.山西某地礦業(yè)開(kāi)發(fā)占用土地狀況遙感監(jiān)測(cè)[J].國(guó)土資源遙感，2018，30（1）：116-120.

[6] 李春雷，劉立聰，張雅莉，等.遙感技術(shù)在河湖“清四亂”中的應(yīng)用方法探究[J].浙江水利科技，2019（4）：74-77.

[7] 盧婷婷.InSAR技術(shù)在地表變化監(jiān)測(cè)中的應(yīng)用[D].西安：長(zhǎng)安大學(xué)，2015：63.

[8] DREZET P M L， QUEGAN S. Satellite-Based Radar Mapping of British Forest Age and Net Ecosystem Exchange Using ERS Tandem Coherence[J].Forest Ecology and Management，2007，238（1-3）：65-80.

[9] 張舞燕，劉清臣，孟秀軍.基于SAR相干系數(shù)圖像的城市邊界提取[J].測(cè)繪與空間地理信息，2014，37（5）：56-59.

[10] 郭唯娜，柯長(zhǎng)青，范宇賓.基于SAR干涉數(shù)據(jù)的東帕米爾高原冰川變化[J].地球信息科學(xué)學(xué)報(bào)，2019，21（11）：1790-1801.

[11] LI Z F， BAO Z， LI H， et al. Image Autocoregistration and InSAR Interferogram Estimation Using Joint Subspace Projection[J].IEEE Transaction on Geoscience and Remote Sensing，2006，44（2）：288-297.

[12] 何曙光，龐蕾，張學(xué)東，等.利用差值法與InSAR相干系數(shù)檢測(cè)城市變化信息[J].遙感信息，2017，32（3）：76-80.

[13] 李新武，郭華東，李震，等.干涉雷達(dá)地表相干特性分析及土地類(lèi)型分類(lèi)[J].高技術(shù)通訊，2002（5）：1-4.

[14] 肖偉山.雙頻雷達(dá)后向散射和相干性分析及土地利用覆蓋分類(lèi)應(yīng)用[D].福州：福州大學(xué)，2010：46.

【責(zé)任編輯 張 帥】