留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

资源三号03星激光测高数据处理与复合测绘应用

李国元 唐新明 周平 陈继溢 刘诏 窦显辉 周晓青 王霞

李国元, 唐新明, 周平, 陈继溢, 刘诏, 窦显辉, 周晓青, 王霞. 资源三号03星激光测高数据处理与复合测绘应用[J]. 红外与激光工程, 2022, 51(5): 20210356. doi: 10.3788/IRLA20210356
引用本文: 李国元, 唐新明, 周平, 陈继溢, 刘诏, 窦显辉, 周晓青, 王霞. 资源三号03星激光测高数据处理与复合测绘应用[J]. 红外与激光工程, 2022, 51(5): 20210356. doi: 10.3788/IRLA20210356
Li Guoyuan, Tang Xinming, Zhou Ping, Chen Jiyi, Liu Zhao, Dou Xianhui, Zhou Xiaoqing, Wang Xia. Laser altimetry data processing and combined surveying application of ZY3-03 satellite[J]. Infrared and Laser Engineering, 2022, 51(5): 20210356. doi: 10.3788/IRLA20210356
Citation: Li Guoyuan, Tang Xinming, Zhou Ping, Chen Jiyi, Liu Zhao, Dou Xianhui, Zhou Xiaoqing, Wang Xia. Laser altimetry data processing and combined surveying application of ZY3-03 satellite[J]. Infrared and Laser Engineering, 2022, 51(5): 20210356. doi: 10.3788/IRLA20210356

资源三号03星激光测高数据处理与复合测绘应用

doi: 10.3788/IRLA20210356
基金项目: 国家自然科学基金(41971425, 41871382)
详细信息
    作者简介:

    李国元,男,副研究员,博士,主要从事激光测高卫星、水资源遥感等方面的研究

    通讯作者: 唐新明,男,研究员,主要从事自然资源陆地卫星总体设计及遥感应用的相关研究。
  • 中图分类号: P237

