李 慧, 楊 濤, 何祺勝, 任偉偉

(1.荒漠與綠洲生態(tài)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 中國(guó)科學(xué)院 新疆生態(tài)與地理研究所, 烏魯木齊830011;2.河海大學(xué) 水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 南京210098;3.河海大學(xué) 地球科學(xué)與工程學(xué)院, 南京210098; 4.中國(guó)科學(xué)院大學(xué), 北京 100101)

新疆天山山區(qū)TRMM衛(wèi)星降水?dāng)?shù)據(jù)的精度檢驗(yàn)和校正方法研究

李 慧1,4, 楊 濤1,2, 何祺勝3, 任偉偉1,4

(1.荒漠與綠洲生態(tài)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 中國(guó)科學(xué)院 新疆生態(tài)與地理研究所, 烏魯木齊830011;2.河海大學(xué) 水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 南京210098;3.河海大學(xué) 地球科學(xué)與工程學(xué)院, 南京210098; 4.中國(guó)科學(xué)院大學(xué), 北京 100101)

在降水資料缺乏的天山山區(qū)，基于重建時(shí)間序列后的NDVI和DEM數(shù)據(jù)，采用CART算法對(duì)TRMM3B43月降水?dāng)?shù)據(jù)進(jìn)行校正。利用研究區(qū)25個(gè)站點(diǎn)實(shí)測(cè)降水量對(duì)校正前的TRMM降水?dāng)?shù)據(jù)和校正后的TRMM降水?dāng)?shù)據(jù)分別進(jìn)行精度檢驗(yàn)。結(jié)果表明：校正前TRMM月降水與站點(diǎn)實(shí)測(cè)降水有很好的一致性，存在顯著的線性相關(guān)關(guān)系，但誤差較大；TRMM降水與實(shí)測(cè)降水的決定系數(shù)(R2)隨時(shí)間尺度的增大而減小，相對(duì)誤差(δ)和均方根誤差(RMSE)則隨之增大，說明TRMM遙感數(shù)據(jù)的精度隨時(shí)間尺度的增加而減??；校正后TRMM降水精度得到了顯著地提高，與實(shí)測(cè)降水在月、季、年尺度上的R2分別為：0.97，0.87，0.83，相比校正前R2提高了10%以上，誤差也有明顯的減小。這說明在天山山區(qū)，使用CART對(duì)TRMM降水?dāng)?shù)據(jù)進(jìn)行校正的方法可行。

天山山區(qū); TRMM; 校正; CART; 精度檢驗(yàn)

降水作為地球水循環(huán)的主要驅(qū)動(dòng)因子，不僅影響著一個(gè)地區(qū)或流域的水分和熱量狀況，也是全球水熱循環(huán)的推動(dòng)力之一[1]。準(zhǔn)確的降水空間分布數(shù)據(jù)，理論上可由高密度分布的站點(diǎn)通過雨量計(jì)測(cè)量的方法采集獲得[2]，但在氣象站分布稀少的西部山區(qū)，可利用的觀測(cè)數(shù)據(jù)有限，所以有關(guān)山區(qū)局地降水分布的研究相對(duì)較少，若能充分利用遙感數(shù)據(jù)產(chǎn)品，則在大面積的空間范圍上為降水的時(shí)空特征研究提供了有利支持[3]。1997年美國(guó)NASA(National Aeronautical and Space Administration)和日本NASDA(National Space Development Agency)聯(lián)合發(fā)射了熱帶降水測(cè)量TRMM(Tropical Rainfall Measuring Mission)衛(wèi)星，它的主要任務(wù)是觀測(cè)熱帶和亞熱帶地區(qū)降水特征及能量交換。TRMM衛(wèi)星可向陸地發(fā)回多種時(shí)間和空間分辨率的探測(cè)數(shù)據(jù)，這為分析降水、閃電、云量等天氣現(xiàn)象提供了最詳細(xì)最直接的信息[4-5]，也使確定降水的空間分布成為可能。TRMM衛(wèi)星降水資料在熱帶濕潤(rùn)地區(qū)有大量應(yīng)用，但要將其用于我國(guó)西北干旱半干旱地區(qū)，還需要進(jìn)一步校正處理[6]，提高該降水?dāng)?shù)據(jù)的精度。

針對(duì)TRMM 數(shù)據(jù)的校正，目前主要方法是傳統(tǒng)的統(tǒng)計(jì)分析方法，如回歸分析等。在回歸分析中，當(dāng)自變量間存在共線性時(shí)，會(huì)影響回歸系數(shù)的正確估計(jì)，致使1個(gè)變量的作用會(huì)掩蓋另1個(gè)的變量作用。決策樹模型可有效避免此類問題，它作為數(shù)據(jù)挖掘中的一種分類技術(shù)，它在根據(jù)各個(gè)屬性變量將數(shù)據(jù)集進(jìn)行分類時(shí)，可以提取一系列的分類規(guī)則，這些規(guī)則中所包含的變量就可能是影響因素。CART是決策樹中比較經(jīng)典的算法，相比其他的決策樹，CART算法的優(yōu)點(diǎn)是它使用二元分支，能夠充分用到全部數(shù)據(jù)；還可以處理孤立點(diǎn)，將孤立點(diǎn)變成單獨(dú)的節(jié)點(diǎn)，這樣不會(huì)影響下一步分裂；并且CART算法能夠處理空缺值，可在一定程度上減小信息的偏差；CART算法應(yīng)用已經(jīng)非常廣泛[7]。

植被是聯(lián)結(jié)土壤、大氣和水分的自然“紐帶”，在全球變化研究中起到“指示器”作用[8]。因此有文獻(xiàn)指出，在分析研究區(qū)域降水時(shí)，將植被作為一個(gè)因子進(jìn)行考慮，并建立了TRMM降水?dāng)?shù)據(jù)與植被指數(shù)NDVI的模型[9-11]。但是已有研究表明，利用MODIS，AVHRR和SPOT4/VEGETATION等中低分辨率的傳感器進(jìn)行觀測(cè)時(shí)，太陽光照角度、觀測(cè)視角以及云的條件隨時(shí)間變化都很大，因此得到的觀測(cè)值是來自地表的雙向反射率信息，包含很多不可預(yù)測(cè)的噪音，嚴(yán)重影響了土地覆被和陸地生態(tài)系統(tǒng)的變化監(jiān)測(cè)結(jié)果[12]。因此，為了真實(shí)地反映出全球氣候變化下植被的動(dòng)態(tài)響應(yīng)，有必要對(duì)NDVI 時(shí)間序列進(jìn)行重建，即通過一系列預(yù)處理方法降低遙感NDVI 合成數(shù)據(jù)中的噪聲水平[13]。簡(jiǎn)單的方法是取一旬之中NDVI 最大值作為這個(gè)期間的NDVI值，由此來減少云的干擾，這種方法叫做最大值合成法(MVC)，然而，如果一旬時(shí)間內(nèi)所有圖像都有云的話，就無法得到該旬無云圖像。因此，僅僅使用最大值合成法是不夠的，還必須在最大值合成后做進(jìn)一步的處理，本文運(yùn)用基于傅立葉變換的HANTS算法來實(shí)現(xiàn)NDVI數(shù)據(jù)時(shí)間序列的重建[14]。

1 研究區(qū)概況及數(shù)據(jù)

1.1 研究區(qū)概況

