王詩(shī)洋， 楊武年， 佘金星

(成都理工大學(xué)　地學(xué)空間信息技術(shù)國(guó)土資源部重點(diǎn)實(shí)驗(yàn)室 ,成都　610059)

?

我國(guó)南海沿岸Landsat影像海岸線提取與變化分析

王詩(shī)洋， 楊武年*， 佘金星

(成都理工大學(xué)地學(xué)空間信息技術(shù)國(guó)土資源部重點(diǎn)實(shí)驗(yàn)室 ,成都610059)

摘要：海岸線是海域變化及陸地變化研究的重要參考依據(jù)。為了精確獲取我國(guó)南海沿岸近十年的變化，這里以TM、OLI數(shù)據(jù)為基礎(chǔ)，利用重復(fù)自組織數(shù)據(jù)分析技術(shù)(ISODATA)自動(dòng)提取海岸線，提取結(jié)果采用區(qū)域統(tǒng)計(jì)分析進(jìn)行精度驗(yàn)證。研究結(jié)果表明，近十年來(lái)，我國(guó)南海沿岸河流湖泊面積呈快速縮減趨勢(shì)。

關(guān)鍵詞：海岸線； 南海； ISODATA； 區(qū)域統(tǒng)計(jì)

0引言

海岸線地理位置隨著時(shí)間的變化而變遷，一定程度地體現(xiàn)了人類(lèi)活動(dòng)的變化以及自然因素的影響，研究海岸線的變遷對(duì)沿岸地區(qū)的發(fā)展十分重要。中國(guó)南海是位于我國(guó)南部的陸緣海，利用遙感影像提取南海沿岸地區(qū)的海岸線并加以解釋判讀，對(duì)分析南海海域近十年變化具有重要意義。利用遙感影像提取海岸線已有很多先例，陳曉英等[1]利用MSS 等數(shù)據(jù)對(duì)海州灣進(jìn)行海岸線的提取，并得出1973年-2003年間海州灣大陸海岸線長(zhǎng)度在波動(dòng)中呈增長(zhǎng)趨勢(shì)、海灣面積呈減少趨勢(shì)；王李娟等[2]應(yīng)用ETM影像進(jìn)行黃河三角洲地區(qū)海岸線的提取，并比較了修復(fù)歸一化水體指數(shù)及Sobel方法的提取結(jié)果，其結(jié)果表明，Sobel較MNDWI提取精度更高一些。這里選用TM以及OLI影像，其中l(wèi)andsat8 OLI數(shù)據(jù)融合后影像分辨率較高。

遙感影像的海岸線提取方法，有目視解譯及計(jì)算機(jī)自動(dòng)提取方法。傳統(tǒng)海岸線目視解譯提取方法，即通過(guò)目視解譯來(lái)提取增強(qiáng)后遙感影像的海岸線，王建步等[3]通過(guò)建立海岸線分類(lèi)體系以及遙感解譯標(biāo)志對(duì)遼河口海岸線進(jìn)行提取，并分析得出1979年-2013年間遼河口海岸線長(zhǎng)度總體呈增長(zhǎng)趨勢(shì)。目視解譯提取方法主觀性太強(qiáng)，數(shù)據(jù)量較大、且處理時(shí)間過(guò)長(zhǎng)，因此應(yīng)多采用計(jì)算機(jī)自動(dòng)提取海岸線的方法。謝明鴻等[4]通過(guò)種子點(diǎn)增長(zhǎng)方法對(duì)雷達(dá)影像進(jìn)行自動(dòng)化的提取，這種計(jì)算機(jī)自動(dòng)提取法主要依靠計(jì)算機(jī)進(jìn)行海岸線提取，對(duì)各類(lèi)型海岸線的提取缺乏多樣性。這里采用非監(jiān)督分類(lèi)方法來(lái)提取海岸線，非監(jiān)督即依據(jù)光譜相似性進(jìn)行集群分類(lèi)，分類(lèi)速度快且客觀性較強(qiáng)。ISODATA(Iterative Self-Organizing Date Analysis Technique)方法是基于迭代的非監(jiān)督分類(lèi)法。采用ISODATA方法對(duì)影像進(jìn)行分類(lèi)能夠獲得矢量結(jié)果數(shù)據(jù)，定義ISODATA分類(lèi)后數(shù)據(jù)，對(duì)定義為海岸線的類(lèi)別即分類(lèi)提取的海岸線依據(jù)不同類(lèi)型遙感解譯標(biāo)志進(jìn)行修改。通過(guò)空間分析進(jìn)行海岸線提取結(jié)果的驗(yàn)證，結(jié)果表明，ISODATA提取的海岸線精度較高。為進(jìn)行南海海岸線變化監(jiān)測(cè)，分析所提取2003年-2014年期間海岸線矢量數(shù)據(jù)，統(tǒng)計(jì)分析海岸線變遷。

