離散Landsat歸一化植被指數(shù)的時間序列融合構(gòu)建

曾德裕 盧曉寧 黃玥 楊柳青 孟成真

摘要：結(jié)合滑動窗的使用，提出了一種基于傅里葉逼近和線性融合降尺度的時間序列構(gòu)建方法。在時間序列上逐像元對離散的Landsat數(shù)據(jù)進行傅里葉逼近，并同分辨率較低的AVHRR傳感器數(shù)據(jù)進行融合，實現(xiàn)了試驗區(qū)Landsat植被指數(shù)的時間序列構(gòu)建，結(jié)合作為實際觀測數(shù)據(jù)的Landsat數(shù)據(jù)對模型結(jié)果的精度進行驗證。結(jié)果表明，模型結(jié)果在色調(diào)與紋理等圖像特征方面與實際觀測數(shù)據(jù)一致，構(gòu)建了白色冰雪覆蓋區(qū)域的NDVI，構(gòu)建結(jié)果與周邊地表具有非常好的空間連續(xù)性，且整幅影像不存在明顯的塊狀效應(yīng);模型結(jié)果與實際觀測數(shù)據(jù)具有較高的相關(guān)性（平均0.869 2）和較低的均方根誤差（平均0.043 5），平均高達75.04%的像元誤差在±0.05之間，高達97.64%的像元誤差在±0.1之間，并且誤差呈良好的正態(tài)分布。模型精度對輸入數(shù)據(jù)量存在一定依賴，同時對混合像元的處理方面尚存在一定的改進空間。

關(guān)鍵詞：Landsat;NDVI;數(shù)據(jù)融合;傅里葉分析

中圖分類號：TP79? ? ? ? ?文獻標識碼：A

文章編號：0439-8114（2019）12-0045-07

DOI：10.14088/j.cnki.issn0439-8114.2019.12.012? ? ? ? ? ?開放科學（資源服務(wù)）標識碼（OSID）：

Abstract： In combination with the use of sliding windows， a time series construction method based on Fourier approximation and linear fusion downscaling is proposed. The Fourier approximation of discrete Landsat data is performed on a time-by-pixel basis and merged with the lower resolution AVHRR sensor data to realize the time series construction of the Landsat vegetation index in experimental area， combined with Landsat as the actual observation data. The data verifies the accuracy of the model results. The results show that： The model results in the hue and texture， image features are consistent with the actual observation data， and the NDVI of the white snow-covered area is constructed. The construction results have very good spatial continuity with the surrounding surface， and there is no obvious blocky effect in the entire image; Model results and actual observed data has a high correlation （mean 0.869 2） and a low root mean square error （average of 0.043 5）， averaging up to 75.04% of the pixels between ±0.05， up to 97.64% of the pixel error between ±0.1， and the error is in a good normal distribution. The accuracy of the model has a certain dependence on the amount of input data， and there is still room for improvement in the processing of mixed pixels.

Key words： Landsat; NDVI; data fusion; Fourier analysis

紅光和紅外波段包含了90%以上的植被信息[1]，被廣泛用于植被指數(shù)的構(gòu)建，并基于此開展氣候、水文、土壤等環(huán)境背景信息的研究和探討，植被指數(shù)在區(qū)域和全球尺度植被生態(tài)系統(tǒng)空間分布和動態(tài)變化研究中得到了廣泛應(yīng)用[2]。長時間序列的植被指數(shù)能反映植被在時間維上的生長變化，對區(qū)域變化研究的價值不言而喻，在作物的長勢監(jiān)測[3]、作物估產(chǎn)[4，5]、作物分類[6，7]等方面均具有重要作用。時間序列長度逾40年的Landsat衛(wèi)星數(shù)據(jù)以其優(yōu)越穩(wěn)定的性能、較高的空間分辨率和超長的時間序列成為區(qū)域植被長期動態(tài)變化研究的首選。但由于云雪等因素的客觀存在及16 d較長的重訪周期，使得Landsat數(shù)據(jù)往往大量缺失，造成數(shù)據(jù)離散程度偏高，大大限制了相關(guān)研究工作的開展。構(gòu)建時間序列Landsat植被指數(shù)數(shù)據(jù)具有重要的現(xiàn)實意義。

