何厚學(xué)

（甘肅鐵道綜合工程勘察院有限公司，甘肅 蘭州 730000）

0.引言

全站儀測(cè)距比同層次測(cè)量?jī)x器測(cè)距更長(zhǎng)，在提升操作便利性的同時(shí)也能夠保證測(cè)量工作效率。但全站儀測(cè)量容易受到測(cè)量環(huán)境干擾，從而對(duì)全站儀測(cè)量精度產(chǎn)生影響。因此，需要對(duì)全站儀測(cè)量誤差識(shí)別，以便于后續(xù)糾正錯(cuò)誤，降低全站儀測(cè)量誤差造成的損失。

國(guó)外對(duì)于全站儀測(cè)量誤差自動(dòng)識(shí)別的研究，尚停留在對(duì)全站儀測(cè)量?jī)x器的維護(hù)上，僅就數(shù)據(jù)觀測(cè)中的誤差及環(huán)境影響做統(tǒng)計(jì)，缺少對(duì)全站儀導(dǎo)線測(cè)量的主要誤差源研究。國(guó)內(nèi)圖像識(shí)別技術(shù)較為成熟，能夠?qū)y(cè)量誤差予以精確識(shí)別分析。正因如此，國(guó)內(nèi)對(duì)全站儀測(cè)量誤差的研究范圍較廣，涉及全站儀損壞誤差、對(duì)中偏心誤差及測(cè)距誤差，其中橫軸與豎軸的垂直誤差為全站儀測(cè)量誤差識(shí)別的主要內(nèi)容，這也是目前針對(duì)全站儀測(cè)量誤差研究的主要方向：利用基坑水平位移，對(duì)全站儀極坐標(biāo)檢測(cè)，同時(shí)分析誤差并合理控制。

文獻(xiàn)[7]中，結(jié)合三角高程測(cè)量的原理對(duì)全站儀測(cè)量闡述，通過(guò)對(duì)工程的測(cè)量確定全站儀三角高程測(cè)量的過(guò)程，對(duì)比普通水準(zhǔn)測(cè)量說(shuō)明全站儀測(cè)量的優(yōu)點(diǎn)，測(cè)量應(yīng)用缺少對(duì)全站儀測(cè)量誤差的計(jì)算。

文獻(xiàn)[8]中，根據(jù)全站儀測(cè)量使用方法總結(jié)全站儀的測(cè)量原理，對(duì)全站儀測(cè)量精度的策略制定，以此提高全站儀煤礦測(cè)量工作效率，雖然全站儀的測(cè)量效率得到提高，但測(cè)量過(guò)程操作規(guī)范度不夠，導(dǎo)致測(cè)量誤差增大。

文獻(xiàn)[9]中，對(duì)全站儀的測(cè)量記錄功能說(shuō)明，對(duì)全站儀測(cè)量工程的精度單位計(jì)算，通過(guò)高斯平面校正投影邊長(zhǎng)，使得礦山測(cè)量工程進(jìn)展流程加快，保證全站儀測(cè)量工程的貫通性，但高斯平面校正投影存在大氣折射誤差。

文獻(xiàn)[10]中，對(duì)全站儀自由設(shè)站的測(cè)量工作分析，通過(guò)高精度全站儀三角高程測(cè)量，設(shè)計(jì)高程測(cè)量誤差傳播公式，根據(jù)公式計(jì)算測(cè)量誤差和大氣折光系數(shù)，以此為基礎(chǔ)提高全站儀水準(zhǔn)測(cè)量精度。但這種方法存在觀測(cè)棱鏡偏移。

本文通過(guò)樣條逼近技術(shù)布設(shè)交叉測(cè)量網(wǎng)絡(luò)，以測(cè)量網(wǎng)絡(luò)為測(cè)量平面建立坐標(biāo)系，計(jì)算全站儀測(cè)量誤差特征向量，規(guī)范測(cè)量操作過(guò)程，在計(jì)算大氣折射誤差的同時(shí)，使得全站儀觀測(cè)棱鏡始終保持在中心位置上，為避免因?qū)蜻呴L(zhǎng)相差過(guò)大引起誤差，需保證全站儀觀測(cè)導(dǎo)線邊長(zhǎng)相等。

