馬洪羽,丁賢榮,葛小平,夏娟

（1 .河海大學(xué)地球科學(xué)與工程學(xué)院,江蘇南京210098；2 .河海大學(xué)水文水資源學(xué)院,江蘇南京210098）

?

輻射沙脊群潮灘地形遙感遙測構(gòu)建

馬洪羽1,丁賢榮2,葛小平2,夏娟1

（1 .河海大學(xué)地球科學(xué)與工程學(xué)院,江蘇南京210098；2 .河海大學(xué)水文水資源學(xué)院,江蘇南京210098）

摘要:潮灘地形是合理開發(fā)利用潮間帶灘涂最重要的指標之一。灘涂水邊線測量法是目前公認的具有可操作性的灘涂地形遙感測量方法。本文以水邊線測量方法為基礎(chǔ),提出一種基于水位測站基線的潮灘地形遙感遙測構(gòu)建方法。以江蘇沿海大面積潮灘為研究區(qū),將該區(qū)域自建的4座水位測站的連線作為水位測站基線,結(jié)合不同潮情下采集的遙感影像提取的灘涂水邊線數(shù)據(jù),通過水位測站基線計算灘涂水邊線高程,從而構(gòu)建出輻射沙脊群潮灘地形。采用實測地形數(shù)據(jù)進行檢驗,結(jié)果表明:構(gòu)建的D E M能夠反應(yīng)當(dāng)前潮灘的基本形態(tài),平面精度達30 m,垂直精度達20 cm。

關(guān)鍵詞:水邊線；水位測站基線；輻射沙脊群；水邊點；水位測站基點；D E M構(gòu)建

1 引言

水邊線的含義是陸地與海水面的瞬時交界線。潮灘上的水邊線通常是隨著水位的升降而處于動態(tài)變化之中,因此具有很好的潮灘高程指示作用?；谒吘€構(gòu)建潮灘高程的方法已經(jīng)成為公認的具有可操作性、效果較好的潮灘地形遙感測量方法。該技術(shù)在測量環(huán)境以及測量范圍方面具有獨特的優(yōu)勢,對于大面積的潮灘地形構(gòu)建能夠滿足精度,當(dāng)前該技術(shù)已經(jīng)有較多的應(yīng)用。M ason等［1］利用地球資源衛(wèi)星的合成孔徑雷達數(shù)據(jù)提取了英國東岸日H u m ber/W ash區(qū)域的水邊線信息,通過復(fù)雜的潮汐-風(fēng)暴水動力模型估算水邊線高程,并據(jù)此構(gòu)建數(shù)字高程模型。該方法的水邊線高程是由水動力學(xué)模型一次性給定,但是由于缺乏實測水位、潮流資料使得模型的精度不高,則高程賦值不夠準確。國內(nèi)學(xué)者鄭宗生等［2］針對小范圍的灘涂,認定其水邊線的高程值是固定的,即認定水邊線就是等高線,然后利用附近的水位站的實測水位高程值或者斷面數(shù)據(jù),對水邊線進行賦值,生成灘涂地形。該方法針對小面積的灘涂地形可用,但是大面積的灘涂,其水邊線并不是等高線,需要更多的實測斷面測量數(shù)據(jù)和水位數(shù)據(jù)進行賦值。2008年鄭宗生等［3］通過水動力模型軟件（Delft3 D）模擬出影像成像時刻的水位值并賦值給對應(yīng)的水邊線,利用GIS空間分析功能對具有高程信息的水邊線進行空間插值,反演出潮灘的D E M。2010年朱昂［4］、2011年康彥彥［5］和李青［6］從納潮盆地地貌特征線入手,利用BP神經(jīng)網(wǎng)絡(luò)、A N U D E M插值等技術(shù)模擬細節(jié)潮灘地形。2013年,陶旭和張東［7］通過分析高光譜影像反射率與實測高程之間的相關(guān)關(guān)系,建立地表反射率與高程的轉(zhuǎn)換關(guān)系模型,實現(xiàn)潮灘大范圍概要地形信息的高光譜定量反演。但大范圍情況下,該方法反演地形要將灘涂分很多區(qū)域,而且分區(qū)沒有明確的標準,根據(jù)光譜反演因子的不同,分區(qū)也不同。有效地高程反演因子很難獲得,需要大量的重復(fù)實驗,較為耗時耗力。目前學(xué)者在該方面研究有較多成果,但并未規(guī)范化、程序化。本文在前人研究的基礎(chǔ)上,針對江蘇沿海岸灘地形變化的規(guī)律性,結(jié)合同時期的遙感影像和實測數(shù)據(jù),研究構(gòu)建大范圍岸灘地形的方法,為后續(xù)深入研究（動力環(huán)境模擬、淤泥質(zhì)潮灘沖淤演變）提供基礎(chǔ)的背景資料,同時為海岸帶土地利用、海岸工程建設(shè)、濕地環(huán)境保護及海岸帶資源的可持續(xù)利用等提供決策支持,理論及現(xiàn)實意義重大。

