一种适用于被动微波资料的广义RI指数频率干扰检测及订正方法技术

本发明专利技术公开一种适用于被动微波资料的广义RI指数频率干扰检测及订正方法，属于大气微波遥感技术领域。包括数据读取、地球表面分类、数据检测、干扰区域订正、信息分析和显示；数据读取对星载微波辐射计观测数据进行读取并进行通道分离处理，地球表面分类是对全球的微波亮温数据按无积雪覆盖陆地表面、冰雪覆盖表面和洋面分类，数据检测是对分组的微波亮温数据计算广义RI指数并据其大小检测干扰信号的分布；干扰区域订正计算受干扰通道的预期亮温值，定量给出订正后亮温的分布；信息分析和显示对上述各部分的输出信息进行判断、标记、处理、显示和存储。具有用一个通用检测方法检测不同地球表面类型，方法简单，计算量小，适用性好等优点。

A Generalized RI Exponential Frequency Interference Detection and Correction Method for Passive Microwave Data

The invention discloses a generalized RI index frequency interference detection and correction method suitable for passive microwave data, which belongs to the field of atmospheric microwave remote sensing technology. It includes data reading, earth surface classification, data detection, interference region correction, information analysis and display; data reading reads the observation data of spaceborne microwave radiometer and carries out channel separation processing. The earth surface classification is to classify the global microwave brightness temperature data according to the land surface, ice-snow surface and ocean surface without snow cover, and data detection is to group microwave brightness. Temperature data calculate the generalized RI index and detect the distribution of interference signal according to its size; interference area correction calculates the expected brightness temperature value of the interference channel and gives the corrected brightness temperature distribution quantitatively; information analysis and display judges, marks, processes, displays and stores the output information of the above parts. It has the advantages of simple method, small calculation and good applicability to detect different types of earth surface with a general detection method.

【技术实现步骤摘要】
一种适用于被动微波资料的广义RI指数频率干扰检测及订正方法
本专利技术涉及一种适用于被动微波资料的广义RI指数频率干扰检测及订正方法，属于大气微波遥感