當前，時間序列植被指數(shù)的構(gòu)建主要集中于對數(shù)據(jù)進行濾波降噪處理，常用的方法包括最大值合成法（MVC）、最佳指數(shù)斜率提取法（BISE）、Savitzky-Golay濾波（S-GF）、時間窗口的線性內(nèi)插算法（TWO）、基于非對稱的高斯函數(shù)擬合法（AGFF）、基于離散傅里葉分析的系列算法等[8]。然而，若將Landsat數(shù)據(jù)在時間序列上的大量無規(guī)律數(shù)據(jù)缺失亦視為一種“噪聲”的話，則上述方法多數(shù)難以應(yīng)對。近年來，通過融合多源數(shù)據(jù)獲得較高空間分辨率影像研究取得較大進展，其中，Gao等[9]綜合考慮不同時期高低分辨率數(shù)據(jù)的時間、空間和光譜差異，提出了時空適應(yīng)反射率融合模型（Spatial and temporal adaptive reflectance fusion model，STARFM），并利用MODIS和Landsat數(shù)據(jù)進行了類似Landsat數(shù)據(jù)的生產(chǎn)，取得了較好的結(jié)果。之后，出現(xiàn)了大量改進算法，以適應(yīng)復(fù)雜異質(zhì)性區(qū)域[10，11]、特定時間[12]以及進行時間序列數(shù)據(jù)生產(chǎn)[13]等的研究需要，新近又出現(xiàn)了與混合像元分解相結(jié)合的混合算法[14，15]等。從使用目的來看，由單純地進行短波反射率數(shù)據(jù)的融合，發(fā)展為用于土壤濕度[16]、植被指數(shù)[17，18]、地表溫度[19]等地表參數(shù)的定量數(shù)據(jù)融合降尺度研究。然而，多數(shù)算法僅適用于年內(nèi)預(yù)測，年際數(shù)據(jù)融合往往存在較大誤差。基于當前算法進行離散Landsat數(shù)據(jù)的時間序列構(gòu)建普遍面臨較大困難。

為此，本研究提出了一種基于傅里葉逼近和線性融合降尺度的Landsat歸一化植被指數(shù)（NDVI）時間序列的構(gòu)建方法。該方法首先在時間序列上逐像元對離散的Landsat數(shù)據(jù)進行傅里葉逼近，隨后，結(jié)合滑動窗的使用，融合分辨率較低的AVHRR傳感器數(shù)據(jù)，并結(jié)合實際的Landsat數(shù)據(jù)對結(jié)果精度進行驗證，目的是實現(xiàn)長時間序列的、相對高空間分辨率的植被指數(shù)的時間序列融合構(gòu)建。

1? 算法及流程

算法主要包括兩個部分，即離散高分辨率數(shù)據(jù)的傅里葉逼近和逼近結(jié)果與粗分辨率數(shù)據(jù)的線性融合，同時，為避免結(jié)果中出現(xiàn)反映粗分辨率像元特征的斑塊效應(yīng)，設(shè)計了一個滑動窗，使數(shù)據(jù)融合過程在窗口內(nèi)進行。鑒于算法主要面向長時間序列數(shù)據(jù)的融合且為線性融合方法，故稱之為長時間序列線性融合模型（Long timeseries linear fousion model，LTLFM）。

1.1? 傅里葉（Fourier）逼近

時間序列NDVI的傅里葉構(gòu)建的實質(zhì)是把時間波譜數(shù)據(jù)分解成許多不同頻率的正弦曲線和余弦曲線，從中選取若干能夠反映時間序列特征的曲線進行疊加，以達到時間序列數(shù)據(jù)的構(gòu)建目的。構(gòu)建過程中，諧波個數(shù)的選擇十分關(guān)鍵，一方面要考慮觀測數(shù)據(jù)（樣本）量，同時需考慮如何能夠更為準確地逼近觀測值。一般而言，增加諧波個數(shù)可獲得更低的擬合偏差，擬合線更靠近觀測值，但同時也會在時間序列曲線上出現(xiàn)更多的起伏，產(chǎn)生無法解釋的“規(guī)律”[20]。因此，為保證對NDVI的周期和擾動的準確模擬和逼近，同時避免產(chǎn)生不可解釋的異常逼近結(jié)果，結(jié)合Landsat數(shù)據(jù)量，選定8個參數(shù)的傅里葉擬合公式，公式[12]如下：