天山山區(qū)位于39°—45.5°N，73°—96°E，地處歐亞大陸腹地，既有中溫帶大陸干旱氣候特點(diǎn)，又有垂直氣候特征，屬于典型的大陸型氣候。天山山區(qū)降水分布不均，時(shí)間和空間差異懸殊，夏季降水量大，而冬季降水量較小，山區(qū)降水量豐富，自西向東呈逐漸減少的趨勢(shì)，北坡的降水明顯多于南坡。天山北坡多年平均降水量可達(dá)500～700 mm，而且在西部個(gè)別迎風(fēng)坡年降水量多達(dá)1 000 mm，這種條件形成了西北干旱區(qū)的“濕島”景觀[15]。有研究顯示，天山山區(qū)的中山帶海拔3 000 m左右是降水量最大的地區(qū)，并且天山山區(qū)孕育了全疆65%的河流，其徑流量占總流量的54%[16]。天山山區(qū)植被結(jié)構(gòu)復(fù)雜，生態(tài)環(huán)境多樣性強(qiáng)，主要以草地覆蓋為主，森林帶上限多處于低處，林相稀疏。因此在分析本研究區(qū)域降水時(shí)，將植被作為一個(gè)因子進(jìn)行考慮。

1.2 數(shù)據(jù)來源

本研究選用的數(shù)據(jù)產(chǎn)品包括TRMM3B43衛(wèi)星雷達(dá)降水?dāng)?shù)據(jù)、植被指數(shù)NDVI數(shù)據(jù)、DEM(數(shù)字高程模型)數(shù)據(jù)和25個(gè)氣象臺(tái)站降水?dāng)?shù)據(jù)。

TRMM降水?dāng)?shù)據(jù)由美國(guó)航空航天局NASA (http：∥trmm.gsfc.nasa.gov/)提供，其中3B43數(shù)據(jù)為綜合降水?dāng)?shù)據(jù)產(chǎn)品，空間分辨率0.25°×0.25°，時(shí)間間隔為1個(gè)月；植被指數(shù)NDVI數(shù)據(jù)由美國(guó)NASA全球監(jiān)測(cè)與模型研究組發(fā)布，空間分辨率為1 km×1 km；DEM數(shù)據(jù)由NASA 提供，空間分辨率為90 m×90 m；天山山區(qū)25個(gè)氣象臺(tái)站的月降水資料由中國(guó)氣象局國(guó)家氣候中心提供，具體位置及海拔信息見表1。以上所有數(shù)據(jù)的時(shí)間均選取2000年1月—2014年12月。

1.3 數(shù)據(jù)預(yù)處理

TRMM3B43數(shù)據(jù)為NC格式，利用MATLAB軟件批處理為文本格式；NDVI數(shù)據(jù)為HDF格式，在ENVI中對(duì)其進(jìn)行拼接、投影、裁剪等處理從而得到研究區(qū)NDVI數(shù)據(jù)，為進(jìn)一步消除云、大氣、太陽高度角等影響，保證NDVI反映的是年地表植被覆蓋狀況[10]，本文用HANTS變換對(duì)NDVI進(jìn)行了處理，得到研究區(qū)重建時(shí)間序列的NDVI數(shù)據(jù)；利用DEM數(shù)據(jù)，基于ArcGIS將其分辨率重采樣為1 km后提取研究區(qū)的坡度和坡向。

表1 天山山區(qū)氣象站點(diǎn)信息

2 研究方法

2.1 校正方法

本文將采用決策樹中CART算法對(duì)TRMM降水?dāng)?shù)據(jù)進(jìn)行校正?；舅枷霝閇17]：從頂結(jié)點(diǎn)開始，CART在每一結(jié)點(diǎn)選取最優(yōu)變量的最優(yōu)分化，形成一顆大樹。然后修剪這一顆大樹獲得一系列子樹，并從中挑選出真誤差率的估計(jì)值為最小的子樹。最后根據(jù)這顆樹對(duì)客體進(jìn)行預(yù)測(cè)。通常，即使學(xué)習(xí)集很大，即學(xué)習(xí)集包含許多客體，并且每一客體有許多變量的測(cè)量值，CART樹的大小仍會(huì)比較適中。因?yàn)樵诤芏嗲闆r下，一小部分變量即包含了絕大部份的預(yù)測(cè)信息，CART樹利用這些變量即可設(shè)法捕獲絕大部分的預(yù)測(cè)信息。

