東北黑土區(qū)典型縣域種植模式遙感識(shí)別與空間格局分析

杜國明，張 瑞，梁常安，胡明宇

（1. 東北農(nóng)業(yè)大學(xué)公共管理與法學(xué)院，哈爾濱 150030；2. 東北農(nóng)業(yè)大學(xué)經(jīng)濟(jì)管理學(xué)院，哈爾濱 150030）

0 引 言

東北黑土區(qū)是全國重要的大豆和玉米主產(chǎn)區(qū)，在保障國家糧食安全方面具有舉足輕重的作用。黑土地雖然土地肥沃、糧食產(chǎn)量較高，但近些年來也面臨著土壤侵蝕加劇、黑土層變薄、土壤養(yǎng)分流失等生態(tài)環(huán)境問題[1]。2020年，習(xí)近平總書記在吉林考察時(shí)提出“保護(hù)好黑土地這一‘耕地中的大熊貓’，留給子孫后代”。作物種植模式對水土保持及黑土資源可持續(xù)利用具有重要影響，已有研究表明，相比于連作，豆米輪作可以明顯改善黑土土壤理化性狀，如提高土壤有機(jī)質(zhì)含量、增加微生物群落多樣性和減少病蟲害影響等[2-5]，進(jìn)而促進(jìn)作物生長發(fā)育，實(shí)現(xiàn)大豆和玉米增產(chǎn)[6-7]，是一種能夠均衡用地養(yǎng)地、協(xié)調(diào)經(jīng)濟(jì)效益和生態(tài)效益的種植模式。但在2007－2016年間，由于國家玉米臨儲(chǔ)和各項(xiàng)支農(nóng)惠農(nóng)補(bǔ)貼政策的實(shí)施，在市場作用下，東北地區(qū)糧食生產(chǎn)結(jié)構(gòu)呈現(xiàn)“玉米增、大豆減”的時(shí)空演變基調(diào)[8]，傳統(tǒng)的“大豆-玉米”輪作種植模式受到嚴(yán)重破壞[9]。為促進(jìn)農(nóng)業(yè)可持續(xù)發(fā)展，國家于2016年出臺(tái)《探索實(shí)行耕地輪作休耕制度試點(diǎn)方案》，方案提出在東北冷涼區(qū)開展輪作試點(diǎn)，并建立了相應(yīng)的輪作補(bǔ)助發(fā)放制度。因此，開展東北黑土區(qū)作物種植模式研究，對于探明區(qū)域輪作水平、優(yōu)化輪作模式以及加強(qiáng)黑土地保護(hù)，推進(jìn)“藏糧于地”戰(zhàn)略措施和質(zhì)量興農(nóng)戰(zhàn)略實(shí)施具有重要意義。

