卫星激光测高严密几何模型研究.docx

《卫星激光测高严密几何模型研究.docx》由会员分享，可在线阅读，更多相关《卫星激光测高严密几何模型研究.docx（13页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、卫星激光测高严密几何模型研究采用星载激光测高仪辅助提高卫星立体影像几何定位精度十分是高程精度，已经得到了航天摄影测量界的重视，计划于2018年发射的高分七号卫星上将同时搭载光学立体相机和激光测高仪。固然，已有相关文献针对美国的ICESat卫星上搭载的地球科学激光测高系统GLAS的几何模型和产品精度做了相关介绍，但对其严密的几何定位模型和精度验证目前还没有系统性的阐述。本文较全面地对激光测高卫星的严密几何模型进行了构建与精度分析，并选择ICESat/GLAS的0级辅助文件，采用严密几何模型重现了2级产品的生产经过。将本文计算的结果与ICESat/GLAS的结果进行了比照分析，其中基于几何模型的高

2、程误差约11cm，平面误差在3cm以内，表明所提出的严密几何模型的正确性，同时采用新发射的资源三号02星的激光测高数据进行了初步处理和验证；相关结论可为国产高分后续卫星的激光测高数据处理提供参考。关键字：星载激光测高；严密几何模型；数据处理；精度验证1.引言激光雷达与成像光谱技术、合成孔径雷达技术一起被列为对地观测系统最核心的信息获取技术。美国NASA在1994和1997年两次将航天激光测高仪SLA安装在航天飞机上，以10Hz的频率、间隔750m和100m的脚印光斑采样方式获取并建立全球控制点库1。2003年，美国还发射了目前为止唯一的一颗对地观测激光测高卫星ICESat，该卫星上搭载了地球科

3、学激光测高系统GLAS，在极地冰盖监测、全球森林生物量估算、陆地高程测量等方面得到广泛应用2,3，而后续的ICESat-2和LIST卫星也正在研发之中4,5。计划中的国产高分七号以及陆地生态系统碳监测卫星上均搭载了激光测高仪，在2016年5月30日发射的资源三号02星上搭载了一个激光试验性载荷，唐新明7等建立了国产资源三号高分辨率光学测绘卫星的严密成像模型，而针对国产的激光测高卫星的定位模型研究基本处于空白。因而，对激光测高卫星的严密几何定位模型进行系统性的研究非常必要。在激光测高数据应用方面，朱长明7等从ICEsat/GLAS激光测高数据中反演出湖泊的动态水位高程，但前提是所选的激光测高足印

4、点要精度可靠；Li8等采用ICESat/GALS数据提出了一种基于机器学习的实时地形分类算法，并采用北京区域的数据开展了试验，且对数据的精度要求较高；胡文敏9等对我国的嫦娥一号激光测高数据开展了轨道穿插点分析和平差方法研究，Wu10、赵双明11等开展了嫦娥探月的激光测高数据与光学影像联合平差研究，但与对地观测还有一定差异；王任享12，李国元13等结合对地观测的卫星激光测高仪辅助光学立体影像平差和立体测图开展了相关研究，证实了激光测高数据辅助提高光学影像测图精度的可行性；李建成14、黄海兰15等曾采用ICESat/GLAS数据开展南极冰盖高程模型研究；Wang2等较系统地对ICESat/GLAS

5、在地学方面的应用作了梳理。而较高的数据精度是应用的前提。在精度分析方面，李鑫16等对激光测高卫星在不同地形条件下的探测能力进行了分析，但没有给出定量的结论；范春波17、朱剑锋18等对ICESat/GLAS脚点定位模型进行了介绍，但对于卫星运行速度引起的光行差改正没有提及，而在600km的轨道高度，激光单向传输约2ms的时间，光行差可达5角秒，对应地面距离约15m；马跃19等对星载激光测高系统高程误差进行了分析，但采用的是简化的定位模型，对硬件安装误差、测距误差等没有深化涉及；GLAS的ATBD20报告对GLAS脚点定位流程及误差做了较为具体的介绍，但该文档对一些细节并没有涉及。文章推导了激光测

