邢喆，奚歌，樊妙，李鵬，李四海

（國家海洋信息中心，天津 300171）

隨著我國經(jīng)濟迅速發(fā)展，陸地不可再生資源的日益枯竭，海洋開發(fā)戰(zhàn)略已成為我國經(jīng)濟可持續(xù)發(fā)展的重要因素。由國家海洋局承擔(dān)的“我國近海海洋綜合調(diào)查與評價”專項突出了海洋開發(fā)戰(zhàn)略的主題，“數(shù)字海洋”作為其中的重要組成部分，得到專項的足夠重視。構(gòu)建“數(shù)字海洋”首先需要建立地理空間數(shù)據(jù)框架，其核心是在陸海垂直基準(zhǔn)統(tǒng)一基礎(chǔ)上，實現(xiàn)陸海地形數(shù)據(jù)的無縫集成。

通過本次專項調(diào)查獲取了大量我國海島海岸帶區(qū)域航空LIDAR 數(shù)據(jù)和近海海底多波束/單波束水深數(shù)據(jù)，然兩者高程基準(zhǔn)并不統(tǒng)一，專項獲取的航空LIDAR 數(shù)據(jù)以大地高作為高程基準(zhǔn)，水深資料則采用深度基準(zhǔn)。以上兩種基準(zhǔn)面作為海陸高程/水深的起算面，不可避免存在陸地高程和海洋深度數(shù)據(jù)不能有效銜接的缺陷，故亟待構(gòu)建航空LIDAR 數(shù)據(jù)和水深數(shù)據(jù)間的轉(zhuǎn)換關(guān)系，將其歸算到統(tǒng)一的高程基準(zhǔn)下，為“數(shù)字海洋”基礎(chǔ)框架構(gòu)建提供可靠的數(shù)據(jù)支撐。

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

1.1 機載LIDAR 數(shù)據(jù)

機載LIDAR 是一種安裝在飛機上的激光探測和測距系統(tǒng)，無需任何或僅需少量地面控制點，獲得高精度、高密度的三維坐標(biāo)數(shù)據(jù)，可用來構(gòu)建目標(biāo)物的三維立體模型（張玉方等，2008）。LIDAR系統(tǒng)具有自動化程度高、受天氣影響小、數(shù)據(jù)生產(chǎn)周期短等特點，為獲取高分辨率地球空間信息，提供了一種全新的技術(shù)手段。LIDAR 數(shù)據(jù)的處理包括航跡重構(gòu)、激光掃描測量數(shù)據(jù)處理、“點云”數(shù)據(jù)處理以及DEM 構(gòu)建等過程，流程如圖1 所示。

圖1 LIDAR 數(shù)據(jù)處理流程

1.2 海底地形資料

自單波束和多波束海底聲吶探測技術(shù)被廣泛應(yīng)用以來，獲取全面、精細的海底水深數(shù)據(jù)變得切實可行。數(shù)據(jù)采集階段結(jié)束以后，原始數(shù)據(jù)首先經(jīng)導(dǎo)航信號編輯、噪聲編輯、吃水改正、聲速改正、潮位改正等前期的數(shù)據(jù)處理過程，獲得大量水深離散點，再用有一定結(jié)構(gòu)關(guān)系的網(wǎng)格點代替離散點，通過數(shù)值插值方法對標(biāo)準(zhǔn)化的海底水深離散數(shù)據(jù)進行內(nèi)插網(wǎng)格化處理，從而構(gòu)建生成海底DEM，進行各種應(yīng)用。數(shù)據(jù)處理流程如圖2 所示。

圖2 海底地形數(shù)據(jù)處理流程

2 海陸地形轉(zhuǎn)換關(guān)系實現(xiàn)

由于上述兩種資料采用基準(zhǔn)面的不同，在垂直方向存在較大差異，無法滿足航空DEM 和海底DEM 拼接的需要，必須建立兩種數(shù)據(jù)源垂直基準(zhǔn)向地理信息表達通用高程基準(zhǔn)的轉(zhuǎn)換數(shù)學(xué)模型。