種植模式是指某一地區(qū)或某一耕地地塊在不同復(fù)種模式下作物的種植順序，主要包括連作模式和輪作模式，其中連作是指在同一地塊連季或連年種植相同作物的種植方式，而輪作則是指在同一地塊上有順序地在季節(jié)間或年際間輪換種植不同作物的種植方式[10]。多熟種植地區(qū)的種植模式是季節(jié)間的，由于復(fù)種模式復(fù)雜，學(xué)者們多從高時(shí)間分辨率遙感數(shù)據(jù)中提取時(shí)間序列植被指數(shù)，并根據(jù)植被指數(shù)的變化規(guī)律判斷作物種植順序[11]，進(jìn)而揭示種植模式的空間格局特征。如張霞等[12]利用全年23個(gè)時(shí)相的MODIS EVI時(shí)間序列數(shù)據(jù)結(jié)合作物物候信息制作了華北平原兩熟制地區(qū)的種植模式分布圖。顧曉鶴等[13]從時(shí)間跨度為1.5 a的9幅遙感影像中提取NDVI時(shí)間序列信息，并采用變化向量分析法識(shí)別了北京順義地區(qū)的耕地輪作模式。許青云等[14]則利用460個(gè)時(shí)相的影像建立時(shí)間跨度為10 a的NDVI 長時(shí)間序列數(shù)據(jù)集和農(nóng)作物季相節(jié)律特征識(shí)別了陜西省各年度的作物種植模式。郭昱杉等[15]基于MODIS 8 d NDVI時(shí)序數(shù)據(jù)識(shí)別出黃河三角洲地區(qū)的主要種植模式為冬小麥-夏玉米。Nguyen等[16]通過重建16 a期的MODIS EVI曲線以識(shí)別種植模式，進(jìn)而劃分出越南中部Dak Lak省的單季和雙季玉米種植區(qū)。不同于多熟地區(qū)，一年一熟地區(qū)的種植模式是年際間的，低復(fù)種指數(shù)下僅需1個(gè)時(shí)相的影像即可實(shí)現(xiàn)作物類型提取，而后通過疊加方式獲取種植模式是學(xué)者們普遍采用的思路。如劉佳等[17]通過疊加2個(gè)相鄰年度的作物空間分布數(shù)據(jù)識(shí)別了黑龍江省海倫市的年際間作物變化信息。于鳳榮等[18]以黑龍江省的友誼農(nóng)場為研究區(qū)，在3 a期作物分類數(shù)據(jù)的基礎(chǔ)上統(tǒng)計(jì)了各年度的作物面積，使用重茬率、迎茬率、重茬相對變化率等數(shù)理統(tǒng)計(jì)指標(biāo)構(gòu)建了作物重迎茬面積監(jiān)測方法?？傮w來看，多熟地區(qū)和一年一熟地區(qū)種植模式的識(shí)別方法是有較大差別的，黑土區(qū)作為典型的一年一熟地區(qū)，部分學(xué)者已對該地區(qū)的種植模式監(jiān)測進(jìn)行了初步探索，但局限于2個(gè)相鄰年度間作物變化信息的識(shí)別和重迎茬面積的監(jiān)測。一方面，種植模式不能簡單地等同于作物變化信息，它是后者的進(jìn)一步分類和綜合，而這個(gè)過程是以豐富的作物變化信息為基礎(chǔ)的；另一方面，種植模式格局是在自然地理因素和社會(huì)經(jīng)濟(jì)因素共同作用下形成的，除面積信息外，其空間分布信息亦是不能忽略的。因此，適當(dāng)增加種植模式監(jiān)測研究的時(shí)間跨度并使用能夠保留種植模式空間位置信息的方法對于揭示黑土區(qū)種植模式的類型及其空間格局是有必要的。

本文提出一種基于地學(xué)信息圖譜方法的作物種植模式監(jiān)測方法，以黑龍江省克東縣作為研究區(qū)域，利用6個(gè)時(shí)相的遙感影像實(shí)現(xiàn)6 a間的作物分布提取，疊加之后進(jìn)行圖譜重構(gòu)以識(shí)別種植模式，最后使用核密度估計(jì)法分析了各類種植模式的空間格局特征，以期為優(yōu)化輪作方案、制定輪作制度提供科學(xué)依據(jù)。

1 研究區(qū)概況

克東縣地處小興安嶺和松嫩平原的過渡地帶，介于126°01′～126°39′E，47°42′～48°17′N之間，隸屬于黑龍江省齊齊哈爾市，東至北安市，南鄰拜泉縣，西與克山縣毗鄰，北與五大連池市相望，如圖1所示。

克東縣下轄4鎮(zhèn)3鄉(xiāng)1農(nóng)場，共包括99個(gè)行政村、4個(gè)農(nóng)場、3個(gè)林場（為便于描述，以下統(tǒng)稱為行政村）?？h域內(nèi)地勢丘陵起伏，緩坡漫崗，海拔高度147～411 m，屬于中溫帶大陸性季風(fēng)氣候，年平均氣溫為2.4 °C，年平均降水量為526.5 mm，降水主要集中出現(xiàn)在7－8月，占全年降水量的49.9%。全縣土地總面積為208 324.29 hm2，其中耕地面積為147 084.53 hm2，墾殖率達(dá)70.60%。主要土壤類型為黑土、草甸土、暗棕壤、沼澤土等，其中，黑土覆蓋面積最大，占耕地總面積的75.61%?？藮|縣位于東北典型黑土區(qū)腹地，土質(zhì)肥沃，適宜農(nóng)業(yè)耕作，農(nóng)業(yè)經(jīng)營主體以農(nóng)戶為主，作物熟制為一年一熟，主要種植大豆和玉米?？藮|縣種植業(yè)較為發(fā)達(dá)，2016年，全縣糧食作物播種面積為12.16萬hm2，全年糧豆薯總產(chǎn)量352 327 t，種植業(yè)產(chǎn)值達(dá)到128 089萬元，占全縣生產(chǎn)總值的30%。

2 數(shù)據(jù)來源與研究方法

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

根據(jù)東北地區(qū)主要農(nóng)作物物候特征和已有研究[19]，本文選擇克東縣8月下旬或9月上旬的2012年Landsat7 ETM+和2013－2017年Landsat8 OLI共計(jì)10幅遙感影像作為主要數(shù)據(jù)源（表1）。除此之外，還收集了克東縣第二次土地調(diào)查年度更新數(shù)據(jù)庫、克東縣行政區(qū)劃矢量數(shù)據(jù)以及縣域30 m分辨率的DEM數(shù)據(jù)等作為輔助數(shù)據(jù)。

表1 遙感影像信息統(tǒng)計(jì)Table 1 Remote sensing image information statistics

本文使用ENVI進(jìn)行作物分類，數(shù)據(jù)處理流程如下：1）使用DEM對影像進(jìn)行幾何精校正，并在進(jìn)行輻射定標(biāo)、大氣校正、圖像融合、圖像鑲嵌等一系列預(yù)處理后，使用耕地矢量數(shù)據(jù)建立掩膜對影像進(jìn)行裁剪；2）選取5、4、3波段進(jìn)行標(biāo)準(zhǔn)假彩色影像合成，根據(jù)各類作物的光譜特征，以目視判別方式分別建立大豆、玉米、水稻以及其他作物（雜糧、薯類、漢麻等經(jīng)濟(jì)作物）的解譯標(biāo)志；3）基于解譯標(biāo)志，建立作物分類訓(xùn)練樣本點(diǎn)集。其中，大豆和玉米種植面積較大，所以各選取80個(gè)樣本點(diǎn)，將其隨機(jī)分為50個(gè)訓(xùn)練樣本、30個(gè)驗(yàn)證樣本。水稻和其他作物種植面積較小，所以各選取50個(gè)樣本點(diǎn)，將其隨機(jī)分為30個(gè)訓(xùn)練樣本、20個(gè)驗(yàn)證樣本。為提高解譯的精度，通過計(jì)算不同農(nóng)作物樣本數(shù)據(jù)之間的可分離度來確定訓(xùn)練樣本和驗(yàn)證樣本的選擇是否合理。以2017年8月31日的影像為例，克東縣大豆、玉米、水稻和其他作物的訓(xùn)練樣本可分度見表2。樣本可分離度的取值范圍為[0，2]，若可分離度大于1.9，表明不同作物樣本之間的可分離性較好，樣本的選擇較為合理；若可分離度小于1.8，則表明不同作物樣本之間的可分離性較差，樣本需要重新選擇；若可分離度小于1則樣本可以進(jìn)行合并。從結(jié)果來看（表2），除了大豆和水稻的可分離度小于1.9之外，其他的農(nóng)作物樣本之間的可分離性均較高，說明選擇的農(nóng)作物樣本較為合理；4）使用監(jiān)督分類方法對2012－2017年克東縣各類作物進(jìn)行分類。

表2 訓(xùn)練樣本可分離度統(tǒng)計(jì)Table 2 Statistics of training samples reparability

2.2 研究方法

2.2.1 地學(xué)信息圖譜

地學(xué)信息圖譜分析方法由陳述彭院士在20世紀(jì)90年代提出，它綜合了景觀圖的簡潔性和數(shù)學(xué)模型的抽象性，是一種經(jīng)典的地理時(shí)空復(fù)合分析方法[20]。圖譜具有要素“圖形”和“譜系”的雙重性質(zhì)，“圖”可以直觀地表達(dá)要素的空間分布特征，而“譜”則可以反映要素的過程變化信息，圖譜結(jié)合即可綜合展示要素的時(shí)空演變規(guī)律。借鑒已有研究[21]，地學(xué)信息圖譜中要素柵格單元的計(jì)算公式為

