岳春宇 周楠 楊舒琪 何紅艷

陸地生態(tài)系統(tǒng)碳監(jiān)測衛(wèi)星植被測量子系統(tǒng)無控定位方法

岳春宇 周楠 楊舒琪 何紅艷

（北京空間機電研究所，北京 100094）

陸地生態(tài)系統(tǒng)碳監(jiān)測衛(wèi)星植被測量子系統(tǒng)包含多波束激光雷達和0°相機。同時開機對地觀測時，多波束激光雷達可提供高精度測距值轉化的廣義高程控制點，參與0°相機平差定位，實現(xiàn)在軌無控測圖。文章基于陸地生態(tài)系統(tǒng)碳監(jiān)測衛(wèi)星植被測量子系統(tǒng)探測體制，建立圖像與激光對地觀測模型及其幾何關系模型，引入公開粗格網(wǎng)數(shù)字高程模型（DEM）約束，實現(xiàn)激光輔助圖像攝影測量平差計算。通過仿真實驗分析，陸地生態(tài)系統(tǒng)碳監(jiān)測衛(wèi)星植被測量子系統(tǒng)無控測圖方法可提升0°相機圖像平面定位精度，仿真數(shù)據(jù)最優(yōu)實驗結果從613.05m提高到15.35m。對參與平差計算的激光控制點數(shù)量和布設位置開展了討論，仿真數(shù)據(jù)一景圖像當采用上中下三行均勻布設激光控制點時，精度最高。

植被測量子系統(tǒng) 無控定位 聯(lián)合平差 陸地生態(tài)系統(tǒng)碳監(jiān)測衛(wèi)星

0 引言

陸地生態(tài)系統(tǒng)碳監(jiān)測衛(wèi)星主要用于陸地生態(tài)系統(tǒng)碳儲量監(jiān)測、生態(tài)資源調查、國家重大生態(tài)工程監(jiān)測評價等領域，并為環(huán)保、測繪、氣象、農(nóng)業(yè)、減災等領域提供業(yè)務支撐和研究服務，提高中國應對全球氣候變化的話語權和主導權。陸地生態(tài)系統(tǒng)碳監(jiān)測衛(wèi)星植被測量子系統(tǒng)由多波束激光雷達、光軸監(jiān)視相機和0°相機組成。當衛(wèi)星處于光照區(qū)且側擺角小于5°的情況下，可同時進行多波束激光雷達和0°相機高分辨率對地觀測。激光雷達采用全波形探測體制接收地表激光回波，0°相機采用星下點對地成像模式獲取高分辨率多光譜圖像數(shù)據(jù)，而光軸監(jiān)視相機同時獲取激光雷達各波束及0°相機的指向信息。通過光軸監(jiān)視中繼匹配，可以得到多波束激光雷達對地觀測腳點在0°相機中的位置，達到同時獲取觀測區(qū)域平面高程三維信息的目的。

傳統(tǒng)的衛(wèi)星對地測繪采用立體相機空中三角測量，實現(xiàn)地表三維模型重構，即數(shù)字表面模型（Digital Surface Model，DSM）+數(shù)字正攝影像（Digital Orthophoto Map，DOM）[1-3]?！百Y源三號”02星搭載了我國首臺星載對地觀測激光雷達，文獻[4]使用激光輔助立體測量，在陜西渭南試驗區(qū)進行試驗，結果表明，可以將“資源三號”02星立體測圖無控高程精度從11.54m提高到1.90m，并建議后續(xù)立體測繪衛(wèi)星可搭載激光雷達業(yè)務化應用?！案叻制咛枴毙l(wèi)星要實現(xiàn)1:10 000衛(wèi)星立體測圖時，依靠攝影測量外方位元素觀測精度，無法滿足高程1.5m和平面5m的精度設計指標。增加激光雷達，觀測數(shù)據(jù)作為廣義高程控制點，利用光學立體測繪定位高程與激光測高的差異建立高程約束，進行聯(lián)合空中三角平差測量，最終實現(xiàn)無控平面定位精度7.2m，高程精度1.2m，達到了1:10 000的測繪指標[5-6]。當前激光輔助立體測繪機制下，采用激光測高數(shù)據(jù)作為高程控制點的方法對光斑大小要求較高，當光斑較大時，只能選擇平坦地區(qū)的光斑，且控制點分布的均勻性比較難實現(xiàn)。對于非測繪衛(wèi)星弱交會條件下的立體測繪，一般利用地面已有數(shù)字高程模型（Digital Elevation Model，DEM）內插連接點高程方法，避免交會角度太小、平差計算不收斂的問題[7-8]。

