文章快速检索  
  高级检索
资源三号02星激光测高精度分析与验证
李国元1,2,3, 唐新明1,2,3     
1. 国家测绘地理信息局卫星测绘应用中心, 北京 100048;
2. 武汉大学资源与环境科学学院, 湖北 武汉 430079;
3. 江苏省地理信息资源开发与利用协同创新中心, 江苏 南京 210023
摘要:资源三号02星搭载了我国首台对地观测的卫星激光测高试验性载荷,对该载荷的精度进行了理论分析,并采用多个区域进行了实际精度验证,同时对其在航天测绘中的应用进行了试验。资源三号02星激光测高仪在平坦地区(坡度≤2°)的理论高程精度为0.85 m、平面精度14.2 m。试验表明,资源三号02星激光测高仪获得的有效测高数据约占23.89%,检校场区域其高程精度为0.89 m,平面精度为14.76 m;华北地区高精度DSM地形数据验证其高程精度为1.09 m,内陆渤海海面上的激光高程精度为0.47 m。将激光足印点作为高程控制点时,在陕西渭南试验区能将资源三号02星立体影像无地面控制的高程精度从11.54 m提高到1.90 m。虽然资源三号02星激光测高仪为试验性载荷,但试验结果证实国产卫星激光测高数据能有效提高立体影像无地面控制的高程精度,在全球测图工程中具有推广应用价值,建议后续立体测图卫星搭载业务化应用的激光测高仪。
关键词:资源三号02星    卫星激光测高    精度分析    全球测图    
Analysis and Validation of ZY-3 02 Satellite Laser Altimetry Data
LI Guoyuan1,2,3, TANG Xinming1,2,3     
1. Satellite Surveying and Mapping Application Center, NASG, Beijing 100048, China;
2. School of Resource and Environmental Sciences, Wuhan University, Wuhan 430079, China;
3. Jiangsu Center for Collaborative Innovation in Geographical Information Resource Development and Application, Nanjing 210023, China
Foundation support: The National Natural Science Foundation of China (No. 41601505) The Special Fund for High Resolution Images Surveying and Mapping Application System of China (No. AH1601-8) The Surveying and Mapping Basic Research Program of National Administration of Surveying, Mapping and Geo-information (Nos. 2017KJ0204 2018KJ0204)
First author: LI Guoyuan(1984—), male, PhD, registered surveyor, majors in satellite laser altimetry data processing and application.E-mail:ligy@sasmac.cn
Corresponding author: TANG Xinming, E-mail:tangxinming99@qq.com
Abstract: ZY-3 02 satellite loaded with Chinese first earth observing satellite laser altimeter, and has been launched successfully on 30th May, 2016. In this paper, the theoretical accuracy of the laser altimeter is analyzed, and several experimental areas are used to verify the actual accuracy. At the same time, the application of the laser altimetry data in the field of space-borne photogrammetry is tested. The laser altimetry theoretical accuracy of ZY-3 02 satellite in the flat area (slope less than 2 degrees) is about 0.85 m and 14.2 m in the elevation and planimetry direction, respectively. The effective laser altimetry data account for about 23.89%, and near the calibration field the elevation accuracy is 0.89 m, and planimetry accuracy is about 14.76 m. Moreover, the verified elevation accuracy is 1.09 m in the North China by high precision DSM terrain data, and laser footprint points accuracy on the surface of the Bohai inland sea is about 0.47 m. When the laser foot print point is used as elevation control point, the elevation accuracy of the ZY-3 02 satellite stereo images in Shaanxi Weinan can be increased from 11.54 m to 1.90 m without GCPs. Although ZY3-02 satellite laser altimeter is just a test, the results proved that the domestic satellite laser altimetry data can effectively improve the stereo images without GCPs, which will be valuable in the global mapping project. It is suggest that operational laser altimeter equip on the next satellite of ZY-3 serials.
Key words: ZY-3 02 satellite     satellite laser altimetry     accuracy analysis     global mapping    

卫星激光测高具备主动获取全球地表及目标三维信息的能力,能为快速获取包括境外地区在内的三维控制点以及立体测图提供服务,同时在极地冰盖测量、植被高度及生物量估测、云高测量、海面高度测量以及全球气候监测等方面都可以发挥重要作用。美国从20世纪70年代就开始逐步发展卫星激光测高,之后在国际上一直独占鳌头,在月球、火星及水星等星体测绘方面大量采用了卫星激光测高技术[1-4]。而对地观测领域,美国在1996年和1997年通过航天飞机SLA(shuttle laser altimeter)先后两次搭载激光测高仪,成功获得南北纬57°范围内的全球高精度控制点信息[5];在2003年成功发射了ICESat(ice,cloud and land elevation satellite)卫星,通过搭载GLAS(geo-sciences laser altimeter system)地球科学激光测高系统开展了极地冰盖监测、海冰高程测量、森林生物量估算、全球陆地高程控制点获取等应用[6-7],在国际上形成广泛的影响。美国还计划在2018年发射ICESat-2/ATLAS(advanced topographic laser altimeter system),并在国际空间站上搭载GEDI(global ecosystems dynamics investigation)激光测高载荷,高精度激光地形测量卫星LIST(LiDAR surface topography)也正在预研之中[8-11]