6、高卫星严密几何定位模型，并对光行差、硬件安装误差做了分析，具体介绍了ICESat/GLAS的激光脚点定位原理和数据处理流程，利用GLAS数据重现了从0级到2级的生产经过，将计算结果与GLAS相应产品做了比照分析，同时利用资源三号02星的首轨激光测高数据的初步处理结果进行了验证。2.严密几何定位模型构建2.1严密几何定位模型激光测高卫星的基本原理是：通过卫星发射激光束经地面反射后由卫星接收，计算激光发射和接收的时间间隔t，光的传播速度为c，则激光单程传输距离p=c*t/2，结合卫星上搭载的GPS和星敏感器获得的卫星位置和姿态信息，即可获得激光脚印点的三维坐标。其严密几何模型如图1所示，其中las

7、erP为激光发射的参考点，GPSP为GPS天线相位中心，BodyO为卫星质心，GroundP为激光地面足印点。为描绘星载激光测高严密几何模型，定义卫星本体坐标系：卫星质心为原点，X轴指向卫星飞行方向，Z轴指向天顶方向，Y轴垂直于卫星轨道平面，与X、Z轴构成右手坐标系。激光发射时指向方向与卫星本体坐标系存在一定的夹角，假设激光指向与本体坐标系Z轴的负向夹角为，在XOY平面上的投影与X轴正向夹角为，如图2所示。若激光的测距值，激光参考点在本体坐标系中的坐标为：TLLxLyLz，则在激光脚点在卫星本体坐标系下的坐标为：sincossinsincosBodyLxLyL1卫星搭载的GPS定位系统能够测定

8、卫星的位置TGPSGPSGPSXYZ，但GPS的相位中心一般与卫星质心不完全重合，两者之间存在一定的偏差，假定GPS的相位中心在卫星本体坐标系下的坐标为：TPDxDyDz。卫星星敏感器相对于卫星本体坐标的安置矩阵为StarBodyR，而星敏感器测量的是星敏本体系到J2000坐标系下的旋转矩阵J2000StarR，所以还需要J2000坐标系到ITRF的坐标变换矩阵2000ITRFJR。2.2光行差引起的位置误差激光测高作为主动遥感方式，其传输途径为卫星到地面距离的两倍，假如卫星高度为600km，则往返一次所需时间约4ms，卫星前进距离约30m。以图3为例，卫星在A点发射一束激光，激光出射方向为A

9、G，卫星运行速度为v，光速为c，则当激光束在天球坐标系下的实际运动方向为AG，激光束到达地面的实际位置为G，此时卫星位置为B点，接收到激光回波时卫星位置为C。假如按发射时刻A点时的卫星姿态和位置计算地面点的位置，则与实际存在位置偏差GG。激光发射方向和实际传输方向的夹角为，由光速和卫星运行速度可得22：sinvc3其中为卫星运行速度与激光出射方向的夹角。为了尽量减少因激光高度角而引起的大气折射延迟，激光出射方向一般为天底点方向，即90。因而，5.2角秒，GG15m。此时，考虑因光行差对高程测量误差的影响，如图4所示，当激光指向角1时，此时因光行差而引起的高程误差为0.212m。当卫星轨道高度降

10、为500km时，因光行差而引起的高程误差虽略有下降，但随着指向角的变大会线性变大。为有效消除光行差引起的位置误差，ICESat卫星采用激光发射时刻的卫星姿态、激光到达地面时卫星的位置，即图3中卫星在A点的姿态和B的位置，解算激光足印在地面的实际位置G。2.3硬件安装误差引起的位置误差由公式2可知，硬件安装误差主要包括激光指向与卫星本体系的两个角度测量误差，卫星本体系与星敏本体系的轴向安装误差，以及激光参考点与卫星本体系原点、GPS天线相位中心与卫星本体系原点的偏心量测量误差。偏心量引起的位置误差基本是线性关系，而且目前偏心量的量测误差基本可达mm级。综合考虑，激光指向与卫星本体系以及卫星本体系

11、与星敏本体系的夹角关系，可简化为激光指向与星敏本体系存在两个夹角。3.ICESat/GLAS处理流程介绍ICESat/GLAS于2003年1月发射成功，2009年10月停止工作，分周期共工作了2465天，采集激光点约20亿个。GLAS数据产品分为0级、1级和2级，其中1级包括1A和1B，1级和2级产品又分为15类，即公开的GLA01GLA1523,24。ICESat/GLAS项目组公布了公式2所建立的严密几何定位模型的简化计算方法20，即采用ANC04惯性坐标系旋转矩阵、ANC08精细定轨数据(POD)25、ANC09精细定姿数据(PAD)22、ANC25时间转换数据以及激光测距值计算脚点三维