1數(shù)據(jù)

1.1研究區(qū)概況及數(shù)據(jù)來(lái)源

研究區(qū)選為南海海域靠近海南省東部沿岸區(qū)域，其沿岸海域自海南省東北部文昌市至海南省東南部陵水黎族自治縣。海南省地處熱帶北緣，雨量充沛，因此多云霧天氣。為研究近年來(lái)南海海域面積變化，選取2003年、2009年云量較小的Landsat5 TM影像以及2014年的Landsat8 OLI影像，三幅影像均采用WGS84坐標(biāo)系統(tǒng)，下載數(shù)據(jù)來(lái)源于Image courtesy of the U.S. geological survey，數(shù)據(jù)參數(shù)見(jiàn)表1。Landsat8于2013年2月11日發(fā)射升空經(jīng)過(guò)100天測(cè)試運(yùn)行后開(kāi)始獲取影像，其中Landsat8所搭載的OLI傳感器，有9個(gè)波段，空間分辨率為30 m，其中pan波段15 m。

表1　衛(wèi)星參數(shù)表

1.2數(shù)據(jù)處理

數(shù)據(jù)處理平臺(tái)為ENVI、ARCGIS。ENVI進(jìn)行數(shù)據(jù)預(yù)處理，即三幅影像的輻射定標(biāo)，以及多光譜影像數(shù)據(jù)的大氣校正。下載的landsat8數(shù)據(jù)包含pan波段，對(duì)Landsat8數(shù)據(jù)進(jìn)行pan波段融合，融合方法采用gram-schmidt pan sharpening方法，得到分辨率15 m的融合圖像，其中兩幅TM影像均采用3、2、1波段組合進(jìn)行海岸線提取。

2提取方法

海岸線提取方法有圖像差分算法、水體指數(shù)法以及邊緣檢測(cè)法等。這里采用ISODATA法以及歸一化水體指數(shù)兩種方法，分別對(duì)南海海域沿岸地區(qū)進(jìn)行海岸線的提取。

2.1非監(jiān)督分類(lèi)

非監(jiān)督分類(lèi)即多光譜圖像的聚類(lèi)分析，是計(jì)算機(jī)自動(dòng)根據(jù)地物光譜相似度特征進(jìn)行聚類(lèi)的方法。非監(jiān)督分類(lèi)無(wú)法對(duì)屬性進(jìn)行判定，分類(lèi)后需要定義類(lèi)別。

ENVI進(jìn)行非監(jiān)督分類(lèi)時(shí)，選用ISODATA[5]方法即重復(fù)自組織數(shù)據(jù)分析。ISODATA方法是在K均值算法的基礎(chǔ)上進(jìn)行的，首先設(shè)定最大類(lèi)別數(shù)量，通過(guò)計(jì)算周?chē)裨辆垲?lèi)中心的距離而劃分類(lèi)別，迭代計(jì)算直至達(dá)到給定條件而停止計(jì)算結(jié)束分類(lèi)，此法降低了人工目視解譯的主觀性誤差。

ISODATA法對(duì)研究影像劃分類(lèi)別以后，得到影像分類(lèi)矢量數(shù)據(jù)，定義類(lèi)別并選擇劃分為海岸帶的類(lèi)，依據(jù)不同類(lèi)型海岸線解譯標(biāo)志對(duì)一些分類(lèi)有偏差的區(qū)域進(jìn)行修改。

2.2歸一化水體指數(shù)

歸一化水體指數(shù)是通過(guò)波段計(jì)算而凸顯水體信息的方法，NDWI公式[6]為式(1)。

NDWI=(Green-NIR)/(Green+NIR)

(1)

其中：Green代表綠波段；NIR代表近紅外波段。

在Landsat8 OLI影像上band3和band5，在Landsat5 TM影像上分別為band2和band4。由于水體反射率在近紅外接近為零，在綠波段可見(jiàn)光范圍內(nèi)反射率較強(qiáng)，因此NDWI影像上的水體呈亮色調(diào)。而植被在近紅外反射率很強(qiáng)，因此NDWI影像中植被為暗色調(diào)，由公式(1)可知水體信息在NDWI影像上被凸顯。采用NDWI進(jìn)行海岸線提取的方法，只能依靠目視解譯提取海岸線，主觀性較強(qiáng)且較依賴(lài)影像的空間分辨率，存在較大的主觀誤差。對(duì)2003年-2009年、2009年-2014年時(shí)限影像處理后分別進(jìn)行海岸線提取等工作，其流程見(jiàn)圖1。

