陳佳清,劉立龍,黎峻宇,劉瑞駒,劉家智
(1.欽州市住房和城鄉(xiāng)建設(shè)測(cè)繪院,廣西 欽州 535000;2.桂林理工大學(xué) 測(cè)繪地理信息學(xué)院,桂林 541006;3.廣西空間信息與測(cè)繪重點(diǎn)實(shí)驗(yàn)室,桂林 541006)
全球衛(wèi)星導(dǎo)航系統(tǒng)(global navigation satellite system,GNSS)已經(jīng)在各行各業(yè)中得到了廣泛的應(yīng)用,國(guó)際GNSS服務(wù)組織(International GNSS Service,IGS)也對(duì)各類產(chǎn)品解算不斷進(jìn)行改進(jìn)和優(yōu)化,精密單點(diǎn)定位技術(shù)(precise point positioning,PPP)研究進(jìn)入了一個(gè)新階段[1]。國(guó)內(nèi)外科研院校、生產(chǎn)部門的學(xué)者相繼對(duì)精密單點(diǎn)定位技術(shù)進(jìn)行了大量理論研究和技術(shù)創(chuàng)新[2-8],精密單點(diǎn)定位數(shù)據(jù)處理軟件開發(fā)技術(shù)越來(lái)越成熟。國(guó)外有加拿大Waypoint公司推出的GrafNav軟件、APPLA-Ni X公司的POSPac AIR、挪威Terra Tec公司的Terra POS軟件、瑞士Leica公司的IPAS等,但是這些都是商用軟件,價(jià)格較為昂貴、不提供源代碼,且數(shù)據(jù)處理過程較為繁雜。國(guó)內(nèi)有武漢大學(xué)研發(fā)的精密單點(diǎn)定位軟件TriP,其源代碼不公開[1]。
RTKLIB是日本東京海洋大學(xué)(Tokyo University of Marine Science and Technology)開發(fā)的一個(gè)用于GNSS精確定位的開源軟件[9],與其他商用軟件相比,其最大的優(yōu)勢(shì)在于免費(fèi),且所有源代碼公開,便于學(xué)習(xí)研究和算法測(cè)試[10]。為此,國(guó)內(nèi)外學(xué)者陸續(xù)對(duì)其性能進(jìn)行了大量的相關(guān)試驗(yàn)及應(yīng)用研究。國(guó)外關(guān)于RTKLIB的研究起步較早,并已成功應(yīng)用在相關(guān)領(lǐng)域:文獻(xiàn)[11]中利用低成本的全球定位系統(tǒng)(global positioning system,GPS)模塊u-blox和天線,借助RTKLIB軟件處理獲得了厘米級(jí)的靜態(tài)定位;文獻(xiàn)[12]評(píng)估了RTKLIB的準(zhǔn)確性,成功達(dá)到了安裝在飛機(jī)上的合成孔徑雷達(dá)(SAR)厘米級(jí)的精度要求;文獻(xiàn)[13]利用RTKLIB軟件處理GNSS數(shù)據(jù)測(cè)定冰川坐標(biāo)及流速。國(guó)內(nèi)方面:文獻(xiàn)[10]研究利用RTKLIB軟件精密單點(diǎn)定位精度達(dá)厘米級(jí);文獻(xiàn)[14]利用RTKLIB軟件分析了不同星歷、鐘差對(duì)定位精度的影響,得出事后星歷和快速星歷單天解3維坐標(biāo)分量得收斂值方向偏差在4 cm左右的結(jié)論;文獻(xiàn)[15]對(duì)比了BERNESE5.2和RTKLIB軟件靜態(tài)單點(diǎn)定位解算精度。已有的文獻(xiàn)在研究RTKLIB靜態(tài)精密單點(diǎn)定位過程中,均采用前向?yàn)V波方法進(jìn)行數(shù)據(jù)處理,而RTKLIB軟件提供了多種解算方式[16],因此對(duì)其應(yīng)用研究尚不夠完善和全面。
本文對(duì)比分析RTKLIB在不同星歷產(chǎn)品、不同濾波方法和不同觀測(cè)時(shí)間下的解算精度,探求采用精密星歷結(jié)合組合濾波方法的更優(yōu)解算策略。
GPS單點(diǎn)定位采用載波相位與偽距無(wú)電離層組合[10](Ionosphere-free LC)觀測(cè)方程[5],觀測(cè)方程簡(jiǎn)化為:
(1)
(2)
RTKLIB軟件精密單點(diǎn)定位中采用的數(shù)據(jù)處理算法為擴(kuò)展卡爾曼濾波(extended Kalman filter,EKF)[16]。運(yùn)用卡爾曼濾波算法,通過tk時(shí)刻系統(tǒng)觀測(cè)向量yk的更新估算未知參數(shù)狀態(tài)向量x及其協(xié)方差矩陣P,其濾波方程為:
(3)
(4)
(5)
(6)
RTKLIB在靜態(tài)精密單點(diǎn)定位中采用EKF對(duì)GNSS測(cè)量信號(hào)、對(duì)流層、電離層等進(jìn)行濾波處理。對(duì)于事后PPP,其濾波類型有前向、后向、組合3種方法[16]。
從RTKLIB源代碼可以看到,在RTKLIB預(yù)處理模塊中采用的是電離殘差(L1-L2 geometry free)周跳探測(cè)法,該方法只對(duì)周跳進(jìn)行探測(cè)不進(jìn)行周跳修復(fù)[10]。為獲得精密單點(diǎn)解算最優(yōu)精度進(jìn)行的誤差改正有:1)電離層時(shí)延改正,為了消除GNSS信號(hào)測(cè)量中的電離層效應(yīng),采用無(wú)電離層組合模型進(jìn)行GNSS數(shù)據(jù)處理;2)對(duì)流層時(shí)延改正,采用精密對(duì)流層模型“Estimate ZTD”估計(jì)天頂總對(duì)流層延遲改正;3)星歷、鐘差改正采用國(guó)際IGS中心提供的事后精密星歷、鐘差文件進(jìn)行改正。
算例采用BJFS、CHAN、SHAO 3個(gè)IGS站1 877 GPS星期第5天的數(shù)據(jù),采樣間隔為30 s。數(shù)據(jù)處理之前首先下載該天3個(gè)IGS站的觀測(cè)文件和導(dǎo)航文件(BJFS18770010.16o、BJFS18770010.16n,CHAN18770010.16o、CHAN18770010.16n,SHAO18770010.16o、SHAO18770010.16n),并從ftp://garner.ucsd.edu/網(wǎng)站下載實(shí)驗(yàn)所需sp3格式的IGS精密星歷數(shù)據(jù):igs18775.sp3、igr18775.sp3、igu18775.sp3以及精密鐘差文件igs18775.clk。
在靜態(tài)單點(diǎn)定位解算之前利用RTKLIB軟件對(duì)數(shù)據(jù)質(zhì)量進(jìn)行分析,必要的預(yù)處理是提高解算精度的前提。繪制3個(gè)IGS站衛(wèi)星分布圖(如圖1所示),高度角、多路徑誤差、信噪比如圖2所示。
從圖1可以看出:3個(gè)IGS站的衛(wèi)星分布狀態(tài)良好,BJFS站衛(wèi)星可視化程度較另外2站稍差,初步推斷可能受建筑物或其他物體遮擋。從圖2可以看出:3個(gè)IGS站的衛(wèi)星高度角基本保持一致;多路徑誤差隨觀測(cè)時(shí)間的增加而變小,BJFS站的多路徑誤差均方根(root mean square,RMS)值為0.394 7 m,CHAN站的多路徑誤差RMS值為0.600 7 m,SHAO站的多路徑誤差RMS值為0.475 7 m,多路徑誤差最大為CHAN站;信噪比隨觀測(cè)時(shí)間的增加而增強(qiáng),CHAN、SHAO 2個(gè)IGS站的信噪比SNR在40 至60 dB之間。
根據(jù)上述IGS各站觀測(cè)數(shù)據(jù)質(zhì)量分析結(jié)果初步得出SHAO站觀測(cè)質(zhì)量最優(yōu)。
經(jīng)RTKLIB質(zhì)量分析后,分別對(duì)3個(gè)IGS站進(jìn)行靜態(tài)單點(diǎn)定位解算。解算過程中統(tǒng)一采用精密鐘差產(chǎn)品,對(duì)比不同星歷產(chǎn)品、不同濾波方式、不同觀測(cè)時(shí)間下的解算精度并進(jìn)行統(tǒng)計(jì)分析。IGS各站點(diǎn)東向(E)、北向(N)、天向(U)收斂偏差值統(tǒng)計(jì)如表1所示。
表1 IGS各站點(diǎn)E、N、U方向收斂偏差值 cm
注:前向、后向收斂值為最終收斂值,組合收斂值為單天區(qū)間內(nèi)收斂值;字體加粗表示E方向最小值,字體加粗加斜表示N方向最小值,字體加粗加下劃線表示U方向最小值。
從表1中可以看出:1)采用igs和igr星歷解算單天IGS各站點(diǎn)E、N、U方向收斂值基本一致,優(yōu)于igu星歷,N方向收斂偏差值小于E、U方向;2)采用組合濾波方法優(yōu)于前向和后向?yàn)V波方法,BJFS、CHAN、SHAO 3個(gè)IGS站的3個(gè)方向收斂偏差值均不大于1 cm;3)采用igs星歷結(jié)合組合濾波方法后,收斂偏差最小值E方向?yàn)?.5 cm、N方向?yàn)?.2 cm、U方向?yàn)?.9 cm;4)SHAO站收斂偏差結(jié)果優(yōu)于另外2個(gè)IGS站。IGS各站點(diǎn)E、N、U方向及RMS值統(tǒng)計(jì)如表2所示。
表2 IGS各站點(diǎn)E、N、U方向及平面RMS值 cm
注:字體加粗表示E方向RMS最小值,字體加粗加斜表示N方向RMS最小值,字體加粗加下劃線表示U方向RMS最小值,字體加粗加雙下劃線表示平面點(diǎn)位中誤差最小值。
從表2中可以看出:1)采用igs星歷組合濾波方法,IGS各站點(diǎn)E、N、U方向RMS值及平面中誤差均小于0.6 cm,RMS最小值N方向?yàn)?.06 cm、E方向?yàn)?.1 cm、U方向?yàn)?.29 cm,點(diǎn)位中誤差最小值為0.29 cm;2)3個(gè)IGS站解算結(jié)果均獲得毫米級(jí)精度,解算最優(yōu)結(jié)果的IGS站點(diǎn)為SHAO站,點(diǎn)位精度為0.29 cm;3)解算結(jié)果表明RTKLIB靜態(tài)單點(diǎn)定精度達(dá)到毫米級(jí)。IGS各站的收斂偏差值(事后精密星歷結(jié)合組合濾波法)如圖3~圖5所示。
由圖3可知:BJFS站的E方向偏差收斂不穩(wěn)定,U方向穩(wěn)定性較E方向差,N方向整體平穩(wěn),無(wú)明顯波動(dòng);E、U方向波動(dòng)幅度在0.9 cm內(nèi),N方向波動(dòng)幅度在0.3 cm內(nèi);E、N、U方向RMS值分別為0.47、0.09及0.33 cm。
由圖4可知:CHAN站的E方向偏差收斂不穩(wěn)定,U方向穩(wěn)定性較E方向差,N方向整體平穩(wěn),無(wú)明顯波動(dòng);E、U方向波動(dòng)幅度在1 cm內(nèi),N方向波動(dòng)幅度在0.3 cm內(nèi);E、N、U方向RMS值分別為0.54、0.11及0.33 cm。
由圖5可知:SHAO站的E方向偏差收斂不穩(wěn)定,U方向穩(wěn)定性較E方向差,N方向整體平穩(wěn),無(wú)明顯波動(dòng);E、U方向波動(dòng)幅度在0.9 cm內(nèi),N方向波動(dòng)幅度在0.2 cm內(nèi);E、N、U方向RMS值分別為0.29、0.06及0.29 cm。
對(duì)比圖3~圖5可知:BJFS、CHAN、SHAO 3個(gè)IGS站的E、N、U方向收斂值偏差基本一致,E方向收斂偏差值不穩(wěn)定,但最大偏差值不超過1 cm;N方向收斂偏差值較為平穩(wěn),在0.4 cm以內(nèi);U方向收斂偏差值波動(dòng)較大,但也不超過1 cm。以SHAO站為例,不同觀測(cè)時(shí)間對(duì)解算精度的影響如圖6所示。
由圖6可知:E、N、U方向精度曲線基本一致,精度隨著觀測(cè)時(shí)間的遞增而提高;N、U方向在觀測(cè)24 h時(shí)精度最高,RMS最小值分別為0.06和0.29 cm,略優(yōu)于文獻(xiàn)[2]中的0.1和0.3 cm;E方向在觀測(cè)時(shí)段14時(shí)至20時(shí)由于受到觀測(cè)粗差的影響,精度隨時(shí)間遞增而減??;受到E方向精度的影響,平面精度出現(xiàn)同樣的影響;E方向及平面在觀測(cè)14 h時(shí)精度最高,RMS最小值分別為0.19和0.23 cm;觀測(cè)4 h左右時(shí)的平面精度達(dá)厘米級(jí),在觀測(cè)條件良好條件下,單天解算精度可達(dá)毫米級(jí)。
綜合上述靜態(tài)單點(diǎn)定位精度分析結(jié)果,對(duì)于高精度單點(diǎn)定位解算可參考以下靜態(tài)精密單點(diǎn)定位更優(yōu)解算策略,如表3所示。
表3 靜態(tài)精密單點(diǎn)定位更優(yōu)解算策略
注:PCO全稱為phase center offset,指接收機(jī)天線相位中心改正;PCV全稱為phase center variation,指相位中心偏差;I08.ATX為改正文件。
RTKLIB軟件靜態(tài)精密單點(diǎn)定位解算結(jié)果表明:igs星歷和igr星歷解算精度相當(dāng),優(yōu)于igu星歷;前向?yàn)V波和后向?yàn)V波方法解算精度略有所差異但并不明顯,組合濾波方法相對(duì)于其他2種方法能夠提高10倍以上精度;解算精度隨觀測(cè)時(shí)間的增加而提高,同時(shí)受到觀測(cè)過程中誤差的影響而出現(xiàn)突變。采用事后精密星歷結(jié)合組合濾波方法,3個(gè)IGS站單天解E、N、U3個(gè)方向偏差收斂值均在1 cm內(nèi),E、N、U方向及平面定位精度達(dá)到毫米級(jí)。其中,SHAO站的靜態(tài)精密定位精度為0.29 cm,優(yōu)于其他2站,驗(yàn)證了RTKLIB數(shù)據(jù)質(zhì)量分析中SHAO站觀測(cè)質(zhì)量最優(yōu)的結(jié)論。在靜態(tài)精密單點(diǎn)定位中,RTKLIB具有解算精度高、效率高、可靠性高的特點(diǎn),對(duì)目前的地殼形變監(jiān)測(cè)及工程應(yīng)用具有重要的現(xiàn)實(shí)意義。
[1] 張小紅,何錫揚(yáng),李星星.TriP軟件非差幾何法精密定軌精度分析[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2010,35(11):1327-1330.
[2] 張小紅,左翔,李盼,等.BDS/GPS精密單點(diǎn)定位收斂時(shí)間與定位精度的比較[J].測(cè)繪學(xué)報(bào),2015,44(3):250-256.
[3] 賴允斌,趙春梅,李子申.不同星歷下實(shí)時(shí)精密單點(diǎn)定位精度分析[J].測(cè)繪通報(bào),2015(8):9-12.
[4] 周磊.GNSS精密單點(diǎn)定位理論及其軟件研制[D].西安:西安科技大學(xué),2015.
[5] 臧楠.BDS/GNSS精密單點(diǎn)定位算法研究[D].西安:長(zhǎng)安大學(xué),2015.
[6] OLIVEIRA P S,MOREL L,F(xiàn)UND F,et al.Modeling tropospheric wet delays with dense and sparse network configurations for PPP-RTK[J].GPS Solutions,2017,21(1):237-250.
[7] GANDOLFI S,TAVASCI L,POLUZZI L.Study on GPS-PPP precision for short observation sessions[J].GPS Solutions,2016:1-10.
[8] EL-MOWAFY A,DEO M,RIZOS C.On biases in precise point positioning with multi-constellation and multi-frequency GNSS data[EB/OL].(2016-12-09)[2017-02-18].https://www.researchgate.net/profile/Ahmed_El-Mowafy/publication/291384900_On_biases_in_precise_point_positioning_with_multi-constellation_and_multi-frequency_GNSS_data/links/5705d3c208ae13eb88b97d16.pdf?origin=publication_detail.
[9] Tokyo University of Marine Science and Technology Laboratory of Communications Engineering.RTKLIB:An open source program package for GNSS positioning[EB/OL].(2016-12-19)[2017-02-18].http://www.rtklib.com/rtklib.htm.
[10] 丁慧君,羅端,盧兵.基于RTKLIB的精密單點(diǎn)定位研究[J].城市勘測(cè),2014,10(5):22-25.
[11] WISNIEWSKI B,BRUNIECKI K,MOSZYNSKI M.Evaluation of RTKLIB’s positioning accuracy using low-cost GNSS receiver and ASG-EUPOS[J].Transnav International Journal on Marine Navigation & Safety of Sea Transportation,2013,7(2):79-85.
[12] VIDEKULL R.Evaluate and develop high-performance GPS navigation using free GPS software[D].G?teborg,Sweden:Chalmers University of Technology,2015.
[13] SUGIYAMA S,SAKAKIBARA D,TSUTAKI S,et al.Glacier dynamics near the calving front of Bowdoin Glacier,northwestern Greenland[J].Journal of Glaciology,2015,61(226):223-232.
[14] 張建龍,徐愛功,張兆南,等.RTKLIB軟件靜態(tài)精密單點(diǎn)定位精度測(cè)試與分析[J].全球定位系統(tǒng),2014,39(1):37-41.
[15] 汪洋.GNSS精密單點(diǎn)定位數(shù)據(jù)解算結(jié)果的比較和分析[C]//中國(guó)衛(wèi)星導(dǎo)航系統(tǒng)管理辦公室學(xué)術(shù)交流中心.第七屆中國(guó)衛(wèi)星導(dǎo)航學(xué)術(shù)年會(huì)論文集.北京:中國(guó)衛(wèi)星導(dǎo)航系統(tǒng)管理辦公室,2016:1-4.
[16] TAKASU 1 T,YASUDA A.RTKLIB ver.2.4.2 manual[EB/OL].(2013-04-29)[2017-02-18].http://www.rtklib.com/prog/manual_2.4.2.pdf.