基于生長氣候分區(qū)的土地利用農(nóng)作物種植面積快速提取研究
————以南京市油菜種植面積提取為例

肖善才

（南京農(nóng)業(yè)大學(xué)土地管理學(xué)院，江蘇·南京 210095）

中國是一個農(nóng)業(yè)生產(chǎn)大國，農(nóng)業(yè)一直是其基礎(chǔ)產(chǎn)業(yè)，及時掌握農(nóng)作物空間分布信息及其種植規(guī)模對監(jiān)測土地利用變化、預(yù)測農(nóng)作物產(chǎn)量及確保國家糧食安全具有重要意義[1,2]。傳統(tǒng)的逐層統(tǒng)計(jì)匯報(bào)和實(shí)地測量方法誤差大、耗時耗力，且缺乏直觀的空間分布信息。而遙感技術(shù)在土地利用、地表覆蓋調(diào)查等方面具有獨(dú)特的優(yōu)勢，近年來在農(nóng)業(yè)生產(chǎn)、水資源管理方面的應(yīng)用越來越普及。

國內(nèi)外關(guān)于作物種植面積監(jiān)測的研究較多[3-7]，王連喜等利用中等分辨率的Landsat8影像數(shù)據(jù)和低等分辨率的MDDIS-NDVI時序數(shù)據(jù)采用決策樹和混合像元分解相結(jié)合的方法提取江蘇省的冬小麥種植面積，精度達(dá)到90%[8]；張健康等利用多時相的TM/ETM+遙感影像數(shù)據(jù)和13幅時間序列的MODIS-EVI遙感影像數(shù)據(jù)，采用基于生態(tài)分類法的監(jiān)督分類與決策樹分類相結(jié)合的人機(jī)交互解譯方法，對黑龍港地區(qū)的主要作物進(jìn)行了解譯，精度達(dá)到了91.3%[9]；張煥雪等利用多時相環(huán)境衛(wèi)星HJ-1A/BCCD的NDVI時序數(shù)據(jù)，采用決策樹算法對黑龍江省北安市紅星農(nóng)場開展了農(nóng)作物分類研究，通過與單時相遙感數(shù)據(jù)分類結(jié)果對比分析，發(fā)現(xiàn)NDVI時序數(shù)據(jù)比僅使用3個時相的HJ-1 A/BCCD數(shù)據(jù)分類精度提高了5.45%[10]；Jakubauskas等研究了利用低空間分辨率NOAA/AVHRR數(shù)據(jù)作為數(shù)據(jù)源，并對耕地和作物空間格局動態(tài)變化進(jìn)行了大尺度范圍研究[11]，這些研究都取得了良好的效果并豐富了遙感技術(shù)在農(nóng)業(yè)生產(chǎn)方面應(yīng)用的深度與廣度。

然而目前多數(shù)研究都是在中大尺度上提取作物種植信息[12-21]，小尺度的作物種植信息提取研究鮮有關(guān)注。并且這些研究大都忽略了即使是在同一時期內(nèi)，研究區(qū)內(nèi)部各處的氣候也是有差異的，這會導(dǎo)致同一種作物生育期的變化不一致，因而直接采用NDVI時序數(shù)據(jù)提取研究區(qū)的作物種植信息可能會對提取結(jié)果的精度產(chǎn)生影響。王凱、王連喜等也指出氣候差異會導(dǎo)致作物的生育期不一致，從而對作物種植分布信息提取產(chǎn)生較大的誤差[2,8]。為解決上述問題，本文提出了一種高精度小尺度的多時相作物種植信息快速提取方法。首先根據(jù)影響農(nóng)作物生長發(fā)育的關(guān)鍵氣候條件，結(jié)合研究區(qū)的背景數(shù)據(jù)，劃分出合理的氣候區(qū)，然后采用多時相的中等分辨率遙感影像建立起各氣候區(qū)的作物種植信息提取模型，并利用高分辨率的遙感影像對提取結(jié)果進(jìn)行進(jìn)一步修正，得到最終的作物種植信息，最后利用南京市2016年油菜種植面積進(jìn)行提取驗(yàn)證。

1 研究方法

1.1 數(shù)據(jù)來源與預(yù)處理

遙感數(shù)據(jù)資料為美國地質(zhì)調(diào)查局（United States Geological Survey，USGS）提供的 2015–2016年 Landsat8陸地成像儀（operational land imager，OLI）影像以及2016年的資源3號影像數(shù)據(jù)。Landsat8攜帶的陸地成像儀（OLI）包括9個波段，其中有8個分辨率為30m的多光譜波段：b1（Coastal，0.43–0.45um）、b2（Blue，0.45–0.51um）、b3（Green，0.53–0.59um）、b4（Red，0.64–0.67um）、b5（NIR，0.85–0.88um）、b6（SWIRI，1.57–1.65um）、b7（SWIRI2，2.11–2.29um）、b9（Cirrus，1.36–1.38um），1個分辨率為15m的全色波段：b8（Pan，0.50–0.68um），成像寬幅為185×185km。其中近紅外波段b5排除了0.825um處水汽吸收特征；全色波段b8范圍收窄，這種方式可以在全色圖像上更好地區(qū)分植被和無植被特征。只需要一景影像就可以覆蓋整個南京市，獲取6景的Landsat8影像的日期分別為2015年10月15日、2016年2月25日、2016年3月12日、2016年3月28日、2016年4月29日、2016年5月19日，通過計(jì)算機(jī)獲取它們的NDVI，并進(jìn)行編號，依次為A、B、C、D、E、F，如圖1所示。

