田 苗 單 捷 盧必慧 黃曉軍

(江蘇省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)信息研究所， 南京 210014)

0 引言

及時準確獲取水稻種植面積及空間分布，對水稻生產(chǎn)管理、農(nóng)業(yè)政策制定及糧食安全分析等具有重要意義。目前水稻面積的遙感提取方法主要有常規(guī)分類法和基于MODIS的時間序列方法。常規(guī)分類法一般使用中、高分辨率的遙感數(shù)據(jù)[1-3]，而此種數(shù)據(jù)重訪周期長，易受天氣影響，不易抓住水稻生長的物候信息，因而，難以進行大范圍的水稻面積業(yè)務(wù)化運行。MODIS數(shù)據(jù)具有較高的時間分辨率和較廣闊的空間覆蓋，而且可以免費獲取，所以適合大區(qū)域水稻面積提取和推廣應(yīng)用[4-5]?；贛ODIS時間序列的方法往往根據(jù)水稻的物候歷，確定水稻識別的移栽期、生長期和收獲期等關(guān)鍵時期，根據(jù)水稻在不同時期表現(xiàn)的遙感特征指數(shù)差異來識別水稻[6-11]。較少有學(xué)者考慮與水稻生產(chǎn)能力密切相關(guān)的季節(jié)積分和生長季振幅兩個物候指標。經(jīng)過眾多學(xué)者對植被物候的研究，開發(fā)了許多判斷植被關(guān)鍵生育期的方法，主要包括植被指數(shù)閾值法[12]、植被指數(shù)中點法[13]、最大變化斜率法[14]、植被指數(shù)滑動移動平均法[15]、通過NDVI時間序列累積生長天數(shù)來計算物候信息[16]、植被物候期的頻率分布與遙感數(shù)據(jù)相結(jié)合法[17]等。目前圖像處理已經(jīng)進入人工智能時代，機器學(xué)習(xí)回歸算法的高度適用性和通用性使其比傳統(tǒng)的擬合函數(shù)能更精確地估計物候趨勢。

本研究以江蘇省為研究區(qū)域，應(yīng)用機器學(xué)習(xí)算法對MODIS-EVI時間序列進行重構(gòu)，并提取水稻物候信息，選取季節(jié)積分和生長季振幅兩個指標，采用閾值法提取水稻種植面積，并應(yīng)用統(tǒng)計數(shù)據(jù)和Landsat8 OLI影像對提取結(jié)果進行精度驗證。

1 研究區(qū)域與數(shù)據(jù)獲取

1.1 研究區(qū)域概況

江蘇省地處亞熱帶向暖溫帶過渡的區(qū)域，其自然條件對水稻生長有利，是全國水稻主產(chǎn)省之一。江蘇省南北縱跨5個緯度(30°45′～35°20′N)，東西橫跨5個經(jīng)度(116°18′～121°57′E)，其中，北部地處我國秈稻種植最北邊緣，南部地處我國兩熟制粳稻種植最南邊緣，東部受海洋氣候影響，西部受陸地氣候影響，形成了境內(nèi)多種稻作生態(tài)區(qū)、多種類型稻作熟制和稻作方式、秈粳交錯和多種熟制品質(zhì)并存的水稻生產(chǎn)格局，在全國稻作生產(chǎn)中具有特殊意義。

1.2 數(shù)據(jù)獲取與預(yù)處理

本文采用2019—2020年地表反射率數(shù)據(jù)(MOD09A1)作為基礎(chǔ)數(shù)據(jù)，數(shù)據(jù)來源于美國國家航空航天局網(wǎng)站(https:∥modis.gsfc.nasa.gov/)，MOD09A1數(shù)據(jù)的空間分辨率為500 m，時間分辨率為8 d。使用MODIS重投影工具(MODIS Re-projection tool, MRT)對數(shù)據(jù)進行重采樣、格式轉(zhuǎn)換、拼接和重投影等預(yù)處理。采用中高空間分辨率的Landsat8 OLI遙感影像數(shù)據(jù)進行水稻信息目視解譯，選取水稻建模樣點及精度評估樣點，建立水稻面積提取模型并對提取結(jié)果進行精度驗證。統(tǒng)計數(shù)據(jù)來自2019年和2020年江蘇省農(nóng)村統(tǒng)計年鑒，包括2019年和2020年江蘇省各地級市水稻種植面積。

江蘇省所在的MODIS影像序列號為h27v05和h28v05。文獻[18]研究證明覆蓋作物關(guān)鍵生育期的時間序列長度即可滿足一定分類精度，因此，本研究選擇的影像時期為2019年和2020年6月2日—11月9日(DOY153～DOY313)，包含水稻移栽期至成熟期。對得到的圖像進行波段運算得到EVI時間序列數(shù)據(jù)集。EVI計算公式為

