CN106643804A - 一种预先确定星载激光测高仪的足印位置的方法 - Google Patents
一种预先确定星载激光测高仪的足印位置的方法 Download PDFInfo
- Publication number
- CN106643804A CN106643804A CN201611253206.0A CN201611253206A CN106643804A CN 106643804 A CN106643804 A CN 106643804A CN 201611253206 A CN201611253206 A CN 201611253206A CN 106643804 A CN106643804 A CN 106643804A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- laser
- sin
- cos
- sensing
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
-
- 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)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Manufacturing & Machinery (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种预先确定星载激光测高仪的足印位置的方法,该方法结合了传统的光学相机成像的原理,并充分顾及了激光测高仪自身的测高机理,足印位置确定结果能够较好地符合卫星真实的飞行过境情况,具有较高的计算精度,并且极大限度地减少不必要的人力物力,节省了大量的试验开销。
Description
技术领域
本发明涉及星载激光测高仪在轨几何检校领域,特别涉及一种预先确定星载激光测高仪的足印位置的方法,应用于星载激光测高仪在轨几何检校过程中的激光地面探测器布设位置的预报。
背景技术
在轨几何检校试验的成功开展是星载激光测高仪获得高精度测高数据的必要条件之一。星载激光测高仪足印位置的预报是在轨几何检校的第一步,也是极为关键的一步。预报的精度不仅关系着整个检校试验外业布设工作的难易程度,甚至决定着试验的成败。
目前,仅有美国ICESat-GLAS(Ice Cloud Elevation Satellite GeoscienceLaser Altimeter System)开展过相关的研究和工程实践,其激光脉冲发射重频为40Hz,激光沿轨相邻两个足印中心间距为170米左右,地面足印平均大小为65米左右,沿轨布设大于170米的探测器阵列,必然可以捕获到光斑。由于其硬件的优势,使其在足印预报方法上并不需要开展过多的研究。但是其高重频是根据卫星的使用目的而设计的,即对冰川的厚度进行高精度测量;而对于光学测绘卫星而言,星上本身携带了高分辨率相机、姿态测控***和轨道确定***等多个载荷,激光测高仪的脉冲出射频率必定要控制在一定的范围内,以保证整星的能源分配问题,且用于符合光学测绘的激光测高数据在实际应用中也不需要过高的激光出射频率。低频率的脉冲出射为星载激光测高仪的在轨几何检校增加了试验难度,最突出地体现在激光足印的预报上。由于国际上暂时没有涉及该方面技术的测绘卫星在轨运行,激光足印预报方面的研究几乎属于空白。
因此,本领域需要一种能够准确预报星载激光测高仪足印位置的方法。
发明内容
为此,本发明提出一种预先确定星载激光测高仪的足印位置的方法,弥补了目前激光测高仪高精度足印位置预报方法的空白,并解决了现有简单预报技术的限制和缺陷导致的一个或多个问题。
本发明另外的优点、目的和特性,一部分将在下面的说明书中得到阐明,而另一部分对于本领域的普通技术人员通过对下面的说明的考察将是明显的或从本发明的实施中学到。通过在文字的说明书和权利要求书及附图中特别地指出的结构可实现和获得本发明目的和优点。
本发明提供了一种预先确定星载激光测高仪的足印位置的方法,其特征在于,所述方法包括以下步骤:
步骤1,计算激光测高仪在卫星本体坐标系中的指向Mb_l,步骤1具体包括以下子步骤:
步骤1.1,预估激光测高仪与平台之间的安装角(αr,βr);
步骤1.2,基于步骤1.1预估的安装角,计算激光测高仪在卫星本体坐标系中的指向Mb_l;
步骤2,将激光测高仪的指向从卫星本体坐标系旋转至轨道坐标系Mo_l;
步骤3,将激光测高仪的指向从轨道坐标系旋转至地球固定坐标系,获得激光测高仪在地球固定坐标系下的指向MT_l;
步骤4,通过计算地球固定坐标系下的激光测高仪的指向与地球椭球的交点,确定星载激光测高仪的足印位置,其中,所述步骤4具体包括以下子步骤:
步骤4.1,根据下述公式计算激光指向与地球椭球相交的地球固定坐标系下的地面点坐标(X1,Y1,Z1):
其中,MT=[X1 Y1 Z1]T,X1、Y1和Z1为待求的激光指向与地球椭球相交的地球固定坐标系下的地面点坐标,μ为待求的比例系数,MT_l为激光测高仪在地球固定坐标系下的指向,P0为卫星激光质心在地球固定坐标系的坐标值,a和b分别为地球椭球长半轴和短半轴长度,h为地面大地高程;
步骤4.2,根据下述公式将地球固定坐标系下的地面点坐标(X1,Y1,Z1)转换成大地经纬度坐标系下的地面点坐标(L,B,H):
其中,e1和e2分别为第一偏心率和第二偏心率,a和b分别为地球椭球长半轴和短半轴长度,X1、Y1、Z1为通过步骤4.1计算出的激光指向与地球椭球相交的地球固定坐标系下的地面点坐标。
本发明还提供了一种计算机程序,其中,当执行所述计算机程序时,执行本发明任一实施例所述的方法。
本发明针对星载激光测高仪足印位置预报领域的空白,提供了一种预先确定星载激光足印位置的方法,该方法结合了传统的光学相机成像的原理,并充分顾及了激光测高仪自身的测高机理,预报结果能够较好地符合卫星真实的飞行过境情况,具有较高的预报精度,并且极大限度地减少不必要的人力物力,节省了大量的试验开销。排除激光指向角、卫星平台姿态预报数据和卫星质心轨道预报数据精度损失的情况下,该方法的星载激光测高仪足印位置预报精度能够控制在5米。
附图说明
图1为根据本发明实施例的确定星载激光测高仪的足印位置的方法的流程图。
具体实施方式
下面参照附图对本发明进行更全面的描述,其中说明本发明的示例性
实施例。
根据本发明的实施例,可以随机抽取星载激光测高仪进行足印位置预报模型的构建;优选的,本发明的实施案例以资源三号02星的激光测高仪的足印位置预报模型构建为例,但并不限于此。
如图1所示,激光测高仪预报模型的构建和光学严密几何成像模型的构建存在同样的过程,涉及到各个坐标系之间的旋转,通过各坐标系的旋转最终将星上指向传输至地面点,本发明所提出的星载激光测高仪足印位置预报模型构建方法具体包括以下步骤:
步骤1,计算激光测高仪在卫星本体坐标系中的指向Mb_l,步骤1具体包括以下子步骤:
步骤1.1,预估激光测高仪与平台之间的安装角(αr,βr);
由于卫星发射升空过程中的剧烈震荡和在轨后太空环境与地面实验室的差异较大,会引起激光测高仪与平台之间的安装关系发生变化,导致激光测高仪的安装角发生变化,因此需要计算真实的激光测高仪与平台之间安装角。
参考资源三号02星激光测高仪设计指标,获取激光测高仪的实验室标定的安装角α=90.86229°,β=0.861985°。以实验室标定安装角为基础,采用基于地形匹配的安装角预估方法计算真实的指向角。这里采用金字塔计算策略,逐步缩小计算区域和遍历间隔,将时间限定在一定的区间内,保证计算的精度。将金字塔分为3层来进行计算,第一层金字塔的区域间隔为2°,遍历间隔为0.1°,第二层金字塔的区域间隔为0.4°,遍历间隔为1′,第三层金字塔的区域间隔为0.02°,遍历间隔为1″,具体的计算公式如下式所示:
其中,(αr,βr)为预估的安装角(其中,αr为激光指向与本体坐标系X轴之间的夹角,βr为激光指向与本体坐标系Z轴之间的夹角),表示从1层至3层逐层计算,αi、βi表示第i层遍历的安装角,表示括号内第i层取最小值时对应的αi、βi,并将其赋值给α0、β0,当i=1首次计算时,α0、β0的值取实验室测量值90.86229°、0.861985°,q(i)为第i层遍历间隔,p(i)为第i层的区域间隔,表示对括号内进行遍历,行遍历序号m,列遍历序号n,E(αi,βi)为αi、βi指向下计算某轨的高程值与已知高程值之间的标准差,具体的高程值计算公式可参照激光测高仪的严格几何测高模型,已知高程采用公开全球范围的AW3D30数据。
以资源三号02星的第382轨激光数据为例,计算过程如下:首先进行第一层金字塔遍历,区域间隔为2°,遍历间隔为0.1°,
经过第一层的遍历,可以获得第一层的预估结果α0、β0分别为89.9°、0.0°,下面进行第二层金字塔的遍历,区域间隔为0.4°,遍历间隔为1′,
经过第二层的遍历,可以获得第二层的预估结果α0、β0分别为89.945°、0.046°,下面进行第三层金字塔的遍历,区域间隔为0.02°,遍历间隔为1″,
通过上式可以求得预估的安装角αr=89.94445°,βr=0.04692°
步骤1.2,基于步骤1.1预估的安装角,根据下述公式计算激光测高仪在卫星本体坐标系中的指向Mb_l:
式中,Mb_l为步骤1待求的激光测高仪在卫星本体坐标系中的指向,αr、βr为步骤1.1求得的激光测高仪的安装角。
步骤1.1求得的激光测高仪安装角预估值代入到上式可以求得预估指向
步骤2,将激光测高仪的指向从卫星本体坐标系旋转至轨道坐标系Mo_l。
根据光束的物理传播路径,传输的第一步是将激光测高仪的指向从卫星本体坐标系旋转至轨道坐标系,步骤2具体包括以下子步骤:
步骤2.1,估计卫星平台的本体坐标系相对于轨道坐标系的姿态值数值单位为弧度,其中,roll、pitch、yaw分别为卫星本体坐标系相对于轨道坐标系的三轴欧拉角:横滚角、俯仰角和偏航角,即,roll为卫星本体坐标系相对于轨道坐标系的横滚角,pitch为卫星本体坐标系相对于轨道坐标系的俯仰角,yaw为卫星本体坐标系相对于轨道坐标系的偏航角。
一般情况下,由于卫星姿控***的高精度实时控制,平台的本体坐标系相对于轨道坐标系的姿态变化较小,因此这里采用了取零值的方法,即三轴姿态均为零
步骤2.2,根据下述公式计算本体坐标系相对于轨道坐标系的旋转矩阵Mo_b;
a1=cos(pitch)*cos(yaw)
a2=-cos(pitch)*sin(yaw)
a3=sin(pitch)
b1=-sin(roll)*sin(pitch)*cos(yaw)+cos(roll)*sin(yaw)
b2=sin(roll)*sin(pitch)*sin(yaw)+cos(roll)*cos(yaw)
b3=sin(roll)*cos(pitch)
c1=-cos(roll)*sin(pitch)*cos(yaw)-sin(roll)*sin(yaw)
c2=cos(roll)*sin(pitch)*sin(yaw)-sin(roll)*cos(yaw)
c3=cos(roll)*cos(pitch)
将步骤2.1得到的卫星平台的本体坐标系相对于轨道坐标系的姿态值(即,)代入到Mo_b,可得本体坐标系相对于轨道坐标系的旋转矩阵
本实施例的步骤2.1优选的采用了的方法,即代入到步骤2.2计算。
步骤2.3,根据下述公式将激光测高仪的指向由本体坐标系转换到轨道坐标系;
Mo_l=Mo_b*Mb_l
其中,Mo_l为待求的激光测高仪在轨道坐标系下的指向,Mo_b为步骤2.2计算的本体坐标系相对于轨道坐标系的旋转矩阵,Mb_l为步骤1.2计算的激光测高仪的指向。
将步骤1.2和步骤2.2计算得到的激光测高仪载荷与平台之间的预估指向与本体坐标系相对于轨道坐标系的旋转矩阵相乘,可以得到激光测高仪的指向由本体坐标系转换到轨道坐标系的向量
步骤3,将激光测高仪指向从轨道坐标系旋转至地球固定坐标系。
根据光束的物理传播路径,传输的第二步是将激光测高仪指向从轨道坐标系旋转至地球固定坐标系,步骤3具体包括以下子步骤:
步骤3.1,计算卫星激光质心在地球固定坐标系的坐标值P0=[X0 Y0 Z0]T和位移速度值V0=[VX0 VY0 VZ0]T;
根据本发明的优选实施例,步骤3.1的具体计算过程如下:
式中,为指定时刻的惯性系卫星位置关于时间的二阶导数即加速度,为待求量t为某一指定时刻,r为指定时刻的惯性系卫星位置,v为指定时刻的惯性系卫星速度,CD为大气阻力系数,CR为太阳光压系数,aiR为作用在卫星轨道径向上的第i段经验加速度,aiT为作用在卫星轨道切向上的第i段经验加速度,aiN为作用在卫星轨道法向上的第i段经验加速度,eR(t)为指定时刻卫星轨道径向方向单位向量,eT(t)为指定时刻卫星轨道切向方向单位向量,eN(t)为指定时刻卫星轨道法向方向单位向量,ti为第i段经验加速度参数开始时刻,ti+1为第i+1段经验加速度开始时刻,也即第i段经验加速度结束时刻。
参考时刻卫星位置速度、大气阻力系数CD、太阳光压系数CR、经验加速度参数aiR、aiT、aiN等来自于预报区间之前的精密轨道拟合。
通过上式可求得新的加速度因此可以利用下式求待估参数(P0,V0),已知初值位置和速度P0',V0',CD,CR及加速度,即
式中,CD和CR作为常值估计。
根据待预报地区,粗略估计卫星过境时间t,利用轨道预报计算公式,结合卫星过境前24小时的长时段轨道数据作为初始值,对过境的轨道参数进行预报。由于该预报一个渐进递推的长时段处理,过程涉及的数据较多,这里仅给出一例作为参考,24小时前的轨道数据初值为:P0'=[-728959.1 6827631.7 473238.5]T和V0'=[1543.3 -367.0 7528.4]T。结合其他可获取的参数,可计算得到过境预报区域的卫星的激光质心在地球固定坐标系的坐标值P0=[-1855244.6 4669501.6 4693461.4]T和位移速度值V0=[-287.4 5397.1 -5468.8]T。
步骤3.2,根据下述公式计算轨道坐标系相对于地球固定坐标系的旋转矩阵MT_o;
其中,P0=[X0 Y0 Z0]T为卫星激光质心在地球固定坐标系的坐标值,V0=[VX0 VY0VZ0]T为卫星激光质心在地球固定坐标系的位移速度值,A1,A2,A3,B1,B2,B3,C1,C2,C3为旋转矩阵MT_o的各个元素。
将步骤3.1得到的过境预报区域的卫星的激光质心在地球固定坐标系的坐标值P和位移速度值V代入到步骤3.2的公式,可以求得轨道坐标系相对于地球固定坐标系的旋转矩阵
步骤3.3,根据下述公式将激光测高仪的指向由轨道坐标系转换到地球固定坐标系;
MT_l=MT_o*Mo_l
式中,Mo_l表示激光测高仪在轨道坐标系下的指向,MT_o表示轨道坐标系相对于地球固定坐标系的旋转矩阵,MT_l表示激光测高仪在地球固定坐标系下的指向。
将步骤2.3和步骤3.2计算得到的激光测高仪的指向由本体坐标系转换到轨道坐标系的旋转矩阵与轨道坐标系相对于地球固定坐标系的旋转矩阵相乘,可以得到激光测高仪的指向由本体坐标系转换到地球固定坐标系的向量
步骤4,通过计算地球固定坐标系下的激光测高仪的指向与地球椭球的交点,确定星载激光测高仪的足印位置。
根据光束的物理传播路径,在激光到达地面后,需要求地球固定坐标系下的激光测高仪的指向与地球椭球的交点,所述交点即本发明要确定的星载激光测高仪的足印位置,即激光探测器的布设位置的坐标。步骤4具体包括以下子步骤:
步骤4.1,根据下述公式计算激光指向与地球椭球相交的地球固定坐标系下的地面点坐标(X1,Y1,Z1):
其中,MT=[X1 Y1 Z1]T,X1、Y1和Z1为待求的激光指向与地球椭球相交的地球固定坐标系下的地面点坐标,μ为待求的比例系数,MT_l为激光测高仪在地球固定坐标系下的指向,P0为卫星激光质心在地球固定坐标系的坐标值,a和b分别为地球椭球长半轴和短半轴长度(a=6378137.0米,b=6356752.3米),h为地面大地高程,可以通过其他辅助高程数据获取,初始值给0。
其中,MT=P0+μ*MT_l可以表示为:或表示为将X1、Y1、Z1代入到公式中,可以得到下述方程
式中,X0、Y0和Z0为步骤3.1计算得到卫星激光质心在地球固定坐标系的坐标值,MT_l为步骤3.3计算得到的值。
解上述一元二次方程,得到比例系数μ,因为方程会有两个解,最终选择绝对值小的u=506437.3为最终解,根据MT=P0+μ*MT_l从而可以计算出MT=[-1718742.3 4325848.34347414.8]T。
步骤4.2,根据下述公式将地球固定坐标系下的地面点坐标(X1,Y1,Z1)转换成大地经纬度坐标系下的地面点坐标(L,B,H);
(X1,Y1,Z1)对于我们来说,没法直接使用,需要转换为我们经常用的经纬度格式(L,B,H),这个也是从天上到地面的必要转换。
其中,e1和e2分别为第一偏心率和第二偏心率,a和b分别为地球椭球长半轴和短半轴长度,X1、Y1、Z1为通过步骤4.1计算出的激光指向与地球椭球相交的地球固定坐标系下的地面点坐标。
根据上述公式,可以求得地面坐标(L,B,H)=(111.66887,43.23643,1079.99)。
根据本发明的优选实施例,为了得到更精确的星载激光测高仪的足印位置,步骤4还可以包括:
步骤4.3,通过迭代计算,获得精确的激光探测器的布设位置;步骤4.3具体包括:将步骤4.2计算得到的地面点坐标(L,B)代入到高程数据中,重新内插得到更精确的高程值H',将计算得到的H'代入到步骤4.1和步骤4.2进行迭代运算,以得到精确的星载激光测高仪的足印位置。
具体来说,将步骤4.2得到的地面点坐标(L,B)代入到高程数据中,重新内插得到更精确的高程值H',内插公式采用双线性内插:
H'=(1-Δx)(1-Δy)H1+(1-Δx)ΔyH2+Δx(1-Δy)H3+ΔxΔyH4
式中,H1、H2、H3和H4为H'附近x、y方向上相邻点的高程值,Δx为H'距H3和H4的x方向距离差,Δy为H'距H1和H2的y方向距离差。
将计算得到的H'代入到步骤4.1和步骤4.2,重复计算,此过程一般迭代三次,最终得到精确的星载激光测高仪的足印位置(L,B,H),即,足印中心点的坐标。
将步骤4.2计算得到的地面点坐标(111.66887,43.23643)代入到我们使用的AW3D30数据(公开数据)中,采用上述内插公式,可以内插得到新的高程值H=951.21;将H=951.21代入到步骤4.1和步骤4.2得到新的地面点坐标(111.66879,43.23638,951.21),将(111.66879,43.23638)代入到AW3D30数据中,可以内插得到新的高程值H=950.14;将H=950.14代入到步骤4.1和步骤4.2得到新的地面点坐标(111.66879,43.23638,950.08),该位置即为待铺设探测器阵的中心点。
本发明结合了传统的光学相机成像的原理,并充分顾及了激光测高仪自身的测高机理,预报结果能够较好地符合卫星真实的飞行过境情况,具有较高的预报精度。
以上内容仅为本发明的较佳实施例,对于本领域的普通技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,本说明书内容不应理解为对本发明的限制。
Claims (9)
1.一种预先确定星载激光测高仪的足印位置的方法,其特征在于,所述方法包括以下步骤:
步骤1,计算激光测高仪在卫星本体坐标系中的指向Mb_l,步骤1具体包括以下子步骤:
步骤1.1,预估激光测高仪与平台之间的安装角(αr,βr);
步骤1.2,基于步骤1.1预估的安装角,计算激光测高仪在卫星本体坐标系中的指向Mb_l;
步骤2,将激光测高仪的指向从卫星本体坐标系旋转至轨道坐标系Mo_l;
步骤3,将激光测高仪的指向从轨道坐标系旋转至地球固定坐标系,获得激光测高仪在地球固定坐标系下的指向MT_l;
步骤4,通过计算地球固定坐标系下的激光测高仪的指向与地球椭球的交点,确定星载激光测高仪的足印位置,其中,所述步骤4具体包括以下子步骤:
步骤4.1,根据下述公式计算激光指向与地球椭球相交的地球固定坐标系下的地面点坐标(X1,Y1,Z1):
其中,MT=[X1 Y1 Z1]T,X1、Y1和Z1为待求的激光指向与地球椭球相交的地球固定坐标系下的地面点坐标,μ为待求的比例系数,MT_l为激光测高仪在地球固定坐标系下的指向,P0为卫星激光质心在地球固定坐标系的坐标值,a和b分别为地球椭球长半轴和短半轴长度,h为地面大地高程;
步骤4.2,根据下述公式将地球固定坐标系下的地面点坐标(X1,Y1,Z1)转换成大地经纬度坐标系下的地面点坐标(L,B,H):
其中,e1和e2分别为第一偏心率和第二偏心率,a和b分别为地球椭球长半轴和短半轴长度,X1、Y1、Z1为通过步骤4.1计算出的激光指向与地球椭球相交的地球固定坐标系下的地面点坐标。
2.根据权利要求1所述的方法,其特征在于,所述步骤4还包括:
步骤4.3,通过迭代计算,获得精确的星载激光足印位置。
3.根据权利要求2所述的方法,其特征在于,所述步骤4.3具体包括:将步骤4.2计算得到的地面点坐标(L,B)代入到高程数据中,重新内插得到更精确的高程值H',将计算得到的H'代入到步骤4.1和步骤4.2进行迭代运算,以得到精确的星载激光足印位置。
4.根据权利要求1所述的方法,其特征在于,在所述步骤1.1中,根据下式预估激光测高仪与平台之间的安装角:
其中,(αr,βr)为预估的安装角,表示从1层至3层逐层计算,αi、βi表示第i层遍历的安装角,表示括号内第i层取最小值时对应的αi、βi,并将其赋值给α0、β0,q(i)为第i层遍历间隔,p(i)为第i层的区域间隔,表示对括号内进行遍历,行遍历序号m,列遍历序号n,E(αi,βi)为αi、βi指向下计算某轨的高程值与已知高程值之间的标准差。
5.根据权利要求1所述的方法,其特征在于,在所述步骤1.2中,根据下述公式计算激光测高仪在卫星本体坐标系中的指向:
式中,Mb_l为步骤1待求的激光测高仪在卫星本体坐标系中的指向,αr、βr为步骤1.1求得的激光测高仪的安装角。
6.根据权利要求1所述的方法,其特征在于,所述步骤2具体包括以下子步骤:
步骤2.1,估计卫星平台的本体坐标系相对于轨道坐标系的姿态值其中,roll、pitch、yaw分别为卫星本体坐标系相对于轨道坐标系的三轴欧拉角(横滚角、俯仰角和偏航角);
步骤2.2,根据下述公式计算本体坐标系相对于轨道坐标系的旋转矩阵Mo_b:
a1=cos(pitch)*cos(yaw)
a2=-cos(pitch)*sin(yaw)
a3=sin(pitch)
b1=-sin(roll)*sin(pitch)*cos(yaw)+cos(roll)*sin(yaw)
b2=sin(roll)*sin(pitch)*sin(yaw)+cos(roll)*cos(yaw)
b3=sin(roll)*cos(pitch)
c1=-cos(roll)*sin(pitch)*cos(yaw)-sin(roll)*sin(yaw)
c2=cos(roll)*sin(pitch)*sin(yaw)-sin(roll)*cos(yaw)
c3=cos(roll)*cos(pitch)
步骤2.3,根据下述公式将激光测高仪的指向由本体坐标系转换到轨道坐标系;
Mo_l=Mo_b*Mb_l
其中,Mo_l为待求的激光测高仪在轨道坐标系下的指向,Mo_b为步骤2.2计算的本体坐标系相对于轨道坐标系的旋转矩阵,Mb_l为步骤1.2计算的激光测高仪的指向。
7.根据权利要求1所述的方法,其特征在于,所述步骤3具体包括以下子步骤:
步骤3.1,计算卫星激光质心在地球固定坐标系的坐标值P0=[X0 Y0 Z0]T和位移速度值V0=[VX0 VY0 VZ0]T;
步骤3.2,根据下述公式计算轨道坐标系相对于地球固定坐标系的旋转矩阵MT_o:
其中,其中,P0=[X0 Y0 Z0]T为卫星激光质心在地球固定坐标系的坐标值,V0=[VX0 VY0VZ0]T为卫星激光质心在地球固定坐标系的位移速度值。
步骤3.3,根据下述公式将激光测高仪的指向由轨道坐标系转换到地球固定坐标系;
MT_l=MT_o*Mo_l
式中,Mo_l表示激光测高仪在轨道坐标系下的指向,MT_o表示轨道坐标系相对于地球固定坐标系的旋转矩阵,MT_l表示激光测高仪在地球固定坐标系下的指向。
8.根据权利要求1-7中的任意一项所述的方法,其特征在于,迭代次数优选的为3次。
9.一种计算机程序,其中,当执行所述计算机程序时,执行前述任一权利要求所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611253206.0A CN106643804B (zh) | 2016-12-30 | 2016-12-30 | 一种预先确定星载激光测高仪的足印位置的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611253206.0A CN106643804B (zh) | 2016-12-30 | 2016-12-30 | 一种预先确定星载激光测高仪的足印位置的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106643804A true CN106643804A (zh) | 2017-05-10 |
CN106643804B CN106643804B (zh) | 2019-11-22 |
Family
ID=58836948
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611253206.0A Active CN106643804B (zh) | 2016-12-30 | 2016-12-30 | 一种预先确定星载激光测高仪的足印位置的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106643804B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108414998A (zh) * | 2018-03-02 | 2018-08-17 | 国家测绘地理信息局卫星测绘应用中心 | 一种卫星激光测高仪回波波形模拟仿真方法及设备 |
CN108828537A (zh) * | 2018-04-04 | 2018-11-16 | 南京理工大学 | 一种激光测高仪综合测试***及其综合测试方法 |
CN110006448A (zh) * | 2019-04-09 | 2019-07-12 | 武汉大学 | 星载激光测高仪指向角***误差在轨标定精度的评估方法 |
CN111505608A (zh) * | 2020-05-06 | 2020-08-07 | 自然资源部国土卫星遥感应用中心 | 一种基于星载激光单片足印影像的激光指向在轨检校方法 |
CN111623802A (zh) * | 2020-07-01 | 2020-09-04 | 中国科学院上海技术物理研究所 | 一种基于足印影像的星载激光测高仪在轨指向检校方法 |
CN113640848A (zh) * | 2021-08-05 | 2021-11-12 | 窦显辉 | 无人机的地面激光足印数据采集方法、***、介质及设备 |
CN116840851A (zh) * | 2023-07-05 | 2023-10-03 | 中国科学院空天信息创新研究院 | 一种星载对地激光测高仪地面探测器布设方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104296729A (zh) * | 2014-09-28 | 2015-01-21 | 中国科学院长春光学精密机械与物理研究所 | 两线阵空间测绘相机与激光高度计共光路*** |
CN104931022A (zh) * | 2015-04-21 | 2015-09-23 | 国家测绘地理信息局卫星测绘应用中心 | 基于星载激光测高数据的卫星影像立体区域网平差方法 |
CN105550639A (zh) * | 2015-12-07 | 2016-05-04 | 国家测绘地理信息局卫星测绘应用中心 | 对地观测激光测高卫星高程控制点自动提取方法和数据处理方法 |
CN105631886A (zh) * | 2015-12-01 | 2016-06-01 | 中国科学院上海技术物理研究所 | 一种基于航空影像的激光光斑与足印相机相对定位方法 |
-
2016
- 2016-12-30 CN CN201611253206.0A patent/CN106643804B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104296729A (zh) * | 2014-09-28 | 2015-01-21 | 中国科学院长春光学精密机械与物理研究所 | 两线阵空间测绘相机与激光高度计共光路*** |
CN104931022A (zh) * | 2015-04-21 | 2015-09-23 | 国家测绘地理信息局卫星测绘应用中心 | 基于星载激光测高数据的卫星影像立体区域网平差方法 |
CN105631886A (zh) * | 2015-12-01 | 2016-06-01 | 中国科学院上海技术物理研究所 | 一种基于航空影像的激光光斑与足印相机相对定位方法 |
CN105550639A (zh) * | 2015-12-07 | 2016-05-04 | 国家测绘地理信息局卫星测绘应用中心 | 对地观测激光测高卫星高程控制点自动提取方法和数据处理方法 |
Non-Patent Citations (3)
Title |
---|
MAGRUDER L A ,ET AL: "Laser pointing angle and time of measurement verification of the ICESat laser altimeter using a ground-based electro-optical detection system", 《JOURNAL OF GEODESY》 * |
唐新明等: "卫星激光测高严密几何模型构建及精度初步验证", 《测绘学报》 * |
赵欣等: "星载激光测高仪大气传输延迟对测距精度的影响", 《红外与激光工程》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108414998A (zh) * | 2018-03-02 | 2018-08-17 | 国家测绘地理信息局卫星测绘应用中心 | 一种卫星激光测高仪回波波形模拟仿真方法及设备 |
CN108414998B (zh) * | 2018-03-02 | 2021-10-15 | 国家测绘地理信息局卫星测绘应用中心 | 一种卫星激光测高仪回波波形模拟仿真方法及设备 |
CN108828537A (zh) * | 2018-04-04 | 2018-11-16 | 南京理工大学 | 一种激光测高仪综合测试***及其综合测试方法 |
CN110006448B (zh) * | 2019-04-09 | 2021-09-14 | 武汉大学 | 星载激光测高仪指向角***误差在轨标定精度的评估方法 |
CN110006448A (zh) * | 2019-04-09 | 2019-07-12 | 武汉大学 | 星载激光测高仪指向角***误差在轨标定精度的评估方法 |
CN111505608A (zh) * | 2020-05-06 | 2020-08-07 | 自然资源部国土卫星遥感应用中心 | 一种基于星载激光单片足印影像的激光指向在轨检校方法 |
CN111505608B (zh) * | 2020-05-06 | 2021-03-12 | 自然资源部国土卫星遥感应用中心 | 一种基于星载激光单片足印影像的激光指向在轨检校方法 |
CN111623802A (zh) * | 2020-07-01 | 2020-09-04 | 中国科学院上海技术物理研究所 | 一种基于足印影像的星载激光测高仪在轨指向检校方法 |
CN111623802B (zh) * | 2020-07-01 | 2023-03-28 | 中国科学院上海技术物理研究所 | 一种基于足印影像的星载激光测高仪在轨指向检校方法 |
CN113640848A (zh) * | 2021-08-05 | 2021-11-12 | 窦显辉 | 无人机的地面激光足印数据采集方法、***、介质及设备 |
CN113640848B (zh) * | 2021-08-05 | 2024-02-13 | 窦显辉 | 无人机的地面激光足印数据采集方法、***、介质及设备 |
CN116840851A (zh) * | 2023-07-05 | 2023-10-03 | 中国科学院空天信息创新研究院 | 一种星载对地激光测高仪地面探测器布设方法 |
CN116840851B (zh) * | 2023-07-05 | 2024-01-12 | 中国科学院空天信息创新研究院 | 一种星载对地激光测高仪地面探测器布设方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106643804B (zh) | 2019-11-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106643804B (zh) | 一种预先确定星载激光测高仪的足印位置的方法 | |
CN110646782B (zh) | 一种基于波形匹配的星载激光在轨指向检校方法 | |
CN106643709B (zh) | 一种海上运载体的组合导航方法及装置 | |
CN106871932A (zh) | 基于金字塔搜索地形匹配的星载激光在轨指向检校方法 | |
CN103674034B (zh) | 多波束测速测距修正的鲁棒导航方法 | |
CN103744052B (zh) | 一种应用于空中目标定位的双星测时差测向方法及装置 | |
Pittet et al. | Spin motion determination of the Envisat satellite through laser ranging measurements from a single pass measured by a single station | |
CN104655152A (zh) | 一种基于联邦滤波的机载分布式pos实时传递对准方法 | |
CN105241456B (zh) | 巡飞弹高精度组合导航方法 | |
CN102116628A (zh) | 一种着陆或附着深空天体探测器的高精度导航方法 | |
CN107300700B (zh) | 敏捷合成孔径雷达卫星聚束模式姿态机动需求计算方法 | |
CN109708663B (zh) | 基于空天飞机sins辅助的星敏感器在线标定方法 | |
Lu et al. | Dual-filter transfer alignment for airborne distributed POS based on PVAM | |
CN111505608B (zh) | 一种基于星载激光单片足印影像的激光指向在轨检校方法 | |
CN103017772A (zh) | 一种基于可观性分析的光学和脉冲星融合自主导航方法 | |
CN103968844A (zh) | 基于低轨平台跟踪测量的大椭圆机动航天器自主导航方法 | |
Qin et al. | An innovative navigation scheme of powered descent phase for Mars pinpoint landing | |
CN104765373B (zh) | 一种星上相对运动状态获取方法 | |
CN107449392A (zh) | 基于惯性测量单元的边坡变形测量方法 | |
Li et al. | Opportunity rover localization and topographic mapping at the landing site of Meridiani Planum, Mars | |
Adams et al. | Passive optical terrain relative navigation using APLNav | |
US20220065587A1 (en) | System and method of hypersonic object tracking | |
Huang et al. | Research on UAV flight performance test method based on dual antenna GPS/ins integrated system | |
CN116698048A (zh) | 一种基于脉冲星/星间测距/陆标的组合导航方法 | |
CN103245948A (zh) | 双区成像合成孔径雷达图像匹配导航的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 100048 No. 28 Lianhuachi West Road, Haidian District, Beijing Applicant after: Ministry of Natural Resources Land Satellite Remote Sensing Application Center Address before: 100048 No. 28 Lianhuachi West Road, Haidian District, Beijing Applicant before: Satellite Surveying and Mapping Application Center, NASG |
|
GR01 | Patent grant | ||
GR01 | Patent grant |