圖1 南京市6景Landsat8-NDVI影像Fig.1 Landsat8-NDVI images of Nanjing at six dates

資源3號衛(wèi)星（ZY-3）是我國2012年發(fā)射成功的第一顆民用高分辨率立體測試衛(wèi)星，至今已在生產(chǎn)中得到了廣泛應(yīng)用。資源三號影像共有5個波段，其中4個為多光譜波段：b1（Blue，0.45–0.52um）、b2（Green，0.52–0.59um），b3（Red，0.63–0.69um），b4（NIR，0.77–0.89um），一個全色光波段：b5（Pan：0.45–0.80um），星下點(diǎn)全色分辨率為2.1m，前、后視22°全色分辨率3.6m，星下點(diǎn)多光譜分辨率為5.8m，星下點(diǎn)全色成像寬幅為50×50km，星下點(diǎn)多光譜成像寬幅為52×52km。由于資源3號影像較難獲取，因此本文只獲取了單時相的2景資源3號影像，日期為2016年5月19日，且未能覆蓋全南京。

本文對遙感數(shù)據(jù)的預(yù)處理主要包括：（1）對所有影像進(jìn)行輻射校正、大氣校正、幾何校正和重采樣；（2）利用Mosaicking工具將2景資源三號數(shù)據(jù)鑲嵌，然后利用行政邊界矢量數(shù)據(jù)對該幅影像和landsat8影像進(jìn)行裁剪；（3）對得到的影像數(shù)據(jù)疊加處理，得到每幅圖像相對應(yīng)的NDVI值，便于分析作物參數(shù)，確定分類閾值。

1.2 氣候分區(qū)的必要性與劃分步驟

氣候條件作為影響農(nóng)業(yè)生產(chǎn)的重要自然條件之一，直接或間接地為農(nóng)業(yè)生物的生長發(fā)育及其產(chǎn)量形成提供必需的物質(zhì)與能量。橘生淮南則為橘，生于淮北則為枳，這句話表明了由于氣候條件的不同，具有一定南北跨度的地域內(nèi)不同地區(qū)的同類作物生長情況必然會有差異。由于氣候條件的差異造成的某種作物生育期的變化則會導(dǎo)致不同地區(qū)的該種作物的NDVI的變化趨勢產(chǎn)生差異，影響作物種植信息提取模型的建立從而影響到提取精度，因此根據(jù)研究區(qū)內(nèi)各處的氣候差異劃分氣候區(qū)十分必要。影響作物生長發(fā)育的主要?dú)夂驐l件有溫度、水分、光照以及二氧化碳和風(fēng)等，其中溫度是影響作物生長發(fā)育期的重要因素，作物只有在達(dá)到一定的積溫后才能發(fā)育成熟，本文重點(diǎn)研究了溫度對作物生長的影響。其主要劃分步驟簡述如下：首先獲取研究區(qū)的氣候數(shù)據(jù)，數(shù)據(jù)來源于中國氣候數(shù)據(jù)共享網(wǎng)的《中國地面累年值日值數(shù)據(jù)集》，得到各氣候站點(diǎn)的溫度數(shù)據(jù)；然后采用克里金方法進(jìn)行插值得到溫度分布圖，結(jié)合研究區(qū)背景數(shù)據(jù)，確立分區(qū)標(biāo)準(zhǔn)，劃分氣候區(qū)；最后以野外調(diào)查數(shù)據(jù)對分區(qū)結(jié)果進(jìn)行驗(yàn)證。

1.3 作物種植信息提取模型

在氣候分區(qū)的基礎(chǔ)上，結(jié)合作物的物候期，分析各區(qū)不同農(nóng)作物的NDVI變化趨勢，選擇出合適時相的遙感影像，提取出區(qū)分不同類型作物的關(guān)鍵值作為閾值，采用人機(jī)交互式的決策樹方法建立起各氣候區(qū)作物種植信息的提取模型從而提取出各區(qū)的作物種植信息，將各氣候區(qū)的結(jié)果相加則可得到整個研究區(qū)的作物種植情況。為了獲取小尺度的研究區(qū)作物種植信息，中等分辨率的遙感影像已經(jīng)不能滿足要求，因此利用高分辨率的影像對提取結(jié)果進(jìn)行校正，修正低分辨率地區(qū)，使其滿足高分辨率，從而獲得小尺度的提取結(jié)果。具體的做法是在中等分辨率和高等分辨率的遙感影像上裁剪出多塊大面積種植作物的研究區(qū)，以相同的方法提取出兩種影像的作物種植信息，對比獲得兩種影像的比例系數(shù)，并以此修正中等分辨率的研究區(qū)影像進(jìn)而得到小尺度的作物種植信息。

