郭獻(xiàn)濤，賈 燕

（南京郵電大學(xué) 地理與生物信息學(xué)院，江蘇 南京210023）

1 引言

地面激光掃描（Terrestrial Laser Scanning，TLS）技術(shù)又稱實(shí)景復(fù)制技術(shù)，它是測(cè)繪領(lǐng)域中繼GPS之后的又一次技術(shù)革命。該技術(shù)通過極坐標(biāo)法掃描測(cè)量激光到對(duì)象表面的距離、水平角及垂直角，可以快速、高精度、免棱鏡式地獲取對(duì)象表面的高密度三維數(shù)據(jù)，得到對(duì)象的離散化表示；但TLS測(cè)量精度較大程度上受儀器系統(tǒng)誤差影響，為保證數(shù)據(jù)質(zhì)量，須在實(shí)際掃描測(cè)量前對(duì)儀器實(shí)施檢校［1］。

目前，關(guān)于TLS儀器的檢校還沒有通用的方法與評(píng)價(jià)體系，主要有分項(xiàng)檢校和整體檢校兩種模式。分項(xiàng)檢校主要采用誤差分析法對(duì)測(cè)距測(cè)角等模塊實(shí)施單獨(dú)檢校，需精確建立儀器的誤差模型，受掃描儀專利設(shè)計(jì)影響，這方面的知識(shí)通常較為有限［2］。此外，該類檢校還需檢?；€等特殊檢校設(shè)施，因此難以普及。整體檢校是將掃描儀所有組件作為一個(gè)整體，直接與標(biāo)準(zhǔn)檢校場(chǎng)比較測(cè)量來獲取系統(tǒng)改正數(shù)。該類檢校無(wú)需精確已知儀器的誤差模型，檢校場(chǎng)建立也不復(fù)雜，目前可通過自檢校實(shí)施整體檢校［3］。

近年來，自檢校已逐漸成為TLS檢校的流行方法，它可以確定掃描儀所有系統(tǒng)誤差及其他系統(tǒng)參數(shù)?；谑褂玫臋z校標(biāo)志類型，TLS自檢?？煞譃閮煞N類型［4］：（1）基于點(diǎn)的自檢校，即利用可提取中心坐標(biāo)的標(biāo)志實(shí)施檢校。Lichti給出了一個(gè)自由網(wǎng)平差法，以同時(shí)估計(jì)外方位參數(shù)、對(duì)象點(diǎn)坐標(biāo)與檢校參數(shù)［5］，檢校了Faro 880、Trim?ble GS200、Trimble GX等掃描儀，并分析了檢校參數(shù)的時(shí)序行為，文獻(xiàn)［6］詳細(xì)描述了上述數(shù)學(xué)模型。Krzysztof對(duì)掃描儀FARO 880實(shí)施了一系列自檢校實(shí)驗(yàn)，較好地識(shí)別了TLS儀器誤差［7］。此外，Soudarissanane等人給出了一個(gè)全景相機(jī)擴(kuò)展傳感器模型作為TLS誤差模型，使用掃描強(qiáng)度圖像代替點(diǎn)云對(duì)激光掃描儀Imager 5003實(shí)施建模和檢校［8］。目前這些研究主要集中于相位式掃描儀，且檢校參數(shù)間存在高相關(guān)性導(dǎo)致參數(shù)估計(jì)不確定性的問題。（2）基于平面的自檢校，即利用平面模板反射坐標(biāo)實(shí)施檢校，該方法由Gielsdorf提出，并被應(yīng)用于掃描儀Imager 5003與PoMeS的實(shí)驗(yàn)室檢校［9］。Chow改進(jìn)了該方法，利用室內(nèi)墻壁、地板和天花板代替專門的檢校面板實(shí)施了掃描儀FARO880的自檢校，改善了自檢校殘差標(biāo)準(zhǔn)差的觀測(cè)精度。Dorninger等人使用平面特征進(jìn)行自校準(zhǔn)，主要專注于現(xiàn)場(chǎng)閉合差估計(jì)［10］。此外，Mahbaz建議采用連續(xù)分段線性改正函數(shù)替代諧波函數(shù)和多項(xiàng)式函數(shù)［11］，該改正函數(shù)不需要對(duì)檢校模型做假設(shè)，但需要更多的參數(shù)。相比于點(diǎn)檢校法，基于平面的自檢校法無(wú)需做大量的單點(diǎn)測(cè)量，但結(jié)果卻與基于點(diǎn)自檢校方法相差較大，使系統(tǒng)的誤差識(shí)別復(fù)雜化。