式中C為研究期間要素類型變化的單元屬性值；A為前一階段的要素單元屬性值；B為后一階段要素單元屬性值。

本文以2012－2017年間的6期作物分類柵格數(shù)據(jù)為基礎(chǔ)，首先在ArcGIS中對其進(jìn)行重分類處理，將表征大豆、玉米、其他作物和水稻的柵格單元屬性值分別編碼為1、2、3、4，而后采用ArcGIS空間分析模塊中的“地圖代數(shù)”功能進(jìn)行空間疊加，即可得到克東縣2012－2017年間的作物種植模式柵格數(shù)據(jù)。具體計(jì)算公式為

式中G為表征研究期間內(nèi)作物種植信息變化特征的圖譜單元編碼值；Y2012、Y2013、Y2014、Y2015、Y2016、Y2017分別為表征2012、2013、2014、2015、2016、2017共6個(gè)年度農(nóng)作物種植信息的圖譜單元編碼值。

2.2.2 核密度估計(jì)法

核密度估計(jì)法是在統(tǒng)計(jì)樣本數(shù)據(jù)的基礎(chǔ)上，通過獲取數(shù)據(jù)的概率分布曲線來擬合已知數(shù)據(jù)點(diǎn)，從而實(shí)現(xiàn)由已知估計(jì)未知的密度函數(shù)，屬于非參數(shù)檢驗(yàn)方法。該方法通過計(jì)算要素在其周圍鄰域中的密度，可以反映點(diǎn)群的空間分布熱點(diǎn)、密度及其趨勢等。本文通過計(jì)算各類種植模式地塊的核密度，以展示種植模式的空間集聚特征。計(jì)算公式為

式中f(x)為核密度值；n為已知點(diǎn)數(shù)量；h為帶寬；k為核函數(shù)；x為估計(jì)點(diǎn)坐標(biāo)值；xi為樣本點(diǎn)坐標(biāo)值。

在核密度估計(jì)中，不同的帶寬h對密度分布結(jié)果有顯著影響，隨著h的增大，空間上點(diǎn)密度的變化更為光滑且概化程度更高；h減小時(shí)，估計(jì)點(diǎn)密度變化會(huì)顯的凹凸不平[22]。本文采用核密度估計(jì)法測度克東縣各類種植模式地塊的空間集聚特征。首先使用ArcGIS中“柵格轉(zhuǎn)點(diǎn)”功能將表征種植模式的柵格轉(zhuǎn)換為點(diǎn)狀要素，而后以此點(diǎn)狀要素為基礎(chǔ)進(jìn)行核密度值的測算。

3 結(jié)果與分析

3.1 作物分類精度評(píng)價(jià)

由于目前國內(nèi)尚無針對東北地區(qū)種植模式的遙感提取結(jié)果作為參考，無法直接評(píng)價(jià)本文種植模式的提取精度。本文以混淆矩陣的形式使用驗(yàn)證樣本對作物分類結(jié)果進(jìn)行精度評(píng)價(jià)，評(píng)價(jià)指標(biāo)包括總體分類精度、Kappa系數(shù)、用戶精度、制圖精度4項(xiàng)，以2017年為例（表3）?？藮|縣2017年作物分類總體精度為92%，Kappa系數(shù)為0.89，大豆、玉米、水稻和其他作物的用戶精度分別為96.30%、93.33%、90.91%、85.71%，這表明分類精度滿足后續(xù)分析要求。

表3 克東縣作物分類混淆矩陣Table 3 Confusion matrix of crop classification in Kedong County

3.2 克東縣2012－2017年作物空間布局

根據(jù)作物分類結(jié)果對克東縣2012－2017年各類作物的種植面積進(jìn)行統(tǒng)計(jì)，如圖2所示。