Laser altimetry data processing and combined surveying application of ZY3-03 satellite

  • 摘要: 资源三号03星是自然资源部主持建造的用于1∶50 000立体测图的陆地遥感业务卫星,该星装备了业务化的激光测高仪,主要用于获取高精度高程控制点。论文针对资源三号03星激光测高数据,研究了标准化测绘处理流程和高程控制点提取方法,在内蒙古苏尼特右旗和江苏苏州开展了精度验证,并选择黑龙江和河北两个实验区开展了复合测绘应用验证。精度验证结果表明,资源三号03星激光点在内蒙古苏尼特右旗平坦区域高程精度为(0.051±0.232) m,在江苏苏州城市建成区的激光点总体精度为(0.414±6.213) m,经高程控制点提取和质量标记后的激光点高程误差为(−0.526±0.624) m,能满足1∶50 000测图高程控制需求。复合测绘应用表明,利用资源三号03星激光高程控制点,立体影像高程精度在黑龙江平坦地区能从5.27 m提高到2.58 m,河北太行山区能从11.25 m提高到4.45 m;无论是平地还是山区,资源三号03星激光高程控制点均能有效提高立体影像的高程精度并满足1∶50 000测图需求。
  • 图  1  资源三号03星激光测高仪阈值探测测距原理图

    Figure  1.  Diagram of ZY3-03 satellite laser altimeter based on the threshold detection and ranging

    图  2  资源三号03星激光测高数据处理流程图

    Figure  2.  The workflow of ZY3-03 satellite laser altimetry data processing

    图  3  资源三号03星激光点位分布示意图

    Figure  3.  Distribution of ZY3-03 satellite laser altimetry points all over the world

    图  4  RTK-GPS实际测量验证示意图

    Figure  4.  Illustration of RTK-GPS for validation in filed surveying

    图  5  实验区复合测绘数据和外业检查点分布示意图

    Figure  5.  Distribution of combined surveying data and field check points in experimental regions

    表  1  资源三号02/03星激光测高仪参数对比表

    Table  1.   Parameters comparison between ZY3-02 and ZY3-03 satellite laser altimeter

    ParameterZY3-02[15]ZY3-03
    Number of beamsOnly oneTwo, one as backup
    Frequency/Hz22
    Central wavelength/nm1 064.41 064.4
    Pulse energy/mJ200150
    Pulse width/ns75
    Divergence angle/μrad15090
    Ground footprint size75 m@500 km45 m@500 km
    Pointing change/μrad≤10≤5
    Receiving aperture/mm220260
    Total transmitted pulse>1×108
    Lifetime8 years
    Power consumption/W≤50≤80
    Weight/kg≤40≤95
    Terrain adaptability/(°)≤2≤5
    Ranging accuracy/m≤1≤0.3
    下载: 导出CSV

    表  2  内蒙古苏尼特右旗资源三号03星激光点实地验证结果表(单位:m)

    Table  2.   Validation result of ZY3-03 laser points on Sunid Right Banner in Inner Mongolia (Unit: m)

    DateIDElevationReference elevationError
    2020-09-08
    86912108 1 007.364 1 007.218 0.146
    86912110 1 052.404 1 052.183 0.221
    2020-09-23 92095212 1 073.461 1 073.316 0.145
    92095214 1 079.237 1 079.174 0.063
    92095216 1 074.751 1 074.405 0.346
    92095220 1 086.253 1 086.500 −0.247
    92095222 1 094.673 1 094.992 −0.319
    Statistical result Max 0.346
    Min 0.051
    RMSE 0.232
    下载: 导出CSV

    表  3  2020年10月10日的1183轨激光测高精度(单位:m)

    Table  3.   Laser elevation accuracy of the No.1183 orbit on 10th Oct, 2020 (Unit: m)

    IDECP flagElevationReference elevationErrorRemarks
    97960914 1 12.453 23.52 −11.067 Roof
    97960916 1 12.054 12.97 −0.916 Near the building
    97960920 2 20.576 14.66 5.916 Road with vegetation
    97960924 1 12.677 13.36 −0.683 Cropland
    97960926 1 10.524 10.69 −0.166 Cropland
    97960928 2 14.608 19.92 −5.312 Building and tree
    97960930 1 11.460 12.37 −0.910 Cropland
    97960932 1 9.756 9.83 −0.074 Paddy field
    97960934 2 30.476 31.19 −0.714 Roof
    97960936 1 9.995 10.74 −0.745 Open space
    97960938 3 14.560 11.75 2.810 Vegetation
    97960940 2 11.699 9.97 1.729 Trees
    97960942 1 10.525 11.15 −0.625 Cropland
    97960944 2 31.953 12.52 19.433 Construction site
    97960950 2 116.315 118.29 −1.975 Building
    97960952 1 30.225 30.31 −0.085 Grassland
    下载: 导出CSV

    表  4  实验区不同平差实验结果统计表

    Table  4.   Statistics result of different adjustment in the experimental regions

    Region
    Group
    Laser elevation control point
    Horizontal control point
    Ground control point
    Check point
    RMS of check point/m
    HorizontalElevation
    Taihang mountainA00011526.9911.25
    B0091065.983.32
    C60011527.434.45
    D6901065.974.31
    Heilongjiang flatA00027015.265.27
    B0092564.992.41
    C910027015.352.58
    D91902565.042.56
    下载: 导出CSV
  • [1] 唐新明, 李国元. 对地观测卫星激光测高数据处理方法与工程实践[M]. 北京: 科学出版社, 2019.

    Tang Xinming, Li Guoyuan. Earth Observing Satellite Laser Altimeter Data Processing Method and Engineer Practice[M]. Beijing: Science Press, 2019. (in Chinese)
    [2] Shu Rong, Huang Genghua, Kong Wei. Development and review of space-based laser altimetry technology [J]. Infrared and Laser Engineering, 2020, 49(11): 20201047. (in Chinese) doi:  10.3788/IRLA20201047
    [3] Schutz B E, Zwally H, Shuman C, et al. Overview of the ICESat mission [J]. Geophysical Research Letters, 2005, 32(21): 97-116.
    [4] Thomas A N, Anthony J M, Thorsten M, et al. The ice, cloud, and land elevation satellite-2 mission: A global geolocated photon product derived from the advanced topographic laser altimeter system [J]. Remote Sensing of Environment, 2019, 233: 111325. doi:  10.1016/j.rse.2019.111325
    [5] 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)
    [6] Tang Xinming, Chen Jiyi, Li Guoyuan, et al. Error analysis and preliminary pointing angle calibration of laser altimeter on Ziyuan-3 02 satellite [J]. Geomatics and Information Science of Wuhan University, 2018, 43(11): 1611-1619. (in Chinese)
    [7] Zhang Guo, Li Shaoning, Huang Wenchao, et al. Geometric 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)
    [8] Tang Xinming, Gao Xiaoming, Cao Haiyi, et al. The China ZY3-03 mission surveying and mapping technology for high-resolution remote sensing satellites [J]. IEEE Geoscience and Remote Sensing Magazine, 2020, 8(3): 8-17. doi:  10.1109/MGRS.2019.2929770
    [9] Li Guoyuan, Yao Jiaqi, Zhao Yiming, et al. Progress and prospect of atmospheric scattering correction for laser altimetry satellite [J]. Infrared and Laser Engineering, 2020, 49(11): 20200234. (in Chinese)
    [10] Yue Chunyu, Zheng Yongchao, Xing Yanqiu, et al. Technical and application development study of space-borne LiDAR in forestry remote sensing [J]. Infrared and Laser Engineering, 2020, 49(11): 20200235. doi:  10.3788/IRLA20200235
    [11] González J H, Bachmann M, Scheiber R, et al. Definition of ICESat selection criteria for their use as height references for TanDEM-X [J]. IEEE Transactions on Geoscience and Remote Sensing, 2010, 48(6): 2750-2757. doi:  10.1109/TGRS.2010.2041355
    [12] Li Guoyuan, Tang Xinming, Zhang Chongyang, et al. Multi-criteria constraint algorithm for selecting ICESat/GLAS data as elevation control points [J]. Journal of Remote Sensing, 2017, 21(1): 96-104. (in Chinese) doi:  10.11834/jrs.20175269
    [13] Li Guoyuan, Tang Xinming, Gao Xiaoming, et al. Integration of ZY3-02 satellite laser altimetry data and stereo images for high accuracy mapping [J]. Photogrammetric Engineering and Remote Sensing, 2018, 84(9): 569-578. doi:  10.14358/PERS.84.9.569
    [14] Fang Yong, Cao Bincai, Gao Li, et al. Development and application of lidar mapping satellite [J]. Infrared and Laser Engineering, 2020, 49(11): 20201044. doi:  10.3788/IRLA20201044
    [15] Li Guoyuan, Gao Xiaoming, Chen Jiyi, et al. Data quality analysis of ZY-3 02 satellite laser altimeter [J]. Journal of Remote Sensing, 2019, 23(6): 1159-1166. (in Chinese)
    [16] Tang Xinming, Li Guoyuan, Gao Xiaoming, et al. The rigorous geometric model of satellite laser altimeter and preliminarily accuracy validation [J]. Acta Geodaetica et Cartographica Sinica, 2016, 45(10): 1182-1191. (in Chinese) doi:  10.11947/j.AGCS.2016.20150357
    [17] Gardner C S. Ranging performance of satellite laser altimeters [J]. IEEE Transactions on Geoscience & Remote Sensing, 1992, 30(5): 1061-1072.
    [18] Li Guoyuan, Tang Xinming, Fan Wenfeng, et al. On -orbit geometric calibration of spaceborne laser altimeter based on ground infrared detector [J]. Infrared and Laser Engineering, 2017, 46(11): 1117004. (in Chinese) doi:  10.3788/IRLA201746.1117004
    [19] GB/T 12341—2008.1∶25000 1∶50000 1∶100000地形图航空摄影测量外业规范[S].

    GB/T 12341—2008. Specifications for aerophotogrammetric field work of 1∶25000 1∶50000 1∶100000 topographic maps [S], 2008. (in Chinese)
  • [1] 朱天豪, 周辉, 石岩, 张千胤.  星载激光测高仪多模式回波参数提取方法(特邀) . 红外与激光工程, 2022, 51(1): 20210836-1-20210836-11. doi: 10.3788/IRLA20210836
    [2] 么嘉棋, 李国元, 陈继溢, 周晓青, 国爱燕, 黄庚华, 唐新明, 艾波.  高分七号卫星激光测高仪光斑质心位置变化分析 . 红外与激光工程, 2021, 50(S2): 20210539-1-20210539-10. doi: 10.3788/IRLA20210539
    [3] 谢俊峰, 刘仁, 王宗伟, 莫凡, 窦显辉.  高分七号星载激光测高仪在轨几何检校与精度评估 . 红外与激光工程, 2021, 50(8): 20200396-1-20200396-11. doi: 10.3788/IRLA20200396
    [4] 汪自军, 张扬, 刘东, 王晓波, 袁金如, 潘超, 赵一鸣, 韩晓爽, 周雨迪, 刘群, 王成.  新型多波束陆-海激光雷达探测卫星技术发展研究 . 红外与激光工程, 2021, 50(7): 20211041-1-20211041-11. doi: 10.3788/IRLA20211041
    [5] 郭金权, 李国元, 左志强, 张宁, 裴亮, 卢刚.  高分七号卫星激光测高仪全波形数据质量及特征分析 . 红外与激光工程, 2020, 49(S2): 20200387-20200387. doi: 10.3788/IRLA20200387
    [6] 左志强, 唐新明, 李国元, 李松.  GF-7星载激光测高仪全波形自适应高斯滤波 . 红外与激光工程, 2020, 49(11): 20200251-1-20200251-11. doi: 10.3788/IRLA20200251
    [7] 程勇, 刘旭, 谭朝勇, 王赛, 韦尚方.  星载激光测高仪固体激光器技术研究与发展 . 红外与激光工程, 2020, 49(11): 20201046-1-20201046-10. doi: 10.3788/IRLA20201046
    [8] 高小明, 李国元, 郭金权, 欧阳斯达, 樊文锋.  发展三极观测激光测高卫星的探讨 . 红外与激光工程, 2020, 49(11): 20200233-1-20200233-7. doi: 10.3788/IRLA20200233
    [9] 舒嵘, 黄庚华, 孔伟.  空间激光测高技术发展及展望 . 红外与激光工程, 2020, 49(11): 20201047-1-20201047-10. doi: 10.3788/IRLA20201047
    [10] 张鑫磊, 邢帅, 徐青, 张国平, 李鹏程, 焦麟, 刘宸博.  ATLAS数据与资源三号02星影像联合区域网平差 . 红外与激光工程, 2020, 49(S2): 20200194-20200194. doi: 10.3788/IRLA20200194
    [11] 刘向锋, 黄庚华, 张志杰, 王凤香, 舒嵘.  高分七号激光测高中全波形回波数据的EMD降噪 . 红外与激光工程, 2020, 49(11): 20200261-1-20200261-10. doi: 10.3788/IRLA20200261
    [12] 张文豪, 李松, 张智宇, 刘芮, 马跃.  利用波形匹配实现卫星激光测高脚点精确定位的方法 . 红外与激光工程, 2018, 47(11): 1117007-1117007(8). doi: 10.3788/IRLA201847.1117007
    [13] 李少辉, 周辉, 倪国强.  基于星载激光测高仪多模式回波的激光测距修正值分析 . 红外与激光工程, 2017, 46(10): 1006001-1006001(8). doi: 10.3788/IRLA201759.1006001
    [14] 杨驰, 胡文怡, 罗敏, 孙莹莹, 周辉.  椭圆高斯足印对星载激光测高仪测距值及其误差的影响 . 红外与激光工程, 2016, 45(7): 717003-0717003(7). doi: 10.3788/IRLA201645.0717003
    [15] 周辉, 李松, 王良训, 郑国兴.  噪声对星载激光测高仪测距误差的影响 . 红外与激光工程, 2015, 44(8): 2256-2261.
    [16] 王良训, 周辉, 李子乐, 刘国根, 王虹, 王雅培.  面向星载激光测高仪的陆地目标响应函数仿真 . 红外与激光工程, 2015, 44(11): 3424-3430.
    [17] 马跃, 阳凡林, 王明伟, 李松, 翁寅侃.  利用GLAS激光测高仪计算格陵兰冰盖高程变化 . 红外与激光工程, 2015, 44(12): 3565-3569.
    [18] 马跃, 阳凡林, 卢秀山, 冯成凯, 李松.  对地观测星载激光测高系统高程误差分析 . 红外与激光工程, 2015, 44(3): 1042-1047.
    [19] 王明明, 朱枫, 郝颖明.  测地自旋图院三维物体局部形状描述符 . 红外与激光工程, 2014, 43(8): 2709-2714.
    [20] 马跃, 李松, 翁寅侃, 周辉.  星载激光测高仪大气干项延迟校正 . 红外与激光工程, 2013, 42(4): 909-914.
  • 加载中