本文基于點(diǎn)自檢校法提出一種脈沖TLS掃描儀檢校方法，以降低參數(shù)估計(jì)的相關(guān)性，保證參數(shù)估計(jì)的可靠性。

2 建模與優(yōu)化設(shè)計(jì)

2.1 自檢校建模

TLS檢校的首要步驟是解決掃描儀幾何建模的問題。由于掃描儀的專利設(shè)計(jì)，用戶通常難以了解其內(nèi)部構(gòu)造，需要對(duì)儀器模型做一定的假設(shè)［12］。TLS在測(cè)量原理（見圖1）上類似于免棱鏡測(cè)量的全站儀［13］，本文擬采用全站儀的誤差模型作為自檢校的基礎(chǔ)模型。

圖1 地面三維脈沖激光掃描測(cè)量的基本原理Fig.1 Basic principle for terrestrial laser scanning

參考全站儀的誤差模型，擬定該模型的主要檢校參數(shù)（Calibration Parameters，CPs）為激光測(cè)距儀的零誤差（加常數(shù)）、水平軸誤差、視準(zhǔn)軸誤差、豎盤指標(biāo)差以及尺度誤差。采用這些參數(shù)實(shí)施檢校后，殘差中可能還會(huì)殘留一些系統(tǒng)趨勢(shì)。對(duì)這些趨勢(shì)進(jìn)行分析以確定其他的CPs。

2.2 檢校網(wǎng)優(yōu)化設(shè)計(jì)

TLS檢校網(wǎng)的最優(yōu)設(shè)計(jì)應(yīng)能保持與其他系統(tǒng)參數(shù)較低相關(guān)性的前提下，提供CPs的可靠估計(jì)［14］，主要分為零類設(shè)計(jì)和一類設(shè)計(jì)兩部分。

2.2.1 零類設(shè)計(jì)

檢校網(wǎng)零類設(shè)計(jì)又稱檢校網(wǎng)基準(zhǔn)設(shè)計(jì)，可基于最小化約束或內(nèi)部約束定義TLS檢校網(wǎng)基準(zhǔn)。常見的內(nèi)部約束是通過對(duì)象點(diǎn)坐標(biāo)（Object Point coordinates，OPs）實(shí)施基準(zhǔn)約束，其約束參數(shù)精度較好，但會(huì)導(dǎo)致一些系統(tǒng)參數(shù)間（主要是外 方 位 參 數(shù)（Exterior Orientation Parameters，EOPs）與CPs間的高相關(guān)性。

最小化約束通常是固定TLS測(cè)站的方位和位置，或固定非共線的3個(gè)對(duì)象點(diǎn)的坐標(biāo)，但該約束易增加一些CPs和OPs間的相關(guān)性［15］。所以，以上約束都不是TLS檢校網(wǎng)基準(zhǔn)設(shè)計(jì)的最佳選擇。本文摒棄此類強(qiáng)制約束，采用最小二乘配置先驗(yàn)隨機(jī)信息柔性定義檢校網(wǎng)基準(zhǔn)。

2.2.2 一類設(shè)計(jì)

檢校網(wǎng)一類設(shè)計(jì)，即檢校網(wǎng)結(jié)構(gòu)圖形設(shè)計(jì)?；谇叭藢?shí)施相機(jī)檢校以及TLS檢校的成果與經(jīng)驗(yàn)，為得到CPs的低相關(guān)高精度估計(jì)，本文擬采用以下檢校網(wǎng)配置：