12、坐标，其流程如图7所示。由发射和接收波形计算激光传输距离的原理如图8所示，图中不规则的曲线表示原始波形，规则曲线表示经最小二乘拟合后的波形，从拟合的波形数据中能够提取出发射与接收波形的时间重心，根据时间重心的间隔与光速即可计算出激光往返的距离26。由于GLAS激光发射波形可近似为高斯脉冲，激光经地外表反射后回波波形可近似看做一次或屡次高斯脉冲的叠加。对于多回波波形数据，采用多个高斯函数叠加进行拟合26。结合全波形数据处理、大气延迟改正以及固体潮改正等内容，卫星激光测高数据几何处理流程可描绘为：1对发射和接收波形进行平滑、去噪，提取波形特征参数，确定发射和接收波形重心对应的时刻Ts、Tr；激光传

13、输时间间隔为:tTrTs，激光到达地面时刻2tTgTs，单向传输距离为2tc，sys，其中sys是硬件引起的激光测距系统误差改正值；2根据激光发射时刻Ts和PAD数据计算在ICRF坐标系下激光指向的单位向量u；3根据激光到达地面的时刻Tg和POD数据计算激光参考点在ICRF坐标系下的坐标向量gr；4根据激光到达地面的时刻Tg，计算ICRF到ITRF的旋转矩阵ITRFICRFR；5利用简化模型()ITRFsICRFgrRru计算激光脚印点在ITRF下的坐标向量sr；6利用sr对应的经纬度、高程以及NCEP(NationalCenterforEnvironmentalPrediction)大气参数

14、和大气延迟改正模型计算激光大气延迟改正值atm；7对测距值进一步进行修正为newatm，同时重复步骤5计算snewr；8对新计算的激光脚印点在ITRF下坐标向量snewr采用潮汐改正模型计算改正量tr，最终准确的激光脚印点在ITRF下的坐标向量为:snewtrr。在上述步骤中，步骤5所对应的为1B级产品GLA05，步骤7所对应的为1B级产品GLA06，步骤8所对应的为2级产品GLA1215。4.精度验证试验本文选取了ICESat/GLAS某一天的辅助数据文件ANC04、ANC08、ANC09、ANC25以及1A级产品GLAS01中的波形参数，采用严密几何模型求解每个激光足印光斑的三维坐标，并与

15、2级产品GLA14中的结果进行比照，其中辅助数据文件及GLAS的1A和2级数据产品从NSIDC网站下载而来。同时对资源三号02星激光试验性载荷的第一轨数据进行了处理与精度验证。4.1基于几何模型的精度比照文章选取了2009年3月24日的ICESat/GLAS激光测高数据，其陆地区域的覆盖情况如图9所示，共约15圈22308个点。采用文章第三部分介绍的数据处理流程，将计算结果与GLA14数据进行比照，统计结果见表1。在计算经过中，为了与GLA14结果中的Topex/Poseidon椭球保持一致，在ITRF坐标系转成经纬度以及最终投影时均采用T/P椭球参数，其参数为：长轴a=6378136.30m

16、，扁率1/f=298.257。单位：cm从表1中能够看出，沿轨和垂轨方向的误差与GLAS的结果相比均到达cm级，其中垂轨方向的多个轨迹数据误差统计的平均值接近于0，这与屡次观测的理论平均误差结果是一致的。此外，垂轨和沿轨两者的误差最大值在57cm，从误差数量级的角度来讲，证实整个基于模型的计算经过是正确的。相比GLAS标称的15cm的高程测量误差，本文中基于模型计算的高程相对较差，固然平均值较小约1.7cm，但最大值超15.0cm。经与GLAS的中间结果比照检查，主要是大气延迟改正和波形拟合时间重心鉴别上存在误差，其中本文计算的大气延迟改正与GLAS给定的改正值差值中误差为1.35cm，但差值

17、最大值为7.64cm。在平坦地形如图10本文提取的波形时间重心与GLAS的基本一致，误差在0.2ns，但在森林或复杂地形区域如图11，则与GLAS的时间重心相差较大，误差达1.5ns。4.2资源三号02星激光试验载荷初步处理精度验证资源三号02星于2016年5月30日成功发射，卫星除载有相比01星分辨率略有提升的三线阵相机和多光谱相机外，还额外搭载了试验性的激光测高载荷，该激光测高载荷的有关参数如表2所示。2016年6月24日，资源三号02星激光试验性载荷开机236秒，获得了第一轨共472个激光数据，其中14个为无效点。对激光指向角进行初步优化后，按本文的处理方法，计算出每个点的三维大地坐标经