图(5) / 表(4)
计量
  • 文章访问数:  23
  • HTML全文浏览量:  4
  • PDF下载量:  6
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-11-10
  • 修回日期:  2021-12-20
  • 刊出日期:  2022-06-08

资源三号03星激光测高数据处理与复合测绘应用

doi: 10.3788/IRLA20210356
    作者简介:

    李国元,男,副研究员,博士,主要从事激光测高卫星、水资源遥感等方面的研究

    通讯作者: 唐新明,男,研究员,主要从事自然资源陆地卫星总体设计及遥感应用的相关研究。
基金项目:  国家自然科学基金(41971425, 41871382)
  • 中图分类号: P237

摘要: 资源三号03星是自然资源部主持建造的用于1∶50 000立体测图的陆地遥感业务卫星,该星装备了业务化的激光测高仪,主要用于获取高精度高程控制点。论文针对资源三号03星激光测高数据,研究了标准化测绘处理流程和高程控制点提取方法,在内蒙古苏尼特右旗和江苏苏州开展了精度验证,并选择黑龙江和河北两个实验区开展了复合测绘应用验证。精度验证结果表明,资源三号03星激光点在内蒙古苏尼特右旗平坦区域高程精度为(0.051±0.232) m,在江苏苏州城市建成区的激光点总体精度为(0.414±6.213) m,经高程控制点提取和质量标记后的激光点高程误差为(−0.526±0.624) m,能满足1∶50 000测图高程控制需求。复合测绘应用表明,利用资源三号03星激光高程控制点,立体影像高程精度在黑龙江平坦地区能从5.27 m提高到2.58 m,河北太行山区能从11.25 m提高到4.45 m;无论是平地还是山区,资源三号03星激光高程控制点均能有效提高立体影像的高程精度并满足1∶50 000测图需求。