3結(jié)果分析

3.1海岸線提取

海岸線的提取分別采用NDWI及ISODATA方法進(jìn)行海岸線提取。在NDWI影像中，根據(jù)海岸線解譯標(biāo)志，對(duì)研究區(qū)的砂質(zhì)海岸、淤泥質(zhì)岸線等進(jìn)行海岸線解譯。砂質(zhì)岸線、人工岸線等不同類(lèi)型海岸線在遙感影像中所顯示的解譯標(biāo)志不同，因此，各類(lèi)型海岸線分別進(jìn)行遙感解譯[7]。研究區(qū)大部分區(qū)域?yàn)樯百|(zhì)岸線，遙感解譯砂質(zhì)岸線較平直，遙感影像中，即位于高亮度干燥砂灘及暗色調(diào)濕潤(rùn)砂灘之間[8]。在遙感影像中，淤泥質(zhì)海岸線多位于疏密植被之間[9]。基巖海岸在遙感影像的灰度值高于水體及淤泥質(zhì)海岸，有植被覆蓋[10]。人工海岸是由水泥等材質(zhì)構(gòu)建，在遙感影像上反射率較高的人工海岸與低反射率海水區(qū)分明顯可以確定其海岸線[11]。

采用ISODATA方法提取海岸線，速度快，能很好地將各類(lèi)別海岸線提取出來(lái)。

圖1　流程圖Fig．1　Flow chart

3.2海岸線提取結(jié)果分析與驗(yàn)證

采用空間分析方法對(duì)ISODATA提取的海岸線進(jìn)行精度驗(yàn)證，空間分析可以通過(guò)研究空間對(duì)象

的空間關(guān)系等來(lái)定量化表述空間事物[12]。對(duì)已提取的海岸線生成向海、向陸緩沖區(qū)，由于向海一側(cè)為海洋，向陸一側(cè)均為陸地植被等，因此在向陸一側(cè)緩沖區(qū)的標(biāo)準(zhǔn)差較大，向海一側(cè)緩沖區(qū)標(biāo)準(zhǔn)差較小。取緩沖區(qū)距離為影像一個(gè)像元，設(shè)定向陸、海緩沖區(qū)分別為區(qū)域A、區(qū)域B，由區(qū)域統(tǒng)計(jì)分析可知區(qū)域A、B的標(biāo)準(zhǔn)差，為便于比較，計(jì)算兩區(qū)域標(biāo)準(zhǔn)差的均值并分別記為MSDA、MSDB，比較結(jié)果見(jiàn)表2。

表2　ISODATA提取海岸線精度統(tǒng)計(jì)

由表2可知，2003年-2014年期間，海岸線向陸一側(cè)緩沖區(qū)與向海一側(cè)緩沖區(qū)標(biāo)準(zhǔn)差比值均大于“1”，表明所提取海岸線兩側(cè)地表差異較大，因此海岸線位置較為精確。

3.3沿岸面積變化分析

采用ISODATA分類(lèi)方法所提取的研究區(qū)海岸線，經(jīng)過(guò)遙感解譯修正后所得的海岸線數(shù)據(jù)，通過(guò)研究四個(gè)主要隨時(shí)間變遷較大的區(qū)域，進(jìn)行面積變化統(tǒng)計(jì)。其中區(qū)域1位于陵水黎族自治縣，海南省東南部，屬于淤泥質(zhì)海岸。圖2為2003年-2014年間，區(qū)域1半封閉湖泊面積變化，由2003年的7.750 km2減少為2009年的7.439 km2，2014年其湖泊面積為7.184 km2。

圖2　區(qū)域1 2003年-2014年面積變化Fig．2　Region1 area changes from 2003 to 2014(a)2003年；(b)2009年；(c)2014年

圖3　區(qū)域2 2003年-2014年面積變化Fig．3　Region2 area changes from 2003 to 2014(a)2003年；(b)2009年；(c)2014年

圖4　區(qū)域3 2003年-2014年面積變化Fig．4　Region3 area changes from 2003 to 2014(a)2003年；(b)2009年；(c)2014年

圖5　區(qū)域4 2003年-2014面積變化Fig．5　Region4 area changes from 2003 to 2014(a)2003年；(b)2009年；(c)2014年

區(qū)域2位于萬(wàn)寧市。屬于人工海岸。圖3為區(qū)域2河流2003年-2014年面積變化，由2003年的5.658 km2，減少為2009年的5.152 km2，2014年其面積縮減到4.294 km2。