輻射沙脊群位于江蘇蘇北海岸外,屬于黃海南陸架海域,范圍自北部的射陽河口至長江口北部的嵩枝港［8］。南北長約200 k m,東西寬將近140 k m,在輻射沙脊群的近岸部分,多數(shù)沙脊在低潮時候出露,成為規(guī)模不等的小沙洲［9］。東海前進波與南黃海旋轉(zhuǎn)波在此相匯,沙脊群以弶港為頂點呈褶扇狀向海輻射,主體由輻射狀的沙脊以及深大潮流通道組成,像一朵巨大的蘭花綻放在江蘇沿海的黃海西岸（圖1）。在其掩蔽的江蘇海岸中部,形成了一系列寬大的潮灘,寬度一般在7～8 k m,最寬處在條子泥、高泥一帶可達40 k m。

圖1 研究區(qū)Eig .1 Research area

輻射沙脊群土地資源潛力巨大,總面積約28 000 k m2,其中出露于水面以上的沙洲與灘地面積3 782 k m2。對于江蘇這樣一個人多地少的沿海省份,輻射沙脊群可以提供巨量的后備土地資源,但海域內(nèi)水動力環(huán)境復(fù)雜,動力地貌特殊,特別是潮強灘寬的潮灘區(qū),灘地低潮出露時間短,以車船為載具的常規(guī)地形測量難以開展。這類潮灘潮差大,水流急,潮水溝擺動異常,導(dǎo)致灘地地形動態(tài)變化極為活躍。近幾十年來,未能完整地進行沿海潮灘常規(guī)地形測量,影響了沿海科學(xué)研究和開發(fā)建設(shè)的深入展開。

因此,對輻射沙脊群潮灘地形進行研究,科學(xué)地認識和評估潛在的空間資源,合理地開發(fā)利用與保護這塊稀世珍寶,將是江蘇沿海開發(fā)戰(zhàn)略性的“東方橋頭堡”。

2 數(shù)據(jù)和方法

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

主要數(shù)據(jù)包括:遙感影像數(shù)據(jù)、水位測站的水位數(shù)據(jù)和實測地形數(shù)據(jù)。

遙感影像數(shù)據(jù)主要來源于我國環(huán)境與災(zāi)害監(jiān)測預(yù)報小衛(wèi)星,共收集16景環(huán)境與災(zāi)害監(jiān)測預(yù)報小衛(wèi)星影像以及3景Landsat-T M影像,分辨率均為30 m,時間跨度為2013年12月至2014年5月（表1）。影像均經(jīng)過篩選,幅面能夠覆蓋整個輻射沙脊群海域并且成像清晰、層次豐富、反差適中,可方便進行幾何糾正和水邊線提取。

水位數(shù)據(jù)主要來源于沿海自建的4個水位測站。自北向南依次為:大豐港、東大港、新條漁港、洋口港。其中,大豐港和洋口港從2010年開始監(jiān)測水位,東大港和新條漁港從2012年10月開始監(jiān)測水位。各水位測站采用1985國家高程基面。水位測站的數(shù)據(jù)采集方式是固定時長的連續(xù)采集方式,每隔5 min采集一次,每次采集時長為3 s。由于遙感影像成像時刻不是固定的,并且不能保證與水位數(shù)據(jù)采集時刻一一對應(yīng)。因此,為了得到任意時刻水位測站的水位信息,特別是成像時刻的水位信息,需要進一步對水位數(shù)據(jù)進行插值計算。常用的獲得瞬時時刻水位數(shù)據(jù)的方法有兩種:一是根據(jù)潮汐的調(diào)和常數(shù)內(nèi)插出成像時刻的水位數(shù)據(jù)；二是根據(jù)已知的整點水位觀測數(shù)據(jù)來建立多項式,采用插值或者擬合的方法來獲取任意時刻的水位值,此種方法簡單快速,無需諸多細節(jié)。本文采用方法二,應(yīng)用M atlab工具中一維插值函數(shù)的線性插值方法計算出影像獲取時刻的水位值,重復(fù)上述計算得到19景影像獲取時刻對應(yīng)的水位數(shù)據(jù)（表1）,從表中可以看出各水位測站間的水位同步不同高。

