技术领域
[0001] 本
发明属于空间
微波遥感技术领域，涉及一种确定卫星波束与地球
位置关系的方法。
背景技术
[0002] 星载微波散射计是一种测量目标后向散射系数σ°的典型雷达系统，主要用于海面
风场测量，它是目前能够同时测量海面风速和风向的重要遥感仪器之一。由σ°反演海面风速和风向时，通过确定卫星波束与地球交点的坐标可以得到σ°对应的测量点经纬度，带入反演模型后可以计算出海面的风速和风向。
[0003] 目前，国外主要采用矢量相乘的方法建立卫星轨道
坐标系和地球坐标系的转换矩阵，进而计算卫星波束与地球的交点坐标。但是相关文献中仅给出了矢量从卫星轨道坐标系至地球坐标系的转换结果，没有详细的转换过程和说明，并且该方法没有说明位置矢量差的计算过程，因此无法在实际计算中应用。
发明内容
[0004] 本发明的技术解决问题是：克服
现有技术的不足，提供了一种卫星波束与地球交点坐标的确定方法，在卫星数据反演时能提高反演数据的可靠度。
[0005] 本发明的技术解决方案是：一种卫星波束与地球交点坐标的确定方法，步骤如下：
[0006] (1)求取 卫星的 一个波 束在卫 星本体 坐标系 中的 单位矢 量rfs，其中φs为卫星的视
角， 为卫星的方位角；所述的卫星本体坐标系用oxsyszs表示，原点o为卫星中心，xs轴指向卫星星体纵轴，zs轴指向卫星运动时面向地球的方向，ys轴按右手
正交定则定义；
[0007] (2)利用卫星的
偏航角ψ、
俯仰角θ、
横滚角γ，计算得到卫星本体坐标系与卫星轨道坐标系的转换矩阵
[0008]
[0009] 所述的卫星轨道坐标系用oxayaza表示，原点o为卫星中心，xa、za轴在卫星轨道平面内，xa轴指向卫星运动方向，za轴指向地心，ya轴按右手正交定则定义；
[0010] (3)利用卫星的经度λ和纬度l，计算得到地理坐标系与地球坐标系的转换矩阵[0011]
[0012] 所述的地理坐标系用oxtytzt表示，原点o为卫星中心，xt、yt轴在地理
水平面内，xt轴指向东，yt轴指向北，zt轴在地球质心与卫星的连线上，按右手正交定则定义；所述的地球坐标系用oexeyeze表示，原点oe为地球质心，ze轴指向地球北极方向，xe轴指向零度子午面和协议地球极赤道的交点，ye轴和ze、xe轴构成右手坐标系；
[0013] (4)利用卫星所在子午面和卫星瞬时轨道面的夹角β以及卫星运动方向与地理水平面的夹角α，计算得到辅助坐标系与地理坐标系的转换矩阵
[0014]
[0015] 所述的辅助坐标系用oxbybzb表示，原点o为卫星中心，yb、zb轴在卫星轨道平面内，yb轴指向卫星运动方向，zb轴与地球质心至卫星连线的方向相同，xb轴按右手正交定则定义；
[0016] (5)通过坐标系转换，将单位矢量rfs转换到地球坐标系中，得到方向矢量rfe，其中
[0017] (6)根据卫星波束起始点在地球坐标系中的坐标[xse，yse，zse]T和波束入射方向矢量rfe，得到卫星波束入射的直线参数方程，该直线参数与地球椭球方程联立，计算得到两组T坐标[xje，yje，zje]，
[0018]
[0019] 式中u为直线参数方程所表示直线的斜率，式中参数a为地球赤道半径、b为地球极半径，两组坐标中离卫星较近的一组坐标即为卫星波束与地球交点的坐标。
[0020] 本发明与现有技术相比的优点在于：本发明方法以一个辅助坐标系为“
桥梁”，通过坐标系转换的方法，将卫星相关坐标系(卫星本体坐标系、卫星轨道坐标系)与导航相关坐标系(地理坐标系、地球坐标系)关联起来，建立卫星本体坐标系与地球坐标系的转换矩阵，最终完成交点坐标的计算。该方法过程简便，易于工程实现。采用本发明方法可以方便的获得卫星波束与地球交点的坐标，从而在卫星数据反演时，提高反演数据的可靠度，使得海面风速和风向的确定更加准确可靠。
附图说明
[0022] 图2为本发明方法中辅助坐标系示意图；
[0023] 图3为本发明方法中辅助坐标系与地理坐标系的关系示意图。
具体实施方式
[0024] 如图1所示，本发明方法以一个辅助坐标系为“桥梁”，将卫星相关坐标系(卫星本体坐标系、卫星轨道坐标系)与导航相关坐标系(地理坐标系、地球坐标系)关联起来，建立卫星本体坐标系与地球坐标系的转换矩阵，最终完成交点坐标的计算。
[0025] 下面先给出本发明步骤中涉及到的坐标系的定义及其之间的转换矩阵。
[0026] 卫星本体坐标系
[0027] 用oxsyszs表示，原点o为卫星中心，xs轴指向卫星星体纵轴，与卫星在轨运动方向相同，zs轴指向卫星运动时面向地球的方向，ys轴按右手正交定则定义。
[0028] 卫星轨道坐标系
[0029] 用oxayaza表示，原点o为卫星中心，xa、za轴在卫星轨道平面内，xa轴指向卫星运动方向，za轴指向地心，ya轴按右手正交定则定义。当卫星的
姿态角为0°时，卫星本体坐标系与卫星轨道坐标系重合。
[0030] 地理坐标系
[0031] 用oxtytzt表示，原点o为卫星中心，xt、yt轴在地理水平面内，xt轴指向东，yt轴指向北，zt轴在地球质心与卫星的连线上，按右手正交定则定义。
[0032] 地球坐标系
[0033] 用oexeyeze表示，原点oe为地球质心，ze轴指向地球北极方向，xe轴指向零度子午面和协议地球极赤道的交点，ye轴和ze、xe轴构成右手坐标系。
[0034] 辅助坐标系
[0035] 用oxbybzb表示，原点o为卫星中心，yb、zb轴在卫星轨道平面内，yb轴指向卫星运动方向，zb轴与地球质心至卫星连线的方向相同，xb轴按右手正交定则定义，如图2所示。该坐标系在卫星轨道面内，
原型来自地理坐标系，当瞬时轨道倾角为90°，且卫星运动方向水平时，辅助坐标系与地理坐标系重合。采用辅助坐标系作为卫星轨道坐标系和地理坐标系的中间坐标系，使卫星轨道坐标系至地理坐标系的转换矩阵建立过程清晰，便于掌握。
[0036] 卫星本体坐标系与卫星轨道坐标系的转换矩阵(也称姿态矩阵)
[0037] 由卫星的姿态信息(偏航角ψ、俯仰角θ、横滚角γ)可以得到卫星本体坐标系与卫星轨道坐标系的转换矩阵：
[0038]
[0039] 地理坐标系与地球坐标系的转换矩阵
[0040] 由卫星的经度λ和纬度l可以得到地理坐标系与地球坐标系的转换矩阵：
[0041]
[0042] 式(1)、式(2)的计算结果见文献：陈哲.捷联惯导系统原理.北京：宇航出版社。
[0043] 卫星波束在卫星本体坐标系下的方向矢量rfs，经过姿态矩阵 转换后得到卫星波束在卫星轨道坐标系下的方向矢量rfa。而方向矢量rfa经过卫星轨道坐标系与地球坐标系的转换矩阵(转换矩阵 )转换后，可以得到卫星波束在地球坐标系下的方向矢量rfe。由方向矢量rfe和卫星位置可以得到波束直线方程，该波束直线方程与地球椭球方程联立求解，即可得到卫星波束与地球交点的坐标。
[0044] 由此可以看出，转换矩阵 的求取是一个关键的环节。转换矩阵 为卫星轨道坐标系(a系)与地球坐标系(e系)之间的转换矩阵，计算过程为：卫星轨道坐标系→辅助坐标系→地理坐标系→地球坐标系。
[0045] (1)根据上面坐标系的定义部分可知，卫星轨道坐标系和辅助坐标系的关系为：xb＝ya，yb＝xa，zb＝-za，由此可以得到这两个坐标系之间的转换矩阵为：
[0046]
[0047] (2)辅助坐标系与地理坐标系的关系示意图如图3所示。两个坐标系通过两个参数关联起来：卫星所在子午面和瞬时轨道面的夹角β，以及卫星运动方向与地理水平面的夹角α。
[0048] 如图3所示，S为当前时刻GPS确定的卫星位置[xse，yse，zse]T，So为S之前的最近时间由GPS确定的卫星位置。由于地球在自转，在惯性空间里，地球坐标系以24小时的周期绕ze轴匀速转动，所以，So时刻至S时刻地球坐标系绕ze轴转动了ωet的角度。定义So时刻的地球坐标系为e′系，S时刻的地球坐标系为e系，e′系绕ze′轴转ωet后与e系重T合。So时刻卫星位置表示为[xs0e′，ys0e′，zs0e′]，则So时刻卫星位置在e系中表示为：
[0049]
[0050] 其中ωe为地球自转
角速度，t为So与S点所对应的GPS时间差。则矢量SoS为：
[0051]
[0052] 将矢量 转换为地理坐标系下的矢量
[0053]
[0054] 矢量 与oxtyt平面的夹角α为：
[0055]
[0056] 矢量 在oxtyt平面的投影与yb轴在同一个方向上，该投影与yt轴的夹角为β。则β角为：
[0057]
[0058] 辅助坐标系绕xb轴转α，绕z′a(辅助坐标系绕xb轴转α角后的中间坐标系的z轴)转β，将会与地理坐标系重合。则辅助坐标系至地理坐标系的转换矩阵 为：
[0059]
[0060] 假设卫星轨道坐标系(a系)与地球坐标系(e系)之间的转换矩阵 那么[0061]
[0062] 卫星一个波束(视角φs、方位角 )在卫星本体坐标系中的单位矢量为：
[0063]
[0064] 该波束在地球坐标系中的方向矢量rfe为：
[0065]
[0066] 已知波束起始点(卫星位置)在地球坐标系中的坐标S＝[xse，yse，zse]T和波束入射方向rfe，可以确定波束入射的直线参数方程为：
[0067] B(u)＝S+rfe·u (13)
[0068] 式中，u为参数。该直线与地球椭球面有两个交点，通过与椭球方程联立可以计算出两个坐标：
[0069]
[0070] 式中参数a为地球赤道半径、b地球极半径。两个坐标中，离卫星较近的即为卫星T波束与地球的交点坐标[xje，yje，zje]。
[0071] 本发明
说明书中未作详细描述的内容属本领域技术人员的公知技术。