（1）使用兩個(gè)或兩個(gè)以上掃描儀位置估計(jì)測(cè)距儀加常數(shù)；（2）使用大量的點(diǎn)狀標(biāo)志提供平差觀測(cè)的高冗余度，確保有效檢測(cè)非模型系統(tǒng)誤差導(dǎo)致的殘差趨勢(shì)；（3）檢校標(biāo)志在垂直方向均勻分布，保證垂直角依賴誤差（如水平軸誤差、視準(zhǔn)軸誤差）的可靠估計(jì)；（4）從同一標(biāo)稱位置實(shí)施正交掃描，以降低EOPs和CPs間的相關(guān)性。

3 方 法

傳統(tǒng)利用參數(shù)化模型或混合最小二乘平差模型實(shí)施TLS自檢校，通常把一部分變量視為觀測(cè)值，其他視為參數(shù)。伴隨觀測(cè)值的是以標(biāo)準(zhǔn)差或方差-協(xié)方差（Variance-Covariance，VCV）矩陣形式存在的先驗(yàn)隨機(jī)信息，而對(duì)于參數(shù)并沒有給出這些信息。

本文提出的方法通過最小二乘配置對(duì)檢校平差中的所有參數(shù)配置先驗(yàn)權(quán)重，把模型中的全部變量作為觀測(cè)值。對(duì)TLS自檢校來說，其主要優(yōu)勢(shì)是易于引入EOPs先驗(yàn)隨機(jī)信息。同時(shí)，通過“隨機(jī)約束”定義檢校網(wǎng)基準(zhǔn)，從而避免了強(qiáng)制約束易導(dǎo)致參數(shù)間高相關(guān)性的不足。

3.1 檢校模型

本文方法的底層數(shù)學(xué)模型是附有檢校參數(shù)與外方位信息的掃描儀坐標(biāo)系統(tǒng)與外部坐標(biāo)系統(tǒng)之間的剛性變換。假設(shè)在檢校場(chǎng)中利用點(diǎn)狀檢校標(biāo)志進(jìn)行檢校，且每站實(shí)施一次掃描，則從第i站掃描測(cè)量第j個(gè)檢校標(biāo)志的條件方程為：

式中：p為掃描儀測(cè)站數(shù)，m為檢校標(biāo)志數(shù)，ΔX(i)表示第i個(gè)掃描站的平移向量，R(ω(i)，?(i)，κ(i))表示坐標(biāo)系統(tǒng)間的旋轉(zhuǎn)矩陣，X je為第j個(gè)檢校標(biāo)志的外部坐標(biāo)向量，x(i)scanj為在第i站上掃描第j個(gè)檢校標(biāo)志校正掃描儀CPs影響后的掃描儀坐標(biāo)向量，ΔX，ΔY，ΔZ，ω，φ，κ為掃描儀外方位參數(shù)。其中：

式中r，φ，θ分別為斜距、水平角與垂直角。誤差項(xiàng) Δr，Δφ，Δθ的 建 模 為 ：Δr=a0，Δφ=b1/cosθj(i)+b2tanθj(i)，Δθ=c0，而a0，b1，b2和c0分別表示測(cè)距儀加常數(shù)、視準(zhǔn)軸誤差、水平軸誤差與豎盤指標(biāo)差。由于檢校場(chǎng)尺寸通常較為有限，尺度參數(shù)估計(jì)較弱，本文檢校模型中沒有包含該誤差參數(shù)。

由于參數(shù)和觀測(cè)隱性相關(guān)，選擇線性化混合平差模型。式（1）線性化可得：

式中：A，B分別為式（1）關(guān)于參數(shù)和掃描觀測(cè)的衍生矩陣，W與W X為閉合差向量，X為參數(shù)向量，其近似值為X0；Xobs，V X分別為參數(shù)“觀測(cè)”值與“觀測(cè)”參數(shù)有關(guān)的殘差向量，uCP為CPs數(shù)量。式（3）和式（4）聯(lián)立得到：

式中：向量X，W X，V X及矩陣A被分解為三部分，分別對(duì)應(yīng)掃描站EOPs(AEOP，XEOP，WEOP，VEOP)、掃 描 儀 CPs(ACP，XCP，WCP，VCP)與 OPs（AOP，XOP，WOP，VOP）。

