岳春宇 何紅艷 鮑云飛 邢坤 周楠

（北京空間機(jī)電研究所，北京 100094）

0 引言

美國 1985年提出采用星載激光測距儀測量兩極的冰面變化，并在 2003年將地球科學(xué)激光高度計（GLAS）作為 ICESAT衛(wèi)星的主要載荷之一發(fā)射升空，實現(xiàn)了星載激光高度計對地觀測。星載激光高度計采取主動的方式精確獲得地面高度信息，可以有效彌補(bǔ)衛(wèi)星光學(xué)影像三維探測能力的不足[1-4]。星載激光高度計為外推定位和數(shù)字地面模型（DSM）獲取提供了新的方式，其測距精度很高，一般為分米級，可以滿足大部分對地觀測任務(wù)。但是其平面精度較低，只有百米級甚至千米級[5]，而平面精度是制約星載激光高度計應(yīng)用的主要障礙，所以研究如何提高星載激光高度計的平面精度具有重要意義。

由于星載激光高度計發(fā)射頻率低于航空或地面激光雷達(dá)，不能獲取密集的點云數(shù)據(jù)，尤其當(dāng)衛(wèi)星高軌運行時，相鄰激光腳點間距為百米左右[6-7]。這也造成了星載激光高度計在軌檢校的困難，嚴(yán)重影響了其定位精度。本文從星載激光高度計幾何定位模型出發(fā)，研究星載激光高度計直接對地定位時各外方位元素的誤差傳播規(guī)律，分析定位誤差的主要來源及其影響，為星載激光高度計設(shè)計和應(yīng)用提供參考。

1 星載激光高度計幾何定位模型

星載激光高度計通過發(fā)射和接收激光脈沖的時間差來推算目標(biāo)與衛(wèi)星之間的距離，再根據(jù)衛(wèi)星的姿態(tài)軌道參數(shù)解算出衛(wèi)星到該目標(biāo)處參考橢球的距離，從而得到目標(biāo)的高程[8-10]。星載激光高度計對地觀測原理見圖1[11-13]。

圖1 激光測高儀幾何定位原理Fig.1 Geometric positioning of space-borne laser altimeter

圖1中S為激光高度計中激光發(fā)射點的空間位置，在空間坐標(biāo)系中坐標(biāo)為 （Xs,Ys,Zs）；P為觀測的地面點，在空間坐標(biāo)系下坐標(biāo)為（X,Y,Z）；L為距離測量值；θ為激光高度計指向角；（φ,ω,κ）為激光高度計外方位角元素。激光高度計對地觀測方程為

式中 （a1,a2,a3,b1,b2,b3,c1,c2,c3）為激光高度計外方位角元素（φ,ω,κ）的方向余弦參數(shù)，求解公式為

2 星載激光高度計定位誤差傳播模型

2.1 定位誤差傳播系數(shù)

式（1）中，星載激光高度計對地定位的誤差源為外方位元素誤差（ΔXS，ΔYS，ΔZS，Δφ，Δω，Δκ）、距離測量誤差ΔL和指向角誤差Δθ。在平差計算中，最終需要解算這8個誤差改正數(shù)，修正觀測的初值，然后再根據(jù)式（1）來對地觀測。以這8個誤差改正數(shù)（ΔXS，ΔYS，ΔZS，Δφ，Δω，Δκ，ΔL，Δθ）為未知數(shù)，用泰勒公式對式（1）線性展開，得到式（3）。

式中 （X0,Y0,Z0）為地面點坐標(biāo)觀測值，作為展開式的近似值。其中各誤差改正數(shù)的誤差傳播系數(shù)為：

2.2 定位誤差傳播分析

在僅考慮星下點測距，不側(cè)擺的情況下，設(shè)定距離測量值與軌道高度相等為500km，即L0=500km，令φ0=ω0=κ0=0，激光高度計指向星下點，即θ0=0。此時各誤差變量的誤差傳播系數(shù)見表1。

表1 激光高度計定位模型誤差傳播系數(shù)Tab.1 Geopositioning error propagation coefficients of laser altimeter system

由表1可以得出，激光高度計對地定位X方向的誤差源與外方位元素中XS坐標(biāo)和俯仰角φ的測量精度有關(guān)，并且主要受φ測量精度的影響，影響程度與軌道高度有很大關(guān)系；激光高度計對地定位Y方向的誤差源與外方位元素中YS坐標(biāo)、側(cè)擺角ω和激光高度計指向角θ的測量精度有關(guān)，并且主要受ω和θ的測量精度的影響，影響程度也與軌道高度有很大關(guān)系；激光高度計對地定位Z方向的誤差源與外方位元素中ZS坐標(biāo)與測距值L的測量精度有關(guān)，并且主要受L測量精度的影響。由于激光測距精度一般較高，且其誤差傳播系數(shù)較小，所以星載激光高度計的高程精度較高。衛(wèi)星的定姿定軌精度雖然很高，但是姿態(tài)角及指向角的誤差傳播系數(shù)很大，是激光測距誤差傳播系數(shù)的幾十萬倍，所以整體上星載激光高度計的平面定位精度要低于高程定位精度。

