孫孟昊，蔡玉林，顧曉鶴，王雪麗，秦鵬

(1.山東科技大學(xué)，山東 青島 266590；2.基礎(chǔ)地理信息與數(shù)字化技術(shù)山東省重點實驗室，山東 青島 266590；3.北京農(nóng)業(yè)信息技術(shù)研究中心，北京 100097)

0 引言

近海岸地區(qū)是人類活動的密集區(qū)域，近海岸的人口密度幾乎是世界平均人口密度的3倍，人類活動的加劇對于沿海地區(qū)海岸線變化具有重要影響[1-2]。海岸線被國際地質(zhì)科學(xué)聯(lián)盟作為“地理指標(biāo)”之一，對海岸線進行準(zhǔn)確提取和監(jiān)測有助于解釋海岸變化過程，對海岸帶土地利用規(guī)劃、海岸侵蝕變化分析、海洋資源優(yōu)化管理等具有重要意義[3]。傳統(tǒng)的人工勘測方法具有周期長、時效性差、工作量大等缺點，應(yīng)用遙感技術(shù)進行海岸線監(jiān)測客觀性強、成本低廉，且具有較強的歷史追溯能力，可實現(xiàn)長時間序列的海岸線時空格局動態(tài)變化監(jiān)測。

近年來國內(nèi)外許多學(xué)者在利用衛(wèi)星影像數(shù)據(jù)快速解譯海岸線、利用多時相遙感影像分析海岸線時間變遷規(guī)律等方面開展了大量的研究。Komar等[4]分別對尼羅河三角洲、美國西北太平洋海岸、加利福尼亞州南部人工海灘以及堰洲島海岸等區(qū)域進行了海岸線的侵蝕與增長監(jiān)測研究；Purkis等[5]采用邊緣檢測和人工編輯相結(jié)合的方法，研究了迭戈加西亞島半個世紀(jì)的海岸線變化；Matghhany[6]利用TOPSAR數(shù)據(jù)通過建立波譜模型結(jié)合航拍數(shù)據(jù)和地面真實數(shù)據(jù)，探測了馬來群島Terengganu海岸線的侵蝕速率，找尋到海岸線侵蝕變化最快的地區(qū)；Ghosh[7]通過遙感(RS)和地理信息系統(tǒng)(GIS)的綜合技術(shù)，以水邊線作為海岸線監(jiān)測1989—2010年孟加拉國哈提亞島海岸線侵蝕和吸積的變化。楊金中等[8]利用長時間序列的Landsat MSS、TM等遙感影像，通過目視解譯結(jié)合地形圖資料對杭州灣地區(qū)海岸線進行提??；朱小鴿等[9]利用神經(jīng)網(wǎng)絡(luò)分類技術(shù)快速提取珠海口人工海岸水邊線；王琳等[10]利用改進的歸一化水體指數(shù)(modified normalized difference water index，MNDWI)提取廈門海域，進而使用紅綠變化檢測法分析其海域變化情況；樊彥國等[11]針對黃河三角洲地區(qū)淤泥質(zhì)海岸地區(qū)，以蝦池、公路等人工地物作為海岸線，采用圖像增強方法處理，提取出黃河口段和刁口段海岸線；常軍等[12]分析海岸砂不同含水量的光譜特征，進而利用多光譜影像提取海岸線；劉東等[13]提出了一種基于樣本選擇和支持向量機的算法，實現(xiàn)了遙感提取海岸線的自動提取、自動跟蹤，大大提高了自動化程度。Xu等[14]以平均大潮線作為海岸線，利用Landsat衛(wèi)星、中巴地球資源衛(wèi)星和環(huán)境一號衛(wèi)星數(shù)據(jù)相結(jié)合，對中國北方海岸線空間分布進行研究，并計算了分形維數(shù)等特征；陳正華等[15]利用四期衛(wèi)星數(shù)據(jù)以人工建筑(岸堤、海塘)作為海岸線，對浙江省大陸海岸線在1986—2009年間的動態(tài)變化進行了分析。