表1 遙感影像信息及成像時刻各水位測站水位值（m）Tab .1 RS images information and each hydrological station water level of corresponding time（m）

實測數(shù)據(jù)用于對構(gòu)建的輻射沙脊群潮灘D E M驗證,以評價D E M的質(zhì)量以及構(gòu)建方法的優(yōu)劣。包括2013年7月實測斷面、沿海集團提供的2012年1月條子泥測點和測繪部門提供的2006年Lidar地形。

2.2 研究方法

2.2.1 水邊線提取方法

圖2 4種水邊線提取方法Eig .2 Eour extraction methods of water boundary

本文對比4種水邊線提取方法,單波段閾值分割法（圖2a）、邊緣檢測算法（圖2b）、多波段譜間關(guān)系法（圖2c）和面向?qū)ο罂臻g特征提?。▓D2d）。前三者的核心是閾值確定,閾值的選取直接決定水體提取的精確性。當(dāng)水體和影像中其他地物的反射率區(qū)分困難時,很難找到一個準確的閾值將水體與其他地物完全的區(qū)分開,從而造成水體提取質(zhì)量較差。因此,通過確定閾值來提取水體信息的方法對地物類型比較豐富、地物在所選波段上的灰度值相近的影像是有一定局限性。應(yīng)用E N VI軟件的面向?qū)ο罂臻g特征提取模塊——Eeature Extraction（簡稱E N VI E X）提取水邊線。該技術(shù)主要是集合臨近像元為對象用來識別感興趣的光譜要素,充分利用高分辨率的全色和多光譜數(shù)據(jù)的空間、紋理以及光譜信息來分割和分類的特點,以高精度的分類結(jié)果或者矢量輸出［10］。此法基于邊緣分割算法,計算快且僅需一個輸入?yún)?shù)就能產(chǎn)生多尺度的分割效果,在操作過程中隨時預(yù)覽影像分割效果,適用于中高分辨率多光譜和全色影像。最終采用較為符合本文所用地形構(gòu)建方法的面向?qū)ο罂臻g特征提取方法提取21景影像的水邊線。

2.2.2 基于水位測站基線的水邊點高程賦值方法

由于輻射沙脊群的空間范圍跨度較大,衛(wèi)星過境時刻的海面不可能被視為一個平面,所以同一條水邊線上的不同位置的高程也是不同的,從遙感圖像中提取出的水邊線也無法直接用來構(gòu)建輻射沙脊群潮灘地形。因此,將水邊線按照影像的分辨率即30 m間隔進行離散,得到水邊點數(shù)據(jù)。

水邊點的高程賦值是地形構(gòu)建的關(guān)鍵環(huán)節(jié),本文采用一種新的賦值方法,通過水位測站基線計算水邊點高程?；€是指在測量、計算或定位中的一種基本參照。水位測站基線是4個水位測站之間的連線,大豐港-東大港,東大港-新條漁港,新條漁港-洋口港（圖3）,以下簡稱基線。該基線建立的原理是,在水動力環(huán)境相同水域可建立多測站間的同步水位基線,即文中的水位測站基線,水動力環(huán)境不同,則不可建立。如輻射沙脊群南北兩側(cè),受條子泥、高泥、東沙等高灘阻隔,低潮階段西洋與黃沙洋、爛沙洋不可建立同步水位基線,而高潮階段沙頂淹沒,則可建立。對于多測站海域,任意兩測站之間可形成一條同步水位基線。該基線兩端水位測站的水位具有同步不同高的性質(zhì),基線上任意一點的同步水位可按照相鄰兩站距離權(quán)重進行插值獲得。

潮位站站網(wǎng)基線是折線,水邊線是與灘涂地形無縫吻合的曲線,二者的平面空間可重合也可不重合。根據(jù)遙感影像空間分辨率大小,水位測站基線同樣按照30 m分辨率離散,生成水位測站基點,以下簡稱基點。根據(jù)海面連續(xù)性特性,在潮位可比海域,按照最近距離水位傳遞原則,水邊點水位分別取自與其最鄰近基點的水位（圖3）。以此類推,可進行系列水邊線的高程賦值。

基點的高程值由水位測站的同步水位插值獲得。本文應(yīng)用M AT L A B工具,對比了4種插值方法（圖4）:線性插值、三次樣條插值、三次多項式插值和最近鄰點插值（linear、spline、cubic、nearest）其調(diào)用格式為:

