黨亞民, 章傳銀，周興華, 薛樹強(qiáng)，許 軍

(1. 中國測繪科學(xué)研究院，北京 100830； 2. 海軍大連艦艇學(xué)院，遼寧 大連 116018； 3. 國家海洋局第一海洋研究所，山東 青島 266061)

海島岸線遙感立體測圖精細(xì)測量方法

黨亞民1, 章傳銀1，周興華3, 薛樹強(qiáng)1，許 軍2

(1. 中國測繪科學(xué)研究院，北京 100830； 2. 海軍大連艦艇學(xué)院，遼寧 大連 116018； 3. 國家海洋局第一海洋研究所，山東 青島 266061)

提出了一種海島岸線遙感測圖精細(xì)測量新方法，該方法直接基于理論定義的海岸線，利用航空影像瞬時(shí)水涯線數(shù)據(jù)在立體測圖環(huán)境中提取瞬時(shí)水位高程；利用海島周邊精密海潮模型和瞬時(shí)水位高程推算海島岸線高程；最后依據(jù)海島岸線高程，采用立體測圖方法測制海島岸線的平面位置。該方法確保了海島岸線成果的唯一性和連續(xù)性，適合大比例尺的大陸海岸線和海島岸線測量。測試結(jié)果顯示，在較高精度海潮模型和海面地形支持下，海島岸線高程精度優(yōu)于0.2 m，可滿足1∶2000測圖要求。

海島岸線；立體測圖；瞬時(shí)水位；海洋潮汐

海島岸線定義為平均大潮高潮時(shí)刻的海陸分界線[1]，即平均大潮高潮線矢量。海島岸線決定了海島的形狀和大小，是海島基礎(chǔ)地理信息中最為重要的要素之一。隨著測圖比例尺的不同，海島岸線的長度和形狀都存在較大差異，具有典型的非線性分形特性。許多情況下，海島岸線沒有明顯的地物特征或痕跡，不能通過實(shí)地定位和遙感測圖方法直接獲取[2]。目前國內(nèi)外海島岸線測量方法主要有痕跡岸線現(xiàn)場測量、遙感影像痕跡岸線判繪和影像水邊線提取3種[3]。受自然岸線的模糊性和多義性等影響，岸線測量成果存在不一致的現(xiàn)象，很難滿足大比例尺海島測圖要求。海島岸線測繪一直以來是一個(gè)國際性技術(shù)難題。

海島岸線具體實(shí)現(xiàn)和確定方法有所不同[4]。我國將平均大潮高潮時(shí)留下的痕跡線作為海島岸線[5]，實(shí)質(zhì)上以測定海水及漂浮物在海岸和島礁上留下的較重痕跡作為岸線的測定手段[3]。一般情況下，陸、島海岸上留下的痕跡線不唯一，也不明顯，有一定寬度，岸線測量中往往依賴于測量人員的工作經(jīng)驗(yàn)。此外，海島痕跡岸線的高程普遍不等，有的相差數(shù)米，可見，海島痕跡岸線作為海島岸線有一定的局限性。

本文采用嚴(yán)格定義的海島岸線代替海島痕跡岸線，測量平均大潮高潮時(shí)海島的海陸分界線作為海島岸線，在充分綜合遙感信息、潮汐特征與海面地形信息的基礎(chǔ)上，采用與陸地一致的海島岸線測繪成果的坐標(biāo)系統(tǒng)和高程基準(zhǔn)，同步解決海島岸線平面位置與海島岸線高程測定的技術(shù)問題。該方法具有快速、高效、高程測定精度高、便于大規(guī)模測繪作業(yè)等特點(diǎn)，解決了基于遙感測圖的海島岸線測繪技術(shù)難題。用于替代傳統(tǒng)的野外登島測繪方法，作業(yè)成本大幅降低，作業(yè)周期明顯縮短，且可避免由于痕跡岸線的模糊性導(dǎo)致的傳統(tǒng)岸線測繪成果難以集成的現(xiàn)象發(fā)生，能有效保證海島岸線測繪成果的唯一性和連續(xù)性。