在上述研究中，由于潮汐數(shù)據(jù)獲取困難，部分學(xué)者忽略潮汐影響，以水邊線作為海岸線；部分學(xué)者為了研究人類活動對海岸的影響，以人工建筑線(蝦池、攔海堤壩、鹽田等)作為海岸線，以上這2種方法都欠妥。人工岸線可以以水邊線作為海岸線，基巖類海岸以陡峭巖石線作為海岸線，而其他類型海岸線提取時都應(yīng)做潮位校正。以本研究為例，潮汐影響最多可達上百米。而以人工建筑線作為海岸線的方法忽略了灘涂的變化，如自然情況下泥沙的淤積等。另有部分學(xué)者采用航拍數(shù)據(jù)和地面實測數(shù)據(jù)監(jiān)測海岸線變化，此方法精度高，但在時間上有一定局限性，且耗費人工；部分學(xué)者采用邊緣檢測等算法對特定地區(qū)的海岸線實現(xiàn)快速提取，此方法提取的仍是水邊線，且此類方法不具有普適性。因此，本文利用中分辨率衛(wèi)星數(shù)據(jù)提取海岸線，對缺少潮汐數(shù)據(jù)的年份，利用潮汐規(guī)律對其潮位校正，以期提高海岸線的遙感提取精度，進而實現(xiàn)長時間序列的海岸線時空格局變化監(jiān)測。

1 數(shù)據(jù)與方法

1.1 研究區(qū)域

研究區(qū)域位于中國山東省青島市，區(qū)域范圍在35°6′N～36°8′N，119°9′E～121°21′E之間。海岸帶北起栲栳灣，南至黃家塘灣。隨著我國沿海地區(qū)經(jīng)濟高度發(fā)展，使脆弱的海岸帶及其近海生態(tài)資源環(huán)境發(fā)生顯著變化[16]。青島是中國北方典型的經(jīng)濟發(fā)展快速的沿海城市，優(yōu)越的地理位置和狹長的海岸帶給青島提供了豐富的資源和條件，同時由于自然環(huán)境的變化和人類對海洋的開發(fā)，近幾十年間海岸線發(fā)生了明顯的變化，也對青島沿海區(qū)域的生態(tài)環(huán)境造成了嚴重的破環(huán)。

1.2 數(shù)據(jù)處理及信息提取

本研究使用的遙感影像數(shù)據(jù)為Landsat衛(wèi)星影像、GDEMV2 30 m分辨率的數(shù)字高程影像以及天地圖的亞米級高空間分辨率歷史影像。研究時間跨度為1979—2019年，時間間隔為10年(表1)。圖像的預(yù)處理包括輻射定標(biāo)、大氣校正、幾何校正、裁剪拼接。

表1 影像數(shù)據(jù)

1.3 海岸線提取

提取海岸線首先要進行水陸分離，徐涵秋[17]提出改進的歸一化水體指數(shù)(MNDWI)，該指數(shù)相比于歸一化水體指數(shù)(normalized difference water index，NDWI)能削弱地表高亮地物(如建筑物)的影響。本研究區(qū)域中存在大量人工海岸，適合采用MNDWI指數(shù)。利用MNDWI指數(shù)結(jié)合閾值分割法實現(xiàn)水陸分離。對濕度較大的非海洋區(qū)域如水塘、魚池等小塊區(qū)域，采用中值濾波的方法消除其影響，最后提取出水邊線。

遙感影像所提取的水邊線為衛(wèi)星過境時間點的瞬時水邊線(圖1)，不是海岸線，潮汐對于水邊線提取存在較大影響，有必要分析潮汐日變化規(guī)律，對瞬時水位線進行潮位校正[18]。

圖1 瞬時水邊線示意圖

1.4 潮位校正