资源三号02星于2016年5月30日成功发射,搭载了国内首台对地观测的试验性激光测高载荷,主要用于测试激光测高仪的功能和性能,探索地表高精度的高程控制点数据获取的可行性,以及采用该数据辅助提高光学卫星影像无地面控制立体测图精度的可能性。未来几年,用于1:10 000比例尺立体测图的国产高分七号卫星以及陆地生态系统碳监测卫星也均将搭载激光测高仪,开展全球高程控制点获取、树高和生物量等森林参数的估算等应用。在资源三号02星之前,由于缺乏相关国产卫星型号及数据,我国在对地观测卫星激光测高领域的发展相对较慢。因此,借助目前资源三号02星试验性激光测高载荷,开展数据处理与应用实践,能有效促进我国相关技术的发展,为后续的高分七号、陆地生态系统碳监测卫星等激光测高载荷提供参考。

在卫星激光测高数据处理方面,文献[12]对不同地形条件下的卫星激光测高精度进行了模拟仿真分析;文献[13-15]对资源三号02星激光测高仪在轨几何检校进行了研究和试验,其中文献[15]提出了一种“两步法”在轨几何检校方法。笔者在文献[16]中对卫星激光测高严密几何模型进行了研究,并对资源三号02星激光测高仪的基本参数做了介绍,对第一轨激光测高数据的精度进行了初步分析;文献[17]对大气折射引起的卫星激光测距延迟进行了分析,经模型修正后的大气折射误差小于2 cm,不是测距的主要误差源;文献[18-19]则重点围绕激光测高数据在林业森林参数提取应用中的全波形数据处理进行了研究,对激光几何定位处理方面则相对较少;文献[20-25]针对卫星激光测高数据与光学影像联合平差处理开展了研究,证实二者可以互补提高测图精度。

文章在第1节对资源三号02星的激光测高理论精度进行了分析,第2节对目前已有的02星激光测高数据进行了多个区域的精度验证,并将激光测高数据与立体影像进行了联合平差试验,最后在第3节进行了总结并提出了一些建议。

1 资源三号02星激光测高理论误差分析

资源三号02星激光测高仪的相关参数在笔者所著的文献[16]中已经做了较详细的说明,针对我国首台对地观测的试验性激光测高仪,结合各类误差源以及可行的改正精度,对其最终的测量误差进行定量分解非常必要,既有利于有针对性地开展数据处理工作,也能指导后续相关卫星的技术改进。国外在GLAS测量误差分解方面做了大量的工作,文献[26]分析了引起GLAS最终的测高误差(约0.15 m)的主要来源包括:器件本身的测距误差(0.10 m)、卫星轨道径向误差(0.05 m)以及激光指向角与地形坡度的综合影响(坡度/°且指向角误差为/″时误差为0.05 m)。在激光指向角接近天底方向时,对于大光斑激光足印点而言,因大气折射引起的平面位置偏差基本可以忽略,而大气折射引起的距离延迟最大可达2.35 m,必须加以改正。清洁大气条件下,采用NCEP(national center for environmental prediction)大气参数,目前大气折射延迟改正算法精度可控制在2 cm左右,不是激光测高的主要误差源。采用IERS(international earth rotation and reference systems service)公布的有关模型,潮汐影响的改正误差也可控制在1 cm。本文将资源三号02星试验性激光测高仪与GLAS的误差指标进行了对比,如表 1所示。在误差分解的过程中,假定了相应的系统误差经过了几何检校或部分误差可以通过数据后处理的方式得到消除。其中激光测距系统本身的系统误差、激光与星敏感器的安装角误差或激光指向角误差的系统误差需要通过在轨几何检校进行确定[13-15]。而约2.35 m的大气折射距离延迟以及40 cm的潮汐误差,经过事后大气折射延迟和潮汐模型修正后,仅残留2 cm和1 cm的随机误差。但云、雾霾等能见度不佳的天气条件下,大气对激光的影响不容忽视,相关的精细化改正算法还有待深入研究。

表 1 资源三号02星与ICESat/GLAS激光定位误差对比表 Tab. 1 Positioning error comparison between the ZY-3 02 and ICESat/GLAS laser altimeter
资源三号02星激光测高仪 ICESat/GLAS
随机误差 系统误差 产生的误差
激光器测距误差 68 cm 68 cm 10 cm
轨道径向误差 5 cm 5 cm 5 cm
安装偏心量误差 1 cm 1 cm 1 cm
安装偏置角误差 1″ 平面约2.5 m;
高程约4.1 cm
姿态/激光指向角误差 0.6″ 精密定姿后,无地面
控制5″
平面约14 m;
高程约23.7 cm
7.5 cm
大气折射 大气压10 mba,
可降水量5 mm
最大2.35 m,改正后
小于2 cm
3°以内,平面偏移
3 cm/°;高程:2 cm
2 cm
潮汐影响 1 cm 最大40 cm,改正后
小于1 cm
1 cm 1 cm
时间同步误差 20 us <1 cm <1 cm
地形起伏 假定坡度为1°且姿态/指向角误差为1″ 4.1 cm 5 cm
中误差 平面:14.2 m
高程:0.85 m(坡度=2°)
平面:4.5 m
高程:0.137 m

