一种星载微波辐射计天线误差校正方法与流程

本发明涉及一种星载微波辐射计天线误差校正方法，属于微波遥感领域。

背景技术：

星载微波辐射计通过被动接收海面微波辐射信号对海洋环境参数如海面温度、风速、大气液态水含量和大气水气含量等进行遥感测量。星载微波辐射计通过天线周期性旋转实现海洋目标的连续扫描观测，在每个扫描周期中，微波辐射计依次对地球场景、高温定标源和低温定标源进行观测。高温定标源利用受到精密温控处于恒定温度的微波吸收体实现，低温定标源由波束指向宇宙冷空背景的天线实现。在每个扫描周期中通过观测两个定标源并结合地面测试参数可确定系统定标参数，进而将地球场景观测计数值转换为亮温。

在星载微波辐射计观测地球场景和冷空时，天线实际接收亮温为天线方向图在整个4π空间区域内对空间分布亮温tb(θ,φ)积分，即

其中tb(θ,φ)为空间任意方向(θ,φ)辐射亮温，α和ω分别表示天线波束观测角和方位角，fn(β,γ)表示在波束指向(α,ω)方向时空间任意点(θ,φ)对应的天线方向图。

根据星载微波辐射计观测几何特性和空间辐射亮温分布特性，可将整个4π空间区域划分为地球区域和冷空区域两部分，天线实际接收亮温表达式可写为：

为天线波束效率归一化因子，ωe和ωcs分别表示地球区域和冷空区域。

按照不同的积分区域，可将上式表示为波束效率因子对辐射亮温加权求和形式：

fe(α,ω)+fcs(α,ω)＝1(3-6)

表示地球平均辐射亮温，tc0表示冷空标称辐射亮温，fe(α,ω)和fcs(α,ω)分别表示在天线波束指向(α,ω)方向时地球区域和冷空区域在天线方向图中所占波束效率因子。

从式(3-3)可以看出，当星载微波辐射计在进行对地观测时，受到天线波束效率因子影响，fe(α,ω)<1，且导致天线实际接收亮温ta小于地球场景真实辐射亮温当星载微波辐射计在每个扫描周期内进行冷空定标观测时，同样因为受到天线波束效率因子影响，fcs(α,ω)<1，且导致天线实际接收亮温大于冷空标称辐射亮温，即ta>tc0，若直接采用冷空标称辐射亮温tc0进行定标，将导致定标出的亮温小于场景真实辐射亮温，带来测量误差。综上所述，在星载微波辐射计进行对地观测和冷空定标观测时有限的波束效率导致地球观测亮温小于场景真实辐射亮温，需要结合星载微波辐射计天线方向图数据和观测几何参数对仪器亮温观测数据进行校正，提高仪器测量精度。

技术实现要素：

本发明解决的技术问题是：克服现有技术的不足，提供了一种星载微波辐射计天线误差校正方法，能够校正有限的天线波束效率引入的定标误差和地球亮温观测误差。本发明通过在卫星轨道坐标系下将空间区域划分为地球区域和冷空区域两部分，在考虑卫星姿态因素影响的情况下，利用旋转矩阵法得到卫星轨道坐标系下空间任意点在天线坐标系下坐标，然后分别在地球区域和冷空区域进行积分得到对应于地球观测和冷空观测的天线波束效率因子，分别对地球观测亮温和冷空定标亮温进行修正，提高了仪器亮温测量精度。

本发明的技术解决方案是：一种星载微波辐射计天线误差校正方法，包括如下步骤：

s1，在卫星轨道坐标系下将整个空间区域划分为冷空区域和地球区域，分别计算当天线波束指向地球区域和冷空区域进行观测时的地球波束效率因子和冷空波束效率因子；

s2，利用s1中所述冷空波束效率因子对冷空定标亮温进行校正，得到冷空定标亮温校正值；

s3，利用s2中所述冷空定标亮温校正值和空间参数对各个天线波束指向下的地球场景观测计数值进行定标，得到地球观测亮温定标值；所述空间参数包括热源观测计数值、冷空观测计数值、热源亮温和通道非线性系数；

s4，利用s1中所述地球波束效率因子对s3中所述地球观测亮温定标值进行校正，得到校正后的地球观测亮温。