由圖2可知，2012－2017年，克東縣大豆與玉米的種植面積之和占總面積的比例均高于94%，水稻和其他作物的種植面積占比均較低且變化相對穩(wěn)定。從大豆和玉米種植面積的年度變化來看，2013－2015年，大豆種植面積由9.20×104hm2下降至5.75×104hm2，玉米種植面積則由4.64×104hm2增加到8.28×104hm2。從2016年開始，玉米的種植面積顯著下降，由2015年的8.28×104hm2下降至2017年的3.74×104hm2，而大豆種植面積持續(xù)升高，由2015年的5.75×104hm2上升到2017年的1.01×105hm2，達(dá)到6 a間的最高值。

由于相鄰年度的作物空間分布情況相似，本文對主要作物種植面積的變化轉(zhuǎn)折點(diǎn)和極值點(diǎn)（2013年、2015年、2017年）的空間分布情況進(jìn)行可視化處理（圖3）?？梢姡蠖?、玉米種植地塊在縣域廣泛分布，6 a間兩者的分布區(qū)域不斷變化且呈“互為進(jìn)退”關(guān)系；水稻在縣域北部呈東西向沿河流谷地帶狀分布；其他作物則在縣域內(nèi)零散分布。

3.3 克東縣種植模式類型

使用地學(xué)信息圖譜方法，將6 a期作物分類柵格數(shù)據(jù)進(jìn)行疊加之后，獲取作物變化信息圖譜。經(jīng)統(tǒng)計(jì)，作物變化信息圖譜單元有2 857種編碼，依據(jù)種植模式的概念以及編碼特征進(jìn)行分類[10]，即圖譜重構(gòu)，大體可以歸并為6類種植模式（表4）。其中，僅由大豆和玉米兩種作物單元疊加產(chǎn)生的編碼有64種，可以發(fā)現(xiàn)：

1）表征某地塊在6 a間至少連續(xù)4 a種植大豆的編碼有8種，將其統(tǒng)一命名為“大豆連作模式”，此類模式意味著地塊的大豆連作問題較為突出；

2）表征某地塊在6 a間至少連續(xù)4 a種植玉米的編碼有8種，將其統(tǒng)一命名為“玉米連作模式”，此類模式則表示地塊的玉米連作障礙現(xiàn)象較為明顯；

3）表征某地塊在6 a間無明顯周期性豆米輪作規(guī)律的編碼有27種，本文將其統(tǒng)一命名為“無序種植模式”，此類模式在一定程度上實(shí)現(xiàn)了作物的交替種植；

4）表征某地塊在6 a間以3 a為周期，進(jìn)行2個(gè)周期的以玉米為基礎(chǔ)茬口、大豆為中軸作物的“玉米-玉米-大豆”或以大豆為基礎(chǔ)茬口、玉米為中軸作物的“大豆-大豆-玉米”的種植模式編碼有6種，將其統(tǒng)一命名為“三年輪作模式”，這是耕地輪作的另一種形式；

5）表征某地塊在6 a間以2 a為周期，至少進(jìn)行2個(gè)周期的“大豆-玉米”的交替種植的編碼有15種，將其統(tǒng)一命名為“兩年輪作模式”，此類模式是東北地區(qū)傳統(tǒng)的豆米輪作模式。

除此之外，共有2 793種編碼中至少含有一個(gè)代碼3或4，這表征某地塊在6 a間至少有1 a種植水稻或其他作物，將其統(tǒng)一命名為“混合種植模式”，其中包含水稻連作、水旱輪作、雜糧輪作等眾多種植模式。

根據(jù)種植模式分類模型，將作物變化信息圖譜中表征大豆連作模式、玉米連作模式、無序種植模式、三年輪作模式、兩年輪作模式、混合種植模式等6類種植模式的地塊進(jìn)行制圖（圖4）。由于混合種植模式編碼類型較為繁雜，且面積占比（21.48%）相對較低，為簡化分析，本文以下內(nèi)容對混合種植模式不做過多討論，僅選擇表4中種植模式二級(jí)類的前五項(xiàng)作為主要研究對象?？藮|縣各類種植模式的面積由大到小分別為無序種植模式、大豆連作模式、兩年輪作模式、玉米連作模式、三年輪作模式，其中，無序種植模式、大豆連作模式以及兩年輪作模式的面積占比之和為83.90%，是克東縣主要的種植模式。