表 1可以看出,如果按资源三号02星给定的硬件测距误差1.0 m(3σ),取0.68 m(1σ)进行分析,在坡度为1°的平坦地区,理论上02星激光足印点的绝对高程精度为0.74 m,平面精度为14.2 m,若坡度为2°,则高程精度为0.85 m。其中激光器本身的测距误差是足印点高程的主要误差源,卫星姿态/激光指向角测量误差是平面误差的主要来源,而指向角与地形起伏的综合影响也不能忽视,尤其在有一定坡度区域,平面误差会引起较大的高程误差,这也是提取激光高程控制点时尽量选取光斑内地形平坦的原因所在[27]。虽然目前资源三号02星激光测高数据相关精度指标与GLAS还有一定距离,但参照国标的有关规定[28],平坦地区1.0 m的高程精度能够满足资源三号用于1:50 000立体测图的高程控制需求,而且该平坦地区特指激光足印光斑内的地形,即如果在山区或高山区,激光足印光斑正好位于一小块平坦区域,该激光点仍能作为山区的高程控制点使用。

2 资源三号02星激光测高精度验证 2.1 激光测高数据精度验证试验

为了保证资源三号02星激光测高数据处理的精度,消除硬件方面的系统性偏差,由国家测绘地理信息局卫星测绘应用中心牵头,联合黑龙江测绘地理信息局、陕西测绘地理信息局、内蒙古自治区测绘地理信息局、武汉大学、中科院安徽光机所、航天五院、北京市遥感信息研究所以及苏尼特右旗人民武装部等多家单位,精心组织了近100余名工程技术人员,于2016年8月中旬到9月上旬,在内蒙古自治区锡林郭勒盟苏尼特右旗相关区域进行了多次在轨几何检校试验,并取得圆满成功[13-14]。通过在轨几何检校,解算出了02星激光测高仪的指向角与测距系统性偏差。本文采用在轨几何检校后的参数、精密轨道参数、激光测距值,结合严密几何模型,并对大气折射延迟、潮汐等误差进行了改正,计算出激光足印点的三维坐标[16]

资源三号02星激光测高仪自2016年6月24日第一次开机,成功获得第一轨试验性激光测高数据以来,截至目前总共获得了44轨数据。笔者通过开发相应的处理软件,已经实现了02星所有激光测高数据的工程化处理,数据分布如图 1所示。由于是试验性载荷,仪器寿命有限,因此目前激光测高数据总量还不够多。02星激光测高仪每轨开机时间约为10 min,获得的点数在800~1000个左右,目前获得的总点数为35 314个。

图 1 资源三号02星44轨激光测高数据分布示意图 Fig. 1 Distribution of ZY-3 02 satellite 44 orbits laser altimetry data

本文针对02星获取的44轨激光测高数据,结合AW3D30(ALOS World 3D 30 m)数据进行了高程精度分析,对不同高差的点数所占比例进行了统计。由于AW3D30的高程精度为5 m[29],因此与其高差小于5 m的02星激光足印点才可能是有效点。总体而言,目前02星仅有的44轨激光数据,与参考AW3D30相比,高差小于2 m的点所占比例为10.37%,总数为3661个;小于5 m的为23.89%,总数8435个,统计结果如表 2所示。

