CN102819019A - 一种卫星波束与地球交点坐标的确定方法 - Google Patents
一种卫星波束与地球交点坐标的确定方法 Download PDFInfo
- Publication number
- CN102819019A CN102819019A CN2012102523300A CN201210252330A CN102819019A CN 102819019 A CN102819019 A CN 102819019A CN 2012102523300 A CN2012102523300 A CN 2012102523300A CN 201210252330 A CN201210252330 A CN 201210252330A CN 102819019 A CN102819019 A CN 102819019A
- Authority
- CN
- China
- Prior art keywords
- satellite
- coordinate system
- cos
- sin
- axle
- 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
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
一种卫星波束与地球交点坐标的确定方法,步骤为:(1)求取卫星波束在卫星本体坐标系中的单位矢量rfs;(2)利用卫星的偏航角、俯仰角、横滚角计算得到卫星本体坐标系与卫星轨道坐标系的转换矩阵(3)利用卫星的经度λ和纬度l计算得到地理坐标系与地球坐标系的转换矩阵(4)利用卫星所在子午面和卫星瞬时轨道面的夹角β以及卫星运动方向与地理水平面的夹角α计算得到辅助坐标系与地理坐标系的转换矩阵(5)通过坐标系转换,将单位矢量rfs转换到地球坐标系中得到rfe;(6)将卫星波束入射的直线参数方程与地球椭球方程联立,计算得到两组坐标[xje,yje,zje]T,其中离卫星较近的一组坐标即为卫星波束与地球交点的坐标。
Description
技术领域
本发明属于空间微波遥感技术领域,涉及一种确定卫星波束与地球位置关系的方法。
背景技术
星载微波散射计是一种测量目标后向散射系数σ°的典型雷达系统,主要用于海面风场测量,它是目前能够同时测量海面风速和风向的重要遥感仪器之一。由σ°反演海面风速和风向时,通过确定卫星波束与地球交点的坐标可以得到σ°对应的测量点经纬度,带入反演模型后可以计算出海面的风速和风向。
目前,国外主要采用矢量相乘的方法建立卫星轨道坐标系和地球坐标系的转换矩阵,进而计算卫星波束与地球的交点坐标。但是相关文献中仅给出了矢量从卫星轨道坐标系至地球坐标系的转换结果,没有详细的转换过程和说明,并且该方法没有说明位置矢量差的计算过程,因此无法在实际计算中应用。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供了一种卫星波束与地球交点坐标的确定方法,在卫星数据反演时能提高反演数据的可靠度。
本发明的技术解决方案是:一种卫星波束与地球交点坐标的确定方法,步骤如下:
(1)求取卫星的一个波束在卫星本体坐标系中的单位矢量rfs,其中φs为卫星的视角,为卫星的方位角;所述的卫星本体坐标系用oxsyszs表示,原点o为卫星中心,xs轴指向卫星星体纵轴,zs轴指向卫星运动时面向地球的方向,ys轴按右手正交定则定义;
所述的卫星轨道坐标系用oxayaza表示,原点o为卫星中心,xa、za轴在卫星轨道平面内,xa轴指向卫星运动方向,za轴指向地心,ya轴按右手正交定则定义;
所述的地理坐标系用oxtytzt表示,原点o为卫星中心,xt、yt轴在地理水平面内,xt轴指向东,yt轴指向北,zt轴在地球质心与卫星的连线上,按右手正交定则定义;所述的地球坐标系用oexeyeze表示,原点oe为地球质心,ze轴指向地球北极方向,xe轴指向零度子午面和协议地球极赤道的交点,ye轴和ze、xe轴构成右手坐标系;
所述的辅助坐标系用oxbybzb表示,原点o为卫星中心,yb、zb轴在卫星轨道平面内,yb轴指向卫星运动方向,zb轴与地球质心至卫星连线的方向相同,xb轴按右手正交定则定义;
(5)通过坐标系转换,将单位矢量rfs转换到地球坐标系中,得到方向矢量rfe, 其中
(6)根据卫星波束起始点在地球坐标系中的坐标[xse,yse,zse]T和波束入射方向矢量rfe,得到卫星波束入射的直线参数方程,该直线参数与地球椭球方程联立,计算得到两组坐标[xje,yje,zje]T,
式中u为直线参数方程所表示直线的斜率,式中参数a为地球赤道半径、b为地球极半径,两组坐标中离卫星较近的一组坐标即为卫星波束与地球交点的坐标。
本发明与现有技术相比的优点在于:本发明方法以一个辅助坐标系为“桥梁”,通过坐标系转换的方法,将卫星相关坐标系(卫星本体坐标系、卫星轨道坐标系)与导航相关坐标系(地理坐标系、地球坐标系)关联起来,建立卫星本体坐标系与地球坐标系的转换矩阵,最终完成交点坐标的计算。该方法过程简便,易于工程实现。采用本发明方法可以方便的获得卫星波束与地球交点的坐标,从而在卫星数据反演时,提高反演数据的可靠度,使得海面风速和风向的确定更加准确可靠。
附图说明
图1为本发明方法的流程框图;
图2为本发明方法中辅助坐标系示意图;
图3为本发明方法中辅助坐标系与地理坐标系的关系示意图。
具体实施方式
如图1所示,本发明方法以一个辅助坐标系为“桥梁”,将卫星相关坐标系(卫星本体坐标系、卫星轨道坐标系)与导航相关坐标系(地理坐标系、地球坐标系)关联起来,建立卫星本体坐标系与地球坐标系的转换矩阵,最终完成交点坐标的计算。
下面先给出本发明步骤中涉及到的坐标系的定义及其之间的转换矩阵。
用oxsyszs表示,原点o为卫星中心,xs轴指向卫星星体纵轴,与卫星在轨运动方向相同,zs轴指向卫星运动时面向地球的方向,ys轴按右手正交定则定义。
用oxayaza表示,原点o为卫星中心,xa、za轴在卫星轨道平面内,xa轴指向卫星运动方向,za轴指向地心,ya轴按右手正交定则定义。当卫星的姿态角为0°时,卫星本体坐标系与卫星轨道坐标系重合。
用oxtytzt表示,原点o为卫星中心,xt、yt轴在地理水平面内,xt轴指向东,yt轴指向北,zt轴在地球质心与卫星的连线上,按右手正交定则定义。
用oexeyeze表示,原点oe为地球质心,ze轴指向地球北极方向,xe轴指向零度子午面和协议地球极赤道的交点,ye轴和ze、xe轴构成右手坐标系。
用oxbybzb表示,原点o为卫星中心,yb、zb轴在卫星轨道平面内,yb轴指向卫星运动方向,zb轴与地球质心至卫星连线的方向相同,xb轴按右手正交定则定义,如图2所示。该坐标系在卫星轨道面内,原型来自地理坐标系,当瞬时轨道倾角为90°,且卫星运动方向水平时,辅助坐标系与地理坐标系重合。采用辅助坐标系作为卫星轨道坐标系和地理坐标系的中间坐标系,使卫星轨道坐标系至地理坐标系的转换矩阵建立过程清晰,便于掌握。
由卫星的姿态信息(偏航角ψ、俯仰角θ、横滚角γ)可以得到卫星本体坐标系与卫星轨道坐标系的转换矩阵:
由卫星的经度λ和纬度l可以得到地理坐标系与地球坐标系的转换矩阵:
式(1)、式(2)的计算结果见文献:陈哲.捷联惯导系统原理.北京:宇航出版社。
卫星波束在卫星本体坐标系下的方向矢量rfs,经过姿态矩阵转换后得到卫星波束在卫星轨道坐标系下的方向矢量rfa。而方向矢量rfa经过卫星轨道坐标系与地球坐标系的转换矩阵(转换矩阵)转换后,可以得到卫星波束在地球坐标系下的方向矢量rfe。由方向矢量rfe和卫星位置可以得到波束直线方程,该波束直线方程与地球椭球方程联立求解,即可得到卫星波束与地球交点的坐标。
(1)根据上面坐标系的定义部分可知,卫星轨道坐标系和辅助坐标系的关系为:xb=ya,yb=xa,zb=-za,由此可以得到这两个坐标系之间的转换矩阵为:
(2)辅助坐标系与地理坐标系的关系示意图如图3所示。两个坐标系通过两个参数关联起来:卫星所在子午面和瞬时轨道面的夹角β,以及卫星运动方向与地理水平面的夹角α。
如图3所示,S为当前时刻GPS确定的卫星位置[xse,yse,zse]T,So为S之前的最近时间由GPS确定的卫星位置。由于地球在自转,在惯性空间里,地球坐标系以24小时的周期绕ze轴匀速转动,所以,So时刻至S时刻地球坐标系绕ze轴转动了ωet的角度。定义So时刻的地球坐标系为e′系,S时刻的地球坐标系为e系,e′系绕ze′轴转ωet后与e系重合。So时刻卫星位置表示为[xs0e′,ys0e′,zs0e′]T,则So时刻卫星位置在e系中表示为:
其中ωe为地球自转角速度,t为So与S点所对应的GPS时间差。则矢量SoS为:
该波束在地球坐标系中的方向矢量rfe为:
已知波束起始点(卫星位置)在地球坐标系中的坐标S=[xse,yse,zse]T和波束入射方向rfe,可以确定波束入射的直线参数方程为:
B(u)=S+rfe·u (13)
式中,u为参数。该直线与地球椭球面有两个交点,通过与椭球方程联立可以计算出两个坐标:
式中参数a为地球赤道半径、b地球极半径。两个坐标中,离卫星较近的即为卫星波束与地球的交点坐标[xje,yje,zje]T。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (1)
1.一种卫星波束与地球交点坐标的确定方法,其特征在于步骤如下:
(1)求取卫星的一个波束在卫星本体坐标系中的单位矢量rfs,其中φs为卫星的视角,为卫星的方位角;所述的卫星本体坐标系用oxsyszs表示,原点o为卫星中心,xs轴指向卫星星体纵轴,zs轴指向卫星运动时面向地球的方向,ys轴按右手正交定则定义;
所述的卫星轨道坐标系用oxayaza表示,原点o为卫星中心,xa、za轴在卫星轨道平面内,xa轴指向卫星运动方向,za轴指向地心,ya轴按右手正交定则定义;
所述的地理坐标系用oxtytzt表示,原点o为卫星中心,xt、yt轴在地理水平面内,xt轴指向东,yt轴指向北,zt轴在地球质心与卫星的连线上,按右手正交定则定义;所述的地球坐标系用oexeyeze表示,原点oe为地球质心,ze轴指向地球北极方向,xe轴指向零度子午面和协议地球极赤道的交点,ye轴和ze、xe轴构成右手坐标系;
所述的辅助坐标系用oxbybzb表示,原点o为卫星中心,yb、zb轴在卫星轨道平面内,yb轴指向卫星运动方向,zb轴与地球质心至卫星连线的方向相同,xb轴按右手正交定则定义;
(6)根据卫星波束起始点在地球坐标系中的坐标[xse,yse,zse]T和波束入射方向矢量rfe,得到卫星波束入射的直线参数方程,该直线参数与地球椭球方程联立,计算得到两组坐标[xje,yje,zje]T,
式中u为直线参数方程所表示直线的斜率,式中参数a为地球赤道半径、b为地球极半径,两组坐标中离卫星较近的一组坐标即为卫星波束与地球交点的坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012102523300A CN102819019B (zh) | 2012-07-20 | 2012-07-20 | 一种卫星波束与地球交点坐标的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012102523300A CN102819019B (zh) | 2012-07-20 | 2012-07-20 | 一种卫星波束与地球交点坐标的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102819019A true CN102819019A (zh) | 2012-12-12 |
CN102819019B CN102819019B (zh) | 2013-11-20 |
Family
ID=47303244
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012102523300A Active CN102819019B (zh) | 2012-07-20 | 2012-07-20 | 一种卫星波束与地球交点坐标的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102819019B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103413006A (zh) * | 2013-08-27 | 2013-11-27 | 中国科学院空间科学与应用研究中心 | 一种空间惯性定向姿态卫星数传天线波束的设计方法 |
CN107831502A (zh) * | 2017-10-17 | 2018-03-23 | 中国科学院电子学研究所 | 卫星指向夹角近似值的计算方法 |
CN110531328A (zh) * | 2019-06-25 | 2019-12-03 | 西安空间无线电技术研究所 | 一种星载环扫雷达波束中心指向误差估计方法 |
CN112269197A (zh) * | 2020-10-12 | 2021-01-26 | 合肥四相至新科技有限公司 | 坐标转换方法及其系统 |
CN113640799A (zh) * | 2021-08-11 | 2021-11-12 | 北京无线电测量研究所 | 一种雷达波束中心照射点确定方法、装置及存储介质 |
CN114200537A (zh) * | 2021-12-13 | 2022-03-18 | 北京环境特性研究所 | 一种卫星下视地球的探测边界确定方法 |
CN115276760A (zh) * | 2022-06-24 | 2022-11-01 | 中国联合网络通信集团有限公司 | 波束中心的位置确定方法、装置及计算机存储介质 |
CN115622594A (zh) * | 2022-10-10 | 2023-01-17 | 中国空间技术研究院 | 一种基于传感器视角的卫星通信点波束覆盖边缘计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6225942B1 (en) * | 1999-07-30 | 2001-05-01 | Litton Systems, Inc. | Registration method for multiple sensor radar |
CN101414003A (zh) * | 2008-11-28 | 2009-04-22 | 北京航空航天大学 | 一种基于星地坐标转换的星载sar图像地理编码方法 |
CN102565799A (zh) * | 2012-01-31 | 2012-07-11 | 北京航空航天大学 | 一种多平台多模式sar回波的统一仿真实现方法 |
CN102565797A (zh) * | 2011-12-21 | 2012-07-11 | 北京航空航天大学 | 一种针对聚束模式星载sar图像的几何校正方法 |
-
2012
- 2012-07-20 CN CN2012102523300A patent/CN102819019B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6225942B1 (en) * | 1999-07-30 | 2001-05-01 | Litton Systems, Inc. | Registration method for multiple sensor radar |
CN101414003A (zh) * | 2008-11-28 | 2009-04-22 | 北京航空航天大学 | 一种基于星地坐标转换的星载sar图像地理编码方法 |
CN102565797A (zh) * | 2011-12-21 | 2012-07-11 | 北京航空航天大学 | 一种针对聚束模式星载sar图像的几何校正方法 |
CN102565799A (zh) * | 2012-01-31 | 2012-07-11 | 北京航空航天大学 | 一种多平台多模式sar回波的统一仿真实现方法 |
Non-Patent Citations (2)
Title |
---|
李德志等: "卫星点波束覆盖区域算法研究", 《载人航天》 * |
陈哲: "《捷联惯导系统原理》", 31 December 1986 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103413006B (zh) * | 2013-08-27 | 2016-06-15 | 中国科学院空间科学与应用研究中心 | 一种空间惯性定向姿态卫星数传天线波束的设计方法 |
CN103413006A (zh) * | 2013-08-27 | 2013-11-27 | 中国科学院空间科学与应用研究中心 | 一种空间惯性定向姿态卫星数传天线波束的设计方法 |
CN107831502B (zh) * | 2017-10-17 | 2023-05-09 | 中国科学院电子学研究所 | 卫星指向夹角近似值的计算方法 |
CN107831502A (zh) * | 2017-10-17 | 2018-03-23 | 中国科学院电子学研究所 | 卫星指向夹角近似值的计算方法 |
CN110531328A (zh) * | 2019-06-25 | 2019-12-03 | 西安空间无线电技术研究所 | 一种星载环扫雷达波束中心指向误差估计方法 |
CN110531328B (zh) * | 2019-06-25 | 2021-09-07 | 西安空间无线电技术研究所 | 一种星载环扫雷达波束中心指向误差估计方法 |
CN112269197A (zh) * | 2020-10-12 | 2021-01-26 | 合肥四相至新科技有限公司 | 坐标转换方法及其系统 |
CN112269197B (zh) * | 2020-10-12 | 2024-04-19 | 合肥四相至新科技有限公司 | 坐标转换方法及其系统 |
CN113640799A (zh) * | 2021-08-11 | 2021-11-12 | 北京无线电测量研究所 | 一种雷达波束中心照射点确定方法、装置及存储介质 |
CN114200537B (zh) * | 2021-12-13 | 2023-12-22 | 北京环境特性研究所 | 一种卫星下视地球的探测边界确定方法 |
CN114200537A (zh) * | 2021-12-13 | 2022-03-18 | 北京环境特性研究所 | 一种卫星下视地球的探测边界确定方法 |
CN115276760A (zh) * | 2022-06-24 | 2022-11-01 | 中国联合网络通信集团有限公司 | 波束中心的位置确定方法、装置及计算机存储介质 |
CN115276760B (zh) * | 2022-06-24 | 2023-06-23 | 中国联合网络通信集团有限公司 | 波束中心的位置确定方法、装置及计算机存储介质 |
CN115622594A (zh) * | 2022-10-10 | 2023-01-17 | 中国空间技术研究院 | 一种基于传感器视角的卫星通信点波束覆盖边缘计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102819019B (zh) | 2013-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102819019B (zh) | 一种卫星波束与地球交点坐标的确定方法 | |
CN102565797B (zh) | 一种针对聚束模式星载sar图像的几何校正方法 | |
CN104848860B (zh) | 一种敏捷卫星成像过程姿态机动规划方法 | |
CN104880191B (zh) | 一种基于太阳矢量的偏振辅助导航方法 | |
CN101339244B (zh) | 一种机载sar图像自动目标定位方法 | |
CN103675773B (zh) | 一种定标器与卫星指向对准的确定方法 | |
CN106124170A (zh) | 一种基于高精度姿态信息的相机光轴指向计算方法 | |
CN103940446B (zh) | 一种基于横坐标系的极区航行船舶捷联惯性导航系统重调方法 | |
CN106197425A (zh) | 基于卫星姿态角的地面目标点位置的计算方法 | |
CN105184002A (zh) | 一种数传天线指向角度的仿真分析方法 | |
CN104880192A (zh) | 一种基于偏振罗盘的载体航向角计算方法 | |
CN103985952B (zh) | 船载a-e-c三轴卫星通信天线极化偏差角实时修正方法 | |
CN104215242A (zh) | 一种基于横向游移坐标系的极区惯性导航方法 | |
CN106292677A (zh) | 基于恒星时角的姿态控制方法和系统 | |
CN104729457A (zh) | 太阳相对近地轨道微小卫星位置的确定方法 | |
CN104749598A (zh) | 一种产生gnss掩星路径的方法 | |
CN103675760A (zh) | 一种星载地球同步轨道合成孔径雷达姿态导引方法 | |
CN103344958B (zh) | 基于星历数据的星载sar高阶多普勒参数估算方法 | |
CN103940429A (zh) | 一种惯性导航系统横坐标系下载体姿态的实时测量方法 | |
CN105021212B (zh) | 一种初始方位信息辅助下潜航器快速传递对准方法 | |
CN109470276A (zh) | 基于零速修正的里程计标定方法与装置 | |
CN103471614A (zh) | 一种基于逆坐标系的极区传递对准方法 | |
CN102707080A (zh) | 一种星敏感器模拟捷联惯导陀螺的方法 | |
CN105045298B (zh) | 一种基于惯导系统量测滞后的动中通天线跟踪控制方法 | |
CN102607560A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |