[发明专利]一种考虑空间分布的窗口滑动GPM数据订正方法

权利要求书

1.一种考虑空间分布的窗口滑动GPM数据订正方法，其特征在于：包括如下步骤，

S1、对卫星栅格数据和地面站点数据进行预处理；

S2、根据预处理后的卫星栅格数据和地面站点数据，设定初始窗口；

S3、对初始窗口内的卫星栅格数据进行订正，移动初始窗口依次对所有的卫星栅格数据进行订正；

S4、对订正结果进行评估；

步骤S3具体包括如下内容，

S31、分别计算初始窗口内地面站点观测降水数据的算术平均值以及所含卫星栅格数据的算术平均值，计算初始窗口内地面站点观测降水数据的算术平均值与所含卫星栅格数据的算术平均值之间的比值；利用比值乘以初始窗口内卫星栅格数据，获取初始窗口内卫星栅格数据的订正结果；

S32、将初始窗口按照1个网格单位的步长进行移动，并按照步骤S31的方法对窗口中的卫星栅格数据依次进行订正；

S33、对于单个网格内存在的多个订正结果，进行算术平均值计算得到整个卫星栅格数据的订正结果。

2.根据权利要求1所述的考虑空间分布的窗口滑动GPM数据订正方法，其特征在于：步骤S1具体包括如下内容，

读取原始卫星.nc文件，将卫星栅格数据处理为矩阵形式，并将卫星栅格数据各网格对应的经纬度坐标分别处理为矩阵形式；

将各地面站点在一定时间尺度下的实测降水数据以及地面站点对应的经纬度坐标分别处理为矩阵形式。

3.根据权利要求2所述的考虑空间分布的窗口滑动GPM数据订正方法，其特征在于：步骤S2具体包括如下内容，

S21、从卫星栅格数据起始网格开始，固定顶点，以1个网格单位作为矩形窗口边长；

S22、判断该矩形窗口内是否同时包含卫星栅格数据和至少一个地面站点数据，若是，则进入S23，否则，将该矩形窗口的边长递增1个网格单位，并重新判断；

S23、选取该边长的矩形窗口作为初始窗口。

4.根据权利要求3所述的考虑空间分布的窗口滑动GPM数据订正方法，其特征在于：步骤S22采用卫星栅格经纬度坐标与地面站点坐标的相关关系判断矩形窗口内是否同时包含卫星栅格数据和至少一个地面站点数据，具体公式为，

|lon_卫星-lon_地面|N

|lat_卫星-lat_地面|T

其中，lon _卫星为卫星栅格数据各网格经度；lon_地面为各地面站点经度；lat_卫星为卫星栅格数据各网格纬度；lat_地面为各地面站点纬度；N、T分别为所选卫星降水数据的空间分辨率。

5.根据权利要求1所述的考虑空间分布的窗口滑动GPM数据订正方法，其特征在于：步骤S4具体为，以地面站点观测降水数据为参照，分别采用相关系数、相对偏差、均方根误差、平均绝对误差、探测率、误报率、临界成功指数、频率偏差和公正先兆评分在内的九项评价指标对订正后的卫星栅格数据进行评估；

相关系数越接近于1、相对偏差越小、均方根误差越小、平均绝对误差越小、探测率越大、误报率越小、临界成功指数越大、频率偏差越小、公正先兆评分越大，则表示卫星栅格数据的订正效果越好。