秦艷芳，葉 琳，陳思明

（閩江學(xué)院海洋學(xué)院，福州 350108）

快速而又準(zhǔn)確地測定海岸線的動態(tài)變化，適時關(guān)注沿海地區(qū)岸線變遷具有十分重要的研究價值和現(xiàn)實意義[1]。廈門市自成立經(jīng)濟(jì)特區(qū)以來，經(jīng)濟(jì)持續(xù)高速發(fā)展，大范圍填海造地，使其海岸線發(fā)生了顯著變化。近年來，作為連接“一帶”和“一路”的主要陸海樞紐城市，廈門市積極建設(shè)與各國家和地區(qū)互聯(lián)互通、經(jīng)濟(jì)貿(mào)易以及人文文化交流的重點地區(qū)，如擴(kuò)張新機(jī)場、推動口岸基礎(chǔ)設(shè)施建設(shè)等[2]。這一系列舉措對廈門市城鎮(zhèn)建設(shè)、產(chǎn)業(yè)布局以及港口經(jīng)濟(jì)發(fā)展均產(chǎn)生重要影響，也對海岸線及沿岸生態(tài)環(huán)境的變化產(chǎn)生深遠(yuǎn)影響。

遙感具有大范圍實時同步獲取不同尺度時空信息的優(yōu)勢，在海岸線調(diào)查、動態(tài)監(jiān)測等方面表現(xiàn)出極大的優(yōu)勢。國內(nèi)外學(xué)者就海岸線遙感提取方法已經(jīng)開展了大量的研究，并提出多種算法。Ryan等對陸地和水體的遙感分類時引入神經(jīng)網(wǎng)絡(luò)算法，完成了海岸線提取[3]。Alesheikh等利用單一波段直方圖閾值和波段比值運算相結(jié)合，提取了烏爾米耶湖的海岸線[4]。劉鵬程利用水平集理論對海岸線進(jìn)行提取，引入海岸邊界區(qū)域特征分割算法[5]。胡亞斌等采用NDWI方法來分割水陸信息，提取山東半島海陽沙灘的砂質(zhì)岸線[6]。方亮等提出了基于均值漂移的海島（礁）岸線快速提取方法[7]。目前，基于遙感的海岸線提取方法已基本實現(xiàn)自動化，并引入了許多新的算法，在一定程度上提高了提取精度，但提取結(jié)果仍然需要人工交互的方式進(jìn)行修訂。

本文采用Landsat系列遙感影像，通過多波段譜間關(guān)系法分別提取1999、2005、2011、2017年廈門市海岸線數(shù)據(jù)，分別從岸線長度和陸地面積變化來分析其變遷特征和主要驅(qū)動因素。

1 研究區(qū)概況

廈門市位于118°04′04″E、24°26′46″N附近，處于福建省的東南端，九龍江的入海口，并與臺灣島隔海相望，是我國東南沿海對外貿(mào)易的主要港口，也是一個風(fēng)景旅游城市。研究區(qū)可分為島內(nèi)和島外兩部分，島內(nèi)地區(qū)包含了思明區(qū)和湖里區(qū)，島外地區(qū)包含了海滄區(qū)、集美區(qū)、同安區(qū)以及翔安區(qū)，陸地總面積約1 600 km2，海域面積約390 km2。

2 數(shù)據(jù)與研究方法

2.1 數(shù)據(jù)來源

本研究以Landsat系列遙感數(shù)據(jù)作為主要數(shù)據(jù)源，分別獲取廈門市1999年9月7日ETM+數(shù)據(jù)、2005年7月13日和2011年9月16日的TM數(shù)據(jù)、2017年10月2日OLI數(shù)據(jù)共4個時相的遙感影像，云量均不超過10%。各期遙感影像均經(jīng)過大氣校正和增強(qiáng)處理，并以2017年OLI影像為準(zhǔn)對其他3期進(jìn)行幾何精校正以保證數(shù)據(jù)可比性。

此外，還使用了2014年1:10萬福建省鄉(xiāng)鎮(zhèn)行政區(qū)劃矢量圖，用于提取廈門市行政區(qū)劃的內(nèi)陸邊界。

