黃宗露
(福州市勘測院有限公司,福建 福州 350000)
20世紀(jì)80年代初,我國通過天文大地網(wǎng)整體平差,建立了1980西安坐標(biāo)系[1],它的表現(xiàn)形式為二維平面坐標(biāo),是參心坐標(biāo)系。隨著科技的進(jìn)步和社會的不斷發(fā)展,2000國家大地坐標(biāo)系應(yīng)運而生。這是我國自主建立,是以全球參考基準(zhǔn)為背景的、全國統(tǒng)一的、協(xié)調(diào)一致的坐標(biāo)系統(tǒng),具有高精度、動態(tài)、適應(yīng)現(xiàn)代化空間技術(shù)發(fā)展趨勢等特點。
福清市國土已有成果大多數(shù)采用1980西安坐標(biāo)系,與國家最新要求的2000國家大地坐標(biāo)系不統(tǒng)一。為了保證坐標(biāo)轉(zhuǎn)換的精度,本文就福清市域范圍內(nèi),對均勻分布的30個控制點采用GPS-C級進(jìn)行觀測,平差后的坐標(biāo)成果用于計算福清區(qū)域內(nèi)1980西安坐標(biāo)系與2000國家大地坐標(biāo)系之間的轉(zhuǎn)換參數(shù),并分析其轉(zhuǎn)換精度和經(jīng)驗總結(jié)。
1.1.1 福州市GNSS綜合服務(wù)系統(tǒng)(FZCORS)
福州市建設(shè)好的GNSS基準(zhǔn)站有7座,等級為GPS-A級,分別為福清站、江陰站、平潭站、永泰站、連江站、軟件園站、長樂站,構(gòu)成了福州市GNSS綜合服務(wù)系統(tǒng),基本覆蓋了整個福州市。所有基站均確定了精確的2000國家大地坐標(biāo)系、1980西安坐標(biāo)系、1954北京坐標(biāo)系、WGS-84坐標(biāo)系成果,該成果已通過驗收。因此,F(xiàn)ZCORS可作為本項目GPS觀測點1980西安坐標(biāo)及2000國家大地坐標(biāo)的起算基準(zhǔn)。
1.1.2 控制點選取
控制點選取參照因素有2個:一是福清國土要求1.5個月左右的時間完成驗收并出具成果,如果重新埋設(shè)點位,工期達(dá)不到業(yè)主的需求;二是利用GPS控制點計算坐標(biāo)轉(zhuǎn)換參數(shù),GPS控制點需均勻分布,覆蓋整個測區(qū),并保證GPS測量成果時效性(或現(xiàn)勢性)。
綜上考慮后,收集了福清市國土已有點位資料,其中C級GPS點有20個,D級GPS點有106個,E級GPS點326個,經(jīng)過內(nèi)業(yè)網(wǎng)形設(shè)計和實地踏勘,選取均勻、穩(wěn)定分布的30個點位進(jìn)行觀測。
1.1.3 其他資料
收集測區(qū)已有的地形圖、影像圖和政務(wù)電子地圖等資料,可以作為網(wǎng)形設(shè)計及點之記使用,也可以充分了解測區(qū)交通線路等情況。
利用GPS技術(shù),按照整體設(shè)計、統(tǒng)一布網(wǎng)的原則,建立和維持GPS控制網(wǎng)[2],提供高精度的二維控制成果用于計算轉(zhuǎn)換參數(shù)。技術(shù)流程如圖1所示。
圖1 技術(shù)流程圖
1.3.1 外業(yè)觀測
按照GPS-C級觀測要求,本次采用Trimble R4-3型GNSS接收機進(jìn)行外業(yè)觀測,作業(yè)模式為靜態(tài)定位,與FZCORS站進(jìn)行聯(lián)測,連續(xù)觀測時間大于8小時。
1.3.2 外業(yè)數(shù)據(jù)整理
(1)天線高的整理歸類:列出并仔細(xì)檢查天線高的信息,包括天線的類型、天線高值、天線高的量測方式等[3]。
(2)RINEX格式的轉(zhuǎn)換:外業(yè)數(shù)據(jù)采集采用天寶GNSS接收機,為了方便數(shù)據(jù)處理,將接收機原始觀測數(shù)據(jù)統(tǒng)一轉(zhuǎn)換為RINEX格式[4],用軟件將觀測數(shù)據(jù)并按時間拆分成2個時段。
(3)數(shù)據(jù)歸類:以時段號建立目錄,目錄名最好為六個字符,如YYDAYS,YY為年的后兩位,DAY為年積日,S為時段號[5]。將同一年積日不同測站的數(shù)據(jù)拷貝至該目錄下,每個測站包含2個文件,一個為觀測數(shù)據(jù)O文件,一個為觀測星歷N文件。
(4)下載所需要的FZCORS站數(shù)據(jù)。
1.3.3 GPS數(shù)據(jù)處理
基線解算采用軟件是Trimble Business Center 5.0,將外業(yè)觀測數(shù)據(jù)和FZCORS站數(shù)據(jù)進(jìn)行基線解算,在基線解算完成后,根據(jù)規(guī)范相關(guān)要求,對各項精度指標(biāo)進(jìn)行評定,當(dāng)指標(biāo)都符合要求后才能進(jìn)行下一步工作。
平差采用的軟件是武漢大學(xué)編制的CosaGPS V5.21,對各觀測網(wǎng)的獨立基線向量及其全協(xié)方差矩陣進(jìn)行無約束平差和約束平差,得到30個觀測點位在1980西安坐標(biāo)系和2000國家大地坐標(biāo)系下的二維網(wǎng)平差成果,對各項指標(biāo)進(jìn)行評定,并與點位已有的坐標(biāo)值進(jìn)行比對,保證平差結(jié)果的準(zhǔn)確性。
1.4.1 轉(zhuǎn)換模型的選擇
目前通常使用的轉(zhuǎn)換模型有布爾沙模型、三維七參數(shù)轉(zhuǎn)換模型、二維七參數(shù)轉(zhuǎn)換模型、二維四參數(shù)轉(zhuǎn)換模型、三維四參數(shù)轉(zhuǎn)換模型等眾多轉(zhuǎn)換模型[6],不同的轉(zhuǎn)換模型其轉(zhuǎn)換流程、適用范圍也不盡相同。本項目將采用二維四參數(shù)轉(zhuǎn)換模型計算轉(zhuǎn)換參數(shù)。轉(zhuǎn)換參數(shù)計算流程如圖2所示。
圖2 參數(shù)計算流程
1.4.2 重合點選取原則
選取具有1980西安坐標(biāo)系和2000國家大地坐標(biāo)系這兩種坐標(biāo)系的控制點作為重合點。重合點選取的基本原則是:等級高、精度高、分布均勻、覆蓋整個轉(zhuǎn)換區(qū)域、局部變形小[7]。
1.4.3 本次計算參數(shù)的軟件
本次計算參數(shù)的軟件采用南方坐標(biāo)轉(zhuǎn)換處理軟件。
1.4.4 計算坐標(biāo)轉(zhuǎn)換參數(shù)
為了計算轉(zhuǎn)換參數(shù)的可靠性和準(zhǔn)確性,將30個重合點代入轉(zhuǎn)換模型,采用最小二乘法計算轉(zhuǎn)換參數(shù),用得到的轉(zhuǎn)換參數(shù)計算重合點坐標(biāo)殘差,剔除殘差大于3倍點位中誤差的重合點,再以剩下的點重新計算轉(zhuǎn)換參數(shù),直到所有參與計算轉(zhuǎn)換參數(shù)的點都滿足精度要求為止,然后利用未參與計算轉(zhuǎn)換參數(shù)的重合點作為外部檢核點進(jìn)行檢核,最終求得轉(zhuǎn)換參數(shù)[8]。本次項目最終從30個控制點選取22個重合點作為計算轉(zhuǎn)換參數(shù),另外8個點位作為外部檢核點。
1.5.1 內(nèi)符合精度
通過22個重合點參與四參數(shù)計算,內(nèi)符合精度見圖3。
圖3 內(nèi)符合精度圖
經(jīng)統(tǒng)計,X方向最大殘差為0.879 mm,Y方向最大殘差為0.668 mm,X方向點位中誤差為±0.365 mm,Y方向點位中誤差為±0.331 mm,平面點位中誤差為±0.493 mm,殘差均小于3倍的點位中誤差,精度良好,符合轉(zhuǎn)換要求。
1.5.2 外符合精度
利用未參與計算的8個重合點參與轉(zhuǎn)換參數(shù)檢核,外符合精度見圖4。
圖4 外符合精度圖
經(jīng)統(tǒng)計,檢核點X方向最大差值為0.384 mm,Y方向最大差值為-0.343 mm,X方向點位中誤差為±0.227 mm,Y方向點位中誤差為±0.239 mm,平面點位中誤差為±0.330 mm,精度良好,符合轉(zhuǎn)換要求。
(1)該項目工期緊,從開始到結(jié)束,只有短短一個多月的時間,利用已有點位進(jìn)行踏勘、選點,從而保證工期順利進(jìn)行。
(2)部分點位位于山上,把點位導(dǎo)入奧維軟件中,方便尋找位置,順便把軌跡記錄下來,為后面靜態(tài)觀測找點做準(zhǔn)備。
(3)由于觀測較長時間,為了防止儀器出現(xiàn)故障,電池老化續(xù)航差等情況,項目組在觀測前將儀器進(jìn)行測試,并購買了外掛電池,確保儀器能連續(xù)記錄8小時以上GPS靜態(tài)數(shù)據(jù),再投入作業(yè)中。
(4)每天觀測點位多、覆蓋范圍廣,對觀測人員、車輛、儀器等安排需做嚴(yán)密的計劃,確保觀測人員在規(guī)定時間內(nèi)到達(dá)點位上。
(5)做好數(shù)據(jù)備份,生產(chǎn)管理中應(yīng)重視過程數(shù)據(jù)的保存、留底。由于項目作業(yè)人員多、作業(yè)范圍廣,可能發(fā)生各種造成數(shù)據(jù)丟失的突發(fā)情況。
(6)內(nèi)業(yè)處理過程中部分指標(biāo)容易超限,需要對各種數(shù)據(jù)進(jìn)行分析。項目組在解算過程中咨詢了開發(fā)軟件的老師,并查閱相關(guān)規(guī)范,保證解算的各項指標(biāo)符合規(guī)范要求。
(7)在求取轉(zhuǎn)換參數(shù)時,將所有的重合點導(dǎo)入軟件中,對于不符合精度和不夠均勻的點位進(jìn)行取舍,選取精度最好的一組重合點位參與計算參數(shù)。
(8)依據(jù)求取的參數(shù),把轉(zhuǎn)換后的數(shù)據(jù)與已有的資料進(jìn)行比對,進(jìn)一步分析其精度,直到滿足所有精度后才投入使用。
福清市2000國家大地坐標(biāo)系坐標(biāo)轉(zhuǎn)換項目的實施既符合國家、福建省國土資源廳等相關(guān)政策要求,又符合福清市城市發(fā)展及業(yè)務(wù)實際需求[9]。從項目情況來看,依據(jù)最新觀測的GPS點求取坐標(biāo)轉(zhuǎn)換參數(shù),保證了現(xiàn)勢性和精度,達(dá)到了毫米級精度的轉(zhuǎn)換效果,能夠滿足實際工作的需求。該項目的順利進(jìn)行,可為同類工程的實施提供很好的思路和參考意義。