根據(jù)表1中的誤差傳播系數(shù)可知，激光高度計的定位精度主要受外方位元素中角元素和指向角的影響，且主要影響平面定位精度。下面就主要針對外方位角元素中φ、ω與θ對激光高度計平面定位精度的影響進(jìn)行分析。

激光高度計的外方位線元素（XS，YS，ZS）由衛(wèi)星平臺軌道參數(shù)計算得到，測定精度可達(dá)到分米級[14]，這里取誤差ΔXS=ΔYS=ΔZS=0.1m。外方位角元素中κ誤差取1″，即Δκ=1″。分別討論Δφ、Δω和Δθ變化對定位精度影響時，一方變化則另兩方暫定誤差為1″。參考GLAS系統(tǒng)激光測距的誤差補(bǔ)償方法[15]，補(bǔ)償大氣延遲、固體潮以及光斑大小等對測距的影響，在平坦地區(qū)，假設(shè)激光高度計指向角誤差Δθ=1″，激光測距誤差ΔL=0.3m。根據(jù)誤差傳播定律，激光高度計的定位誤差為：

根據(jù)式（7），φ、ω與θ精度分別對激光高度計平面定位精度的影響見表2。

在表 2中，激光高度計對外方位角元素的觀測誤差變化比較敏感，當(dāng)直接定位時，Δφ、Δω或 Δθ誤差在1″以上時，只能滿足1∶5萬以上的小比例尺測圖要求。隨著俯仰角誤差Δφ增加，激光高度計X方向定位誤差增大，而Y方向定位誤差不變，平面定位誤差增加；隨著側(cè)擺角誤差Δω或Δθ增加，激光高度計Y方向定位誤差增大，而X方向定位誤差不變，平面定位誤差增加。由于Δω與Δθ的誤差傳播系數(shù)絕對值相同，所以兩者分別對平面精度的影響相同。由于Y方向定位誤差受Δω與Δθ共同影響，所以激光高度計Y方向定位誤差稍微大于X方向定位誤差。

由于φ、ω是姿態(tài)參數(shù)與衛(wèi)星平臺有關(guān)，假設(shè)其測量精度一致，即Δφ=Δω。φ、ω與θ精度共同對激光高度計平面定位精度的影響見表3。

表2 φ、ω與θ誤差分別對平面精度的影響Tab.2 Influence of φ, ω and θ on ground positioning plane accuracy

表3 φ、ω與θ精度共同對平面精度的影響Tab.3 Comprehensive influences of φ, ω and θ on ground positioning plane accuracy

在表3中，當(dāng)Δφ、Δω和Δθ增大時，平面定位誤差也增大；當(dāng)Δφ、Δω和Δθ中任一個超過3″時，定位誤差已經(jīng)接近10m。在使用激光高度計作為控制數(shù)據(jù)進(jìn)行1∶5萬測圖時，姿態(tài)控制需要在1″以內(nèi)。

3 結(jié)束語

本文從星載激光高度計幾何定位模型出發(fā)，推導(dǎo)了星載激光高度計幾何定位各誤差源的誤差傳播系數(shù)，分析了星載激光高度計幾何定位誤差傳播規(guī)律。根據(jù)誤差傳播系數(shù)可以得出下述結(jié)論：

1）星載激光高度計定位精度主要受俯仰角φ、側(cè)擺角ω和激光高度計指向角θ的觀測精度影響，影響的程度與軌道高度正相關(guān)；

2）偏航角κ對星載激光高度計定位精度沒有影響；

3）星載激光高度計定位平面精度和高程精度誤差源不同，平面精度中X方向定位精度主要受俯仰角φ的觀測精度影響，軌道高度越高，影響越大；Y方向定位精度主要受側(cè)擺角ω和激光高度計指向角θ的測量精度影響，軌道高度越高，影響越大；高程精度主要與測距精度有關(guān)；

4）在星載激光高度計幾何定位模型誤差傳播系數(shù)中，影響平面精度的誤差源多于影響高程精度的誤差源，且影響平面精度的誤差源對應(yīng)的誤差傳播系數(shù)遠(yuǎn)大于影響高程精度的誤差源的傳播系數(shù)。在當(dāng)前衛(wèi)星定軌精度較高的情況下，星載激光高度計高程精度幾乎僅受測距精度影響。因為影響某項指標(biāo)的各項誤差源的精度乘以其誤差傳播系數(shù)，得到的積取平方和，將得到的和取平方根，等于該項指標(biāo)最終的誤差值，所以星載激光高度計的平面精度低于高程精度。

通過仿真分析，可知星載激光高度計幾何定位對外方位元素和指向角精度要求較高，直接定位精度較低，不能滿足大比例尺測圖要求；X方向定位精度高于Y方向定位精度。

本文僅從星載激光高度計幾何定位模型出發(fā)，沒有考慮星載激光高度計足印光斑及回波信號等因素對定位精度的影響，后續(xù)研究應(yīng)在幾何定位模型的基礎(chǔ)上加入激光器本身光學(xué)性質(zhì)對定位影響的修正。

（

）

