秦 菡 曾維特 李華彥 楊新發(fā) 張同同 鄭建宜 沈愛斯

(1. 海南省海洋地質資源與環(huán)境重點實驗室, 海南 ?？?570206;2. 海南省海洋地質調查研究院, 海南 ?？?570206)

0 引言

海岸帶是海陸相互作用的交集地帶,也是生態(tài)脆弱災害多發(fā)的敏感區(qū)域[1],海岸帶環(huán)境的演化直接影響了人類的生活質量和生存發(fā)展,科學認識海岸帶地質環(huán)境的形成和演化規(guī)律及其人類活動對海岸帶的影響,是實現(xiàn)人地協(xié)調、可持續(xù)發(fā)展的基礎和前提。目前,海南島砂質岸線的主要地質環(huán)境問題為:岸坡變陡、岸灘變窄,岸灘砂質細化、黑化,岸線侵蝕后退等。為了進一步做好海岸線地質資源開發(fā)利用與保護、海岸帶生態(tài)環(huán)境保護與修復、近岸重大工程基礎設施地質安全和海岸帶地質災害防治等工作,需要準確及時的空間地理信息作為支持。

與傳統(tǒng)的依靠全站儀、水準儀、全球定位系統(tǒng)(global positioning system,GPS)等測量工具的全野外大比例尺地形圖測圖方式相比,無人機(unmanned aerial vehicle,UAV)攝影測量[2-3]作為載人遙感和衛(wèi)星遙感的補充手段,具有周期短、效率高、成本低、操控簡單、受地面狀況和天氣情況影響較小、數(shù)據(jù)更新較為靈活簡便等特點,同時配合全球導航衛(wèi)星系統(tǒng)(global navigation satellite system,GNSS)、慣性測量單元(inertial measurement unit,IMU)等技術來獲取測區(qū)影像數(shù)據(jù),經(jīng)一系列內業(yè)處理工作[3]后,獲取地物的空間地理信息數(shù)據(jù),從而廣泛應用于小范圍大比例尺地形圖測量[4],在海岸帶測量中成為不可或缺的技術手段[5]。

1 實驗區(qū)概況與作業(yè)內容

1.1 實驗區(qū)概況

受到豐富的地表水和強烈的化學風化作用,海南島的海岸線每年均會受到不同程度的侵襲,而?？谑形骱０短幱诤？谑幸?guī)劃的濱海旅游發(fā)展軸線和城市發(fā)展軸線之間,區(qū)域資源和區(qū)域位置突出。因此本文將位于海口西海岸的國際會展中心(110°11′21″E、20°3′41″N附近)周邊區(qū)域作為實驗區(qū),重點關注沙灘被侵蝕的動態(tài)變化,為西海岸岸線治理等工程提供數(shù)據(jù)指導。

1.2 技術路線

無人機攝影測量以無人機為飛行平臺,搭載數(shù)碼相機等傳感器設備,在目標區(qū)域進行航攝飛行,獲取測區(qū)內影像數(shù)據(jù)及每張相片對應的定位定姿系統(tǒng)(position and orientation system, POS)點數(shù)據(jù),再對數(shù)據(jù)進行空中三角測量等內業(yè)處理,從而獲得數(shù)字表面模型(digital surface model,DSM)、數(shù)字正射影像(digital orthophoto map,DOM)、數(shù)字線劃圖(digital line graphic,DLG)等基礎測繪產品,為海岸帶地質環(huán)境調查提供基礎地理資料[6]。技術路線見圖1。

圖1 海岸帶地質環(huán)境調查航測技術路線

1.3 野外作業(yè)內容

本文所涉及的岸段測量范圍僅為作業(yè)期間海岸帶地區(qū)露出水面的陸地部分,即岸線向陸地延伸200 m的陸地部分;使用大疆Phantom 4實時動態(tài)載波相位差分技術(real time kinematic,RTK)多旋翼無人機配合單鏡頭相機及中海達全球導航衛(wèi)星系統(tǒng)接收機(Ubase)硬件進行航攝。

重點岸段附近的地形復雜多變,水域面積較大,岸線彎曲程度不同且現(xiàn)有的地物特征較少,故需制定針對性的作業(yè)流程[7]。

在正式飛行前,需充分整理測區(qū)相關航攝資料、已知控制點點之記、地形圖等資料,根據(jù)測區(qū)地形地貌、測圖比例尺等進行航攝分區(qū)、飛行架次、相控布設點、航向重疊和旁向重疊以及曝光間隔、航線間隔的設定。

1.3.1布設控制點

控制點布設可采用區(qū)域網(wǎng)五點法布設控制點[8],即測區(qū)四角4個控制點,中心1個控制點,其中四角控制點處于兩個航測相鄰區(qū)域時,像控點盡量布設在航線重疊結合處,鄰區(qū)盡量公用。

