一种确定波谱仪数据降雨标志的方法技术

本发明专利技术公开了一种确定波谱仪数据降雨标志的方法，利用TRMM卫星PR雷达9°入射角观测数据，基于9°入射角海面后向散射系数(NRCS)不随海面风、浪等动力环境因素变化的特性，通过分析无雨和降雨海面NRCS差异，确定了9°入射角海面NRCS降雨阈值，该阈值可用于确定波谱仪数据的降雨标志。确定降雨标志可用于剔除波谱仪部分受降雨干扰的无效数据，对保证波谱仪关于海面动力环境遥感产品的质量，起着非常重要的作用。同时具有简单实现的优势，适合于实际应用。

【技术实现步骤摘要】

本专利技术属于海洋微波遥感
，具体是一种确定波谱仪数据降雨标志的方法。
技术介绍
降雨是一种常见的大气现象，海面的降雨会引起波谱仪海面NRCS的变化，由于目前的波谱仪动力环境遥感模型仅适用无雨海面，所以海面降雨会干扰波谱仪对海面风、浪等动力环境的遥感观测。海面降雨标志的确定，对于剔除波谱仪部分受降雨干扰的无效数据，保证遥感产品的质量，具有非常重要的现实意义和科研价值。波谱仪是一种小入射角星载雷达，入射角范围在0°-11°，主要用于测量海面的风、浪等动力环境参数，全球第一颗波谱仪卫星将搭载中法海洋卫星(CFOSAT)，预计于2018年发射(Ren et al.,2015)。目前，对于小入射角雷达，只有TRMM卫星PR雷达能够实现海面降雨量的遥感。PR雷达作为一种气象雷达，是通过观测大气的反射率，利用Z-R关系(反射率-降雨量关系)来确定降雨量，进而确定降雨标志(Iguchi et al.2000；Li et al.2002；Seto and Iguchi 2007)。而对于动力环境遥感卫星，如波谱仪，散射计和合成孔径雷达(SAR)等，都是基于观测的海面NRCS来遥感海面动力环境参数，并未配置测量大气发射率的功能模块。因此，波谱仪难以应用与PR雷达相同的方法来确定波谱仪数据降雨标志。由于波谱仪不能观测大气发射率，所以只能从海面NRCS的角度来确定降雨标志。在降雨条件下，波谱仪观测到的海面NRCS，是海面风、浪等动力环境因素和降雨的共同作用结果，而实际海面的动力环境始终处于变化之中，很难根据经验用一个预先设定的海面NRCS阈值来区分是否降雨。因此，当需要同时考虑降雨和动力环境因素对海面NRCS影响时，如何应用简单的海面NRCS降雨阈值，来确定波谱仪数据降雨标志，是一个非常困难的问题。针对上述问题，解决的关键在于如何消除动力环境因素对海面NRCS的影响。根据小入射角海面NRCS特性研究(Hesany et al.,2000；Tran,et al.,
2007)，9°-11°入射角观测的海面NRCS，基本不随海面动力环境变化，这主要是因为在该入射角范围内，海面NRCS是准镜面散射和布拉格散射的共同作用结果，准镜面散射随着粗糙度增大而减小，而布拉格散射随着粗糙度增大而增大，从而导致不同动力环境下海面NRCS整体的不变性。因此，若选取该入射角范围内观测数据来确定波谱仪数据降雨标志，则可以摆脱海面动力环境因素对NRCS的影响，保证海面NRCS变化量仅与降雨量相关，从而有利于确定海面NRCS降雨阈值。然而迄今为止，仍然没有研究将其用于确定波谱仪数据降雨标志。
技术实现思路
本专利技术的目的是提供了一种确定波谱仪数据降雨标志的方法。本专利技术通过以下技术方案实现：确定波谱仪数据降雨标志的方法包括以下步骤：S01.收集研究区域TRMM卫星PR雷达2A25标准产品数据；S02.从步骤S01获得的数据中选取海面观测数据，选取方法：提取产品中的“TRMM 2A25Method Flag”因子，当设置值为0时，表示海面观测数据；S03.从海面观测数据中分别选取无雨和降雨海面数据，选取方法：提取产品中的“Rain Flag”因子，当设置值为0时，表示无雨数据；当设置值为1时，表示降雨数据；S04.从无雨海面数据，提取9°入射角无雨海面NRCS数据，所述的9°入射角是指8.5°-9.5°的入射角范围；S05.从ECMWF模式预报数据中，提取与9°入射角无雨海面NRCS数据时间和空间上同步的动力环境参数，S06.从降雨海面数据，分别提取9°入射角降雨海面NRCS数据和同步的平均降雨量数据；S07.从步骤S01和S05获得的数据中，选取典型动力环境参数下的9°入射角无雨海面NRCS，典型海面风速设定为1-13m/s，典型海面有效波高设定为0.3-3.5m；S08.从步骤S06获得的数据中，选取有效平均降雨量下的9°入射角降雨海面NRCS，有效平均降雨量设定为大于1mm/hr；S09.利用概率统计方法估计典型动力环境参数下9°入射角无雨海面NRCS的
95％置信区间；S10.利用概率统计方法估计有效平均降雨量下9°入射角降雨海面NRCS的95％置信区间；S11.计算9°入射角无雨海面NRCS置信区间下限和9°入射角降雨海面NRCS置信区间上限的平均值，作为9°入射角海面NRCS降雨阈值，将9°入射角海面NRCS降雨阈值用于研究区域波谱仪数据确定降雨标志，当波谱仪9°入射角海面NRCS大于降雨阈值时，标志为无雨海面；当波谱仪9°入射角海面NRCS小于降雨阈值时，标志为降雨海面。进一步的，所述步骤S05中的动力环境参数包括动力环境参数指海面风速和海面有效波高。进一步的，所述步骤S05中同步时间窗口为10分钟，同步空间窗口为5km。本专利技术的有益效果是：(1)本专利技术提供了一种确定波谱仪数据降雨标志的方法，利用TRMM卫星PR雷达9°入射角观测数据，基于9°入射角海面NRCS不随海面风、浪等动力环境因素变化的特性，通过比较无雨和降雨海面NRCS差异，确定了9°入射角海面NRCS降雨阈值，该阈值可用于确定波谱仪数据的降雨标志；(2)确定波谱仪数据的降雨标志，可用于剔除波谱仪在降雨条件下的受干扰的无效数据，以保证动力环境遥感产品的质量；(3)摆脱了小入射角雷达需要从大气反射率来遥感降雨量的惯性思维，为解决今后建立降雨条件下动力环境遥感模型问题奠定了研究基础。附图说明图1是确定波谱仪数据降雨标志方法的技术路线图；图2(a)是提取的TRMM卫星PR雷达9°入射角无雨海面NRCS数据随典型海面风速变化关系；图2(b)是提取的TRMM卫星PR雷达9°入射角无雨海面NRCS数据随典型海面有效波高变化关系；图3是提取的TRMM卫星PR雷达9°入射角降雨海面NRCS数据随有效平均降雨量变化关系；图4是提取的TRMM卫星PR雷达9°入射角海面NRCS数据在不同降雨阈值下得到的漏警率(图4a)和虚警率(图4b)。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及具体实例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实例仅仅用以解释本专利技术，并不用于限定本专利技术。此外，下面所描述的本专利技术各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。如图1所示，一种确定波谱仪数据降雨标志的方法，包括以下步骤：一种确定波谱仪数据降雨标志的方法，其特征在于其包括以下步骤：S01.收集研究区域TRMM卫星PR雷达2A25标准产品数据；S02.从步骤S01获得的数据中选取海面观测数据，选取方法：提取产品中的“TRMM 2A25Method Flag”因子，当设置值为0时，表示海面观测数据；S03.从海面观测数据中分别选取无雨和降雨海面数据，选取方法：提取产品中的“Rain Flag”因子，当设置值为0时，表示无雨数据；当设置值为1时，表示降雨数据；S04.从无雨海面数据，提取9°入射角无雨海面NRCS数据，9°入射角包括8.5-9.5°；S05.从ECMWF模式预报数据中，提取与9°入射角无雨海面NRCS数据时间和空间上同步的动力环境参数，动力环境参数指海面风速和海面有效波高，同步时间窗口为10分钟，同步空间窗口为5本文档来自技高网...