(1)

式中ρnir——近紅外通道反射率

ρred——紅色通道反射率

ρblue——藍色通道反射率

L——土壤調(diào)節(jié)參數(shù)，取1

2 研究方法

2.1 水稻物候提取

基于遙感手段監(jiān)測植被物候涉及兩個重要步驟，一是植被指數(shù)的重構(gòu)，二是植被物候參數(shù)提取方法的選擇。目前，圖像處理已進入了人工智能時代，機器學(xué)習(xí)算法已成為時間序列數(shù)據(jù)處理的有力工具[19-21]。機器學(xué)習(xí)回歸算法(Machine learning regression algorithm, MLRAs)的高度適應(yīng)性和通用性使其比傳統(tǒng)的擬合函數(shù)能更精確地估計物候趨勢[22]。高斯過程回歸(Gaussian process regression，GPR)[23]作為一種最優(yōu)機器學(xué)習(xí)回歸算法在成功重建植被指數(shù)和提取可靠的物候指標方面具有潛力。本研究應(yīng)用GPR算法對時間序列進行重構(gòu)，在此基礎(chǔ)上提取物候信息，提取方法為：生長季開始日期/生長季結(jié)束日期在曲線的左/右部分達到曲線上升/衰減部分的季節(jié)振幅的30%時被識別。物候指標包括：生長季振幅、最大值、最大值日期、生長季開始日期(Start of season, SOS)、生長季結(jié)束日期(End of season, EOS)、季節(jié)積分、生長季長度。

在生長發(fā)育期識別的各個指標中，季節(jié)積分和生長季振幅是兩個非常重要的指標，不僅顯示了水稻的生育特點，還能夠反映水稻的生長狀況[24]。經(jīng)過反復(fù)實驗，本研究選用這兩個指標監(jiān)測江蘇省水稻種植面積。

2.2 水稻種植面積提取

本研究采用閾值法開展水稻面積提取。閾值法是圖像分割中的經(jīng)典方法，它利用圖像中要提取的目標與背景在灰度上的差異，通過設(shè)置閾值來把像素級分成若干類，從而實現(xiàn)目標的提取。最優(yōu)閾值選擇方法通常有人工經(jīng)驗選擇法、直方圖法和最大類間方差法[25]。本研究閾值確定方法選用直方圖法，結(jié)合Landsat8選取的水稻樣點和統(tǒng)計面積數(shù)據(jù)，分別確定江蘇省13個地級市季節(jié)積分和生長季振幅的閾值。確定步驟為：

(1)以高分辨率遙感影像(Landsat8 OLI)為基礎(chǔ)，以目視解譯的方式在每個地級市選擇若干個水稻樣點(圖1)，并提取各個水稻樣點對應(yīng)的MODIS-EVI時間序列。

圖1 2019年水稻樣點分布圖Fig.1 Distribution map of rice sampling points in 2019

(2)應(yīng)用MODIS-EVI時間序列提取每個水稻樣點的物候信息，得到每個水稻樣點的季節(jié)積分和生長季振幅指標。

(3)以地級市作為分區(qū)，應(yīng)用每個區(qū)域的若干個樣點，做該區(qū)域水稻的季節(jié)積分頻數(shù)分布圖和生長季振幅頻數(shù)分布圖。

(4)根據(jù)頻數(shù)分布圖，初步確定季節(jié)積分和生長季振幅的閾值，并將該值用于提取該區(qū)域水稻種植面積。

(5)將步驟(4)提取的水稻種植面積與該區(qū)域水稻種植面積的統(tǒng)計數(shù)據(jù)對比分析，根據(jù)頻數(shù)分布重新調(diào)整兩個指標的閾值，使水稻種植面積的總體誤差控制在±10%以內(nèi)，并將該閾值作為該區(qū)域的水稻提取閾值。

(6)應(yīng)用2019年的數(shù)據(jù)，經(jīng)過步驟(1)～(5)，確定各個地級市水稻提取的閾值，應(yīng)用該閾值提取2020年的水稻種植面積，并對2020年的提取結(jié)果進行精度評價。

水稻面積提取精度評價以兩種方式進行：與統(tǒng)計面積數(shù)據(jù)進行對比分析；通過高分辨率影像選取水稻和非水稻樣點，建立混合矩陣對提取結(jié)果進行定量分析。

3 結(jié)果分析

3.1 水稻物候信息提取

本研究以2019年的數(shù)據(jù)為建模數(shù)據(jù)，2020年的數(shù)據(jù)為精度驗證數(shù)據(jù)，采用2019年中、高空間分辨率的Landsat8 OLI遙感影像數(shù)據(jù)進行水稻信息目視解譯，在江蘇省內(nèi)選取了432個水稻樣點(圖1)，應(yīng)用GPR方法將2019年水稻樣點時間序列(6月2日—11月9日)進行重構(gòu)，重構(gòu)后的時間序列時間間隔為1 d(圖2a)。圖2a顯示GPR模型可以很好地將時間序列中的異常值和缺失值進行重構(gòu)，使曲線更加平滑合理，在此基礎(chǔ)上提取水稻物候信息，物候提取參數(shù)生長季振幅設(shè)置為30%(圖2b)。

圖2 單點時間序列重構(gòu)及水稻生育期提取Fig.2 Single-point time series reconstruction and extraction of rice growth period

圖2b顯示了水稻物候提取參數(shù)，其中橫坐標為日期，第1天為時間序列開始日期，即6月2日，通過與江蘇省水稻物候數(shù)據(jù)對比，該模型提取的SOS和EOS分別對應(yīng)江蘇水稻的分蘗期和成熟期。按以上方法提取所有水稻樣點的物候信息，并分區(qū)域整理所有樣點的季節(jié)積分和生長季振幅。

在水稻樣點物候信息提取的基礎(chǔ)上，逐像素提取2019年和2020年水稻物候信息，得到2019年和2020年生長季振幅和季節(jié)積分分布(圖3)，從圖3中可以看出，生長季振幅和季節(jié)積分的分布具有很好的區(qū)域性，水體和城市所在地區(qū)的值都較低，能夠很好的與周圍地物分開，植被覆蓋區(qū)顯示也較清晰，有利于水稻物候信息提取。

圖3 江蘇省2019年和2020年生長季振幅和季節(jié)積分Fig.3 Amplitude and area total of Jiangsu Province in 2019 and 2020

3.2 水稻種植面積提取

分析2019年432個水稻樣點的生長季振幅和季節(jié)積分，通過直方圖和水稻面積的統(tǒng)計數(shù)據(jù)確定13個市兩個指標的閾值(圖4)，圖4顯示，生長季振幅和季節(jié)積分閾值的分布也具有區(qū)域相似性，如南京市、蘇州市、無錫市、鎮(zhèn)江市、常州市，兩個指標的閾值均相似，淮安市和宿遷市也具有相似性，連云港市、南通市和鹽城市具有相似性等。

圖4 江蘇省13個市的生長季振幅和季節(jié)積分閾值Fig.4 Amplitude and area total thresholds of 13 cities in Jiangsu Province

將閾值相似的區(qū)域合并(圖5)，圖5充分體現(xiàn)了江蘇省水稻種植的特點：東部受海洋氣候的影響(連云港市、鹽城市和南通市)，其水稻種植具有一定的相似性，西部受陸地氣候的影響，形成了多種稻作生態(tài)區(qū)，南北縱跨5個緯度，南北水稻種植方式和格局有較大差異。徐州市作為江蘇省最北端的地級市與蘇南5市的差異較大，蘇中水稻種植情況也具有一定差異性，徐州市、連云港市和鹽城市等蘇北地區(qū)也具有一定的相似性(圖4)。因此，圖4和圖5從側(cè)面反映了模型建模的合理性。

圖5 閾值合并后的分區(qū)Fig.5 Partition map after threshold merging

根據(jù)2019年數(shù)據(jù)確定13個市的生長季振幅和季節(jié)積分的閾值，分別提取2020年13個市的水稻種植面積及空間分布，合并得到江蘇省2020年水稻種植面積及空間分布(圖6)。從圖中可以看出，水稻的分布具有很好的區(qū)域性，能夠較好地將水體、城市和山區(qū)等剔除，符合江蘇省水稻種植的分布情況。

3.3 精度驗證