陸地生態(tài)系統(tǒng)碳監(jiān)測衛(wèi)星植被測量子系統(tǒng)多波束激光雷達與0°相機同時對地探測成像，多波束激光雷達以40Hz工作頻率提供地面高密度測距格網(wǎng)，可以作為控制點與0°相機圖像聯(lián)合平差定位，實現(xiàn)無控測圖。本文基于陸地生態(tài)系統(tǒng)碳監(jiān)測衛(wèi)星植被測量子系統(tǒng)探測方法，通過建立圖像幾何成像模型和激光對地觀測模型之間的關系，構建激光輔助面陣攝影測量平差模型。最后經(jīng)過數(shù)據(jù)仿真計算，對無地面控制情況下的測圖精度開展實驗驗證。

1 植被測量子系統(tǒng)定位平差模型

1.1 遙感圖像的幾何定位平差模型

衛(wèi)星遙感圖像的成像幾何模型用式（1）表示[9-11]

通過少數(shù)地面控制點，根據(jù)式（1），構建光束法區(qū)域網(wǎng)平差模型，如式（2）所示[12-14]

根據(jù)光束法區(qū)域網(wǎng)平差原理，通過構建上述誤差方程，可以實現(xiàn)遙感圖像檢校和幾何校正以及DOM制作[15-16]。

1.2 激光雷達對地定位模型

激光雷達對地定位幾何模型如式（3）所示[17-21]

陸地生態(tài)系統(tǒng)碳監(jiān)測衛(wèi)星光軸監(jiān)視相機能夠記錄激光出光時的光斑圖像，從而獲得激光出射時的指向角。光軸監(jiān)視相機中激光出光指向與衛(wèi)星坐標系之間的關系如圖1所示。

圖1 光軸監(jiān)視相機坐標系示意圖

式（3）與式（4）聯(lián)立，即構建了激光雷達與光軸監(jiān)視相機的對地定位三維關聯(lián)模型。

1.3 植被測量子系統(tǒng)區(qū)域網(wǎng)平差模型

陸地生態(tài)系統(tǒng)碳監(jiān)測衛(wèi)星植被測量子系統(tǒng)5波束激光雷達與0°相機同時工作。0°相機地面分辨率2m，幅寬20km，5波束激光雷達每波束每秒發(fā)射40個脈沖，獲取與0°相機同視場下均勻分布的測距值，并由光軸監(jiān)視相機記錄光斑的地面位置。配準光軸監(jiān)視相機圖像與0°相機圖像，將激光光斑位置轉刺到0°相機圖像上，作為控制點。聯(lián)立式（1）、（3），建立連接點的激光輔助面陣相機定位方程

在式（2）基礎上，建立激光輔助面陣相機無控制圖，針對上述誤差改正數(shù)的聯(lián)合平差模型

由于激光雷達與0°相機姿軌精度一致且均由衛(wèi)星平臺確定，直接聯(lián)合平差對平面定位精度提升有限，一般需要先對激光雷達幾何檢校后使用[23]。本文為了實現(xiàn)無控測圖，可先將5波束激光雷達生成的激光控制點格網(wǎng)與粗格網(wǎng)DEM數(shù)據(jù)匹配，修正激光控制點平面定位精度。經(jīng)過修正后的激光控制點，可以作為平高控制點，進一步提升系統(tǒng)定位精度。

2 植被測量子系統(tǒng)無控平差實驗與分析

2.1 實驗流程

本文方法算法流程如圖2所示。由多波束激光雷達對地探測幾何模型獲取激光地面三維坐標，與公開粗格網(wǎng)DEM匹配提升平面定位精度，作為激光控制點，與0°相機圖像構建激光輔助自檢校聯(lián)合平差模型，平差計算0°相機及多波束激光雷達幾何模型的內外方位元素改正數(shù)，同時獲取地面點三維坐標。

圖2 激光輔助0°相機無控定位方法

具體流程為：

1）分別對0°相機和多波束激光雷達構建幾何定位模型；

2）根據(jù)多波束激光雷達幾何定位模型計算激光控制點坐標；

3）激光控制點構成DEM與公開粗格網(wǎng)DEM匹配，精化激光控制點平面坐標；

4）分別根據(jù)0°相機和多波束激光雷達幾何定位模型構建其定位誤差方程；

5）以激光控制點為控制點，構建0°相機和多波束激光雷達自檢校聯(lián)合平差模型，解算0°相機和多波束激光雷達定位誤差方程改正數(shù)及地面加密點三維坐標；

6）使用檢查點驗證自檢校平差后的圖像幾何定位精度。

2.2 實驗及分析

采用某區(qū)域“資源三號”03星下視2.1m分辨率圖像和仿真多波束激光雷達數(shù)據(jù)進行無控平差實驗。為了仿真0°相機圖像大小，將“資源三號”03星下視圖像裁剪為20km×20km范圍，激光雷達控制點布設方案采用中間布設、上中下布設和5×5均勻布設三種模式，圖3中三角形代表選用的激光控制點。

圖3 試驗區(qū)及激光控制點布設方案

仿真實驗參數(shù)見表1。

表1 仿真實驗參數(shù)

不實際平差計算，使用該仿真實驗參數(shù)，由共線條件方程誤差傳播公式預估得到的定位精度結果為方向529.61m，方向265.39m，平面592.39m。

采用試驗區(qū)50個檢查點對平差結果進行精度評價。

（1）實驗1

圖像直接無控平差結果見表2。

表2 0°相機圖像直接定位精度

從表2可以看出，仿真圖像直接在沒有任何控制點約束的情況下，平面定位精度為613.05m，與根據(jù)仿真參數(shù)精度預估的定位精度結果相仿，證明平差仿真計算結果正確。

（2）實驗2

采用圖3（b）方案平差結果見表3。

表3 0°相機圖像與中間5波束激光聯(lián)合定位平差精度

加入5波束激光雷達獲取的激光控制點后，平面定位精度有了顯著提升。仿真誤差定軌精度方向優(yōu)于方向，所以在沒有控制約束的情況下，定位精度方向應優(yōu)于方向，與表3中計算結果一致。由于激光控制點與粗格網(wǎng)DEM匹配后，修正了激光控制點的平面定位精度方向到相同水平，所以在表3中平差精度與方向已沒有明顯差異，且得到了整體提升，平面總體定位精度與無控制點相比提升近10倍。

（3）實驗3

采用圖3（c）方案平差結果見表4。

表4 0°相機圖像與上中下15波束激光聯(lián)合定位平差精度

加入少量均勻分布的激光控制點后，平差結果平面精度有了大幅提升，得益于控制點位分布對圖像邊緣和中心都能夠有效控制，且平面定位精度已達到公開粗格網(wǎng)DEM平面精度，基本達到了修正后激光控制點的平面絕對定位精度。

（4）實驗4

采用圖3（c）方案平差結果見表5。

表5 0°相機圖像與5×5均勻分布25波束激光聯(lián)合定位平差精度

在上中下均勻分布的基礎上再加密到5×5共25個激光控制點后，平面定位精度比實驗（3）下降了一倍，但是方向沒有明顯下降。這是由于多波束激光雷達外方位元素一致，獲取的激光控制點手段具有強相關性，超過一定數(shù)量造成過度參數(shù)化[24]。而方向的定位精度主要由激光測距值來控制，激光控制點的激光測距值相對于平面定位精度受到外方位元素的影響小，更加可靠穩(wěn)定，這也是“資源三號”系列衛(wèi)星采用激光測距值作為廣義高程控制點的原因。