式中，x為朱利安日期（Julian date，JD，儒略日）;T為每年的天數(shù)（T=365.25）;a0用于對NDVI的整體值進行估計;a1、a2、a3和b1、b2、b3分別用于對NDVI的年內(nèi)和年際變化進行估計;c1用于對NDVI隨時間變化的線性趨勢進行估計，該部分的引入使上述模型一定程度上區(qū)別于一般的諧波分析模型。

1.2? 數(shù)據(jù)融合

根據(jù)線性混合理論，同時考慮不同傳感器之間的光譜響應(yīng)差異、波段劃分以及角度和過境時間等差異所造成的系統(tǒng)性誤差，則有：

式中，Ct為低分辨率像元反射率;a為調(diào)整系數(shù);FFti為高分辨率均質(zhì)像元表面反射率;Ati是面積比例;i表示高分辨率像元的空間索引（位置）;？著k為系統(tǒng)性誤差。假設(shè)各像元面積比例相等，則相應(yīng)地有：

式中，F(xiàn)ti為融合結(jié)果;N為與粗分辨率像元相對應(yīng)的高分辨率像元數(shù)。

1.3? 滑動窗的設(shè)計

線性光譜混合模型具有構(gòu)模簡單、物理含義明確的優(yōu)點，理論上也有較好的科學性，但線性混合模型適用于本質(zhì)上屬于或者基本屬于線性混合的地物，以及在大尺度上可以認為是線性混合的地物[21，22]。因此，線性融合會無可避免地在結(jié)果中產(chǎn)生反映低分辨率像元特征的塊狀效應(yīng)。滑動窗的設(shè)計能夠從表面上、視覺效果上消除這種塊狀效應(yīng)，實質(zhì)是對光譜非線性混合的一種近似。

本研究設(shè)計了一個與低分辨率像元大小相同的滑動窗，在窗口進行線性降尺度。該窗口每次滑動一個細分辨率像元并進行一次降尺度運算，將結(jié)果累積到相應(yīng)高分辨率像元，并求平均，得到最終結(jié)果。一個低分辨率像元含n個高分辨率像元的情況下，每個高分辨率像元參與n×n次線性降尺度運算。

滑動窗設(shè)計的關(guān)鍵在于窗口值的計算。如圖1所示，該窗口最多涉及4個粗分辨率像元，相應(yīng)地窗口與粗分辨率的交叉部分最多可為4部分，窗口值W=Part1+Part2+Part3+Part4，各部分根據(jù)其像元值所占比例和所對應(yīng)的粗分辨率像元值進行計算，公式如下：

式中，Lpi為第i部分的細分辨率像元值（之和）;Lfi為第i部分對應(yīng)粗分辨率像元所覆蓋的細分辨率像元值（之和）;Gi為該部分所對應(yīng)的粗分辨率像元值。當該部分所在粗分辨率像元對應(yīng)的細分辨率像元有效像元數(shù)不足80%時，則該部分的值不再根據(jù)像元值比例計算，而是根據(jù)像元數(shù)比例，結(jié)合所對應(yīng)的粗分辨率像元值進行計算。滑動過程中，還可能僅包含一或兩個交叉部分，計算原理相同。

2? 試驗區(qū)及數(shù)據(jù)處理

2.1? 試驗區(qū)概況

試驗區(qū)位于青藏高原東緣，覆蓋該區(qū)域的Landsat衛(wèi)星行列號為H131V37?？偯娣e1.849 4萬km2，土地利用類型以草地為主（15 045.06 km2，占試驗區(qū)總面積的81.35%），其次是林地（908.50 km2，占試驗區(qū)總面積的4.91%），并有少量耕地（92.19 km2，占試驗區(qū)總面積的0.50%）。海拔高度位于3 081～4 371 m，地勢西南高東北低，在平坦低洼的東北部形成大片沼澤濕地（面積2 041 km2，占試驗區(qū)總面積的11.04%）。該區(qū)為典型高寒濕地生態(tài)系統(tǒng)，主要包括若爾蓋濕地核心區(qū)及其南部、西部及西南部山地區(qū)域。黃河及黑河、白河兩條支流貫穿其中。氣候寒冷濕潤，東北部地勢平坦，河曲發(fā)育，湖泊眾多，排水不暢，泥炭沼澤廣泛發(fā)育，年平均氣溫在0 ℃左右，多年平均降水量500～600 mm，蒸發(fā)量小于降水量，地表經(jīng)常處于過濕狀態(tài)。