在資料采集及資料成果在垂直方向上主要涉及以下基準(zhǔn)面（許家琨等，2011），如圖3 所示。

圖3 垂直方向上各基準(zhǔn)面關(guān)系圖

圖3 中，H 為LIDAR 觀測數(shù)據(jù)大地高；Hr 為正常高；ζ 為高程異常值；X 為驗潮站水準(zhǔn)點與1985 國家高程的垂直距離，即驗潮站水準(zhǔn)點在1985 國家高程基準(zhǔn)中的高程；ΔH 為驗潮站水準(zhǔn)點與水尺驗潮零點的垂直距離；MSL 為平均海面到水尺零點的垂直距離；Δh 為當(dāng)?shù)仄骄Ｃ媾c1985國家高程基準(zhǔn)的垂直距離，即海面地形；h 為從深度基準(zhǔn)面起算的成果水深；h1 為海底在1985 國家高程基準(zhǔn)中的高程負值，即海陸拼接后使用統(tǒng)一高程基準(zhǔn)的海底地形高程。

由圖3 知，以下關(guān)系式成立：

通過式（1） 可將LIDAR 測得的大地高歸算到正常高，通過式（2） 可將水深成果從理論深度基準(zhǔn)面歸算至1985 國家高程，至此構(gòu)建了區(qū)域高程基準(zhǔn)和深度基準(zhǔn)的相互轉(zhuǎn)換數(shù)學(xué)模型。

由式（1）、（2） 可知，實現(xiàn)此轉(zhuǎn)換關(guān)鍵在于如何準(zhǔn)確計算ζ、Δh 及L。

（1） ζ 的計算方法：采用GPS 觀測和水準(zhǔn)聯(lián)測方法獲取測區(qū)內(nèi)的大地高和正常高，通過高程擬合，建立兩者之間的轉(zhuǎn)換模型。

（2） L 的計算方法：我國深度基準(zhǔn)面采用由13 個分潮（Sa，Ssa，Q1，O1，P1，K1，N2，M2，S2，K2，M4，MS4，M6） 確定的理論最低潮面（簡稱TLT）。根據(jù)《海道測量規(guī)范·GB12327—1998》TLT 計算公式計算求得。也可取海圖（及潮汐表）采用值。

（3） Δh 的計算方法：采用幾何水準(zhǔn)測量的方法，按照國家三、四等水準(zhǔn)測量要求，直接聯(lián)測1985 國家高程的水準(zhǔn)點至驗潮站主要水準(zhǔn)點（或工作水準(zhǔn)點） 的高差x，根據(jù)下式計算得到該驗潮站平均海面的高程Δh：

式中，ΔH 通過常規(guī)測量方法可得。

3 應(yīng)用研究

在上海海岸帶區(qū)域，選取12 個GPS 控制點，其中包括3 個上海市E 級GPS 控制網(wǎng)點。以該3個控制網(wǎng)點為基準(zhǔn)點進行聯(lián)測，測量其余9 個聯(lián)測點WGS84 坐標(biāo)（B，L，H），同時按照四等水準(zhǔn)要求進行水準(zhǔn)聯(lián)測，獲得其正常高。同時，在測區(qū)范圍內(nèi)均勻采集了10 個校核點的WGS84 坐標(biāo)，按四等水準(zhǔn)測量施測其正常高。

在上海沿岸布設(shè)5 個臨時水位站，點位布設(shè)如圖4 所示。采用直立式直讀水尺，人工觀測潮水位，潮位觀測漲、落潮均按每半小時觀測一次，為在其周邊2 m 等深線以淺海域單波束測量提供水位改正數(shù)，將測深成果歸算至深度基準(zhǔn)，水深測量測線布設(shè)間隔為100 m。