表4 克東縣種植模式分類模型Table 4 Classification model of cropping patterns in Kedong County

3.4 克東縣種植模式的空間分布特征

本文采用核密度估計(jì)法測度克東縣各類種植模式地塊的空間分布。在帶寬的選擇上，以表征大豆連作模式的點(diǎn)狀要素為例，本文結(jié)合核密度分析工具中自動(dòng)生成的1 600 m帶寬，分別設(shè)置1 000、2 000、3 000、4 000、5 000 m帶寬進(jìn)行對比研究，據(jù)此比選最優(yōu)帶寬。由于要體現(xiàn)縣域尺度上種植模式地塊的整體分布趨勢，因此選擇較大的帶寬，經(jīng)過多次試驗(yàn)，確定最優(yōu)帶寬為4 km。在此帶寬下，既能較為清晰地辨別出各類種植模式地塊的密度中心，也能較完整地覆蓋各類種植模式地塊的分布區(qū)域（圖5）。

大豆連作模式地塊主要分布于縣域西部和北部，包括寶泉鎮(zhèn)、金城鄉(xiāng)以及千豐鎮(zhèn)；玉米連作模式地塊在縣域內(nèi)呈現(xiàn)東北-西南走向的帶狀分布特征，其中以玉崗鎮(zhèn)最為明顯；無序種植模式地塊在縣域中部、東部和南部呈現(xiàn)集中連片分布的形態(tài)，在西部的金城鄉(xiāng)和北部的寶泉鎮(zhèn)亦有一定集聚分布；三年輪作模式地塊和兩年輪作模式地塊在縣域內(nèi)均呈現(xiàn)“局部集聚、全域分散”的分布特征，同時(shí)，兩者在趙光農(nóng)場北部的分布趨勢較為一致。將各類種植模式的分布區(qū)域進(jìn)行對比，可以發(fā)現(xiàn)，無序種植模式地塊與大豆連作模式地塊存在一定的空間互補(bǔ)關(guān)系，與其他3類種植模式地塊的分布區(qū)域卻多有重疊，加之無序種植模式為縣域內(nèi)面積占比最大的種植模式。因此，在縣域中部、東部和南部，無序種植模式可以視為克東縣種植模式的“基底”；對于面積占比為第二位的大豆連作模式地塊，與玉米連作模式地塊的主要分布區(qū)域呈現(xiàn)“互斥但不互補(bǔ)”的空間分布關(guān)系，與三年輪作模式地塊和兩年輪作模式地塊則亦有重疊。因此，在縣域西部和北部，大豆連作模式可以視為克東縣種植模式的“基底”。

3.5 克東縣種植模式的空間結(jié)構(gòu)特征

為進(jìn)一步分析各行政村的種植模式類型及比例，本文選取各行政村內(nèi)面積占比為前3位的種植模式，并按其所占比例由高到低進(jìn)行排列組合形成該村的種植模式結(jié)構(gòu)，以展示行政村尺度種植模式的空間結(jié)構(gòu)特征。通過計(jì)算各行政村前3位種植模式的面積之和占村域內(nèi)所有種植模式面積的比例，得出克東縣所有行政村上述比例的總體平均值為88.29%，這證明此方法是適用的。為便于表達(dá)和制圖，將大豆連作模式、玉米連作模式、無序種植模式、三年輪作模式、兩年輪作模式分別賦予代碼“S、M、D、T、W”，最終將106個(gè)行政村歸并為10類種植模式結(jié)構(gòu)（表5）。

表5 克東縣種植模式結(jié)構(gòu)分區(qū)Table 5 Zoning of the cropping pattern structure in Kedong County