1 立體測圖環(huán)境中海島岸線測量方法

立體測圖環(huán)境中，影像水邊線的平面和高程已知，其高程精度也與相應(yīng)比例尺地物點(diǎn)的高程精度相當(dāng)。假設(shè)影像水邊線為一近似等高線[6]，若能采用潮汐和水位推算方法求得水邊線與海島岸線的高差，就可由水邊線高程得到海島岸線的高程。這樣，在立體測圖環(huán)境中，可依據(jù)海島岸線高程，采樣跟蹤出海島岸線的平面坐標(biāo)[7]。立體測圖環(huán)境中海島岸線跟蹤技術(shù)流程如圖1所示。

圖1 立體測圖環(huán)境中海島岸線跟蹤技術(shù)流程

1.1 影像水邊線的選擇與高程提取

影像水邊線高程將用于計(jì)算海島岸線高程的起算基準(zhǔn)，因此要求所選擇的水邊線高程不符值盡量小些，即滿足水邊線為近似等高線的條件。

影像水邊線高程及其不符值一般具有如下性質(zhì)：

(1) 曝光時(shí)刻海況越好，水邊線不符值越小。

(2) 附近有浪花的影像水邊線，其高程不符值相對(duì)較大。

(3) 海島背風(fēng)面的影像水邊線，其高程不符值相對(duì)較小。

(4) 對(duì)于有風(fēng)浪影響的影像，在迎風(fēng)面的開闊地帶、岸線平滑地方或能阻擋風(fēng)浪的海灣處，其水邊線的高程不符值相對(duì)較小。

(5) 海岸地形坡度越小，立體測圖環(huán)境中水邊線高程提取的準(zhǔn)確度越高。

(6) 當(dāng)海況條件較好時(shí)，曝光時(shí)間間隔短(如小于5 min)的不同影像，在一定空間范圍內(nèi)(如小于10 km)，其水邊線的高程值應(yīng)比較接近。

1.2 海島岸線高程推算

當(dāng)水邊線為近似等高線時(shí)，利用提取后的水邊線高程計(jì)算水邊線平均高程。

利用潮汐模型計(jì)算平均大潮高潮位，依據(jù)影像曝光時(shí)刻，結(jié)合周邊水位觀測數(shù)據(jù)推算水邊線瞬時(shí)水位；將平均大潮高潮位與水邊線瞬時(shí)水位相減，求得海島岸線與水邊線的高差；由該高差值加上水邊線平均高程，求得海島岸線的高程。

平均大潮高潮位精度僅受潮汐模型精度影響。但是水邊線瞬時(shí)水位精度不僅受潮汐模型影響，還受余水位和海況條件影響[8]。因此，海島岸線與水邊線高差計(jì)算時(shí)，最好有曝光時(shí)刻前后該海區(qū)的同步水位觀測數(shù)據(jù)支持。

1.3 立體測圖環(huán)境中海島岸線采樣跟蹤

在立體測圖環(huán)境中，依據(jù)計(jì)算的海島岸線高程，按類似等高線跟蹤方法采樣跟蹤出海島岸線。

在海島岸線采樣跟蹤過程中，應(yīng)參考痕跡岸線、水邊線的空間形態(tài)及附近地貌特征，最大限度地保持海島岸線的逼真信息。當(dāng)跟蹤出的海島岸線與實(shí)際地形要素發(fā)生矛盾時(shí)，可在一定的高程限差范圍內(nèi)，調(diào)整岸線平面位置，消除矛盾。

通常情況下，海岸地形坡度越大，立體環(huán)境下采樣跟蹤海島岸線的平面位置越準(zhǔn)確。因此，實(shí)際作業(yè)過程中，當(dāng)遇到地形坡度很小或海岸地形復(fù)雜的地方，應(yīng)仔細(xì)采樣跟蹤，以提高岸線的準(zhǔn)確性。

海島平均水位線采樣跟蹤也可以參考上述流程實(shí)施。