基于外業(yè)觀測成果，采用二次曲面擬合方法，求解出陸上各觀測點的高程異常值，并通過檢校點進行了精度評定，結(jié)果如圖5 所示。

由圖5 可見，高程擬合的殘差基本在厘米級，中誤差為0.05 m，滿足專項海岸帶調(diào)查DEM 數(shù)據(jù)處理1 ∶10 000 比例尺精度要求（國家海洋局908專項辦公室，2006）。

對臨時水位站的水位觀測值進行了潮位計算及校核，并繪制成相應(yīng)的水位過程線作合理性分析，根據(jù)測量瞬時水深和相應(yīng)潮位計算出測點高程，經(jīng)解算求得Δh = 0.023 m；依據(jù)《2011年潮汐表》中，得L=0.202 m，再由式（2） 即可建立深度基準(zhǔn)面和1985 國家高程基準(zhǔn)之間的關(guān)系。

圖4 上海市海岸帶地區(qū)外業(yè)觀測點位布設(shè)圖

圖5 擬合高程與水準(zhǔn)高程值對比

利用測區(qū)LIDAR 資料和單波束觀測結(jié)果分別構(gòu)建了陸、海DEM，并基于轉(zhuǎn)換關(guān)系進行了拼接，圖6（a） 顯示了基準(zhǔn)轉(zhuǎn)換前陸海拼接模型（大地高和深度基準(zhǔn)面），可以明顯地看出陸海拼接處有一縫隙，地形剖面變化值顯著；圖6（b） 顯示了統(tǒng)一到1985 國家高程基準(zhǔn)后的海陸拼接模型，可以看處，地形剖面變化值平緩，海陸接邊處實現(xiàn)了無縫拼接。

4 結(jié)論與展望

圖6 基準(zhǔn)轉(zhuǎn)換前后海陸拼接模型對比

本文在分析研究海陸地形轉(zhuǎn)換構(gòu)建方法的基礎(chǔ)上，以上海市海岸帶為例，采用高精度GPS 測量和水準(zhǔn)測量，確定了區(qū)域大地高和正常高之間準(zhǔn)確的轉(zhuǎn)換關(guān)系，并結(jié)合潮位數(shù)據(jù)，應(yīng)用水深觀測數(shù)據(jù)，建立起深度基準(zhǔn)面與1985 國家高程基準(zhǔn)之間的轉(zhuǎn)換模型，實現(xiàn)了海陸高程數(shù)據(jù)的無縫拼接，研究成果可為“數(shù)字海洋”構(gòu)建、海洋公益性數(shù)據(jù)服務(wù)、海洋環(huán)境監(jiān)測與預(yù)警等提供基礎(chǔ)技術(shù)支撐。但本文數(shù)據(jù)源多依靠外業(yè)觀測手段獲取，對于資料稀缺或難以施測區(qū)域的應(yīng)用有待進一步深入研究。

GB12327-1998 海道測量規(guī)范：31-32.

鞏淑楠，陳云，徐敏，2010.測繪與空間地理信息，33（5）：165-167．

國家海洋局908 專項辦公室，2006.海島海岸帶航空遙感調(diào)查技術(shù)規(guī)程.北京：海洋出版社.

韓凌云，楊英力，2003.地形圖與海圖拼接的矛盾問題及其處理.海洋測繪，23（4）：34-37．

史照良，曹敏，2007.基于LIDAR 技術(shù)的海島礁_灘涂測繪研究.測繪通報，5：49-53．

許家琨，申家雙，繆世偉，等，2011.海洋測繪垂直基準(zhǔn)的建立與轉(zhuǎn)換.海洋測繪，31（1）：4-8.

陽凡林，李家彪，吳自銀，等，2008.淺水多波束勘測數(shù)據(jù)精細處理方法.測繪學(xué)報，37（1）：444-457.

張玉方，程新文，歐陽平，等，2008.機載LIDAR 數(shù)據(jù)處理及其應(yīng)用綜述．工程地球物理學(xué)報，5（1）：119-124.