1.基于樣條逼近的全站儀測(cè)量誤差自動(dòng)識(shí)別方法設(shè)計(jì)

1.1 計(jì)算全站儀測(cè)量誤差特征向量

對(duì)測(cè)量誤差自動(dòng)識(shí)別前，需要識(shí)別誤差特征向量，因此要對(duì)全站儀測(cè)量誤差特征向量計(jì)算。（如圖1所示）在測(cè)量物體表面設(shè)置數(shù)個(gè)固定測(cè)點(diǎn)。測(cè)點(diǎn)分別為（M1，M2，…，Mn），n≥4，設(shè)置固定測(cè)量的垂足點(diǎn)為（M1′，M2′，…，Mn′），由此求得固定測(cè)點(diǎn)到垂足點(diǎn)的距離為（M1M1′，M2M2′，…，MnMn′），設(shè)測(cè)量物體表面到垂足點(diǎn)的距離為r，由此得到該距離條件為r＝。計(jì)算距離條件為r的全站儀固定測(cè)距向量，得到測(cè)距向量的集合為，…，。因此，計(jì)算得軸線向量的垂直條件，如公式（1）所示：

圖1 全站儀測(cè)量誤差特征向量數(shù)據(jù)

據(jù)此，得到全站儀測(cè)量固定的共線條件，為軸線上四個(gè)垂足的距離相等，即。

測(cè)量固定共線條件中常見(jiàn)線性函數(shù)，設(shè)固定測(cè)點(diǎn)的觀測(cè)值誤差為（m1，m2，…，mn）。獨(dú)立觀測(cè)值與觀測(cè)函數(shù)值并不相等，由此得到誤差mz在函數(shù)z中計(jì)算公式，如公式（2）所示：

式（2）中，（k1，k2，k3，…，kn）為對(duì)應(yīng)固定測(cè)點(diǎn)誤差值（m1，m2，…，mn）的倍率系數(shù)，設(shè)全站儀獨(dú)立觀測(cè)函數(shù)z的表達(dá)式如公式（3）所示：

式（3）中，（x1，x2，…，xn）為函數(shù)總的處理觀測(cè)值，獨(dú)立觀測(cè)函數(shù)對(duì)應(yīng)的誤差值為（m1，m2，…，mk），利用真差關(guān)系處理對(duì)應(yīng)觀測(cè)值與函數(shù)值，得到函數(shù)的全微分表達(dá)式，如公式（4）所示：

對(duì)獨(dú)立觀測(cè)函數(shù)全微分表達(dá)式做真差處理，得到真差替代處理后的關(guān)系表達(dá)式如公式（5）所示：

式（5）中，獨(dú)立觀測(cè)函數(shù)對(duì)xi的偏導(dǎo)數(shù)可以表示為（i＝1，2，…，n），確定觀測(cè)值與函數(shù)式，將偏導(dǎo)數(shù)整理為常數(shù)的形式，將誤差計(jì)算公式與偏導(dǎo)數(shù)建立關(guān)系，得到公式（6）所示：

利用公式中的函數(shù)中誤差，計(jì)算全站儀測(cè)量誤差特征向量，表達(dá)式如公式（7）所示：

根據(jù)全站儀測(cè)量誤差的特征向量，分別確定對(duì)應(yīng)觀測(cè)值的誤差值，根據(jù)全站儀測(cè)量誤差特征向量，設(shè)計(jì)全站儀測(cè)量誤差自動(dòng)識(shí)別算法。

1.2.設(shè)計(jì)全站儀測(cè)量誤差自動(dòng)識(shí)別算法

利用全站儀測(cè)量誤差特征向量，制定全站儀測(cè)量誤差特征標(biāo)簽。收集全站儀測(cè)量誤差特征向量標(biāo)簽中的信息，對(duì)類(lèi)似的特征向量信息予以存儲(chǔ)，以此作為方向設(shè)計(jì)測(cè)量誤差自動(dòng)識(shí)別算法。