1.4 提取結(jié)果的驗(yàn)證方法

為了驗(yàn)證該方法的合理性，采用空間和定量兩種方法對提取結(jié)果進(jìn)行驗(yàn)證?？臻g驗(yàn)證方法利用作物感興趣區(qū)與上述方法提取出的作物種植分布信息進(jìn)行對比，確定提取精度；定量驗(yàn)證方法采用統(tǒng)計(jì)數(shù)據(jù)與提取的作物種植面積進(jìn)行比較，從而確定提取精度。

2 實(shí)例應(yīng)用研究

2.1 研究區(qū)概況

2.2 氣候分區(qū)

2.2.1 氣候條件的選取

影響油菜生長發(fā)育的主要因素是溫度和水分，根據(jù)中國地面累年值年值數(shù)據(jù)集查得南京市不同地區(qū)的年降雨量并無明顯差異，但積溫差異較大，因此選取積溫為核心氣候指標(biāo)，劃分氣候區(qū)域，南京市氣候站點(diǎn)油菜各生育期平均溫度見表1，從北到南依次為六合、浦口、南京主城區(qū)、溧水、高淳。由表1可以看出，由北向南各生育期的平均溫度逐漸升高，積溫也呈現(xiàn)相似的趨勢，且六合與高淳的積溫相差最大。

表1 南京市氣候站點(diǎn)油菜各生育期平均溫度（℃）Table 1 Average temperature of rape growth period in Nanjing meteorological sites (℃)

2.2.2 插值處理

克里金法是1951年由南非地質(zhì)學(xué)家 Krige 提出，并因此而命名。在1962年法國地理學(xué)家Matheron給出了Kriging的一般公式，其公式為：

式中，s為不同位置的點(diǎn)，可以認(rèn)為是用經(jīng)緯度表示的空間坐標(biāo)；z(s)為s處的變量值，其可分解確定趨勢值μ(s)和自相關(guān)隨機(jī)誤差?(s)。通過對這個公式進(jìn)行變化，可以生成克里金方法的不同類型。

在地統(tǒng)計(jì)分析中，克里金法是建立在平穩(wěn)假設(shè)的基礎(chǔ)上，這種假設(shè)在一定程度上要求所有數(shù)據(jù)集具有相同的變異性。另外，一些克里金插值（如普通克里金法、簡單克里金法和泛克里金法等）都假設(shè)數(shù)據(jù)服從正態(tài)分布。如果數(shù)據(jù)不服從正態(tài)分布，需要進(jìn)行一定的數(shù)據(jù)變換，使其服從正態(tài)分布。因此，在地統(tǒng)計(jì)分析前，需要對選取的氣候指標(biāo)數(shù)據(jù)檢驗(yàn)其分布特征，剔除異常值后，進(jìn)行插值得到氣候區(qū)。首先利用正態(tài)QQ圖檢驗(yàn)17個站點(diǎn)的積溫分布特征如圖2所示。

圖2 氣候站點(diǎn)積溫正態(tài)QQ圖Fig.2 Accumulated temperature normal QQ figure of meteorological sites

從圖2中可以發(fā)現(xiàn)，除右上角的一個站點(diǎn)—高淳的數(shù)據(jù)外基本都在一條直線上，即這些站點(diǎn)都服從正態(tài)分布，采用普通克里金法對南京市周圍的16個氣候站點(diǎn)的積溫?cái)?shù)據(jù)進(jìn)行插值。普通克里金法是區(qū)域化變量的線性估計(jì)，它假設(shè)數(shù)據(jù)變化成正態(tài)分布，認(rèn)為區(qū)域化變量Z的期望值是未知的常量。插值的過程類似于加權(quán)滑動平均，權(quán)重值的確定來自于空間數(shù)據(jù)分析，插值后得到南京市氣候區(qū)分布圖，由北向南劃分為北區(qū)、中區(qū)和南區(qū)。劃分的依據(jù)即積溫閾值是根據(jù)南京的農(nóng)作物種植格局來確定的，南京市南部如高淳、溧水，北部如六合是傳統(tǒng)的農(nóng)業(yè)種植中心，區(qū)內(nèi)油菜較其他行政區(qū)都呈大面積分布，而在中部如浦口、江寧等主要是城市和鄉(xiāng)鎮(zhèn)，油菜呈小面積或零星分布。依據(jù)該種情況劃分了氣候區(qū)，如圖3所示。

圖3 南京市氣候區(qū)分布圖Fig.3 Distribution of climate area of Nanjing

