陳慧欣， 陳超， 張自力， 汪李彥， 梁錦濤

(1.浙江海洋大學(xué)海洋科學(xué)與技術(shù)學(xué)院，舟山 316022； 2.浙江省生態(tài)環(huán)境監(jiān)測中心(浙江省生態(tài)環(huán)境監(jiān)測預(yù)警及質(zhì)控研究重點(diǎn)實(shí)驗(yàn)室),杭州 310012)

0 引言

潮間帶是陸地生態(tài)系統(tǒng)和海洋生態(tài)系統(tǒng)之間的緩沖區(qū)域，是濱海濕地的重要組成部分，具有較高的生態(tài)經(jīng)濟(jì)價(jià)值，通常位于高潮線與低潮線之間。在陸地環(huán)境和海水的雙重作用下，潮間帶濕地有著極其豐富的生物資源[1]，在維持海岸線穩(wěn)定、促進(jìn)碳匯、調(diào)節(jié)氣候和維持生態(tài)平衡等方面發(fā)揮極高的生態(tài)價(jià)值[2-5]。同時(shí)也在水產(chǎn)養(yǎng)殖、供應(yīng)土地資源和旅游資源等方面創(chuàng)造高經(jīng)濟(jì)效益[6-7]。然而，受全球氣候變化、海平面上升以及人類活動的影響[8]，潮間帶面積萎縮，生物多樣性下降，濕地環(huán)境面臨嚴(yán)峻威脅。及時(shí)準(zhǔn)確地掌握潮間帶濕地的分布狀況和演變規(guī)律，有利于實(shí)現(xiàn)潮間帶濕地的科學(xué)管理，促進(jìn)海岸帶地區(qū)的可持續(xù)發(fā)展[9]。

遙感技術(shù)的發(fā)展為獲取潮間帶信息提供方便快捷的手段，更多高空間分辨率的衛(wèi)星影像在監(jiān)測潮間帶變化過程中提供數(shù)據(jù)支撐。Google Earth Engine(GEE)是一個(gè)行星尺度的地理空間信息分析平臺，為地球觀測研究提供了巨大的潛力[10-11]。GEE平臺可以訪問和計(jì)算Landsat/Sentinel 等公開獲取的遙感影像數(shù)據(jù)和其他輔助數(shù)據(jù)，使得快速處理和分析海量的地理空間數(shù)據(jù)成為可能，為進(jìn)行大尺度長時(shí)間序列的分析和應(yīng)用提供技術(shù)支撐[12]。目前在土地利用變化[13-14]、濱海濕地分類[15]、監(jiān)測森林類型[16-17]、水體信息的提取和動態(tài)變化[18-19]、紅樹林制圖[20]以及監(jiān)測潮間帶變化[21-22]等領(lǐng)域都得到廣泛應(yīng)用。

潮間帶是一個(gè)復(fù)雜的地理體，遙感監(jiān)測潮間帶要將衛(wèi)星影像與地學(xué)知識相結(jié)合，提取遙感數(shù)據(jù)中所包含的光譜空間信息等，獲取潮間帶地形的動態(tài)變化信息。國內(nèi)外學(xué)者利用遙感技術(shù)對潮間帶調(diào)查和制圖開展了大量研究。從小尺度如監(jiān)測葡萄牙里斯本的塔古斯河口[23]、廣西防城港市[24]，到中大尺度如澳大利亞潮間帶[25]、中國沿海地區(qū)[26-27]、東亞濱海地區(qū)[28]，以及對全球尺度的潮間帶地區(qū)進(jìn)行制圖和監(jiān)測[29]。在研究方法方面，部分學(xué)者利用潮汐高程模型和基于地形的方法來劃定潮間帶的空間范圍和灘涂[30-32]，隨機(jī)森林算法[33]和支持向量機(jī)算法[34]等基于訓(xùn)練樣本的機(jī)器學(xué)習(xí)方法也得到廣泛應(yīng)用。此外，還有通過選取潮灘的普遍閾值來進(jìn)行灘涂制圖的基于知識的決策樹方法[21]?？偟膩碚f，提取大陸海岸潮間帶信息的方法依賴潮位數(shù)據(jù)、訓(xùn)練樣本或者人工干預(yù)閾值等輔助數(shù)據(jù)的使用，且實(shí)現(xiàn)過程較為復(fù)雜。同時(shí)面向群島區(qū)域的研究較少，缺乏對海島潮灘資源較為深入和系統(tǒng)的研究，需要一種應(yīng)用于群島區(qū)域潮間帶信息的自動、快速提取方法。