本文基于MATLAB平臺(tái)，使用CART算法校正時(shí)，輸入的變量為：TRMM數(shù)據(jù)、坡度、坡向、經(jīng)度、緯度、海拔和重建時(shí)間序列后的NDVI。期望輸出為實(shí)測(cè)降水量，實(shí)際輸出即為校正后的TRMM降水量。

2.2 精度檢驗(yàn)

通過25個(gè)氣象站點(diǎn)的實(shí)測(cè)降水?dāng)?shù)據(jù)對(duì)TRMM和校正后的TRMM數(shù)據(jù)進(jìn)行檢驗(yàn)，采用以下3個(gè)指標(biāo)驗(yàn)證：

決定系數(shù)(R2)

相對(duì)誤差(δ)

均方根誤差(RMSE)

3 結(jié)果與分析

以25個(gè)站點(diǎn)的實(shí)測(cè)降水?dāng)?shù)據(jù)為準(zhǔn)，分別從月、季、年時(shí)間尺度分析TRMM降水?dāng)?shù)據(jù)精度，并通過CART算法對(duì)TRMM降水?dāng)?shù)據(jù)進(jìn)行校正，比較校正前后的精度差異。

3.1 校正前精度檢驗(yàn)

(1) 月尺度。將天山山區(qū)2000—2014年25個(gè)站點(diǎn)的實(shí)測(cè)降水?dāng)?shù)據(jù)作為降水量的“真實(shí)值”，以各雨量站15 a月降水量為自變量，以與其對(duì)應(yīng)的TRMM3B43月降水量為因變量做散點(diǎn)圖(圖1)進(jìn)行一元線性回歸分析，并將兩組數(shù)據(jù)做成多年平均月降水的時(shí)間序列圖(圖2)。

圖1給出了TRMM 降水?dāng)?shù)據(jù)和地面雨量站2000—2014年各月降水量的散點(diǎn)分布一元線性回歸方程。經(jīng)過檢驗(yàn)，R2=0.87，δ=22.40%，RMSE=3.89 mm，可以看出TRMM月降水?dāng)?shù)據(jù)與實(shí)測(cè)數(shù)據(jù)之間存在顯著的線性相關(guān)性，但誤差較大。

圖1TRMM降水?dāng)?shù)據(jù)與實(shí)測(cè)月降水?dāng)?shù)據(jù)散點(diǎn)圖

圖2中，2000—2014年TRMM月降水?dāng)?shù)據(jù)與站點(diǎn)實(shí)測(cè)的多年平均月降水?dāng)?shù)據(jù)的總體變化趨勢(shì)是一致的。經(jīng)計(jì)算，TRMM降水?dāng)?shù)據(jù)普遍高于實(shí)測(cè)降水?dāng)?shù)據(jù)，尤其是降水量較少的1—3月、11月、12月，TRMM降水?dāng)?shù)據(jù)比實(shí)測(cè)降水?dāng)?shù)據(jù)分別高：3.47 mm，3.14 mm，1.91 mm，3.73 mm和3.54 mm，而在降水量相對(duì)較多的4月，5月，9月，10月，TRMM降水略低于實(shí)測(cè)降水。

圖2TRMM降水?dāng)?shù)據(jù)與實(shí)測(cè)月降水?dāng)?shù)據(jù)時(shí)間序列

(2) 季尺度。將TRMM 降水?dāng)?shù)據(jù)與站點(diǎn)實(shí)測(cè)降水的數(shù)據(jù)按季節(jié)劃分(春季：3—5月，夏季：6—8月，秋季：9—11月，冬季：12—2月)，整理為季平均降水?dāng)?shù)據(jù)，并對(duì)各季節(jié)的降水?dāng)?shù)據(jù)分別做散點(diǎn)圖(圖3)和時(shí)間序列圖(圖4)。

圖3TRMM降水?dāng)?shù)據(jù)與實(shí)測(cè)季降水?dāng)?shù)據(jù)散點(diǎn)圖

