一种确定波谱仪数据降雨标志的方法
【专利摘要】本发明公开了一种确定波谱仪数据降雨标志的方法，利用TRMM卫星PR雷达9°入射角观测数据，基于9°入射角海面后向散射系数(NRCS)不随海面风、浪等动力环境因素变化的特性，通过分析无雨和降雨海面NRCS差异，确定了9°入射角海面NRCS降雨阈值，该阈值可用于确定波谱仪数据的降雨标志。确定降雨标志可用于剔除波谱仪部分受降雨干扰的无效数据，对保证波谱仪关于海面动力环境遥感产品的质量，起着非常重要的作用。同时具有简单实现的优势，适合于实际应用。
【专利说明】
一种确定波谱仪数据降雨标志的方法
技术领域
[0001]本发明属于海洋微波遥感技术领域，具体是一种确定波谱仪数据降雨标志的方法。
【背景技术】
[0002]降雨是一种常见的大气现象，海面的降雨会引起波谱仪海面NRCS的变化，由于目前的波谱仪动力环境遥感模型仅适用无雨海面，所以海面降雨会干扰波谱仪对海面风、浪等动力环境的遥感观测。海面降雨标志的确定，对于剔除波谱仪部分受降雨干扰的无效数据，保证遥感产品的质量，具有非常重要的现实意义和科研价值。
[0003]波谱仪是一种小入射角星载雷达，入射角范围在0°_11°，主要用于测量海面的风、浪等动力环境参数，全球第一颗波谱仪卫星将搭载中法海洋卫星(CFOSAT)，预计于2018年发射(Ren et al.，2015)。目前，对于小入射角雷达，只有TRMM卫星PR雷达能够实现海面降雨量的遥感。PR雷达作为一种气象雷达，是通过观测大气的反射率，利用Z-R关系(反射率-降雨量关系)来确定降雨量，进而确定降雨标志(Iguchi et al.2000；Li et al.2002；Setoand Iguchi 2007)。而对于动力环境遥感卫星，如波谱仪，散射计和合成孔径雷达(SAR)等，都是基于观测的海面NRCS来遥感海面动力环境参数，并未配置测量大气发射率的功能模块。因此，波谱仪难以应用与PR雷达相同的方法来确定波谱仪数据降雨标志。
[0004]由于波谱仪不能观测大气发射率，所以只能从海面NRCS的角度来确定降雨标志。在降雨条件下，波谱仪观测到的海面NRCS，是海面风、浪等动力环境因素和降雨的共同作用结果，而实际海面的动力环境始终处于变化之中，很难根据经验用一个预先设定的海面NRCS阈值来区分是否降雨。因此，当需要同时考虑降雨和动力环境因素对海面NRCS影响时，如何应用简单的海面NRCS降雨阈值，来确定波谱仪数据降雨标志，是一个非常困难的问题。
[0005]针对上述问题，解决的关键在于如何消除动力环境因素对海面NRCS的影响。根据小入射角海面NRCS特性研究(Hesany et al.,2000；Tran,et al.，2007)，9°-11°入射角观测的海面NRCS，基本不随海面动力环境变化，这主要是因为在该入射角范围内，海面NRCS是准镜面散射和布拉格散射的共同作用结果，准镜面散射随着粗糙度增大而减小，而布拉格散射随着粗糙度增大而增大，从而导致不同动力环境下海面NRCS整体的不变性。因此，若选取该入射角范围内观测数据来确定波谱仪数据降雨标志，则可以摆脱海面动力环境因素对NRCS的影响，保证海面NRCS变化量仅与降雨量相关，从而有利于确定海面NRCS降雨阈值。然而迄今为止，仍然没有研究将其用于确定波谱仪数据降雨标志。
【发明内容】
[0006]本发明的目的是提供了一种确定波谱仪数据降雨标志的方法。
[0007]本发明通过以下技术方案实现:
[0008]确定波谱仪数据降雨标志的方法包括以下步骤:
[0009]SOl.收集研究区域TRMM卫星PR雷达2A25标准产品数据；
[0010]S02.从步骤SOl获得的数据中选取海面观测数据，选取方法:提取产品中的“TRMM2A25Method Flag”因子，当设置值为O时，表示海面观测数据；
[0011]S03.从海面观测数据中分别选取无雨和降雨海面数据，选取方法:提取产品中的“Rain Flag”因子，当设置值为O时，表示无雨数据;当设置值为I时，表示降雨数据；
[0012]S04.从无雨海面数据，提取9°入射角无雨海面NRCS数据，所述的9°入射角是指8.5°-9.5°的入射角范围；
[0013]S05.从ECMWF模式预报数据中，提取与9°入射角无雨海面NRCS数据时间和空间上同步的动力环境参数，
[0014]S06.从降雨海面数据，分别提取9°入射角降雨海面NRCS数据和同步的平均降雨量数据；
[0015]S07.从步骤SOl和S05获得的数据中，选取典型动力环境参数下的9°入射角无雨海面NRCS，典型海面风速设定为l-13m/s，典型海面有效波高设定为0.3-3.5m;
[0016]S08.从步骤S06获得的数据中，选取有效平均降雨量下的9°入射角降雨海面NRCS，有效平均降雨量设定为大于lmm/hr;
[0017]S09.利用概率统计方法估计典型动力环境参数下9°入射角无雨海面NRCS的95%置信区间；
[0018]S1.利用概率统计方法估计有效平均降雨量下9°入射角降雨海面NRCS的95 %置信区间；
[0019]Sll.计算9°入射角无雨海面NRCS置信区间下限和9°入射角降雨海面NRCS置信区间上限的平均值，作为9°入射角海面NRCS降雨阈值，将9°入射角海面NRCS降雨阈值用于研究区域波谱仪数据确定降雨标志，当波谱仪9°入射角海面NRCS大于降雨阈值时，标志为无雨海面;当波谱仪9°入射角海面NRCS小于降雨阈值时，标志为降雨海面。
[0020]进一步的，所述步骤S05中的动力环境参数包括动力环境参数指海面风速和海面有效波尚。
[0021 ] 进一步的，所述步骤S05中同步时间窗口为10分钟，同步空间窗口为5km。
[0022]本发明的有益效果是:
[0023](I)本发明提供了一种确定波谱仪数据降雨标志的方法，利用TRMM卫星PR雷达9°入射角观测数据，基于9°入射角海面NRCS不随海面风、浪等动力环境因素变化的特性，通过比较无雨和降雨海面NRCS差异，确定了9°入射角海面NRCS降雨阈值，该阈值可用于确定波谱仪数据的降雨标志；(2)确定波谱仪数据的降雨标志，可用于剔除波谱仪在降雨条件下的受干扰的无效数据，以保证动力环境遥感产品的质量；(3)摆脱了小入射角雷达需要从大气反射率来遥感降雨量的惯性思维，为解决今后建立降雨条件下动力环境遥感模型问题奠定了研究基础。
【附图说明】
[0024]图1是确定波谱仪数据降雨标志方法的技术路线图；
[0025]图2(a)是提取的TRMM卫星PR雷达9°入射角无雨海面NRCS数据随典型海面风速变化关系；
[0026]图2(b)是提取的TRMM卫星PR雷达9°入射角无雨海面NRCS数据随典型海面有效波尚变化关系；
[0027]图3是提取的TRMM卫星PR雷达9°入射角降雨海面NRCS数据随有效平均降雨量变化关系;
[0028]图4是提取的TRMM卫星PR雷达9°入射角海面NRCS数据在不同降雨阈值下得到的漏警率(图4a)和虚警率(图4b)。
【具体实施方式】
[0029]为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及具体实例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
[0030]如图1所示，一种确定波谱仪数据降雨标志的方法，包括以下步骤:
[0031]—种确定波谱仪数据降雨标志的方法，其特征在于其包括以下步骤:
[0032]SOl.收集研究区域TRMM卫星PR雷达2A25标准产品数据；
[0033]S02.从步骤SOl获得的数据中选取海面观测数据，选取方法:提取产品中的“TRMM2A25Method Flag”因子，当设置值为O时，表示海面观测数据；
[0034]S03.从海面观测数据中分别选取无雨和降雨海面数据，选取方法:提取产品中的“Rain Flag”因子，当设置值为O时，表示无雨数据;当设置值为I时，表示降雨数据；
[0035]S04.从无雨海面数据，提取9°入射角无雨海面NRCS数据，9°入射角包括8.5-9.5° ;
[0036]S05.从ECMWF模式预报数据中，提取与9°入射角无雨海面NRCS数据时间和空间上同步的动力环境参数，动力环境参数指海面风速和海面有效波高，同步时间窗口为10分钟，同步空间窗口为5km;
[0037]S06.从降雨海面数据，分别提取9°入射角降雨海面NRCS数据和同步的平均降雨量数据；
[0038]S07.从步骤SOl和S05获得的数据中，选取典型动力环境参数下的9°入射角无雨海面NRCS，典型海面风速设定为l-13m/s，典型海面有效波高设定为0.3-3.5m。
[0039]S08.从步骤S06获得的数据中，选取有效平均降雨量下的9°入射角降雨海面NRCS，有效平均降雨量设定为大于lmm/hr。
[0040]S09.利用概率统计方法估计典型动力环境参数下9°入射角无雨海面NRCS的95%置信区间。
[0041]S10.利用概率统计方法估计有效平均降雨量下9°入射角降雨海面NRCS的95%置信区间。
[0042]Sll.计算9°入射角无雨海面NRCS置信区间下限和9°入射角降雨海面NRCS置信区间上限的平均值，作为9°入射角海面NRCS降雨阈值，当9°入射角海面NRCS大于降雨阈值时，标志为无雨海面;当9°入射角海面NRCS小于降雨阈值时，标志为降雨海面。
[0043]S12.将9°入射角海面NRCS降雨阈值用于研究区域波谱仪数据确定降雨标志。
[0044]下面以中国南海2010年01月期间TRMM卫星PR雷达200轨数据为例，具体步骤如下:
[0045]SOl.收集中国南海2010年01月期间TRMM卫星PR雷达200轨2A25标准产品数据；
[0046]S02.从步骤SOl获得的数据中选取海面观测数据，选取方法:提取产品中的“TRMM2A25Method Flag”因子，当设置值为O时，表示海面观测数据；
[0047]S03.从海面观测数据中分别选取无雨和降雨海面数据，选取方法:提取产品中的“Rain Flag”因子，当设置值为O时，表示无雨数据;当设置值为I时，表示降雨数据；
[0048]S04.从无雨海面数据，提取9°入射角无雨海面NRCS数据，9°入射角包括8.5-9.5° ;
[0049]S05.从ECMWF模式预报数据中，提取与9°入射角无雨海面NRCS数据时间和空间上同步的动力环境参数，动力环境参数指海面风速和海面有效波高，同步时间窗口为10分钟，同步空间窗口为5km;
[0050]S06.从降雨海面数据，分别提取9°入射角降雨海面NRCS数据和同步的平均降雨量数据；
[0051]S07.从步骤SOl和S05获得的数据中，选取典型动力环境参数下的9°入射角无雨海面NRCS，典型海面风速设定为l-13m/s，典型海面有效波高设定为0.3-3.5m，典型动力环境代表海面常见环境，基于典型动力环境确定出的方法具有良好的适用性。如图2(a)和(b)所示，9°入射角无雨海面NRCS基本不随典型海面风速和有效波高等动力环境参数变化，表明确定降雨标志的过程中，可以摆脱海面动力环境因素对NRCS的影响。
[0052]S08.从步骤S06获得的数据中，选取有效平均降雨量下的9°入射角降雨海面NRCS，有效平均降雨量设定为大于lmm/hr，因为小于lmm/hr平均降雨量对海面NRCS影响甚小，难以对波谱仪数据海面动力环境遥感产品产生干扰，因此这里仅选取大于lmm/hr平均降雨量下的海面NRCS。如图3所示，9°入射角降雨海面NRCS随有效平均降雨量降低，表明无雨和降雨海面NRCS存在差异。
[0053]S09.利用概率统计方法估计典型动力环境参数下9°入射角无雨海面NRCS的95%置信区间，估计的置信区间约为[5.5，7.5]。
[0054]S10.利用概率统计方法估计有效平均降雨量下9°入射角降雨海面NRCS的95%置信区间，估计的置信区间约为[1.0，7.0]。
[0055]Sll.计算9°入射角无雨海面NRCS置信区间下限和9°入射角降雨海面NRCS置信区间上限的平均值，作为9°入射角海面NRCS降雨阈值，计算的降雨阈值约为6.2dB，当9°入射角海面NRCS大于降雨阈值时，标志为无雨海面；当9°入射角海面NRCS小于降雨阈值时，标志为降雨海面。图4是提取的TRMM卫星PR雷达9°入射角海面NRCS数据在不同降雨阈值下得到的漏警率和虚警率。可以看出，当降雨阈值为6.2dB时对应的漏警率约为8.1% (图4a)，虚警率约为6.7% (图4b)，表明能有效区分无雨和降雨海面。
[0056]将9°入射角海面NRCS降雨阈值用于研究区域波谱仪数据确定降雨标志，可用于剔除波谱仪在降雨条件下的受干扰的无效数据，以保证动力环境遥感产品的质量。
【主权项】
1.一种确定波谱仪数据降雨标志的方法，其特征在于其包括以下步骤: SO1.收集研究区域TRMM卫星PR雷达2A25标准产品数据； S02.从步骤SOI获得的数据中选取海面观测数据，选取方法:提取产品中的“TRMM2A25Method Flag”因子，当设置值为O时，表示海面观测数据； 503.从海面观测数据中分别选取无雨和降雨海面数据，选取方法:提取产品中的“RainFlag”因子，当设置值为O时，表示无雨数据;当设置值为I时，表示降雨数据； 504.从无雨海面数据，提取9°入射角无雨海面NRCS数据，所述的9°入射角是指8.5° -9.5°的入射角范围； 505.从ECMWF模式预报数据中，提取与9°入射角无雨海面NRCS数据时间和空间上同步的动力环境参数， 506.从降雨海面数据，分别提取9°入射角降雨海面NRCS数据和同步的平均降雨量数据； 507.从步骤SOl和S05获得的数据中，选取典型动力环境参数下的9°入射角无雨海面NRCS，典型海面风速设定为l-13m/s，典型海面有效波高设定为0.3-3.5m; 508.从步骤S06获得的数据中，选取有效平均降雨量下的9°入射角降雨海面NRCS，有效平均降雨量设定为大于lmm/hr; 509.利用概率统计方法估计典型动力环境参数下9°入射角无雨海面NRCS的95%置信区间； 510.利用概率统计方法估计有效平均降雨量下9°入射角降雨海面NRCS的95%置信区间； 511.计算9°入射角无雨海面NRCS置信区间下限和9°入射角降雨海面NRCS置信区间上限的平均值，作为9°入射角海面NRCS降雨阈值，将9°入射角海面NRCS降雨阈值用于研究区域波谱仪数据确定降雨标志，当波谱仪9°入射角海面NRCS大于降雨阈值时，标志为无雨海面;当波谱仪9°入射角海面NRCS小于降雨阈值时，标志为降雨海面。2.如权利要求1所述的一种确定波谱仪数据降雨标志的方法，其特征在于所述步骤S05中的动力环境参数包括动力环境参数指海面风速和海面有效波高。3.如权利要求1所述的一种确定波谱仪数据降雨标志的方法，其特征在于所述步骤S05中同步时间窗口为10分钟，同步空间窗口为5km。
【文档编号】G01S13/95GK105938194SQ201610556806
【公开日】2016年9月14日
【申请日】2016年7月12日
【发明人】任林, 杨劲松, 郑罡, 王隽
【申请人】国家海洋局第二海洋研究所