海水在天體(主要是太陽和月亮)的引潮力作用下產(chǎn)生潮汐現(xiàn)象，潮汐現(xiàn)象具有規(guī)律性，古代科學(xué)家就已發(fā)現(xiàn)潮汐現(xiàn)象和月亮的圓缺有關(guān)[19]。農(nóng)歷歷法就是以月相變化為依據(jù)的，實際上古人按照農(nóng)歷歷法總結(jié)了潮汐規(guī)律，如農(nóng)歷每月初一和十五出現(xiàn)大潮，每月上下弦出現(xiàn)小潮[20]?，F(xiàn)代學(xué)者丁奕、趙敬瑛等分別結(jié)合錢塘江和大遼河的潮位資料得出了兩地的潮汐規(guī)律[21-22]。為了研究潮汐規(guī)律做潮位校正的可行性，本文對2017—2019年農(nóng)歷六月十八潮汐數(shù)據(jù)(http：//www.chinaports.com)進行了數(shù)學(xué)分析。

青島屬于半日潮型即一日有2次高潮和低潮，如表2獲取2組高潮和低潮的時刻和潮高數(shù)據(jù)，對高低潮時刻和潮高做相關(guān)性分析，相關(guān)系數(shù)均達到0.99以上。用每個時刻潮位數(shù)據(jù)擬合潮位曲線，發(fā)現(xiàn)不同年分潮位曲線具有相同的走勢(圖2)。計算3條曲線同一時刻任意2條之間的潮位高度差，計算所有潮位高度差平均值為28 cm。分析歷年農(nóng)歷同月同日的潮位數(shù)據(jù)，同樣具有高度相似性，因此可作為潮汐數(shù)據(jù)缺失年份的潮位修正依據(jù)。

表2 2組低潮和高潮數(shù)據(jù)

圖2 擬合2017—2019年同一天的潮位曲線

本研究以日最高潮位線作為海岸線，采用外牽引法[23]對淤泥質(zhì)海岸和砂質(zhì)海岸做潮位校正。根據(jù)衛(wèi)星過境時潮位和日最高潮位計算出潮位差，再由潮位差和海岸平均坡度計算外牽引距離(圖3)，計算公式如下：

L=(H2-H1)/Tan(α)

(1)

式中：α為海岸平均坡度；H1為瞬時潮位高度；H2為最高潮位高度；h為潮位差；L為外牽引距離。

圖3 計算外牽引距離示意圖

由DEM圖像生成坡度圖(圖4)，計算從2009年水邊線向陸地推進80 m的范圍內(nèi)平均坡度為2.1°。

對數(shù)據(jù)缺失年份采用有數(shù)據(jù)年份的農(nóng)歷同月同日潮汐數(shù)據(jù)校正，由上述2017—2019年農(nóng)歷同日同時刻平均潮位高差為28 cm，根據(jù)三角函數(shù)關(guān)系計算此誤差在水平方向上為7.6 m，而衛(wèi)星分辨率為30 m，遠遠小于一個像元寬度，誤差在合理范圍內(nèi)。因此提取海岸線時，利用潮汐規(guī)律可以對缺失數(shù)據(jù)的年份做潮位修正。

圖4 研究區(qū)域DEM圖

研究區(qū)域內(nèi)海岸類型有人工海岸、基巖海岸、砂質(zhì)海岸和淤泥質(zhì)海岸。人工海岸和基巖海岸水陸分明、坡度大而陡，潮汐影響小不需要做潮位校正。僅對砂質(zhì)海岸和淤泥質(zhì)海岸做潮位校正，計算潮位校正外牽引距離結(jié)果如表4所示。

表4 潮位數(shù)據(jù)以及計算結(jié)果

1.5 精度驗證

對衛(wèi)星影像中模糊區(qū)域使用較高分辨率天地圖臨近日期的歷史影像進行目視檢驗。以2019年提取結(jié)果為例，圖5左邊是不需潮位校正的人工海岸提取結(jié)果；右邊是砂質(zhì)海岸潮位校正前后對比，其中黃色為瞬時水邊線，紅色為校正后的海岸線。

在天地圖臨近歷史影像中觀察缺少潮位數(shù)據(jù)年份的潮位校正情況，也具有較高精度。選取的是天地圖2010年1月16日歷史影像，與研究使用的Landsat2009年影像相差41天，如圖5黃色為瞬時水邊線，紅色為海岸線。2010年1月16日為農(nóng)歷臘月初二，Landsat影像獲取日期2009年12月6日，為農(nóng)歷十月二十，從潮汐數(shù)據(jù)得知2018年農(nóng)歷十月二十和臘月初二的最高潮位的潮位差為0.3 m，說明可以用天地圖2010年1月16日影像檢驗。如圖6所示海岸線位置恰好在積雪位置附近，積雪線所處的位置為最高潮位線，說明提取結(jié)果精度較高。

圖5 不需潮位校正的人工海岸提取結(jié)果及砂質(zhì)海岸潮位校正前后對比

圖6 潮位數(shù)據(jù)缺失年份潮位校正結(jié)果

2 結(jié)果與分析

2.1 青島海岸線變化分析

海岸線提取最終結(jié)果如圖7所示。提取結(jié)果可以客觀地反應(yīng)出青島海岸線近40年的變遷情況。由圖8可看出，1979—2019年間海岸線整體呈向海洋方向推進的趨勢，總體向海洋推進158 km2。從栲栳灣到橫門灣大部分屬于基巖海岸，岸線穩(wěn)定；鰲山灣由于該區(qū)域海水養(yǎng)殖業(yè)發(fā)展和青島國際博覽中心經(jīng)濟圈的形成，使得該地區(qū)岸線向海洋推進；鰲山灣到膠州灣之間有嶗山風(fēng)景區(qū)，多為基巖海岸，岸線穩(wěn)定；膠州灣變化最為明顯，由于灣內(nèi)人工造陸、圍海養(yǎng)殖、港口建設(shè)、城市擴建、鹽田開發(fā)，導(dǎo)致灣內(nèi)面積不斷縮小，岸線由自然曲折的自然岸線演變?yōu)槠街钡娜斯ぐ毒€；早期唐島灣由于自然泥沙淤積和發(fā)展養(yǎng)殖業(yè)，海岸線有較明顯的向海洋發(fā)展趨勢，近十年來唐島灣發(fā)展旅游業(yè)岸線變化也十分明顯；黃家塘灣至陳家貢灣由于發(fā)展交通業(yè)，大力擴建和修建港口，導(dǎo)致灣內(nèi)面積減少。