其平差先驗(yàn)VCV矩陣計(jì)算如下：

式中：C為VCV矩陣，Q為輔因子矩陣，P為權(quán)重矩陣。在平差開始，假設(shè)先驗(yàn)方差因子σ20=1.0，并假設(shè)式（6）中系統(tǒng)參數(shù)各不相關(guān)。VCV的子矩陣Cobs為：

式中：σr，σφ，σθ分別表示測(cè)距、水平方向與垂直角的先驗(yàn)信息（標(biāo)準(zhǔn)差），參考儀器說明書設(shè)定；σbeam表示激光光斑內(nèi)測(cè)距位置的不確定性，約等于波束寬度的1/4。

X的法方程系統(tǒng)表示為：

式中：

式中：

則解變?yōu)椋?/p>

為消除設(shè)計(jì)矩陣基準(zhǔn)秩虧及應(yīng)對(duì)傳統(tǒng)約束導(dǎo)致的高相關(guān)性，代替引入OPs作為觀測(cè)值，對(duì)其坐標(biāo)引入內(nèi)部約束。該模型中隨機(jī)內(nèi)部約束為：

式中：VC是約束下的殘差向量，D為內(nèi)部約束矩陣。本文估計(jì)的不是OPs，而是其近似值的改正（XOP）。將式（5）轉(zhuǎn)變?yōu)橄铝芯仃嚪匠蹋?/p>

其約束平差的先驗(yàn)VCV矩陣為：

式中：CC，QC和PC分別為約束下“觀測(cè)”的VCV矩陣、輔因子矩陣與權(quán)重矩陣。具體實(shí)施中，矩陣CC和QC可設(shè)置為單位矩陣乘以一個(gè)極小數(shù)，矩陣PC中對(duì)角元素趨向無(wú)窮大。同式（12）求解內(nèi)部約束的平差解。

3.2 方差分量估計(jì)與粗差檢測(cè)

在TLS自檢校過程中需要利用后驗(yàn)方差分量估計(jì)觀測(cè)精度，以評(píng)估檢校平差中使用的先驗(yàn)觀測(cè)權(quán)重是否合理。此外，為了成功實(shí)施后驗(yàn)方差分量估計(jì)須剔除粗差，而要有效檢驗(yàn)粗差，觀測(cè)VCV矩陣應(yīng)正確無(wú)誤，因此，需同時(shí)實(shí)施方差分量估計(jì)和粗差檢驗(yàn)。

本文采用Forstner法同時(shí)實(shí)施方差分量估計(jì)和魯棒估計(jì)［16］。首先利用檢校模型中計(jì)算的方差-協(xié)方差矩陣實(shí)施最小二乘平差，然后通過以下權(quán)重基于平差殘差計(jì)算新的等價(jià)權(quán)矩陣，即有：

式中σl i和v i分別為第i站觀測(cè)的先驗(yàn)標(biāo)準(zhǔn)差和殘差。根據(jù)以上加權(quán)法則，具有較大殘差的觀測(cè)被視為異常值，使用殘差指數(shù)函數(shù)降低其權(quán)重，持續(xù)迭代直至收斂。

最后，估計(jì)觀測(cè)方差分量，并利用F測(cè)試檢驗(yàn)它是否與先驗(yàn)值存在顯著差異，若不顯著，則接受平差結(jié)果；若顯著，則計(jì)算新的后驗(yàn)權(quán)矩陣，并重復(fù)上述檢校流程。

4 實(shí)驗(yàn)

為了測(cè)試本文方法的性能以及找到CPs之間以及CPs和其他系統(tǒng)參數(shù)間低相關(guān)性的最優(yōu)檢校網(wǎng)設(shè)計(jì)，利用掃描儀Leica Scan Station（具有直接地理定位功能）在一個(gè)12 m×9 m×3 m的室內(nèi)檢校場(chǎng)中實(shí)施自檢校實(shí)驗(yàn)（見圖2）。

圖2 檢校場(chǎng)設(shè)計(jì)平面圖Fig.2 Schematic diagram of self-calibration field