18、度、纬度和大地高，其地理位置跨蒙古国以及我国的内蒙古、甘肃和青海地区,如图12所示，高程剖面如图13所示。同时收集了该区域由资源三号和高精度控制数据生产的精度优于3.0m、格网大小为15m的DSM数据作为参考6，评价激光点的绝对高程精度。初步分析发现部分激光点因云、霾及地形等原因导致高程值较大，如图13中最右侧的一些数据点，其高程值到达20000m，存在明显误差，基本能确定是受大气中影响，造成测距值偏小，引起高程偏大。结合实验区的DSM数据的高程范围，实验中以2500米为阈值，将解算出来的激光点高程值大于2500m的去掉不介入计算，对剩下的177个点，利用激光点的经纬度从DSM中内插的高程作为

19、参考值，对激光点的高程误差进行统计分析，高程误差分布直方图见图14，其统计结果如表3所示。5.分析与结论文章第二部分对星载激光测高卫星的严密几何定位模型进行系统构建，同时对光行差引起的平面和高程误差进行了分析。在600km轨道高度，因激光光行差而引起的平面位置误差约15m，在定位时必须考虑。而由图4可知，为了保证激光测高的高程精度，激光指向角必须接近0，若为了使高程误差优于15cm，则激光指向角应不大于1，这在图5和6中也得到进一步的验证。在让激光指向角尽量小的同时，指向角的测量误差d应优于1.5，这与ICESat/GLAS的要求是一致的，而另一个指向角及其测量误差要求则能够相对降低，但为了控

20、制平面误差，角也应该尽量小，将激光指向尽量控制在卫星本体坐标系的OXZ平面范围内且接近Z轴的负向。通过采用目前唯一的对地观测激光测高卫星ICESat/GLAS的实际数据对严密几何定位模型的精度进行验证，与GLAS的二级产品GLA14的结果比照，其精度在平面方向与GLAS结果最大偏差在57cm，且多轨数据在垂轨方向的均值接近于0。高程方向的偏差中误差为11.3cm，均值约1.7cm，但最大偏差超过19cm，主要由于波形时间重心拟合值与GLAS的值存在一定偏差，且部分区域大气改正值大于GLAS的结果。对资源三号02星第一轨激光数据进行初步处理，同时采用地面DSM数据进行了验证，初步结果表明激光点的

21、高程中误差为5.73m，离激光测距标称的精度还有一定距离。其原因，一方面是由于DSM内插的高程与实际高程有一定的偏差；另一方面是由于激光器的各项硬件参数十分是激光指向角未经过准确在轨检校，与实验室测定的结果有一定的差距，需进一步开展准确标定工作。总体实验结果表明，文中第2部分的模型是正确的，只是在复杂地形及云霾等大气条件下还需进一步的改良和讨论。6.总结与瞻望文章对激光测高卫星严密几何定位模型进行了较为系统地描绘，同时对其他文献没有涉及的光行差以及硬件安装误差引起的位置误差进行了分析，相应的分析均建立在严密几何定位模型而非简化模型基础上。分析结果表明，在600km的轨道高度，由于卫星相对运动而

22、产生的激光光行差可引起地面约15m的平面位置偏差，在定位经过中必须加以考虑。此外激光指向角和其测量误差d需重点关注，为了使高程误差优于15cm，激光指向角应不大于1，d应优于1.5。文章结合ICESat/GLAS数据以及资源三号02星第一轨激光数据对几何模型进行了分析验证，结果表明该几何模型是正确的，对于ICESat/GLAS基于几何模型的计算精度与官方公布基本一致，能到达ICESat标称的0.15m的测高精度；对于国产的资源三号02星激光试验性载荷的第一轨数据，经该模型初步处理，剔除明显的错误点，采用参考DSM进行评价，其高程中误差为5.73m，与标称的精度还有一定差距，还有待继续深化分析研究。当前，对地观测卫星发展突飞猛进，搭载了激光测高仪的国产高分七号以及陆地生态系统碳监测卫星也在发射计划中，对激光测高卫星的研究即将进入一个蓬勃发展期，等待本文的研究能对推动国产激光测高卫星的发展做出一点奉献。