整合測(cè)量誤差特征向量到同一個(gè)訓(xùn)練集中，并設(shè)置訓(xùn)練集的標(biāo)簽。（如圖2所示）將全站儀測(cè)量誤差特征向量轉(zhuǎn)換為空間直角三維向量。以方差陣的形式設(shè)計(jì)每個(gè)角度方向的向量表達(dá)式B＝［β α］T，其中，T為全站儀測(cè)量訓(xùn)練集合維度，得到協(xié)方差的空間直角三維向量計(jì)算，如公式（8）所示：

圖2 坐標(biāo)轉(zhuǎn)換為空間直角三維向量

利用判斷全站儀測(cè)量誤差協(xié)方差的空間直角三維向量，確定全站儀測(cè)量誤差分量識(shí)別內(nèi)容，設(shè)全站儀測(cè)量誤差出現(xiàn)的次數(shù)為n，針對(duì)測(cè)量誤差識(shí)別特征點(diǎn)xi重復(fù)測(cè)量，并對(duì)比測(cè)量標(biāo)準(zhǔn)偏差s（xi），由此得到全站儀測(cè)量誤差的識(shí)別平均值為s（x）。設(shè)計(jì)全站儀測(cè)量誤差的識(shí)別不確定度的計(jì)算，如公式（9）所示：

式（9）中，sj為全站儀測(cè)量誤差的標(biāo)準(zhǔn)平均值；m為總測(cè)量次數(shù)，根據(jù)不確定度的來(lái)源識(shí)別合成測(cè)量誤差的輸出量，設(shè)Y為待測(cè)目標(biāo)的測(cè)量誤差結(jié)果，整合識(shí)別的測(cè)量誤差輸入量N＝［X1，X2，…，XN］，根據(jù)測(cè)量誤差識(shí)別輸入量與輸出量之間的映射關(guān)系，得到線性函數(shù)Y＝f（X1，X2，…，XN），利用線性函數(shù)整理全站儀測(cè)量誤差自動(dòng)識(shí)別輸出量與輸入量，如公式（10）所示：

式（10）中，N為線性函數(shù)整理全站儀測(cè)量誤差自動(dòng)識(shí)別輸入量個(gè)數(shù)，XN為識(shí)別測(cè)量誤差輸入量的末端值，據(jù)此得到全站儀測(cè)量誤差的自動(dòng)識(shí)別矢量公式，如公式（11）所示：

根據(jù)自動(dòng)識(shí)別矢量公式中的單元結(jié)構(gòu)參數(shù)，判斷空間姿態(tài)的識(shí)別參數(shù)分量。對(duì)比分量中的分量常數(shù)及單元方向向量，確定全站儀測(cè)量誤差自動(dòng)識(shí)別的內(nèi)容。根據(jù)全站儀測(cè)量自動(dòng)識(shí)別內(nèi)容，結(jié)合樣條逼近技術(shù)構(gòu)建自動(dòng)識(shí)別模型，通過(guò)全站儀測(cè)量誤差的自動(dòng)識(shí)別相關(guān)性，對(duì)全站儀空間分量的測(cè)量誤差識(shí)別。

1.3 構(gòu)建樣條逼近全站儀測(cè)量誤差自動(dòng)識(shí)別模型

根據(jù)全站儀測(cè)量誤差自動(dòng)識(shí)別算法，分配樣條逼近全站儀測(cè)量誤差自動(dòng)識(shí)別的分量，據(jù)此給出相互獨(dú)立的測(cè)量誤差自動(dòng)識(shí)別內(nèi)容，通過(guò)輸入映射核函數(shù)空間特征，實(shí)現(xiàn)非線性間隔分類(lèi)。

針對(duì)全站儀測(cè)量誤差自動(dòng)識(shí)別算法，設(shè)計(jì)測(cè)量誤差損失函數(shù)，根據(jù)選擇的不平衡測(cè)量參數(shù)，區(qū)分小類(lèi)樣本測(cè)量誤差識(shí)別與訓(xùn)練集測(cè)量誤差識(shí)別。構(gòu)建全站儀測(cè)量最優(yōu)參數(shù)集合，并據(jù)此設(shè)計(jì)全站儀測(cè)量誤差識(shí)別增長(zhǎng)序列g(shù)＝2－15，2－13，…，23，設(shè)為全站儀測(cè)量誤差自動(dòng)識(shí)別修正值，則識(shí)別誤差觀測(cè)點(diǎn)間的距離如公式（12）所示：