表 2 资源三号02星44轨激光测高数据高程精度统计 Tab. 2 Elevation error statistical result of ZY-3 02 satellite 44 orbits laser altimetry data
序号 轨道号 时间 总数 无效 |dh|<2000 m |dh|<5 m |dh|<2 m |dh|<5 m/(%) |dh|<2 m/(%)
1 382 2016-06-24 472 14 319 268 144 56.78 30.51
2 656 2016-07-12 670 67 178 108 55 16.12 8.21
3 900 2016-07-28 911 43 38 15 4 1.65 0.44
4 913 2016-07-29 710 71 75 20 10 2.82 1.41
5 914 2016-07-29 706 44 450 308 148 43.63 20.96
6 915 2016-07-29 810 87 166 103 59 12.72 7.28
7 929 2016-07-03 706 48 402 256 150 36.26 21.25
8 930 2016-07-03 825 145 316 197 101 23.88 12.24
9 944 2016-07-31 862 47 569 320 102 37.12 11.83
10 945 2016-07-31 762 77 290 174 96 22.83 12.6
11 959 2016-08-01 619 31 376 174 69 28.11 11.15
12 960 2016-08-01 315 15 247 196 108 62.22 34.29
13 984 2016-08-02 813 335 421 289 115 35.55 14.15
14 990 2016-08-03 709 36 439 291 160 41.04 22.57
15 991 2016-08-03 906 71 541 326 153 35.98 16.89
16 1004 2016-08-04 458 68 165 66 31 14.41 6.77
17 1005 2016-08-04 825 50 421 304 135 36.85 16.36
18 1006 2016-08-04 885 127 522 354 133 40 15.03
19 1021 2016-08-05 819 104 458 247 128 30.16 15.63
20 1035 2016-08-06 970 89 317 141 53 14.54 5.46
21 1036 2016-08-06 527 76 237 165 88 31.31 16.69
22 1050 2016-08-07 589 12 82 57 27 9.68 4.58
23 1051 2016-08-07 746 76 434 250 142 33.51 19.03
24 1052 2016-08-07 946 75 259 190 58 20.08 6.13
25 1067 2016-08-08 926 85 149 89 38 9.61 4.1
26 1081 2016-08-09 890 69 331 171 69 19.21 7.75
27 1082 2016-08-09 918 121 344 189 91 20.59 9.91
28 1083 2016-08-09 905 3 842 689 272 76.13 30.06
29 1111 2016-08-11 1058 34 455 203 64 19.19 6.05
30 1112 2016-08-11 668 102 136 23 13 3.44 1.95
31 1143 2016-08-13 1048 118 413 239 92 22.81 8.78
32 1157 2016-08-14 1062 114 470 237 71 22.32 6.69
33 1172 2016-08-15 1070 72 463 243 95 22.71 8.88
34 1188 2016-08-16 699 67 128 61 27 8.73 3.86
35 1218 2016-08-18 825 101 138 43 12 5.21 1.45
36 1219 2016-08-18 933 112 435 183 57 19.61 6.11
37 1220 2016-08-18 784 13 694 293 75 37.37 9.57
38 1233 2016-08-19 900 120 364 119 38 13.22 4.22
39 1248 2016-08-02 951 76 324 137 61 14.41 6.41
40 1263 2016-08-21 760 58 231 145 70 19.08 9.21
41 1264 2016-08-21 716 91 167 72 43 10.06 6.01
42 1309 2016-08-24 826 99 321 134 55 16.22 6.66
43 1385 2016-08-29 811 85 490 239 84 29.47 10.36
44 1476 2016-09-04 1003 131 195 107 65 10.67 6.48
合计 35 314 4378 14 812 8435 3661
所占比例/(%) 12 41.90 23.89 10.37

为了验证02星激光测高精度,采用精检校后计算的激光足印点平面位置,在检校区附近平坦地区利用高精度RTK-GPS测量了8个激光足印点的绝对高程坐标,对比了检校前、基于地形匹配的粗检校以及基于红外探测器的精检校的高程误差[14],统计结果如表 3所示。显然,检校前资源三号02星激光测高数据存在较大的系统偏差,“两步法”在轨几何检校法[15]后能大幅提高绝对高程精度,该区域的激光高程验证精度为0.89 m,部分点高程误差小于0.5 m。

表 3 检校区附近部分激光足印点高程误差统计 Tab. 3 Elevation error statistical of laser altimetry footprint points near the field calibration region
点号 与RTK-GPS高程偏差/m
检校前 粗检校 精检校后
311 90.45 2.52 0.28
310 90.19 2.25 0.55
317 90.18 2.26 0.28
316 90.23 2.29 0.32
313 91.51 3.58 1.60
312 90.32 2.38 0.41
309 90.61 2.65 0.67
308 90.99 3.03 1.05
统计值 最大误差 91.51 3.58 1.60
最小误差 90.18 2.25 0.28
中误差 90.56 2.66 0.89
均值 90.56 2.62 0.65
标准偏差 0.471 0.47 0.47

由于大光斑激光的平面精度验证比较困难,因此采用的是在检校过程中被激光脉冲击中的地面探测器作为检查点。探测器的位置采用RTK-GPS可以精确测量获得,按式(1)计算得到激光足印点中心的真实位置,统计检校后经数据处理获得的激光足印点平面位置的偏差

(1)

式中,E(i, j)被击中的探测的能量值;(i, j)代表探测器所在格网的行列号;(ic, jc)为按探测器的位置计算出的激光足印中心点的相对位置,根据格网间隔以及RTK-GPS测量得到的探测器阵列的绝对位置,可内插出激光足印点的平面绝对位置坐标。

资源三号02星共进行了5次检校试验,但只有3次成功[13],共有3个激光足印点可用,分别为编号A6809、E2818、W1909,其平面位置误差统计如表 4所示。经检校和数据处理后平面精度优于15.0 m的水平,本文采用式(1)后结果比文献[14]有一定进步。由于平面精度检验的样本有限,其精度的可靠性以及稳定性还有待深入地分析验证,但对于大光斑激光足印点作为高程控制点应用而言,目前15.0 m的水平精度在足印点光斑为平坦地形时可以满足应用需求[27]

表 4 检校前后被探测到的激光足印点平面位置误差统计 Tab. 4 Planimetry error statical result of the captured laser footprint points
点号 平面位置偏差/m
直接定位 粗检校 精检校后
A6809 8 066.29 95.46 12.46
E2818 8 047.11 86.34 16.61
W1909 8 065.10 84.87 14.91
中误差 8 059.505 89.013 14.76

