呂艷兵
(忻州市水利勘測設計院,山西 忻州 034000)
水利工程確權(quán)劃界,是依法對河道及水利工程劃定管理和保護范圍,是有效推進河長制工作的重要手段。傳統(tǒng)測量技術(shù)成本高、效率低,已無法滿足現(xiàn)代工程管理需要。無人機航測具有成本低,效率高等優(yōu)勢。將無人機航測技術(shù)應用到水利工程勘測中,在保證勘測成果質(zhì)量的同時,可提高勘測作業(yè)效率。
工程平面坐標系統(tǒng)采用CGCS 2000,中央子午線為111°,高斯投影為3°帶;高程采用1985 國家黃海高程;航測數(shù)字化地形圖成圖比例尺為1∶2000,等高距為2 m;測圖范圍包括運行區(qū)、工程區(qū)(從劃界范圍的周邊外延300 m)、庫區(qū)周邊(從校核水位外延300 m);數(shù)字線劃圖(DLG)為AUTOCAD 2004 的DWG 數(shù)據(jù)格式;數(shù)字正射影像圖(DOM)為TIF 數(shù)據(jù)格式(0.07 m分辨率)。
本次測量的D 級控制網(wǎng),選取了14 點組成D 級網(wǎng)進行觀測,并按點名進行編號(FD 01—13),點間距為4~6 km。GPS 網(wǎng)使用三臺南方GPS 接收機(儀器標稱精度:靜態(tài)定位±5 mm+2 ppm),三臺華測GPS 接收機(儀器標稱精度:靜態(tài)定位±5 mm + 2ppm)進行觀測,采用三點推進邊連式構(gòu)網(wǎng)法同步觀測,觀測時,有效衛(wèi)星數(shù)應不少于4 顆,衛(wèi)星高度角最小為15°,數(shù)據(jù)間隔為15 s,時段長度大于1 h,從天線的三面分別測量天線高三次,取平均值。
實測GPS 控制點15 個,其中已知控制點2 個,未知點13 個,組成同步環(huán)66 個,異步環(huán)38 個。獨立基線20 條,其中必要基線14 條,多余基線36 條,平均重復設站數(shù)為1.6/站,符合《全球定位系統(tǒng)(GPS)測量規(guī)范》規(guī)定1.6/站[1],詳見表1。
表1 GPS 觀測網(wǎng)的概述
外業(yè)觀測后所有數(shù)據(jù)轉(zhuǎn)化為Rinex 數(shù)據(jù),用計算機對外業(yè)數(shù)據(jù)進行檢查,檢查得中誤差小于40 mm,方差比大于3。此次觀測共形成52 條合格基線,基線最短邊長349 m,最大邊長7670.4 m,平均邊長4 047 m,詳見表2。
表2 GPS網(wǎng)的基線
本次觀測共形成66 個同步環(huán),同步環(huán)最優(yōu)閉合差1.768 mm,最弱閉合差33.018 mm。所有同步環(huán)坐標分量相對閉合差和全長相對閉合差均小于閉合限差。
本次觀測形成的異步環(huán)38 個。異步環(huán)最優(yōu)閉合差6.443 mm,最弱閉合差36.577 mm。網(wǎng)形質(zhì)量良好,構(gòu)網(wǎng)基線合理?;€質(zhì)量統(tǒng)計見表3。
對GPS 網(wǎng)三維無約束平差后,點位中誤差為±0.0033 cm、最弱點為FD 01?;€相對精度最弱邊YSN 1392141C-FD 13141BN、長度979.796 m,相對精度1/92146,表明控制網(wǎng)的內(nèi)符合精度較高。
表3 基線質(zhì)量統(tǒng)計
CGCS 2000 坐標系下無約束平差,中央子午線為111°,高斯投影為3°帶,固定平面起算點YSN 進行無約束平差,見表4。
表4 GPS 網(wǎng)精度統(tǒng)計
原高程基準為大沽高程,本次測量采用1985 高程基準,以三等水準聯(lián)測水庫高程基準BM 02 與國家Ⅱ等水準點靜汾15 后,1985 高程比大沽高程高1.518 m,D 級控制網(wǎng)的高程采用原四等水準網(wǎng)成果。
測區(qū)布設像控點89 個,檢查點85 個,像控點及檢查點采用GPS RTK 測量。像控點均勻布設于水庫周邊。利用數(shù)字空三進行平高加密。像控點打小木樁和地面齊平并且做好L 形白粉標志,像控點一律采用順序號前冠英文字母XK 方法編號,整個測區(qū)沒有重號。檢查點采用順序號前冠英文字母JC 方法編號。
在D 級GNSS 點的基礎上,利用SXCORS RTK 架設三腳架三次初始化測量像控點,84 坐標系與GCCS 2000 坐標系的轉(zhuǎn)換關(guān)系的獲取方法,采用在測區(qū)現(xiàn)場,通過點校正的方法獲取。轉(zhuǎn)換參數(shù)的求解,采用13 個D 級點,并且能覆蓋整個測區(qū),通過合適的數(shù)學模型及可靠性檢驗,對多點組合方式進行計算和優(yōu)選。通過RTK 測量像控點坐標時,其轉(zhuǎn)換殘差應不大于圖上±0.1 mm。RTK 測量像控點高程,擬合殘差為等高距的1/10[2]。全部像控點又經(jīng)過山西省測繪工程院轉(zhuǎn)換GC 2000 坐標系和1985 國家高程基準,較差在允許范圍。
采用大鵬無人機系統(tǒng)(CW-10)搭載嚴密檢校過的高分辨率數(shù)碼相機Sony ILCE7R,完成所有航攝工作。按航片地面平均分辨率為0.07 m,相對航高為513 m 飛行。結(jié)合測區(qū)地形高差不應大于1/6 航高,分為3 個攝區(qū),共飛行8 架次,完成航測面積為102 km2。攝影測量時間選擇,應確保影像能真實地顯現(xiàn)地面細部,太陽高度角應大于30°,陰影應小于2 倍時進行[3-4]。
2.4.1 GPS 數(shù)據(jù)采集
為提高搜索空間,無人機飛行時,差分GPS 天線置于飛機前方;相機拍下的每張照片,都和飛控系統(tǒng)記錄的坐標點一一對應,飛行結(jié)束后,如果從飛控系統(tǒng)下載的POS 數(shù)據(jù)和照片的數(shù)量對應,則表明GPS數(shù)據(jù)采集正確。
2.4.2 照片數(shù)據(jù)質(zhì)量
影像質(zhì)量特別強調(diào)影像清晰,反差適中,顏色飽和,色彩鮮明,色調(diào)一致。有較豐富的層次、能辨別與地面分辨率相適應的細小地物影像,滿足外業(yè)全要素精確調(diào)繪和室內(nèi)判讀的要求。
2.4.3 飛行質(zhì)量
無人機飛行選在風速不大且穩(wěn)定,日照充足的時段,飛行質(zhì)量較高。相片的航測范圍覆蓋了整個測區(qū),無漏測區(qū)域。航測時航向重疊率大于等于80%;旁向重疊率大于等于75%。經(jīng)檢查,所有航測數(shù)據(jù)完整清晰,滿足精度和規(guī)范的要求。
為了正射影像精度小于4 個像素,DOM 和DLG成圖精度可達1∶2 000 的要求,選用Pix4D 進行數(shù)據(jù)處理。
影像預處理包括增強圖像、相片變形校正、編輯等工序,預處理后的影像更加鮮明,有利于后期空三加密。
在測區(qū)已知點架設GPS 接收機,為靜態(tài)采集模式。參數(shù)設置為:在無人機起飛前直至降落持續(xù)采樣,間隔為1 s。飛行時公共衛(wèi)星數(shù)多于4 顆,且解算全過程顯示綠色,說明數(shù)據(jù)質(zhì)量較好,解算得到WGS 84 坐標下的POS 坐標數(shù)據(jù)。
為達到1∶2 000 的精度要求,空三加密采用Pix4D進行,且刺點的誤差最大不能超過3 個像素[5]。
在武漢點云科技有限公司開發(fā)的SKYPHOTO-MAP3D 軟件下,基于數(shù)字正射影像DOM 和數(shù)字高程模型DEM 生成實景三維模型DSM,提取檢查點坐標,與檢查點RTK 實測坐標比較,進行DSM 質(zhì)量評價。經(jīng)計算,檢查點平面中誤差±0.015 cm,高程中誤差0.045 cm。
數(shù)據(jù)采集在武漢點云科技有限公司開發(fā)的SKYPHOTO-MAP3D 軟件下,基于數(shù)字正射影像DOM和數(shù)字高程模型DEM 生成實景三維模型DSM,點狀地物在地物的定位點上采集,線狀地物在地物中心線采集,偏差不能超過圖上0.1 mm。
DLG 成圖軟件使用《南方CASS10.1 地形圖系統(tǒng)》,按《地形圖圖式》標準進行修改直至滿足成圖要求。
生成DEM 后,可進行DOM 的制作。生成DOM后,要對其拉絲、劃痕、變形及重疊進行檢查,并用PHOTOSHOP 軟件進行調(diào)色和修飾。最后在圖幅上加地名、圖名等打印成圖。