式（12）中，K為樣條逼近測(cè)量誤差的修正集合，即K＝［K1，K2，K3］，轉(zhuǎn)換得到識(shí)別誤差修正集合的常數(shù)矩陣，利用修正集合中的常數(shù)項(xiàng)，對(duì)誤差觀點(diǎn)間距離計(jì)算識(shí)別向量的方差陣予以表達(dá)，如公式（13）所示：

設(shè)全站儀測(cè)量誤差自動(dòng)識(shí)別中的測(cè)點(diǎn)坐標(biāo)為T(mén)＝［x0，y0，z0］，根據(jù)誤差識(shí)別測(cè)點(diǎn)坐標(biāo)表達(dá)方差陣方向向量，利用測(cè)得的全站儀向量方差陣對(duì)誤差特征予以表達(dá)。設(shè)定特定的全站儀測(cè)量誤差自動(dòng)識(shí)別閾值，以閾值為限確定樣條逼近的全站儀測(cè)量誤差的識(shí)別結(jié)果，判斷全站儀測(cè)量誤差的識(shí)別擬合結(jié)果，對(duì)向量中的擬合維度判斷。利用曲線逼近的辨別臨界滑動(dòng)面，判斷全站儀測(cè)量誤差自動(dòng)識(shí)別的特征誤差擬合，以此為基礎(chǔ)對(duì)樣條曲線重構(gòu)，壓縮并控制觀測(cè)點(diǎn)。利用漸進(jìn)迭代逼近的方法，對(duì)測(cè)量擴(kuò)展曲線識(shí)別。以此完成對(duì)樣條逼近全站儀測(cè)量誤差自動(dòng)識(shí)別模型的構(gòu)建。

利用構(gòu)建的樣條逼近全站儀測(cè)量誤差自動(dòng)識(shí)別模型，確定全站儀測(cè)量誤差自動(dòng)識(shí)別的向量維度，根據(jù)維度中測(cè)量誤差的特征，對(duì)應(yīng)自動(dòng)識(shí)別的測(cè)量誤差指標(biāo)，在確定模型的F1-measure小于常規(guī)值時(shí)，將測(cè)量誤差的向量維度導(dǎo)入自動(dòng)識(shí)別資料庫(kù)中。根據(jù)資料庫(kù)中全站儀測(cè)量誤差維度，確定自動(dòng)識(shí)別的填補(bǔ)項(xiàng)。并針對(duì)全站儀測(cè)量誤差自動(dòng)識(shí)別訓(xùn)練集，對(duì)全站儀測(cè)量誤差特征分類(lèi)，以確定自動(dòng)識(shí)別對(duì)應(yīng)的特征向量?jī)?nèi)容。對(duì)比特征向量?jī)?nèi)容與預(yù)定閾值大小，完成對(duì)全站儀測(cè)量誤差的自動(dòng)識(shí)別。

1.4 完成全站儀測(cè)量誤差自動(dòng)識(shí)別

利用全站儀測(cè)量誤差自動(dòng)識(shí)別模型，對(duì)全站儀測(cè)量誤差自動(dòng)識(shí)別。構(gòu)建不同種類(lèi)的全站儀測(cè)量誤差識(shí)別標(biāo)準(zhǔn)，考慮到全站儀的測(cè)量視場(chǎng)范圍較大，要確定角度儀器的修正范圍，設(shè)定位瞄準(zhǔn)的相關(guān)重復(fù)性誤差系數(shù)為0.4，全站儀的自準(zhǔn)分辨率在以0.01＂下，設(shè)全站儀測(cè)量誤差自動(dòng)識(shí)別中線坐標(biāo)軌跡的計(jì)算，如公式（14）所示：

式（14）中，軌跡檢測(cè)誤差的測(cè)量值為g，則水平傾角在軌跡中線以外的測(cè)量誤差識(shí)別，需要以全站儀測(cè)量誤差獨(dú)立坐標(biāo)原點(diǎn)為中心，做坐標(biāo)轉(zhuǎn)換處理，如公式（15）所示：