为进一步验证02星激光点的绝对高程精度,选取了华北某地高精度的DSM数据作为参考,该DSM为1:2000基础地理信息成果图,格网大小为2.0 m,平面精度1.0 m,高程精度平地0.4 m,丘陵地区0.7 m,平面坐标系为WGS 84,高程基准为WGS 84大地高。2016年9月4日资源三号02星第1476轨正好经过该区域,激光测高仪开机获得了部分有效激光测高数据,如图 2所示。

图 2 华北某地高精度DSM及资源三号02星激光足印点示意图 Fig. 2 Distribution of ZY-3 02 satellite laser footprint points and referenced DSM in north China

该区域离内蒙古的外业检校区有几百千米的距离,且激光获取时间与检校时间也有一定间隔,在时间和空间上有一定代表性,能反映出检校后的激光测高精度水平。区域内共有15个02星激光足印点落在参考DSM的范围内,其中有2个属于硬件记录无效点,对余下13个点采用AW3D30以5 m作为阈值进行有效点筛选,最终保留8个点,利用1:2000高精度DSM参考地形数据进行了精度评价。由于1:2000 DSM格网大小为2.0 m,而02星激光为大光斑足印,因此高程精度统计时采用了两种统计方式:①根据激光足印点的平面坐标直接从参考DSM中内插出高程;②根据激光足印点平面位置及光斑大小统计光斑范围内所有的DSM高程值的平均值,其中根据外场检校结论,足印光斑大小取120 m。最终统计结果如表 5所示。从表 5中可以看出,经处理后的02星绝对高程精度可以达到1.09 m的水平,少部分激光点的高程误差小于0.5 m。

表 5 ZY-3 02星第1476轨部分激光足印点高程统计表 Tab. 5 Elevation error statistical result of ZY-3 02 satellite laser footprint points on orbit 1476
点编号 经度/(°E) 纬度/(°N) 大地高/m 与AW3D30
高差/m
与1:2000 DSM高差/m
直接内插 取平均
494 117.558 690 39.355 50 -5.432 -1.99 -1.60 -1.19
500 117.504 948 39.167 580 -3.550 -1.68 -0.21 -0.08
501 117.496 017 39.136 256 -3.997 -0.98 0.31 0.53
502 117.487 096 39.104 932 -7.737 -4.45 -1.86 -1.78
503 117.478 156 39.073 607 -4.430 -5.38 1.33 1.04
504 117.469 238 39.042 278 -4.866 -3.65 -0.24 0.33
505 117.460 337 39.010 948 -6.947 -2.35 -1.19 -1.18
508 117.433 636 38.916 982 -3.268 -1.96 1.15 1.45
统计结果 均值 -2.81 -0.29 -0.11
中误差 3.15 1.16 1.09
最大误差 -5.38 -1.86 -1.78
最小误差 -0.98 -0.21 -0.08

对第1035轨激光测高数据,选取海面上中间的一段激光足印点,如图 3所示,采用EGM2008(earth gravitational model 2008)地球重力场模型数据[30],计算激光足印点大地水准面高,其高程剖面如图 4所示。经统计,该轨数据在海面上的激光足印点高程均值为-0.23 m,中误差为0.47 m,说明其内部相对精度比较高,而最大与最小高程相差1.32 m可能是海浪的影响,同时还包括EGM2008模型的本身误差。

图 3 第1035轨激光足印点海面足印分布图 Fig. 3 Distribution of orbit 1035 laser footprint points on the sea

图 4 第1035轨海面上激光足印点高程剖面图 Fig. 4 The profile of orbit 1035 laser footprint points on the sea

2.2 激光测高数据与立体影像联合平差试验

针对资源三号02星上同时搭载的立体相机和试验性激光测高仪,选取位于西部渭南地区作为试验区开展联合平差试验。具体的平差算法,笔者在文献[20-21]中均进行了阐述,试验中采用基于有理函数模型RFM(rational function model)加像方平移的方式对RPC(rational polynomial coefficients)参数进行补偿,如式(2)和(3)所示

(2)
(3)

式中,(φn, λn, Hn)为归一化的物方三维坐标;(cn, rn)为归一化的像方坐标;(c0, r0)为像方坐标的平移参数;(φ0, λ0, H0)为物方三维坐标的平移参数;(cs, rs)为像方坐标的缩放系数;(φs, λs, Hs)为物方三维坐标的缩放系数;p1ijkp2ijkp3ijkp4ijk(i=1, 2, 3;j=1, 2, 3;k=1, 2, 3)代表标准的RPC参数;(Δc, Δr)表示像方补偿数值,本文中主要采用平移变换模型。

选取的试验区的地理范围为:[33.97°N,35.31°N],[108.74°E,110.07°E],区域中间是平原和丘陵,其他部分是山区,整个区域高差在1000 m左右。资源三号02星8月9日的第1081轨和8月14日的第1157轨经过该区域,同步获得了激光和影像数据,分别有36个和38个硬件记录为有效的激光点落在试验区。试验中选取了27个外业高精度GPS点作为检查点,综合考虑影像上的刺点误差,检查点的绝对平面精度优于0.5 m,高程精度优于0.2 m。在与参考的AW3D30高差小于2.0 m的11个激光足印点中,选择了6个点作为高程控制点,剩下的5个由于点位难以确定或地形存在明显的不确定性而被舍弃,激光高程控制点与检查点分布如图 5(a)所示。考虑到AW3D30本身存在一定误差,而保留下的高差小于2.0 m的激光足印点太少,进一步放宽高差限制阈值,将与AW3D30相比高差小于5.0 m的激光足印点也进行了保留,从中选取了16个点作为高程控制点,此时的激光高程控制点与检查点分布如图 5(b)所示。