從圖3中可以看出，南京市的積溫由南向北逐漸減少，這說明江南與江北地區(qū)的油菜成熟時間有差異，且由南向北，油菜成熟時間越來越晚。這符合野外實(shí)地調(diào)研過程中所觀察到的現(xiàn)象——在5月20日至5月24日為期4天的野外調(diào)查中，我們發(fā)現(xiàn)江南的油菜例如高淳和溧水地區(qū)油菜成熟的時間較早，某些地區(qū)的油菜甚至已經(jīng)完成了收割工作，而江北地區(qū)的油菜例如六合和浦口地區(qū)等普遍沒有成熟，距離收割工作還需一段時間。氣候區(qū)的劃分與野外調(diào)研的結(jié)果基本一致，南區(qū)積溫最高，油菜成熟最早，基本已完成收割；中區(qū)積溫居中，油菜成熟比南區(qū)晚，收割工作基本還未開始；北區(qū)積溫最低，油菜距離成熟還需一段時間，因此，溫度是影響南京市油菜成熟期的重要因素，以積溫為氣候指標(biāo)為南京市劃分氣候區(qū)合理。

2.3 油菜種植信息提取模型建立、結(jié)果提取與校正

2.3.1 油菜提取模型的建立

南京市越冬農(nóng)作物只有冬小麥和油菜，而油菜和冬小麥具有相同的物候期，都是在10月下旬至次年5月上旬；所以，提取油菜種植面積的關(guān)鍵是區(qū)分油菜和冬小麥。表2為南京市油菜與冬小麥的物候期。

表2 南京市油菜與冬小麥物候期Table 2 Rape and winter wheat phenological period in Nanjing

油菜和冬小麥在11–12月處于緩慢生長階段，NDVI曲線平緩上升；1–2月處于越冬期，油菜和冬小麥停止生長，NDVI變化很??；3–4月返青后，油菜和冬小麥均快速生長，NDVI曲線快速上升；4–5月成熟收獲期，NDVI曲線開始下降。兩者相比較，在3月中旬之后，油菜與冬小麥的NDVI曲線有相反的變化，在此期間油菜進(jìn)入盛花期，所以NDVI曲線有所下降；而冬小麥在拔節(jié)期繼續(xù)生長，其NDVI曲線則繼續(xù)上升，因此可根據(jù)這種差異區(qū)分出油菜和小麥。而且同時間段內(nèi)的其他農(nóng)作物（如早稻、棉花及玉米）尚未移栽，所以3、4月份是識別冬小麥和油菜的最佳時相。然而由于氣候條件、種植時間等因素的不同，3、4月份的遙感影像中油菜和冬小麥的NDVI值會出現(xiàn)相同的情況，單從一景影像很難區(qū)分油菜與冬小麥，因此利用多時相遙感影像進(jìn)行作物種植信息提取。

Landsat-8影像的選擇也遵循上述規(guī)律，2015年10月15日油菜和小麥還未播種，此時NDVI值很低；2016年的2月25日至3月12日是油菜和小麥返青時期，NDVI值迅速上升；2016年3月28日至4月29日小麥在拔節(jié)期繼續(xù)生長，NDVI繼續(xù)上升，而油菜進(jìn)入盛花期，NDVI降低；2016年5月19日油菜已經(jīng)基本成熟或已收割完成，NDVI值很低。

利用野外調(diào)查采取的油菜和小麥樣點(diǎn)結(jié)合油菜與小麥的物候期及它們的NDVI值變化趨勢，采用決策樹方法建立起各個氣候區(qū)的油菜面積提取模型。各個氣候區(qū)的部分油菜樣點(diǎn)的6個時相的NDVI值如表3所示。

表3 氣候區(qū)油菜樣點(diǎn)6個時相的NDVI值Table 3 NDVI value of rape sample in climate area at 6 dates

決策樹各分支約束條件的閾值即野外樣點(diǎn)對應(yīng)的NDVI值，通常以某一類別的最小值為閾值，爭取到大部分相似的像元，然后再輔以其他約束條件進(jìn)行篩選和刪除。根據(jù)NDVI曲線、NDVI值和決策樹模型閾值選取原則，得到每個子研究區(qū)的決策樹識別模型如下（式中b1、b2、b3、b4、b5、b6 代表各個時相的遙感影像的NDVI值）：

（1）北區(qū)油菜識別模型：對北區(qū)的油菜NDVI值的變化趨勢研究后發(fā)現(xiàn)，小麥的NDVI值在第4個時相和第5個時相時總有一個大于0.68，油菜的NDVI值最高也沒有超過0.68。那么，北區(qū)油菜識別模型為：

b1＜b2 and b2＜b3 and b3＜b4 and b5＞b6

b1＜0.2 and b6＜0.2 and b4＜0.68 and b5＜0.68

（2）中區(qū)油菜識別模型：研究中區(qū)的油菜NDVI值的變化趨勢后可以知道，研究區(qū)內(nèi)各個樣點(diǎn)的油菜NDVI值的變化趨勢基本相同，第3時相到第4時相油菜的NDVI值大于0.5，第4時相的NDVI值小于0.65，而小麥大都大于0.65，第5時相油菜的NDVI值都大于0.58，則中區(qū)油菜識別模型為：

b1＜b2 and b2＜b3 and b5＞b6

b1＜0.2 and b6＜0.2 and b3＞0.5 and b4＜0.65 and b5＞0.58

（3）南區(qū)油菜識別模型：同理，研究南區(qū)油菜的NDVI值的變化趨勢后發(fā)現(xiàn)，第3時相到第4時相NDVI值的降低是區(qū)分油菜和小麥的最明顯標(biāo)志，此時高淳區(qū)內(nèi)油菜的NDVI值處于上升階段，且小麥的NDVI值在第四個時相時總是大于0.65。因此，南區(qū)油菜識別模型為：

b1＜b2 and b2＜b3 and b3＞b4

b1＜0.2 and b2＞0.35 and b4＜0.65

2.3.2 提取結(jié)果

對提取結(jié)果研究分析后發(fā)現(xiàn)，多時相的影像很好地區(qū)分了油菜與小麥，其中C、D、E三個時相是區(qū)分油菜與小麥的關(guān)鍵。由表3可以看出，區(qū)分效果最好的是南區(qū)，只利用了4個時相的影像就很好的區(qū)分開油菜與小麥，其中C時相的影像是區(qū)分油菜與小麥的關(guān)鍵，C時相之后，油菜的NDVI值呈下降趨勢，而小麥的NDVI值呈上升趨勢；中區(qū)和南區(qū)的區(qū)分效果較南區(qū)較差，都需要用6個時相的影像才能較好地區(qū)分油菜與小麥，其中D、E兩個時相的影像是區(qū)分油菜與小麥的關(guān)鍵，中區(qū)和北區(qū)的油菜NDVI值變化趨勢各異，因此只能設(shè)置合理的閾值來區(qū)分油菜與小麥。這可能是由于中區(qū)和北區(qū)油菜種植較為分散，部分油菜通常是伴隨著大面積的小麥種植的，因而各處的NDVI值變化趨勢各異，而南區(qū)油菜是呈大面積種植，小麥種植較油菜相對較少，所以區(qū)分效果較中區(qū)和北區(qū)更佳。

利用Landsat8數(shù)據(jù)采用上述方法提取的南京市油菜種植面積為17.28千公頃。

2.3.3 高分辨率遙感影像校正結(jié)果

以資源3號衛(wèi)星數(shù)據(jù)來校正低空間分辨率遙感影像像元，利用南京市行政區(qū)劃矢量圖裁剪后，得到部分研究區(qū)的影像，選取該研究區(qū)中國油菜呈大面積種植的區(qū)域，采用決策樹分類方法對這些區(qū)域內(nèi)的油菜面積進(jìn)行提取，得到5m小尺度油菜面積。裁剪出與這些區(qū)域相同大小的中等精度的Landsat8影像數(shù)據(jù)，運(yùn)用相同的決策樹方法得到該地區(qū)的中等空間分辨率的油菜面積，從而得到該地區(qū)的校正系數(shù)，利用該系數(shù)對基于氣候區(qū)建立的油菜面積提取模型提取出的結(jié)果進(jìn)行校正，得到5m小尺度的南京市油菜種植信息。

X為校正系數(shù)，S1為高分辨率影像的作物面積，S2為低分辨率影像的作物面積。

在利用Landsat8數(shù)據(jù)提取的油菜面積結(jié)果上，利用資源3號進(jìn)行修正后的油菜種植面積為17.80千公頃。

2.3.4 提取精度驗(yàn)證

（1）定量精度分析：根據(jù)南京市農(nóng)委的統(tǒng)計(jì)資料顯示2016年南京市油菜種植面積約為20千公頃，則利用高分辨率的資源3號影像修正前和修正后的提取精度分別為86.3%和89%，這表明采用上述方法能夠提取出高精度小尺度的油菜。

（2）空間精度分析：利用野外調(diào)查選取的樣點(diǎn)制作的油菜感興趣區(qū)與分類結(jié)果進(jìn)行對比得到油菜的提取精度，如表4所示。

表4 基于氣候區(qū)的油菜提取精度Table 4 Rape extraction precision based on climate area