2 影像水邊線高程與海面瞬時(shí)水位推算

遙感影像在拍攝瞬間，海面與海島相截處會(huì)在影像上顯示出一條截線，這條截線稱為影像水邊線，水邊線表征了海面瞬時(shí)水位與海島的空間關(guān)系。

海洋中涉及的水位是指海面整體的垂直升降，海風(fēng)等引起的涌浪不計(jì)入水位垂直變化中，在水位觀測中或水位預(yù)處理時(shí)應(yīng)盡量消除或減弱波浪的影響。瞬時(shí)水位主要分為兩個(gè)部分：主體是引潮力在特定海底地形和海岸形狀等因素制約下引起的海面升降，通常稱為天文潮位；氣壓、風(fēng)等氣候、氣象作用引起的水位變化，其中周期性部分以氣象分潮形式歸入天文潮位，而剩余的短期無周期性部分，通常稱為余水位或異常水位，其主要激發(fā)因素是短期氣象變化，在一定范圍內(nèi)具有較強(qiáng)的空間相關(guān)性。因此，瞬時(shí)水位變化h(t)可表示為[8]

h(t)=MSL+T(t)+R(t)+Δ(t)

(1)

式中，MSL為當(dāng)?shù)亻L期平均海面在驗(yàn)潮零點(diǎn)上的垂直高度；T(t)為天文潮位；Δ(t)為觀測誤差；R(t)為余水位。

天文潮位T(t)是諸多分潮對(duì)引潮力響應(yīng)的疊加，當(dāng)?shù)亻L期平均海面可看作是其平衡位置，相對(duì)于平衡位置的天文潮位可表示為

(2)

式中，m為分潮的個(gè)數(shù)；Hi、gi分別稱為分潮i的振幅和遲角，構(gòu)成分潮i的調(diào)和常數(shù)，在潮汐調(diào)和分析時(shí)為未知的待求量，而用于潮汐預(yù)報(bào)時(shí)作為已知量；vi(t)為分潮的天文幅角。

振幅H和遲角g反映了海洋對(duì)一頻率外力的響應(yīng)，可近似為常數(shù)。振幅H和遲角g稱為實(shí)際潮汐分潮的調(diào)和常數(shù)。

遙感影像立體測圖時(shí)，可測出影像水邊線的平面位置矢量和高程值，因此，可以得到影像拍攝時(shí)刻海面瞬時(shí)水位的高程。顯然，水位高程的高程基準(zhǔn)與測圖高程基準(zhǔn)一致。

3 平均大潮高潮面與海島岸線高程推算

平均大潮高潮面與海島地面相截，其截線就是理論定義的海島岸線。平均大潮高潮面與海島岸線高程推算步驟如下：

3.1 平均大潮高潮面的選擇

平均大潮高潮面是一種特征潮位面，用其相對(duì)于當(dāng)?shù)仄骄Ｃ娴母叨葋肀硎?，該高度值稱為凈空基準(zhǔn)值Q，其大小與當(dāng)?shù)爻毕膹?qiáng)弱即潮差的大小有著密切的關(guān)系。在海洋潮汐學(xué)中，平均大潮高潮面是指半日潮大潮期間高潮位的平均值，即定義限制于半日潮類型海域，這與大潮概念只存在于半日潮類型有關(guān)。在日潮類型海域，回歸潮與半日潮類型的大潮具有相似的極值意義，但日潮類型海域回歸潮時(shí)的日潮不等現(xiàn)象十分明顯，故在日潮類型海域可采用平均回歸潮高高潮面。對(duì)于混合潮類型海域，依據(jù)不規(guī)則半日潮類型與不規(guī)則日潮類型對(duì)應(yīng)選擇。因此，平均大潮高潮面的定義依據(jù)潮汐類型，而本質(zhì)上代表的特征潮位面并不相同，這與美國海岸線定義中的平均高潮面相似，其在日潮類型海域?qū)嵸|(zhì)是平均高高潮面[9-10]。