图 5 激光高程控制点与检查点分布示意图 Fig. 5 Distribution of ZY-3 02 laser elevation control points and check points

在实际平差试验中,共进行了如下几种处理,最终统计结果如表 6所示。

表 6 资源三号02星激光与影像联合处理精度统计表 Tab. 6 Statistical result of ZY-3 02 laser and stereo images combined processing
平差类型 控制点 检查点 中误差/m
东方向 北方向 平面 高程
无地面控制 0 27 22.424 12.350 25.599 11.542
27平高 27 0 3.353 2.949 4.466 1.629
5平高 5 22 3.857 3.152 4.981 1.930
5平0高 5 22 3.328 2.547 4.192 10.135
0平6高 6 27 19.787 11.672 22.973 2.936
0平16高 16 27 19.050 11.537 22.271 1.902

(1) 完全无地面控制,统计27个检查点的误差;

(2) 将27个点全部作为平高控制点,没有检查点,统计最终控制点的残差;

(3) 选择四角和中间共5个点作为平高控制点,统计余下的22个检查点误差;

(4) 选择四角和中间共5个点作为平面控制点,统计余下的22个检查点误差;

(5) 选择与AW3D30高差小于2.0 m的6个激光足印点作为高程控制点,统计27个检查点误差;

(6) 选择与AW3D30高差小于5.0 m的16个激光足印点作为高程控制点,统计27个检查点误差。

表 6可以看出,在有激光高程控制点参与的情况下,影像无地面控制的高程精度能有较大提升。基于有理函数模型,在16个激光高程控制点下,本试验区资源三号02星影像无地面控制点高程测量精度从11.542 m提高到1.902 m,与采用5个全野外控制点的结果相当,提升效果非常明显。

3 总结与展望

本文针对资源三号02星试验性激光测高载荷,开展了精度分析与多个区域的验证,同时对国产激光测高数据提高卫星影像无地面控制高程精度的应用开展了初步试验,能得到如下几个结论:

(1) 作为我国首台对地观测的卫星激光测高试验性载荷,资源三号02星激光测高仪能获得部分有效数据,以AW3D30高程作为参考,02星激光有效数据约占23.89%。

(2) 在平坦地区(坡度≤2°),资源三号02星激光测高仪的理论高程精度为0.85 m、平面精度14.2 m。激光器硬件本身测距精度以及高精度姿态测量精度的不足,是制约国产卫星激光测高几何精度的主要因素。华北地区高精度DSM地形数据验证02星激光足印点高程精度为1.09 m;内陆渤海海面上的激光高程精度为0.47 m;在检校场区域高程验证精度为0.89 m,平面精度为14.76 m。

(3) 将激光足印点作为高程控制点时,在陕西渭南试验区能将资源三号02星立体影像无地面控制的高程精度从11.54 m提高到1.90 m,影像的高程精度提升效果非常明显。

资源三号02星激光测高仪为试验性载荷,暂时无法业务化运行工作,但卫星激光测高数据在全球测图工程中具有重要的应用价值,国家测绘地理信息局在十三五规划中已明确提出“全球地理信息资源开发”是五大核心任务之一,因此建议在资源三号后续卫星中能配备业务化应用的激光测高仪。此外,为了在全球气候变暖、森林碳储量监测、极地冰盖测量以及境外高精度地形测绘等方面掌握主动权,我国需要大力发展专门的激光测高卫星,要瞄准下一代的单光子激光测高卫星,提前布局、积极开展相关的指标论证与关键技术攻关,努力构建我国自主的激光测高卫星体系。

致谢: 特别感谢黑龙江测绘地理信息局、陕西测绘地理信息局、内蒙古自治区测绘地理信息局、武汉大学、中科院安徽光机所、航天五院、北京市遥感信息研究所以及苏尼特右旗人民武装部等资源三号02星激光测高仪外业检校参研单位,感谢高小明、付兴科、谢俊峰、陈继溢、窦显辉、张悦、朱广彬等对论文试验数据提供的帮助。