式（15）中，R為全站儀測(cè)量誤差旋轉(zhuǎn)矩陣，據(jù)測(cè)判斷全站儀測(cè)量誤差自動(dòng)識(shí)別后的內(nèi)積結(jié)果Z＝XiXc＋YiYc＋ZiZc。根據(jù)測(cè)量誤差的自動(dòng)識(shí)別內(nèi)積結(jié)果設(shè)定測(cè)量誤差外積處理內(nèi)容，針對(duì)自動(dòng)識(shí)別的內(nèi)積設(shè)定結(jié)果計(jì)算測(cè)量誤差自動(dòng)識(shí)別的不確定度。

將全站儀測(cè)量誤差自動(dòng)識(shí)別的算法轉(zhuǎn)換為矩陣乘法，根據(jù)全站儀視準(zhǔn)軸的空間位姿確定不確定度的單元?jiǎng)討B(tài)位姿，確定全站儀測(cè)量誤差初始向量水平軸方向。得到樣條逼近的全站儀測(cè)量誤差自動(dòng)識(shí)別不確定度的公式，如公式（16）所示：

式（16）中，?為全站儀測(cè)量誤差的包含因子；θ為測(cè)量誤差擴(kuò)展角度，利用樣條逼近的全站儀測(cè)量誤差自動(dòng)識(shí)別的不確定度，對(duì)樣條逼近的全站儀測(cè)量誤差計(jì)算，以此確定測(cè)量誤差識(shí)別的矢量大小。根據(jù)置信概率重新構(gòu)建自動(dòng)識(shí)別誤差矢量矩陣，得到自動(dòng)識(shí)別測(cè)量誤差的動(dòng)態(tài)坐標(biāo)（如圖3所示）：

圖3 自動(dòng)識(shí)別測(cè)量誤差的動(dòng)態(tài)坐標(biāo)

圖3中，按照點(diǎn)號(hào)的順序，對(duì)全站儀測(cè)量誤差自動(dòng)識(shí)別，識(shí)別結(jié)果自動(dòng)生成折線，根據(jù)測(cè)量結(jié)果的偏差，對(duì)全站儀的測(cè)量誤差計(jì)算。至此，完成對(duì)基于樣條逼近的全站儀測(cè)量誤差自動(dòng)識(shí)別方法設(shè)計(jì)。

2.實(shí)驗(yàn)分析

為提高誤差識(shí)別的速率，設(shè)計(jì)對(duì)比實(shí)驗(yàn)，對(duì)比數(shù)據(jù)分析識(shí)別法、二值化檢測(cè)識(shí)別法、基于樣條逼近的全站儀測(cè)量誤差自動(dòng)識(shí)別方法的測(cè)量誤差識(shí)別SVN（SVN全名Subversion，即版本控制系統(tǒng)）檢出量。

2.1.準(zhǔn)備過(guò)程

應(yīng)用全站儀對(duì)某目標(biāo)進(jìn)行測(cè)量，對(duì)全站儀測(cè)量現(xiàn)場(chǎng)中自動(dòng)識(shí)別需要的基礎(chǔ)數(shù)值（如圖4所示）進(jìn)行測(cè)量。

圖4 全站儀自動(dòng)識(shí)別基礎(chǔ)數(shù)據(jù)采集示意圖

處理圖中的基礎(chǔ)數(shù)據(jù)，根據(jù)觀測(cè)目標(biāo)到全站儀的距離D，確定觀測(cè)目標(biāo)的實(shí)際測(cè)量測(cè)距，E為全站儀到觀測(cè)目標(biāo)的實(shí)際距離；yp為觀測(cè)目標(biāo)到測(cè)量平面y的垂足；xp為觀測(cè)目標(biāo)到測(cè)量平面x的垂足；v和Hz分別為全站儀測(cè)量角度和觀測(cè)目標(biāo)；0為全站儀測(cè)量所在的坐標(biāo)原點(diǎn)；P為觀測(cè)目標(biāo)所在的點(diǎn)；th為觀測(cè)目標(biāo)距離最短平面的距離；ih為全站儀觀測(cè)平面到地面的距離。