舟山群島四面環(huán)海，云雨天氣較多，潮位處于不斷變化之中，海島地區(qū)陸地面積狹小，潮灘零星分布，難以在最低潮和最高潮時(shí)捕捉到高質(zhì)量的遙感影像，進(jìn)而無法準(zhǔn)確提取潮間帶信息。針對此問題，本研究選擇舟山群島作為研究區(qū)，基于GEE云平臺，以Landsat8時(shí)間序列影像為主要數(shù)據(jù)源，結(jié)合舟山群島海岸帶地理特征選擇適合的遙感指數(shù)，利用最大光譜指數(shù)合成算法(maximum spectral index composite，MSIC)和最大類間方差法(Otsu algorithm，OTSU)構(gòu)建多層自動決策樹分類模型，并在此基礎(chǔ)上提出基于融合數(shù)字高程模型(digital elevation model，DEM)的決策樹算法。將提取結(jié)果與谷歌地球高空間分辨率影像進(jìn)行精度評價(jià)，驗(yàn)證本研究方法的可行性及穩(wěn)定性，為自動、快速實(shí)現(xiàn)潮間帶信息提取提供可行的方法和思路，為海岸帶資源的開發(fā)與保護(hù)提供技術(shù)支持。

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

1.1 研究區(qū)域概況

舟山群島隸屬于浙江省，地處中國東部黃金海岸線與長江黃金水道的交匯處，背靠長江三角洲廣闊經(jīng)濟(jì)腹地，研究區(qū)位置示意圖如圖1所示。舟山群島擁有2 085個(gè)島嶼和270多km深水海岸線，是中國第一大群島和重要港口城市，下轄定海、普陀2區(qū)和岱山、嵊泗2縣，2021年常住人口為116.5萬人。舟山群島四面環(huán)海，屬亞熱帶季風(fēng)氣候，冬暖夏涼，溫和濕潤，年平均氣溫為16℃。該研究區(qū)全年降水量顯著，常年降水量在927～1 620 mm之間，夏秋之際易受熱帶風(fēng)暴侵襲。舟山群島屬于低山丘陵地貌類型，主要由平原和丘陵組成。舟山群島植被屬中亞熱帶常綠闊葉林，以闊葉喬木和灌叢混交林為主[35]。

圖1 研究區(qū)位置示意圖Fig.1 Location of study area

1.2 數(shù)據(jù)獲取與處理

Landsat是美國國家航空航天局(national aero-nautics and space administration，NASA)發(fā)射的陸地衛(wèi)星系列，從1972年以來已經(jīng)發(fā)射了9顆，積累了大范圍、長時(shí)間序列的地表遙感影像數(shù)據(jù)，被廣泛應(yīng)用于對地觀測研究中。Landsat8衛(wèi)星于2013年2月發(fā)射，包括陸地成像儀(operational land imager，OLI)和熱紅外傳感器(thermal infrared sensor，TIRS)2種載荷，重訪周期為16 d[36-37]。根據(jù)目前數(shù)據(jù)源狀況和研究需求，以2021年的衛(wèi)星影像合成為一期遙感數(shù)據(jù)進(jìn)行分析，篩選2021年間云量小于25%、成像效果良好的Landsat8 OLI衛(wèi)星影像，利用質(zhì)量評估波段(quality assessment band，QA)去除云和陰影等觀測質(zhì)量較差的像素。