在圖3中春、夏、秋、冬四季的R2分別為0.71，0.80，0.67，0.61，經(jīng)計(jì)算，δ分別為13.23%，10.57%，13.54%和62.62%，RMSE分別為6.08 mm，9.19 mm，5.77 mm和10.82 mm。說明夏季相關(guān)性最好，冬季相關(guān)性最差且誤差較大。

從圖5中可以看出，春夏秋冬四季的TRMM降水與實(shí)測(cè)降水變化趨勢(shì)基本保持一致，其中夏季和冬季降水的變化趨勢(shì)的一致性優(yōu)于春季和秋季，但夏季和冬季TRMM降水明顯高于實(shí)測(cè)降水，尤其冬季差值較大，平均高10.25 mm。

(3) 年尺度。將TRMM 降水?dāng)?shù)據(jù)與站點(diǎn)實(shí)測(cè)降水的數(shù)據(jù)統(tǒng)計(jì)為年降水?dāng)?shù)據(jù)，分別做散點(diǎn)圖(圖5)和時(shí)間序列圖(圖6)，在年時(shí)間尺度上對(duì)TRMM降水的精度進(jìn)行對(duì)比。

圖5中，TRMM多年平均降水?dāng)?shù)據(jù)與實(shí)測(cè)值的R2=0.70，δ=44.23%，RMSE=26.5 mm，可以看出，TRMM年降水?dāng)?shù)據(jù)與雨量站觀測(cè)數(shù)據(jù)之間存在線性相關(guān)性，但誤差較大。從圖6中可以看出，2000—2014年，TRMM年均降水?dāng)?shù)據(jù)與實(shí)測(cè)降水量總體變化趨勢(shì)還是一致的，但同一年份降水量差異較大，除2001年TRMM年均降水比站點(diǎn)實(shí)測(cè)降水低9.3 mm外，其余年份均為TRMM降水?dāng)?shù)據(jù)高，且平均高24.73 mm。

圖4TRMM降水?dāng)?shù)據(jù)與實(shí)測(cè)季降水?dāng)?shù)據(jù)時(shí)間序列

從研究區(qū)的月、季、年尺度上可以看出，TRMM降水與站點(diǎn)實(shí)測(cè)的降水都有很好的一致性，可以大致反映天山山區(qū)降水的分布特征。但大多數(shù)情況下，尤其是降水較少的月份，TRMM降水產(chǎn)品高估了實(shí)際降水，并且TRMM降水?dāng)?shù)據(jù)與實(shí)測(cè)降水?dāng)?shù)據(jù)的決定系數(shù)(R2)隨時(shí)間尺度的增大而減小，相對(duì)誤差(δ)和均方根誤差(RMSE)則隨之增大，說明TRMM遙感數(shù)據(jù)的精度隨時(shí)間尺度的增加而減小。

3.2 校正后精度檢驗(yàn)

以25個(gè)站點(diǎn)的實(shí)測(cè)降水?dāng)?shù)據(jù)為準(zhǔn)，分別從月、季、年3個(gè)時(shí)間尺度對(duì)校正后的TRMM數(shù)據(jù)進(jìn)行精度檢驗(yàn)。

(1) 月尺度。提取每個(gè)站點(diǎn)對(duì)應(yīng)位置校正后的TRMM月降水?dāng)?shù)據(jù)與實(shí)測(cè)月降水量分別構(gòu)建一元線性回歸模型(圖7)，并對(duì)兩組數(shù)據(jù)做時(shí)間序列圖(圖8)。

圖5TRMM降水?dāng)?shù)據(jù)與實(shí)測(cè)年降水?dāng)?shù)據(jù)散點(diǎn)圖

在校正后的TRMM降水與實(shí)測(cè)月降水?dāng)?shù)據(jù)的一元線性回歸分析(圖7)中，R2=0.97，δ=9.47%，RMSE=1.85 mm。可以看出，校正后的TRMM降水?dāng)?shù)據(jù)與實(shí)測(cè)月降水?dāng)?shù)據(jù)擬合度很高，誤差較小。

圖6TRMM降水?dāng)?shù)據(jù)與實(shí)測(cè)年降水?dāng)?shù)據(jù)時(shí)間序列

圖7校正后的TRMM降水?dāng)?shù)據(jù)與實(shí)測(cè)月降水?dāng)?shù)據(jù)散點(diǎn)圖

圖8中，校正后的TRMM月降水?dāng)?shù)據(jù)與站點(diǎn)實(shí)測(cè)多年平均月降水?dāng)?shù)據(jù)的總體變化趨勢(shì)是一致的。經(jīng)計(jì)算，校正后的TRMM月降水?dāng)?shù)據(jù)與實(shí)測(cè)降水?dāng)?shù)據(jù)的多年月平均差值為1.32 mm，在降水量較少、TRMM數(shù)據(jù)與實(shí)測(cè)降水?dāng)?shù)據(jù)誤差較大的1—3月、11月，12月分別差0.74 mm，0.77 mm，0.99 mm，1.55 mm和0.97 mm。

圖8校正后的TRMM降水?dāng)?shù)據(jù)與實(shí)測(cè)月降水?dāng)?shù)據(jù)時(shí)間序列

(2) 季尺度。將兩組校正后的TRMM降水?dāng)?shù)據(jù)與站點(diǎn)實(shí)測(cè)降水?dāng)?shù)據(jù)整理為季降水?dāng)?shù)據(jù)，并對(duì)各季節(jié)的降水?dāng)?shù)據(jù)做散點(diǎn)圖(圖9)和時(shí)間序列圖(圖10)。

在圖9中，春夏秋冬4個(gè)季節(jié)，校正后的TRMM降水?dāng)?shù)據(jù)與實(shí)測(cè)降水?dāng)?shù)據(jù)的R2分別為0.84，0.94，0.79，0.74，δ分別為6.87%，2.47%，11.02%，11.89%，RMSE分別為3.03 mm，2.58 mm，4.4 mm，2.41 mm。

圖10顯示，校正后的TRMM季降水?dāng)?shù)據(jù)更接近實(shí)測(cè)值，冬季的一致性更明顯。經(jīng)過計(jì)算可知，春夏秋冬4個(gè)季節(jié)校正后的TRMM降水與實(shí)測(cè)降水?dāng)?shù)據(jù)的平均差值分別為：2.67 mm，1.91 mm，3.78 mm和1.88 mm。

圖9校正后的TRMM降水?dāng)?shù)據(jù)與實(shí)測(cè)季降水?dāng)?shù)據(jù)散點(diǎn)圖

圖10校正后的TRMM降水?dāng)?shù)據(jù)與實(shí)測(cè)季降水?dāng)?shù)據(jù)時(shí)間序列

(3) 年尺度。將校正后的TRMM降水?dāng)?shù)據(jù)與站點(diǎn)實(shí)測(cè)降水?dāng)?shù)據(jù)統(tǒng)計(jì)為年降水?dāng)?shù)據(jù)，分別做散點(diǎn)圖(圖11)和時(shí)間序列圖(圖12)。

圖11校正后的TRMM降水?dāng)?shù)據(jù)與實(shí)測(cè)年降水?dāng)?shù)據(jù)散點(diǎn)圖

圖12校正后的TRMM降水?dāng)?shù)據(jù)與實(shí)測(cè)年降水?dāng)?shù)據(jù)時(shí)間序列

圖11中，校正后的TRMM降水?dāng)?shù)據(jù)與實(shí)測(cè)年降水?dāng)?shù)據(jù)的R2=0.83，相關(guān)性相對(duì)于月和季尺度較小。δ=4.16%，RMSE=9.1 mm。從圖12可以看出，校正后的TRMM年降水與站點(diǎn)實(shí)測(cè)降水有較好的一致性，但仍存在一定誤差。計(jì)算得知，校正后的TRMM年降水?dāng)?shù)據(jù)比實(shí)測(cè)降水?dāng)?shù)據(jù)平均高6.23 mm，其中2006年，2011—2014年誤差較大，分別高13.28 mm，10.25 mm，18.33 mm，13.06 mm和24.94 mm。