2.2? 數(shù)據(jù)源及數(shù)據(jù)處理

2.2.1? Landsat數(shù)據(jù)及其處理? Landsat數(shù)據(jù)下載自地理空間數(shù)據(jù)云（http：//www.gscloud.cn/），覆蓋該區(qū)域的影像行列號為H131V37，空間分辨率為30 m。篩選出其中有一定面積的無云覆蓋區(qū)域的影像，共計47景，數(shù)據(jù)日期如下（格式：年+年內(nèi)序）：

1990189、1991176、1993245、1994184、1994216、1994248、1995107、2000265、2000281、2000313、2001187、2001203、2001283、2002014、2002046、2002350、2004260、2004292、2004308、2005006、2005038、2005086、2005150、2006217、2006281、2007140、2007204、2007220、2007268、2008095、2008127、2008175、2008191、2008239、2008287、2008303、2008335、2008351、2009049、2009097、2009177、2009209、2009225、2009337、2010212、2010260、2011199。

對數(shù)據(jù)進行輻射定標、大氣校正、云檢測及標記等處理，并計算歸一化植被指數(shù)，為節(jié)約計算時間，使用像元聚合方式（Pixel aggregate）將其重采樣到與MODIS數(shù)據(jù)相同的250 m分辨率。歸一化植被指數(shù)的計算公式為NDVI=（NIR-R）/（NIR+R），即近紅外與紅光波段之差同近紅外與紅光波段之和的比值。

2.2.2? AVHRR數(shù)據(jù)及其處理? AVHRR是搭載在NOAA氣象衛(wèi)星上的傳感器，目前有多種AVHRR數(shù)據(jù)集，在AVHRR NDVI數(shù)據(jù)集中，GIMMS NDVI最能反映植被動態(tài)變化[23]，也是目前數(shù)據(jù)系列最長、數(shù)據(jù)最完整且應(yīng)用最廣泛的AVHRR時間序列數(shù)據(jù)集[24]，因此本研究選擇使用該數(shù)據(jù)集。該數(shù)據(jù)集來自美國宇航局艾姆斯研究中心生態(tài)預(yù)測實驗室數(shù)據(jù)網(wǎng)站（http：//ecocast.arc.nasa.gov/），為GIMMS NDVI 3G數(shù)據(jù)，時長1981—2013年，為半月（15 d）數(shù)據(jù)，全年共24期，累計780景，空間分辨率為8 km。

對其處理過程包括定義投影、格式轉(zhuǎn)換、數(shù)據(jù)裁剪等操作，同時，為消除云與冰雪等的影響，使用Savitzky-Golay濾波算法對其進行濾波重構(gòu)，得到高質(zhì)量的GIMMS NDVI數(shù)據(jù)。

2.2.3? MODIS數(shù)據(jù)及其處理? MODIS數(shù)據(jù)主要用來從時間序列上對模型結(jié)果進行驗證。使用的是16 d合成、空間分辨率為250 m的MOD13Q1產(chǎn)品，來自美國宇航局LAADS（Level 1 and atmosphere archive and distribution system）數(shù)據(jù)網(wǎng)站（https：//ladsweb.nascom.nasa.gov/）。

數(shù)據(jù)的處理過程在MRT中進行，首先對其進行投影轉(zhuǎn)換和格式轉(zhuǎn)換，并裁剪出試驗區(qū)。同樣地，為消除云與冰雪等的影響，使用Savitzky-Golay濾波算法對其進行濾波重構(gòu)。

2.2.4? 數(shù)據(jù)一致性處理? 由于傳感器及數(shù)據(jù)本身性質(zhì)不同（Landsat為瞬時觀測數(shù)據(jù)，GIMMS為半月合成數(shù)據(jù)），Landsat和GIMMS NDVI難免存在時空上的差異。考慮到Landsat具有更高的數(shù)據(jù)質(zhì)量，故以Landsat為優(yōu)化目標，作為實際觀測數(shù)據(jù)，在GIMMS空間分辨率尺度，基于各像元時間序列數(shù)據(jù)，建立GIMMS與Landsat數(shù)據(jù)之間的線性關(guān)系：Landsat=a+b×GIMMS，并基于回歸系數(shù)，對GIMMS數(shù)據(jù)進行校正，得到相對一致的多源數(shù)據(jù)。

3? 結(jié)果與分析

使用8參數(shù)（式1中的a0、a1、a2、a3、b1、b2、b3、c1這8個參數(shù)）的傅里葉公式對時間序列上離散的Landsat NDVI進行擬合，并根據(jù)擬合參數(shù)公式，以半月步長構(gòu)建了1981—2013年時間序列上試驗區(qū)Landsat NDVI，之后同GIMMS NDVI數(shù)據(jù)進行融合，得到最終的LTLFM模型結(jié)果。隨機選取多期實際Landsat NDVI數(shù)據(jù)，分別從圖像特征、差值直方圖、散點圖等對模型結(jié)果進行空間上的驗證，并結(jié)合所有Landsat數(shù)據(jù)，從時間序列上對模型結(jié)果進行驗證。

3.1? LTLFM模型結(jié)果的空間特征

不同季節(jié)各隨機選取了一景結(jié)果與實際Landsat NDVI數(shù)據(jù)以圖像形式進行對比，結(jié)果見圖2。模型結(jié)果中，河流、山脈等地物紋理特征明顯，且與實際觀測數(shù)據(jù)一致，色調(diào)均勻，與實際觀測數(shù)據(jù)較為接近;同時，構(gòu)建了白色冰雪覆蓋區(qū)域的NDVI，且原冰雪覆蓋區(qū)的NDVI模擬結(jié)果與周邊地表具有非常好的空間連續(xù)性，整幅影像不存在明顯的塊狀效應(yīng)。視覺上分季節(jié)來看，夏秋季略好于冬春季?？傮w上反映了較高的精度。

為更廣泛和定量地考察模型結(jié)果的精度，除上述4期數(shù)據(jù)之外，分別在各季節(jié)補充選取兩期共計12期數(shù)據(jù)，對結(jié)果與實際觀測數(shù)據(jù)的散點圖及差值直方圖進行分析。由圖3和圖4可知，各期LTLFM模型結(jié)果與實際觀測數(shù)據(jù)的相關(guān)系數(shù)絕大多數(shù)（83%）在0.80以上，最高0.920 4，最低0.789 3，平均為0.869 2，反映出較好的相關(guān)性;斜率在0.623 9～1.065 7，平均為0.866 1，50%以上在1.0±0.1，反映出模擬數(shù)據(jù)與實際觀測數(shù)據(jù)空間上有較好的一致性，同時存在一定程度的（總體<20%）低值高估或高值低估現(xiàn)象。

由圖4可知，根據(jù)LTLFM模型模擬出的數(shù)據(jù)與實際觀測數(shù)據(jù)也存在一定的季節(jié)性。相對而言，夏秋季一致性較好，平均相關(guān)系數(shù)R2分別為0.893 9和0.886 8，冬春季則分別為0.844 3和0.851 9。同時平均斜率也更接近1，其中，夏秋季分別為1.002 5和0.887 3，而冬春季則分別為0.810 7和0.859 3，秋季略好于春季而差于夏季。實際上這種情況的出現(xiàn)主要與輸入數(shù)據(jù)及其處理過程有關(guān)：①春夏秋冬各季輸入數(shù)據(jù)分別為7、22、13和5期，在相關(guān)時間點上形成控制，降低了傅里葉擬合結(jié)果在相應(yīng)時段的不確定性;②試驗區(qū)位于青藏高原東緣，冬春季部分區(qū)域積雪積冰較多，有效觀測數(shù)據(jù)相對較少，從而導(dǎo)致這些區(qū)域NDVI的傅里葉擬合精度降低，這表明要獲得更好的模擬結(jié)果仍有賴于更多的有效輸入數(shù)據(jù)，且只有達到一定的數(shù)據(jù)量，才能夠保證結(jié)果的精度;③斜率的季節(jié)性特征反映在單期圖像上，相對于實際情況，存在低值的線性高估或高值的線性低估現(xiàn)象，而根據(jù)散點圖來看，主要是低值的高估，這種情況的存在很大程度上是由于冰雪、云及其陰影的不徹底去除和標記所帶來的影響，如能有效加強這一環(huán)節(jié)，將使模型結(jié)果本身的精度及驗證的結(jié)果精度大大提高。