SRTM(shuttle Radar topography mission)DEM是由美國NASA和國防部國家測繪局聯(lián)合測量得到，獲取N60°～S56°之間的雷達(dá)影像數(shù)據(jù)，覆蓋地球80%以上的陸地表面。按精度可分為SRTM1和SRTM3這2種數(shù)據(jù)，分別對應(yīng)30 m和90 m空間分辨率精度的數(shù)據(jù)[38-39]。本文使用的DEM數(shù)據(jù)是SRTM 30 m空間分辨率數(shù)據(jù)，對舟山群島的高程和坡度信息進(jìn)行提取分析。

2 研究方法

研究利用2021年Landsat8影像集合構(gòu)建自動、快速提取潮間帶的方法流程，基于遙感指數(shù)特征提取水體信息，包括以下4個(gè)步驟： ①基于GEE云平臺，通過云量篩選和QA波段去云構(gòu)建Landsat8影像堆棧； ②深入分析舟山群島潮間帶濕地的時(shí)間序列變化特征，選取適合的光譜指數(shù)，基于 MSIC算法和統(tǒng)計(jì)分析，利用OTSU算法構(gòu)建多層決策樹自動分類模型，從Landsat8數(shù)據(jù)集合中提取潮間帶初始信息； ③融合DEM數(shù)據(jù)，去除養(yǎng)殖塘和湖泊等內(nèi)陸水體的干擾，提取潮間帶最終信息； ④利用谷歌地球的高空間分辨率影像對提取的潮間帶范圍進(jìn)行驗(yàn)證。具體技術(shù)路線如圖2所示。

圖2 潮間帶信息提取流程Fig.2 Workflow of intertidal information extraction

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

對Landsat8衛(wèi)星遙感數(shù)據(jù)進(jìn)行預(yù)處理，去除研究區(qū)內(nèi)2021年度Landsat8影像集合中的云像素，流程如下： 首先，在 GEE 中調(diào)用2021年的所有覆蓋研究區(qū)域的Landsat8(USGS Landsat8 Surface Reflectance Tier 1)產(chǎn)品，建立研究區(qū)初始密集時(shí)間序列影像集合，該產(chǎn)品包含已經(jīng)過大氣校正和輻射校正的地表反射率數(shù)據(jù)； 然后，篩選云量小于25%的影像，再利用 QA波段去除云和陰影等觀測質(zhì)量較差的像素，以獲得高質(zhì)量觀測數(shù)據(jù)。

2.2 潮間帶粗提取

MSIC算法是指基于時(shí)間序列Landsat8影像集合中選取一個(gè)光譜指數(shù)作為逐像素排序函數(shù)，根據(jù)集合中具有所選光譜指數(shù)的最大值來設(shè)置合成影像中的每個(gè)像素。在GEE平臺的API中提供了“image Collection. quality Mosaic()”的函數(shù)，可以使用該函數(shù)調(diào)用MSIC算法。

在本研究中，將高潮與低潮之間的淤泥質(zhì)或沙質(zhì)無植被海岸帶定義為潮間帶。首先，基于像素合成影像的MSIC算法，利用所選光譜指數(shù)最大值的像素來合成高質(zhì)量影像，在一定程度上可以避免云噪聲像素對影像質(zhì)量的影響，解決濱海地區(qū)的影像缺失問題，并且可以緩解因潮汐運(yùn)動帶來的潮位變化的不確定性； 其次，通過OTSU算法對改進(jìn)的歸一化差異水體指數(shù)(modified normalized difference water index，mNDWI)最大值合成圖像進(jìn)行二值分割，自動生成二元分類結(jié)果，即水體和非水體，獲得2021年度最大地表水覆蓋； 應(yīng)用直方圖統(tǒng)計(jì)對自動水體提取指數(shù)(automated water extraction index，AWEI)最小值合成圖像進(jìn)行閾值分割，有效去除非水體像元； 最后，應(yīng)用歸一化差異植被指數(shù)(normalized difference vegetation index，NDVI)的二值分割圖像，進(jìn)一步識別錯分地物并去除，完成初始潮間帶信息提取。

2.2.1 基于遙感指數(shù)的最大和最小水面合成