整理好基礎(chǔ)數(shù)據(jù)后，移動(dòng)觀測(cè)目標(biāo)，并對(duì)上述基礎(chǔ)數(shù)值進(jìn)行移動(dòng)后的測(cè)量，針對(duì)全站儀移動(dòng)線路計(jì)算，得到全站儀測(cè)量誤差自動(dòng)識(shí)別的線路偏差結(jié)果（如圖5所示）：

圖5 全站儀測(cè)量誤差自動(dòng)識(shí)別的線路偏差結(jié)果

分析圖中線路偏差結(jié)果可知：水平誤差與垂直誤差的線路偏差結(jié)果相近。由此可知線路測(cè)量選點(diǎn)操作與誤差檢出的線路偏差波動(dòng)集中在3.75mm～7.50mm之間，為了保證線路偏差結(jié)果更接近測(cè)量誤差標(biāo)準(zhǔn)，選擇角度中誤差在0.4°～0.6°的線路偏差區(qū)間作為實(shí)驗(yàn)測(cè)量?jī)?nèi)容。線路偏差區(qū)間中的測(cè)量誤差數(shù)據(jù)分為10組。通過(guò)整理全站儀測(cè)量的自動(dòng)識(shí)別測(cè)量結(jié)果，根據(jù)測(cè)量選點(diǎn)確定偏差波動(dòng)的數(shù)據(jù)特征，針對(duì)上述數(shù)據(jù)特征做自動(dòng)識(shí)別監(jiān)測(cè)。針對(duì)全站儀測(cè)量誤差自動(dòng)識(shí)別中的水平誤差與垂直誤差做區(qū)分處理，分別計(jì)算線路偏差與角度中誤差，由此整理得到十組測(cè)量誤差識(shí)別SVN檢出量。

2.2 分析結(jié)果

整理得到全站儀測(cè)量識(shí)別誤差SVN檢出量（如圖6所示）：

圖6 測(cè)量誤差自動(dòng)識(shí)別前后的誤差SVN檢出量

由圖6可知：數(shù)據(jù)分析識(shí)別法的全站儀測(cè)量識(shí)別誤差SVN檢出量在0.60以下，由于檢出量受測(cè)量次數(shù)干擾逐漸趨于平穩(wěn)，因此在后續(xù)的7～10組測(cè)量檢出無(wú)明顯變化，對(duì)后續(xù)測(cè)量結(jié)果的誤差檢出始終保持不變。可見(jiàn)應(yīng)用數(shù)據(jù)分析識(shí)別法對(duì)全站測(cè)量誤差的檢出效果并不明顯。

應(yīng)用二值化檢測(cè)識(shí)別法的全站儀測(cè)量識(shí)別誤差SVN檢出量在0.50以下，測(cè)量檢出率在達(dá)到第7組時(shí)趨于平穩(wěn)，整理測(cè)量識(shí)別檢出率低于數(shù)據(jù)分析識(shí)別法。

應(yīng)用基于樣條逼近的全站測(cè)量誤差自動(dòng)識(shí)別方法，全站測(cè)量誤差自動(dòng)識(shí)別前后的誤差SVN檢出量在前幾組保持穩(wěn)定，處于0.60左右，在第5組測(cè)量檢出結(jié)果后，誤差SVN檢出量明顯上升，在第10組達(dá)到檢出最高為0.88。因此，基于樣條逼近的全站儀測(cè)量誤差自動(dòng)識(shí)別方法更優(yōu)異。

3.結(jié)束語(yǔ)

通過(guò)本文研究，提高了全站儀測(cè)量誤差自動(dòng)識(shí)別的檢出量。構(gòu)建樣條逼近全站儀測(cè)量誤差自動(dòng)識(shí)別模型，為全站儀測(cè)量誤差自動(dòng)識(shí)別提供基礎(chǔ)數(shù)據(jù)。今后應(yīng)繼續(xù)研究提高全站儀測(cè)量精度，制定全站儀的測(cè)量識(shí)別指標(biāo)，以此為基礎(chǔ)建立基坑水平位移測(cè)量模型，利用測(cè)繪的空間地理信息提高全站儀的利用率。