根據(jù)實際測量過程中存在位置難以到達、像主點落水等問題,在滿足空中三角測量以及立體測圖的要求下,根據(jù)岸段形狀、植被分布等情況進行控制點布設,對于經(jīng)濟發(fā)達的海岸帶地區(qū),選取相對固定且易于準確量測和定位的地方或較易布設與容易到達的地方均勻布設像控點;對于較難到達或測區(qū)內部分區(qū)域難以到達以及周圍地物不能明顯區(qū)分的海岸帶地區(qū),特征點較少,使用半徑0.15 m的圓形像控布全野外布點測量,像控點位置應做到大小適宜、易于辨認、不易遮擋,并位于高潮線以上;同時能最大限度均勻控制測繪面積,本文共布設5個控制點及11個檢查點,如圖2所示。

圖2 控制點與檢查點示意圖

1.3.2無人機數(shù)據(jù)采集

按規(guī)定,無人機飛行需進行空域申請[9],但本次工作飛行高度為100 m,滿足“120 m以下非機場軍方、政府、監(jiān)獄等禁飛區(qū)不用申請空域”的條件,因此無須做空域協(xié)調。現(xiàn)場作業(yè)前,對飛行平臺、航攝儀、飛控系統(tǒng)等進行檢查和維護;在適宜位置架設Ubase接收機測站采集靜態(tài)數(shù)據(jù);采用垂直起降方式起飛,按預設航線飛行過程中時刻監(jiān)視飛行高度、速度、飛行姿態(tài)等;起飛點與降落點在同一位置,起降點附近無遮擋,平整開闊,GNSS衛(wèi)星信號好。

為控制數(shù)據(jù)采集速度,同時兼顧無人機航拍的承載能力,航線設置為航向重疊率80%,旁向重疊率70%,飛行高度為100 m,飛行速度為7.5 m/s,距地面最大安全高度為120 m。

為控制實驗區(qū)內的相片質量,航線邊距在測區(qū)邊線的基礎上外擴50 m;采用定距模式拍照,照片比例為4∶3,云臺角度為-90°;采用中央平均的測光模式,防止過曝或者欠曝。

無人機降落后,要及時導出無人機數(shù)據(jù)(影像數(shù)據(jù)、機載端數(shù)據(jù)、拍照點時刻文件)及基站數(shù)據(jù)(基站靜態(tài)文件、基站對應的基站坐標)并檢查,若有不合格數(shù)據(jù),應依據(jù)補飛方案進行補飛。

1.4 數(shù)據(jù)處理內容

1.4.1后差分處理

在后差分處理(post processed kinematic,PPK)軟件中導入無人機及基站數(shù)據(jù),采用自主高精度的PPK算法向前和向后解算,采用高精度的插值算法,獲取目標橢球BLH及ENZ坐標格式的POS信息,使曝光瞬間的位置更可靠。

1.4.2空中三角測量

光線束法區(qū)域網(wǎng)空中三角測量是以投影中心點、像點以及對應的地面點為一條直線成立的共線方程,該方法以單張相片為解算單元,借助相片之間的重疊部分和控制點把所有相片的光束連成一個區(qū)域進行整體平差,解算相片中各點的三維坐標[10]。結合解算后的POS文件、控制點文件,經(jīng)過影像特征點提取、同名點匹配等步驟實現(xiàn)自由網(wǎng)空三測量。由于自由網(wǎng)空三測量后需要人機交互進行控制點配準工作,所以需要進行二次甚至多次全自動空三加密優(yōu)化[11-12],直至空三測量的精度滿足要求。

1.4.3立體模型建立

根據(jù)空三加密測量成果進行模型分塊,并在此基礎上,全自動提取測區(qū)密集點云,以便進行點云分類、影像圖制作、可視化及分析等;生成不規(guī)則三角網(wǎng),同時利用鏡頭獲取的紋理信息,自動貼圖,得到測區(qū)三維網(wǎng)格模型。

1.4.4數(shù)字成果生成

采用立體測圖生產數(shù)字成果可免去大量的外業(yè)實測及調繪工作,大大提高了生產的工作效率。本文實驗在完成立體建模后,自動生成實驗區(qū)可互操作的光柵圖層—正射影像DOM及數(shù)字表面模型DSM[13]。

DLG數(shù)字成果借助于立體測圖軟件中的三維模塊,將osgb格式的三維網(wǎng)格模型轉換成DSM表面模型,借助模型的空間尺寸信息,利用三維實景模型的空間量測功能,直接進行地形、地物的量算及采集工作。

2 精度評價

2.1 空三精度評價

P01、P02、P03點為沙灘控制點,P04、P05點為水泥路控制點,各點的中誤差見表1。

表1 控制點的平面位置與高程中誤差 單位:m