各期數(shù)據(jù)的均方根誤差（sx）在0.025 4～0.064 0，平均0.043 5，其中，平均高達75.04%的像元在±0.05之間，高達97.64%的像元誤差在±0.1之間，并且誤差呈良好的正態(tài)分布特征（圖5）。模型結(jié)果與實際觀測數(shù)據(jù)像元數(shù)峰值對應(yīng)NDVI往往存在一定偏差，主要由如下兩種原因?qū)е拢孩偎惴ㄒ暩叻直媛氏虻头直媛蔔DVI的轉(zhuǎn)換為線性變換，對高低分辨率數(shù)據(jù)之間總體系統(tǒng)誤差的估計尚存在一定不足;②Landsat為瞬時觀測數(shù)據(jù)，重訪周期為16 d，而GIMMS為最大值合成產(chǎn)品，每半月一期，這種時間上的不完全匹配也是產(chǎn)生誤差的重要原因。因此，由算法及數(shù)據(jù)源本身的時空匹配問題而造成的高低分辨率數(shù)據(jù)之間的系統(tǒng)誤差的估計問題，尚存在一定改進余地，理論上模擬結(jié)果的精度應(yīng)高于實際驗證精度。

3.2? LTLFM模型結(jié)果的時間序列特征

構(gòu)建長時間序列NDVI數(shù)據(jù)集的核心目的就在于研究其在時間上的變化趨勢，反映其在不同時刻的氣候、水文、生態(tài)等環(huán)境背景信息。因此，模擬結(jié)果能否在時間序列上準確反映NDVI實際的變化趨勢，是判斷研究方法是否成功達到預(yù)期的關(guān)鍵。由于圖幅限制，選取自2002年始的Landsat數(shù)據(jù)和時間序列MODIS NDVI數(shù)據(jù)，對所構(gòu)建的時間序列NDVI數(shù)據(jù)集進行驗證（其中Landsat數(shù)據(jù)截至2011年，MODIS數(shù)據(jù)截至2014年）。

圖6分別考察了耕地、林地和草地3種試驗區(qū)主要土地利用類型NDVI模擬結(jié)果同實際觀測數(shù)據(jù)（Landsat、MODIS）在時間序列上的一致性。結(jié)果顯示，模擬數(shù)據(jù)對Landsat數(shù)據(jù)有較好的逼近效果，尤其是耕地和草地，平均絕對值誤差分別僅為0.045 8和0.043 8，而對試驗區(qū)林地的逼近效果相對較差，但其平均絕對值誤差亦僅為0.058 1。相較于MODIS數(shù)據(jù)，各土地利用類型NDVI模擬值系統(tǒng)性偏低，多年平均而言，耕地偏低0.047 8（11.63%），林地偏低0.088 1（17.17%），草地偏低0.010 3（2.48%），但對變化趨勢的反映情況總體較好，尤其是對生長季NDVI的模擬。

試驗區(qū)不同土地利用類型NDVI模擬結(jié)果的上述差異由多種原因?qū)е拢菏紫仁腔旌舷裨獑栴}的存在。試驗區(qū)位于青藏高原東緣，草地和耕地的分布相對集中，而林地的分布則相對離散和破碎，使得在250 m和8 km空間分辨率尺度上均存在大量混合像元，增加了線性融合降尺度模型的不確定性;其次，林地本身較多地集中分布于易形成積雪覆蓋、海拔相對較高的山地區(qū)域，使得有效的Landsat像元本身較少，這也是造成林地NDVI模擬結(jié)果偏差較大的原因。這表明模型雖已具備一定精度，但對混合像元的處理方面仍有一定的提升空間，同時，模型精度對有效輸入數(shù)據(jù)量存在依賴關(guān)系，有效數(shù)據(jù)量越大，模型精度越高。

4? 小結(jié)

結(jié)合滑動窗的使用，提出了一種基于傅里葉逼近和線性融合降尺度的Landsat NDVI時間序列的構(gòu)建方法（LTLFM），在時間序列上逐像元對離散的Landsat數(shù)據(jù)進行傅里葉逼近，并同分辨率較低的AVHRR傳感器數(shù)據(jù)進行了融合，結(jié)合實際的Landsat數(shù)據(jù)對結(jié)果精度進行驗證，反映了較好的結(jié)果。