圖7 海岸線提取結(jié)果(1979—2019年)

圖8 海岸線提取結(jié)果

2.2 不同時間段典型區(qū)域海岸線變化分析

1)在1979—1989年間，海岸線總體向海洋推進6 km2。在此期間岸線整體穩(wěn)定，岸線變化主要由于泥沙淤積等自然因素的引起。膠州灣東岸小部分區(qū)

域經(jīng)過人工改造，導(dǎo)致了原來曲折自然岸線部分轉(zhuǎn)化為平直的人工岸線(圖9(a))。

2)在1989—1999年期間，海岸線共向海洋推進30 km2。膠州灣是變化面積最大區(qū)域(圖9(b))，在此期間膠州灣由于少量泥沙淤積、人工造陸發(fā)展旅游業(yè)以及修建沿海公路和攔海堤等城市工程，海岸線共向海洋方向推進23 km2。尤其膠州灣東部岸線基本已全部演變?yōu)槠街钡娜斯ぐ毒€，岸線趨于穩(wěn)定。唐島灣變化也較為突出，由于泥沙淤積海岸線向海洋推進2 km2。

3)在1999—2009年期間，總共向海洋推進78 km2，是40年間面積變化最大的階段。變化最突出區(qū)域仍然是膠州灣，此期間膠州灣開發(fā)西南部(黃島)。十年間黃島區(qū)變化劇烈，該區(qū)域填海造陸、修建港口幾乎完全被改造為人工海岸，海岸線大幅度向海洋方向推進變化極為明顯(圖9(c))。