該實(shí)驗(yàn)場(chǎng)包含144個(gè)均勻分布于檢校場(chǎng)空間（墻壁、地板和天花板上）的點(diǎn)狀檢校標(biāo)志，其垂直分布范圍為-45°～83°，幾乎覆蓋了整個(gè)垂直維掃描視場(chǎng)，分別采用2站和4站掃描設(shè)置實(shí)施兩個(gè)系列的模擬自檢校實(shí)驗(yàn)（具體設(shè)計(jì)見表1）。

表1 自檢校模擬設(shè)計(jì)Tab.1 Specific indexes of self-calibration test

實(shí)驗(yàn)中采用的掃描儀觀測(cè)先驗(yàn)標(biāo)準(zhǔn)差從儀器說明書中獲取，誤差項(xiàng)σbeams忽略不計(jì)，掃描儀CPs先驗(yàn)標(biāo)準(zhǔn)差參考前人研究成果設(shè)置（見表2）；OPs先驗(yàn)標(biāo)準(zhǔn)差基于假設(shè)“其中一站掃描的檢校標(biāo)志坐標(biāo)可作為OPs的先驗(yàn)值”；EOPs先驗(yàn)標(biāo)準(zhǔn)差參考掃描儀可實(shí)現(xiàn)的對(duì)中整平精度；掃描方位先驗(yàn)精度可設(shè)為一較大值，以實(shí)現(xiàn)該參數(shù)的完全自由。最后，加入隨機(jī)噪聲（對(duì)應(yīng)于各個(gè)觀測(cè)的標(biāo)準(zhǔn)差）到掃描儀觀測(cè)、掃描儀中心坐標(biāo)和“傾”角中。

表2 自檢校中模擬CPsTab.2 Simulation of CPs in self-calibration test

根據(jù)表1中的實(shí)驗(yàn)設(shè)計(jì)，分別實(shí)施兩個(gè)系列的檢校實(shí)驗(yàn)。在第一個(gè)系列實(shí)驗(yàn)中不加入任何已知點(diǎn)，在第二個(gè)實(shí)驗(yàn)系列中加入3個(gè)位置分布良好的已知檢校標(biāo)志（XYZ坐標(biāo)先驗(yàn)標(biāo)準(zhǔn)差約為1 mm），以研究已知點(diǎn)對(duì)模擬結(jié)果的影響。

5 結(jié)果與分析

5.1 實(shí)驗(yàn)結(jié)果

檢校實(shí)驗(yàn)結(jié)果見表3與表4（篇幅所限，僅給出部分掃描站檢校數(shù)據(jù)）以及圖3和圖4（圖中無(wú)已知點(diǎn)參與的結(jié)果以黑色表示，包含3個(gè)已知點(diǎn)的模擬結(jié)果以灰色表示）。表3和表4中給出了估計(jì)的CPs以及其標(biāo)準(zhǔn)差信息、EOPs和OPs的相關(guān)系數(shù)（最大絕對(duì)值）以及平均標(biāo)準(zhǔn)偏差信息。由于掃描儀自檢校中主要關(guān)注CPs的精確估計(jì)，所以表中僅給出了CPs和其他系統(tǒng)參數(shù)間的相關(guān)系數(shù)。這里用ρ表示相關(guān)系數(shù)。

圖3 自檢校模擬得到的CPs估計(jì)和標(biāo)準(zhǔn)差Fig.3 Estimation and standard deviations of CPs in self-calibration simulation

圖4 自檢校模擬中不同參數(shù)間最大（絕對(duì)值）相關(guān)系數(shù)Fig.4 Maximal correlation coefficients（absolute value）between different parameters

表3 無(wú)已知點(diǎn)參與的自檢校模擬結(jié)果Tab.3 Simulation results of self-calibration without points （°）

表4 包含3個(gè)已知點(diǎn)的自檢校模擬結(jié)果Tab.4 Simulation results of self-calibration with 3 points

5.2 結(jié)果分析