3.2 凈空基準(zhǔn)值的統(tǒng)計(jì)算法

由平均大潮高潮面定義可知，其計(jì)算應(yīng)采用實(shí)測水位數(shù)據(jù)或預(yù)報(bào)潮位的統(tǒng)計(jì)算法。算法的關(guān)鍵是推算或判斷大潮或回歸潮出現(xiàn)日期。推算是指首先由潮汐類型決定采用大潮或回歸潮，進(jìn)而分別依朔望(或陰歷)或月球赤緯推算日期。判斷是指直接由潮差大小判斷大潮或回歸潮出現(xiàn)日期，不判斷潮汐類型也不計(jì)算月相與月球赤緯，可作為所有潮汐類型統(tǒng)一的統(tǒng)計(jì)算法。每次大潮或回歸潮應(yīng)選取前后共3天的高潮或高高潮，取19年的平均值，該平均值與當(dāng)?shù)亻L期平均海面在垂直方向上的差距即為凈空基準(zhǔn)值。

3.3 海島岸線高程推算

平均大潮高潮面相對(duì)于平均海面的高度、影像拍攝時(shí)刻瞬時(shí)水位相對(duì)于平均海面的高度，以及平均大潮高潮面相對(duì)于瞬時(shí)水位的高差，都可由精密海潮模型推算。將立體測圖過程中得到的影像拍攝時(shí)刻瞬時(shí)水位高程與該高差求和，就得到平均大潮高潮面的高程，即海島岸線高程。顯然，這里得到海島岸線高程的高程基準(zhǔn)也與測圖高程基準(zhǔn)一致。

4 測試分析與結(jié)論

本文采用航空遙感數(shù)據(jù)對(duì)某海島及周邊區(qū)域進(jìn)行1∶2000比例尺測圖(如圖2所示)，在立體測圖環(huán)境中按上述方法測定影像拍攝時(shí)刻海島瞬時(shí)水位線的高程為1.69 m。

由本區(qū)域海洋潮汐模型、海面地形模型推算出影像拍攝時(shí)刻海島瞬時(shí)水位線的高程為1.78 m，平均大潮高潮面高程(海島岸線高程)為2.01 m，進(jìn)而得到海島岸線與瞬時(shí)水位線的高差為2.01-1.78=0.23 m。

由立體測圖環(huán)境中瞬時(shí)水位線的高程和推算求得的海島岸線與瞬時(shí)水位線高差求得，立體測圖環(huán)境中海島岸線高程值為1.69+0.23=1.92 m。

在立體測圖環(huán)境中，由測圖海島岸線高程值1.92 m，按該方法跟蹤出海島岸線平面位置矢量(如圖2所示)。

為測試遙感測圖海島岸線測量性能，選擇海島背風(fēng)處海岸，采用GPS RTK技術(shù)實(shí)測痕跡岸線點(diǎn)400余個(gè)，將痕跡岸線點(diǎn)大地高減去大地水準(zhǔn)面高(16.12 m)，以2.5倍中誤差為限差，迭代求取平均值，求得實(shí)測痕跡岸線高程為2.07 m。與推算的海島岸線高程相差2.07-2.01=0.06 m，與立體測圖中海島岸線高程相差2.07-1.92=0.15 m。

綜合以上面分析結(jié)果，可以看出：

(1) 海島痕跡岸線在現(xiàn)場難以辨認(rèn)，RTK測量方式獲取的痕跡岸線不光滑，細(xì)節(jié)表達(dá)不清，尤其在狹長的海灣及礁石眾多的地方，誤差較大。

圖2 立體測圖獲取的海島岸線與RTK測量的痕跡岸線對(duì)比

(2) 立體測圖環(huán)境中測量的海島岸線具有連續(xù)、精細(xì)的地理特征，在較高精度海潮模型和海面地形支持下，高程精度能優(yōu)于0.2 m。

(3) 測試結(jié)果顯示，實(shí)測痕跡岸線高程、岸線推算高程與立體測圖環(huán)境下海島岸線高程互差不大于0.2 m。