首先，應(yīng)用2020年江蘇省13個地級市水稻種植面積的統(tǒng)計數(shù)據(jù)與MODIS提取數(shù)據(jù)進行對比分析(表1)。表1顯示2020年提取的全省水稻總面積的相對誤差為-6.10%。在13個市的面積統(tǒng)計中，連云港、淮安、宿遷、徐州、鹽城、揚州、泰州、南通和蘇州9個市的相對誤差都在±10%以內(nèi)，南京、鎮(zhèn)江、常州、無錫4個市的相對誤差較大。經(jīng)分析，可能原因有：①在建模過程中，蘇南4市選取的樣點較少，閾值的確定不夠準確。②蘇南4市水稻種植面積較少，在絕對誤差相同的情況下，相對誤差較高。③蘇南地形復(fù)雜，水稻種植情況不夠穩(wěn)定，每年的種植情況變化較大(從統(tǒng)計面積上看，2020年蘇南4市的水稻種植面積較2019年減少了近10%)，且水稻種植的空間分布情況變化也較大，因此，應(yīng)用2019年的數(shù)據(jù)建立的模型，不適用于2020年的水稻面積提取。④在本研究中使用的MODIS影像空間分辨率為500 m，即一個像元內(nèi)存在像元“同物異譜”和“同譜異物”的現(xiàn)象，對面積較小、地塊破碎的地區(qū)提取精度不夠，從而影響水稻的提取精度。除蘇南4市外，江蘇省大部分區(qū)域水稻種植面積隨時間變化具有一定的穩(wěn)定性，可以用該模型進行水稻面積的提取。在后續(xù)研究過程中，還可以不斷地加入歷史數(shù)據(jù)或最新數(shù)據(jù)，不斷更新調(diào)整閾值，從而更準確地把握每個區(qū)域的水稻種植特點，使參數(shù)具有穩(wěn)定性、代表性和可推廣性。

表1 2019年和2020年江蘇省水稻面積提取結(jié)果與農(nóng)業(yè)統(tǒng)計數(shù)據(jù)對比Tab.1 Comparison of rice area extraction results and agricultural statistics in Jiangsu Province in 2019 and 2020

其次，僅僅依靠統(tǒng)計面積的比較進行精度驗證，在理論上存在一定程度的不可靠性，因為相同的面積可能會有不同的形式組合，用高分辨率的遙感影像驗證大尺度(低分辨率)提取精度，這種方法的優(yōu)點在于省時省力效率高，彌補缺乏野外驗證的不足。本研究利用Landsat8和目視解譯方法選取各個地類的樣點共3 423個，通過混淆矩陣對提取結(jié)果進行定量分析(表2)，分析結(jié)果表明，水稻的制圖精度為92.90%，用戶精度為89.09%。計算得水稻提取的總體精度為92.55%，Kappa系數(shù)為0.846 3。

表2 混淆矩陣精度評價Tab.2 Confusion matrix accuracy evaluation

將非水稻樣點分為水體、建筑、其他植被及其他4類，分析各個類別的提取精度(表3)。其中水體、建筑及其他類型由于與水稻差異較大，精度分別為99.52%、98.11%和92.33%，而其他植被由于與水稻具有相似性，因此錯分為水稻的概率較大，其精度為87.73%。

表3 非水稻類型精度評價Tab.3 Accuracy evaluation for non-rice types

在精度評價方面，盡管已經(jīng)開展了大量工作，但是對這種光譜指數(shù)時間序列水稻面積提取方法的精度評價還不夠充分，本研究主要是針對統(tǒng)計年鑒統(tǒng)計數(shù)據(jù)進行的行政區(qū)尺度水稻面積評估，雖然本研究也應(yīng)用高空間分辨率影像，提取樣點進行了精度驗證，但仍不能全面完整地體現(xiàn)面積提取方法的空間分布精度，對于水稻面積總量漏分和誤分的空間分布的精度評價工作還有待進一步深入研究。

4 結(jié)論

(1)基于機器學(xué)習(xí)的GPR模型能夠較好地重構(gòu)MODIS-EVI時間序列數(shù)據(jù)集的異常點和缺失值，為后續(xù)水稻物候期的提取提供了較好的基礎(chǔ)。根據(jù)提取的水稻物候信息，提取水稻種植面積，精度較好，與統(tǒng)計數(shù)據(jù)相比，大部分地級市的水稻提取精度都在90%以上。與高分辨率影像數(shù)據(jù)相比，水稻提取的總體精度為92.55%，Kappa系數(shù)為0.846 3，水稻的制圖精度為92.90%，水稻的用戶精度為89.09%。由于其他植被與水稻生長具有一定的相似性，其他植被類型的錯分率較高，其精度為87.73%，水體、建筑和其他類型的錯分率都較低，精度都在90%以上。

(2)江蘇省水稻種植方式、種植格局有較大差異，本研究分13個地級市確定水稻物候指標的閾值，通過閾值相似性合并，很好地體現(xiàn)了江蘇省水稻種植的特點。該結(jié)論一方面說明了本研究的合理性，另一方面說明對大區(qū)域水稻種植面積的分析與建模，需要分區(qū)域進行。

(3)本研究僅應(yīng)用2019年的數(shù)據(jù)作為建模數(shù)據(jù)，將模型用于2020年水稻面積提取，并取得了較好的結(jié)果。