【技术保护点】
一种确定波谱仪数据降雨标志的方法，其特征在于其包括以下步骤：S01.收集研究区域TRMM卫星PR雷达2A25标准产品数据；S02.从步骤S01获得的数据中选取海面观测数据，选取方法：提取产品中的“TRMM 2A25Method Flag”因子，当设置值为0时，表示海面观测数据；S03.从海面观测数据中分别选取无雨和降雨海面数据，选取方法：提取产品中的“Rain Flag”因子，当设置值为0时，表示无雨数据；当设置值为1时，表示降雨数据；S04.从无雨海面数据，提取9°入射角无雨海面NRCS数据，所述的9°入射角是指8.5°‑9.5°的入射角范围；S05.从ECMWF模式预报数据中，提取与9°入射角无雨海面NRCS数据时间和空间上同步的动力环境参数，S06.从降雨海面数据，分别提取9°入射角降雨海面NRCS数据和同步的平均降雨量数据；S07.从步骤S01和S05获得的数据中，选取典型动力环境参数下的9°入射角无雨海面NRCS，典型海面风速设定为1‑13m/s，典型海面有效波高设定为0.3‑3.5m；S08.从步骤S06获得的数据中，选取有效平均降雨量下的9°入射角降雨海面NRCS，有效平均降雨量设定为大于1mm/hr；S09.利用概率统计方法估计典型动力环境参数下9°入射角无雨海面NRCS的95％置信区间；S10.利用概率统计方法估计有效平均降雨量下9°入射角降雨海面NRCS的95％置信区间；S11.计算9°入射角无雨海面NRCS置信区间下限和9°入射角降雨海面NRCS置信区间上限的平均值，作为9°入射角海面NRCS降雨阈值，将9°入射角海面NRCS降雨阈值用于研究区域波谱仪数据确定降雨标志，当波谱仪9°入射角海面NRCS大于降雨阈值时，标志为无雨海面；当波谱仪9°入射角海面NRCS小于降雨阈值时，标志为降雨海面。...