[1] 李磊, 鄭永超, 彭鳳超, 等.地形測繪激光成像雷達(dá)技術(shù)研究[J].紅外與激光工程, 2006, 35（增刊）: 294- 298.

LI Lei，ZHENG Yongchao，PENG Fengchao, et al.Research of Three-dimension Imaging Lidar on Land Topography[J].Infrared and Laser Engineering, 2006, 35（Supplement）: 294-298.（in Chinese）

[2] 岳春宇, 鄭永超, 陶宇亮.星載激光測高儀輔助衛(wèi)星攝影測量淺析[J].航天返回與遙感, 2013, 34（4）: 71- 76.

YUE Chunyu, ZHENG Yongchao, TAO Yuliang.Study on Space-borne Laser Altimeter Supported Satellite Photogrammetry[J].Spacecraft Recovery & Remete Sensing, 2013, 34（4）: 71-76.（in Chinese）

[3] Wehr A, Lohr U.Airborne Laser Scanning-an Introduction and Overview[J].ISPRS Journal of Photogrammetry and Remote Sensing, 1999, 54（2）: 68-82.

[4] Baltsavias E P.Airborne Laser Scanning: Existing Systems and Firms and Other Resourses[J].ISPRS Journal of Photogrammetry and Remote Sensing, 1999, 54（2）: 164-198.

[5] 李然, 王成, 蘇國中, 等.星載激光雷達(dá)的發(fā)展和應(yīng)用[J].科技導(dǎo)報, 2007, 25（14）: 58-63.

LI Ran, WANG Cheng, SU Guozhong, et al.Development and Applications of Spaceborne LiDAR[J].Science & Technology Review, 2007, 25（14）: 58-63.（in Chinese）

[6] Sungkoo Bae, Charles Web, Bob Schutz.GLAS PAD Calibration Using Laser Reference Sensor Data[C].In AIAA/AAS Astrodynamics Specialist Conference and Exhibit, Providence, RI, USA, 2004:1-10.

[7] 鄂棟臣, 徐瑩, 張小紅.ICESat衛(wèi)星及其在南極 Dome A 地區(qū)的應(yīng)用[J].武漢大學(xué)學(xué)報: 信息科學(xué)版, 2007, 32（12）:1139-1142.

E Dongchen, XU Ying, ZHANG Xiaohong.ICESat’s Performance and Its Application in Dome A Area in Antarctica[J].Geomatics and Information Science of Wuhan University, 2007, 32（12）: 1139-1142.（in Chinese）

[8] Garvin J B, Bufton J, Blair J B, et al.Observations of the Earth’s Topograph from the Shuttle Laser Altimeter （SLA）: Laser-pulse Echo-recovery Measurements of Terrestrial Surfaces[J].Physics and Chemistry of the Earth, 1998, 23（9）:1053-1068.

[9] Abshire J B, Sun Xiaoli, Riris H, et al.Geoscience Laser Altimeter System （GLAS） on the ICESat Mission: On-orbit Measurement Performance[J].Geophysical Research Letters, 2005, 32（21）: 1-4.

[10] Axelsson P.Processing of Laser Scanner Data-algorithms and Applications[J].ISPRS Journal of Photogrammetry and Remote Sensing, 1999, 54（2）: 138-147.

[11] 文漢江, 程鵬飛.ICESAT/GLAS激光測高原理及其應(yīng)用[J].測繪科學(xué), 2005, 30（5）: 33-35.

WEN Hanjiang, CHENG Pengfei.Introducton to Principle of ICESAT/GLAS Laser Altimetry and its Applications[J].Science of Surveying and Mapping, 2005, 30（5）: 33-35.（in Chinese）

[12] Zwally H J.ICESat’s Laser Measurements of Polar Ice, Atmosphere, Ocean and Land[J].Journal of Geomatics, 2002, 34:405-445.

[13] 余俊鵬, 孫世君.衛(wèi)星攝影測量觀測技術(shù)發(fā)展的若干思路[C].第二十三屆全國空間探測學(xué)術(shù)交流會, 廈門, 2010: 1-6.

YU Junpeng, SUN Shijun.Some Strategies about Satellite Photogrammetry[C].The 23rd National Conference on Space Explorer, Xiamen, 2010: 1-6.（in Chinese）

[14] 余俊鵬, 孫世君, 毛建杰.衛(wèi)星遙感影像外方位元素的誤差傳播研究[J].航天返回與遙感, 2011, 32（1）: 18-23.

YU Junpeng, SUN Shijun, MAO Jianjie.Study on Error Propagation of Exterior Orientation Elements of Satellite Remote Sensing Imagery[J].Spacecraft Recovery & Remete Sensing, 2011, 32（1）: 18-23.（in Chinese）

[15] 范春波, 李建成, 王丹, 等.ICESAT/GLAS激光腳點定位及誤差分析[J].大地測量與地球動力學(xué), 2007, 27（1）: 104-106.

FAN Chunbo, LI Jiancheng, WANG Dan, et al.ICEsat/GLAS Laser Footprint Geolocation and Error Analysis[J].Journal of Geodesy and Geodynamic, 2007, 27（1）: 104-106.（in Chinese）