CN106289156B - 一种卫星以任意姿态成像时获取摄影点太阳高度角的方法 - Google Patents
一种卫星以任意姿态成像时获取摄影点太阳高度角的方法 Download PDFInfo
- Publication number
- CN106289156B CN106289156B CN201610862773.XA CN201610862773A CN106289156B CN 106289156 B CN106289156 B CN 106289156B CN 201610862773 A CN201610862773 A CN 201610862773A CN 106289156 B CN106289156 B CN 106289156B
- Authority
- CN
- China
- Prior art keywords
- satellite
- coordinate system
- angle
- point
- attitude
- 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
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000003384 imaging method Methods 0.000 claims abstract description 32
- 239000011159 matrix material Substances 0.000 claims abstract description 19
- 238000009434 installation Methods 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims description 15
- 230000003287 optical effect Effects 0.000 description 12
- 230000001174 ascending effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000002310 reflectometry Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 210000001747 pupil Anatomy 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C1/00—Measuring angles
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Photovoltaic Devices (AREA)
- Stereoscopic And Panoramic Photography (AREA)
Abstract
本发明涉及一种卫星以任意姿态成像时获取摄影点太阳高度角的方法,步骤为:(1)根据卫星机动成像时的姿态矩阵和相机在卫星本体系的安装矩阵,计算相机本体坐标系相对于轨道坐标系的等效滚动角和俯仰角;(2)基于步骤(1)的结果和卫星轨道高度计算地面摄影点至星下点大弧段对应的地心角;(3)基于步骤(1)的结果计算摄影点至星下点大弧段与星下线前进方向的球面夹角;(4)根据前两步的结果,结合实时的轨道根数计算地面摄影点的赤经和赤纬;(5)根据步骤(4)的结果,结合姿轨控分系统自身可用的太阳方向矢量计算摄影点的太阳高度角。本发明方法可用于任意卫星姿态和任意相机安装方位条件下进行遥感成像时准确计算地面摄影点的太阳高度角。
Description
技术领域
本发明属于光学遥感卫星总体设计领域,涉及一种卫星以任意姿态成像时通过姿轨控分系统实时获取地面摄影点太阳高度角的方法。
背景技术
光学遥感卫星已经广泛应用于国土资源管理、环境监测、土地利用、城市规划和应急救灾等领域,它一般通过搭载一台或多台高分辨率的对地观测空间相机来获取遥感图像。空间相机成像时的入瞳辐亮度主要受光照条件、地物目标的反射率等因素影响。由于空间相机多采用推扫摄影方式,通常单次十几分钟左右的摄影过程覆盖的地物目标相距上千公里,横跨较大的纬度范围,这就会导致每个时刻的光照条件差别很大。为了获取层次分明、亮度统一的高质量遥感图像,可通过在成像过程中根据地面摄影点的太阳高度角和反射率对相机参数进行实时在轨调整来实现。对于线阵CCD相机来说,可调整的参数为增益;对于TDICCD相机来说,可调整的参数为增益和积分级数。
地面摄影点太阳高度角是空间相机参数在轨调整的重要依据,其准确性直接影响到相机参数调整的正确与否。太阳高度角的大小主要由地面摄影点与太阳的相对位置关系决定。现有的太阳高度角计算方法一般是基于地面摄影点的纬度、成像时刻的太阳赤纬角和太阳时角,依据三角函数关系公式计算获取(张云松,冯钟葵,石丹.卫星观测方位对遥感成像的影响.遥感学报,2007,11(4):433-438;武星星,刘金国,周怀得等.基于光照条件的空间相机增益在轨自动调整.光学学报,2014,34(03):032800-1-7)。这种方法的缺点是太阳赤纬角和太阳时角都是依据经验公式计算,未考虑太阳位置的实时变化,计算结果与实际值会存在偏差,导致太阳高度角计算结果的准确度不高。
此外,上述太阳高度角计算方法一般只考虑了星下点成像的情况,星下点太阳高度角的计算相对容易;而对于卫星在完成侧摆姿态机动后成像、俯仰姿态机动后成像或者侧摆-俯仰双轴姿态机动后成像的情况,以及空间相机的安装不满足相机光轴与卫星偏航轴平行的情况,尚缺少有效的获取地面摄影点太阳高度角的方法。
发明内容
本发明解决的技术问题是:克服现有技术的不足,针对光学遥感卫星在任意飞行姿态和任意相机安装方位条件下对地成像的复杂情况,提供一种通过姿轨控分系统实时获取地面摄影点太阳高度角的方法,可用于给空间相机在轨参数调整提供准确的依据。
本发明的技术解决方案是:一种卫星以任意姿态成像时获取摄影点太阳高度角的方法,包括如下步骤:
(1)根据卫星成像时的姿态矩阵A和相机在卫星本体坐标系的安装矩阵M,选用1-2-3欧拉角转序计算相机本体坐标系相对于轨道坐标系的等效滚动角和俯仰角θ,
θ=arcsin(A'31),
其中,A'是相机本体坐标系相对于轨道坐标系的姿态矩阵,A'=MA;(2)根据步骤(1)的计算结果,计算摄影点至星下点大弧段对应的地心角α,
其中,H为轨道高度,R为地球半径;
(3)根据步骤(1)的计算结果,计算摄影点至星下点大弧段与星下线前进方向的球面夹角β,
(4)根据步骤(2)和步骤(3)的计算结果和当前成像时刻的轨道根数计算地面摄影点的赤经RA和赤纬Dec,Ω为升交点赤经,i是卫星的轨道倾角,u=ω+f是当前时刻卫星离升交点的角矩,其中ω和f分别是当前时刻的近地点幅角和真近点角;当侧摆角时:
Dec=arcsin(sin∠B′OCsini),
其中,∠B′OC=u+sign(θ)*α,sign(·)为符号函数;
当侧摆角时:
其中,
cos∠B′OC=cosαcosu-sinαcosβsinu;
(5)根据步骤(4)的计算结果,结合姿轨控分系统各时刻可用的地心赤道惯性坐标系下的单位太阳方向矢量SI计算当前成像时刻摄影点的太阳高度角SunElevAng,
SunElevAng=90°-arccos(SI·NI),
NI=[cos(Dec)cos(RA)cos(Dec)sin(RA)sin(Dec)]T,
其中,NI为地心赤道惯性坐标系下当前时刻摄影点水平面的单位法向量(方向为指向球面外侧)。
本发明与现有技术相比的优点在于:现有的光学遥感卫星太阳高度角的获取方法中,往往只考虑了星下点成像且空间相机的光轴与卫星偏航轴平行的情况,并且在计算公式中要使用依据经验公式获取的太阳赤纬角和太阳时角,致使计算出来的太阳高度角存在偏差,从而导致不能给空间相机参数调整提供准确的输入。本发明方法提高了太阳高度角的计算精度,可满足光学遥感卫星多种场合的使用需求。首先在求解空间相机的姿态矩阵时,直接将卫星的姿态矩阵和相机的安装矩阵作为输入,使得本发明方法不仅适用于侧摆机动、俯仰机动以及侧摆和俯仰双轴机动成像的情况,也适用于相机光轴与卫星偏航轴不平行的情况。在太阳高度角的计算过程中,考虑了摄影点位置和太阳位置两者的实时变化,基于轨道根数并借助空间几何关系实时计算出地面摄影点的赤经和赤纬,结合姿轨控分系统现有的太阳方向矢量,仅在赤道惯性坐标系下就能准确的计算得到太阳高度角,无需经过任何坐标系之间的转换。此外,由于卫星单次成像任务中姿态角和轨道高度基本保持不变,在成像过程中本发明方法的前面三个步骤仅需在成像起始时刻执行一次,后续各时刻只执行步骤(4)和步骤(5),根据实时轨道根数更新摄影点的赤经、赤纬并推算出太阳高度角即可,其计算复杂度较小,适合星上在轨应用。本发明方法适用于任意卫星姿态和任意相机安装方位的情况,覆盖了目前光学遥感卫星的绝大多数情况,可为空间相机在轨参数调整提供可靠的输入。
附图说明
图1为本发明方法的流程框图;
图2为光学遥感卫星成像几何模型示意图;
图3为卫星姿态机动后成像星下点与地面摄影点位置关系图。
具体实施方式
如图1所示,为本发明方法的流程框图,(1)根据卫星机动成像时的姿态矩阵和相机在卫星本体系的安装矩阵,计算相机本体坐标系相对于轨道坐标系的等效滚动角和俯仰角;(2)基于步骤(1)的结果和卫星轨道高度计算地面摄影点至星下点大弧段对应的地心角;(3)基于步骤(1)的结果计算摄影点至星下点大弧段与星下线前进方向的球面夹角;(4)根据前两步的结果,结合实时的轨道根数计算地面摄影点的赤经和赤纬;(5)根据步骤(4)的结果,结合姿轨控分系统自身可用的太阳方向矢量计算摄影点的太阳高度角。本发明方法可用于任意卫星姿态和任意相机安装方位条件下进行遥感成像时准确计算地面摄影点的太阳高度角。
具体实现步骤如下:
(1)相机本体坐标系相对于轨道坐标系的等效侧摆角与俯仰角计算
本发明根据光学遥感卫星成像原理,可以建立如图2所示的成像几何模型:卫星运行在轨道高度为H的轨道上,卫星的轨道坐标系为S-XoYoZo,S为卫星的质心,Zo轴指向地心,Xo轴指向飞行方向,Yo由右手定则确定,星下点为B,地面摄影点为B′。假设相机在卫星本体坐标系的安装矩阵为M,卫星本体坐标系初始时刻与轨道坐标系重合,当卫星成像时的姿态矩阵为A。则相机相对于轨道坐标系的姿态矩阵为:
A'=MA
如图2所示,选用1-2-3欧拉角转序,即将轨道坐标系先绕Xo轴旋转角(滚动角),再绕经过第一次旋转后的Yo轴旋转θ角(俯仰角),再绕经过两次旋转后的Zo轴旋转φ角(偏航角)。这里的滚动角、俯仰角和偏航角不是卫星本体姿态角,而是相机本体坐标系相对于轨道坐标系的等效姿态角。则有:
由于这两个角度可由下式计算:
θ=arcsin(A'31)
(2)摄影点至星下点大弧段的地心角计算
在图2的三角形SOB′中,
因此,摄影点至星下点大弧段对应的地心角
其中,H为轨道高度,R为地球半径。
(3)摄影点至星下点大弧段与星下线前进方向的球面夹角计算
根据相机等效侧摆角和俯仰角θ,可以计算得到摄影点至星下点大弧段与星下线前进方向的球面夹角β,β的取值范围是[0,π],
(4)当前成像时刻摄影点的赤经和赤纬计算
如图3所示,O-XIYIZI是地心赤道惯性坐标系。在当前成像时刻,B点为星下点,B′点为地面摄影点,C点为升交点。实时轨道根数为:升交点赤经Ω,卫星的轨道倾角i,卫星离升交点的角矩u=ω+f,其中ω和f分别是近地点幅角和真近点角。卫星在轨工作期间,姿轨控分系统会通过地面站定期注入精轨数据或者通过引入GPS轨道数据的方式来确保当前轨道数据的准确性。
a)如果侧摆角摄影点B′位于BC圆弧上
∠B′OC=u+sign(θ)*α
其中,sign(·)为符号函数。
在直角球面三角形CB′D′中,由球面三角形的正弦定理和余弦定理可得到当前时刻摄影点B′的赤经RA和赤纬Dec:
Dec=arcsin(sin∠B′OCsini)
b)如果侧摆角摄影点B′位于BC圆弧以外
此时,球面角∠B′BC=π-β。在球面三角形B′BC中,由球面三角形正弦定理和余弦定理可以推出:
cos∠B′OC=cosαcosu-sinαcosβsinu
在球面三角形B′CD′中,由球面三角形正弦定理和余弦定理可以推出:
将上面两式展开化简得到:
进一步可得到当前时刻摄影点B′的赤经RA和赤纬Dec:
在单次成像任务期间,等效滚动角轨道倾角i、地心角α和球面夹角β可以认为是常值,它们对应的三角函数值仅需在成像起始时刻计算一次,因此计算量比较小。
(5)当前成像时刻摄影点的太阳高度角计算
太阳高度角的定义为太阳方向矢量与地面摄影点水平面的夹角。
考虑到太阳敏感器是光学遥感卫星控制分系统标配的姿态敏感器之一,在使用太阳敏感器时,姿轨控分系统需要根据太阳星历实时计算地心赤道惯性坐标系下的单位太阳方向矢量SI。因此,每一时刻SI的大小都是可供直接使用的。
在地心赤道惯性坐标系下,当前时刻摄影点B′点水平面的单位法向量(方向为指向球面外侧)可以表示为:
NI=[cos(Dec)cos(RA) cos(Dec)sin(RA) sin(Dec)]T
当前时刻摄影点B′点的太阳高度角为:
SunElevAng=90°-arccos(SI·NI)。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (6)
1.一种卫星以任意姿态成像时获取摄影点太阳高度角的方法,其特征在于包括以下步骤:
(1)根据卫星成像时的姿态矩阵A和相机在卫星本体坐标系的安装矩阵M,计算相机本体坐标系相对于轨道坐标系的等效滚动角和俯仰角θ;
(2)根据步骤(1)的等效滚动角和俯仰角θ,结合卫星轨道高度计算摄影点至星下点大弧段对应的地心角α;
(3)根据步骤(1)的计算结果,计算摄影点至星下点大弧段与星下线前进方向的球面夹角β;
(4)根据步骤(2)的所述地心角α、步骤(3)所述球面夹角β和当前成像时刻的轨道根数,计算地面摄影点的赤经RA和赤纬Dec;
(5)根据步骤(4)的计算结果,结合姿轨控分系统各时刻可用的地心赤道惯性坐标系下的单位太阳方向矢量SI计算当前成像时刻摄影点的太阳高度角SunElevAng。
2.根据权利要求1所述的一种卫星以任意姿态成像时获取摄影点太阳高度角的方法,其特征在于:所述步骤(1)中,选用1-2-3欧拉角转序计算相机本体坐标系相对于轨道坐标系的等效滚动角和俯仰角θ,
θ=arcsin(A'31),
其中,A'是相机本体坐标系相对于轨道坐标系的姿态矩阵,A'=MA,A'32:3-2转序时相机本体坐标系相对于轨道坐标系的姿态矩阵,A'33:3-3转序时相机本体坐标系相对于轨道坐标系的姿态矩阵,A'31:3-1转序时相机本体坐标系相对于轨道坐标系的姿态矩阵。
3.根据权利要求1所述的一种卫星以任意姿态成像时获取摄影点太阳高度角的方法,其特征在于:所述步骤(2)中,所述地心角α计算如下:
其中,H为轨道高度,R为地球半径。
4.根据权利要求1所述的一种卫星以任意姿态成像时获取摄影点太阳高度角的方法,其特征在于:所述步骤(3)中,所述球面夹角β计算如下:
5.根据权利要求1所述的一种卫星以任意姿态成像时获取摄影点太阳高度角的方法,其特征在于:所述步骤(4)中,赤经RA和赤纬Dec计算如下:
Ω为升交点赤经,i是卫星的轨道倾角,u=ω+f是当前时刻卫星离升交点的角距,其中ω和f分别是当前时刻的近地点幅角和真近点角;
当时:
Dec=arcsin(sin∠B′OCsini),
其中,∠B′OC=u+sign(θ)*α,sign()为符号函数;
当时:
其中,
cos∠B′OC=cosαcosu-sinαcosβsinu。
6.根据权利要求1所述的一种卫星以任意姿态成像时获取摄影点太阳高度角的方法,其特征在于:所述步骤(5)中,太阳高度角SunElevAng计算如下:
SunElevAng=90°-arccos(SI·NI),
NI=[cos(Dec)cos(RA) cos(Dec)sin(RA) sin(Dec)]T,
其中,NI为地心赤道惯性坐标系下当前时刻摄影点水平面的单位法向量,即方向为指向球面外侧。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610862773.XA CN106289156B (zh) | 2016-09-28 | 2016-09-28 | 一种卫星以任意姿态成像时获取摄影点太阳高度角的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610862773.XA CN106289156B (zh) | 2016-09-28 | 2016-09-28 | 一种卫星以任意姿态成像时获取摄影点太阳高度角的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106289156A CN106289156A (zh) | 2017-01-04 |
CN106289156B true CN106289156B (zh) | 2019-11-29 |
Family
ID=57715761
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610862773.XA Active CN106289156B (zh) | 2016-09-28 | 2016-09-28 | 一种卫星以任意姿态成像时获取摄影点太阳高度角的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106289156B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108803642B (zh) * | 2018-06-14 | 2021-04-13 | 上海卫星工程研究所 | 光学成像卫星相机的太阳保护姿控关联设计方法 |
CN110220492B (zh) * | 2019-06-06 | 2020-09-29 | 北京千乘探索科技有限公司 | 一种卫星星下点太阳高度角仿真分析方法及系统 |
CN112857306B (zh) * | 2020-12-31 | 2022-12-13 | 航天东方红卫星有限公司 | 一种视频卫星任意视向点的连续太阳高度角确定方法 |
CN114187351B (zh) * | 2021-12-09 | 2022-07-26 | 北京劢亚科技有限公司 | 一种应用于卫星的图像采集方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10221110A (ja) * | 1997-02-10 | 1998-08-21 | Nec Eng Ltd | 太陽角度検出機能付き突き抜け走査型地球センサ |
EP1106505A2 (en) * | 1999-12-10 | 2001-06-13 | Nec Corporation | Attitude angle sensor correcting apparatus for an artificial satellite |
CN104848860A (zh) * | 2015-05-19 | 2015-08-19 | 北京空间飞行器总体设计部 | 一种敏捷卫星成像过程姿态机动规划方法 |
CN105953803A (zh) * | 2016-04-25 | 2016-09-21 | 上海航天控制技术研究所 | 数字式太阳敏感器测量坐标系与棱镜坐标系偏差测量方法 |
-
2016
- 2016-09-28 CN CN201610862773.XA patent/CN106289156B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10221110A (ja) * | 1997-02-10 | 1998-08-21 | Nec Eng Ltd | 太陽角度検出機能付き突き抜け走査型地球センサ |
EP1106505A2 (en) * | 1999-12-10 | 2001-06-13 | Nec Corporation | Attitude angle sensor correcting apparatus for an artificial satellite |
CN104848860A (zh) * | 2015-05-19 | 2015-08-19 | 北京空间飞行器总体设计部 | 一种敏捷卫星成像过程姿态机动规划方法 |
CN105953803A (zh) * | 2016-04-25 | 2016-09-21 | 上海航天控制技术研究所 | 数字式太阳敏感器测量坐标系与棱镜坐标系偏差测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106289156A (zh) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106124170A (zh) | 一种基于高精度姿态信息的相机光轴指向计算方法 | |
CN103983254B (zh) | 一种新型敏捷卫星机动中成像方法 | |
US5072396A (en) | Navigation systems | |
CN106289156B (zh) | 一种卫星以任意姿态成像时获取摄影点太阳高度角的方法 | |
US5809457A (en) | Inertial pointing and positioning system | |
CN107450582B (zh) | 一种基于星上实时规划的相控阵数传引导控制方法 | |
CN113720360B (zh) | 一种敏捷光学遥感卫星星/地相机夹角在轨标定方法 | |
CN108663052B (zh) | 一种星上自主空间非合作目标相对导航相机指向控制方法 | |
CN107600464A (zh) | 利用星敏感器信息的飞轮控制捕获太阳及对日定向方法 | |
CN108710379A (zh) | 静止卫星成像偏航导引角计算方法 | |
CN101825475A (zh) | 空间光学遥感传感器的像移补偿方法 | |
CN105160125B (zh) | 一种星敏感器四元数的仿真分析方法 | |
CN107525492B (zh) | 一种适用于敏捷对地观测卫星的偏流角仿真分析方法 | |
CN105115477A (zh) | 对地面点目标推扫成像任务参数的星载求解方法 | |
CN103389099A (zh) | 基于x射线脉冲星的航天器姿态和位置测量系统及方法 | |
CN114679541B (zh) | 一种星上运动目标跟踪方法 | |
CN109178345B (zh) | 一种用于空中动目标跟踪的云台指向及星体姿态协同控制方法 | |
CN111381256A (zh) | 主动遥感卫星天线相位中心偏移误差计算的方法和系统 | |
CN103941740B (zh) | 一种考虑地球椭率的多轴机动成像卫星偏航姿态控制方法 | |
CN103955138A (zh) | 一种基于增量式偏流角的动中成像卫星姿态控制方法 | |
CN109506662B (zh) | 一种小天体着陆初始对准方法、其相对导航基准确定方法及装置 | |
CN102279001A (zh) | 一种星载相机的相移补偿方法 | |
CN109708667A (zh) | 一种基于激光陀螺的双动态目标跟踪引导方法 | |
CN111121765B (zh) | 共用遥感卫星平台的对地灵巧观测与对天geo目标监视方法 | |
CN108489483B (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 |