近年來，水體指數(shù)法因使用簡單和計(jì)算耗時(shí)少，已成為常用的水體信息提取方法。目前常用的水體指數(shù)主要有以下幾種[40]： McFeeters[41]提出基于綠光波段和近紅外波段的歸一化差異水體指數(shù)(normalized difference water index，NDWI)； Rogers等[42]提出基于紅光波段和短波紅外波段的NDWI； Xu[43]提出的基于綠光波段和中紅外波段的mNDWI； Feyisa等[44]提出的AWEI。

根據(jù)Ji等[45]和Jia等[30]的研究結(jié)果，在提取水體信息時(shí)，mNDWI比NDWI取得更好的效果，且具有穩(wěn)定的閾值。所以在提取潮間帶方面，mNDWI增強(qiáng)了堤壩、道路和池塘外緣等人工岸線與水體的差異，有效地抑制養(yǎng)殖池等其他類型背景水體的影響，因此本研究選擇 mNDWI 用于提取潮間帶水體信息。另外，AWEI是通過波段間的相減、相加以及給波段賦予不同的系數(shù)來實(shí)現(xiàn)水體像元和非水體像元最大程度的分離，通過設(shè)定合適的閾值，可以有效去除非水體像元。NDVI能夠有效識別植被信息[46]，彌補(bǔ) mNDWI 和 AWEI 水體指數(shù)的不足，成為本文中潮間帶信息提取的輔助指數(shù)。其公式分別為：

(1)

AWEInsh=4(GREEN-SWIR1)-(0.25NIR+2.75SWIR1)，

(2)

(3)

式中RED，GREEN，NIR和SWIR1分別為紅光波段、綠光波段、近紅外波段和短波紅外1波段反射率值。

接下來的幾天里，我把目光牢牢地鎖定在了她們身上。我發(fā)現(xiàn)，但凡是班級排隊(duì)，諸如上下午大課間、體育課、信息課以及放學(xué)，她們都會刻意保持距離，儼然熟悉的陌生人。我找了個(gè)由頭，調(diào)整了排隊(duì)站位順序，當(dāng)然主要還是為了將她們隔開，試圖讓距離淡化她們之間的矛盾。

2.2.2 基于OTSU算法的水陸分離

OTSU算法是大津展之于1979年提出的基于圖像的統(tǒng)計(jì)屬性選取最佳閾值的非參數(shù)方法，適用于決策樹分類方法。按照圖像的灰度特性，基于最小二乘法原理將圖像分割成背景和前景2類，計(jì)算其類間方差最大、類內(nèi)方差最小的值為最優(yōu)分割閾值，進(jìn)行圖像的自動二值分類得到2類像素最大的區(qū)分度。原理如下： 將圖像l(x,y)分為前景和背景，分割閾值記作T，采用遍歷的方法得到使類間方差g最大的閾值T，即為分割閾值[47-48]。公式為：

g=ω0ω1(μ0-μ1)2，

(4)

式中：μ為圖像的總平均灰度；g為類間方差；ω0為前景圖像像素?cái)?shù)占整幅圖像的比例；μ0為前景圖像平均灰度；ω1為背景圖像像素?cái)?shù)占整幅圖像的比例；μ1為背景圖像平均灰度。

2.3 潮間帶精提取

DEM廣泛應(yīng)用于地學(xué)的各個(gè)領(lǐng)域，如地質(zhì)地貌學(xué)、地圖制圖和水文分析等各項(xiàng)研究，為全球的地理信息提取和分析提供數(shù)據(jù)支持。海水養(yǎng)殖是舟山市漁業(yè)的特色，養(yǎng)殖空間向淺海擴(kuò)展，全市養(yǎng)殖塘數(shù)量較多，且和潮灘特征相似，在遙感指數(shù)的識別上不易區(qū)分。應(yīng)用DEM數(shù)據(jù)的高程和坡度信息對舟山群島潮間帶粗提取結(jié)果進(jìn)行處理分析，可有效去除養(yǎng)殖塘等內(nèi)陸水體的干擾，進(jìn)一步區(qū)分灘涂和水體，提取潮間帶最終信息。

