CN104614730B - 一种利用月面激光反射阵坐标标校月球高程模型的方法 - Google Patents
一种利用月面激光反射阵坐标标校月球高程模型的方法 Download PDFInfo
- Publication number
- CN104614730B CN104614730B CN201510027624.7A CN201510027624A CN104614730B CN 104614730 B CN104614730 B CN 104614730B CN 201510027624 A CN201510027624 A CN 201510027624A CN 104614730 B CN104614730 B CN 104614730B
- Authority
- CN
- China
- Prior art keywords
- elevation
- lunar
- elevation model
- laser
- laser reflection
- 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
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 title claims abstract description 17
- 230000009897 systematic effect Effects 0.000 claims abstract description 15
- 238000012937 correction Methods 0.000 claims abstract description 7
- 238000005259 measurement Methods 0.000 claims description 15
- 238000003491 array Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 4
- 239000000523 sample Substances 0.000 description 2
- 241001673112 Citrus clementina Species 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/06—Systems determining position data of a target
- G01S17/08—Systems determining position data of a target for measuring distance only
-
- 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)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开一种利用月面激光反射阵坐标标校月球高程模型的方法包括:利用激光测高数据建立月球高程模型;计算月面激光反射阵精确高程值;计算月球高程模型的系统性偏差均值;对月球高程模型进行偏差修正。采用本发明的技术方案,可以有效消除高程模型的系统性偏差,精确获取月球高程模型。
Description
技术领域
本发明属于月球行星科学技术领域,尤其涉及一种利用月面激光反射阵坐标标校月球高程模型的方法。
技术背景
月球高程模型描述的是月球表面地形相对于某参考正球(一般取为半径为1738km的正球)的高度变化,反映的是月球表面的起伏状况,亦称为月球地形模型。目前月球高程模型均是通过月球探测器所携带的卫星激光高度计获取的测高数据反演得到,因此获取到的模型精度不仅与探测器轨道的特征相关,还与高度计自身的属性相关。我国嫦娥一号卫星自2007年10月份至2008年7月份共获取7,247,000个有效激光测高月面高程点,激光测高数据可以有效反映月球的地形地貌,利用嫦娥一号激光测高数据建立了月球地形模型CLTM-s01,日本学者指出该模型存在系统性偏差。事实上,不止国内的月球地形模型存在该问题,较为常见的月球地形模型如Clementine,ULCN均有类似问题。研究表明,激光高度计在数据预处理中,首先需要对测距值进行晶振引起的系统误差修正,如嫦娥卫星携带的高度计是通过测距值乘以修正因子0.9645得到的卫星至月球表面距离,如果修正因子存在千分之一的误差会引起高程百米的误差,该误差无法通过轨道或是激光测高数据进行消除。
发明内容
本发明要解决的技术问题是,提供提出一种利用月面激光反射阵坐标标校月球高程模型的方法,可以有效消除高程模型的系统性偏差,精确获取月球高程模型。
为解决上述问题,本发明采用如下的技术方案:
一种利用月面激光反射阵坐标标校月球高程模型的方法包括以下步骤:
步骤S1、利用激光测高数据建立月球高程模型
设格网内有n个观测点,则
其中,hk为曲面给出的高程值,和分别为激光测高数据点经度与格网中心点经度之差、激光测高数据点纬度与格网中心点纬度之差,N为曲面拟合次数,N≤3,ci,j为曲面多项式系数,i,j分别为经度差的指数、纬度差的指数,且i+j≤N,
设c=[c00,c10,…,cij,…,c0N]T,H=[h1,h2,…,hn]T,则,
c=(ATA)-1ATH
其中,参数c中的C00分量为格网中心点的高程值,所有格网中心点的高程值构成格网化的月球高程模型;
步骤S2、计算月面激光反射阵精确高程值
计算第i′个月面激光反射阵的精确高程值为:其中,为第i′个月面激光反射阵的精确坐标;
步骤S3、计算月球高程模型的系统性偏差均值
计算月球高程模型的系统性偏差均值S为:其中,n′为激光反射阵个数,hi′为第i′个月面激光反射阵所处经纬度处的高程值,为激光反射阵的精确坐标所处的经纬度;
步骤S4、对月球高程模型进行偏差修正
将月球高程模型所有格网中心点的高程值均减去偏差均值S,得到修正后的高程值,修正后的格网中心高程值为c00-S。
本发明利用月球高程模型计算月面激光反射阵的精确坐标所处经纬度的高程值;然后比较激光反射阵高程与精确高程的差异,得到系统性偏差均值;最后利用系统性偏差均值对月球高程模型进行偏差修正,采用本发明的技术方案,可以有效消除高程模型的系统性偏差,精确获取月球高程模型。
附图说明
图1为激光测高原理示意图;
图2为格网中心点与相邻节点的关系示意图;
图3为利用CE-1激光测距直接计算得到的高程模型与LOLA-720的差异示意图;
图4为利用本发明获取的高程模型与LOLA-720模型的差异示意图。
具体实施方式
本发明提供一种利用月面激光反射阵坐标标校月球高程模型的方法,包括以下步骤:
步骤S1、利用激光测高数据建立月球高程模型
如图1所示,给出了月球卫星激光测高、卫星轨道和姿态数据计算月面高程值的主要原理。探测器携带的激光高度计发射一组激光束至月表,经月表反射后回波再返回至高度计,高度计通过记录激光束接收信号与发射信号的时间差获取探测器至月表的距离。设为某观测时刻激光高度计的观测矢量,由激光高度计测距值和姿态数据确定;为观测时刻卫星在月心坐标系下的位置矢量,由卫星轨道给出;R为月球参考半径。利用和得到观测时刻光斑中心点在月心坐标系中的位置矢量利用和月球参考半径R的关系可得到月面光斑中心点对应的月面高程h和经纬度
其中,(x,y,z)为在月心坐标系下的坐标分量。
对所有获取的激光测高数据进行如上处理,获取月表离散的高程值,但是上面获取的高程值在空间的分布是不均匀的。通常月面高程模型为按经纬度均分分布的高程值,这就需要平滑处理。根据高程模型分辨率要求设置格网间隔,用该格网中的所有高程数据拟合产生网格中心点的高程数据。格网点高程采用最小二乘曲面拟合的方式给出,对格网内n个观测点,可建立以下曲面方程:
其中,hk为曲面给出的高程值,和分别为激光测高数据点经度与格网中心点经度之差、激光测高数据点纬度与格网中心点纬度之差,N为曲面拟合次数,N≤3,ci,j为曲面多项式系数,i,j分别为经度差的指数、纬度差的指数,且i+j≤N。
令c=[c00,c10,…,cij,…,c0N]T,H=[h1,h2,…,hn]T,
则公式(1)可表示成,
H=Ac (2)
根据最小二乘原理,公式(2)的解为,
c=(ATA)-1ATH (3)
其中,参数c中的c00分量为格网中心点的高程值,所有格网中心点的高程值构成了格网化的月球高程模型。
在具体拟合时,根据格网点中激光测高数据个数n来确定拟合采用曲面次数,如表1所示:
表1拟合方法设置
观测数据个数 | 拟合方法 |
n<3 | 直接利用高程数据的平均值 |
3≤n<6 | 利用平面拟合,N=1 |
6≤n<10 | 二次曲面拟合,N=2 |
n≥10 | 三次曲面拟合,N=3 |
步骤S2、计算月面激光反射阵精确高程值
计算第i′个月面激光反射阵的精确高程值为:
其中,为第i′个月面激光反射阵的精确坐标,由表2给出。
表2月面激光反射阵精确坐标
步骤S3、计算月球高程模型的系统性偏差均值
利用第步骤S1中计算得到的月球高程模型,使用公式(1)插值计算第i′个月面激光反射阵所处经纬度处的高程值hi′,其中,为激光反射阵的精确坐标所处的经纬度。
如图2所示,根据激光反射阵的经纬度确定月球高程模型中与相邻的4个节点这4个节点处高程值分别为 则月面激光反射阵高程值hi′为:
其中,Hj′,j′=1,2,3,4为相邻4个节点的高程值,
W1=α2β2(9-6α-6β+4αβ)
W2=(1-α)2β2(9-6(1-α)-6β+4(1-α)β)
W3=(1-α)2(1-β)2(9-6(1-α)-6(1-β)+4(1-α)(1-β))
W4=α2(1-β)2(9-6α-6(1-β)+4α(1-β))
其中,α,β为月面激光反射阵坐标的函数,
比较计算激光反射阵的高程值与精确高程值的差异,统计月球高程模型的系统性偏差均值S,偏差均值S则为月球高程模型的系统偏差标校值,
其中,n′为月面激光反射阵个数。
步骤S4、对月球高程模型进行偏差修正
将月球高程模型所有格网中心点的高程值均减去偏差均值S,得到修正后的高程值,修正后的格网中心高程值为c00-S。
利用本发明方法对嫦娥一号卫星获取的月球高程模型进行标校,并使用美国最新发布的LOLA-720模型作为基准,用于模型精度的比较验证。图3为未利用该方法建立的高程模型与LOLA-720模型的差异。图4为利用本发明方法标校的高程模型与LOLA-720模型的差异,消除了整体的系统性偏差。比较图3-4,利用本发明方法进行高程模型标校,可以有效去除激光测高设备系统偏差引入的误差,获取精确的高程模型。
说明书中描述的只是该发明的具体实施方式。虽然结合附图描述了本发明的实施方式,但是本领域内熟练的技术人员可以在所附权利要求的范围内做出各种变形或修改。
Claims (2)
1.一种利用月面激光反射阵坐标标校月球高程模型的方法,其特征在于,包括以下步骤:
步骤S1、利用激光测高数据建立月球高程模型
设格网内有n个观测点,则
其中,hk为曲面给出的高程值,和分别为激光测高数据点经度与格网中心点经度之差、激光测高数据点纬度与格网中心点纬度之差,N为曲面拟合次数,N≤3, cij为曲面多项式系数,i,j分别为经度差的指数、纬度差的指数,且i+j≤N,c为曲面多项式;
设c=[c00,c10,…,cij,…,c0N]T,H=[h1,h2,…,hn]T,则,
c=(ATA)-1ATH
其中,系数c00分量为格网中心点的高程值,所有格网中心点的高程值构成格网化的月球高程模型;
步骤S2、计算月面激光反射阵精确高程值
计算第i′个月面激光反射阵的精确高程值为:其中,为第i′个月面激光反射阵的精确坐标,R为月球参考半径;
步骤S3、计算月球高程模型的系统性偏差均值
计算月球高程模型的系统性偏差均值S为:其中,n′为激光反射阵个数,hi′为第i′个月面激光反射阵所处经纬度处的高程值,为激光反射阵的精确坐标所处的经纬度;
步骤S4、对月球高程模型进行偏差修正
将月球高程模型所有格网中心点的高程值均减去偏差均值S,得到修正后的高程值,修正后的格网中心点的高程值为c00-S。
2.如权利要求1所述的一种利用月面激光反射阵坐标标校月球高程模型的方法,其特征在于,步骤S3中,根据激光反射阵的经纬度确定月球高程模型中与相邻的4个节点所述4个节点处高程值分别为H1 月面激光反射阵高程值hi′为:其中,Hj′,j′=1,2,3,4为相邻4个节点的高程值,
其中,α,β为月面激光反射阵坐标的函数,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510027624.7A CN104614730B (zh) | 2015-01-20 | 2015-01-20 | 一种利用月面激光反射阵坐标标校月球高程模型的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510027624.7A CN104614730B (zh) | 2015-01-20 | 2015-01-20 | 一种利用月面激光反射阵坐标标校月球高程模型的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104614730A CN104614730A (zh) | 2015-05-13 |
CN104614730B true CN104614730B (zh) | 2018-02-02 |
Family
ID=53149264
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510027624.7A Active CN104614730B (zh) | 2015-01-20 | 2015-01-20 | 一种利用月面激光反射阵坐标标校月球高程模型的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104614730B (zh) |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101344391B (zh) * | 2008-07-18 | 2011-05-11 | 北京工业大学 | 基于全功能太阳罗盘的月球车位姿自主确定方法 |
CN101950436A (zh) * | 2010-09-29 | 2011-01-19 | 中国科学院国家天文台 | 利用激光高度计数据制作数字高程模型的方法 |
CN102313533B (zh) * | 2011-07-28 | 2012-12-19 | 东南大学 | 月球激光高度计在轨探测数据误差分析的方法 |
PT2641833T (pt) * | 2012-03-20 | 2018-01-03 | Airbus Defence & Space Ltd | Um método de ocultação solar |
CN102645229B (zh) * | 2012-04-06 | 2014-06-25 | 东南大学 | 星载激光高度计在轨月球高程探测不确定度评定方法 |
-
2015
- 2015-01-20 CN CN201510027624.7A patent/CN104614730B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104614730A (zh) | 2015-05-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110275184B (zh) | 一种gnss掩星电离层残差修正方法、系统、设备及存储介质 | |
CN103926589B (zh) | 星载激光测高系统固体地表目标平面和高程精度检测方法 | |
CN113189620B (zh) | 一种gnss掩星临近空间气候数据反演方法及系统 | |
CN106052717B (zh) | 一种利用精密弹道的电波折射修正效果标定方法 | |
CN106547002B (zh) | 一种提高rtk定位精度的选星方法 | |
CN110672031B (zh) | 一种同时以点与面特征约束的三维激光扫描的检校方法 | |
CN109959898B (zh) | 一种座底式水声被动定位基阵自校准方法 | |
CN107179533A (zh) | 一种多参数机载LiDAR系统误差自动检校方法 | |
CN112013874B (zh) | 一种基于激光足印先验坐标的星载激光测高仪在轨标定方法 | |
Shi et al. | Calibrating the scale of the NRLMSISE00 model during solar maximum using the two line elements dataset | |
CN105783944B (zh) | 太阳敏感器标定方法及系统 | |
CN111025243A (zh) | 一种基于参数优化的大气折射误差实时修正方法 | |
CN113074703A (zh) | 多源卫星高度计的数据处理方法、装置、设备及存储介质 | |
CN111060139A (zh) | 星载激光测高仪无场几何定标方法及系统 | |
Zeitlhöfler et al. | Station-dependent satellite laser ranging measurement corrections for TOPEX/Poseidon | |
Tang | Mathematical methods for camera self-calibration in photogrammetry and computer vision | |
CN103727962B (zh) | 大视场红外光电经纬仪精度标定方法 | |
CN104614730B (zh) | 一种利用月面激光反射阵坐标标校月球高程模型的方法 | |
CN109100705A (zh) | 星载激光测高仪在轨标定模型中权矩阵的确定方法 | |
CN109764892B (zh) | 近红外星等的校准方法 | |
CN114722455A (zh) | 一种联合全站仪和激光跟踪仪的三维工程控制网构建方法 | |
CN111044489B (zh) | 一种基于多波长测量获得大气折射率高度分布剖面的方法 | |
CN109506645B (zh) | 一种星敏感器安装矩阵地面精测方法 | |
JP2014048075A (ja) | Gnssによる位置計測装置 | |
WO2020178984A1 (ja) | 情報処理装置および測位補強情報送信方法 |
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 |