English Abstract

    • 基础地理信息数据是国家重要的基础性、战略性信息资源,卫星测绘是支撑国民经济社会建设高质量发展的重要手段。随着激光技术的快速发展,卫星激光测高已经成为卫星测绘的重要手段,特别是在全球高程控制点获取、境外地区高精度立体测绘方面发挥了越来越重要的作用[1-2]。美国先后通过航天飞机SLA-1/2、ICESat-1/2获取了全球大量的高程控制点数据,支撑了其全球战略实施及对外政治军事行动[3-4]。我国先后在资源三号02星、高分七号01星上装备了激光测高仪,通过激光与影像复合测绘的方式分别满足1:50 000、1:10 000无地面控制点的立体测图高程精度要求。其中资源三号02星上实现了我国对地观测领域卫星激光测高“从无到有”的重要突破,为卫星激光测高的快速发展奠定了重要技术基础[5-7]

      资源三号03星基本继承资源三号02星[8],其中激光测高仪升级为业务化载荷,主要用于获取高精度的高程控制点,然而受大气[9]、地形[10]、地物等因素影响,并不是每个激光点均能用作高程控制,自动筛选高精度、高可靠的激光高程控制点非常关键。参考文献[11-12]对ICESat/GLAS的高程控制点提取方法进行了较系统的研究,但资源三号03星激光测高仪由于受体积、口径、功耗等条件限制,没有采用全波形记录模式,仅记录了回波离散阈值,这对高精度高程控制点提取提出了一定挑战。将激光高程控制点参与立体影像区域网平差处理,是激光与影像复合测绘的重要步骤,能有效提高立体影像在无地面控制点条件的高程测量精度,在境外地区测图有非常重要的应用价值[13-14]

      文章介绍了资源三号03星的基本参数,并将03星的激光测高仪与02星进行了对比,研究了激光测高数据处理及高程控制点提取方法,对03星立体影像与激光复合测绘进行了应用分析。

    • 资源三号是我国首个高分辨率民用立体测图卫星星座,01星、02星分别于2012年、2016年成功发射,资源三号03星于2020年7月25日在太原卫星发射中心成功发射。除继承资源三号02星几何分辨率分别为2.1 m正视、2.7 m前后视的三线阵相机和5.8 m多光谱相机载荷外,资源三号03星还实现了国产星敏感器以及激光测高仪的业务化,同时将卫星设计寿命从原来的5年提高为8年。卫星运行于回归周期为59天的太阳同步回归圆轨道,通过三星组网,重访周期从3天缩短到1天,保证我国高分辨率立体测绘数据的长期稳定获取,形成全球领先的业务化立体观测能力[8]

    • 资源三号03星激光测高仪在02星试验性载荷的基础上做了适当的优化改进,如表1所示。具体包括:设计寿命从试验性载荷变为8年,同时增加了一个备份激光器,激光发散角从150 μrad缩小为90 μrad,对应地面光斑大小约45 m,采用前后沿鉴别的方式,具有多阈值探测能力。

      表 1  资源三号02/03星激光测高仪参数对比表

      Table 1.  Parameters comparison between ZY3-02 and ZY3-03 satellite laser altimeter

      ParameterZY3-02[15]ZY3-03
      Number of beamsOnly oneTwo, one as backup
      Frequency/Hz22
      Central wavelength/nm1 064.41 064.4
      Pulse energy/mJ200150
      Pulse width/ns75
      Divergence angle/μrad15090
      Ground footprint size75 m@500 km45 m@500 km
      Pointing change/μrad≤10≤5
      Receiving aperture/mm220260
      Total transmitted pulse>1×108
      Lifetime8 years
      Power consumption/W≤50≤80
      Weight/kg≤40≤95
      Terrain adaptability/(°)≤2≤5
      Ranging accuracy/m≤1≤0.3

      资源三号03星激光测高仪的测距原理如图1所示。受硬件资源等限制,03星激光测高仪的硬件误差主要包括发射波形时刻鉴别误差、回波时刻鉴别误差、晶振计时误差和TDC (Time Distance Calculator)时间间隔测量误差,平地理论测距精度优于0.3 m。

      图  1  资源三号03星激光测高仪阈值探测测距原理图

      Figure 1.  Diagram of ZY3-03 satellite laser altimeter based on the threshold detection and ranging

      资源三号03星激光测高仪回波有四个通道,包括三个宽带通道和一个窄带通道。其中宽带通道又分为三个子通道,分别为高阈值TH1子通道,低阈值TH2子通道和噪声监测阈值TH4子通道。窄带通道为TH3通道。噪声监测阈值TH4子通道类似于资源三号02星设计状态,仅记录门内第一个回波的前沿,输出测距值,同时可以输出距离内的过阈值脉冲个数,用于在轨检测噪声,指导阈值设置;其他TH1,TH2和TH3三个通道均记录门内回波脉冲的前沿和后沿,并进行输出,每个通道均具有记录四次回波脉冲的能力。

      根据主波时刻和回波时刻计算激光测高仪的飞行时间为:

      $${T_{tof}} = {T_{echo}} - {T_{tx}} - \Delta {T_{THi}}$$ (1)

      式中:$\Delta {T_{THi}}$为各通道延迟常数。

      对大气延迟影响$\Delta R$等进行改正后,激光的测距值为:

      $$R = \frac{c}{2}{T_{tof}}{\rm{ - }}\Delta R$$ (2)
    • 对地观测卫星激光测高数据在获取过程中,受到激光指向、卫星轨道和姿态、大气、地球物理潮汐、大光斑内复杂地形地物等多种因素的影响,标准化的测绘处理是确保资源三号03星高精度激光测高的重要前提。

    • 图2所示,根据卫星激光测高严密几何模型进行激光地面足印几何定位,结合卫星上搭载的GPS定位仪器和星敏感器获得的卫星位置和姿态信息,解算激光足印精确的三维坐标[16]。卫星激光测高严密几何定位模型如下式所示:

      图  2  资源三号03星激光测高数据处理流程图

      Figure 2.  The workflow of ZY3-03 satellite laser altimetry data processing

      $$\begin{split}& \left( {\begin{array}{*{20}{c}} {{X_P}} \\ {{Y_P}} \\ {{Z_P}} \end{array}} \right) = \left( {\begin{array}{*{20}{c}} {{X_S}} \\ {{Y_S}} \\ {{Z_S}} \end{array}} \right){\rm{ + }}\;{{R}}_{J2000}^{WGS84}{{R}}_{Body}^{J2000}\\& \left[ {\left( {\begin{array}{*{20}{c}} {Lx} \\ {Ly} \\ {L{\textit{z}}} \end{array}} \right) - \left( {\begin{array}{*{20}{c}} {Dx} \\ {Dy} \\ {D{\textit{z}}} \end{array}} \right) + \left. {\left( {\begin{array}{*{20}{c}} {\rho \sin \theta \cos \alpha } \\ {\rho \sin \theta \sin \alpha } \\ { - \rho \cos \theta } \end{array}} \right)} \right]} \right.\end{split}$$ (3)

      式中:${\left[ {{X_P}\;{Y_P}\;{Z_P}} \right]^{\rm{T}}}$为激光足印点在WGS84坐标系下的三维坐标;${\left[ {{X_S}\;{Y_S}\;{Z_S}} \right]^{\rm{T}}}$为该激光足印点对应时刻卫星在WGS84坐标系的三维坐标;${{R}}_{Body}^{J2000}$为星敏感器测定的激光足印点对应时刻卫星本体坐标系与J2000坐标系之间的旋转矩阵;${{R}}_{WGS84}^{J2000}$为激光足印点对应时刻WGS84坐标系与J2000坐标系之间的旋转矩阵;${\left[ {Dx\;Dy\;Dz} \right]^{\rm{T}}}$为GPS的相位中心在卫星本体坐标系下的坐标;${\left[ {Lx\;Ly\;Lz} \right]^{\rm{T} }}$为激光参考中心在卫星本体坐标系下的坐标;$\; \rho $为激光所测距离值;$\theta $为激光指向与本体坐标系Z轴的负向夹角;$\alpha $为激光指向在XOY平面上的投影与X轴正向夹角。激光指向角也可以采用与卫星本体系三个轴的三个夹角进行描述,两者属于等价关系[1]

      截止到2021年4月30日,资源三号03星激光发射次数超过125万次,有效点数约72万个,占比约57.6%,全球分布情况如图3所示。资源三号03星激光测高标准产品目前已在自然资源部国土卫星遥感应用中心实现业务化处理。

      图  3  资源三号03星激光点位分布示意图

      Figure 3.  Distribution of ZY3-03 satellite laser altimetry points all over the world

    • 高程控制点提取是复合测绘应用的重要基础,具备全波形采样的激光高程控制点提取方法已得到有效应用[11,17],但资源三号03星不具备全波形采样功能,文中在借鉴资源三号02星试验性激光测高载荷处理经验的基础上[5],结合资源三号03星激光测高仪阈值探测的特点,针对性地提出了高程控制点ECP (Elevation control Point)提取及质量标记的方法,具体流程如下:

      (1)确保大气、潮汐、精密轨道和姿态等参数的有效性,对激光测高标准产品的完备性进行质量控制;

      (2)利用全球30 m的地表覆盖数据GLC30,将落在水面的激光点标记为10;

      (3)结合阈值探测和地形地物影响,结合波形阈值计算的脉宽,对没有落水的激光点,设定回波脉宽阈值为${\tau _s}$≤20 ns时,默认光斑内坡度小于2°,此时激光点标记为1;若20 ns≤${\tau _s}$≤40 ns,默认坡度介于2°~5°,质量标记为2,其余的标记为3;

      (4)进一步地,对标记为1和2的激光点,借鉴资源三号02星激光点的提取思想,采用全球30 m格网的AW3D30产品,计算激光点高程与光斑范围内AW3D30平均高的差值dh。如果高差|dh|≤2 m,则原标记不变;若高差2 m<|dh|≤5 m,则原标记为1的点变为2;若|dh|>5 m,则原标记为1和2的点均变为3。

      经筛选提取后,激光点标记为1的能用作高程控制点使用,标记为2的点仅推荐在标记为1的点较少的情况下参考使用,标记为3的误差较大不推荐使用,标记为10的代表该点落水。

    • 利用地面红外探测器是实现激光点绝对几何精度验证的最直观方法[18],但该方法成本较大。文中在实现激光点标准化测绘处理后,选择内蒙古苏尼特右旗和江苏苏州两个验证区,结合外业RTK-GPS和高精度机载LiDAR点云数据,对资源三号03星激光点精度进行验证和分析。

      (1)在轨定标区附近验证结果

      在内蒙古苏尼特右旗外场定标区附近,针对2020年9月8日、9月23日的资源三号03星实际激光落点位置,选择地形平坦、没有植被和地物影响的7个激光点,采用高精度的RTK-GPS测量方式,对每个激光点以约5 m的间隔测量其附近地面点的三维坐标,形成RTK-GPS点阵,如图4所示。

      图  4  RTK-GPS实际测量验证示意图

      Figure 4.  Illustration of RTK-GPS for validation in filed surveying

      以激光点平面位置为中心,外扩25 m统计光斑范围内RTK-GPS的平均高作为参考,统计激光点的高程误差,结果如表2所示。该区域激光测高精度较高,高程精度为(0.051±0.232) m。

      表 2  内蒙古苏尼特右旗资源三号03星激光点实地验证结果表(单位:m)

      Table 2.  Validation result of ZY3-03 laser points on Sunid Right Banner in Inner Mongolia (Unit: m)

      DateIDElevationReference elevationError
      2020-09-08
      86912108 1 007.364 1 007.218 0.146
      86912110 1 052.404 1 052.183 0.221
      2020-09-23 92095212 1 073.461 1 073.316 0.145
      92095214 1 079.237 1 079.174 0.063
      92095216 1 074.751 1 074.405 0.346
      92095220 1 086.253 1 086.500 −0.247
      92095222 1 094.673 1 094.992 −0.319
      Statistical result Max 0.346
      Min 0.051
      RMSE 0.232

      (2)附加ECP标记的江苏苏州验证结果

      2020年10月10日资源三号03星第1183轨经过江苏苏州,有16个激光点位于验证区内。参考数据为高精度的机载LiDAR点云数据,获取时间为2018年,高程精度优于0.15 m,高程基准统一为WGS84大地高。针对进行了高程控制点ECP标记的资源三号03星激光测高标准产品,以激光点平面位置为中心,外扩25 m统计光斑范围内LiDAR点云数据平均高,计算其与激光点高程误差,结果如表3所示。

      表 3  2020年10月10日的1183轨激光测高精度(单位:m)

      Table 3.  Laser elevation accuracy of the No.1183 orbit on 10th Oct, 2020 (Unit: m)

      IDECP flagElevationReference elevationErrorRemarks
      97960914 1 12.453 23.52 −11.067 Roof
      97960916 1 12.054 12.97 −0.916 Near the building
      97960920 2 20.576 14.66 5.916 Road with vegetation
      97960924 1 12.677 13.36 −0.683 Cropland
      97960926 1 10.524 10.69 −0.166 Cropland
      97960928 2 14.608 19.92 −5.312 Building and tree
      97960930 1 11.460 12.37 −0.910 Cropland
      97960932 1 9.756 9.83 −0.074 Paddy field
      97960934 2 30.476 31.19 −0.714 Roof
      97960936 1 9.995 10.74 −0.745 Open space
      97960938 3 14.560 11.75 2.810 Vegetation
      97960940 2 11.699 9.97 1.729 Trees
      97960942 1 10.525 11.15 −0.625 Cropland
      97960944 2 31.953 12.52 19.433 Construction site
      97960950 2 116.315 118.29 −1.975 Building
      97960952 1 30.225 30.31 −0.085 Grassland

      表3中可以看出,在农田、草地、裸地等区域,03星激光点高程精度较高,在居民区、林地等有建筑、树木的区域,因光斑内地形地物复杂,激光点的高程精度相对较差,这些区域的激光点应该尽量避免使用。该验证区16个资源三号03星激光点的总体高程精度为(0.414±6.213) m,经过高程控制点自动提取后,控制点标记为1的有9个,占比约52.9%,除其中点号为97960914的激光点因位于居民区的建筑物附近,高程误差达到11 m外,剩下的8个点高程误差为(−0.526±0.624) m,满足1:50 000立体测图高程控制点精度要求[19]

    • 为了提升资源三号03星三线阵立体影像的高程精度,利用资源三号03星光学立体影像和激光测高数据同平台获取的特点,开展了激光测高点与光学立体影像条带模式的联合测绘处理,评价利用激光测高点提升立体影像高程精度是否满足1∶50 000测图精度要求。选择了黑龙江试验区、河北太行山两个实验区,其中黑龙江实验区位于东北平原,面积约3万 km2,整个区域的地形起伏从海拔120~550 m,其中南部是平地地形,北部是丘陵和山地地形。河北太行山实验区位于河北省太行山地区,整个区域的地形起伏从海拔40~220 m,属于山地地形。

      2020年9月4日和11月4日资源三号03星第663轨和1552轨分别经过河北太行山、黑龙江试验区,同步获取了三线阵立体影像和激光点数据。为了验证试验精度,分别采用分布于试验区内的RTK-GPS测量的115个和270个高精度外业点作为影像检查点、并取其中的9个作为平面或平高控制点,所有点的大地基准为WGS84坐标系,高程基准为WGS84椭球高,平面和高程精度均优于0.2 m,点位分布如图5所示。

      图  5  实验区复合测绘数据和外业检查点分布示意图

      Figure 5.  Distribution of combined surveying data and field check points in experimental regions

      在区域网平差中,自动匹配的连接点中误差优于0.3 pixel。如表4所示,为了有效验证和对比激光与影像复合测绘的实验效果,每个试验区均开展了如下4组区域网平差实验:

      表 4  实验区不同平差实验结果统计表

      Table 4.  Statistics result of different adjustment in the experimental regions

      Region
      Group
      Laser elevation control point
      Horizontal control point
      Ground control point
      Check point
      RMS of check point/m
      HorizontalElevation
      Taihang mountainA00011526.9911.25
      B0091065.983.32
      C60011527.434.45
      D6901065.974.31
      Heilongjiang flatA00027015.265.27
      B0092564.992.41
      C910027015.352.58
      D91902565.042.56

      A: 无任何控制条件下光学立体影像自由网平差,验证立体影像的原始几何精度;

      B: 选取少量RTK-GPS点作为平高控制点,针对立体影像开展区域网平差,验证常规控制条件下立体影像的几何精度。

      C: 利用激光点作为高程控制点,开展立体影像与激光点联合区域网平差,验证利用激光测高数据提升光学立体影像的效果;

      D: 利用激光点作为高程控制,选取少量RTK-GPS点作为平面控制点,开展联合区域网平差,验证引入平面控制时,利用激光测高数据提升光学立体影像的效果。

    • 我国1∶50 000立体测图成果的平面精度要求为平地优于25 m、山地优于37.5 m,高程精度要求为平地优于3 m、山地优于5 m。

      对比实验组A和B可以发现,两个验证区在完全无控制的情况下,平面精度能达标,但高程无法直接满足要求,利用外业控制点后均能满足相关测图规范要求。实验组C结果表明,单独采用激光点作为高程控制,在没有其他控制数据的情况下,对立体影像的平面精度几乎没有影响,但高程精度有较大提升,其中河北太行山区从11.25 m提高到4.45 m,黑龙江地区从5.27 m提高到2.58 m,均能满足1∶50 000测图要求。

      对比实验组B和C,黑龙江地区基于激光点和高精度地面控制点的高程精度结果分别为2.58 m和2.41 m,河北太行山区两者的高程精度提升对比结果为4.45 m和3.32 m,说明平原地区激光点的效果更明显,基本可相当于外业控制点的水平。山区则因目前的阈值探测模式及地形限制,点的分布如图5(a)所示可能有一定欠缺,导致结果比外业控制点稍差,但是高程精度也已经满足我国1∶50 000山地区域高程精度要求。

      实验组D的结果则表明,在激光与光学立体影像联合区域网平差中加入适量的平面控制点后,在不影响激光点对立体影像高程精度提升效果的同时,可进一步提升试验影像的平面精度,使之完全满足我国1∶50 000比例尺立体测图精度要求。

    • 资源三号03星作为我国1∶50 000比例尺的民用高分辨率立体测图业务卫星,将在01/02星的基础上,继续重点保障新型基础测绘、全球测图、地理国情监测、实景三维中国等重大测绘工程建设。论文对资源三号03星的激光测高数据处理及精度验证进行了分析,同时对激光与立体影像复合测绘应用进行了探讨,结果表明:

      (1) 资源三号03星激光测高仪基本继承02星,虽然受卫星平台及硬件水平限制,但寿命指标以及精度水平仍有较大改善,从试验性载荷变为业务化载荷,从单阈值探测变为多阈值探测, 内蒙古定标场附近激光点高程验证精度为(0.051±0.232) m,优于02星的激光测高水平;

      (2)基于阈值约束的激光高程控制点提取方法能有效标记激光点质量,江苏苏州城市建成区的激光点总体精度为(0.414±6.213) m,经高程控制点提取和质量标记后的激光点高程误差为(−0.526±0.624) m,能满足1∶50 000立体测图高程控制精度要求;

      (3)无论平原、丘陵或山地地形区域,通过激光测高点与立体影像联合区域网平差处理后的立体影像,其高程精度均可显著提升至满足我国1∶50 000比例尺立体测图高程精度的水平,这表明资源三号03星的激光测高数据可以适用于全地形区域的立体影像高程精度提升。

      资源三号卫星星座已成为我国自然资源卫星遥感观测体系的重要组成部分,是自然资源季度遥感监测的主要数据源,资源三号03星成功交付以及三星顺利组网运行,将进一步提升1∶50 000比例尺立体测图及1∶2.5万比例尺地理信息更新能力,进一步提高山水林田湖草等自然资源调查监测能力,进一步有力支撑国家治理能力和治理体系现代化。

参考文献 (19)

目录

    /

    返回文章
    返回