圖3 水位測站基線圖Eig .3 H ydrological station baseline

式中,d代表的是已知數(shù)據(jù)點的橫坐標向量,即由北至南水位測站間的距離,d必須單調(diào)；h代表已知數(shù)據(jù)點的縱坐標向量,即對應(yīng)水位測站的同步水位；D為未知數(shù)據(jù)點的橫坐標數(shù)據(jù)（或數(shù)據(jù)組）,D不能超出d的范圍,而method為可選插值方法［11］。從插值結(jié)果來看,線性插值的正確性最高,數(shù)值無外拋或者內(nèi)凹現(xiàn)象；又根據(jù)水面連續(xù)性和水位“跳繩效應(yīng)”特征（圖5）,選擇線性插值最為合理。因此,采用線性插值方法計算出21個時刻對應(yīng)的水位測站基點的高程值。

水位測站基點和水邊點均是常規(guī)三維空間點,采用最鄰近算法將基點的高程值賦給水邊點。最鄰近算法是指計算基點和水邊點之間互為最近距離的點, 圖3中水位傳遞示意線是該算法的基本表達,依據(jù)最近距離將基點的高程值賦給對應(yīng)水邊點。二者在海面上東西向分布,其之間最大距離為9 k m,采用最鄰近算法賦值在精度允許范圍內(nèi)。此外,海水面的連續(xù)性特征,使得空間上越靠近,屬性越相似,相關(guān)性越高,即兩點之間距離越小,水位越相近。從而得到具有高程信息的水邊點。

圖4 計算水位測站基點高程的4種插值方法Eig .4 Eour interpolation methods of hydrological station baseline points elevation calculating

圖5 水位的跳繩效應(yīng)Eig .5 Rope skipping effect of water level

3 成果分析

3.1 潮灘DE M構(gòu)建

插值是利用有限數(shù)目的樣本點來估計未知樣本點值的一種方法,這種估值方法可用于生成高程、降雨量、化學(xué)污染程度、噪聲等級、湖泊水質(zhì)等級等連續(xù)表面［12］。在獲取了潮灘水邊線以及完成給水邊點賦值后,進一步探討構(gòu)建輻射沙脊群潮灘地形。

插值的前提是空間地物之間具有一定的空間相關(guān)性,距離越近的地物相關(guān)性越高,其值也就越接近。由點數(shù)據(jù)插值生成柵格的方法有很多種,本文討論了自然鄰近插值法（Natural Neighbor）、反距離權(quán)重插值法（Inverse Distance W eighting,ID W）、克里格插值以及規(guī)則樣條函數(shù)插值4種常用方法。以泥螺珩為實驗區(qū),將已賦值的水邊點的高程數(shù)據(jù)作為輸入,內(nèi)插生成潮灘地形。與2006年實測Lidar數(shù)據(jù)作差,獲取誤差分布圖（圖6）,結(jié)果表明Natural Neighbor生成的地形較為真實,精度最好。

此外,從4種方法的誤差指標對比來看（表2）:除樣條函數(shù)插值法,其他插值方法的誤差平均值都很小。雖然ID W法和Natural Neighbors法兩者的各精度指標相近,但從D E M的表面光滑性來看,Natural Neighbors方法更優(yōu)。因此,本文采用Natural Neighbors方法構(gòu)建D E M ,通過調(diào)用ArcGIS中3 D Analyst模塊的Natural Neighbor Interpolation命令,輸入21個時刻水邊點的高程數(shù)據(jù),cell size設(shè)置為30 m×30m,以保證整體數(shù)據(jù)分辨率的一致性,插值構(gòu)建出輻射沙脊群潮灘D E M（圖7）。

圖6 誤差直方圖Eig .6 Error histogram

表2 4種地形插值方法的誤差指標表（單位:m）Tab .2 Four error indicators of terrain interpolation methods（Unit :Meter）

3.2 精度驗證

2012年條子泥測點數(shù)據(jù)是由沿海集團提供, 1∶1萬測圖比例,廢黃河基面,測區(qū)面積128.01 k m2,主要分布在西大港、東大港區(qū)域的灘涂（圖8）。境,潮水溝變化速度快也是實地測量的重要困難之一。

圖7 輻射沙脊群潮灘地形反演結(jié)果Eig .7 Terrain of radial sand ridges tidalflat