1）模型結(jié)果在色調(diào)與紋理等圖像特征方面與實際觀測數(shù)據(jù)一致，構(gòu)建了白色冰雪覆蓋區(qū)域的NDVI，構(gòu)建結(jié)果與周邊地表具有非常好的空間連續(xù)性，且整幅影像不存在明顯的塊狀效應(yīng)。

2）模型結(jié)果與實際觀測數(shù)據(jù)具有較高的相關(guān)性（平均0.869 2）和較低的均方根誤差（平均0.043 5），平均高達75.04%的像元誤差在±0.05之間，高達97.64%的像元誤差在±0.1之間，并且誤差呈良好的正態(tài)分布。

3）模型結(jié)果精度呈現(xiàn)出一定的季節(jié)特征，即夏秋季略高于冬春季，主要與有效輸入數(shù)據(jù)量有關(guān)，有效輸入數(shù)據(jù)量越大，模型精度越高;同時在不同土地利用類型方面的表現(xiàn)存在一定差異，即耕地和草地優(yōu)于林地，反映了模型在對混合像元的處理方面尚存在一定的改進空間。

4）與傅里葉逼近方法相比，基于LTLFM模型具有更加清晰的紋理和細部特征等圖像特征，在時間序列上更接近真實，體現(xiàn)出該方法較原方法的優(yōu)越性。這對土地覆蓋的變化、地表植被生態(tài)系統(tǒng)的物候變化等的研究檢測具有重要意義。

要獲得更高精度的植被指數(shù)時間序列構(gòu)建結(jié)果，可從模型和數(shù)據(jù)源兩個方面進行改進：

1）模型改進。傅里葉公式選定了固定個數(shù)的變量，尚不夠靈活，在實際觀測樣本足夠的情況下，可考慮依據(jù)規(guī)則使用變量可變的傅里葉公式實現(xiàn)時間序列NDVI的分段、滑動擬合。線性融合降尺度模型處理均質(zhì)像元效果較好，處理混合像元效果較差，可考慮使用更復(fù)雜的融合模型，但伴隨而來的對土地利用數(shù)據(jù)的需求及其所帶來的不確定性問題，也是需要認真考慮的。

2）數(shù)據(jù)源問題。首先，要獲得更好的模擬結(jié)果有賴于更多的有效輸入數(shù)據(jù)，未來可考慮使用更多其他數(shù)據(jù)，同Landsat一起，提供更多時刻的空間異質(zhì)性信息，使傅里葉逼近結(jié)果更準確。其次，應(yīng)選擇時空匹配較好的高低分辨率數(shù)據(jù)源，或者在數(shù)據(jù)融合之前，更好地解決多源數(shù)據(jù)之間合成與瞬時、過境時間不一致等帶來的時空匹配問題。

參考文獻：

[1] BARET F，GUYOT G，MAJOR D J. TSAVI：A vegetation index which minimizes soil brightness effects on LAI and APAR estimation[A].Geoscience and Remote Sensing Symposium，1989.IGARSS'89.12th Canadian symposium on remote sensing[C].IEEE，1989.1355-1358.

[2] 史? 軍，孫國武，崔林麗.中國東部NDVI時空變化及對東亞夏季風響應(yīng)[J].高原氣象，2008（S1）：170-174.

[3] 鄒文濤，吳炳方，張? 淼，等.農(nóng)作物長勢綜合監(jiān)測——以印度為例[J].遙感學報，2015，19（4）：539-549.

[4] 劉? 佳，王利民，楊福剛，等.基于HJ時間序列數(shù)據(jù)的農(nóng)作物種植面積估算[J].農(nóng)業(yè)工程學報，2015，31（3）：199-206.

[5] 劉煥軍，孟令華，張新樂，等.基于時間序列Landsat影像的棉花估產(chǎn)模型[J].農(nóng)業(yè)工程學報，2015，31（17）：215-220.

[6] 李鑫川，徐新剛，王紀華，等.基于時間序列環(huán)境衛(wèi)星影像的作物分類識別[J].農(nóng)業(yè)工程學報，2013，29（2）：169-176.

