发明名称 一种地球同步轨道卫星微波图像分辨率的增强方法
摘要 本发明公开一种地球同步轨道卫星微波图像分辨率的增强方法,涉及微波遥感领域,包括以下步骤:确定卫星微波遥感系统参数;设置目标区域参数;通过过采样获得地球同步轨道卫星对地观测实测数据,利用重采样获得完整的微波遥感图像;使用投影球方法处理微波遥感图像,去除采样过程中因地球曲率引入的图像信息误差;对坐标投影后的图像进行维纳滤波,提高图像空间分辨率;对维纳滤波处理后的图像进行逆傅里叶变换和坐标变换,获得分辨率增强后的图像。利用本发明可以得到分辨率更高的地球同步轨道卫星微波遥感图像,特别适用于位于远离星下点地区、探测跨度较大的图像,实现效率更优,具有更强的适应性和兼容性。
申请公布号 CN103810679A 申请公布日期 2014.05.21
申请号 CN201410056715.9 申请日期 2014.02.19
申请人 北京航空航天大学 发明人 刘大伟;刘凯;吕长春;魏川奇;苗俊刚
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 周长琪
主权项 1.一种地球同步轨道卫星微波遥感图像分辨率的增强方法,其特征在于:通过以下步骤实现:步骤一:地球同步轨道卫星微波遥感系统参数与目标区域地理位置参数;步骤二:获取微波遥感图像;a、获得在地球坐标系中以地球表面点C(x<sub>0</sub>,y<sub>0</sub>)为瞬时视场中心的采样矩阵;使用高斯函数近似采样矩阵<img file="FDA0000467343550000011.GIF" wi="95" he="78" />并利用空间分辨率:<maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>G</mi><mo>&OverBar;</mo></mover><mi>C</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msup><mrow><mn>2</mn><mi>&pi;D</mi></mrow><mn>2</mn></msup></mfrac><mi>exp</mi><mo>{</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>[</mo><mfrac><msup><mrow><mo>(</mo><msub><mrow><mi>x</mi><mo>-</mo><mi>x</mi></mrow><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><msup><mi>D</mi><mn>2</mn></msup></mfrac><mo>+</mo><mfrac><msup><mrow><mo>(</mo><msub><mrow><mi>y</mi><mo>-</mo><mi>y</mi></mrow><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><msup><mi>D</mi><mn>2</mn></msup></mfrac><mo>]</mo><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(1)中,D为空间分辨率;x、y为瞬时视场内的点在地球坐标系中的坐标;x<sub>0</sub>、y<sub>0</sub>为点C在地球坐标系中的坐标;b、利用采样矩阵<img file="FDA0000467343550000013.GIF" wi="68" he="73" />对探测地区进行采样;通过采样矩阵与地面相应位置地球表面亮温乘积的积分来描述搭载在卫星上的天线探测地球表面亮温的过程,即:<maths num="0002"><![CDATA[<math><mrow><msub><mi>t</mi><mi>A</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><msub><mrow><mo>&Integral;</mo><mo>&Integral;</mo></mrow><mi>C</mi></msub><msub><mover><mi>G</mi><mo>&OverBar;</mo></mover><mi>C</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><msub><mi>t</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mi>dxdy</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,t<sub>A</sub>为采样后的地表亮温,t<sub>B</sub>为地表真实亮温;将式(2)等价表示成<img file="FDA0000467343550000015.GIF" wi="197" he="77" />与t<sub>B</sub>(x,y)的卷积,即:<maths num="0003"><![CDATA[<math><mrow><msub><mi>t</mi><mi>A</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><msub><mover><mi>G</mi><mo>&OverBar;</mo></mover><mi>C</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&CircleTimes;</mo><msub><mi>t</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>从而,获得地球同步轨道卫星测得的原始探测图像;c、对原始探测图像进行重采样,生成微波遥感图像;步骤三:建立投影球坐标系;投影球坐标系以卫星为球心,投影球坐标系中各点由投影球坐标系的方位角和仰角<img file="FDA0000467343550000017.GIF" wi="149" he="65" />来定义,并将卫星与地球连线中点坐标定义为(0°,0°);步骤四:将微波遥感图像投影在投影球坐标系中,得到:<img file="FDA0000467343550000018.GIF" wi="1305" he="96" />式(4)中,t<sub>A</sub>(x<sub>0</sub>,y<sub>0</sub>)<img file="FDA0000467343550000019.GIF" wi="209" he="83" />和<img file="FDA00004673435500000110.GIF" wi="174" he="81" />分别是t<sub>A</sub>(x<sub>0</sub>,y<sub>0</sub>)、<img file="FDA00004673435500000111.GIF" wi="196" he="83" />和t<sub>B</sub>(x,y)投影到投影球坐标系后的表示形式;将式(4)等价为<img file="FDA00004673435500000112.GIF" wi="190" he="78" />与<img file="FDA00004673435500000113.GIF" wi="174" he="82" />的卷积,得到:<img file="FDA00004673435500000114.GIF" wi="1288" he="82" />式(5)中,<img file="FDA00004673435500000115.GIF" wi="166" he="76" />为探测噪声;步骤四:对投影后的微波遥感图像进行去卷积处理,重建得到频域中的高分辨率微波遥感图像;步骤三中式(5)进行傅里叶变换,在频域得到:T<sub>A</sub>(u,v)=H(u,v)T<sub>B</sub>(u,v)+N(u,v)  (6)其中,T<sub>A</sub>(u,v),H(u,v),T<sub>B</sub>(u,v)和N(u,v)分别是<img file="FDA00004673435500000117.GIF" wi="639" he="83" />和<img file="FDA0000467343550000021.GIF" wi="162" he="81" />的在频域中对应的表示形式;根据式(7)得到维纳滤波算子W(u,v)为:<maths num="0004"><![CDATA[<math><mrow><mi>W</mi><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>H</mi><msup><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><mo>*</mo></msup></mrow><mrow><msup><mrow><mo>|</mo><mi>H</mi><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><mfrac><msup><mi>NE&Delta;T</mi><mn>2</mn></msup><mrow><msub><mi>P</mi><mi>TATA</mi></msub><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow></mrow></mfrac></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(8)中,NEΔT<sup>2</sup>为辐射计噪声;P<sub>TATA</sub>(u,v)为频域中微波遥感图像T<sub>A</sub>(u,v)的功率谱密度。通过W(u,v)重建得到频域中的高分辨率微波遥感图像:T<sub>C</sub>(u,v)=W(u,v)T<sub>A</sub>(u,v)   (9)步骤五:对频域中的高分辨率微波遥感图像变换到投影球坐标系,再将高分辨率地球表面亮温图像变换到地球表面的地球坐标系,得到最终的分辨率增强结果。
地址 100191 北京市海淀区学院路37号