由2012年測點內(nèi)插生成DE M,與同期影像疊加對比,二者吻合,證明數(shù)據(jù)的可靠性。2014年DE M與2012年DE M作差值運算,獲得該區(qū)域沖淤空間分布和地形誤差統(tǒng)計（圖9,10）。其平均值- 0.28 m,最大值4.53 m,最小值- 6.27 m,誤差小于- 2 m的數(shù)據(jù)較多,使誤差分布不服從正態(tài)分布。因此,應(yīng)用A RC M AP中R ASTER C A LC U L ATE工具計算沖刷高度在- 2 m以下區(qū)域（圖11）,可以看出西大港和東大港的潮水溝擺動造成灘涂沖刷劇烈,尤其是東大港潮水溝。對該區(qū)域沖淤進行定量統(tǒng)計（表3）,其中沖刷是差值小于0區(qū)域,淤積是差值大于0區(qū)域。潮水溝的擺動是造成沖刷和淤積量大的主要原因,也符合該地區(qū)水動力環(huán)

圖8 2012年1月條子泥1∶1萬實測點圖Eig .8 January 2012 Tiaozini 1∶10 000 measuring points

圖9 條子泥灘涂2014 - 2012年沖淤空間分布Eig .9 The spatial distribution of erosion of 2014 - 2012 in Tiaozinitidalflat

圖10 地形誤差統(tǒng)計圖Eig .10 Statisticalfigure of terrain error

圖11 沖刷高度在- 2 m以下區(qū)域Eig .11 Scour altitude below - 2 m in the area

表3 2012 - 2014年沖淤量統(tǒng)計Tab .3 Quantity of statistics of scouring and siltingin 2012 - 2014

需要指出的是,水位測站基線法構(gòu)建的輻射沙脊群D E M與2012年D E M的總平均差實際為0.08 m。因為,前者采用1985國家高程基準,后者則是廢黃河基面,二者在江蘇中部輻射沙脊群海岸相差0.2 m,即廢黃河基面低于黃海85基面0.2 m。因此,水位測站基線法構(gòu)建的輻射沙脊群D E M減去2012年D E M的總平均差應(yīng)該是0.2 m - 0.28 m = - 0.08 m。由此說明,灘涂體積總體穩(wěn)定,潮水溝地區(qū)沖淤變化較大,2014年D E M與所測量的局部地區(qū)相比略有沖刷趨勢。

圖12 2014 - 2006年輻射沙脊群地形沖淤分布Eig .12 The spatial distnbution of erosion of 2014 - 2006

2006年Lidar地形數(shù)據(jù)由測繪部門提供,1∶1萬測圖比例,分辨率是5 m×5 m,黃海1985基面,有效測區(qū)面積286.399 km2。由于2006年實測數(shù)據(jù)與2014年地形露灘的形態(tài)不同,因此,通過空間掩膜運算提取二者有效值的區(qū)域,計算差值,獲得2014 - 2006年地形沖淤分布及沖淤統(tǒng)計（圖12,14）。沖淤統(tǒng)計結(jié)果:平均值2.16 m,最大值12.27 m,最小值- 6.71 m,在5～10 m間出現(xiàn)異常誤差值。因此,提取淤積高度在5～10 m區(qū)域（圖13）,是較2006年灘涂明顯增長區(qū)域,由圖13可以看出,東大港和西大港潮水溝的擺動仍舊是造成沖刷的主要地帶,東大港潮水溝明顯向南擺動,條子泥灘涂南部淤積大于北部,局部小區(qū)域略有沖刷,對于工程實施有利。總體上從2006 - 2014年平均淤積2.16 m,平均年淤積率0.27 m。

圖13 沖刷高度在5 m以上區(qū)域Eig .13 Scour altitude above 5 m in the area

圖14 2014 - 2006年地形誤差沖淤統(tǒng)計圖Eig .14 Statisticalfigure of erosion of 2014 - 2006

為檢驗本文所構(gòu)建地形,只考慮兩D E M公共區(qū)域,獲取誤差統(tǒng)計（圖15）,平均值0.03 m,最大值11.62 m,最小值- 6.71 m,標準差1.46 m。由此說明,公共區(qū)域灘涂穩(wěn)定性較高,地形基本穩(wěn)定。同時,二者相減的誤差平均值為3 cm,符合地形測量精度。對該區(qū)域沖淤進行定量統(tǒng)計（表4）,對比表3可知,局部沖淤變化很大的情況下,隨著范圍的擴大,沖淤被2013年條子泥實測斷面數(shù)據(jù),黃海1985高程基面,共12條測量斷面,其中3條位于本文地形構(gòu)建區(qū)域,每隔3.5 m一個測點,共12 058個測點（圖18）。平均化,但總體為淤長型灘涂。并且沖淤變化較大的原因是水動力引起的潮水溝的擺動。