從表4可以看出，南區(qū)的油菜提取精度達(dá)到了90%以上，這與實(shí)際情況相符，較北區(qū)和中區(qū)，南區(qū)的油菜呈大面積分布，特別是在固城湖和石臼湖附近分布很集中，純凈像元較多，所以提取精度能夠達(dá)到90%；中區(qū)的油菜提取精度為為85.10%，是3個氣候區(qū)中提取精度最低的，這是因?yàn)橹袇^(qū)內(nèi)低山，丘陵，崗地較南北區(qū)較多，境內(nèi)有老山、黃龍山、湯山等，故油菜分布較為分散，因而提取精度較低；北區(qū)的油菜提取精度為87.84%，研究后發(fā)現(xiàn)，北區(qū)北部為丘陵山崗地區(qū)，油菜分布分散，而中南部為河谷平原、崗地區(qū)，油菜分布較為集中，因而提取精度只有87.84%。采用上述方法提取的各區(qū)的油菜面積精度都在85%以上，且總體精度達(dá)到了89.40%，說明了本方法能夠很好的提取油菜信息。

2.4 方法創(chuàng)新性驗(yàn)證

為了能夠更好地體現(xiàn)本方法的創(chuàng)新性，本文在不劃分氣候區(qū)的基礎(chǔ)上，以南京市為研究區(qū)域利用多時相遙感影像，提取南京市油菜種植面積，進(jìn)而與文章提取的新方法比較優(yōu)勢性。

通過對各氣候區(qū)的油菜NDVI值變化趨勢的研究后，綜合每個區(qū)的NDVI值的變化趨勢，確定整個南京地區(qū)的油菜識別模型：

b2＜b3 and b1＜b2 and b5＞b6

b2＞0.35 and b3＞0.45 and b4＞0.48 and b5＞0.45

b1＜0.2 and b3＜0.65

b4＜0.7 and b5＜0.7and b6＜0.2

（1）定量精度分析：以校正系數(shù)來修正上述方法提取出的結(jié)果，得到南京市油菜種植面積為24.20千公頃，較南京市真實(shí)油菜播種面積多出了4.20千公頃，比文章提取的新方法提取出的油菜面積多了6.40千公頃。

（2）空間精度分析：利用ENVI軟件結(jié)合野外調(diào)研資料制作的油菜驗(yàn)證數(shù)據(jù)與分類結(jié)果進(jìn)行對比得到油菜的提取精度，其對照像元4147個，正確分類像元3579個，分類精度為86.3%。由此可知，以南京市為研究區(qū)提取油菜面積精度為86.3%，低于表4中以氣候區(qū)提取油菜面積精度，后者比前者高了3.1%。

由定量和空間兩種分析方法可知，本文所提出的方法較已有的不考慮研究區(qū)內(nèi)部氣候差異，直接以研究區(qū)作為整體建立作物種植信息提取模型來提取作物分布情況的方法具有創(chuàng)新性。

3 結(jié)論與討論

獲取高精度的作物種植規(guī)模一直都是作物種植信息提取研究中的重點(diǎn)，本文根據(jù)研究區(qū)內(nèi)部氣候差異，劃分氣候區(qū)，然后利用中等分辨率的多時相遙感影像建立作物種植面積提取模型，并利用高分辨率的影像對提取結(jié)果修正，最后通過南京市2016年油菜種植面積進(jìn)行提取檢驗(yàn)。與統(tǒng)計(jì)數(shù)據(jù)相比，總體精度達(dá)到89%，表明該方法能夠很好地提取出作物種植信息。為了更好地表明本方法的優(yōu)勢，將本研究提出的新方法與當(dāng)前多數(shù)研究以研究區(qū)為整體提取作物分布信息進(jìn)行比較[22-24]，結(jié)果表明運(yùn)用本文方法提取的結(jié)果優(yōu)于后者，證明本文提出的新方法更具優(yōu)勢性。由于具有一定南北跨度的地區(qū)之間氣候必然會有差異，而這樣的地區(qū)在中國很多見，因此本方法又具有普遍適用性。

另外，與以往較為復(fù)雜的混合像元分解法[8,12]相比——以最小二乘法為基礎(chǔ)，建立分解模型，本文中提出了以高等空間分辨的影像像元來校正中等空間分辨率的像元的方法，并在實(shí)例研究中取得了較好的效果，提取精度提高了2.7%。在未使用多景高分辨率遙感影像的情況下，只利用一景中分辨率的影像就獲得了高分辨率高精度的南京市作物種植分布信息，證明該種方法切實(shí)可行。

Wardlow B.D.等利用12個月16d一期的MODIS數(shù)據(jù)對美國中部大平原的作物信息進(jìn)行了提取[15]，Pringle等人也利用了MODIS數(shù)據(jù)對土地覆蓋進(jìn)行分類并對特定地物進(jìn)行了面積的提取[24]。雖然MODIS數(shù)據(jù)時間分辨率高，但空間分辨率低，多針對于大尺度區(qū)域的作物種植結(jié)構(gòu)提取，面向小尺度研究區(qū)域則容易出現(xiàn)較大誤差[25]，且所需要的景數(shù)較多，而如果都利用高空間分辨率的影像來提取小尺度的作物信息，一是成本太高，可能需要多景影像才能覆蓋一個研究區(qū)，二是合適高空間分辨率影像難獲取，例如影像云量較多又或者作物關(guān)鍵生育期時相的影像沒有等。本文的方法中只需要少量的可免費(fèi)獲取中等空間分辨率遙感影像外加幾景高空間分辨率的遙感影像結(jié)合作物物候期就能夠提取出小尺度的作物種植信息，具有經(jīng)濟(jì)高效性。