进一步地，所述校正后的地球观测亮温为其中，αe和ωe分别为在卫星本体坐标系下天线波束指向地球时的观测角和方位角，ta(αe,ωe)为地球观测亮温定标值，tc0为冷空标称辐射亮温，fee(αe,ωe)和fecs(αe,ωe)分别为天线波束指向地球区域时的地球波束效率因子和冷空波束效率因子，且满足fee(αe,ωe)+fecs(αe,ωe)＝1。

进一步地，所述地球观测亮温定标值为其中，cs(αe,ωe)为地球场景观测计数值，和为系数，tccor为冷空定标亮温校正值，ch为热源观测计数值，cc为冷空观测计数值，th为热源亮温，u为通道非线性系数。

进一步地，所述冷空定标亮温校正值为其中，为地球平均辐射亮温，tc0为冷空标称辐射亮温，fce和fccs分别为天线波束指向冷空区域时的地球波束效率因子和冷空波束效率因子，且满足fce+fccs＝1。

进一步地，所述计算冷空波束效率因子的方法为：在天线波束指向冷空区域进行观测时，在天线坐标系下按冷空区域对天线方向图进行积分，得到所述冷空波束效率因子为0≤θ≤θe,和θe<θ≤π,且fce+fccs＝1；其中，fn(β,γ)为归一化天线方向图，n为天线方向图积分归一化因子，re为地球平均半径，he为卫星轨道高度。

进一步地，所述计算地球波束效率因子的方法为：在天线波束指向地球区域进行观测时，在天线坐标系下按地球区域对天线方向图进行积分，得到所述地球波束效率因子为0≤θ≤θe,和θe<θ≤π,且fee(αe,ωe)+fecs(αe,ωe)＝1；其中，fn(β,γ)为归一化天线方向图，n为天线方向图积分归一化因子，re为地球平均半径，he为卫星轨道高度。

优选地，将所述卫星轨道坐标系转换为所述天线坐标系的方法为：取卫星轨道坐标系任意点则该点坐标转换为天线坐标系坐标为β和γ分别为在天线坐标系中的俯仰角和方位角，x＝sinθcosθ，z＝cosθ。

本发明与现有技术相比的优点在于：

(1)本发明采用旋转矩阵法实现卫星轨道坐标系下空间任意点坐标至天线坐标系中坐标转换，直接利用根据星载微波辐射计天线波束观测角和方位角计算出的坐标转换矩阵进行卫星轨道坐标系到天线坐标系的转换。相比于传统采用建立辅助坐标系和在空间多次投影实现坐标转换的方法相比，本发明采用的坐标转换方法更为简单和易于实现。

(2)本发明在卫星轨道坐标系到天线坐标系转换过程中考虑了星载微波辐射计在轨扫描观测过程中卫星姿态角(俯仰角、偏航角和横滚角)变化影响，修正了由卫星姿态角变化带来的地球和冷空波束效率因子计算误差，提高了星载微波辐射计定标精度和地球观测亮温修正精度，可进一步提升星载微波辐射计数据应用效果。

附图说明

图1为微波辐射计定标原理框图；

图2为本发明方法实施流程框图。

具体实施方式

如图1所示为微波辐射计定标原理框图。

本发明方法星载微波辐射计天线误差校正步骤如下：

(1)波束效率因子计算

在卫星轨道坐标系下将整个空间区域划分为冷空和地球两部分，在天线波束指向不同区域进行观测时(地球观测和冷空定标观测)，可按照冷空和地球两部分区域对天线方向图进行积分，得到不同波束指向位置处地球和冷空区域波束效率因子。设卫星轨道高度为he，地球平均半径为re，则地球与冷空区域边界对应的天底角为：

在卫星轨道坐标系下地球区域和冷空区域可表示为表4-1中所示角度范围：

表4-1地球区域和冷空区域角度范围

当星载微波辐射计天线波束指向冷空进行冷空定标观测时，设天线波束主轴za在卫星本体坐标下与卫星本体坐标系zs轴夹角为αcs，za在xsosys平面内投影与xs夹角为ωcs，在进行冷空定标观测时αcs和ωcs通常为固定角度，则校正系数为：

fce+fccs＝1(4-4)

当星载微波辐射计波束指向地球观测地球场景目标时，设波束主轴za在卫星本体坐标下与zs轴夹角为αe，za在xsosys平面内投影与xs夹角为ωe，则校正系数为：