表4 2006 - 2014年沖淤量統(tǒng)計Tab .4 Quantity of statistics of scouring and silting during 2006 to 2014

圖15 2006 - 2014年地形誤差統(tǒng)計圖Eig .15 Statisticalfigure of terrain error of 2014 - 2006

應(yīng)用A R C M A P空間分析模塊的E X T R A C T VA L U ES T O P OIN TS工具,提取與實測斷面點位置對應(yīng)的輻射沙脊群D E M高程值,計算二者高程差,得到誤差統(tǒng)計（圖17）,服從正態(tài)分布。圖18是兩斷面高程值的折線圖,可以看出灘涂高程值整體變化趨勢整體一致,局部潮水溝地區(qū)變化較大。一方面原因是水動力造成潮水溝的擺動、沖淤；另一方面是本文構(gòu)建的D E M分辨率是30 m,而實測斷面分辨率是4 m,構(gòu)建的地形相對較粗略,點序1420 - 1807之間就是此情況引起的誤差較大。

應(yīng)用select by attributes工具提取誤差值大于1 m的點,圖中藍色高亮點,可以看出同樣是由于潮水溝的沖淤造成。尤其是條魚港區(qū)域的淤積十分明顯,對應(yīng)表中點序7354 - 8902。

圖16 2013年數(shù)據(jù)測點信息Eig .16 The measured points information in 2013

圖17 2014 - 2013年誤差統(tǒng)計圖Eig .17 Statisticalfigure of error of 2014 - 2013

圖18 兩斷面高程值的折線圖Eig .18 Two section elevation value of line chart

4 結(jié)論

以水邊線法為基礎(chǔ),結(jié)合水位測站基線賦值方法,構(gòu)建輻射沙脊群潮灘地形的模擬方法在理論上可行,精度可信,結(jié)果符合潮灘地形基本規(guī)律。該方法實現(xiàn)水邊線從2 D到3 D的轉(zhuǎn)換,除大比例尺工程測量的精度要求外,基本滿足對潮灘演變、沖淤計算等基礎(chǔ)研究的地形精度要求,可作為主要數(shù)據(jù)資料使用。D E M結(jié)果表明該方法是建立在遙感影像和水位數(shù)據(jù)相結(jié)合現(xiàn)代信息技術(shù)基礎(chǔ)上的新型測量方法,對于灘涂濕地動態(tài)多變地形具有快速、高效的測量效果,此方法具有可操作性,制圖成果具有實用性。但根據(jù)現(xiàn)有影像很難獲取最低潮影像；影像分辨率30 m,無法滿足大比例尺精度；該賦值方法存在一定的缺陷,由于受輻射沙脊群沙脊的影響,使該區(qū)域不能視為開暢水域,導(dǎo)致亮月沙、東沙東部灘涂賦值不夠準確,但從驗證結(jié)果來看,并未導(dǎo)致結(jié)果不正確,說明只需進一步優(yōu)化,達到更高精度。這些將是本文進一步研究的方向。該方法的優(yōu)點是大范圍灘涂地形的構(gòu)建可以程序化、客觀化,真正意義上達到快速高效。

參考文獻:

［1］ M ason D C,Davenport I J,Elather R A,et al. Cover A digital elevation model oftheinter-tidal areas ofthe W ash,England,produced by the waterline method［J］.International Journal of Remote Sensing,1998,19（8）:1455 - 1460 .

［2］ 鄭宗生,周云軒,蔣雪中,等.崇明東灘水邊線信息提取與潮灘D E M的建立［J］.遙感技術(shù)與應(yīng)用,2007,22（1）:35 - 38,94 . Zheng Zongsheng,Zhou Yunxuan,Jiang Xuezhong,et al. W aterline extraction and D E M reconstruction in Chongming Dongtan［J］. Remote Sensing Technology and Application,2007,22（1）:35 - 38,94 .

［3］ 鄭宗生,周云軒,劉志國,等.基于水動力模型及遙感水邊線方法的潮灘高程反演［J］.長江流域資源與環(huán)境,2008,17（5）:756 - 760 . Zheng Zongsheng,Zhou Yunxuan,Liu Zhiguo,et al. D E M reconstruction based on hydrodynamicmodel and waterline method［J］. Resources and Environ mentin the Yangtze Basin,2008,17（5）:756 - 760 .

［4］ 朱昂.潮水溝LiD A R地形神經(jīng)網(wǎng)絡(luò)插補方法［D］.南京:河海大學(xué),2010:73 . Zhu Ang . A method ofinterpolation for the Voidsin Tidal Channel's D E M obtained by LiD A R base on neural network［D］. Nanjing:H ohai U niversity,2010:73 .

［5］ 康彥彥.基于納潮盆地地貌特征線的潮灘地形模擬研究［D］.南京:河海大學(xué),2011:69 . Kang Yanyan .Simulation of tidalflat terrain based on geomorphologic feature line of tidal basin［D］.Nanjing:H ohai U niversity,2011:69 .

［6］ 李青.潮灘數(shù)字地形神經(jīng)網(wǎng)絡(luò)建模研究［D］.南京:河海大學(xué),2011:78 . Li Qing . Neural network modeling research of tidalflat digital terrain［D］. Nanjing:H ohai U niversity,2011:78 .

［7］ 王穎.黃海陸架輻射沙脊群［M］.北京:中國環(huán)境科學(xué)出版社,2002:29 - 228 . W ang Ying . The Radiated Sand Ridge of Yellow Sea Shelf［M］. Beijing:China Environ mental Science Press,2002:29 - 228 .

［8］ 任美鍔.江蘇省海岸帶和海涂資源綜合調(diào)查報告［M］.北京:海洋出版社,1985:60 - 236 . Ren M ei'e .Jiangsu Coastal Zones and Tideland Resources Comprehensive Survey Report［M］. Beijing:China Ocean Press,1985:60 - 236 .

［9］ 鄧書斌. E N VI遙感圖像處理方法［M］.北京:科學(xué)出版社,2010:212 - 220 . Deng Shubin . E N VI Remote Sensing Image Processing M ethods［M］. Beijing:Science Press,2010:212 - 220 .

［10］ 趙靜.數(shù)學(xué)建模與數(shù)學(xué)實驗［M］.第三版.北京:高等教育出版社,2008:306 - 307 . Zhao Jing . M athematical M odeling and M athematical Experiment［M］.3rd edition . Beijing:Higher Education Press,2008:306 - 307 .

［11］ 湯國安,楊昕. A R C GIS地理信息系統(tǒng)空間分析實驗教程［M］.北京:科學(xué)出版社,2006:309 . Tang Guo'an,Yang Xin . A R C GIS Geographic Information System Spatial Analysis Experiment Tutorial［M］. Beijing:Science Press,2006:309 .

［12］ Lohani B .Construction of a digital elevation model ofthe H olderness Coast using the waterline method and airborne thematic mapper data［J］.International Journal of Remote Sensing,1999,20（3）:593 - 607 .

［13］ M ason D C,Davenport I J,Elather R A,et al. A sensitivity analysis of the waterline method of constructing a digital elevation modelforintertidal areas in E RS SA R scene of eastern England［J］. Estuarine,Coastal and Shelf Science,2001,53（6）:759 - 778 .

［14］ Lohani B M D . Application of airborne scanning laser altimetry to the study oftidalchannel geomorphology［J］.ISPRS Journal of Photogram metry and Remote Sensing,2001,56（2）:100 - 120 .

［15］ Ryu Joo-H yung,W on Joong-Sun,Min Kyung Duck . W aterline extraction from Landsat T M data in a tidalflat A case study in Gomso Bay,Korea ［J］. Romote Sensing of Environ ment,2002,83:442 - 456 .

［16］ 韓震,惲才興,戴志軍,等.淤泥質(zhì)潮灘高程及沖淤變化遙感定量反演方法研究——以長江口崇明東灘為例［J］.海洋湖沼通報,2009（1）:12 - 18 . Han Zhen,Yun Caixing,Dai Zhijun,et al. Remote sensing quantitativeinversion research ofthe silttidalflat elevation and deposition and erosion—Chongming Dongtan as the example［J］. Transactions of Oceanology and Limnology,2009（1）:12 - 18 .

［17］ 張長寬.江蘇省近海海洋環(huán)境資源基本現(xiàn)狀［M］.北京:海洋出版社,2014:101 - 102 . Zhang Changkuan . Basic Status Quo of Offshore M arine Environ ment and Resourcesin Jiangsu Province［M］. Beijing:China Ocean Press,2014: 101 - 102 .

［18］ 鄭宗生,周云軒,沈芳,等.基于D T M的水邊線遙感信息提取方法［J］.國土資源遙感,2007（2）:56 - 59 . Zheng Zongsheng,Zhou Yunxuan,Shen Eang,et al. W aterline extraction from remotely sensedimages with D T M［J］. Remote Sensing for Land & Resource,2007（2）:56 - 59 .