但是，本文提出的方法也有一些需要改進(jìn)的地方：在建立氣候區(qū)時，只考慮了溫度因素，而沒考慮其他對作物有影響的氣候條件，它們對氣候區(qū)的劃分也可能會有一定的影響；沒有考慮農(nóng)民播種時間的影響，這可能會對提取結(jié)果有一定的影響。

參考文獻(xiàn)(References)

[1]王松林,張佳華,劉學(xué)鋒. 基于MODIS多時相的江蘇啟東市油菜種植面積提取[J]. 遙感技術(shù)與應(yīng)用,2015,30(5):946-951.Wang S L, Zhang J H, Liu X F. Extraction of rape planting areas based on multi-temporal MODIS data in Qidong, Jiangsu province[J].Remote Sensing Technology and Application,2015,30(5):946-951.

[2]王凱,張佳華. 基于MODIS數(shù)據(jù)的湖北省油菜種植分布信息提取[J]. 國土資源遙感,2015,27(3):65-70.Wang K, Zhang J H. Extraction of rape seed cropping distribution information in Hubei province based on MODIS images[J].Remote Sensing for Land and Resources, 2015,27(3):65-70.

[3]陳思寧,趙艷霞,申雙和. 基于波譜分析技術(shù)的遙感作物分類方法[J]. 農(nóng)業(yè)工程學(xué)報(bào),2012,28(5):154-160.Chen S N, Zhao Y X, Shen S H. Crop classification by remote sensing based on spectral analysis[J].Transactions of the Chinese Society of Agricultural Engineering, 2012,28(5):154-160.

[4]李丹丹. 基于中高空間分辨率遙感影像的油菜種植面積信息提取研究[D]. 北京:中國農(nóng)業(yè)科學(xué)院,2011.Li D D. Crop planted area extraction based on middle and high resolution remote sensing images[D]. Beijing: Chinese Academy of Agricultural Sciences, 2011.

[5]楊嘯. 基于時序NDVI的湖北省植被覆蓋動態(tài)變化監(jiān)測分析[J].長江流域資源與環(huán)境,2013,22(2):226-231.Yang X. Vegetation cover dynamic monitoring of Hubei province based on time-series NDVI[J].Resource and Environment in the Yangtze Basin, 2013,22(2):226-231.

[6]Liu T T, Liu C. Study on extraction of crop information using time-series MODIS data in the Chao Phraya Basin of Thailand[J].Advances in Space Research, 2010,45(6):775-784.

[7]李芬. 資源三號衛(wèi)星數(shù)據(jù)在土地利用遙感監(jiān)測中的應(yīng)用研究[D]. 長春:吉林大學(xué),2013.Li F. Extraction of vegetation information based on resources satellite No.3 (ZY-3) data[D]. Changchun: Jilin University, 2013.

[8]王連喜,徐勝男,李琪. 基于決策樹和混合像元分解的江蘇省冬小麥種植面積提取[J]. 農(nóng)業(yè)工程學(xué)報(bào),2016,32(5):182-187.Wang L X, Xu S N, Li Q. Extraction of winter wheat planted area in Jiangsu province using decision tree and mixed-pixel methods[J].Transactions of the Chinese Society of Agricultural Engineering,2016,32(5):182-187.

[9]張健康,程彥培,張發(fā)旺. 基于多時相遙感影像的作物種植信息提取[J]. 農(nóng)業(yè)工程學(xué)報(bào),2012,28(2):134-141.Zhang J K, Cheng Y P, Zhang F W. Crops planting information extraction based on multi-temporal remote sensing images[J].Transactions of the Chinese Society of Agricultural Engineering,2012,28(2):134-141.

[10]張煥雪,曹新,李強(qiáng)子. 基于多時相環(huán)境星NDVI時間序列的農(nóng)作物分類研究[J]. 遙感技術(shù)與應(yīng)用,2015,30(2):304-311.Zhang H X, Cao X, Li Q Z. Research on crop identification using multi-temporal NDVI HJ images[J].Remote Sensing Technology and Application, 2015,30(2):304-311.

[11]Jakubauskas M E, Legates D R, Kastens J H. Crop identification using harmonic analysis of time-series AVHRR NDVI data[J].Computer and Electronics in Agriculture, 2002,37(1):127-139.

[12]Lunetta R S, Shao Y, Ediriwickrema J, et al. Monitoring agriculture cropping patterns across the Laurentian Great Lakes Basin using MODIS-NDVI data[J].International Journal of Applied Earth Observation and Geoinformation, 2010,(12):81-88.

[13]Wardlow B D, Egbert S L, Kastens J H. Analysis of time-series MODIS 250m vegetation index for crop classification in the US Central Great Plains[J].Remote sensing of Environment,2007(108):290-310.