圖9 典型變化區(qū)域

唐島灣在此期間由于灣內(nèi)清理淤泥，海岸線向陸地方向移動。從天地圖的1999年和2009年歷史影像中可以看出唐島灣變化情況(圖10)。同時在此期間灣內(nèi)修建了沿海公路、欄海堤和一些公園，岸線成分以人工岸線為主。其他區(qū)域岸線穩(wěn)定沒有發(fā)生明顯的岸線侵蝕和吸積。

圖10 1999年、2009年唐島灣對比圖

4)在2009—2019年間，海岸線向海洋推進44 km2。由于膠州灣和唐島灣之前的沿海公路和攔海堤的建設(shè)，岸線趨于穩(wěn)定。該區(qū)域內(nèi)主要變化是膠州灣北部水產(chǎn)養(yǎng)殖區(qū)的發(fā)展和唐島灣西部人工造陸修建了星光島旅游區(qū)(圖9(d))。

在此期間變化最明顯區(qū)域由中部城市化程度高的膠州灣地區(qū)，轉(zhuǎn)移到了南部的黃家塘灣至陳家貢灣岸段。應(yīng)國家交通部和山東省人民政府要求，在此區(qū)域填海造陸大力擴建和修造港口，岸線大幅向海洋推進(圖9(e))。

3 結(jié)束語

本文利用數(shù)學(xué)分析的方法對不同年份的農(nóng)歷同日期潮汐數(shù)據(jù)對比分析，發(fā)現(xiàn)其農(nóng)歷同日期日潮汐高度走勢具有高度相似性。將這一規(guī)律應(yīng)用在海岸線潮位校正上，提出了利用現(xiàn)有潮汐數(shù)據(jù)校正農(nóng)歷同日期的往期遙感影像的方法。利用此方法對青島海岸線進行提取并分析其變化，結(jié)果表明在1979—2019年間青島海岸線整體向陸地推進158 km2。變化最大的地區(qū)為膠州灣，變化最劇烈的時間段為1999—2009年。從整體變化趨勢看，膠州灣地區(qū)修建了沿海公路和攔海堤，海岸線已趨于穩(wěn)定。近十年來青島南部的黃家塘灣至陳家貢灣岸段，因投資發(fā)展交通業(yè)使海岸線變化劇烈。未來交通業(yè)可能為其帶來巨大經(jīng)濟效益，加快其城市化進程，因此未來將會是一個變化活躍的區(qū)域。青島北部岸線有大量基巖海岸，其他類型海岸功能多為養(yǎng)殖業(yè)、旅游業(yè)，未來將繼續(xù)保持穩(wěn)定，海岸線仍以較慢的速度向海洋方向推進。

前人對水邊線潮位校正距離是根據(jù)岸灘坡度和驗潮站潮高數(shù)據(jù)的三角函數(shù)關(guān)系推算而得[24]。岸灘坡度一般通過DEM影像或地形數(shù)據(jù)獲得；也有學(xué)者根據(jù)不同時相的2幅影像提取的水邊線間距，結(jié)合驗潮站2個時間的潮位差計算灘涂坡度[25-26]，總之都需要有此時刻潮汐數(shù)據(jù)支持。對于久遠年份的海岸線提取，潮高數(shù)據(jù)有無便成了一個決定性因素。本文針對潮位數(shù)據(jù)缺失這一難題，提出的基于潮汐規(guī)律修正的海岸線長時序遙感監(jiān)測方法，提取結(jié)果誤差遠小于中等分辨率衛(wèi)星的像元寬度。此方法操作簡單可行、精度高，解決了無潮汐數(shù)據(jù)年份的潮位校正問題。另外本文使用MNDWI指數(shù)提取水邊線，下一步有待利用神經(jīng)網(wǎng)絡(luò)和邊緣檢測算子的方法，提高水邊線提取的自動化程度，提高效率減少人工勞動。