星载激光测高仪在轨检校方法及系统的制作方法
【技术领域】
[0001] 本发明涉及星载激光测高仪在轨标定技术领域，特别涉及一种基于自然地表测距 残差的星载激光测高仪在轨检校技术方案，用于星载激光测高系统的在轨标定。
【背景技术】
[0002] 星载激光测高仪是一种主动式的测量设备，具备波形记录和分析功能的激光测高 仪通过采集激光测高仪的原始回波波形并提取波形的统计参量，结合测高仪的系统参数及 位置姿态信息，能够有效地获取被测目标的高度分布及地表信息。
[0003] 星载激光测高仪主要由三个部分组成，高精度定姿系统，高精度定位系统及测距 系统。其测量原理如图1所示，其中采用Satellite表示卫星，Footprint表示脚点，X ICRF、 YICRF、ZieRF分别为国际天球参考框架的X、Y、Z轴。实质上是利用矢量加法获得激光脚点的 位置坐标，即通过测距系统测定卫星与目标表面的距离（P )，结合高精度定姿系统测得的 激光指向角单位向量（u)得到测距矢量，进一步结合高精度定位系统测得的系统位置矢量 (r)，生成激光脚点位置矢量（R)。
[0004] 为了满足测高仪的科学测量目的，对测高系统的数据产品精度有着严格的要求。 激光脚点三维坐标数据是星载激光测高系统的核心数据产品，其平面及高程精度受到系统 误差的影响。由测高系统脚点定位原理可知，系统误差来自测高系统的高精度定姿系统，高 精度定位系统及测距系统三个部分及三个部分系统误差的组合，其中主要系统误差为系统 装配误差及测距误差。因此，要提高脚点定位精度，就必须对脚点定位过程中的系统误差进 行有效的剔除，即进行在轨检校。
[0005] 现有的检校方法有三种，一种是基于海洋扫描的测距残差分析方法，其原理是通 过卫星平台在平静海面的姿态机动，对测高系统的测距观测值进行残差（测距模型距离解 算值与测高系统距离实测值的差值）分析，达到检校指向角误差及测距误差的目的。该方 法能对指向角误差进行统计分析，但需要用到星上数据及地表的先验知识。第二种是现场 实时检校方法，其原理是通过独立的测量方式得到地表激光光斑的"坐标真值"，与测高系 统测得的激光点坐标值进行比较，从而达到对指向角误差进行检校的目的。该方法测量光 斑"坐标真值"的方式又分为两种，一种是在选定的标定场沿卫星地表轨迹布设探测器的方 式，另一种是用机载相机对激光光斑成像的方式，两种探测方式各有利弊，前者不受实施时 间及环境的限制，但需要对探测器进行设计及维护；后者能记录激光光斑能量分布，为回波 波形参数提供参考，但需精确设定飞机过境检校场的时间，且实施时间限制在月亮处于地 平线以下的夜间，避免噪声对相机成像的干扰。第三种是基于自然地物表面的检校方法，其 原理是基于脚点定位方程及地表模型这两个空间关系，即脚点定位方程确定的脚点坐标满 足地表模型方程。该方法对检校场的要求不再局限于平坦表面，而是自然地物表面，避免了 载体平台的姿态机动，但需要已知地表的先验模型。
【发明内容】
[0006] 针对现有技术存在的不足，本发明提出一种基于自然地物表面测距残差分析的星 载激光测高仪系统误差的在轨检校方法。
[0007] 本发明的技术方案提供一种星载激光测高系统在轨检校方法，包括以下步骤：
[0008] 步骤1、根据已知地表先验模型，建立星载激光测高仪的测距模型，实现如下，
[0009] 所述地表先验模型采用面元先验模型，设面元先验模型为，
[0010] z = ax+by+c
[0011] 其中，（X，y，Z)为平面上点的坐标，a、b分别为面元法向量的X及Y坐标，c为平面 在Z轴上的截距；
[0012] 设在站心坐标系下卫星位置坐标为（Xs，Ys，Z s)，激光测量参考坐标系由光学平台 坐标系绕其X轴旋转角度r及Y轴旋转角度p得到，建立测距模型如下，
[0013] R = (Zs-aXs-bYs-c) (cos (r) cos (p) +asin (p) cos (r) -bsin (r)) ^Rbias
[0014] 其中，R为卫星激光测高系统的测距值，Rblas为测距误差；
[0015] 步骤2、基于测距模型，建立系统误差与测距残差的关系如下，
[0017] 其中，R1为第i次测量时的R值，S R1为第i次测量时的测距残差，dRblas为测距 误差，dr为指向角横滚方向误差，dp为指向角俯仰方向误差， ei为第i次测量时的随机误 差；
[0018] 步骤3、利用实测距离值与测距模型计算距离值生成测距残差并对系统误差进行 检校，实现如下，
[0019] 设进行了共n次测量，i = 1，2,…，n，将系统误差与测距残差的关系写成矩阵形式 如下，
[0021] 根据上式解算得到测距误差、装配误差横滚方向及俯仰方向系统误差。
[0022] 而且，步骤3采用最小二乘方法实现解算。
[0023] 本发明相应提供一种星载激光测高系统在轨检校系统，包括以下模块：
[0024] 测距模型建立模块，用于根据已知地表先验模型，建立星载激光测高仪的测距模 型，实现如下，
[0025] 所述地表先验模型采用面元先验模型，设面元先验模型为，
[0026] z = ax+by+c
[0027] 其中，（x，y，z)为平面上点的坐标，a、b分别为面元法向量的X及Y坐标，c为平面 在Z轴上的截距；
[0028] 设在站心坐标系下卫星位置坐标为（Xs，Ys，Z s)，激光测量参考坐标系由光学平台 坐标系绕其X轴旋转角度r及Y轴旋转角度P得到，建立测距模型如下，
[0029] R = (Zs-aXs-bYs-c) (cos (r) cos (p) +asin (p) cos (r) -bsin (r)) ^Rbias
[0030] 其中，R为卫星激光测高系统的测距值，Rblas为测距误差；
[0031] 残差关系构建模块，用于基于测距模型，建立系统误差与测距残差的关系如下，
[0033] 其中，R1为第i次测量时的R值，S R1为第i次测量时的测距残差，dRblas为测距 误差，dr为指向角横滚方向误差，dp为指向角俯仰方向误差， ei为第i次测量时的随机误 差；
[0034] 校验模块，用于根据测距残差并对系统误差进行检校，实现如下，
[0035] 设进行了共n次测量，i = 1，2,…，n，将系统误差与测距残差的关系写成矩阵形式 如下，
[0037] 根据上式解算得到测距误差、装配误差横滚方向及俯仰方向系统误差。
[0038] 而且，校验模块采用最小二乘方法实现解算。
[0039] 与现有技术相比，本发明相较于海洋扫描的测距残差分析方法，都利用了残差分 析进行检校，但检校场不再局限于平静的海面，且避免了姿态机动；相较于自然地物表面的 检校方法，都利用了自然地物表面作为检校场，但建立的误差模型更为直观的说明了装配 误差对测距值的影响。本发明技术方案能够很好地检校系统误差，使得估计偏差RSS值要 优于10%。。
【附图说明】
[0040] 图1现有技术中星载激光测高仪的星载激光测量原理；
[0041] 图2本发明实施例的站心坐标系下的激光测距示意图。
【具体实施方式】
[0042] 为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施案例，并参 照附图，对本发明进一步详细说明。
[0043] 本发明提供一种基于自然地物表面测距残差分析的星载激光测高仪系统误差的 在轨检校方法，利用自然地表为在轨检校场，基于自然地表的测距模型，使用测距残差分析 方法，对测高系统的系统误差，如装配误差及测距误差进行在轨检校。具体实施时可采用计 算机软件技术实现自动运行流程。实施例流程包括以下步骤：
[0044] 步骤1、在已知地表先验模型的前提下，建立星载激光测高仪的测距模型。
[0045] 测高系统激光器装配在卫星光学平台上，以光学平台坐标系（包括X^、Y^、Zi3bf 轴）为参考坐标系，理论上其指向为天底方向，由于装配问题、发射过程震动及卫星运行时 环境温度变化等因素影响，会使得激光指向偏离天底方向，从而使激光指向存在误差。此时 激光测量参考坐标系（包括轴）可由光学平台坐标系绕Xi3bf轴旋转角度r (顺 时针为正）及Yi3bf轴旋转角度p得到，此时激光偏离天底指向角为P，如图2所示，其中， X TRF、YTRF、ZTRF分别为站心坐标系下的X、Y、Z轴，Surface Model表示地表模型。由几何关系 可得激光指向角0与横滚角r及俯仰角p有如下关系式。
[0046] cos ( 0 ) = cos (r) cos (p) (1)
[0047] 在站心坐标系下，假设激光脚点坐标为（Xf，Yf，Z f)，卫星位置坐标为（Xs，Ys，Zs)，卫 星在站心坐标系下相对参考面的轨道高度为H，则：
[0048] Xf= X s+Rcos (r) sin (p) (2