CN105115477B - 对地面点目标推扫成像任务参数的星载求解方法 - Google Patents
对地面点目标推扫成像任务参数的星载求解方法 Download PDFInfo
- Publication number
- CN105115477B CN105115477B CN201510447888.8A CN201510447888A CN105115477B CN 105115477 B CN105115477 B CN 105115477B CN 201510447888 A CN201510447888 A CN 201510447888A CN 105115477 B CN105115477 B CN 105115477B
- Authority
- CN
- China
- Prior art keywords
- satellite
- represent
- coordinate system
- ground
- target
- 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
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Navigation (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供了一种对地面点目标推扫成像任务参数的星载求解方法,包括:步骤1:建立坐标系,并定义卫星、地面目标、地心三者间的矢量;步骤2:建立包含摆镜的光学系统等效光路;步骤3:获取卫星星下点太阳高度角信息、卫星位置以及地面目标信息;步骤4:求解出卫星和地面目标在瞬时惯性坐标系下的运动状态表征参数;步骤5:计算出地面目标在瞬时惯性坐标系下的坐标和变化率随时间变化的表达式;步骤6:计算出地面目标在卫星本体坐标系下的坐标和变化率随时间变化的表达式;步骤7:建立约束方程,并求解出成像任务参数。本发明在获取地面点目标坐标后,在轨卫星能够自主、实时、快速、精准的调整摆镜指向目标,以确保可对目标成像。
Description
技术领域
本发明涉及航天空间遥感成像任务技术领域,具体地,涉及一种对地面点目标推扫成像任务参数的星载求解方法。
背景技术
推扫成像相机采用线阵延时积分作为接收器,通过多次曝光解决传统面阵相机通光量不足的问题。由于相机视场角的限制,要对地面特定目标成像,根据星地相对关系,往往需要通过姿态或相机摆镜,实现对相机指向的调整。使得调整后的视场,随着卫星飞行,将在某时刻恰好能够拍摄到特定地面目标。
传统的对地遥感卫星要完成对特定区域成像,通常是通过星地联合控制完成的,在地面预先获取目标位置信息,并由地面解算成像任务参数后,如卫星姿态、成像载荷指向调整、成像时刻等,生成程控作业,在卫星过境时上注星上。由于星地链路无法做到全时互通,因此传统做法无法应对新发现的目标。因此一种星上实时自主获知地面目标坐标,并能够在轨自主解算成像任务参数的卫星需求应运而生,本发明为这种需求提供了一种星载解决方法。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种对地面点目标推扫成像任务参数的星载求解方法。
步骤1:建立坐标系,所述坐标系包括瞬时惯性坐标系、卫星轨道坐标系、卫星本体坐标系、虚拟相机坐标系、虚像平面坐标系,并定义卫星、地面目标、地心三者间的矢量;
步骤2:利用建立的坐标系辅助建立包含摆镜的光学系统等效光路;
步骤3:获取卫星星下点太阳高度角,并判断满足成像的光照条件后,在包括卫星位置、卫星速度、地面目标位置以及系统时间在内的WGS84坐标系下获取卫星和地面目标的各项信息;
步骤4:利用卫星星下点太阳高度角信息,WGS84坐标系下卫星和地面目标的各项信息求解出卫星和地面目标在瞬时惯性坐标系下的运动状态表征参数,所述运动状态表征参数包括:卫星位置、卫星速度、卫星轨道角动量、卫星的倾角、升交点赤经、轨道角速度、实时轨道幅角、卫星轨道坐标系的转换矩阵及变化率以及卫星轨道坐标系到卫星本体坐标系的转换矩阵;
步骤5:根据卫星和地面目标在WGS84坐标系下的运动状态表征参数计算出地面目标在瞬时惯性坐标系下的坐标和变化率随时间变化的表达式;
步骤6:根据卫星和地面目标在瞬时惯性坐标系下的运动状态表征参数、地面目标坐标随时间变化参数以及地面目标变化率随时间变化参数计算出地面目标在卫星本体坐标系下的坐标和变化率随时间变化的表达式;
步骤7:根据卫星本体坐标系下的坐标和变化率随时间变化参数建立约束方程,并求解出成像任务参数,所述成像任务参数包括:相机指向调整参数和成像开始时刻。
优选地,所述步骤1包括:
步骤1.1:建立坐标系,所述坐标系包括瞬时惯性坐标系、卫星轨道坐标系、卫星本体坐标系、虚拟相机坐标系、虚像平面坐标系;
-瞬时惯性坐标系Oexiyizi,Oexi轴在地球赤道平面内,指向GPS时刻的WGS84坐标系X轴方向;Oezi轴垂直于地球赤道平面,与地球自转角速度矢量方向一致;Oeyi与Oexi、Oezi轴垂直,且构成右手坐标系;所述坐标系简称I系;
-卫星轨道坐标系Osxoyozo,Os表示原点,即卫星质心;Osxo轴在卫星轨道平面内,指向卫星的运动方向;Oszo轴指向地心;Osyo与Osxo、Oszo轴垂直,且构成右手坐标系;所述坐标系简称O系;
-卫星本体坐标系Osxbybzb,Os表示原点,即卫星质心;Osxb轴沿卫星飞行方向,Oszb轴与Osyb轴、Osxb轴构成右手坐标系;当卫星处于标称零姿态运动时,Osxb、Osyb、Oszb轴分别与Osxo、Osyo、Oszo轴平行,且方向相同;所述坐标系简称b系;
-虚拟相机坐标系Ocxcyczc,Oc表示原点,即摆镜虚相机的摄像中心;Ocxcyc平面为物镜平面;Oczc轴沿光轴方向,摆镜无摆角时Ocxc、Ocyc轴与轨道系重合,构成右手坐标系;所述坐标系简称c系;
-虚像平面坐标系Opxpyp,Op表示原点,即相机像平面光轴中心处;相机像平面法线与Oczc轴平行,且Opxp、Opyp轴分别与Ocxc、Ocyc平行,且方向相同;所述坐标系简称p系;
步骤1.2:定义卫星、地面目标、地心三者间的矢量;
-Ret用于描述地面目标,表示从I系原点指向地面目标的矢量;
-Res用于描述卫星位置,表示从I系原点指向卫星质心的矢量;
-Rst用于描述地面目标相对卫星位置,表示从卫星质心指向地面目标的矢量;
-Rct用于描述地面目标相对相机位置,表示从c系原点指向地面目标的矢量;
-表示在I系中的位置矢量Ret;
-表示在轨道系O系中的位置矢量Res。
优选地,所述步骤2包括:安装在卫星对地面上的相机光轴经过摆镜的反射后照向地面目标位置;当摆镜偏转角时,相当于相机偏转了角,则通过摆镜就能够实现对地面目标的推扫。
优选地,所述步骤3包括:获取卫星星下点太阳高度角信息,并在WGS84坐标系下获取卫星和地面目标的各项信息,包括:卫星位置、卫星速度、地面目标位置以及系统时间;此时瞬时惯性坐标系与对应的WGS84坐标系重合,每一个GPS测量节拍,都以当时WGS84坐标系的指向确定重合的瞬时惯性坐标系。
优选地,所述步骤4包括:
步骤4.1:判断星下点太阳高度角是否满足成像的光照条件,若满足则进入步骤4.2,若不满足则结束所述对地面点目标推扫成像任务参数的星载求解方法流程;
步骤4.2:利用卫星上的GPS测得的当前时刻的卫星位置和速度,结合地球自转角速度,将短时间的卫星运动简化为等速圆周运动,则推导出瞬时惯性坐标系、卫星轨道坐标系、卫星本体坐标系的转换关系;
步骤4.2.1:GPS测得的当前时刻在WGS84坐标系下的卫星位置和速度,计算公式如下:
式中:RWGS84表示GPS测得的当前时刻在WGS84坐标系下的卫星位置矢量,XRGPS表示RWGS84的X轴分量,YRGPS表示RWGS84的Y轴分量,ZRGPS表示RWGS84的Z轴分量;VWGS84表示GPS测得的当前时刻在WGS84坐标系下的卫星速度矢量,XVGPS表示VWGS84X轴分量,YVGPS表示VWGS84Y轴分量,ZVGPS表示VWGS84Z轴分量;
步骤4.2.2:根据WGS84坐标系下的卫星位置和速度求解在I系下,同时考虑地球自转引起的牵连速度时的卫星位置和速度,计算公式如下:
式中:TGPS表示卫星位置速度对应的时刻,分别表示I系下卫星位置和速度矢量,ωe表示地球自转角速率;
步骤4.2.3:根据I系下的卫星位置和速度求解卫星轨道角动量,计算公式如下:
式中:H表示卫星轨道角动量,HX表示卫星轨道角动量X轴分量,HY表示卫星轨道角动量Y轴分量,HZ表示卫星轨道角动量Z轴分量;
步骤4.2.4:根据卫星轨道角动量求解I系下卫星的倾角、升交点赤经、轨道角速度,计算公式如下:
xi=[1 0 0]T,zi=xi×N
zi(3)≥0时,Ω=arccos(N·xi),zi(3)<0时,Ω=-arccos(N·xi)
其中:
式中:i表示I系下卫星的轨道倾角、Ω表示升交点赤经、ω表示轨道角速度:xi表示惯性系X轴矢量,zi表示惯性系Z轴矢量,zi(3)表示zi在Z轴分量,N表示轨道节线单位矢量;
步骤4.2.5:根据I系下卫星的轨道角动量、倾角、轨道角速度求解TGPS时刻的轨道幅角u0,轨道幅角u,计算公式如下:
ZRGPS≥0时,ZRGPS<0时,
式中,t表示从起始时刻开始的时间差;
步骤4.2.6:根据I系下TGPS时刻的轨道幅角u0,轨道幅角u求解I系到卫星轨道坐标系的转换矩阵及变化率,计算公式如下:
式中:Aoi表示惯性系到轨道系转移矩阵,表示惯性系到轨道系转移矩阵变化率,ωb表示轨道角速率在本体系分量,ω表示轨道角速度;
步骤4.2.7:求解卫星轨道坐标系到卫星本体坐标系的转换矩阵,计算公式如下:
式中:Abo表示轨道系到问题系转移矩阵,θ表示俯仰姿态角,表示滚动自他角,ψ表示偏航姿态角。
优选地,所述步骤5包括:根据卫星和地面目标在在WGS84坐标系下的运动状态表征参数计算出地面目标在瞬时惯性坐标系下的坐标和变化率随时间变化的表达式,计算公式如下:
式中:表示目标初始位置,表示目标位置矢量惯性系分量,表示目标位置矢量变化率,ωe表示地球自转角速率。
优选地,所述步骤6包括:根据卫星和地面目标在瞬时惯性坐标系下的运动状态表征参数、地面目标坐标随时间变化参数以及地面目标变化率随时间变化参数计算出地面目标在卫星本体坐标系下的坐标和变化率随时间变化的表达式;计算公式如下:
式中:表示卫星到目标矢量在卫星本体系分量,X(t)表示X轴分量,Y(t)表示Y轴分量,Z(t)表示Z轴分量;表示卫星到目标矢量变化率在卫星本体系分量,X’(t)表示X轴分量,Y’(t)表示Y轴分量,Z’(t)表示Z轴分量。
优选地,所述步骤7包括:根据卫星本体坐标系下的坐标和变化率随时间变化参数建立约束方程,并求解出成像任务参数;
步骤7.1:假设相机后视角为θ1,建立约束方程如下:
步骤7.2:将约束方程转换为F(t)=0的函数,函数表达式如下:
步骤7.3:求解函数F(t)的零点;
步骤7.3.1:对函数F(t)进行求导运算,得到如下的方程:
步骤7.3.2:利用牛顿迭代求解法计算时间变量t,计算公式如下:
步骤7.3.3:令t0的初始值为0,则t秒之后能够得到c系下地面目标相对卫星位置矢量Rst c;
步骤7.4:求解卫星成像任务参数,包括:摆镜摆角和相机成像时间TShoot;
若需要提前n秒成像,则相机成像时间为:
TShoot=TGPS+t-n,
相机光轴需绕滚动轴调整量为:
若反映到摆镜摆角上,则摆镜摆角为:
步骤7.5:判断是否满足成像条件,若满足则进入步骤7.6,若不满足则结束所述对地面点目标推扫成像任务参数的星载求解方法流程,其中所述成像条件为:在摆角最大摆动范围内,且t大于星上成像前的最小准备时间;
步骤7.6:星上准备好后,提前驱动相机的摆镜,并在TShoot时刻开始成像。
与现有技术相比,本发明具有如下的有益效果:
1、本发明提供的方法是一种星载解决方法,在获取地面目标后,在轨自主解算成像任务参数,使得成像更加精准。
2、本发明提供的方法通过求解星上相机指向调整参数、成像开始时刻等任务执行参数,可应对在轨实时获得地面目标坐标后的快速、自主成像任务。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明提供的多目标连续成像原理示意图;
图2为本发明提供的相机像平面坐标系与虚拟相机坐标系关系示意图;
图3为本发明提供的瞬时惯性坐标系、卫星轨道坐标系、卫星本体坐标系的矢量关系示意图;
图4为本发明提供的含摆镜光学系统光路示意图;
图5为本发明提供的镜面反射后的虚拟光轴示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进。这些都属于本发明的保护范围。
根据本发明提供的对地面点目标推扫成像任务参数的星载求解方法,包括如下步骤:
步骤1:建立坐标系,所述坐标系包括瞬时惯性坐标系、卫星轨道坐标系、卫星本体坐标系、虚拟相机坐标系、虚像平面坐标系,并定义卫星、地面目标、地心三者间的矢量;
步骤2:利用建立的坐标系辅助建立包含摆镜的光学系统等效光路;
步骤3:获取卫星星下点太阳高度角,并判断满足成像的光照条件后,在包括卫星位置、卫星速度、地面目标位置以及系统时间在内的WGS84坐标系下获取卫星和地面目标的各项信息;
步骤4:利用卫星星下点太阳高度角信息,WGS84坐标系下卫星和地面目标的各项信息求解出卫星和地面目标在瞬时惯性坐标系下的运动状态表征参数,所述运动状态表征参数包括:卫星位置、卫星速度、卫星轨道角动量、卫星的倾角、升交点赤经、轨道角速度、实时轨道幅角、卫星轨道坐标系的转换矩阵及变化率以及卫星轨道坐标系到卫星本体坐标系的转换矩阵;
步骤5:根据卫星和地面目标在WGS84坐标系下的运动状态表征参数计算出地面目标在瞬时惯性坐标系下的坐标和变化率随时间变化的表达式;
步骤6:根据卫星和地面目标在瞬时惯性坐标系下的运动状态表征参数、地面目标坐标随时间变化参数以及地面目标变化率随时间变化参数计算出地面目标在卫星本体坐标系下的坐标和变化率随时间变化的表达式;
步骤7:根据卫星本体坐标系下的坐标和变化率随时间变化参数建立约束方程,并求解出成像任务参数,所述成像任务参数包括:相机指向调整参数和成像开始时刻。
具体地,如图1中所示,对于线阵推扫相机,其推扫方向一般与卫星飞行方向相同,视场的扩展通常仅依靠绕卫星滚动轴的姿态偏置或依靠滚动轴的摆镜摆动完成。OXYZ为卫星轨道坐标系,标称状态下,卫星三轴姿态为0,相机对地视场较窄,视轴正对地。在t0时刻,获知地面目标Target1的坐标,星上综合可成像光照约束、卫星当前位置速度参数、目标坐标、成像前准备时间多少等,解算出在Δt时间后,侧摆角为时,相机推扫阵列在地面投影的推扫线中心,恰好遇到地面点目标。
优选地,所述步骤1包括:
步骤1.1:建立坐标系,所述坐标系包括瞬时惯性坐标系、卫星轨道坐标系、卫星本体坐标系、虚拟相机坐标系、虚像平面坐标系;
-瞬时惯性坐标系Oexiyizi,Oexi轴在地球赤道平面内,指向GPS时刻的WGS84坐标系X轴方向;Oezi轴垂直于地球赤道平面,与地球自转角速度矢量方向一致;Oeyi与Oexi、Oezi轴垂直,且构成右手坐标系;所述坐标系简称I系;
-卫星轨道坐标系Osxoyozo,Os表示原点,即卫星质心;Osxo轴在卫星轨道平面内,指向卫星的运动方向;Oszo轴指向地心;Osyo与Osxo、Oszo轴垂直,且构成右手坐标系;所述坐标系简称O系;
-卫星本体坐标系Osxbybzb,Os表示原点,即卫星质心;Osxb轴沿卫星飞行方向,Oszb轴与Osyb轴、Osxb轴构成右手坐标系;当卫星处于标称零姿态运动时,Osxb、Osyb、Oszb轴分别与Osxo、Osyo、Oszo轴平行,且方向相同;所述坐标系简称b系;
-虚拟相机坐标系Ocxcyczc,Oc表示原点,即摆镜虚相机的摄像中心;Ocxcyc平面为物镜平面;Oczc轴沿光轴方向,摆镜无摆角时Ocxc、Ocyc轴与轨道系重合,构成右手坐标系;所述坐标系简称c系;
-虚像平面坐标系Opxpyp,Op表示原点,即相机像平面光轴中心处;相机像平面法线与Oczc轴平行,且Opxp、Opyp轴分别与Ocxc、Ocyc平行,且方向相同;所述坐标系简称p系;
具体地,相机像平面坐标系与虚拟相机坐标系的关系如图2所示。
步骤1.2:定义卫星、地面目标、地心三者间的矢量;
-Ret用于描述地面目标,表示从I系原点指向地面目标的矢量;
-Res用于描述卫星位置,表示从I系原点指向卫星质心的矢量;
-Rst用于描述地面目标相对卫星位置,表示从卫星质心指向地面目标的矢量;
-Rct用于描述地面目标相对相机位置,表示从c系原点指向地面目标的矢量;
-表示在I系中的位置矢量Ret;
-表示在轨道系O系中的位置矢量Res。
具体地,地面目标T、卫星质心Os与摆镜虚相机的摄像中心Oc在空间中的位置与矢量关系如图3所示。
具体地,在步骤1中定义了一种有利于解决问题的瞬时惯性坐标系,在GPS给出卫星位置和速度的测量时刻,瞬时惯性坐标系与对应的WGS84坐标系重合。每一个GPS测量节拍,都对以当时的WGS84坐标系空间的指向确定重合的瞬时惯性坐标系。在当节拍内,瞬时惯性坐标系系与WGS84坐标系的区别仅为WGS84坐标系相对瞬时惯性坐标系系有一个绕Z轴的地球自转角速度。
优选地,所述步骤2包括:当卫星处于标称零姿态时,即卫星本体坐标系OsZb指向地心,OsXb为飞行正方向,OsYb指向轨道面法相;安装在卫星对地面上的相机光轴经过摆镜的反射后照向地面目标位置;当摆镜偏转角时,相当于相机偏转了角,则通过摆镜就能够实现对地面目标的推扫。
具体地,如图4和图5所示,图中给出了光学相机随摆镜驱动后,等效光路图,用于后续分析的基础。带摆镜的TDICCD相机等同于对地相机,此相机可绕本体X轴以二倍摆镜摆角偏置,以扩展对地视场。
优选地,所述步骤3包括:获取卫星星下点太阳高度角信息,并在WGS84坐标系下获取卫星和地面目标的各项信息,包括:卫星位置、卫星速度、地面目标位置以及系统时间;此时瞬时惯性坐标系与对应的WGS84坐标系重合,每一个GPS测量节拍,都以当时WGS84坐标系的指向确定重合的瞬时惯性坐标系。
优选地,所述步骤4包括:
步骤4.1:判断星下点太阳高度角是否满足成像的光照条件,若满足则进入步骤4.2,若不满足则退出程序;
步骤4.2:利用卫星上的GPS测得的当前时刻的卫星位置和速度,结合地球自转角速度,将短时间的卫星运动简化为等速圆周运动,则推导出瞬时惯性坐标系、卫星轨道坐标系、卫星本体坐标系的转换关系;
步骤4.2.1:GPS测得的当前时刻在WGS84坐标系下的卫星位置和速度,计算公式如下:
式中:RWGS84表示GPS测得的当前时刻在WGS84坐标系下的卫星位置矢量,XRGPS表示RWGS84的X轴分量,YRGPS表示RWGS84的Y轴分量,ZRGPS表示RWGS84的Z轴分量;VWGS84表示GPS测得的当前时刻在WGS84坐标系下的卫星速度矢量,XVGPS表示VWGS84X轴分量,YVGPS表示VWGS84Y轴分量,ZVGPS表示VWGS84Z轴分量;
步骤4.2.2:根据WGS84坐标系下的卫星位置和速度求解在I系下,同时考虑地球自转引起的牵连速度时的卫星位置和速度,计算公式如下:
式中:TGPS表示卫星位置速度对应的时刻,分别表示I系下卫星位置和速度矢量,ωe表示地球自转角速率;
步骤4.2.3:根据I系下的卫星位置和速度求解卫星轨道角动量,计算公式如下:
式中:H表示卫星轨道角动量,HX表示卫星轨道角动量X轴分量,HY表示卫星轨道角动量Y轴分量,HZ表示卫星轨道角动量Z轴分量;
步骤4.2.4:根据卫星轨道角动量求解I系下卫星的倾角、升交点赤经、轨道角速度,计算公式如下:
xi=[1 0 0]T,zi=xi×N
zi(3)≥0时,Ω=arccos(N·xi),zi(3)<0时,Ω=-arccos(N·xi)
其中:
式中:i表示I系下卫星的轨道倾角、Ω表示升交点赤经、ω表示轨道角速度:xi表示惯性系X轴矢量,zi表示惯性系Z轴矢量,zi(3)表示zi在Z轴分量,N表示轨道节线单位矢量;
步骤4.2.5:根据I系下卫星的轨道角动量、倾角、轨道角速度求解TGPS时刻的轨道幅角u0,轨道幅角u,计算公式如下:
ZRGPS≥0时,ZRGPS<0时,
式中,t表示从起始时刻开始的时间差;
步骤4.2.6:根据I系下TGPS时刻的轨道幅角u0,轨道幅角u求解I系到卫星轨道坐标系的转换矩阵及变化率,计算公式如下:
式中:Aoi表示惯性系到轨道系转移矩阵,表示惯性系到轨道系转移矩阵变化率,ωb表示轨道角速率在本体系分量,ω表示轨道角速度;
步骤4.2.7:求解卫星轨道坐标系到卫星本体坐标系的转换矩阵,计算公式如下:
式中:Abo表示轨道系到问题系转移矩阵,θ表示俯仰姿态角,表示滚动自他角,ψ表示偏航姿态角。
具体地,步骤4中将卫星未来一段时间内的运行,简化为近圆轨道,认为短时间内轨道角速度不变,且高度不变,因此可以在满足计算精度要求的基础上,无需进行复杂的轨道六根数递推,把卫星运动简化为标准等角速率圆周运动,简化星上算法计算量。
优选地,所述步骤5包括:根据卫星和地面目标在在WGS84坐标系下的运动状态表征参数计算出地面目标在瞬时惯性坐标系下的坐标和变化率随时间变化的表达式,计算公式如下:
式中:表示目标初始位置,表示目标位置矢量惯性系分量,表示目标位置矢量变化率,ωe表示地球自转角速率。
优选地,所述步骤6包括:根据卫星和地面目标在瞬时惯性坐标系下的运动状态表征参数、地面目标坐标随时间变化参数以及地面目标变化率随时间变化参数计算出地面目标在卫星本体坐标系下的坐标和变化率随时间变化的表达式;计算公式如下:
式中:表示卫星到目标矢量在卫星本体系分量,X(t)表示X轴分量,Y(t)表示Y轴分量,Z(t)表示Z轴分量;表示卫星到目标矢量变化率在卫星本体系分量,X’(t)表示X轴分量,Y’(t)表示Y轴分量,Z’(t)表示Z轴分量。
优选地,所述步骤7包括:根据卫星本体坐标系下的坐标和变化率随时间变化参数建立约束方程,并求解出成像任务参数;
具体地,相机的后视角,代表成像的CCD阵列有效区域的中心,与相机Z轴夹角为θ1,假定未来某时刻,随着卫星飞行,卫星CCD阵列用于成像区域的中心在地面投影与目标点重合,是用于求解成像任务参数的主要约束条件。
步骤7.1:假设相机后视角为θ1,建立约束方程如下:
步骤7.2:将约束方程转换为F(t)=0的函数,函数表达式如下:
步骤7.3:求解函数F(t)的零点;
步骤7.3.1:对函数F(t)进行求导运算,得到如下的方程:
步骤7.3.2:利用牛顿迭代求解法计算时间变量t,计算公式如下:
步骤7.3.3:令t0的初始值为0,则t秒之后能够得到c系下地面目标相对卫星位置矢量Rst c;
步骤7.4:求解卫星成像任务参数,包括:摆镜摆角和相机成像时间TShoot;
若需要提前n秒成像,则相机成像时间为:
TShoot=TGPS+t-n,
相机光轴需绕滚动轴调整量为:
若反映到摆镜摆角上,则摆镜摆角为:
步骤7.5:判断是否满足成像条件,若满足则进入步骤7.6,若不满足则结束所述对地面点目标推扫成像任务参数的星载求解方法流程,其中所述成像条件为:在摆角最大摆动范围内,且t大于星上成像前的最小准备时间;
步骤7.6:星上准备好后,提前驱动相机的摆镜,并在TShoot时刻开始成像。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。
Claims (8)
1.一种对地面点目标推扫成像任务参数的星载求解方法,其特征在于,包括如下步骤:
步骤1:建立坐标系,所述坐标系包括瞬时惯性坐标系、卫星轨道坐标系、卫星本体坐标系、虚拟相机坐标系、虚像平面坐标系,并定义卫星、地面目标、地心三者间的矢量;
步骤2:利用建立的坐标系辅助建立包含摆镜的光学系统等效光路;
步骤3:获取卫星星下点太阳高度角,并判断满足成像的光照条件后,在包括卫星位置、卫星速度、地面目标位置以及系统时间在内的WGS84坐标系下获取卫星和地面目标的各项信息;
步骤4:利用卫星星下点太阳高度角信息,WGS84坐标系下卫星和地面目标的各项信息求解出卫星和地面目标在瞬时惯性坐标系下的运动状态表征参数,所述运动状态表征参数包括:卫星位置、卫星速度、卫星轨道角动量、卫星的倾角、升交点赤经、轨道角速度、实时轨道幅角、卫星轨道坐标系的转换矩阵及变化率以及卫星轨道坐标系到卫星本体坐标系的转换矩阵;
步骤5:根据卫星和地面目标在WGS84坐标系下的运动状态表征参数计算出地面目标在瞬时惯性坐标系下的坐标和变化率随时间变化的表达式;
步骤6:根据卫星和地面目标在瞬时惯性坐标系下的运动状态表征参数、地面目标坐标随时间变化参数以及地面目标变化率随时间变化参数计算出地面目标在卫星本体坐标系下的坐标和变化率随时间变化的表达式;
步骤7:根据卫星本体坐标系下的坐标和变化率随时间变化参数建立约束方程,并求解出成像任务参数,所述成像任务参数包括:相机指向调整参数和成像开始时刻。
2.根据权利要求1所述的对地面点目标推扫成像任务参数的星载求解方法,其特征在于,所述步骤1包括:
步骤1.1:建立坐标系,所述坐标系包括瞬时惯性坐标系、卫星轨道坐标系、卫星本体坐标系、虚拟相机坐标系、虚像平面坐标系;
-瞬时惯性坐标系Oexiyizi,Oexi轴在地球赤道平面内,指向GPS时刻的WGS84坐标系X轴方向;Oezi轴垂直于地球赤道平面,与地球自转角速度矢量方向一致;Oeyi与Oexi、Oezi轴垂直,且构成右手坐标系;所述坐标系简称I系;
-卫星轨道坐标系Osxoyozo,Os表示原点,即卫星质心;Osxo轴在卫星轨道平面内,指向卫星的运动方向;Oszo轴指向地心;Osyo与Osxo、Oszo轴垂直,且构成右手坐标系;所述坐标系简称O系;
-卫星本体坐标系Osxbybzb,Os表示原点,即卫星质心;Osxb轴沿卫星飞行方向,Oszb轴与Osyb轴、Osxb轴构成右手坐标系;当卫星处于标称零姿态运动时,Osxb、Osyb、Oszb轴分别与Osxo、Osyo、Oszo轴平行,且方向相同;所述坐标系简称b系;
-虚拟相机坐标系Ocxcyczc,Oc表示原点,即摆镜虚相机的摄像中心;Ocxcyc平面为物镜平面;Oczc轴沿光轴方向,摆镜无摆角时Ocxc、Ocyc轴与轨道系重合,构成右手坐标系;所述坐标系简称c系;
-虚像平面坐标系Opxpyp,Op表示原点,即相机像平面光轴中心处;相机像平面法线与Oczc轴平行,且Opxp、Opyp轴分别与Ocxc、Ocyc平行,且方向相同;所述坐标系简称p系;
步骤1.2:定义卫星、地面目标、地心三者间的矢量;
-Ret用于描述地面目标,表示从I系原点指向地面目标的矢量;
-Res用于描述卫星位置,表示从I系原点指向卫星质心的矢量;
-Rst用于描述地面目标相对卫星位置,表示从卫星质心指向地面目标的矢量;
-Rct用于描述地面目标相对相机位置,表示从c系原点指向地面目标的矢量;
表示在I系中的位置矢量Ret;
表示在轨道系O系中的位置矢量Res。
3.根据权利要求1所述的对地面点目标推扫成像任务参数的星载求解方法,其特征在于,所述步骤2包括:安装在卫星对地面上的相机光轴经过摆镜的反射后照向地面目标位置;当摆镜偏转角时,相当于相机偏转了角,则通过摆镜就能够实现对地面目标的推扫。
4.根据权利要求1所述的对地面点目标推扫成像任务参数的星载求解方法,其特征在于,所述步骤3包括:获取卫星星下点太阳高度角信息,并在WGS84坐标系下获取卫星和地面目标的各项信息,包括:卫星位置、卫星速度、地面目标位置以及系统时间;此时瞬时惯性坐标系与对应的WGS84坐标系重合,每一个GPS测量节拍,都以当时WGS84坐标系的指向确定重合的瞬时惯性坐标系。
5.根据权利要求1所述的对地面点目标推扫成像任务参数的星载求解方法,其特征在于,所述步骤4包括:
步骤4.1:判断星下点太阳高度角是否满足成像的光照条件,若满足则进入步骤4.2,若不满足则结束所述对地面点目标推扫成像任务参数的星载求解方法流程;
步骤4.2:利用卫星上的GPS测得的当前时刻的卫星位置和速度,结合地球自转角速度,将短时间的卫星运动简化为等速圆周运动,则推导出瞬时惯性坐标系、卫星轨道坐标系、卫星本体坐标系的转换关系;
步骤4.2.1:GPS测得的当前时刻在WGS84坐标系下的卫星位置和速度,计算公式如下:
式中:RWGS84表示GPS测得的当前时刻在WGS84坐标系下的卫星位置矢量,XRGPS表示RWGS84的X轴分量,YRGPS表示RWGS84的Y轴分量,ZRGPS表示RWGS84的Z轴分量;VWGS84表示GPS测得的当前时刻在WGS84坐标系下的卫星速度矢量,XVGPS表示VWGS84X轴分量,YVGPS表示VWGS84Y轴分量,ZVGPS表示VWGS84Z轴分量;
步骤4.2.2:根据WGS84坐标系下的卫星位置和速度求解在I系下,同时考虑地球自转引起的牵连速度时的卫星位置和速度,计算公式如下:
式中:TGPS表示卫星位置速度对应的时刻,分别表示I系下卫星位置和速度矢量,ωe表示地球自转角速率;
步骤4.2.3:根据I系下的卫星位置和速度求解卫星轨道角动量,计算公式如下:
式中:H表示卫星轨道角动量,HX表示卫星轨道角动量X轴分量,HY表示卫星轨道角动量Y轴分量,HZ表示卫星轨道角动量Z轴分量;
步骤4.2.4:根据卫星轨道角动量求解I系下卫星的倾角、升交点赤经、轨道角速度,计算公式如下:
xi=[100]T,zi=xi×N
zi(3)≥0时,Ω=arccos(N·xi),zi(3)<0时,Ω=-arccos(N·xi)
其中:
式中:i表示I系下卫星的轨道倾角、Ω表示升交点赤经、ω表示轨道角速度:xi表示惯性系X轴矢量,zi表示惯性系Z轴矢量,zi(3)表示zi在Z轴分量,N表示轨道节线单位矢量;
步骤4.2.5:根据I系下卫星的轨道角动量、倾角、轨道角速度求解TGPS时刻的轨道幅角u0,轨道幅角u,计算公式如下:
ZRGPS≥0时,ZRGPS<0时,
u=u0+ωt,
式中,t表示从起始时刻开始的时间差;
步骤4.2.6:根据I系下TGPS时刻的轨道幅角u0,轨道幅角u求解I系到卫星轨道坐标系的转换矩阵及变化率,计算公式如下:
式中:Aoi表示惯性系到轨道系转移矩阵,表示惯性系到轨道系转移矩阵变化率,ωb表示轨道角速率在本体系分量,ω表示轨道角速度;
步骤4.2.7:求解卫星轨道坐标系到卫星本体坐标系的转换矩阵,计算公式如下:
式中:Abo表示轨道系到问题系转移矩阵,θ表示俯仰姿态角,表示滚动自他角,ψ表示偏航姿态角。
6.根据权利要求1所述的对地面点目标推扫成像任务参数的星载求解方法,其特征在于,所述步骤5包括:根据卫星和地面目标在WGS84坐标系下的运动状态表征参数计算出地面目标在瞬时惯性坐标系下的坐标和变化率随时间变化的表达式,计算公式如下:
式中:表示目标初始位置,表示目标位置矢量惯性系分量,表示目标位置矢量变化率,ωe表示地球自转角速率。
7.根据权利要求1所述的对地面点目标推扫成像任务参数的星载求解方法,其特征在于,所述步骤6包括:根据卫星和地面目标在瞬时惯性坐标系下的运动状态表征参数、地面目标坐标随时间变化参数以及地面目标变化率随时间变化参数计算出地面目标在卫星本体坐标系下的坐标和变化率随时间变化的表达式;计算公式如下:
式中:表示卫星到目标矢量在卫星本体系分量,X(t)表示X轴分量,Y(t)表示Y轴分量,Z(t)表示Z轴分量;表示卫星到目标矢量变化率在卫星本体系分量,X’(t)表示X轴分量关于时间t的导数,Y’(t)表示Y轴分量关于时间t的导数,Z’(t)表示Z轴分量关于时间t的导数,Abo表示轨道系到问题系转移矩阵,Aoi表示惯性系到轨道系转移矩阵,表示惯性系到轨道系转移矩阵变化率,表示在轨道系中的位置矢量,表示目标位置矢量惯性系分量,表示目标位置矢量变化率。
8.根据权利要求2所述的对地面点目标推扫成像任务参数的星载求解方法,其特征在于,所述步骤7包括:根据卫星本体坐标系下的坐标和变化率随时间变化参数建立约束方程,并求解出成像任务参数;
步骤7.1:假设相机后视角为θ1,建立约束方程如下:
式中:X(t)表示X轴分量,Y(t)表示Y轴分量,Z(t)表示Z轴分量;
步骤7.2:将约束方程转换为F(t)=0的函数,函数表达式如下:
步骤7.3:求解函数F(t)的零点;
步骤7.3.1:对函数F(t)进行求导运算,得到如下的方程:
式中:X’(t)表示X轴分量关于时间t的导数,Y’(t)表示Y轴分量关于时间t的导数,Z’(t)表示Z轴分量关于时间t的导数;
步骤7.3.2:利用牛顿迭代求解法计算时间变量t,计算公式如下:
式中:tk表示第k时刻,tk+1表示第k+1时刻,F(tk)表示第k时刻的函数,F′(tk)表示F(tk)关于时间的导数;
步骤7.3.3:令t0的初始值为0,则t秒之后能够得到c系下地面目标相对卫星位置矢量Rst c;
步骤7.4:求解卫星成像任务参数,包括:摆镜摆角和相机成像时间TShoot;
若需要提前n秒成像,则相机成像时间为:
TShoot=TGPS+t-n,
式中:TGPS表示卫星位置速度对应的时刻;
相机光轴需绕滚动轴调整量为:
若反映到摆镜摆角上,则摆镜摆角为:
步骤7.5:判断是否满足成像条件,若满足则进入步骤7.6,若不满足则结束所述对地面点目标推扫成像任务参数的星载求解方法流程,其中所述成像条件为:在摆角最大摆动范围内,且t大于星上成像前的最小准备时间;
步骤7.6:星上准备好后,提前驱动相机的摆镜,并在TShoot时刻开始成像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510447888.8A CN105115477B (zh) | 2015-07-27 | 2015-07-27 | 对地面点目标推扫成像任务参数的星载求解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510447888.8A CN105115477B (zh) | 2015-07-27 | 2015-07-27 | 对地面点目标推扫成像任务参数的星载求解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105115477A CN105115477A (zh) | 2015-12-02 |
CN105115477B true CN105115477B (zh) | 2017-12-15 |
Family
ID=54663522
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510447888.8A Active CN105115477B (zh) | 2015-07-27 | 2015-07-27 | 对地面点目标推扫成像任务参数的星载求解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105115477B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110220492A (zh) * | 2019-06-06 | 2019-09-10 | 北京千乘探索科技有限公司 | 一种卫星星下点太阳高度角仿真分析方法及系统 |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106197434A (zh) * | 2016-06-30 | 2016-12-07 | 中国电子科技集团公司第五十四研究所 | 基于地面目标点位置的卫星姿态角的计算方法 |
CN107704424B (zh) * | 2017-09-22 | 2020-09-18 | 北京空间飞行器总体设计部 | 基于被动推扫成像模式下距离与所需时间的精确计算方法 |
CN109783514B (zh) * | 2018-12-26 | 2021-03-19 | 航天恒星科技有限公司 | 光学遥感卫星对地面目标的观测时间窗口快速计算方法 |
CN109974692B (zh) * | 2019-03-21 | 2021-08-10 | 北京控制工程研究所 | 一种基于中微子信号的隐蔽环境天文定位系统及方法 |
CN110162069B (zh) * | 2019-05-10 | 2020-09-08 | 北京航空航天大学 | 一种近地轨道航天器阳光反射凝视期望姿态解析求解方法 |
CN111090941A (zh) * | 2019-12-17 | 2020-05-01 | 哈尔滨工业大学 | 一种基于多目标优化算法的航天器最优Lambert轨道交会方法 |
CN111158020B (zh) * | 2020-01-06 | 2022-03-15 | 中国科学院微小卫星创新研究院 | 一种用于卫星的星载实时云判系统及方法 |
CN111949922B (zh) * | 2020-07-31 | 2023-11-28 | 上海卫星工程研究所 | 适用于对地探测任务多时间窗的星上快速计算的方法及系统 |
CN115014294A (zh) * | 2022-05-18 | 2022-09-06 | 河海大学 | 一种东西向卫星的河道敏捷成像方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103679711A (zh) * | 2013-11-29 | 2014-03-26 | 航天恒星科技有限公司 | 一种遥感卫星线阵推扫光学相机在轨外方位参数标定方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6271877B1 (en) * | 1999-06-25 | 2001-08-07 | Astrovision, Inc. | Direct broadcast imaging satellite system apparatus and method for providing real-time, continuous monitoring of earth from geostationary earth orbit |
-
2015
- 2015-07-27 CN CN201510447888.8A patent/CN105115477B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103679711A (zh) * | 2013-11-29 | 2014-03-26 | 航天恒星科技有限公司 | 一种遥感卫星线阵推扫光学相机在轨外方位参数标定方法 |
Non-Patent Citations (1)
Title |
---|
星载TDI-CCD推扫相机的偏流角计算与补偿;袁孝康;《上海航天》;20061231(第6期);第10-13页 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110220492A (zh) * | 2019-06-06 | 2019-09-10 | 北京千乘探索科技有限公司 | 一种卫星星下点太阳高度角仿真分析方法及系统 |
CN110220492B (zh) * | 2019-06-06 | 2020-09-29 | 北京千乘探索科技有限公司 | 一种卫星星下点太阳高度角仿真分析方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105115477A (zh) | 2015-12-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105115477B (zh) | 对地面点目标推扫成像任务参数的星载求解方法 | |
CN103674063B (zh) | 一种光学遥感相机在轨几何定标方法 | |
CN104567819B (zh) | 一种星载相机全视场偏流角确定与补偿方法 | |
JP6698670B2 (ja) | コリドー捕捉 | |
US7400950B2 (en) | Optical sensing system and system for stabilizing machine-controllable vehicles | |
CN105203101B (zh) | 一种基于目标天体星历修正的深空探测器捕获段天文导航方法 | |
CN104462776B (zh) | 一种低轨道地球观测卫星对月球绝对辐射定标方法 | |
CN105698764B (zh) | 一种光学遥感卫星影像时变系统误差建模补偿方法及系统 | |
CN107741229A (zh) | 一种光电/雷达/惯性组合的舰载机着舰导引方法 | |
CN105444781B (zh) | 星载自主引导成像地面验证方法 | |
CN106871932A (zh) | 基于金字塔搜索地形匹配的星载激光在轨指向检校方法 | |
CN105444780B (zh) | 星载摆扫式光学相机图像定位验证系统及处理方法 | |
CN105160125B (zh) | 一种星敏感器四元数的仿真分析方法 | |
CN103542981B (zh) | 一种双目视觉测量转动惯量的方法 | |
CN104764449B (zh) | 一种基于星历修正的捕获段深空探测器自主天文导航方法 | |
CN105651261B (zh) | 基于小卫星平台的双面阵立体测绘系统 | |
CN108344396A (zh) | 一种敏捷卫星斜条带成像模式姿态计算方法 | |
CN105606128A (zh) | 一种星载激光高度计外场检校方法 | |
CN108663043A (zh) | 基于单个相机辅助的分布式pos主子节点相对位姿测量方法 | |
CN110435930A (zh) | 低轨光学卫星匀降速推扫姿态规划方法 | |
CN102279001B (zh) | 一种星载相机的相移补偿方法 | |
CN105486315B (zh) | 遥感卫星对月绝对定标姿态调整方法 | |
CN112461204B (zh) | 卫星对动态飞行目标多视角成像联合计算航行高度的方法 | |
CN106289156B (zh) | 一种卫星以任意姿态成像时获取摄影点太阳高度角的方法 | |
CN104144304A (zh) | 一种高分辨率相机不同视场积分时间确定方法 |
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 |