2.2 海岸線的遙感提取方法

本文采用多波段譜間關(guān)系法進(jìn)行海岸線提取，該方法由周成虎等提出[8]。通過分析水體與其他地物在TM影像不同波段上的灰度曲線圖，筆者發(fā)現(xiàn)，水體具有綠波段加紅波段大于近紅外波段加中紅外波段的特征，與其他類型均不相同，由此提出了提取水體的譜間關(guān)系法模型：

式 中，TM2、TM3、TM4和 TM5分 別 對 應(yīng) TM、ETM+影像中的2、3、4、5波段以及OLI影像的3、4、5、6波段；x是能夠區(qū)別水體與陸地的閾值（取值約為0）。

將基于閾值分離的水陸二值圖像通過柵格轉(zhuǎn)矢量、面轉(zhuǎn)線即可自動提取海陸邊界，再經(jīng)平滑處理可消除柵格轉(zhuǎn)矢量特有的鋸齒特征，最后通過人工交互的方法對各期海岸線進(jìn)行修正。

2.3 海岸線變遷特征分析

2.3.1 長度變化

按廈門市行政區(qū)劃，將研究區(qū)海岸線分為海滄區(qū)、集美區(qū)、同安區(qū)、翔安區(qū)、湖里區(qū)和思明區(qū)6段，通過各時期分段長度變化來分析海岸線長度變化的情況。

為了能夠更加客觀、準(zhǔn)確地對海岸線進(jìn)行時空變化分析對比，采用某一時間段內(nèi)海岸線長度的年均變化來表示海岸線的變化速率。

2.3.2 陸地面積變化

分別將4期海岸線與廈門市行政區(qū)劃的內(nèi)陸邊界線進(jìn)行合并，形成整個廈門市各期行政區(qū)劃的陸地邊界，再由線轉(zhuǎn)面并計算面積，從而獲取4期廈門市的陸地總面積和3個時期陸地面積的增大或減少情況。

3 結(jié)果分析

4期廈門市海岸線如圖1所示，海岸線主要變化的區(qū)域集中在海滄區(qū)、同安區(qū)、湖里區(qū)和翔安區(qū)，特別是2011-2017年，大嶝島陸地面積向外擴(kuò)張非常明顯。

圖1 1999-2017年廈門市地區(qū)海岸線空間分布

3.1 長度變化

1999-2017年間廈門市地區(qū)海岸線總長度由235.16 km增至264.98 km，其中在2005-2011年和2011-2017年，分別增長了11.05 km、21.8 km，年變化速率分別為1.84 km/a和3.63 km/a，如表1所示。

結(jié)合圖1，1999-2005年，海岸線總長度減少3.03 km，主要體現(xiàn)在海滄區(qū)和湖里區(qū)的岸線曲折程度明顯下降；2005-2011年海岸線向海域方向明顯移動，曲折度增大，主要變化集中在翔安區(qū)；2011-2017年，除翔安區(qū)大嶝島海岸線大幅度向海域方向移動外，其他地區(qū)變化不明顯。

表1 海岸線長度變化及年平均變化速率

下面分別從各分區(qū)進(jìn)行分析。

翔安區(qū)2017年海岸線總長度為103.74 km，較1999年增長20.06 km。增長最顯著的時期為2011-2017年，共增長14.6 km，年變化速率為2.43 km/a，主要表現(xiàn)為大嶝島在此期間急劇擴(kuò)張。2005-2011年，海岸線長度也增加了9 km，岸線曲折度有所上升，東南部岸線明顯向海域方向移動。此外，1999-2005年局部地區(qū)海岸線由曲折變平直，總長度下降3.54 km。

湖里區(qū)2017年海岸線總長度為39.91 km，較1999年增長8.11 km。1999-2005年、2005-2011年的變化較為顯著，長度增長率分別為22.14%、6.23%，年變化速率分別為1.17 km/a、0.40 km/a。其間變化主要集中表現(xiàn)在高殿村、廈門高崎機(jī)場以及鐘宅村附近沿海方向的陸地不斷擴(kuò)建。

同安區(qū)2017年海岸線總長度為19.51 km，較1999年增長3.04 km。1999-2005年和2005-2011年，海岸線長度不增反減，主要與海岸線曲折度下降有關(guān)。而2011-2017年，海岸線曲折度有所上升，長度也增長了5.72 km，主要集中在大埔附近沿海地區(qū)的填海造地工程。

海滄區(qū)2017年海岸線總長度為48.36 km，較1999年減少了0.91 km。雖然總長度變化不大，但在1999-2005年有明顯降低，主要原因是鐘山村、石塘村等地圍填海工程使得岸線的曲折度大幅下降；而在2011-2017年，海岸線長度又有所上升。

此外，集美區(qū)和思明區(qū)的岸線沒有明顯變化。

3.2 陸地面積變化

1999-2017年，廈門市陸地面積從1 564.93 km2增至1 600.95 km2，總體增加36.02 km2。如表2、圖2所示，1999-2005年，陸地面積增加了9.05 km2，主要分布在海滄區(qū)和湖里區(qū)；2005-2011年，陸地面積增加了15.60 km2，全地區(qū)分布總體較為平均；2011-2017年，陸地面積增加了11.36 km2，主要集中在翔安區(qū)。

表2 1999-2017年陸地面積變化統(tǒng)計

圖2 局部地區(qū)陸地變化

從各分區(qū)來看，湖里區(qū)的新增陸地較為顯著。2017年陸地面積為69.58 km2，較1999年增長了7.72 km2，其中1999-2005年和2005-2011年間分別增加了 4.46 km2、3.69 km2。

此外，海滄區(qū)和翔安區(qū)陸地面積也迅速擴(kuò)張，主要原因是翔安區(qū)的大嶝島陸地面積增加了將近一倍，前浯、澳頭等附近海岸也有明顯擴(kuò)展。1999-2017年海滄區(qū)的陸地面積從170.1 km2增至177.97 km2，增長率為4.63%；翔安區(qū)從350.09 km2增至365.89 km2，增長率為4.51%。在此過程中，海滄區(qū)的陸地面積持續(xù)增加，主要集中在2011年之前；而翔安區(qū)在2005年之前陸地面積變化不明顯，之后開始逐漸增加，特別是在2011-2017年共增加了11.21 km2，增長率為3.16%。

4 結(jié)論

本文選用4期Landsat遙感數(shù)據(jù)，對廈門市近20年海岸線變遷特征進(jìn)行分析，得出以下主要結(jié)論。

（1）1999-2017年，廈門市地區(qū)海岸線增長趨勢十分顯著，特別是在2005-2011年以及2011-2017年2個時間段里分別增長了11.05 km、21.8 km，其年平均變化速率分別為1.84 km/a和3.63 km/a。從行政分區(qū)上來看，顯著表現(xiàn)為2005-2017年在翔安區(qū)的增長和1999-2005年在湖里區(qū)的增長。

（2）廈門市陸地面積總體增加了36.02 km2，集中表現(xiàn)在翔安區(qū)、海滄區(qū)和湖里區(qū)。城鎮(zhèn)建設(shè)等人類活動是廈門市海岸線發(fā)生顯著變化的重要原因。1999年以來，在海滄區(qū)、翔安區(qū)和湖里區(qū)，沿海城鎮(zhèn)建設(shè)用地不斷擴(kuò)張，港口進(jìn)一步擴(kuò)建，翔安區(qū)東南側(cè)水產(chǎn)養(yǎng)殖區(qū)域也有明顯擴(kuò)大，特別是翔安區(qū)大嶝島圍填海工程用于大嶝航空城基地建設(shè)，造成的海岸線變化最為顯著。

（3）近年來，為響應(yīng)“一帶一路”倡議，廈門市經(jīng)濟(jì)建設(shè)的步伐進(jìn)一步擴(kuò)大，在經(jīng)濟(jì)文化方面的對外交流愈發(fā)頻繁，這也增加了廈門市在客運和貨運方面的供給壓力，使得政府積極擴(kuò)建船埠、港口和機(jī)場，從而對廈門市海岸線變遷產(chǎn)生了重要影響。