技术介绍
微波观测资料可以提供各种天气条件下的陆表和大气信息，星载被动微波传感器在数值天气预报模式及资料同化中占有愈来愈重要的地位。然而，低频微波通道被各种主动和被动遥感探测占据，例如通信卫星、天气和军用雷达、GPS信号，移动电话等，这使得卫星微波辐射计接收的信息除了来自真实地表的热辐射信号，还叠加了来自主动传感器的信号或表面反射的辐射信号，称为无线电频率干扰。由于强的无线电频率干扰信号很容易淹没相对较弱的地表的热发射辐射，无线电频率干扰会使局地观测亮温高于正常值范围，从而对卫星观测数据造成污染，进而导致较大的反演误差。因此，无线电频率干扰在主动及被动微波遥感探测领域已成为越来越严重的问题而且，未来的传感器为了实现特定的观测目标，通道频率仍需要使用这些不受保护的波段，所以很多学者对星载微波辐射计的无线电频率干扰识别做了大量的研究，提出了一系列识别方法。提高无线电频率干扰识别的准确性，对评价微波数据反演地表参数的精度，以及提高星载微波资料在陆面过程模式及资料同化中的利用率等诸多方面的研究，都具有非常重要的意义。以往的研究开发了各种各样的检测星载微波辐射计资料中无线电频率干扰信号的方法。Li等(2004)最初提出了用频谱差法来检测无线电频率干扰，确定了AMSR-E在C和X波段的观测值中无线电频率干扰信号的强度和分布范围，随后进一步提出了用主成分分析法来分析陆地区域的无线电频率干扰分布特征；Njoku等(2005)提出了统计平均值和标准差分析法，得出的研究结果是AMSR-E在6.925GHz和10.67GHz通道受无线电频率干扰影响的区域分别处于不同的地理位置；Li等(2006)接着又提出了可以用多通道回归算法来识别WindSat数据中分布于洋面上的无线电频率干扰信号；Lacava等(2016)用多时相法分析了AMSR-EC波段的陆面无线电频率干扰；Zou等(2012)提出用主成分分析法可以检测MWRI陆地表面的无线电频率干扰分布；Zhao等(2013)改进了主成分分析法，提出用双主成分分析法分析WindSat资料在积雪覆盖地区的无线电频率干扰分布；官莉等(2014)提出可用简化的主成分分析法对欧洲陆地区域的AMSR-E无线电频率干扰进行检测。而对于多样的地表覆盖状况，不同检测方法的效果对不同的地表状况有局限性，有的甚至完全不适用于另一种地表状况。因此，需要一种适用于被动微波资料频率干扰的广义RI指数检测方法，从而在晴空/有云大气状况下，针对多种地球表面类型，对多种星载微波传感器资料中的无线电频率干扰进行检测。
技术实现思路
本专利技术要解决的技术问题是针对现有技术只能检测某一种传感器或某一两种表面类型的无线电频率干扰，对不同的地表状况有局限性的缺点，提出一种适用于被动微波资料的广义RI指数频率干扰检测及订正方法，通过计算广义RI指数，获得地球表面的无线电频率干扰信号的分布和强度，并订正受干扰微波观测资料，克服检测的局限性。为解决上述技术问题，本专利技术提供一种适用于被动微波资料的广义RI指数频率干扰检测及订正方法，包括数据读取、地球表面分类、数据检测、干扰区域订正、信息分析和显示五个步骤；数据读取负责对星载微波辐射计观测数据进行读取，并进行通道分离处理；地球表面分类是对全球的微波亮温数据进行分组，按无积雪覆盖陆地表面、冰雪覆盖表面和海洋表面分类；数据检测负责对分组的微波亮温数据计算广义RI指数，根据RI指数的数值大小检测出地球表面的无线电频率干扰信号的分布，并确定干扰信号的强度；干扰区域订正负责计算受干扰通道的预期亮温值，定量给出订正后亮温的分布；信息分析和显示负责对上述各部分的输出信息进行判断、标记、处理、显示和存储。具体步骤包括：(1)数据读取：从星载微波辐射计观测的HDF数据格式文件中读取所有通道的亮温观测值，进行通道分离后，定量显示每个通道观测的亮温空间分布；分别判断6.9GHz水平和垂直极化通道、7.3GHz水平和垂直极化通道、10.7GHz水平和垂直极化通道、18.7GHz水平和垂直极化通道的观测亮温数值是否存在异常高值的视场，若没有则初步判断该通道观测没有被干扰，若有异常高值区则初步判断该通道观测有可能被污染；(2)地球表面分类：对全球的微波亮温数据进行分组，首先用海陆分布掩码，如LandSeaMask数据，对微波辐射计亮温数据进行分组，分成陆地表面和海洋表面两大类；再用以下两个有散射敏感性的谱差指数(Zouetal，2012)来识别冰或雪表面，把陆地表面分成无积雪覆盖陆地表面和冰雪覆盖陆地表面，把海洋表面分成无冰雪覆盖海洋表面和冰雪覆盖海洋表面：TB89H-TB18H≥10K，TB89V-TB18V≥10K；式中，H代表水平极化，V代表垂直极化，TB89H、TB18H分别代表89.0GHz和18.7GHz水平极化通道的亮温值，TB89V、TB18V分别代表89.0GHz和18.7GHz垂直极化通道的亮温值，满足该式条件的为冰雪覆盖陆地表面或海洋表面；然后，针对不同地球表面类型(冰雪覆盖地表和无冰雪陆面、洋面和冰雪表面)，把全球微波亮温观测值分组；(3)数据检测：首先，通过下式计算每组中微波亮温值的散射指数SI，将SI数值大于10K的相应像元判定为降水像元予以剔除，即剔除每组中包含自然云、降水等气象要素的像素点：式中，V表示垂直极化，系数a、b、c和d由无云情况下的SSM/I全球观测数据集导出，Tb19v、Tb22v、Tb85v分别代表19GHz、22GHz、85GHz垂直极化通道的微波亮温值；然后，分别在无积雪覆盖陆地表面、冰雪覆盖陆地表面、无冰雪覆盖海洋表面和冰雪覆盖海洋表面的每一组数据范围内，计算所有其他通道相对于感兴趣通道的经验拟合系数；再按下式计算给定通道与其他辐射计通道的预期亮温值的偏差，并将该偏差定义为广义RI指数ΔTb[i]：式中,i代表受到无线电频率干扰的通道的通道号，j代表除i通道以外所有通道的通道号，a0[i]是常数项，aj[i]是应用于每个通道j以计算i通道亮温值的线性拟合系数；由于发现二次项不会显著改善标准偏差，因而可以采用线性拟合方式，将广义RI指数ΔTb[i]写为上式；式中其他通道的线性组合倾向于检测出感兴趣的信道的亮温受到的无线电频率干扰，目前的星载微波辐射计资料的无线电频率干扰主要出现在C波段和X波段通道，以Kelvins为单位表示的ΔTb[i]的大的正值通常对应于来自无线电频率干扰的非自然发射；最后，根据广义RI指数的数值大小，量化频率干扰信号的强度，并检测出地球表面的无线电频率干扰信号的分布；(4)干扰区域订正：利用自然地表发射辐射在各通道观测之间的高度相关性，通过线性拟合关系订正受到无线电频率干扰的通道的观测值，通过下式计算受干扰通道的预期亮温值Tb[i]，对i通道中受无线电频率干扰区域进行订正，定量给出订正后亮温的分布；式中,i代表受到无线电频率干扰的通道的通道号，j代表所有通道的通道号，a0[i]是常数项，aj[i]是应用于每个通道j以计算i通道亮温值的线性拟合系数，对应于感兴趣的信道的系数本文档来自技高网...