由上述實驗分析可知，本文方法使用激光控制點可有效提升0°相機圖像平面定位精度，僅在圖像中央布設激光控制點不如采用少量激光控制點均勻布設在圖像邊緣，每景圖像采用5×3均勻分布的激光控制點，參與激光與圖像聯(lián)合平差，可實現(xiàn)與參考公開粗格網(wǎng)DEM一致的平面定位精度。采用5×5密集分布激光控制點時，則精度有所下降。對于陸地生態(tài)系統(tǒng)碳監(jiān)測衛(wèi)星植被測量子系統(tǒng)高重頻激光雷達來說，可間隔采樣，避免參數(shù)強相關性影響定位精度。

3 結束語

陸地生態(tài)系統(tǒng)碳監(jiān)測衛(wèi)星植被測量子系統(tǒng)無控測圖方法是利用激光測高數(shù)據(jù)和公開粗格網(wǎng)修正后獲取的激光控制點，提供給0°相機圖像和多波束激光雷達定位聯(lián)合平差，實現(xiàn)無控定位測圖。通過仿真實驗分析，以公開粗格網(wǎng)DEM為參考，測圖平面定位精度可達到15.35m。另外根據(jù)多波束激光雷達對地觀測設計，對激光控制點選取不同分布方案，根據(jù)平差結果分析可知，應采用少量均勻分布在圖像邊緣的激光控制點，才能獲取最優(yōu)定位精度。對于激光控制點的分布及數(shù)量對平差結果的影響后續(xù)應進一步討論分析。

[1] 孫承志, 唐新明, 翟亮. 我國測繪衛(wèi)星的發(fā)展思路和應用展望[J]. 測繪科學, 2009, 34(2): 5-7.

SUN Chengzhi, TANG Xinming, ZHAI Liang. Development Ideas and Application Prospects of Surveying and Mapping Satellite in China[J]. Science of Surveying and Mapping, 2009, 34(2): 5-7. (in Chinese)

[2] 李德仁. 我國第一顆民用三線陣立體測圖衛(wèi)星——資源3號測繪衛(wèi)星[J]. 測繪學報, 2012, 41(3): 317-322.