2.4 精度評價(jià)

混淆矩陣也稱誤差矩陣，通過將每個(gè)實(shí)測像元的位置和類別與分類結(jié)果中的相應(yīng)位置和類別相比較計(jì)算，可以評價(jià)某一類別地物是否被正確分類。具體可分為總體分類精度(overall accuracy，OA)、制圖精度(producer’s accuracy，PA)和用戶精度(user’s accuracy，UA)。OA指被正確分類的像元數(shù)與所有真實(shí)數(shù)據(jù)像元總數(shù)的比值。PA指被正確分類的某一類別像元數(shù)占該類別真實(shí)數(shù)據(jù)像元數(shù)的比例。UA指被正確分類的某一類別像元數(shù)占分類結(jié)果中該類所有像元數(shù)的比例。衡量分類結(jié)果與真實(shí)數(shù)據(jù)吻合度的一個(gè)指標(biāo)是Kappa系數(shù)，Kappa系數(shù)既考慮了被正確分類的像元，也考慮了各種漏分或錯分的像元，其范圍為[0，1]。所以當(dāng)OA，PA，UA和Kappa系數(shù)越大時(shí)，說明結(jié)果精度越高[49]。

3 結(jié)果與分析

3.1 潮間帶信息精度評價(jià)

本研究采用混淆矩陣法進(jìn)行精度評價(jià)，結(jié)合研究期內(nèi)的Landsat8影像信息并利用谷歌地球的高空間分辨率影像隨機(jī)生成樣本點(diǎn)對提取的潮間帶范圍進(jìn)行驗(yàn)證。利用目視解譯判別樣本點(diǎn)屬性，作為區(qū)分綠色植被和地表水體等內(nèi)陸地物與潮灘的依據(jù)，經(jīng)過谷歌地球影像檢查確保每個(gè)樣本點(diǎn)的地物類別屬性一致，沒有清晰的土地覆蓋信息的樣本點(diǎn)被剔除。最后，一共生成305個(gè)樣本點(diǎn)，其中灘涂85個(gè)，非灘涂(包括綠色植被和地表水體)220個(gè)，在研究區(qū)內(nèi)均勻分布，得到的混淆矩陣如表1所示。結(jié)果表明，2021年灘涂的UA為93.0%，PA為91.9%和OA為97.7%均較高，Kappa系數(shù)為0.95，舟山市的潮間帶提取結(jié)果與驗(yàn)證數(shù)據(jù)對比取得良好的一致性，實(shí)現(xiàn)了舟山群島潮間帶的快速自動提取。

表1 混淆矩陣及精度分析Tab.1 Confusion matrix and precision analysis

3.2 潮間帶信息空間分布特征

2021年舟山群島潮灘空間分布如圖3所示。統(tǒng)計(jì)計(jì)算得出，2021年舟山群島潮灘面積為35.19 km2。與現(xiàn)有的潮間帶提取方法對比，該方法不僅減少潮汐運(yùn)動引起的潮位不確定性和單時(shí)相遙感影像對潮間帶信息提取帶來的影響，并且不依賴于訓(xùn)練樣本、人工閾值的選取等，廣泛應(yīng)用的潛力較大。同時(shí)GEE云平臺的海量地理空間數(shù)據(jù)和強(qiáng)大的云計(jì)算能力為提取潮間帶信息提供強(qiáng)有力的支持，為用戶提供實(shí)時(shí)的計(jì)算模式，與傳統(tǒng)方法相比大大減少計(jì)算運(yùn)行的時(shí)間，有效提高潮間帶信息提取的時(shí)效性。

圖3 舟山群島潮間帶提取結(jié)果Fig.3 Extraction result of intertidal zone in Zhoushan Islands