4 結(jié) 論

(1) 校正前，TRMM降水與站點(diǎn)實(shí)測(cè)降水有很好的一致性，TRMM月降水?dāng)?shù)據(jù)與實(shí)測(cè)值的決定系數(shù)(R2)為0.87，存在顯著的線性相關(guān)關(guān)系，說明TRMM數(shù)據(jù)可以反映出研究區(qū)大致的降水情況。相對(duì)誤差(δ)和均方根誤差(RMSE)分別為22.40%和3.89 mm，說明TRMM降水與實(shí)測(cè)降水之間的誤差較大。

(2) TRMM 數(shù)據(jù)與實(shí)測(cè)降水在月、季、年尺度上的決定系數(shù)(R2)分別為：0.87，0.75，0.70，隨時(shí)間尺度的增大而減小，相對(duì)誤差(δ)和均方根誤差(RMSE)則隨之增大，說明TRMM遙感數(shù)據(jù)的精度隨時(shí)間尺度的增加而減小。這與王曉杰研究天山山區(qū)TRMM降水得到的結(jié)論一致。

(3) 校正后，降水量精度得到了顯著地提高。校正后的降水量與實(shí)測(cè)值在月、季、年尺度上的決定系數(shù)(R2)分別為：0.97，0.87，0.83，相對(duì)誤差(δ)分別為：9.47%，8.06%和4.16%，均方根誤差(RMSE)分別為：1.85 mm，3.11 mm和9.1 mm。相比校正前R2提高了10%以上，誤差也有明顯的減小。說明基于重建時(shí)間序列的NDVI數(shù)據(jù)和利用DEM數(shù)據(jù)提取坡度、坡向?qū)RMM數(shù)據(jù)進(jìn)行校正的方法可行。

[1] 劉奇,傅云飛.基于TRMM/TMI的亞洲夏季降水研究[J].中國(guó)科學(xué)D輯:地球科學(xué),2007,37(1):111-122.

[2] 朱會(huì)義,賈紹鳳.降雨信息空間插值的不確定性分析[J].地理科學(xué)進(jìn)展,2004,23(2):34-42.

[3] 吳雪嬌,楊梅學(xué),吳洪波,等. TRMM多衛(wèi)星降水?dāng)?shù)據(jù)在黑河流域的驗(yàn)證與應(yīng)用[J].冰川凍土,2013,35(2):310-319.

[4] 穆振俠,姜卉芳,劉豐.基于TRMM/TMI與實(shí)測(cè)站點(diǎn)的降水垂直分布差異性探討[J].干旱區(qū)研究,2010,27(4):515-521.

[5] 白愛娟,方建剛,張科翔. TRMM衛(wèi)星資料對(duì)陜西及周邊地區(qū)夏季降水的探測(cè)[J].災(zāi)害學(xué),2008,23(2):41-45.

[6] 季漩,羅毅. TRMM降水?dāng)?shù)據(jù)在中天山區(qū)域的精度評(píng)估分析[J].干旱區(qū)地理,2013,26(2):253-262.

[7] 張松林. CART-分類與回歸樹方法介紹[J].火山地質(zhì)與礦產(chǎn),1997,18(1):67-75.

[8] 劉少華,嚴(yán)登華,史曉亮,等.中國(guó)植被NDVI與氣候因子的年際變化及相關(guān)性研究[J].干旱區(qū)地理,2014,37(3):480-489.

[9] Immerzeel W W, Pellicciotti F, Shrestha A B. Glaciers as a proxy to quantify the spatial distribution of precipitation in the Hunza basin[J]. Mountain Research and development, 2012,32(1):30-38.

[10] 張曉.天山中段衛(wèi)星降水?dāng)?shù)據(jù)的降尺度研究及時(shí)空特征分析[D].蘭州:西北師范大學(xué),2013.

