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

说明书

本发明公开了一种考虑空间分布的窗口滑动GPM数据订正方法，包括如下步骤，S1、对卫星栅格数据和地面站点数据进行预处理；S2、根据预处理后的卫星栅格数据和地面站点数据，设定初始窗口；S3、对初始窗口内的卫星栅格数据进行订正，移动初始窗口依次对所有的卫星栅格数据进行订正；S4、对订正结果进行评估。优点是：能够考虑卫星降水数据的时空分布特征，将地面站点实测降水数据作为参考，对卫星栅格数据进行局部订正，降低订正误差，使得订正结果具有相对较高的精度。

技术领域

本发明涉及卫星降水数据订正技术领域，尤其涉及一种考虑空间分布的窗口滑动GPM数据订正方法。

背景技术

现有的技术是卫星降水数据融合技术，采用的方法可以划分为全局和局部两大类。其中，全局融合方法有平均偏差修正法、线性回归法、双核平滑法等；局部融合方法主要包括协同克里金、地理加权回归、贝叶斯融合等。但是对于精度本身不高、存在误差的卫星降水数据而言，直接进行融合，可能不仅无法提高产品质量，反倒会引入更大的误差，难以保证结果的可靠性。此外，常见的卫星降水数据融合技术通常没有考虑降水在空间分布上以及时间的差异性，从而导致融合结果会产生很大的不确定性。由此可见，现有技术存在一定的局限性，卫星数据的合理高效订正就显得格外重要。

发明内容

本发明的目的在于提供一种考虑空间分布的窗口滑动GPM数据订正方法，从而解决现有技术中存在的前述问题。

为了实现上述目的，本发明采用的技术方案如下：

一种考虑空间分布的窗口滑动GPM数据订正方法，包括如下步骤，

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

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

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

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

优选的，步骤S1具体包括如下内容，

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

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

优选的，步骤S2具体包括如下内容，

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

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

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

优选的，步骤S22采用卫星栅格经纬度坐标与地面站点坐标的相关关系判断矩形窗口内是否同时包含卫星栅格数据和至少一个地面站点数据，具体公式为，

|lon_卫星-lon_地面|N

|lat_卫星-lat_地面|T

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

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