【技术特征摘要】
1.一种确定波谱仪数据降雨标志的方法，其特征在于其包括以下步骤：S01.收集研究区域TRMM卫星PR雷达2A25标准产品数据；S02.从步骤S01获得的数据中选取海面观测数据，选取方法：提取产品中的“TRMM 2A25Method Flag”因子，当设置值为0时，表示海面观测数据；S03.从海面观测数据中分别选取无雨和降雨海面数据，选取方法：提取产品中的“Rain Flag”因子，当设置值为0时，表示无雨数据；当设置值为1时，表示降雨数据；S04.从无雨海面数据，提取9°入射角无雨海面NRCS数据，所述的9°入射角是指8.5°-9.5°的入射角范围；S05.从ECMWF模式预报数据中，提取与9°入射角无雨海面NRCS数据时间和空间上同步的动力环境参数，S06.从降雨海面数据，分别提取9°入射角降雨海面NRCS数据和同步的平均降雨量数据；S07.从步骤S01和S05获得的数据中，选取典型动力环境参数下的9°入射角无雨海面NRCS，典型海面风速设定为1-13m/s，典型海面有效波高设定为0.3-...

【专利技术属性】
技术研发人员：任林，杨劲松，郑罡，王隽，
申请(专利权)人：国家海洋局第二海洋研究所，
类型：发明
国别省市：浙江;33