區(qū)域3位于萬(wàn)寧市，屬于人工海岸，周?chē)嗑用竦?。圖4為區(qū)域3半封閉湖泊的2003年-2014年面積變化。區(qū)域3的面積由43.278 km2增長(zhǎng)到為44.778 km2，2014年其面積又減少到42.497 km2。因此區(qū)域3湖泊面積2003年到2014年減少了0.781 km2。

區(qū)域4位于瓊海市，屬于人工海岸線，周邊多居民地及植被。圖5區(qū)域4的面積由2003年的6.352 km2減少為2009年的6.007 km2，2014年其面積減少為4.791 km2。

圖6　2003年-2014年NDWI提取海岸線變遷Fig．6　Changes of coastline extracted by NDWI 　　　from 2003 to 2014

3.4海岸線變遷分析

從NDWI提取海岸線圖6可以看出，研究區(qū)南?？拷Ｄ鲜|南部地區(qū)的海域近十年來(lái)，研究區(qū)海岸變遷主要在沿岸河流湖泊處及一些灣區(qū)。

采用ISODATA方法提取海岸線，同NDWI提取的海岸線相比，如圖7所示，岸線位置以及海岸線變化基本相同，但NDWI通過(guò)目視解譯提取海岸線存在無(wú)法忽視的主觀性誤差。

圖7　2003年-2014年ISODATA分類(lèi)　　　提取海岸線變遷Fig．7　Changes of coastline extracted by ISODATA 　　　from 2003 to 2014

4結(jié)論

南海是我國(guó)南部的邊緣海，為研究近年來(lái)南海海域變化，選用2003年-2009年，2009年-2014年時(shí)限的影像，通過(guò)海岸線的提取，對(duì)我國(guó)南海海域沿岸地區(qū)進(jìn)行變化監(jiān)測(cè)。采用ISODATA分類(lèi)方法自動(dòng)提取海岸線，并綜合遙感解譯知識(shí)完成海岸線提取工作，與傳統(tǒng)的主觀性目視解譯提取的海岸線相比，ISODATA方法客觀性強(qiáng)，且提取速度快所用時(shí)間少、精度較高，經(jīng)驗(yàn)證采用ISODATA方法所得分類(lèi)結(jié)果較準(zhǔn)確。

數(shù)據(jù)分析可得2003年-2014年間，變遷主要發(fā)生在沿岸湖泊以及入海河流等區(qū)域，且大部分封閉湖泊河流基本呈縮減狀態(tài)。人類(lèi)活動(dòng)的影響，自然環(huán)境的改變，全球變暖氣候變化等多種因素都造成了海岸線變遷。

基于研究所得南海海域沿岸地區(qū)十年來(lái)變化趨勢(shì)，可對(duì)南海沿岸地區(qū)未來(lái)的發(fā)展與規(guī)劃提供一種參考。

參考文獻(xiàn)：

[1]陳曉英，張杰，馬毅.近40年來(lái)海州灣海岸線時(shí)空變化分析[J].海洋科學(xué)進(jìn)展，2014，32(3):324-334.

CHEN X Y,ZHANG J,MA Y.Analysis of the spatial and temporal changes of the coastline in the Haizhou bay during the past 40 years[J].Advances In Marine Science,2014,32(3):324-334.(In Chinese)

[2]王李娟，牛錚，趙德剛.基于ETM的遙感影像的海岸線提取與驗(yàn)證研究[J].遙感技術(shù)與應(yīng)用，2010,25(2):235-239.

WANG L J,NIU Z,ZHAO D G.The study of coastline extraction and validation using ETM remote sensing image[J].Remote Sensing Technology and Application,2010,25(2):235-239. (In Chinese)

[3]王建步，張杰，陳景云，等，近30余年遼河口海岸線遙感變遷分析[J].海洋環(huán)境科學(xué).2015，34(1):86-92.

WANG J B,ZHANG J,CHEN J Y,et al.Analysis of the Liaohe Estuary coastline changes basing on the remote sensing image in the past thirty years[J].Marine Environmental Science,2015,34(1):86-92.(In Chinese)

[4]謝明鴻，張亞飛，付琨.基于種子點(diǎn)增長(zhǎng)的SAR圖像海岸線自動(dòng)提取算法[J].中國(guó)科學(xué)院研究生院學(xué)報(bào),2007，2(1):93-98.

XIE M H,ZHANG Y F,FU K.Algorithm of detection coastline from SAR images based on seeds growing[J].Journal of the Graduate School of the Chinese Academy of Sciences,2007,2(1):93-98. (In Chinese)