[14]梁益同,萬君. 基于HJ-1A/B-CCD影像的湖北省冬小麥和油菜分布信息的提取方法[J].中國農(nóng)業(yè)氣象,2012,33(4):573-578.Liang Y T, Wan J. Application of HJ-1A/B-CCD images extracting the distribution of winter wheat and rape in Hubei province[J].Chinese Journal of Agrometeorology, 2012,33(4):573-578.

[15]權(quán)文婷,王釗. 冬小麥種植面積遙感提取方法研究[J]. 國土資源遙感,2013,25(4):8-15.Quan W T, Wang Z. Researches on the extraction of winter wheat planting area using remote sensing method[J].Remote Sensing for Land and Resources, 2013,25(4):8-15.

[16]楊紅衛(wèi),童小華. 中高分辨率遙感影像在農(nóng)業(yè)中的應(yīng)用現(xiàn)狀[J].農(nóng)業(yè)工程學(xué)報(bào),2012,28(24):138-149.Yang H W, Tong X H. Application status of middle and high resolution remote sensing images in agriculture[J].Transactions of the Chinese Society of Agricultural Engineering,2012,28(24):138-149.

[17]宋明明,鄭文龍,卞國棟. 基于多時相Landsat8數(shù)據(jù)的南方丘陵區(qū)典型地物信息提取—以秦淮河流域?yàn)槔齕J]. 科學(xué)技術(shù)與工程,2016,16(18):211-216.Song M M, Zheng W L, Bian G D. Typical southern hillside area information extraction based on multiple time Landsat8 data[J].Science Technology and Engineering, 2016,16(18):211-216.

[18]韓婷婷,習(xí)曉環(huán),王成. 基于決策樹方法的云南省森林分類研究[J]. 遙感技術(shù)與應(yīng)用,2014,29(5):744-751.Han T T, Xi X H, Wang C. Study on forest classification in Yunnan baesd on decision tree algorithm[J].Remote Sensing Technology and Application, 2014,29(5):744-751.

[19]苗翠翠. 基于NDVI時間序列的水稻面積提取研究[D]. 南京:南京林業(yè)大學(xué),2010.Miao C C. Rice area extraction based on NDVI time seris[D].Nanjing: Nanjing Forest University, 2010.

[20]郝衛(wèi)平,梅旭榮,蔡學(xué)良. 基于多時相遙感影像的東北三省作物分布信息提取[J]. 農(nóng)業(yè)工程學(xué)報(bào),2011,27(1):201-207.Hao W P, Mei X R, Cai X L. Crop planting extraction based on multi-temporal remote sensing data in northeast China[J].Transactions of the Chinese Society of Agricultural Engineering,2011,27(1):201-207.

[21]徐涵秋,唐菲. 新一代Landsat系列衛(wèi)星:Landsat8遙感影像新增特征及其生態(tài)環(huán)境意義[J]. 生態(tài)學(xué)報(bào),2013,33(11):3249-3257.Xu H Q, Tang F. Analysis of new characteristics of the first Landsat 8 image and their eco-environmental significance[J].Acta Ecologica Sinica, 2013,33(11):3249-3257.

[22]Pringle M J, Denham R J, Devadas R. Identification of cropping activity in central and sourthern Queensland, Australia, with the aid of MODIS MOD13Q1 imagery[J].International Journal of Applied Earth Observation and Geoinformation, 2012,19:276-285.

[23]汪擴(kuò)軍. 農(nóng)業(yè)氣候業(yè)務(wù)[M]. 北京:氣象出版社,2013:183-188.Wang K J. Agricultural meteorological operations[M]. Beijing:Meteorological Press, 2013:183-188.

[24]Oliver M A. Kringing: A method of interpolation for geographical information system[J].International Journal of Geographic Information Systems, 1990,49(4):313-332.

[25]李軍龍,張劍,張叢. 氣候要素空間插值方法的比較分析[J]. 草業(yè)科學(xué),2006,23(8):6-11.Li J L, Zhang J, Zhang C. Analysis and compare the spatial interpolation methods for climate factor[J].Pratacultural Science,2006,23(8):6-11.

[26]阿力木江,朱秀芳,努爾麥麥提. 基于時序歸一化植被指數(shù)的撂荒地調(diào)查研究[J]. 上海國土資源,2016,37(4):34-37.Alimujiang, Zhu X F, Nuermaimaiti. Studying abandoned farmland using temporal NDVI data[J].Shanghai Land & Resources,2016,37(4):34-37.

[27]涂曉揚(yáng). 福建省植被指數(shù)與地形要素多尺度空間格局分析[J]. 上海國土資源,2016,37(2):9-12,20.Tu X Y. Multi-scale features of spatial patterns of the normalised difference vegetation index and topographic factors in Fujian province[J].Shanghai Land & Resources, 2016,37(2):9-12,20.