【技术保护点】
1.一种适用于被动微波资料的广义RI指数频率干扰检测及订正方法，其特征在于：包括数据读取、地球表面分类、数据检测、干扰区域订正、信息分析和显示五个步骤；数据读取是对星载微波辐射计观测数据进行读取，并进行通道分离处理；地球表面分类是对全球的微波亮温数据进行分组，按无积雪覆盖陆地表面、冰雪覆盖表面和海洋表面分类；数据检测是对分组的微波亮温数据计算广义RI指数，根据RI指数的数值大小检测出地球表面的无线电频率干扰信号的分布，并确定干扰信号的强度；干扰区域订正是计算受干扰通道的预期亮温值，定量给出订正后亮温的分布；信息分析和显示是对上述各部分的输出信息进行判断、标记、处理、显示和存储。

【技术特征摘要】
1.一种适用于被动微波资料的广义RI指数频率干扰检测及订正方法，其特征在于：包括数据读取、地球表面分类、数据检测、干扰区域订正、信息分析和显示五个步骤；数据读取是对星载微波辐射计观测数据进行读取，并进行通道分离处理；地球表面分类是对全球的微波亮温数据进行分组，按无积雪覆盖陆地表面、冰雪覆盖表面和海洋表面分类；数据检测是对分组的微波亮温数据计算广义RI指数，根据RI指数的数值大小检测出地球表面的无线电频率干扰信号的分布，并确定干扰信号的强度；干扰区域订正是计算受干扰通道的预期亮温值，定量给出订正后亮温的分布；信息分析和显示是对上述各部分的输出信息进行判断、标记、处理、显示和存储。2.根据权利要求1所述的适用于被动微波资料的广义RI指数频率干扰检测及订正方法，其特征在于：所述检测及订正方法的具体步骤包括：(1)数据读取：从星载微波辐射计观测的HDF数据格式文件中读取所有通道的亮温观测值，进行通道分离后，定量显示每个通道观测的亮温空间分布；分别判断6.9GHz水平和垂直极化通道、7.3GHz水平和垂直极化通道、10.7GHz水平和垂直极化通道、18.7GHz水平和垂直极化通道的观测亮温数值是否存在异常高值的视场，若没有则初步判断该通道观测没有被干扰，若有异常高值区则初步判断该通道观测有可能被污染；(2)地球表面分类：对全球的微波亮温数据进行分组，首先用海陆分布掩码对微波辐射计亮温数据进行分组，分成陆地表面和海洋表面两大类；再用以下两个有散射敏感性的谱差指数来识别冰或雪表面，把陆地表面分成无积雪覆盖陆地表面和冰雪覆盖陆地表面，把海洋面分成无冰雪覆盖海洋表面和冰雪覆盖海洋表面：丁B89H-TB18H≥10K，TB89V-TB18V≥10K；式中，H代表水平极化，V代表垂直极化，TB89H、TB18H分别代表89.0GHz和18.7GHz水平极化通道的亮温值，TB89V、TB18V分别代表89.0GHz和18.7GHz垂直极化通道的亮温值，满足该式条件的为冰雪覆盖陆地表面或海洋表面；然后，针对不同地球表面类型，把全球微波亮温观测值分组；(3)数据检测：首先，通过下式计算每组中微波亮温...

【专利技术属性】
技术研发人员：吴莹，钱博，王振会，
申请(专利权)人：南京信息工程大学，
类型：发明
国别省市：江苏,32