由表3～表4和圖3分析可知，掃描儀加常數(shù)、水平軸誤差和豎盤指標(biāo)差等CPs在所有模擬設(shè)計(jì)下都得到了高精度估計(jì)，其估計(jì)標(biāo)準(zhǔn)差約為其估計(jì)值的1/10，且估計(jì)值與模擬值吻合較好；在2站掃描設(shè)計(jì)下，視準(zhǔn)軸誤差b1與掃描儀方位κ高度相關(guān)（ρ約為-0.8），導(dǎo)致這兩個(gè)參數(shù)估計(jì)結(jié)果較差。雖然通過多次掃描，b1的標(biāo)準(zhǔn)差有所下降，但并沒有降低b1和κ間的相關(guān)性，同時(shí)采取不同測(cè)站高度掃描未發(fā)現(xiàn)對(duì)參數(shù)估計(jì)有顯著影響。除b1與κ間的高相關(guān)性外，所有CP-EOP間的相關(guān)性都相對(duì)較低，且在檢校中影響不顯著。采用4站掃描設(shè)計(jì)顯著改善了b1的估計(jì)結(jié)果（其估計(jì)標(biāo)準(zhǔn)差比2站配置下改善約2倍），且ρ下降到0.5。

分析圖4可知，掃描儀中心位置的已知有助于降低a0和ΔX及ΔY間的相關(guān)性（ρ最高約為0.3）。a0和ΔY間的相關(guān)性高于它和ΔX間的相關(guān)性，這可能與外部坐標(biāo)系統(tǒng)的Y軸沿檢校場(chǎng)長(zhǎng)尺寸方向有關(guān)；此外，CP-OP相關(guān)性非常低，觀測(cè)到的最大相關(guān)性（0.4～0.5）存在于加常數(shù)a0和X，Y標(biāo)志點(diǎn)坐標(biāo)間。無(wú)已知點(diǎn)檢校的最大相關(guān)系數(shù)為0.28，使用3個(gè)已知點(diǎn)檢校的最大相關(guān)系數(shù)為0.39，以上相關(guān)性對(duì)TLS自檢校而言并不算高；而OPs估計(jì)的平均標(biāo)準(zhǔn)差處于亞毫米水平，由于加常數(shù)a0和標(biāo)志點(diǎn)X，Y軸向間的相關(guān)性相對(duì)較高，OPs平面標(biāo)準(zhǔn)差稍差于高程標(biāo)準(zhǔn)差。

最后結(jié)合圖3可知，檢校網(wǎng)中已知點(diǎn)的存在對(duì)CPs估計(jì)無(wú)顯著影響，但對(duì)參數(shù)相關(guān)性有較大影響，如b1和κ間的相關(guān)性以及a0與ΔX，ΔY間的相關(guān)性。相對(duì)于無(wú)已知點(diǎn)檢校實(shí)驗(yàn)，其相關(guān)性變大，CP-OP間的相關(guān)性也略有增加，而EOPs和OPs的標(biāo)準(zhǔn)差卻略有下降。

6 結(jié)論

本文以全站儀誤差模型為基礎(chǔ)提出了一種脈沖TLS儀器系統(tǒng)誤差自檢校方法，通過對(duì)估計(jì)參數(shù)配置先驗(yàn)權(quán)重信息以及對(duì)檢校網(wǎng)進(jìn)行優(yōu)化設(shè)計(jì)，降低了脈沖TLS自檢校平差中大多數(shù)CPs和其他系統(tǒng)參數(shù)間（尤其是與外方位參數(shù)EOPs間）的相關(guān)性，從而提高了估計(jì)參數(shù)的可靠性。檢校實(shí)驗(yàn)結(jié)果表明，4站掃描檢校配置的設(shè)計(jì)結(jié)果顯著優(yōu)于2站掃描設(shè)計(jì)，而掃描儀設(shè)站的高度差異、同一測(cè)站多次掃描、正交掃描等措施并未對(duì)參數(shù)估計(jì)產(chǎn)生顯著影響。此外，為可靠估計(jì)掃描儀加常數(shù)，需以約1 mm的高精度測(cè)定掃描儀站坐標(biāo)，由此顯示出采用隨機(jī)內(nèi)部約束方法定義檢校網(wǎng)基準(zhǔn)的優(yōu)勢(shì)。