為便于進(jìn)一步分析，選取行政村數(shù)量占比為前5位的種植模式結(jié)構(gòu)（代碼為D-S-W、S-D-W、D-M-W、D-W-S、D-W-M）作為縣域的主要種植模式結(jié)構(gòu)類型，并根據(jù)相似性原則將數(shù)量較少的S-D-T歸入S-D-W，D-S-M和D-S-T歸入D-S-W中，D-M-S和D-M-T歸入D-M-W中，可視化結(jié)果見圖6。由圖6可知，“無序種植-大豆連作-兩年輪作”種植模式結(jié)構(gòu)（代碼為D-S-W）呈西北-東南向在縣域內(nèi)連片分布，包含37個(gè)行政村，是縣域內(nèi)占比最大的種植模式結(jié)構(gòu)類型。該區(qū)域年際間的作物種植隨意性較高，且大豆連作障礙現(xiàn)象相對明顯，輪作以兩年輪作模式為主但在該區(qū)不占優(yōu)勢，這與克東縣整體種植模式結(jié)構(gòu)是一致的；“大豆連作-無序種植-兩年輪作”種植模式結(jié)構(gòu)（代碼為S-D-W）零散分布于縣域西部和北部，在西南部則呈集聚分布特征，包含22個(gè)行政村?？傮w來看，該地區(qū)大豆連作問題比較突出，農(nóng)戶種植行為亦缺乏引導(dǎo)，輪作政策實(shí)施效果較差；“無序種植-玉米連作-兩年輪作”種植模式結(jié)構(gòu)（代碼為D-M-W）在縣域內(nèi)呈東北-西南向帶狀分布，包含21個(gè)行政村。該區(qū)域作物種植隨意性亦較高，玉米連作障礙現(xiàn)象相對明顯，兩年輪作模式有一定分布但不占優(yōu)勢；“無序種植-兩年輪作-大豆連作”種植模式結(jié)構(gòu)（代碼為D-W-S）在縣域中部呈南北向離散分布，包含16個(gè)行政村。該區(qū)域的作物種植缺乏引導(dǎo)，但兩年輪作模式的推廣效果相對較好，大豆連作現(xiàn)象占較小比例；“無序種植-兩年輪作-玉米連作”種植模式結(jié)構(gòu)（代碼為D-W-M）在縣域中部、南部和東部零散分布，涉及10個(gè)行政村。該區(qū)域作物種植無序特征相對明顯，兩年輪作效果較為突出，但玉米連作障礙現(xiàn)象同樣存在。

4 討 論

農(nóng)戶的作物種植行為受農(nóng)產(chǎn)品價(jià)格政策、市場價(jià)格、輪作種植等因素影響較大[9]。2013－2015年，克東縣大豆種植面積下降，玉米種植面積上升，這可能是受到大豆進(jìn)口價(jià)格持續(xù)走低的影響，加之玉米單產(chǎn)及均價(jià)均高于大豆，導(dǎo)致農(nóng)戶大豆種植意愿較低。從2015年開始，國家為調(diào)整糧食生產(chǎn)結(jié)構(gòu)，鼓勵(lì)農(nóng)民種豆積極性，在東北三省下調(diào)臨儲(chǔ)玉米掛牌收購價(jià)格，并調(diào)減“鐮刀灣”地區(qū)的玉米種植面積，造成玉米的種植面積顯著縮減，大豆種植面積快速攀升，調(diào)控政策富有成效。

克東縣無序種植模式占比最大，這可能是由于農(nóng)戶作為“理性經(jīng)濟(jì)人”，其作物種植決策行為易受農(nóng)產(chǎn)品經(jīng)濟(jì)收益影響而不斷調(diào)整；由于該地農(nóng)戶已具備成熟的大豆種植技術(shù)和相對完善的機(jī)械設(shè)施，種植玉米反而會(huì)導(dǎo)致收益下降，故大豆連作現(xiàn)象較為嚴(yán)重；兩年周期的豆米輪作是當(dāng)?shù)貍鹘y(tǒng)的輪作方式，農(nóng)戶出于保護(hù)耕地的目的，使得兩年輪作模式的面積占比相對較高。

與劉佳等[17]的研究相比，本文在獲取作物變化信息的基礎(chǔ)上進(jìn)一步識(shí)別出了種植模式。與于鳳榮等[18]的研究相比，本文除了分析種植模式的面積，還對種植模式的空間分布特征進(jìn)行了測度。地學(xué)信息圖譜作為一種經(jīng)典的地學(xué)分析方法論，在土地利用變化、景觀格局演變、生態(tài)信息表達(dá)等方面具有廣泛的應(yīng)用[23]，但是在農(nóng)業(yè)領(lǐng)域的應(yīng)用研究相對較少。同時(shí)，種植模式是一種典型的時(shí)空復(fù)合體[24]，這與地學(xué)信息圖譜可以綜合表達(dá)空間位置和變化過程信息的內(nèi)涵是一致的。因此，使用圖譜方法進(jìn)行種植模式的識(shí)別是可行的。