舟山群島位于杭州灣外緣、錢塘江入?？?，風(fēng)浪和潮汐運(yùn)動對舟山海域的泥沙運(yùn)動和含沙量具有重要作用。舟山群島海岸線漫長曲折有利于形成可供開發(fā)利用的潮灘資源，且所形成的潮灘資源具有明顯的海島特征： 以淤泥質(zhì)灘涂為主，質(zhì)地軟弱； 為港灣型潮灘，多分布于港灣和島嶼之間，地勢平坦，具有優(yōu)越的開發(fā)利用條件，發(fā)展?jié)摿^大。灘涂資源分布面積廣，但分布不均勻，規(guī)模性差，潮灘的單塊面積較小，在一些陸域面積較大的海島沿岸地區(qū)呈塊狀分布，如岱山島(圖4(a))、舟山本島(圖4(b))以及六橫島(圖4(c))等。岱山島海岸坡度平緩，潮灘面積廣闊，多為淤泥質(zhì)灘涂，可促進(jìn)當(dāng)?shù)貫I海旅游的發(fā)展。舟山本島周圍的潮灘分布大多呈狹長型帶狀，外側(cè)與海相連接，內(nèi)側(cè)相接海塘或者壁岸，適用于圍墾造地和養(yǎng)殖業(yè)發(fā)展。六橫島擁有綿長的深水海岸線，灘涂連片發(fā)展，風(fēng)浪小以及氣候溫暖濕潤有利于養(yǎng)殖灘涂的發(fā)展。

(a) 岱山島(圖3a) (b) 舟山本島(圖3b) (c) 六橫島 (圖3c)圖4 舟山群島潮間帶局部提取結(jié)果Fig.4 Local extraction results of intertidal zone in Zhoushan Islands

4 結(jié)論與討論

為了減少潮汐運(yùn)動和單時(shí)相遙感影像對潮間帶信息提取帶來的影響，利用GEE云平臺提供的Landsat8 時(shí)序影像，提出基于融合DEM數(shù)據(jù)的決策樹算法，實(shí)現(xiàn)了舟山群島潮間帶信息的準(zhǔn)確、快速、自動提取。主要結(jié)論如下：

1)利用遙感指數(shù)和決策樹算法進(jìn)行潮間帶粗提取，為潮間帶信息的準(zhǔn)確提取奠定基礎(chǔ)。融合DEM數(shù)據(jù)的決策樹算法進(jìn)一步區(qū)分灘涂和水體，有效去除養(yǎng)殖塘等內(nèi)陸水體的干擾，提高潮間帶提取的準(zhǔn)確性，實(shí)現(xiàn)了舟山群島潮間帶信息的快速自動提取。

2)采用混淆矩陣法對提取的潮灘結(jié)果進(jìn)行精度評價(jià)，得到的舟山群島潮間帶提取的總體分類精度為97.7%，Kappa系數(shù)為0.95，具有良好的提取結(jié)果。

3)2021年舟山群島潮間帶面積為35.19 km2，舟山群島潮灘資源以淤泥質(zhì)灘涂為主，分布廣，在舟山本島、岱山島以及六橫島等陸域面積較大的海島沿岸地區(qū)呈塊狀分布。

在本研究中，利用GEE云平臺，基于融合DEM數(shù)據(jù)的決策樹算法是快速和自動進(jìn)行潮間帶信息提取的有效方法。GEE云平臺的海量地理空間數(shù)據(jù)和強(qiáng)大的云計(jì)算能力極大程度地降低數(shù)據(jù)的傳輸時(shí)間，大大減少處理數(shù)據(jù)所耗費(fèi)的人力物力，提高潮間帶信息提取的效率。融合DEM數(shù)據(jù)的決策樹算法能夠有效去除內(nèi)陸河流和湖泊等水體的影響，不依賴選取大量的訓(xùn)練樣本以及人工選取閾值?；贕EE云平臺最大程度地利用覆蓋研究區(qū)的所有Landsat影像，取得良好的提取結(jié)果。在更廣泛區(qū)域開展模型驗(yàn)證，并利用海量衛(wèi)星遙感數(shù)據(jù)進(jìn)行潮間帶長時(shí)間序列分析是下一步需要關(guān)注的研究內(nèi)容。