5個控制點中,平面位置中誤差最高為0.029 m,高程中誤差最高為0.045 m;整體平面位置中誤差為0.021 m,高程中誤差為0.024 m。

看出11個檢查點的位置,其中6個點在實驗區(qū)域內,4個點在實驗區(qū)邊界線外擴50 m范圍內,1個點位于實驗區(qū)擴充區(qū)域外,各點的中誤差見表2。

表2 檢查點的平面位置與高程中誤差 單位:m

如表2所示,11個檢查點中,平面位置中誤差最高為0.059 m,高程中誤差最高為0.149 m;整體平面位置中誤差為0.031 m,高程中誤差為0.065 m。根據(jù)控制點及檢查點中誤差可知,空三精度較高,滿足《數(shù)字航空攝影測量空中三角測量規(guī)范》的1∶1 000地形圖測圖要求。

2.2 影像精度評價

使用RTK實測地物特征點與影像圖上的同名像點進行對比分析來檢查影像的平面位置及高程精度,如表3所示。

表3 影像精度分析驗證 單位:m

通過16個均勻布設的RTK實測點與DOM、DSM中同名像點做對比,由上表可以看出,平面位置中誤差為0.044 m,高程中誤差為0.048 m,滿足規(guī)范要求。

2.3 數(shù)字線劃圖精度評價

使用RTK實測地物特征點與1∶1 000地形圖上的同名點進行對比分析來檢查影像的平面位置及高程精度(地形圖上同名點的高程可通過三維模型直接量取),如表4所示。

表4 DLG精度分析驗證 單位:m

通過計算得到平面位置中誤差為0.041 m,高程中誤差為0.034 m,查找規(guī)范可知1∶1 000比例尺平面、丘陵地的點位中誤差為0.6 m,平地的高程中誤差為0.2 m,對比得到0.041<0.6,0.034<0.2,因此滿足規(guī)范要求。

3 數(shù)字成果精度分析

將影像和數(shù)字線劃圖的平面位置中誤差及高程中誤差作為精度評定指標,可發(fā)現(xiàn)數(shù)字成果均能夠滿足規(guī)范及作業(yè)要求。

使用影像和數(shù)字線劃圖的點位誤差及高程誤差進行精度分析,以此來確定最優(yōu)控制點及檢查點的布設方式及位置,為大面積海岸帶地質環(huán)境調查提供理論依據(jù),如圖3所示。

由圖3可知,C01RTK點的平面位置誤差均大于0.12 m;C11號RTK點的高程誤差變化在0.07～0.11 m之間。

(a)影像平面位置精度統(tǒng)計折線圖

(b)影像高程精度統(tǒng)計折線圖

(c)DLG點位精度統(tǒng)計折線圖

(d)DLG高程精度統(tǒng)計折線圖

根據(jù)實際測量環(huán)境分析,檢查點C01布設在水泥路上的雨水篦子邊角處,控制點與周邊地物顏色變化不明顯,導致人機交互誤差較大,因此該點處平面精度誤差較大,且該點處位于水泥路上面,因此高程精度影響不大;而檢查點C11點位于會展中心高層建筑物的中心位置,周邊區(qū)域有一個高程垂直變化的現(xiàn)象,對該點位高程精度有較大影響。

實驗區(qū)外擴區(qū)域內外的4個檢查點C02、C03、C04、C09的影像平面位置中誤差為0.039 m,高程中誤差為0.044 m;DLG平面位置中誤差為0.032 m,高程中誤差為0.017 m,在規(guī)劃測區(qū)以外但無人機可以拍攝到的區(qū)域,數(shù)字成果的精度均滿足規(guī)范要求。

由此可驗證,布設像控點時要避免有高層建筑,且要與周邊地物明顯區(qū)分;若布設的像控點得當合理,航攝范圍線外擴50 m范圍內的數(shù)字成果亦可滿足精度要求,在數(shù)字成果分幅、裁切時可保留。

4 結束語

實驗將影像同名像點、DLG同名點與實測RTK測量點做對比,雖個別結果受周邊地形與人機交互的影響,點位誤差相差較大,但仍在規(guī)范范圍內,點位精度和高程精度均滿足生產測量要求,因此無人機低空遙感技術能夠完全滿足海岸帶地質環(huán)境調查的需要,所獲數(shù)據(jù)形式多元化,可用于測區(qū)的全方位分析,并且能夠大幅度減少外業(yè)工作量,提高工作效率。

本文在數(shù)據(jù)處理過程中發(fā)現(xiàn),點云數(shù)據(jù)未做到精細分類,房屋、樹木等對DEM制作影響較大,影像反映的是地表高程數(shù)據(jù),而并未最真實地進行地面地形數(shù)字化模擬,因此無人機航測技術在海岸帶地質環(huán)境調查過程中如何實現(xiàn)點云分類,將是下一步工作的重點。