為實(shí)現(xiàn)土壤改良、牧業(yè)發(fā)展、生態(tài)保護(hù)、增加糧油供給等目標(biāo)，農(nóng)業(yè)部、中央農(nóng)辦等于2016年發(fā)布《探索實(shí)行耕地輪作休耕制度試點(diǎn)方案》，提出在東北冷涼區(qū)推廣“一主四輔”輪作模式，但未明確各類模式的具體輪作區(qū)域。因此，立足東北地區(qū)不同區(qū)域的資源條件和生態(tài)特點(diǎn)，探索各地區(qū)適宜的輪作作物、輪作格局、輪作周期并構(gòu)建多目標(biāo)協(xié)同的輪作區(qū)劃應(yīng)當(dāng)是下一步研究的重點(diǎn)[3]。

作物空間分布格局包含種植結(jié)構(gòu)、種植面積等信息，是區(qū)域種植模式形成的基礎(chǔ)[25]。反之，在區(qū)域理想的輪作模式框架下，如何在空間維度上合理安排作物布局，在時(shí)間維度上協(xié)調(diào)作物種植順序，進(jìn)而推動(dòng)區(qū)域種植業(yè)結(jié)構(gòu)調(diào)整和優(yōu)化東北地區(qū)專業(yè)作物帶，對于確保國家糧食安全意義重大。

對于一年一熟地區(qū)種植模式的后續(xù)研究應(yīng)注意選取適宜的時(shí)間尺度。過小的時(shí)間尺度難以識(shí)別出周期性的輪作規(guī)律，如2 a的時(shí)間尺度僅能體現(xiàn)作物變化信息；同時(shí)，農(nóng)戶的作物種植決策行為亦在不斷調(diào)整中，多周期的輪作是難以實(shí)現(xiàn)的，過大的時(shí)間尺度下，種植行為的無序特征將愈加明顯?？傮w來看，時(shí)間尺度在3～6 a之間選取是較為合適的。

5 結(jié) 論

本文利用Landsat 8 OLI遙感數(shù)據(jù)提取克東縣作物類型，并結(jié)合地學(xué)信息圖譜方法識(shí)別種植模式，進(jìn)而分析了種植模式的類型、面積、空間集聚性和空間結(jié)構(gòu)等空間格局特征。主要結(jié)論如下：

1）2012－2017年，克東縣大豆和玉米種植總面積占比均高于94%，且兩者的種植面積在6 a間分別呈現(xiàn)“先減后增”和“先增后減”的變化趨勢。水稻和其他作物的種植面積占比較低且變化相對穩(wěn)定。

2）克東縣各類作物種植模式按面積由大到小排序?yàn)闊o序種植模式、大豆連作模式、兩年輪作模式、玉米連作模式、三年輪作模式。其中，前3位種植模式面積占比之和為83.90%，是克東縣主要的種植模式。

3）大豆連作模式在縣域西部和北部占有明顯優(yōu)勢，無序種植模式在縣域中部、東部和南部則是主要的種植模式。玉米連作模式呈東北-西南向帶狀分布，三年輪作模式和兩年輪作模式呈“局部集聚、全局分散”的分布形態(tài)。

4）行政村尺度的種植模式結(jié)構(gòu)大體可以分為5種類型，其中以“無序種植-大豆連作-兩年輪作”類型為主，涉及的行政村廣泛分布于在縣域西北-東南方向；“大豆連作-無序種植-兩年輪作”和“無序種植-玉米連作-兩年輪作”種植模式結(jié)構(gòu)次之，兩者占比相當(dāng)，前者在縣域東部零散分布，后者呈東北-西南向帶狀分布；“無序種植-兩年輪作-大豆連作”和“無序種植-兩年輪作-玉米連作”兩類涉及的行政村零散分布于縣域東部和東南部。