fee(αe,ωe)+fecs(αe,ωe)＝1(4-7)

由于空间任意点(θ,φ)与天线方向图fn(β,γ)存在一一对应关系，因此在星载微波辐射计天线波束指向空间任意方向时均能计算得到在该波束指向下地球区域和冷空区域波束效率因子，进而通过对空间各区域辐射亮温值加权求和得到天线实际接收亮温。

卫星轨道坐标系下空间任意一点在特定波束指向(指向地球和冷空)下对应在天线坐标系下的角度(β,γ)可按照如下步骤计算：

a)定义卫星轨道坐标系ooxoyozo，坐标系原点为卫星中心，xo、zo轴在卫星轨道平面内，zo轴指向地心，xo轴指向卫星运动方向，yo轴按右手正交定则定义。在卫星轨道坐标系下空间任意点方向表示为则沿着从卫星轨道坐标系原点到该点方向的单位矢量为：

b)定义卫星本体坐标系osxsyszs，原点为卫星中心，xs轴指向卫星星体纵轴，zs轴位于第一象限，指向卫星运动时面向地球的方向，ys轴按右手正交定则定义。设卫星偏航、横滚和俯仰三个姿态角度分别为yaw、roll和pitch，可以得到卫星本体坐标系到卫星轨道坐标系转换矩阵：

mos＝mpitch·(mroll·myaw)(4-9)

在卫星本体坐标系下矢量为：

c)定义天线坐标系oaxayaza，原点oa为卫星中心，za为波束主轴方向，设从oa到地心的单位长度矢量为则xa轴满足右手法则。设波束主轴za在卫星本体坐标下与zs轴夹角为α，za在xsosys平面内投影与xs夹角为ω，则从卫星本体坐标系到天线坐标系转换矩阵为：

d)根据卫星本体坐标系和天线坐标系转换关系可将转换至天线坐标系：

根据坐标可计算出空间任意点(θ,φ)对应在天线坐标系下的角度：

根据β、γ和天线方向图数据采用查找表方法可得到空间任意点对应的天线方向图数据fn(β,γ)，最后可计算出不同波束指向下地球区域和冷空区域对应的波束效率因子。

(2)冷空定标亮温校正

根据式(4-2)和(4-3)在计算得到天线波束效率因子fce和fccs后，对冷空定标亮温进行校正：

为地球平均辐射亮温，tc0为冷空标称辐射亮温。

(3)地球观测亮温定标

利用冷空定标亮温校正值tccor﹑热源观测计数值ch、冷空观测计数值cc﹑热源亮温th和通道非线性系数u对在各个天线波束指向下得到的地球场景观测计数值cs(αe,ωe)进行定标，得到地球观测亮温定标值ta(αe,ωe)。

(4)地球观测亮温校正

根据式(4-5)和式(4-6)计算出的校正系数可对(3)步骤中地球观测亮温定标值进行校正，校正在进行地球场景观测时冷空辐射带来的误差，校正后的地球观测亮温为：

下面结合附图进一步说明。

如图2所示，本专利提出的冷空定标亮温修正方法实施过程。

(1)在卫星轨道坐标系下将整个空间区域划分为地球区域和冷空区域两部分，采用旋转矩阵法得到卫星轨道坐标系下空间任意点坐标在天线坐标系中对应坐标，当天线波束指向地球区域时，分别在地球区域和冷空区域对天线方向图进行积分得到用于校正地球场景观测亮温误差的天线波束效率因子。当天线波束指向冷空区域时，分别在地球和冷空区域对天线方向图进行积分得到用于校正冷空定标亮温的天线波束效率因子。

(2)冷空定标亮温校正

利用冷空定标观测时计算出的天线波束效率因子对冷空定标亮温进行校正，校正星载微波辐射计在进行冷空定标观测时地球场景辐射对冷空定标观测带来的误差。

(3)地球场景观测亮温定标

利用冷空定标亮温修正值、热源温度、冷空定标观测计数值、热源观测计数值和通道非线性系数对地球场景观测计数值进行定标，得到地球场景亮温定标值。

(4)利用天线波束指向地球区域时计算出的天线波束效率因子对(3)中的地球场景亮温定标值进行校正，校正星载微波辐射计观测地球场景时天线有限波束效率带来的亮温测量误差。

本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。