[11] Jia S, Zhu W, LüA, et al. A statistical spatial downscaling algorithm of TRMM precipitation based on NDVI and DEM in the Qaidam Basin of China[J]. Remote Sensing of Environment, 2011,115(12):3069-3079.

[12] Cihlar J, Ly H, Li Z, et al. Multitemporal, multichannel AVHRR data sets for land biosphere studies： artifacts and corrections[J]. Remote Sensing of Environment, 1997,60(1):35-57.

[13] 顧娟,李新,黃春林. NDVI時(shí)間序列數(shù)據(jù)集重建方法述評(píng)[J].遙感技術(shù)與應(yīng)用,2006,21(4):391-395.

[14] 王丹,姜小光,唐伶俐,等.利用時(shí)間序列傅里葉分析重構(gòu)無云NDVI圖像[J].國(guó)土資源遙感,2005,64(2):29-32.

[15] 藍(lán)永超,吳素芬,鐘英君,等.近50年來新疆天山山區(qū)水循環(huán)要素的變化特征與趨勢(shì)[J].山地學(xué)報(bào),2007,25(2):177-183.

[16] 藍(lán)永超,吳素芬,韓萍,等.全球變暖情境下天山山區(qū)水循環(huán)要素變化的研究[J].干旱區(qū)資源與環(huán)境,2008,22(6):99-104.

[17] Rozante J R, Cavalcanti I F A. Regional Eta model experiments:SALLJEX and MCS development[J]. Journal of Geophysical Research:Atmospheres, 2008,113(D17):1-20.

ResearchonAccuracyValidationandCalibrationMethodsofTRMMinTianshanMountainsofXinjiang

LI Hui1,4, YANG Tao1,2, HE Qisheng3, REN Weiwei1,4

(1.StateKeyLaboratoryofOasisEcologyandDesertEnvironment,XinjiangInstituteofEcologyandGeography,CAS,Urumqi830011,China; 2.StateKeyLaboratoryofHydrologyWaterResourcesandHydraulicsEngineering,HohaiUniversity,Nanjing210098,China; 3.EarthScienceandEngineering,HohaiUniversity,Nanjing210098,China; 4.UniversityofChineseAcademyofSciences,Beijing100049,China)

Based on the NDVI data after time series reconstruction, we used DEM data to extract the geographical factors, such as slope and aspect, using the method of CART, and mainly calibrated the TRMM3B43 monthly precipitation data in the Tianshan Mountains lacking of observed precipitation data. We used the 25 stations precipitation data to test the precision of TRMM precipitation data and TRMM precipitation data after calibration, respectively. The results showed that TRMM monthly precipitation had a good consistency with observed data，there was a significant linear correlation relationship, but had a the significant deviation. The coefficient of determination (R2) between TRMM precipitation and observed precipitation decreased with time, the relative error (δ) and root mean square error (RMSE) increased, illustrating that the precision of remote sensing data of TRMM decreased with time; the accuracy of TRMM precipitation after calibration had been significantly improved,R2were 0.97, 0.87 and 0.83, which was between TRMM precipitation and observed precipitation on different temporal scales of month, season and year, and increased by more than 10% compared to the previous correction, deviation had decreased obviously. The results proved that the calibration method of CART to the TRMM precipitation data in Tianshan Mountains was feasible.

Tianshan Mountains; TRMM; calibration; CART; accuracy validation

K903

A

1005-3409(2017)06-0327-07

2016-03-03

2016-04-23

國(guó)家自然科學(xué)基金面上項(xiàng)目(41371051);中科院重點(diǎn)部署課題(KZZD-EW-12-1)；荒漠與綠洲生態(tài)國(guó)家重點(diǎn)實(shí)驗(yàn)室基金(G2014-02-04)

李慧(1991—),女,陜西榆林人,碩士研究生,研究方向?yàn)樗呐c水資源.E-mail:527919855@qq.com