[7] 李登科，郭? 鈮.基于植被生長規(guī)律的陜西省植被遙感分類[J].高原氣象，2008，27（1）：215-221.

[8] 李? 儒，張? 霞，劉? 波，等.遙感時間序列數(shù)據(jù)濾波重建算法發(fā)展綜述[J].遙感學報，2009，13（2）：335-341.

[9] GAO F，MASEK J，SCHWALLER M，et al. On the blending of the Landsat and MODIS surface reflectance：Predicting daily Landsat surface reflectance[J].IEEE transactions on geoscience & remote sensing，2006，44（8）：2207-2218.

[10] ZHU X，CHEN J，GAO F，et al. An enhanced spatial and temporal adaptive reflectance fusion model for complex heterogeneous regions[J].Remote sensing of environment，2010，114：2610-2623.

[11] HWANG T，SONG C，BOLSTAD P V，et al. Downscaling real-time vegetation dynamics by fusing multi-temporal MODIS and Landsat NDVI in topographically complex terrain[J].Remote sensing of environment，2011，115：2499-2512.

[12] ZHU Z，WOODCOCK C E，HOLDEN C，et al. Generating synthetic Landsat images based on all available Landsat data：Predicting Landsat surface reflectance at any given time[J].Remote sensing of environment，2015，162：67-83.

[13] BINDHU V M，NARASIMHAN B. Development of a spatio-temporal disaggregation method （DisNDVI） for generating a time series of fine resolution NDVI images[J].ISPRS journal of photogrammetry and remote sensing，2015，101：57-68.

[14] GEVAERT C M，GARC？魱A-HARO F J. A comparison of STARFM and an unmixing-based algorithm for Landsat and MODIS data fusion[J].Remote sensing of environment，2015，156：34-44.

[15] 謝登峰，張錦水，孫佩軍，等.結(jié)合像元分解和 STARFM模型的遙感數(shù)據(jù)融合[J].遙感學報，2016，20（1）：62-72.

[16] ZHAO W，LI A. A downscaling method for improving the spatial resolution of AMSR-E derived soil moisture product based on MSG-SEVIRI data[J].Remote sensing，2013，5（12）：6790-6811.

[17] 蒙繼華，吳炳方，杜? 鑫，等.高時空分辨率NDVI數(shù)據(jù)集構(gòu)建方法[J].遙感學報，2011，15（1）：52-59.

[18] BUSETTO L，MERONI M，COLOMBO R. Combining medium and coarse spatial resolution satellite data to improve the estimation of sub-pixel NDVI time series[J].Remote sensing of environment，2008，112：118-131.

[19] ZHAN W F，CHEN Y H，ZHOU J，et al. Spatial simulation of urban heat island intensity based on the support vector machine technique：A case study in Beijing[J].Acta geodaetica et cartographica sinica，2011，40（1）：96-103.

[20] 林忠輝，莫興國. NDVI時間序列諧波分析與地表物候信息獲取[J].農(nóng)業(yè)工程學報，2006，22（12）：138-144.

[21] KWON H，NASRABADI N M. Kernel orthogonal subspace projection for hyperspectral signal classification[J].IEEE transactions on geoscience and remote sensing，2005，43（12）：2952-2962.

[22] 趙遼英，張? 凱，厲小潤.高光譜圖像目標檢測的核信號空間正交投影法[J].遙感學報，2011，15（1）：13-28.

[23] BECK H E，MCVICAR T R，VAN DIJK A I J M，et al. Global evaluation of four AVHRR-NDVI data sets：Inter-comparison and assessment against Landsat imagery[J].Remote sensing of environment，2011，115：2547-2563.

[24] 馬明國，宋? 怡，王旭峰，等.AVHRR，VEGETATION和MODIS時間系列遙感數(shù)據(jù)產(chǎn)品現(xiàn)狀與應(yīng)用研究進展[J].遙感技術(shù)與應(yīng)用，2012，27（5）：663-670.

收稿日期：2018-09-11

基金項目：四川省教育廳重點（自然科學）項目（17ZA0075）;成都市科技局科技惠民技術(shù)研發(fā)項目（2016-HM01-00392-SF）;國家自然

作者簡介：曾德裕（1994-），男，湖南臨湘人，在讀碩士研究生，研究方向為資源環(huán)境遙感，（電話）18708121831（電子信箱）673526122@qq.com。