基于卫星遥感与回归克里格的地面日降雨量预测方法
【专利摘要】本发明公开了一种基于卫星遥感与回归克里格的地面日降雨量预测方法。首先利用卫星遥感快速获取数据，根据降雨量分级，建立各级降雨量的地面观测值、TRMM、DEM及地理位置间的回归关系，得到各级回归估计值和回归残差。其次分析各级回归残差的空间集聚度，对回归残差进行趋势去除，并进行回归残差的克里格插值，得到1km的各级回归残差空间分布特征。然后分别将各级回归估计值与回归残差相加，得到1km的各级降雨的地面预测值。最终将各级降雨的地面预测值合并，得到1km的日降雨量预测值。本发明能准确地预测地面降雨量的时空分布特征，提高了地面日降雨的预测精度，又提高了预测的空间分辨率，解决水利部门进行地面降雨预报的关键问题。
【专利说明】基于卫星遥感与回归克里格的地面日降雨量预测方法【技术领域】
[0001]本发明涉及一种地面降雨预测方法，尤其是涉及一种基于卫星遥感与回归克里格的地面日降雨量预测方法。
技术背景
[0002]降水时空变化较为剧烈，准确预测降水对水文、气象及灾害预报等具有重要意义。目前常用的降水数据主要通过地面观测站和卫星测雨等方法获取。地面观测站作为一种常规的测量降水手段，具有应用广泛、精度较高且技术成熟的特点，但是地面观测站布网密度满足相关研究工作的质量需要，且地面观测站观测的降水量仅代表地面观测站及周边一定范围的降水状况，在表述大面积降水分布特征时可能引起偏差。卫星测雨遥感技术能够提供较高时空分辨率的降水数据，覆盖空间范围更广，克服了地面雨量站和测雨雷达的局限，为监测全球降水变化提供有力的数据支持。
[0003]近年来，随着气象卫星技术的发展，全球尺度高时空分辨率的测雨卫星应运而生，其中利用热带降雨测量卫星(Tropical Rainfall Measuring Mission, TRMM)上搭载的微波成像仪(TRMM Microwave Imager, TMI)及测雨雷达(Precipitation Radar, PR)验证降雨已经取得了长足的进步。TRMM降雨卫星是第一颗专门用于定量观测热带、亚热带降雨的气象卫星，能够以较高的时空分辨率，提供覆盖全球50° S~50° N以内的区域的降雨数据。
[0004]国内外学者在利用TRMM卫星数据进行降水分析和降水资料验证等方面的工作上做出了大量贡献，Islam 等(Islam M N, Uyeda H.Use of TRMM in determiningthe climatic characteristics of rainfall over Bangladesh.Remote Sensing ofEnvironment, 2007,108:264— 276.)利用TRMM卫星数据在孟加拉国的对比结果精度达到96%,表明使用TRMM卫星反演降水能够很好的反映区域的降水状况。Huffman等(HuffmanG J, Robert F A, Dav`id T B.The TRMM Multisatellite Precipitation Analysis (TMPA):Quas1-global, multiyear,combined-sensor precipitation estimates at fine scales.Journal of Hydrometeorology, 2006, 8(1):38-55.)在澳大利亚进行了 TRMM 卫星数据与气象观测数据、雷达观测数据的一致性检测，结果表明TRMM卫星数据与观测数据在日、月尺度上均具有较高的线性相关特征。杨传国等(杨传国，余钟波，林朝晖.基于TRMM卫星雷达降雨的流域陆面水文过程.水科学进展，2009，20 (4):461-466.)利用TRMM卫星数据研究流域尺度陆面水文过程，结果表明利用TRMM卫星数据模拟的结果与利用地面观测数据模拟的结果相似。白爱娟等(白爱娟，方建刚，张科翔.TRMM卫星资料对陕西及周边地区夏季降水的探测.灾害学，2008，23(2):41-45.)采用TRMM卫星数据在陕西及其周边地区与地面气象站观测资料进行了对比，发现TRMM卫星数据对区域夏季降水过程的具有良好的探测能力。
[0005]以上研究表明TRMM卫星数据在预测降雨方面具有很强的适应性和可行性。但是，TRMM卫星的原始分辨率较低(空间分辨率为0.25。，约30km)，在预测区域尺度降雨方面具有局限性和误差，因此需要针对TRMM数据进行空间尺度的转换，以得到较高空间分辨率的降雨量预测值。Jia 等(Jia, S.，Zhu, ff., L Q , A., &Yan, T.(2011).Astatistical spatialdownscaling algorithm of TRMM precipitation based on NDVI and DEM in the QaidamBasin of China.Remote Sensing of Environment, 115，3069-3079.)利用 NDVI 和 DEM数据，对柴达木盆地的TRMM降雨数据进行降尺度研究，发现降尺度后的r2从0.72上升到
0.96，取得较好成果。马金辉等(马金辉，屈创，张海莜，夏燕秋.2001-2010年石羊河流域上游TRMM降水资料的降尺度研究.地理科学进展，2013, 32 (9): 1423-1432.)利用Ikm DEM数据，采用回归方程结合残差的Thiessen插值方法，将TRMM年降水量的空间分辨率提高到lkm。经检验发现，降尺度后的数据精度有一定程度的提高。
[0006]作为一种间接估测降雨的工具，TRMM卫星不仅存在空间尺度上的局限性，而且在应用到不同的地形及降雨条件时，同样具有很大程度上的不确定性及误差。Tian等(Tian, Y.D., and Peters-Lidard, C.D., 2007, Systematic anomalies over inland waterbodies in satellite-based precipitation estimates.Geophysical Research Letters, 34, L14403, do1:10.1029/2007GL030787.)发现TRMM卫星数据应用到海洋及南美洲等地时，不确定性为40 - 60%,而应用到高纬度地区、海岸带、水体等地区时，不确定性达到了100 - 140%oAs_Syakur 等(As-Syakur, A.R., Tanaka, T., Prasetia, R., Swardika, 1.K., andKasa, 1.ff., 2011, Comparison of TRMM multisatellite precipitation analysis (TMPA)products and daily-monthly gauge data over Bal1.1nternational Journal of RemoteSensing, 32，pp.8969-8982.)通过分析巴厘岛地面观测站点数据与TRMM卫星数据之间的关系，发现在干季时，TRMM卫星数据与地面感测数据具有较好的相关性，而这种相关性在湿季较差。
[0007]因此，需要考虑地形、经纬度和距海岸线位置等因素，针对不同的降雨量级别，分别建立不同的预测模型，从而得到更准确地降雨量预测值。这种新的思路需要引进DEM、经纬度信息等，针对不同级别降雨量，对TRMM数据的误差进行分别分析，这就必须引入地统计及空间分析相关方法。
[0008]回归克里格(regression kring)方法是一种结合了回归模型和残差空间插值的混合空间建模方法，以常规的多元回归分析技术和相关的辅助预测数据(如DEM)来分离趋势项，然后对残差进行克里格插值，最后将多元回归预测的趋势项和残差的克里格估计值相加进行估测，回归克里格方法中的趋势项是辅助预测数据的函数。该方法既分析主要影响因素也考虑随机因素，既模拟其空间分布趋势也模拟不确定性，不仅能提高模型的空间精度，预测精度也要优于其他模型。Lloyd等(Lloyd, C.D.2005.Assessing the effectof integrating elevation data into the estimation of monthly precipitation inGreat Britain.Journal of Hydrology.308.pp.128-150.)将回归克里格方法用于地面站点观测降雨量的制图；Yemefack 等(Yemefack, M., D.G.Rossiter, and R.Njomgang.2005.Mult1-scale characterization of soil variability within an agriculturallandscape mosaic system in southern Cameroon.Geoderma.125.pp.117-143.)利用回归克里格方法进行土壤属性预测。
【发明内容】
[0009]本发明的目的在于提供一种基于卫星遥感与回归克里格的地面日降雨量预测方法，结合降雨分级，分析不同模型中残差的空间趋势，建立不同的预测模型来估算地面降雨量的空间分布特征。
[0010]本发明采用的技术方案的步骤如下:
[0011]步骤(I)数据获取:获取TRMM气象卫星遥感影像数据及ASTERGDEM卫星遥感影像数据，同时收集地面观测站点的日降雨量观测值；
[0012]步骤(2)数据预处理:将步骤(I)获取的TRMM气象卫星遥感影像数据的时间分辨率处理为与地面观测站点的日降雨量观测值的时间分辨率相等，同时将ASTER⑶EM卫星遥感影像数据的空间分辨率重采样为1km，并将数据以文件形式存放；
[0013]步骤(3)卫星遥感影像数据提取:根据地面观测站点的空间位置信息，分别提取步骤(2 )中相应的TRMM气象卫星遥感影像数据及ASTER⑶EM卫星遥感影像数据的栅格值，从而得到每个地面观测站点对应的卫星遥感影像数据的栅格值；
[0014]步骤(4)日降雨量数据的分位数分级:根据步骤(I)收集的地面观测站点的日降雨量观测值及步骤(3)中获取的卫星遥感影像数据的栅格值，进行日降雨量数据的分位数分级；
[0015]步骤(5)各级线性回归关系建立:根据步骤(4)的分位数分级，建立地面观测站点的日降雨量观测值、对应TRMM气象卫星遥感影像数据的栅格值、ASTER⑶EM卫星遥感影像数据的栅格值以及经纬度坐标之间的线性回归关系，得到各级降雨量的回归估计值与回归
残差;
[0016]步骤(6)各级回归残差分析:分析步骤(5)中得到的各级降雨量的回归残差的空间自相关性及空间集聚性；
[0017]步骤(7)各级回归残差处理:对步骤(6)中具有空间集聚性的各级回归残差进行趋势去除，然后根据ASTER⑶EM卫星遥感影像数据对各级回归残差进行Block克里格插值，得到Ikm空间分辨率的各级回归残差的空间分布特征图；
[0018]步骤(8)各级降雨量预测:将根据步骤(5)得到的各级降雨量的回归估计值与步骤(7)得到的Ikm空间分辨率的各级回归残差相加，得到Ikm空间分辨率的各级降雨量预测值；
[0019]步骤(9)地面日降雨量预测:将根据步骤(8)所得的各级降雨量预测值合并，得到地面Ikm空间分辨率的日降雨量预测值；
[0020]步骤(10)地面日降雨量预测精度评价:利用交叉检验的方法对步骤(9)中的地面日降雨量预测值进行预测精度验证，交叉检验选用均方根误差、平均误差以及偏差作为评价参数。
[0021]步骤(I)中，所述的TRMM气象卫星遥感影像数据空间分辨率为0.25° X0.25°，时间分辨率为3小时，所述的ASTER⑶EM卫星遥感影像数据的空间分辨率为30m。
[0022]步骤(2 )中，所述的时间分辨率为I天。
[0023]步骤(3)中，所述的两种卫星遥感影像数据的栅格值的提取采用的是点对点的方法，即分别提取最靠近地面站点的TRMM气象卫星遥感影像数据及ASTER⑶EM卫星遥感影像数据的栅格值。
[0024]步骤(4)中，所述的日降雨量数据的分位数分级方法为四分位法。[0025]步骤(5)中，所述的线性回归关系建立采用的是多元线性回归模型，其具体公式如下，
[0026]Z (X0) = α + β J (X0) + β 2D (X0) + β 3Lon (X0) + β 4Lat (X0) + ε (χ0)
[0027]其中，Z(Xtl)为地面观测站点Xtl的降雨回归估计值，α、β 1、β 2、β 3、β 4均是回归系数，T(X0)是地面观测站点Xtl对应的TRMM卫星估测值，D(X0)是地面观测站点Xtl对应的ASTER⑶EM卫星遥感影像数据的栅格值，Lon(Xtl)为xQ的经度坐标值，Lat (xQ)为xQ的纬度坐标值，ε (Xtl)为回归残差。
[0028]步骤(6)中，所述的回归残差的空间自相关性采用半方差模型来判断，
【权利要求】
1.一种基于卫星遥感与回归克里格的地面日降雨量预测方法，其特征在于，包括以下步骤: 步骤(1)数据获取:获取TRMM气象卫星遥感影像数据及ASTERGDEM卫星遥感影像数据，同时收集地面观测站点的日降雨量观测值； 步骤(2)数据预处理:将步骤(1)获取的TRMM气象卫星遥感影像数据的时间分辨率处理为与地面观测站点的日降雨量观测值的时间分辨率相等，同时将ASTER⑶EM卫星遥感影像数据的空间分辨率重采样为1km，并将数据以文件形式存放； 步骤(3)卫星遥感影像数据提取:根据地面观测站点的空间位置信息，分别提取步骤(2)中相应的TRMM气象卫星遥感影像数据及ASTER⑶EM卫星遥感影像数据的栅格值，从而得到每个地面观测站点对应的卫星遥感影像数据的栅格值； 步骤(4)日降雨量数据的分位数分级:根据步骤(1)收集的地面观测站点的日降雨量观测值及步骤(3)中获取的卫星遥感影像数据的栅格值，进行日降雨量数据的分位数分级； 步骤(5)各级线性回归关系建立:根据步骤(4)的分位数分级，建立地面观测站点的日降雨量观测值、对应TRMM气象卫星遥感影像数据的栅格值、ASTER⑶EM卫星遥感影像数据的栅格值以及经纬度坐标之间的线性回归关系，得到各级降雨量的回归估计值与回归残差; 步骤(6)各级回归残差分析:分析步骤(5)中得到的各级降雨量的回归残差的空间自相关性及空间集聚性； 步骤(7)各级回归残 差处理:对步骤(6)中具有空间集聚性的各级回归残差进行趋势去除，然后根据ASTER⑶EM卫星遥感影像数据对各级回归残差进行Block克里格插值，得到Ikm空间分辨率的各级回归残差的空间分布特征图； 步骤(8)各级降雨量预测:将根据步骤(5)得到的各级降雨量的回归估计值与步骤(7)得到的Ikm空间分辨率的各级回归残差相加，得到Ikm空间分辨率的各级降雨量预测值；步骤(9)地面日降雨量预测:将根据步骤(8)所得的各级降雨量预测值合并，得到地面Ikm空间分辨率的日降雨量预测值； 步骤(10)地面日降雨量预测精度评价:利用交叉检验的方法对步骤(9)中的地面日降雨量预测值进行预测精度验证，交叉检验选用均方根误差、平均误差以及偏差作为评价参数。
2.根据权利要求1所述的基于卫星遥感与回归克里格的地面日降雨量预测方法，其特征在于，步骤(1)中，所述的TRMM气象卫星遥感影像数据空间分辨率为0.25° X0.25°，时间分辨率为3小时，所述的ASTER⑶EM卫星遥感影像数据的空间分辨率为30m。
3.根据权利要求1所述的基于卫星遥感与回归克里格的地面日降雨量预测方法，其特征在于，步骤(2)中，所述的时间分辨率为I天。
4.根据权利要求1所述的基于卫星遥感与回归克里格的地面日降雨量预测方法，其特征在于，步骤(3)中，所述的两种卫星遥感影像数据的栅格值的提取采用的是点对点的方法，即分别提取最靠近地面站 点的TRMM气象卫星遥感影像数据及ASTER⑶EM卫星遥感影像数据的栅格值。
5.根据权利要求1所述的基于卫星遥感与回归克里格的地面日降雨量预测方法，其特征在于，步骤(4)中，所述的日降雨量数据的分位数分级方法为四分位法。
6.根据权利要求1所述的基于卫星遥感与回归克里格的地面日降雨量预测方法，其特征在于，步骤(5)中，所述的线性回归关系建立采用的是多元线性回归模型，其具体公式如下，
7.根据权利要求1所述的基于卫星遥感与回归克里格的地面日降雨量预测方法，其特征在于，步骤(6)中，所述的回归残差的空间自相关性采用半方差模型来判断，
8.根据权利要求1所述的基于卫星遥感与回归克里格的地面日降雨量预测方法，其特征在于，步骤(6)中，所述的回归残差的空间聚集性利用Moran’ I系数来判断，
9.根据权利要求1所述的基于卫星遥感与回归克里格的地面日降雨量预测方法，其特征在于，步骤(7)中，所述的各级回归残差进行趋势去除，是采用二阶趋势去除法去除。
10.根据权利要求1所述的基于卫星遥感与回归克里格的地面日降雨量预测方法，其特征在于，步骤(8)中，所述的各级降雨量预测采用以下公式计算，
Z(xQ) = a +^1T(X0) + β2Ο(χ0) + β? Lon(X0) + fi4Lat(x0) + ?(χ0) 其中，是各级降雨量的预测值，?、/5,%我、爲、歲均是回归估计系数，Τ(χ0)是地面观测站点Xtl对应的TRMM卫星估测值，D(x0)是地面观测站点X0对应的ASTER⑶EM卫星遥感影像数据的栅格值，Lon(X0)为Xtl的经度坐标值，Lat(Xci)为Χ(ι的纬度坐标值，啟.?)为回归残差的估计值。
【文档编号】G06F19/00GK103810376SQ201410021364
【公开日】2014年5月21日 申请日期:2014年1月17日 优先权日:2014年1月17日
【发明者】史舟, 滕洪芬, 马自强, 张健 申请人:浙江大学