[1] 國家質(zhì)量技術(shù)監(jiān)督局.海洋學(xué)術(shù)語·海洋地質(zhì)學(xué)：GB/T 18190—2000[S]．北京: 中國標(biāo)準(zhǔn)出版社，2000．

[2] 黨亞民，程鵬飛，章傳銀，等.海島礁測繪技術(shù)與方法[M].北京：測繪出版社, 2012.

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

[4] LI R，MA R，DI K.Digital Tide-coordinated Shoreline [J]. Journal of Marine Geodesy，2002，25(1)：27-36.

[5] 許家琨，劉雁春，許希啟，等.平均大潮高潮面的科學(xué)定位和現(xiàn)實(shí)描述[J].海洋測繪，2007，27(6)：19-24.

[6] 申家雙，郭海濤，宋瑞麗，等.海岸帶水邊線等高條件控制下的航空影像外部定向[J].海洋測繪，2011，31(4)：22-25.

[7] 王之卓.攝影測量原理[M].武漢：武漢大學(xué)出版社，2007.

[8] 侯世喜，黃辰虎，陸秀平.基于余水位配置的海洋潮汐推算研究[J].海洋測繪，2005，25(6)：29-33.

[9] 暴景陽，劉雁春，晁定波，等.中國沿岸主要驗(yàn)潮站海圖深度基準(zhǔn)面的計(jì)算與分析[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版)，2006，31(3)：224-228.

[10] 暴景陽，章傳銀.關(guān)于海洋垂直基準(zhǔn)的討論[J].測繪通報(bào)，2001(6)：10-11.

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

[12] 謝秀琴. 基于遙感圖像的海岸線提取方法研究[J]. 福建地質(zhì), 2012, 31(1):60-66.

[13] 馬小峰. 海岸線衛(wèi)星遙感提取方法研究[D].大連：大連海事大學(xué), 2007.

[14] LUQUE I F, TORRES F J A, TORRES M A A, et al. A New Robust and Accurate Method to Extract Tide-coordinated Shorelines from Coastal Elevation Models[J]. Journal of Coastal Research, 2005, 28(3):683-699.

ShorelineSurveyingMethodBasedontheStereoImagingandMapping

DANG Yamin1，ZHANG Chuanyin1，ZHOU Xinghua3，XUE Shuqiang1，XU Jun2

(1. Chinese Academy of Surveying and Mapping, Beijing 100830,China； 2. PLA Dalian Naval Academy, Dalian 116018，China； 3. The First Institute of Oceanography,State Oceanic Administration, Qingdao 266061，China)

A new method for accurately surveying and mapping the shoreline of an island is proposed where the theoretically shoreline defined by the high water (MHHW) is adopted as the shoreline indicator of the island. This method utilizes the height of the instantaneous shoreline extracted from the aerial images to obtain the height of the shoreline by combining the precise ocean tide model and the height of the instantaneous shoreline, and then the horizontal position of the shoreline can be tracked by the stereo imaging and mapping method with the obtained height of MHHW. The proposed method can keep and ensure the continuity and uniqueness of the shoreline relative to the field measurements for shoreline mapping, and this method can be widely applied to large-scale shoreline mapping of mainland and islands.Experimental test indicates that the height accuracy can reach 0.2 m and it can be applied to the island mapping of the scale 1∶2000.

island shoreline; stereo mapping；instantaneous water level；ocean tide

黨亞民，章傳銀，周興華，等.海島岸線遙感立體測圖精細(xì)測量方法[J].測繪通報(bào)，2017(11)：47-50.

10.13474/j.cnki.11-2246.2017.0346.

P229

A

0494-0911(2017)11-0047-04

2017-03-09

黨亞民(1965—)，男，研究員，主要從事GNSS精密定位、大地測測量基準(zhǔn)、海島礁測繪和地球動(dòng)力學(xué)領(lǐng)域的研究工作。E-mail: dangym@casm.ac.cn