LI Deren. China’s First Civilian Three-line-array Stereo Mapping Satellite: ZY-3[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(3): 317-322. (in Chinese)

[3] 曹海翊, 戴君, 張新偉, 等. “高分七號”高精度光學立體測繪衛(wèi)星實現(xiàn)途徑研究[J]. 航天返回與遙感, 2020, 41(2): 17-28.

CAO Haiyi, DAI Jun, ZHANG Xinwei, et al. Study on the Development Approach of GF-7 High Precision Optical Stereo Mapping Satellite[J]. Spacecraft Recovery & Remote Sensing, 2020, 41(2): 17-28. (in Chinese)

[4] 李國元, 唐新明. 資源三號(02)星激光測高精度分析與驗證[J]. 測繪學報, 2017, 46(12): 1939-1949.

LI Guoyuan, TANG Xinming. Analysis and Validation of ZY-3 02 Satellite Laser Altimetry Data[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(12): 1939-1949. (in Chinese)

[5] 曹海翊, 劉付強, 趙晨光, 等. 高分辨率立體測繪衛(wèi)星技術研究[J]. 遙感學報, 2021, 25(7): 1400-1410.

CAO Haiyi, LIU Fuqiang, ZHAO Chenguang, et al. The Study of High Resolution Stereo Mapping Satellite[J]. National Remote Sensing Bulletin, 2021, 25(7): 1400-1410. (in Chinese)

[6] 張新偉, 賀濤, 趙晨光, 等. 高分七號衛(wèi)星測繪體制與性能評估[J]. 航天器工程, 2020, 29(3): 1-11.

ZHANG Xinwei, HE Tao, ZHAO Chenguang, et al. Demonstration of Surveying and Mapping System and Performance Evaluation of GF-7 Satellite[J]. Spacecraft Engineering, 2020, 29(3): 1-11. (in Chinese)

[7] 曹寧. 以DEM為高程輔助的弱交會衛(wèi)星影像的區(qū)域網(wǎng)平差[D]. 蘭州: 蘭州交通大學, 2014.

CAO Ning. DEM-aided Block Adjustment for Satellite Images with Weak Convergence Geometry[D]. Lanzhou: Lanzhou Jiaotong University, 2014. (in Chinese)

[8] 曹輝, 陶鵬杰, 李海鴻, 等. DEM約束的衛(wèi)星影像定位法[J]. 測繪學報, 2020, 49(1): 79-91.

CAO Hui, TAO Pengjie, LI Haihong, et al. Using DEM as Full Controls in Block Adjustment of Satellite Imagery[J]. Acta Geodaetica et Cartographica Sinica, 2020, 49(1): 79-91. (in Chinese)

[9] POLI D. General Model for Airborne and Spaceborne Linear Array Sensors[C]//ISPRS Commission I Symposium—— Integrating Remote Sensing at the Global, Regional and Local Scale, November 10-15, 2002, Denver, CO (USA). ISPRS, 2002: 177-182.

[10] KIM T, KIM H, RHEE S, et al. Investigation of Physical Sensor Models for Modeling SPOT 3 Orbits[J]. The Photogrammetric Record, 22(119): 257-273.

[11] 余俊鵬. 高分辨率衛(wèi)星遙感影像的精確幾何定位[D]. 武漢: 武漢大學, 2009.

YU Junpeng. Accurate Geometric Positioning of High Resolution Satellite Remote Sensing Imagery[D]. Wuhan: Wuhan University, 2009. (in Chinese)

[12] 賈博, 姜挺, 張銳, 等. 基于定向片模型的SPOT-5遙感影像自檢校光束法平差[J]. 測繪科學, 2014, 39(9): 3-10.

JIA Bo, JIANG Ting, ZHANG Rui, et al. Self-calibration Bundle Adjustment of SPOT-5 Remote-sensing Imagery Based on Orientation Image Model[J]. Science of Surveying and Mapping, 2014, 39(9): 3-10. (in Chinese)

[13] TOUTIN T. Spatiotriangulation with Multisensor VIR/SAR Images[J]. IEEE Transactions on Geoscience and Remote Sensing, 2004, 42(10): 2096-2103.

[14] TOUTIN T. Generation of DSMs from SPOT-5 In-track HRS and Across-track HRG Stereo Data using Spatiotriangulation and Autocalibration[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2006, 60(3): 170-181.

[15] 張過. 缺少控制點的高分辨率衛(wèi)星遙感影像幾何糾正[D]. 武漢: 武漢大學, 2005.

ZHANG Guo. Rectification for High Resolution Remote Sensing Image Under Lack of Ground Control Points[J]. Wuhan: Wuhan University, 2005. (in Chinese)

[16] 張過, 厲芳婷, 江萬壽, 等. 推掃式光學衛(wèi)星影像系統(tǒng)幾何校正產(chǎn)品的3維幾何模型及定向算法研究[J]. 測繪學報, 2010, 39(1): 34-38.

ZHANG Guo, LI Fangting, JIANG Wanshou, et al. Study of Three-dimensional Geometric Model and Orientation Algorithms for Systemic Geometric Correction Product of Pushbroom Optical Satellite Image[J]. Acta Geodaetica et Cartographica Sinica, 2010, 39(1): 34-38. (in Chinese)

[17] XIE J, TANG X, MO F, et al. In-Orbit Geometric Calibration and Experimental Verification of the ZY3-02 Laser Altimeter[J]. The Photogrammetric Record, 2018, 33(163): 341-362.

[18] 張過, 李少寧, 黃文超, 等. 資源三號02星對地激光測高系統(tǒng)幾何檢校及驗證[J]. 武漢大學學報(信息科學版), 2017, 42(11): 1589-1596.

ZHANG Guo, LI Shaoning, HUANG Wenchao, et al. Calibration and Validation of ZY3-02 Satellite Laser Altimeter System[J]. Geomatics and Information Science of Wuhan University, 2017, 42(11): 1589-1596. (in Chinese)

[19] LUTHCKE S B, STONEKING E. Spaceborne Laser-altimeter-pointing Bias Calibration from Range Residual Analysis[J]. Journal of Spacecraft and Rockets, 2000, 37(3): 374-384.

[20] MAGRUDER L A, WEBB C E, URBAN T J, et al. ICESat Altimetry Data Product Verification at White Sands Space Harbor[J]. IEEE Transactions on Geosciences and Remote Sensing, 2007, 45(1): 147-155.

[21] 田世強. 國產(chǎn)星載激光雷達在軌幾何檢校與驗證技術研究[D]. 西安: 長安大學, 2017.

TIAN Shiqiang. Geometric Calibration and Validation of Domestic Space-Borne Laser Altimeter[D]. Xi’an: Chang’an University, 2017. (in Chinese)

[22] 岳春宇, 鄭永超, 陶宇亮. 星載激光測高儀輔助衛(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 & Remote Sensing, 2013, 34(4): 71-76. (in Chinese)

[23] 唐新明, 劉昌儒, 張恒, 等. 高分七號衛(wèi)星立體影像與激光測高數(shù)據(jù)聯(lián)合區(qū)域網(wǎng)平差[J]. 武漢大學學報(信息科學版), 2021, 46(10): 1423-1430.

TANG Xinming, LIU Changru, ZHANG Heng, et a1. GF-7 Satellite Stereo Images Block Adjustment Assisted with Laser Altimetry Data[J]. Geomatics and Information Science of Wuhan University, 2021, 46(10): 1423-1430. (in Chinese)

[24] PI Y, LI X, YANG B. Global Iterative Geometric Calibration of a Linear Optical Satellite Based on Sparse GCPs[J]. IEEE Transactions on Geoscience and Remote Sensing, 2020, 58(1): 436-446.

Positioning without Ground Control Points of Vegetation Measurement System on Terrestrial Ecosystem Carbon Inventory Satellite

YUE Chunyu ZHOU Nan YANG Shuqi HE Hongyan

（Beijing Institute of Space Mechanics & Electricity, Beijing 100094, China）

The Vegetation Measurement System on Terrestrial Ecosystem Carbon Inventory Satellite includes a multi-beam LiDAR and a 0° camera, the high accuracy distance obtained by multi-beam LiDAR can be used as generalized elevation control point for the 0°camera image positioning adjustment without control points, when they are working at the same time in-orbit. In this paper, the positioning combined adjustment model is based on the geometric model of multi-beam LiDAR and the 0° camera, with the LiDAR control points corrected by public coarse grid DEM. The location accuracy of 0° camera is improved in the simulation experiment from 613.05m to 15.35m in the best case. And the number and the distribution of LiDAR control points are also analysed, when the LiDAR control points are evenly set in the upper, middle, and lower three lines of one scene image, the location accuracy is the highest.

vegetation measurement system; positioning without ground control points; combined adjustment; Terrestrial Ecosystem Carbon Inventory Satellite

TP751

A

1009-8518(2022)06-0050-09

10.3969/j.issn.1009-8518.2022.06.005

2022-10-17

國家重大科技專項工程

岳春宇, 周楠, 楊舒琪, 等. 陸地生態(tài)系統(tǒng)碳監(jiān)測衛(wèi)星植被測量子系統(tǒng)無控定位方法[J]. 航天返回與遙感, 2022, 43(6): 50-58.

YUE Chunyu, ZHOU Nan, YANG Shuqi, et al. Positioning without Ground Control Points of Vegetation Measurement System on Terrestrial Ecosystem Carbon Inventory Satellite[J]. Spacecraft Recovery & Remote Sensing, 2022, 43(6): 50-58. (in Chinese)

岳春宇，男，1983年生，2012年獲武漢大學攝影測量與遙感專業(yè)博士學位，高級工程師。主要研究方向航天攝影測量與遙感。E-mail：ycy1893@163.com。

（編輯：龐冰）