参考文献
[1] SMITH D E, ZUBER M T, NEUMANN G A, et al. Topography of the Moon from the Clementine Lidar[J]. Journal of Geophysical Research Atmospheres, 1997, 102(E1): 1591–1611. DOI:10.1029/96JE02940
[2] CAVANAUGH J F, SMITH J C, SUN Xiaoli, et al. The Mercury Laser Altimeter Instrument for the MESSENGER Mission[J]. Space Science Reviews, 2007, 131(1-4): 451–479. DOI:10.1007/s11214-007-9273-4
[3] SMITH D E, ZUBER M T, FREY H V, et al. Mars Orbiter Laser Altimeter:Experiment Summary after the First Year of Global Mapping of Mars[J]. Journal of Geophysical Research:Planets, 2001, 106(E10): 23689–23722. DOI:10.1029/2000JE001364
[4] SMITH D E, ZUBER M T, JACKSON G B, et al. The Lunar Orbiter Laser Altimeter Investigation on the Lunar Reconnaissance Orbiter Mission[J]. Space Science Reviews, 2010, 150(1-4): 209–241. DOI:10.1007/s11214-009-9512-y
[5] GARVIN J, BUFTON J, BLAIR J, et al. Observations of the Earth's Topography from the Shuttle Laser Altimeter(SLA):Laser-pulse Echo-Recovery Measurements of Terrestrial Surfaces[J]. Physics and Chemistry of the Earth, 1998, 23(9-10): 1053–1068. DOI:10.1016/S0079-1946(98)00145-1
[6] SCHUTZ B E, ZWALLY H J, SHUMAN C A, et al. Overview of the ICESat Mission[J]. Geophysical Research Letters, 2005, 32(21): 213–219.
[7] WANG Xianwei, CHENG Xiao, GONG Peng, et al. Earth Science Applications of ICESat/GLAS:A Review[J]. International Journal of Remote Sensing, 2011, 32(23): 8837–8864. DOI:10.1080/01431161.2010.547533
[8] ABDALATI W, ZWALLY H J, BINDSCHADLER R, et al. The ICESat-2 Laser Altimetry Mission[J]. Proceedings of the IEEE, 2010, 98(5): 735–751. DOI:10.1109/JPROC.2009.2034765
[9] YU A W, KRAINAK M A, STEPHEN M A, et al.Spaceflight Laser Development for Future Remote Sensing Applications[C]//Proceedings of SPIE 8182, Lidar Technologies, Techniques, and Measurements for Atmospheric Remote Sensing:VⅡ.Prague, Czech Republic:SPIE, 2011:818-204.
[10] COYLE D B, STYSLEY P R, POULIOS D, et al.Laser Transmitter Development for NASA's Global Ecosystem Dynamics Investigation (GEDI) Lidar[C]//Proceedings of LiDAR Remote Sensing for Environmental Monitoring:XV.San Diego, California:SPIE, 2015:1021-1028.
[11] 崔成玲, 李国元, 黄朝围, 等. 国内外星载激光测高系统发展现状及趋势[C]//中国测绘地理信息年会优秀青年论文集. 南昌: [s. n. ], 2015: 288-295.
CUI Chengling, LI Guoyuan, HUANG Chaowei, et al.The Status and Trend of Satellite Laser Altimeter at Home and Abroad[C]//Proceedings of the Best Papers of the Annual of Chinese Surveying Mapping and Geo-information.Nanchang:[s.n.], 2015:288-295.
[12] 李鑫, 廖鹤, 赵美玲, 等. 激光测绘卫星对不同地表形貌探测能力分析[J]. 测绘学报, 2014, 43(12): 1238–1244.
LI Xin, LIAO He, ZHAO Meiling, et al. Research on LiDAR Surveying Satellite Detection Capacity for Different Terrains[J]. Acta Geodaetica et Cartographica Sinica, 2014, 43(12): 1238–1244. DOI:10.13485/j.cnki.11-2089.2014.0188
[13] 唐新明, 谢俊峰, 付兴科, 等. 资源三号02星激光测高仪在轨几何检校与试验验证[J]. 测绘学报, 2017, 46(6): 714–723.
TANG Xinming, XIE Junfeng, FU Xingke, et al. ZY3-02 Laser Altimeter On-orbit Geometrical Calibration and Test[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(6): 714–723. DOI:10.11947/j.AGCS.2017.20160597
[14] 李国元, 唐新明, 陈继溢, 等. 基于地面红外探测器的星载激光测高仪在轨几何定标[C]//第三届新型光电探测技术及其应用研讨会会议论文集. 西安: [s. n. ], 2016.
LI Guoyuan, TANG Xinming, CHEN Jiyi, et al.On-orbit Geometric Calibration of Satellite Laser Altimeter Using Ground-based IR Detectors[C]//Proceedings of the Third Seminar on Novel Dectection Technology and Application.Xi'an:[s.n.], 2016.
[15] 李国元. 对地观测卫星激光测高数据处理方法与工程实践[D]. 武汉: 武汉大学, 2017.
LI Guoyuan.Earth Observing Satellite Laser Altimeter Data Processing Method and Engineer Application[D].Wuhan:Wuhan University, 2017.
[16] 唐新明, 李国元, 高小明, 等. 卫星激光测高严密几何模型构建及精度初步验证[J]. 测绘学报, 2016, 45(10): 1182–1191.
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. DOI:10.11947/j.AGCS.2016.20150357
[17] 李松, 马跃, 周辉, 等. 星载激光测高系统大气折射延迟的影响分析[J]. 应用光学, 2014, 35(3): 520–524.
LI Song, MA Yue, ZHOU Hui, et al. Analysis on Atmospheric Refraction Delay for Satellite Laser Altimeter System[J]. Journal of Applied Optics, 2014, 35(3): 520–524.
[18] 李增元, 刘清旺, 庞勇. 激光雷达森林参数反演研究进展[J]. 遥感学报, 2016, 20(5): 1138–1150.
LI Zengyuan, LIU Qingwang, PANG Yong. Review on Forest Parameters Inversion Using LiDAR[J]. Journal of Remote Sensing, 2016, 20(5): 1138–1150.
[19] 王成, 习晓环, 骆社周, 等. 星载激光雷达数据处理与应用[M]. 北京: 科学出版社, 2015.
WANG Cheng, XI Xiaohuan, LUO Shezhou, et al. Space-borne LiDAR Data Processing and Its Application[M]. Beijing: Science Press, 2015.
[20] 李国元, 唐新明, 王华斌, 等. GLAS激光测高数据辅助的资源三号三线阵区域网平差研究[C]//第三届高分辨率对地观测学术年会优秀论文集. 北京: [s. n. ], 2015.
LI Guoyuan, TANG Xinming, WANG Huabin, et al.Research on the ZY-3 Block Adjustment Supported by the GLAS Laser Altimetry Data[C]//Proceedings of Best Papers of the 3rd Annual Conference on High Resolution Earth Observation.Beijing:[s.n.], 2015.
[21] LI Guoyuan, TANG Xinming, GAO Xiaoming, et al. ZY-3 Block Adjustment Supported by GLAS Laser Altimetry Data[J]. The Photogrammetric Record, 2016, 31(153): 88–107. DOI:10.1111/phor.2016.31.issue-153
[22] LI Guoyuan, TANG Xinming, GAO Xiaoming, et al.Improve the ZY-3 Height Accuracy Using ICESat/GLAS Laser Altimeter Data[C]//International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Volume XLI-B1.Prague, Czech Republic:ISPRS, 2016:37-42.
[23] 胡文敏, 邸凯昌, 岳宗玉, 等. 嫦娥一号激光高度计数据交叉点分析与平差处理[J]. 测绘学报, 2013, 42(2): 218–224.
HU Wenmin, DI Kaichang, YUE Zongyu, et al. Crossover Analysis and Adjustment for Chang'E-1 Laser Altimeter Data[J]. Acta Geodaetica et Cartographica Sinica, 2013, 42(2): 218–224.
[24] 赵双明, 冉晓雅, 付建红, 等. CE-1立体相机与激光高度计数据联合平差[J]. 测绘学报, 2014, 43(12): 1224–1229.
ZHAO Shuangming, RAN Xiaoya, FU Jianhong, et al. Combined Adjustment of CE-1 Stereo Camera Image and Laser Altimeter Data[J]. Acta Geodaetica et Cartographica Sinica, 2014, 43(12): 1224–1229. DOI:10.13485/j.cnki.11-2089.2014.0178
[25] 王任享, 王建荣. 二线阵CCD卫星影像联合激光测距数据光束法平差技术[J]. 测绘科学技术学报, 2014, 31(1): 1–4.
WANG Renxiang, WANG Jianrong. Technology of Bundle Adjustment Using Two-line-array CCD Satellite Image Combined Laser Ranging Data[J]. Journal of Geomatics Science and Technology, 2014, 31(1): 1–4.
[26] SCHUTZ B E, URBAN T J.The GLAS Algorithm Theoretical Basis Document for Laser Footprint Location(Geolocation) and Surface Profiles[R].[s.l.]:NASA, 2014.
[27] 李国元, 唐新明, 张重阳, 等. 多准则约束的ICESat/GLAS高程控制点筛选[J]. 遥感学报, 2017, 21(1): 96–104.
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.
[28] 国家测绘局. GB/T 12341-2008 1: 25 000 1: 50 000 1: 100 000地形图航空摄影测量外业规范[S]. 北京: 中国标准出版社, 2008.
State Bureau of Surveying and Mapping.GB/T 12341-2008 Specifications for Aerophotogrammetric Field Work of 1:25 000 1:50 000 1:100 000 Topographic Maps[S].Beijing:China Standard Press, 2008.
[29] TAKAKU J, TADONO T, TSUTSUI K.Generation of High Resolution Global DSM from ALOS PRISM[C]//Proceedings of International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences:XL-4.Suzhou, China:ISPRS, 2014:243-248.
[30] PAVLIS N K, HOLMES S A, KENYON S C, et al.An Earth Gravitational Model to Degree 2160:EGM 2008[C]//Proceedings of the 2008 General Assembly of the European Geosciences Union.Vienna, Austria:[s.n.], 2008.
http://dx.doi.org/10.11947/j.AGCS.2017.20170174
中国科学技术协会主管、中国测绘地理信息学会主办。
0

文章信息

李国元,唐新明
LI Guoyuan, TANG Xinming
资源三号02星激光测高精度分析与验证
Analysis and Validation of ZY-3 02 Satellite Laser Altimetry Data
测绘学报,2017,46(12):1939-1949
Acta Geodaetica et Cartographica Sinica, 2017, 46(12): 1939-1949
http://dx.doi.org/10.11947/j.AGCS.2017.20170174

文章历史

收稿日期:2017-04-11
修回日期:2017-08-21

相关文章

工作空间