杜向鋒,張興福,段 杰
(1.廣東工貿(mào)職業(yè)技術(shù)學(xué)院,廣州 510510;2.廣東工業(yè)大學(xué),廣州 510006;3.廣東省核工業(yè)地質(zhì)局測(cè)繪院,廣州 510800)
利用全球定位系統(tǒng)(global positioning system,GPS)定位技術(shù)可以快速獲得地面點(diǎn)高精度的三維坐標(biāo)——大地緯度B、大地經(jīng)度L和大地高H,通過(guò)平差或坐標(biāo)轉(zhuǎn)換的方法可以快速計(jì)算出該點(diǎn)國(guó)家坐標(biāo)系或地方坐標(biāo)系的高斯平面坐標(biāo),其精度可達(dá)到10-6~10-9量級(jí)[1],高程方面因?yàn)槲覈?guó)采用的是正常高系統(tǒng),故而GPS測(cè)量獲得高精度的大地高需要通過(guò)高程轉(zhuǎn)換才能應(yīng)用于生產(chǎn),目前GPS高程轉(zhuǎn)換理論方法研究已經(jīng)成為大地測(cè)量領(lǐng)域的一個(gè)研究熱點(diǎn),一般說(shuō)來(lái)GPS高程轉(zhuǎn)換的方法有以下4類:模型擬合法、坐標(biāo)轉(zhuǎn)換法、神經(jīng)網(wǎng)絡(luò)法和借助區(qū)域似大地水準(zhǔn)面模型法[2],目前國(guó)內(nèi)許多城市和地區(qū)正在紛紛建立或者已經(jīng)建立起本區(qū)域的高精度的似大地水準(zhǔn)面模型[3-6]。
由于歷史原因,2008年以前韶關(guān)市區(qū)不同地區(qū)之間、同一地區(qū)不同時(shí)期的基礎(chǔ)測(cè)繪數(shù)據(jù)所采用的坐標(biāo)體系互不相同,這一歷史遺留問(wèn)題嚴(yán)重影響了韶關(guān)市國(guó)土資源的統(tǒng)一管理,制約了信息化工作的正常開展,同時(shí)也給社會(huì)各應(yīng)用部門帶來(lái)諸多不便。為了保證測(cè)量工作能真正走在城市建設(shè)各項(xiàng)工作的前面,有效發(fā)揮測(cè)繪部門的職能,及時(shí)地為城市規(guī)劃、建設(shè)和管理提供準(zhǔn)確的測(cè)繪基礎(chǔ)資料,滿足快速發(fā)展的經(jīng)濟(jì)建設(shè)給測(cè)量部門提出的新要求,韶關(guān)市啟動(dòng)了韶關(guān)市基礎(chǔ)控制網(wǎng)改造項(xiàng)目,希望建立起統(tǒng)一的、長(zhǎng)期的、現(xiàn)代化的三維空間基準(zhǔn)框架,并以此為契機(jī)建立韶關(guān)市區(qū)高精度似大地水準(zhǔn)面,項(xiàng)目實(shí)施按照項(xiàng)目?jī)?nèi)容共分兩階段實(shí)施,第一階段是進(jìn)行基礎(chǔ)控制網(wǎng)改造工作,第二階段利用基礎(chǔ)控制網(wǎng)改造的成果計(jì)算似大地水準(zhǔn)面模型并對(duì)模型精度進(jìn)行評(píng)定,本文將詳細(xì)介紹項(xiàng)目實(shí)施的具體過(guò)程和方法。
本次D級(jí)GPS控制網(wǎng)改造共施測(cè)D級(jí)點(diǎn)47個(gè),其中原有網(wǎng)點(diǎn)7個(gè),新增網(wǎng)點(diǎn)40個(gè),聯(lián)測(cè)C級(jí)GPS點(diǎn)6個(gè),外業(yè)觀測(cè)共使用了7臺(tái)Trimble GPS接收機(jī),其中2臺(tái)Trimble 5700雙頻接收機(jī)、2臺(tái)Trimble R8雙頻接收機(jī)和3臺(tái)Trimble SPS781雙頻接收機(jī),各種儀器在使用前均進(jìn)行過(guò)檢驗(yàn),具有可靠的穩(wěn)定性,狀態(tài)良好,觀測(cè)過(guò)程嚴(yán)格按照相關(guān)規(guī)范要求進(jìn)行,為獲得各網(wǎng)點(diǎn)高精度的大地高數(shù)據(jù),外業(yè)觀測(cè)時(shí)段長(zhǎng)度均不少于2h,數(shù)據(jù)處理先利用PINNACLE軟件進(jìn)行基線解算,再利用同濟(jì)大學(xué)的TGPPSw6.0軟件進(jìn)行網(wǎng)平差,計(jì)算結(jié)果表明觀測(cè)數(shù)據(jù)質(zhì)量良好,起算點(diǎn)匹配,各項(xiàng)指標(biāo)均滿足相關(guān)規(guī)范要求,成果合格,其中三維約束平差后大地高精度優(yōu)于1cm。
三等水準(zhǔn)網(wǎng)改造共直接施測(cè)新埋石水準(zhǔn)點(diǎn)35座,直接聯(lián)測(cè)GPS-D級(jí)點(diǎn)23座,聯(lián)測(cè)測(cè)區(qū)范圍內(nèi)的已知I等級(jí)水準(zhǔn)點(diǎn)7個(gè),其它水準(zhǔn)點(diǎn)14個(gè),觀測(cè)水準(zhǔn)路線往返里程共約675km,水準(zhǔn)網(wǎng)計(jì)算前對(duì)各測(cè)段高差觀測(cè)值分別進(jìn)行了尺長(zhǎng)改正和水準(zhǔn)面不平行改正,水準(zhǔn)網(wǎng)計(jì)算采用南方平差易2005進(jìn)行了嚴(yán)密平差,計(jì)算結(jié)果表明各項(xiàng)精度指標(biāo)符合規(guī)范要求,成果合格,其中每千米高差中誤差4.25mm,平均高程中誤差6.50mm。
與此同時(shí)為了似大地水準(zhǔn)面精化需要,部分分布在樓頂無(wú)法直接聯(lián)測(cè)水準(zhǔn)的GPS-D級(jí)點(diǎn),先將高程采用三等水準(zhǔn)引測(cè)到樓下設(shè)置的臨時(shí)水準(zhǔn)點(diǎn)上,再利用水準(zhǔn)儀懸吊鋼尺的方法將高程引測(cè)到樓上GPS點(diǎn),為確保引測(cè)精度,采用多測(cè)回樓上樓下同時(shí)觀測(cè)并往返測(cè)的觀測(cè)方法進(jìn)行施測(cè),數(shù)據(jù)處理顧及了尺長(zhǎng)、溫度和拉力等各項(xiàng)改正,最終獲得26個(gè)D級(jí)GPS點(diǎn)高程值,其中有3個(gè)點(diǎn)為檢查點(diǎn),同時(shí)進(jìn)行了水準(zhǔn)聯(lián)測(cè)和水準(zhǔn)儀懸吊鋼尺聯(lián)測(cè),兩者獲得結(jié)果的互差均優(yōu)于5mm。
3.1.1 重力數(shù)據(jù)
項(xiàng)目組在精化區(qū)域共收集到了107個(gè)點(diǎn)地面重力觀測(cè)數(shù)據(jù),數(shù)據(jù)分布如圖1所示。
圖1 重力數(shù)據(jù)分布圖
3.1.2 地形數(shù)據(jù)
數(shù)字地面模型作為局部似大地水準(zhǔn)面的短波部分的重要影響,在區(qū)域似大地水準(zhǔn)面精化中發(fā)揮重要作用,利用SRTM作為基礎(chǔ)地面模型[7],并進(jìn)行局部改正后作為地形數(shù)據(jù),范圍為112~115°E,23~26°N,分辨率為3″(90m)。地形數(shù)據(jù)如圖2所示。
圖2 數(shù)字地面模型
3.1.3 GPS/水準(zhǔn)數(shù)據(jù)
控制網(wǎng)改造共獲得的GPS/水準(zhǔn)數(shù)據(jù)41個(gè),其分布如圖3所示。
3.2.1 總體思路
基于計(jì)算區(qū)域有一定密度的地面重力數(shù)據(jù),以及離散分布的GPS水準(zhǔn)高程異常數(shù)據(jù),該區(qū)域似大地水準(zhǔn)面精化計(jì)算的總體思路是首先基于莫洛金斯基(Molodesky)邊值問(wèn)題理論[8],利用地面重力數(shù)據(jù)和SRTM數(shù)字地面模型計(jì)算重力似大地水準(zhǔn)面;其次利用GPS水準(zhǔn)數(shù)據(jù)擬合GPS水準(zhǔn)幾何似大大地水準(zhǔn)面,最后利用超定邊值問(wèn)題方法或最小二乘配置方法將重力似大地水準(zhǔn)面和GPS水準(zhǔn)幾何大地水準(zhǔn)面進(jìn)行融合得到最終的GPS水準(zhǔn)重力似大地水準(zhǔn)面。技術(shù)方案如圖4所示,其中參考重力場(chǎng)元采用高精度的EGM2008重力場(chǎng)模型計(jì)算。
圖3 用于水準(zhǔn)面精化的GPS/水準(zhǔn)數(shù)據(jù)分布
圖4 似大地水準(zhǔn)面計(jì)算技術(shù)方案
3.2.2 計(jì)算過(guò)程
(1)計(jì)算重力似大地水準(zhǔn)面,地面重力數(shù)據(jù)處理的目標(biāo)是由離散的地面重力異常數(shù)據(jù)計(jì)算適當(dāng)分辨率的地面平均空間異常,處理內(nèi)容主要包括離散重力異常計(jì)算、平面布格改正、精密地形處理、潮汐基準(zhǔn)轉(zhuǎn)換、橢球校正以及地面重力格網(wǎng)化計(jì)算。當(dāng)計(jì)算區(qū)域面積小于10萬(wàn)km2時(shí),對(duì)于厘米級(jí)似大地水準(zhǔn)面精化目的,由于存在GPS水準(zhǔn)實(shí)測(cè)高程異常的控制數(shù)據(jù),可以省略橢球校正和潮汐基準(zhǔn)轉(zhuǎn)換兩個(gè)步驟,然后采用球面一維FFT算法利用地面平均空間異常與參考重力場(chǎng)的差異(剩余空間異?;驍_動(dòng)重力)按斯托克斯(Stokes)或Hotine積分公式計(jì)算(零階)剩余高程異常,接著進(jìn)行高程異常的地形改正計(jì)算[9],在線性Molodesky理論中,高程異常地形改正是指Molodesky I、Ⅱ階項(xiàng)對(duì)高程異常的貢獻(xiàn)。本項(xiàng)目的高程異常的地形改正計(jì)算方法是用重力異常的地形改正代替Molodesky I階項(xiàng),進(jìn)行地面Stokes積分求得高程異常的地形改正。其中直接計(jì)算Molodesky I、Ⅱ階項(xiàng)及其對(duì)高程異常的影響作為檢核算法,重力地面高程異常(重力似大地水準(zhǔn)面高)等于模型高程異常、零階剩余高程異常和高程異常地形改正三者之和,從而確定出重力似大地水準(zhǔn)面。
(2)計(jì)算GPS水準(zhǔn)幾何似大大地水準(zhǔn)面,利用離散的實(shí)測(cè)GPS/水準(zhǔn)計(jì)算出的高程異常值,通過(guò)空間插值的方式獲得計(jì)算區(qū)域的GPS水準(zhǔn)幾何似大地水準(zhǔn)面模型。
(3)計(jì)算GPS水準(zhǔn)重力似大地水準(zhǔn)面
由于GPS衛(wèi)星定位、高程基準(zhǔn)與全球重力場(chǎng)參考系統(tǒng)不一致,以及水準(zhǔn)高差系統(tǒng)偏差等影響,GPS水準(zhǔn)實(shí)測(cè)高程異常通常與重力地面高程異常存在一定差別,即使采用擬合方法提取了系統(tǒng)偏差信息后,還會(huì)存在殘余的空間噪聲。為有效地消除高程異常的不一致性,需要進(jìn)行GPS水準(zhǔn)重力融合計(jì)算,采用GPS水準(zhǔn)與地面重力聯(lián)合平差的方法進(jìn)行計(jì)算,以GPS水準(zhǔn)高程異常殘差和平均場(chǎng)元改正數(shù)之間的重力場(chǎng)積分公式(Stokes或Hotine積分)為約束條件,以平均場(chǎng)元為觀測(cè)量,按條件平差法計(jì)算格網(wǎng)重力異?;驍_動(dòng)重力的改正數(shù)[10-11], 并 進(jìn) 行 精 度 評(píng) 定, 再 對(duì) 改 正 數(shù) 按Stokes或Hotine積分,求得高程異常的改正值,并最終建立韶關(guān)市區(qū)GPS水準(zhǔn)重力似大地水準(zhǔn)面,該模型分別提供30″×30″分辨率(相當(dāng)于1km×1km)似大地水準(zhǔn)面模型,見圖5。
圖5 GPS水準(zhǔn)重力似大地水準(zhǔn)面
圖6 內(nèi)符合精度檢核結(jié)果
內(nèi)符合精度就是利用計(jì)算點(diǎn)來(lái)檢驗(yàn)區(qū)域似大地水準(zhǔn)面模型的精度,內(nèi)符合精度評(píng)定過(guò)程是基于用于計(jì)算的41個(gè)GPS水準(zhǔn)點(diǎn),將其經(jīng)緯度和大地高輸入到精化后的似大地水準(zhǔn)面模型,從而得到正常高,將其給定的正常高進(jìn)行比較,最小偏離值0.1mm,最大偏離值10.2mm,平均值0.0mm,標(biāo)準(zhǔn)差為±3.9mm,見圖6。
為了進(jìn)一步量化該計(jì)算方法的可靠性以及精化后的似大地水準(zhǔn)面模型的精度,同時(shí)采用了外符合方法來(lái)加以驗(yàn)證。所謂的外符合精度驗(yàn)證即利用未參與似大地水準(zhǔn)面精化的GPS/水準(zhǔn)點(diǎn)來(lái)驗(yàn)證模型精度。本驗(yàn)證采用均勻分布于計(jì)算區(qū)域內(nèi)的其余5個(gè)GPS/水準(zhǔn)點(diǎn)作為檢核點(diǎn)對(duì)其它36個(gè)點(diǎn)精化的似大地水準(zhǔn)面模型進(jìn)行進(jìn)行外符合精度驗(yàn)證,驗(yàn)證結(jié)果見表1,其標(biāo)準(zhǔn)差為±0.70cm,綜上內(nèi)符合精度和外符合精度驗(yàn)證的結(jié)果表明:韶關(guān)市區(qū)似大地水準(zhǔn)面模型的精度達(dá)到1cm精度水平。
表1 外符合精度檢驗(yàn)
本文以韶關(guān)市基礎(chǔ)控制網(wǎng)改造項(xiàng)目為契機(jī),綜合利用測(cè)區(qū)一定密度的地面重力數(shù)據(jù),數(shù)字地面模型,以及基礎(chǔ)控制網(wǎng)改造獲得的GPS/水準(zhǔn)數(shù)據(jù)建立了韶關(guān)市區(qū)似大地水準(zhǔn)面模型,通過(guò)對(duì)模型的內(nèi)符合和外符合精度的檢驗(yàn)表明,該模型的精度在1cm精度水平。該模型投入生產(chǎn)后可替代低效率低等級(jí)精度的水準(zhǔn)測(cè)量或三角高程測(cè)量,大大降低企業(yè)生產(chǎn)成本,降低外業(yè)勞動(dòng)強(qiáng)度,提高測(cè)繪服務(wù)的信息化水平。
[1] 鄂棟臣,詹必偉,姜衛(wèi)平,等.應(yīng)用 GAMIT/GLOBK軟件進(jìn)行高精度的GPS數(shù)據(jù)處理[J].極地研究,2005,17(3):173-182.
[2] 何美琳,文鴻雁,潘元進(jìn),等.GPS高程擬合的方法比較[J].測(cè)繪科學(xué),2013,38(3):63-65.
[3] 唐秋祥,伍百發(fā).省級(jí)似大地水準(zhǔn)面的局部再精化[J].測(cè)繪科學(xué),2013,38(4):7-9.
[4] 唐秋祥,伍百發(fā).湖南省似大地水準(zhǔn)面的精度檢測(cè)和特征分析[J].測(cè)繪科學(xué),2013,38(3):10-13.
[5] 羅天文,覃鋒.GPS擬合高程建立和精化山區(qū)大地水準(zhǔn)面應(yīng)用探討[J].測(cè)繪科學(xué),2011,36(3):91-94.
[6] 陳俊勇,李建成,寧津生,等.全國(guó)及部分省市地區(qū)高精度高分辨率似大地水準(zhǔn)面的研究和實(shí)施[J].測(cè)繪通報(bào),2005(5):1-5.
[7] Consortium for Spatial Information.SRTM 90mDigital Elevation Database v4.1[EB/OL].[2014-02-20].http://www.cgiar-csi.org/data/srtm-90m-digital-elevation-database-v4-1.
[8] 張利明,李斐,章傳銀.GPS/重力邊值問(wèn)題實(shí)用公式推導(dǎo)及分析[J].地球物理學(xué)進(jìn)展,2008(6):1746-1750.
[9] GUO Jun-yi.The Basis on Physical Geodesy[M].Wuhan:Publishing House of Wuhan Science and Technical University of Surveying and Mapping,1994.
[10] 章傳銀,丁劍,晁定波,等.GPS水準(zhǔn)和重力數(shù)據(jù)聯(lián)合平差技術(shù)及其性能分析[J].測(cè)繪科學(xué),2006,31(增刊):17-19,16.
[11] 黃鵬,李靜濤.組合法精化膠州市大地水準(zhǔn)面[J].測(cè)繪科學(xué),2013,38(3):43-46.