[5]趙英時(shí).遙感應(yīng)用分析原理與方法[M].北京：科學(xué)出版社,2003.

ZHAO Y S.The principle and method of analysis of remote sensing[M].Application.Beijing:Science Press,2003. (In Chinese)

[6]S K MCFEETERS.The use of the normalized difference water index (NDWI) in the delineation of open water features[J].International Journal of Remote Sensing,1996,17(7):1425-1432.

[7]馬小峰，趙冬至，邢小罡，等.海岸線衛(wèi)星遙感提取方法研究[J].海洋環(huán)境科學(xué),2007,26(2):185-189.

MA X F,ZHAO D Z,XING X G,et al.Means of withdrawing coastline by remote sensing[J].Marine Environmental Science,2007,26(2):185-189. (In Chinese)

[8]張霞，莊智，張旭凱，等.秦皇島市海岸線遙感提取及變化檢測(cè)[J].遙感技術(shù)與應(yīng)用，2014,29(4):625-630.

ZHANG X,ZHUANG Z,ZHANG X K,et al.Coastline extraction and change monitoring by remote sensing technology in Qinhuangdao City[J].Remote Sensing Technology And Application,2014,29(4):625-630. (In Chinese)

[9]孫偉富，馬毅，張杰，等.不同類(lèi)型海岸線遙感解譯標(biāo)志建立和提取方法研究[J].測(cè)繪通報(bào),2011(3):41-44.

SUN W F,MA Y,ZHANG J,et al.Study of remote sensing interpretation keys and extraction technique of different types of shoreline[J].Bulletin of Surveying and Mapping,2011(3):41-44.(In Chinese)

[10]孫欽幫，蘇媛媛，馬軍，等. 長(zhǎng)興島海岸線變化遙感動(dòng)態(tài)監(jiān)測(cè)及分形特征[J].海洋環(huán)境科學(xué),2011，30(3)：389-393.

SUN Q B,SU Y Y,MA J,et al. Analysis of shoreline changes based on remote sensing in Changxing Island,Liaoning Province[J].Marine Environmental Science,2011,30(3):389-393.(In Chinese)

[11]馬小峰，趙冬至，張豐收，等.海岸線衛(wèi)星遙感提取方法研究進(jìn)展[J].遙感技術(shù)與應(yīng)用,2007,22(4):575-580.

MA X F,ZHAO D Z,ZHANG F S,et al.An overview of means of withdrawing coastline by remote sensing[J].Remote Sensing Technology And Application,2007,22(4):575-580. (In Chinese)

[12]張剛，楊昕，湯國(guó)安.GIS軟件的空間分析功能比較[J].南京師范大學(xué)學(xué)報(bào)：工程技術(shù)版,2013,13(2):41-47.

ZHANG G,YANG X,TANG G A.A comparison of spatial analysis function of GIS software[J].Journal of Nanjing Normal University：Engineering and Technology Edition,2013,13(2):41-47.(In Chinese)

Change analysis and extraction of the coastline in coastal area of south China sea basing on Landsat images

WANG Shi-yang, YANG Wu-nian*, SHE Jing-xing

(Key Laboratory of Geoscience Spatial Information Technology,Ministry of Land and Resources of the P.R.China, Chengdu University of technology,Chengdu610059,China)

Abstract:The coastline is an important reference for the research of sea changes and land changes. In order to obtain the coast changes of the south China sea accurately in recent 10 years, this paper extracted the coastline automatically basing on the TM, OLI data by iterative self-organizing date analysis technique (ISODATA). The precision verification of the extraction results was carried out by zonal statistics as table. The results indicates that in recent 10 years, the area of rivers and lakes along the south China sea coast showed a rapid reduction trend.

Key words:coastline; south China sea; ISODATA; zonal statistics as table

中圖分類(lèi)號(hào)：TP 79

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

DOI：10.3969/j.issn.1001-1749.2016.01.21

文章編號(hào)：1001-1749(2016)01-0139-06

作者簡(jiǎn)介：王詩(shī)洋(1991-),女,碩士,研究方向?yàn)橘Y源與環(huán)境遙感，E-mail:sophieatsdkd@163.com。*通信作者：楊武年(1954-),男,博導(dǎo),主要從事3S技術(shù)及地學(xué)應(yīng)用科研與教學(xué)，E-mail:ywn@cdut.edu.cn。

基金項(xiàng)目:國(guó)家自然科學(xué)基金資助項(xiàng)目(41372340，41071265);高等學(xué)校博士學(xué)科點(diǎn)專(zhuān)項(xiàng)科研基金課題(20105122110006)

收稿日期：2015-05-27改回日期：2015-06-23