［19］ 韓震,郭永飛,李睿,等.長江口淤泥質(zhì)潮灘環(huán)形水邊線信息提取方法研究［J］.國土資源遙感,2010（4）:64 - 66 . Han Zhen,Guo Yongfei,Li Rui,et al. Research on the method for ring watersidelineinformation extraction from mudflatin the Yangtze River estuary［J］. Remote Sensing for Land & Resource,2010（4）:64 - 66 .

［20］ Chen L C . A unified solution for digitalterrain modeland orthoimage generation from SP O T stereopairs［J］.IEEE Transactions on Geoscience and Remote Sensing,1993,31（6）:1243 - 1252 .

［21］ 張鷹,丁賢榮,王文.水深遙感與潮灘地形沖淤變化分析［J］.港口工程,1998（2）:26 - 30 . Zhang Ying,Ding Xianrong,W ang W en . Analysis remote sensing to water depth and change in scour and fillin tidal flat［J］. Port Engineering, 1998（2）:26 - 30 .

［22］ 錢一婧,張鷹,李洪靈,等.常用邊緣檢測算法在遙感影像水邊線提取比較［J］.人民長江,2008,39（13）:95 - 97 . Qian Yiqian,Zhang Ying,Li H ongling,et al. Com monly used image edge detection algorith m in remote sensing to extract water edge［J］. Yangtze River,2008,39（13）:95 - 97 .

［23］ 劉永學(xué),李滿春,張忍順.江蘇輻射沙洲水邊線自動提取方法研究［J］.海洋科學(xué),2004,28（6）:42 - 45 . Liu Yongxue,Li M anchun,Zhang Renshun .Jiangsu radiation ahoal water edge automatic extraction method research［J］. M arine Science,2004,28 （6）:42 - 45 .

中圖分類號:P737.2

文獻標志碼:A

文章編號:0253-4193（2016）03-0111-12

收稿日期:2015-02-12；

修訂日期:2015-05-14。

基金項目:國家科技支撐計劃（2012B A B03B01）；國家海洋公益性行業(yè)科研專項（201005006）；江蘇省基礎(chǔ)研究計劃（自然科學(xué)基金）（B K2012414）。

作者簡介:馬洪羽（1990—）,女,吉林省遼源市人,主要從事海洋遙感與地理信息系統(tǒng)研究。E-mail:hhu_m hy @ 163 .com

馬洪羽,丁賢榮,葛小平,等.輻射沙脊群潮灘地形遙感遙測構(gòu)建［J］.海洋學(xué)報,2016,38（3）:111 - 122,doi:10.3969/j.issn .0253-4193.2016.03.011

M a H ongyu,Ding Xianrong,Ge Xiaoping,et al. Remote sensing and remote measuring approach to construct tidalflat terrain in the radial sand ridges［J］. Haiyang Xuebao,2016,38（3）:111 - 122,doi:10 .3969/j.issn .0253-4193.2016.03.011

Remote sensing and remote measuring approach to construct tidal flatterrain in the radial sand ridges

M a H ongyu1,Ding Xianrong2,Ge Xiaoping2,Xia Juan3
（1 . HohaiUniversity Schoolof Earth Scienceand Engineering,Nanjing 210098,China；2 . HohaiUniversity Collegeof Hydrology and Water Resources,Nanjing 210098,China；3 . Hohai University Schoolof Earth Science and Engineering,Nanjing 210098,China）

Abstract:Tidalflat terrain is one of the mostimportantindexes for us to exploitintertidal mudflatlegitimately .Its measure method that using water boundary of tidal flat is now widely approved as a feasible method by remote sensing . Based on this method,the paper provided a new way named hydrological station baseline to constructtidal flat terrain . The way is established in a large tidalflat of radial sand ridges in coastal area of Jiangsu Province byconnecting 4 hydrological stationsin a row,used the values of waterlevel,getting from the hydrological stations,to calculate water boundaries' elevation which were extracted from remote sensing images with multi tide situations . Accordingly,completed the construction of the tidalflat terrain of radial sand ridges . The accuracy was confirmed by actual measurement terrain . Results show that the D E M can reflect the current basic form of tidalflat and the position accuracy approach 30 m and the elevation accuracy approach 20 cm .

Key words:water boundary；hydrological station baseline；radial sand ridges；water boundary points；hydrological station baseline points；D E M construction