CN105842679B - 一种国产卫星激光高度计在轨几何标定方法及*** - Google Patents
一种国产卫星激光高度计在轨几何标定方法及*** Download PDFInfo
- Publication number
- CN105842679B CN105842679B CN201610255083.8A CN201610255083A CN105842679B CN 105842679 B CN105842679 B CN 105842679B CN 201610255083 A CN201610255083 A CN 201610255083A CN 105842679 B CN105842679 B CN 105842679B
- Authority
- CN
- China
- Prior art keywords
- mtd
- msub
- mrow
- mtr
- mover
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/497—Means for monitoring or calibrating
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明公开了一种国产卫星激光高度计在轨几何标定方法及***，所述方法通过卫星过顶前利用机载激光雷达或者GPS测量标定场的地形，采用激光雷达点云数据或者地面GPS测量数据拟合标定场地表面几何模型。将卫星激光高度计的激光点位于某一已知平面上作为约束条件，建立卫星激光高度计严密几何标定模型，利用最小二乘原理，求解卫星激光高度计的几何标定参数。本方法利用卫星过境时地面自然地形，可以很好地解决由于卫星激光高度计安装偏差以及激光测距误差造成的激光光斑定位精度不高这一问题，提高了国产卫星激光高度计的测量精度，具有较强的适用性和实用价值。
Description
技术领域
本发明属于遥感对地观测技术领域，特别涉及一种国产卫星激光高度计在轨几何标定方法及***。
背景技术
卫星激光高度计是在卫星平台上安装激光高度计以实现对地观测的一种主动遥感技术，作为深空探测最重要的载荷之一，具有激光脚印大小可控、激光光束指向精确以及高精度测量等优点，几乎所有的航天工程都搭载这种传感器。然而，激光高度计对地观测的测量精度受到很多因素影响，其中激光测距和激光光束指向会直接造成激光光斑位置产生偏移，激光高度计测距主要受到大气、测时等因素的影响，而安装误差会造成激光光束指向偏离天底点。
为了提高卫星激光高度计在轨测量精度，我们必须在过境时对卫星激光高度计进行在轨检校，以获取标定参数值，进而确定它的***误差，恢复卫星激光高度计照射到地面激光光斑的正确位置。
目前国际上还没有为卫星激光高度计专门制定检校标准、方法和流程。卫星激光高度计目前主要检校方法是基于红外信号探测器，即在检校场布设一系列规则红外发射器，捕捉卫星过境时的激光信号，根据激光光斑和红外发射器的相对位置关系确定激光光斑的真实坐标，将真实坐标值与定位模型计算得到的坐标进行对比检核。激光高度计在轨几何标定不但人工成本非常高，其可靠性也低，卫星过境时检校失误造成的经济损失巨大。
发明内容
本发明目的在于解决现有技术不足，提供一种国产卫星激光高度计在轨几何标定方法及***，将卫星激光高度计照射到地表面的地形特征作为检校约束条件，采用最小二乘方法计算***几何标定参数，克服了卫星激光高度计在轨检校成本较高等难题。
本发明通过以下技术方案实现：
一种卫星激光高度计在轨几何标定方法，所述方法包括以下步骤：
步骤1.1，针对检校场平面，利用机载激光雷达或者GPS测量，获取高精度激光点云或者GPS测量点，利用这些离散点构建数字表面模型，建立检校区域的平面方程。所述步骤1.1包括以下具体步骤：
步骤1.1.1，标定场地形数据采集。标定场地形数据采集可以采取机载激光雷达测量和GPS野外测量两种数据采集方式，选择较为晴朗的天气，在卫星过境前测量标定场地形，采集标定场地面数据。在进行地形数据采集时，需要保证点位分布均匀，且同一个平面至少包含三个以上的点；
步骤1.1.2，建立检校场内每个平面的平面方程，通过地面采集点平面拟合，计算各个平面方程的系数，标定场平面方程模型如下：
z＝aix+biy+ci
其中，i代表检校场内的平面标号，ai、bi、ci分别是第i个平面的平面方程系数，(x,y,z)代表第i个平面内所有的地面采集数据点集，具体形式为：
上式中，n代表某个平面内的所有点数量，将数据点集(x,y,z)代入标定场平面模型，计算标定场内各个平面方程的系数
步骤1.2，建立卫星激光高度计检校模型。在回波波形信号滤波处理后，根据激光高斯波形信号的峰值位置，确定激光发射点到地面光斑的起始距离。具体地说，激光测距模型可以表示为
式中为激光光斑在WGS-84坐标系下的地面坐标；是激光坐标系原点在WGS-84坐标系的坐标，可以从星上下传的轨道文件中获得；R为激光测距值；r为激光测距***误差；Tbody2WGS-84是卫星平台坐标系到WGS-84坐标系之间的转换矩阵，描述卫星运行时的姿态，由星上下传的姿态文件获得；Tlaser2body是卫星激光高度计载荷在卫星本体平台下的安装矩阵，描述载荷与卫星平台的相对关系，其大小会在卫星发射前精确测量；θx和θy分别是激光光束沿X轴和Y轴的偏转参数；θx、θy和r是卫星激光高度计在轨几何标定的待定三个参数。
代入平面方程
zs′＝axs′+bys′+c
其中，a、b、c可以从平面方程的系数P中确定；
步骤1.3，卫星激光高度计检校参数解算；将激光高度计测距模型代入平面方程展开后，得到
ax0+by0-z0+c+(R+r)·Tbody2WGS-84·Tlaser2body·(cosθxcosθy-asinθy+bsinθxcosθy)＝0令
f(r,θx,θy)＝ax0+by0-z0+c+(R+r)·Tbody2WGS-84·Tlaser2body·(cosθxcosθy-asinθy+bsinθxcosθy)将f进行泰勒函数一阶展开后，得到
其中，f0(r,θx,θy)是f函数的近似值，f函数用矩阵表示为
Ax+B＝0
其中，
采用迭代方式对未知数进行解算，即
式中分别为第k+1次和第k次参数估计值，为第k次解算的参数估计值的改正数，当满足相邻两次检校前后激光光斑之间的距离小于某一参考值时，迭代停止，用公式描述为
上式中τ为迭代停止的阈值，可以在***检校前根据经验参数来确定。
相应地，本发明提供一种卫星激光高度计在轨几何标定***，其特征在于，所述方法包括以下模块：数据采集模块、平面方程建模模块、激光高度计几何标定模块、数据误差纠正模块；其中
数据采集模块，采集激光测距数据、地面GPS数据以及激光雷达数据；
平面方程建模模块，利用GPS数据或者激光雷达数据生成各个平面方程，并获取每个平面的参数；
激光高度计几何标定模块，利用平面方程和激光测距模型联合解算，获取激光测距误差改正值和激光光束偏转角度；
数据误差纠正模块，将计算得到的激光标定参数代入激光测距方程，获取激光光斑纠正后的位置。
本方法利用卫星过境时地面自然地形，采用卫星激光高度计照射到地表面的地形特征作为检校约束条件，采用最小二乘方法计算***几何标定参数，可以很好地解决由于卫星激光高度计安装偏差以及激光测距误差，提高了国产卫星激光高度计的测量精度，克服了卫星激光高度计在轨检校成本较高等难题，具有较强的适用性和实用价值。
附图说明
图1是根据本发明实施例的一种国产卫星激光高度计在轨几何标定流程图；
图2是根据本发明实施例的卫星激光高度计测量原理示意图；
图3是根据本发明实施例的国产卫星激光高度计在轨几何标定***功能模块图。
具体实施方式
下面结合具体实施案例，参照附图对本发明作进一步详细的描述。
如图1所示，本发明提供一种国产卫星激光高度计在轨几何标定方法，其包括以下步骤：
步骤1.1，针对检校场平面，利用机载激光雷达或者GPS测量，获取高精度激光点云或者GPS测量点，利用这些离散点构建数字表面模型，建立检校区域的平面方程。所述步骤1.1包括以具体下步骤：
步骤1.1.1，标定场地形数据采集。标定场地形数据采集可以采取机载激光雷达测量和GPS野外测量两种数据采集方式，选择较为晴朗的天气，在卫星过境前测量标定场地形，采集标定场地面数据。在进行地形数据采集时，需要保证点位分布均匀，且同一个平面至少包含三个以上的点；
步骤1.1.2，建立检校场内每个平面的平面方程，通过地面采集点平面拟合，计算各个平面方程的系数，标定场平面方程模型如下：
z＝aix+biy+ci
其中，i代表检校场内的平面标号，ai,bi,ci分别是第i个平面的平面方程系数，(x,y,z)代表第i个平面内所有的地面采集数据点集，具体形式为：
上式中，n代表某个平面内的所有点数量，将数据点集(x,y,z)代入标定场平面模型，计算标定场内各个平面方程的系数
步骤1.2，建立卫星激光高度计检校模型。在回波波形信号滤波处理后，根据激光高斯波形信号的峰值位置，确定激光发射点到地面光斑的起始距离。具体地说，激光测距模型可以表示为
式中为激光光斑在WGS-84坐标系下的地面位置坐标；
是激光坐标系原点在WGS-84坐标系的坐标，可以从星上下传的轨道文件中获得；如图2所示，R为激光测距值；r为激光测距***误差；Tbody2WGS-84是卫星平台坐标系到WGS-84坐标系之间的转换矩阵，描述卫星运行时的姿态，由星上下传的姿态文件获得；Tlaser2body是卫星激光高度计载荷在卫星本体平台下的安装矩阵，描述载荷与卫星平台的相对关系，其大小会在卫星发射前精确测量；θx和θy分别是激光光束沿X轴和Y轴的偏转参数；θx、θy和r是卫星激光高度计在轨几何标定的待定三个参数。
代入平面方程
zs′＝axs′+bys′+c
其中，a、b、c可以从平面方程的系数P中确定；
步骤1.3，卫星激光高度计检校参数解算；将激光高度计测距模型代入平面方程展开后，得到
ax0+by0-z0+c+(R+r)·Tbody2WGS-84·Tlaser2body·(cosθxcosθy-asinθy+bsinθxcosθy)＝0令
f(r,θx,θy)＝ax0+by0-z0+c+(R+r)·Tbody2WGS-84·Tlaser2body·(cosθxcosθy-asinθy+bsinθxcosθy)将f进行泰勒函数一阶展开后，得到
其中，f0(r,θx,θy)是f函数的近似值，f函数用矩阵表示为
Ax+B＝0
其中，
采用迭代方式对未知数进行解算，即
式中分别为第k+1次和第k次参数估计值，为第k次解算的参数估计值的改正数，当满足相邻两次检校前后激光光斑之间的距离小于某一参考值时，迭代停止，用公式描述为
上式中τ为迭代停止的阈值，可以在***检校前根据经验参数来确定。
如图3所示，本发明提供一种卫星激光高度计在轨几何标定***，包括以下模块：数据采集模块、平面方程建模模块、激光高度计几何标定模块、数据误差纠正模块。所述数据采集模块主要是采集激光测距数据、地面GPS数据或者激光雷达数据；所述平面方程建模模块主要是利用GPS数据或者激光雷达数据生成各个平面方程，并获取每个平面的参数；所述激光高度计几何标定模块主要是利用平面方程和激光测距模型联合解算，获取激光测距误差改正值和激光光束偏转角度；所述数据误差纠正模块主要是将计算得到的激光标定参数代入激光测距方程，获取激光光斑纠正后的位置。
Claims (2)
1.一种卫星激光高度计在轨几何标定方法，其特征在于包括以下步骤：
步骤1，针对检校场平面，利用机载激光雷达或者GPS测量，获取高精度激光点云或者GPS测量点，利用这些离散点构建数字表面模型，建立检校区域的平面方程；所述步骤1包括以下具体步骤：
步骤1.1，标定场地形数据采集，标定场地形数据采集采取机载激光雷达测量和GPS野外测量两种数据采集方式，选择较为晴朗的天气，在卫星过境前测量标定场地形，采集标定场地面数据；在进行地形数据采集时，需要保证点位分布均匀，且同一个平面至少包含三个以上的点；
步骤1.2，建立检校场内每个平面的平面方程，采集地面点进行平面拟合，计算各个平面方程的系数，标定场平面方程模型如下：
z＝aix+biy+ci
其中，i代表检校场内的平面标号，ai,bi,ci分别是第i个平面的平面方程系数，(x,y,z)代表该平面内所有的地面采集数据点集，具体形式为：
<mrow>
<mo>&lsqb;</mo>
<mtable>
<mtr>
<mtd>
<mi>x</mi>
</mtd>
<mtd>
<mi>y</mi>
</mtd>
<mtd>
<mi>z</mi>
</mtd>
</mtr>
</mtable>
<mo>&rsqb;</mo>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mi>n</mi>
</msub>
</mtd>
<mtd>
<msub>
<mi>y</mi>
<mi>n</mi>
</msub>
</mtd>
<mtd>
<msub>
<mi>z</mi>
<mi>n</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
上式中，n代表某个平面内的所有点数量，将该数据点集代入标定场平面模型，计算标定场内各个平面方程的系数
<mrow>
<mi>P</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
</mtd>
<mtd>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
</mtd>
<mtd>
<msub>
<mi>c</mi>
<mi>i</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
步骤2，建立卫星激光高度计检校模型，在回波波形信号滤波处理后，根据激光高斯波形信号的峰值位置，确定激光发射点到地面光斑的起始距离；具体地说，激光测距模型表示为
<mrow>
<msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<msup>
<mi>S</mi>
<mo>&prime;</mo>
</msup>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>y</mi>
<msup>
<mi>S</mi>
<mo>&prime;</mo>
</msup>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>z</mi>
<msup>
<mi>S</mi>
<mo>&prime;</mo>
</msup>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mi>W</mi>
<mi>G</mi>
<mi>S</mi>
<mo>-</mo>
<mn>84</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mi>O</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>y</mi>
<mi>O</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>z</mi>
<mi>O</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mi>W</mi>
<mi>G</mi>
<mi>S</mi>
<mo>-</mo>
<mn>84</mn>
</mrow>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mi>R</mi>
<mo>+</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>b</mi>
<mi>o</mi>
<mi>d</mi>
<mi>y</mi>
<mn>2</mn>
<mi>W</mi>
<mi>G</mi>
<mi>S</mi>
<mo>-</mo>
<mn>84</mn>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>l</mi>
<mi>a</mi>
<mi>s</mi>
<mi>e</mi>
<mi>r</mi>
<mn>2</mn>
<mi>b</mi>
<mi>o</mi>
<mi>d</mi>
<mi>y</mi>
</mrow>
</msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>sin&theta;</mi>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>sin&theta;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>cos&theta;</mi>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>cos&theta;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>cos&theta;</mi>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中为激光光斑在WGS-84坐标系下的地面位置坐标；
是激光坐标系原点在WGS-84坐标系的坐标，从星上下传的轨道文件中获得；R为激光测距值；r为激光测距***误差；Tbody2WGS-84是卫星平台坐标系到WGS-84坐标系之间的转换矩阵，描述卫星运行时的姿态，由星上下传的姿态文件获得；Tlaser2body是卫星激光高度计载荷在卫星本体平台下的安装矩阵，描述载荷与卫星平台的相对关系，其大小会在卫星发射前精确测量；θx和θy分别是激光光束沿X轴和Y轴的偏转参数；θx、θy和r是卫星激光高度计在轨几何标定的待定三个参数；
将上式代入平面方程
zs′＝axs′+bys′+c
其中，a、b、c从平面方程的系数P中确定；
步骤3，卫星激光高度计检校参数解算；将激光高度计测距模型代入平面方程展开后，得到
ax0+by0-z0+c+(R+r)·Tbody2WGS-84·Tlaser2body·(cosθxcosθy-asinθy+bsinθxcosθy)＝0令
f(r,θx,θy)＝ax0+by0-z0+c+(R+r)·Tbody2WGS-84·Tlaser2body·(cosθxcosθy-asinθy+bsinθxcosθy)
将f进行泰勒函数一阶展开后，得到
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>&theta;</mi>
<mi>x</mi>
</msub>
<mo>,</mo>
<msub>
<mi>&theta;</mi>
<mi>y</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>&theta;</mi>
<mi>x</mi>
</msub>
<mo>,</mo>
<msub>
<mi>&theta;</mi>
<mi>y</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>f</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mi>d</mi>
<mi>r</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>f</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>x</mi>
</msub>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<msub>
<mi>d&theta;</mi>
<mi>x</mi>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>f</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>y</mi>
</msub>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<msub>
<mi>d&theta;</mi>
<mi>y</mi>
</msub>
</mrow>
其中，f0(r,θx,θy)是f函数的近似值，f函数用矩阵表示为
Ax+B＝0
其中，
<mrow>
<mi>A</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>f</mi>
<mn>1</mn>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>f</mi>
<mn>1</mn>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>x</mi>
</msub>
</mrow>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>f</mi>
<mn>1</mn>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>y</mi>
</msub>
</mrow>
</mfrac>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>f</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>f</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>x</mi>
</msub>
</mrow>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>f</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>y</mi>
</msub>
</mrow>
</mfrac>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>f</mi>
<mi>n</mi>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>r</mi>
</mrow>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>f</mi>
<mi>n</mi>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>x</mi>
</msub>
</mrow>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>f</mi>
<mi>n</mi>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&theta;</mi>
<mi>y</mi>
</msub>
</mrow>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mi>B</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>f</mi>
<mn>0</mn>
<mn>1</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>&theta;</mi>
<mi>x</mi>
</msub>
<mo>,</mo>
<msub>
<mi>&theta;</mi>
<mi>y</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>f</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>&theta;</mi>
<mi>x</mi>
</msub>
<mo>,</mo>
<msub>
<mi>&theta;</mi>
<mi>y</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>f</mi>
<mn>0</mn>
<mi>n</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>&theta;</mi>
<mi>x</mi>
</msub>
<mo>,</mo>
<msub>
<mi>&theta;</mi>
<mi>y</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mi>x</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>d</mi>
<mi>r</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>d</mi>
<msub>
<mi>&theta;</mi>
<mi>x</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>d&theta;</mi>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>A</mi>
<mi>T</mi>
</msup>
<mi>A</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>A</mi>
<mi>T</mi>
</msup>
<mi>B</mi>
</mrow>
采用迭代方式对未知数进行解算，即
<mrow>
<msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mover>
<mi>r</mi>
<mo>^</mo>
</mover>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mi>x</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mi>y</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mover>
<mi>r</mi>
<mo>^</mo>
</mover>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mi>x</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mi>y</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>d</mi>
<mi>r</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>d</mi>
<msub>
<mi>&theta;</mi>
<mi>x</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>d&theta;</mi>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
式中分别为第k+1次和第k次参数估计值，其中为第k次解算的参数估计值的改正数，当满足相邻两次检校前后激光光斑之间的距离小于某一参考值时，迭代停止，用公式描述为
<mrow>
<mrow>
<mo>||</mo>
<mrow>
<msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mover>
<mi>r</mi>
<mo>^</mo>
</mover>
<mi>sin</mi>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mi>r</mi>
<mo>^</mo>
</mover>
<mi>sin</mi>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mi>x</mi>
</msub>
<mi>cos</mi>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mover>
<mi>r</mi>
<mo>^</mo>
</mover>
<mi>cos</mi>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mi>x</mi>
</msub>
<mi>cos</mi>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mover>
<mi>r</mi>
<mo>^</mo>
</mover>
<mi>sin</mi>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mi>r</mi>
<mo>^</mo>
</mover>
<mi>sin</mi>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mi>x</mi>
</msub>
<mi>cos</mi>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mover>
<mi>r</mi>
<mo>^</mo>
</mover>
<mi>cos</mi>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mi>x</mi>
</msub>
<mi>cos</mi>
<msub>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>k</mi>
</msub>
</mrow>
<mo>||</mo>
</mrow>
<mo><</mo>
<mi>&tau;</mi>
</mrow>
上式中τ为迭代停止的阈值，在***检校前根据经验参数来确定。
2.一种用于实现权利要求1所述的卫星激光高度计在轨几何标定方法的卫星激光高度计在轨几何标定***，其特征在于，所述在轨几何标定***包括数据采集模块、平面方程建模模块、激光高度计几何标定模块、数据误差纠正模块；
数据采集模块采集激光测距数据、地面GPS数据或者激光雷达数据；平面方程建模模块利用GPS数据或者激光雷达数据生成各个平面方程，并获取每个平面的参数；激光高度计几何标定模块利用平面方程和激光测距模型联合解算，获取激光测距误差改正值和激光光束偏转角度；数据误差纠正模块将计算得到的激光标定参数代入激光测距方程，获取激光光斑纠正后的位置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610255083.8A CN105842679B (zh) | 2016-04-22 | 2016-04-22 | 一种国产卫星激光高度计在轨几何标定方法及*** |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610255083.8A CN105842679B (zh) | 2016-04-22 | 2016-04-22 | 一种国产卫星激光高度计在轨几何标定方法及*** |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105842679A CN105842679A (zh) | 2016-08-10 |
| CN105842679B true CN105842679B (zh) | 2018-02-16 |
Family
ID=56589960
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610255083.8A Active CN105842679B (zh) | 2016-04-22 | 2016-04-22 | 一种国产卫星激光高度计在轨几何标定方法及*** |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105842679B (zh) |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106780386B (zh) * | 2016-12-16 | 2019-11-26 | 武汉理工大学 | 一种三维激光扫描变形提取可靠性评价方法 |
| US10509120B2 (en) * | 2017-02-16 | 2019-12-17 | GM Global Technology Operations LLC | Lidar-radar relative pose calibration |
| CN106840212A (zh) * | 2017-04-12 | 2017-06-13 | 国家测绘地理信息局卫星测绘应用中心 | 基于地面激光光斑质心位置的星载激光在轨几何检校方法 |
| CN106960468B (zh) * | 2017-04-12 | 2019-07-09 | 武汉理工大学 | 一种三维激光扫描点云精度评价方法 |
| CN106871932B (zh) * | 2017-04-20 | 2019-12-24 | 自然资源部国土卫星遥感应用中心 | 基于金字塔搜索地形匹配的星载激光在轨指向检校方法 |
| CN107179534B (zh) * | 2017-06-29 | 2020-05-01 | 北京北科天绘科技有限公司 | 一种激光雷达参数自动标定的方法、装置及激光雷达 |
| CN107678013B (zh) * | 2017-09-14 | 2021-08-20 | 同济大学 | 远程激光雷达标定***与方法 |
| CN108226907B (zh) * | 2017-12-11 | 2020-05-19 | 武汉万集信息技术有限公司 | 用于激光测距设备的测距校准方法和装置 |
| CN108469624B (zh) * | 2018-03-26 | 2021-06-25 | 西安电子科技大学 | 基于激光探测器阵列超分辨的卫星定标方法 |
| CN110308436B (zh) * | 2018-03-27 | 2023-07-21 | 深圳砺剑天眼科技有限公司 | 一种多线激光扫描仪的激光光轴标校方法及*** |
| CN111366070B (zh) * | 2018-12-25 | 2022-04-29 | 苏州笛卡测试技术有限公司 | 一种复合式线激光测量***多轴空间坐标系标定方法 |
| CN111060139A (zh) * | 2019-12-31 | 2020-04-24 | 武汉大学 | 星载激光测高仪无场几何定标方法及*** |
| CN111366913A (zh) * | 2020-03-24 | 2020-07-03 | 上海卫星工程研究所 | 星载激光雷达光轴指向测量误差的标定方法 |
| CN111623802B (zh) * | 2020-07-01 | 2023-03-28 | 中国科学院上海技术物理研究所 | 一种基于足印影像的星载激光测高仪在轨指向检校方法 |
| CN116840851B (zh) * | 2023-07-05 | 2024-01-12 | 中国科学院空天信息创新研究院 | 一种星载对地激光测高仪地面探测器布设方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN2852150Y (zh) * | 2005-09-15 | 2006-12-27 | 中国科学院上海光学精密机械研究所 | 瑞利测风激光雷达的标定装置 |
| CN101566692B (zh) * | 2009-05-26 | 2011-04-13 | 吉林大学 | 利用卫星遥感数据中的云影信息检测云高的方法 |
| CN103364766B (zh) * | 2012-04-01 | 2015-02-04 | 中国科学院电子学研究所 | 星载InSAR***的外定标方法 |
| US9504425B2 (en) * | 2013-12-16 | 2016-11-29 | Verily Life Sciences Llc | Method of location coordination via wireless protocol between multiple devices |
| CN105068065B (zh) * | 2015-07-29 | 2017-06-27 | 武汉大学 | 星载激光测高仪在轨检校方法及*** |
| CN105512678A (zh) * | 2015-12-01 | 2016-04-20 | 中国科学院上海技术物理研究所 | 一种卫星激光高度计足印相机影像精确分类的方法 |
-
2016
- 2016-04-22 CN CN201610255083.8A patent/CN105842679B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105842679A (zh) | 2016-08-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105842679B (zh) | 一种国产卫星激光高度计在轨几何标定方法及*** | |
| CN106125069B (zh) | 一种基于指向角残差的星载激光测高仪指向角***误差标定方法 | |
| CN106291488B (zh) | 一种雷达标定误差校正方法 | |
| US20200103530A1 (en) | Method for extracting elevation control point with assistance of satellite laser altimetry data | |
| CN103217688B (zh) | 一种基于不规则三角网机载激光雷达点云平差计算方法 | |
| CN106990401B (zh) | 基于全波形机载激光雷达数据二类高程误差修正方法 | |
| CN102506824B (zh) | 一种城市低空无人机***生成数字正射影像图的方法 | |
| CN105068065B (zh) | 星载激光测高仪在轨检校方法及*** | |
| CN104061932B (zh) | 一种利用引力矢量和梯度张量进行导航定位的方法 | |
| CN108061889A (zh) | Ais与雷达角度***偏差的关联方法 | |
| CN102707726B (zh) | 一种无人机目标定位方法 | |
| CN103954953B (zh) | 一种基于数据驱动的机载激光雷达盲源误差补偿方法 | |
| CN106970398A (zh) | 顾及卫星遮挡条件的卫星可见性分析及星历预报方法 | |
| CN107179533A (zh) | 一种多参数机载LiDAR***误差自动检校方法 | |
| CN107886531A (zh) | 一种基于激光测距以及物方匹配的虚拟控制点获取方法 | |
| CN101839976B (zh) | 利用北斗一号终端进行雷达***误差估计的方法 | |
| CN105573318A (zh) | 基于概率分析的环境构建方法 | |
| CN102176003B (zh) | 一种机载激光雷达航测参数的优化设计方法 | |
| CN106382917A (zh) | 一种室内环境三维空间信息连续采集方法 | |
| CN105931200A (zh) | 一种小面阵星载tdi ccd相机的快速几何精校正方法 | |
| CN104655135A (zh) | 一种基于地标识别的飞行器视觉导航方法 | |
| CN110440762B (zh) | 一种多旋翼无人机山区航测成图的网格化像控点布设方法 | |
| CN105606128A (zh) | 一种星载激光高度计外场检校方法 | |
| CN110081905A (zh) | 一种基于单站光电经纬仪的光波大气折射误差计算方法 | |
| CN107991676A (zh) | 星载单航过InSAR***对流层误差校正方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |