基于局域空間自相關的吉林省耕地質量空間分布分析

周巖

摘要將耕地質量評價的自然等、利用等和經濟等作為耕地質量圖斑的空間屬性，應用局域空間自相關分析方法，對耕地質量空間分布格局進行分析，探討吉林省耕地質量的空間分布特征，并對耕地質量評價結果的合理性進行探討。結果表明，吉林省耕地質量整體上表現出顯著的聚集性；HH型和LL型表現為正相關的聚集分布；HL型和LH型為負相關的零星分布。最終通過對各種類型的分布特征分析耕地質量評價結果的合理性，為耕地質量評價更新和后續管理提供支持。

關鍵詞耕地質量；局域空間自相關；空間分布格局

中圖分類號F301.2文獻標識碼A文章編號0517-6611（2017）32-0200-04

Spatial Distribution of Cropland Quality in Jilin Province Based on Local Spatial Autocorrelation

ZHOU Yan

（Land Consolidation and Rehabilitation Center of Jilin Province， Changchun， Jilin 130061）

AbstractSpatial attributes of cropland patches were set as natural grade， utilization grade and economic grade of cropland quality assessment. With local spatial autocorrelation method， the spatial distribution of cropland quality was analyzed and the characteristic of spatial distribution and the rationality of results of cropland quality assessment were discussed. The results showed that significant aggregation was one of the characteristics of cropland quality of Jilin Province. HH and LL types were positive correlation and the distribution characteristics were aggregated. HL and LH types were negative correlation and the distribution characteristics were scattered. Through analyzing the rationality of the results of cropland quality assessment， the method and the results of this study can be used to support the update assessment and management of cropland quality in the future.

Key wordsCropland quality；Local spatial autocorrelation；Spatial distribution pattern

耕地質量是農業生產的最重要保障，保持與逐步提高耕地質量是耕地管理中的重要內容之一[1]。開展耕地質量等別調查與評定，全面掌握和科學量化耕地質量等別狀況，是深化國土資源管理與改革，促進國土資源科學化、規范化管理的重要基礎工作。從耕地質量評價的內容來看，主要的研究方法是立足于耕地質量與立地條件分析的耕地質量評價[2-5]。以美國的農用地保護的LESA（land evaluation and site assessment）體系思想為指導，構建耕地質量與立地條件綜合分析體系，采取多因素綜合對比分析法，開展耕地建設的適宜性評價，并劃分耕地質量級別。隨著遙感和地理信息技術的發展，在GIS和RS技術的支持下，構建了耕地質量評價模型和評價指標體系，對區域耕地質量等別進行了廣泛的研究，并對耕地質量的空間分布進行了探討[6-9]。由此可見，耕地質量評價的研究內容、指標體系、評價方法等已經基本完善，已經極大地豐富了耕地質量研究。在此基礎上，經過綜合研究與多年探索，國土資源部于2012年頒布了《農用地質量分等規程》（GB/T 28407—2012），指導全國的耕地質量等別評價研究工作，并且已經取得了全國性的耕地質量等別成果。但是，該類成果僅對耕地的數量、質量進行了研究，關于耕地質量的空間分布特征等研究明顯不足。

筆者應用GIS的空間統計、空間分析技術，采用空間自相關分析的方法，對耕地質量的空間分布進行分析，探討耕地質量的空間差異性和結構性特征，為耕地保護保護、基本農田劃定等提供依據，為耕地的差異化管理提供幫助。同時，通過對耕地質量的空間分布特征進行分析，對耕地質量評價結果空間分布的合理性進行判斷，為耕地質量更新評價提供支持。

1數據來源與研究方法

1.1研究區概況

吉林省地處中國東北中部，東北亞地理中心。地跨121°38′～131°19′E、40°50′～46°19′N；吉林省位于日本、俄羅斯、朝鮮、韓國、蒙古與中國東北部組成的東北亞幾何中心地帶。北接黑龍江省，南鄰遼寧省，西接內蒙古自治區，東與俄羅斯聯邦接壤，東南部與朝鮮隔江相望。

面積19.1萬km2，人口2 729.82萬，轄1個自治州、8個地級市、57個縣（市、區）、3個自治縣。吉林省土地資源主要劃分為東部長白山地區和西部松嫩平原區，在結構上表現為以農用地為主體，“三多一少” 的特點：一是農用地多。農用地（包括耕地、林地、草地、農田水利用地、養殖水面等）約1 639.32萬hm2，占土地總面積的86%，高出全國平均水平17百分點。二是耕地多。耕地總量占土地總面積的30%，在全國排第9位。三是耕地后備資源多。全省共有104.6萬hm2未利用地可作為后備資源，其中宜農地30.59萬hm2；四是建設用地少（包括居民點及工礦用地，交通運輸用地和水利設施用地）。總面積約106.53萬hm2，占土地總面積的5.5%。布局上表現為區域差異顯著，大致呈現東林、中農、西牧的格局：東部以林地為主，占東部總面積的81.0%；中部以耕地為主，占中部總面積的61.7%；西部草地、濕地和荒地較多，占西部總面積的36.6%。吉林省還是世界聞名的黑土帶，黑土面積約110.0萬hm2，黑土耕地約83.2萬hm2，占全省耕地面積的15.0%，黑土區糧食產量占全省一半以上。

吉林省是我國的農業大省，作為國家重要的商品糧基地和糧食主產區，吉林省承擔著維護國家糧食安全的重要責任。農業穩則全局穩，穩定的農業生產離不開良好的土地資源，全面掌握吉林省的耕地數量、質量和利用水平，對我國農業可持續發展、糧食安全和社會穩定尤為重要。

1.2數據來源與處理

吉林省土地整治中心按照《國土資源部辦公廳關于部署開展2015年全國耕地質量等別調查評價與監測工作的通知》（國土資廳發〔2015〕17號）的部署，全面開展耕地質量等別評價工作，依據《耕地質量等別年度更新評價技術手冊》有關技術標準，取得了全省耕地質量的等別成果。

由于該耕地質量等別成果是基于吉林省年度變更調查數據庫完成的，耕地質量等別評價單元為1∶1萬耕地圖斑，數據量大。為了便于分析，分別對自然等別、利用等別和經濟等別進行融合，并對離群值進行剔除，將其融合到最鄰近的耕地圖斑中。

1.3研究方法

局域空間自相關是檢驗具有空間位置的某要素的觀測值與周邊地區該觀測值的空間差異程度，即分析其與相鄰空間單元之間空間變量的相關程度。局域空間自相關分析的內容包括被分析對象的空間位置關系以及要素的屬性值[10]。該研究采用ArcGIS 10.2的局域空間自相關指數（Local Morans I）進行分析[10-11]。對于空間單元i，其局域空間自相關指數表達形式為：

Ii=ZijwijZj

Zi=（xi-）/1nni=1（xi-）2

式中，Ii為每個空間單元的莫蘭指數（Moran）；Zi為第i個空間單元的屬性觀測值的標準化值；為空間單元的屬性觀測值的平均值；wij為空間單元之間的權重。根據莫蘭指數和其對應的Z得分，借助Moran散點圖，將計算結果劃分為4種局部空間自相關類型，分別表示為：HH型（高值單元被其他高值單元所包圍），表示該耕地圖斑與周圍圖斑的耕地質量都較高，空間聚集性較強；LH型（低值單元被其他高值單元所包圍），表示該圖斑與周圍圖斑的質量有較強的異質性，空間聚集性不強；LL型（低值單元被其他低值單元所包圍），表示該耕地圖斑與周圍圖斑的耕地質量都較低，空間聚集性較強；HL型（高值單元被其他低值單元所包圍），表示該圖斑與周圍圖斑的質量有較強的異質性，空間聚集性不強。通過對耕地質量圖斑的分布類型進行分析，判斷耕地質量評價結果空間分布的合理性。

空間權重是進行局域空間自相關分析的基礎，用于度量空間單元之間在位置上相互依賴的程度。該研究中，采用空間距離方式進行空間關系的概念化。空間權重矩陣的計算采用反距離權重的方法計算獲得。

2結果與分析

2.1吉林省耕地質量空間分布

吉林省耕地質量等別最后劃分為10個自然等、7個利用等和7個經濟等。從表1可以看出，自然等別中9等、10等、11等地面積最大，為205.17萬、229.29萬、110.16萬hm2，分別占耕地總面積29.29%、32.74%、15.73%；15等地面積較小，為39.61 hm2。從圖1可以看出，松嫩平原區自然質量較優，該區內有面積比例最多的9等、10等和11等地，共占該區面積的83.67%。而12等和13等地則大部分分布在長白山地區，占該區面積的20.06%。長白山地區10等、11等和12等地分布最多，共占該區面積的77.20%，長白山地區沒有6等和7等地、8等和15等地在2區均則呈零星分布。

吉林省耕地利用等為6～14等共9個等別。其中10等、11等和12等地面積最大，為190.25萬、148.93萬、155.88萬hm2，分別占耕地總面積27.17%、21.26%、22.26%；松嫩平原區耕地質量較優，該區內有面積比例最多的9等、10等和11等地，共占該區面積的75.86%。而12等和13等地則大部分分布在長白山地區，占總面積的20.06%。長白山地區10等、11等和12等地分布最多，共占該區面積的77.20%，長白山地區沒有6等和7等地，松嫩平原區沒有15等地，8等和15等地則呈零星分布。

吉林省耕地經濟等為8～14等共8個等別，各等別的面積分布相對較平衡，14等和8等地面積較少，僅有35.15萬hm2，占全省耕地總面積的5.00%；12等地面積為174.13萬hm2，占全省耕地總面積的24.86%。除12等地主要分布在長白山地區外，其他等別的耕地主要分布在松嫩平原區。

從整體空間分布來看，耕地主要分布在中西部平原區，東部山區較少。自然等的較高等別主要分布在公主嶺市、雙遼市、長嶺縣、農安縣、前郭爾羅斯蒙古族自治縣、乾安縣、德惠市、舒蘭市；自然等的較低等別主要分布在東部的延邊朝鮮族自治州、白山市和通化市的大部分地區，西部的松原市部分耕地的自然等別也較低。其他地區的自然等別為中等狀況（圖1a）。

利用等的較高等別主要分布在梨樹縣、公主嶺市、農安縣、長春市、乾安縣、舒蘭市和東遼縣和長嶺縣的部分區域；較低的利用等等別主要分布在東部的延邊朝鮮族自治州、通化市、白山市和西部的通榆縣、洮南市和鎮賚縣。其他地區的利用等別為中等狀況（圖1b）。經濟等的較高等別范圍較少，主要分布在梨樹縣、長春市、農安縣、榆樹市。其他區域的利用等別均較低（圖1c）。

安徽農業科學2017年

2.2吉林省耕地質量局部空間自相關分析

應用ArcGIS 10.2軟件分別以耕地質量自然等、利用等和經濟等等別作為空間屬性，進行局域自相關分析。在95%的置信度下，有一

半的耕地質量自然等在空間上表現為顯著聚集性（50.97%）。圖2a顯示，耕地質量要素之間屬于正相關類型的HH型和LL型分別占總耕地質量圖斑的22.46%和26.52%，累積占到總耕地質量圖斑的48.98%。從空間分布格局來看，LL型（低值與低值聚集分布）的耕地質量圖斑主要分布在東部長白山地區，主要有汪清縣、琿春市、延吉市、龍井市、安圖縣，撫松縣、白山市、長白縣、柳河縣、梅河口市和東豐縣。在敦化市和蛟河市的部分地區有少量分布。HH型（高值與高值聚集分布）的耕地質量圖斑主要分布在公主嶺市、四平市、雙遼市、長嶺縣、前郭爾羅斯蒙古族自治縣、乾安縣、大安市、舒蘭市、吉林市以及鎮賚縣、磐石市和洮南市的部分區域（圖2a）。耕地質量要素之間屬于負相關類型的HL和LH型所占比例非常小，分別僅占總耕地質量圖斑的0.97%和1.02%。從空間分布格局來看，HL型（高值分布于低值區域內）主要呈零散方式分布于東遼縣、梅河口市的部分區域。LH型（低值分布于高值區域內）也呈零散方式分布于蛟河市、吉林市和伊通滿族自治縣的部分區域（圖2a）。

耕地質量利用等在空間上也表現為顯著的聚集性（53.28%，95%的置信區間）。圖2b顯示，呈正相關的HH型和LL型分別占總耕地質量圖斑的21.76%和30.08%。從空間分布格局來看，LL型的耕地質量圖斑主要分布在東部的延邊州地區、白山市地區和通化市地區，其次為東豐縣和鎮賚縣北部部分。HH型的耕地質量圖斑主要分布在中部的梨樹縣、公主嶺市、長春市、德惠市、九臺市、榆樹市、扶余市、遷安市、磐石市和舒蘭市、東遼縣大部、永吉縣、前郭爾羅斯蒙古族自治縣北部和四平市的局部區域。耕地質量圖斑利用等呈負相關的HL型和LH型所占比例非常小，分別為0.35%和1.10%。從其空間分布格局來看，HL型東遼縣西部部分地區和臨江市的部分地區。LH型主要分布在樺甸市、蛟河市、雙陽區、和洮南市的小部分區域，分布范圍零散，所占面積小（圖2b）。

耕地質量經濟等在空間上呈正相關的聚集性稍大于利用等和自然等（57.11%，95%的置信區間）。從圖2c可以看出，HH型主要分布在東部的延邊朝鮮族自治州、白山市和通化市，東豐縣、樺甸市南部和舒蘭市的大部分區域，占耕地質量圖斑的23.88%。LL型主要分布在梨樹縣、長春市、九臺市、德惠市、榆樹市、扶余縣、前郭爾羅斯蒙古族自治縣、長嶺縣、大安市、白城市、永吉縣和吉林市，以及雙遼市、站來西安、磐石市、蛟河市的部分區域，占耕地總面積的31.42%。呈負相關的HL型和LH型分別占耕地質量圖斑總面積的0.67%和1.13%。其中HL型主要分布在柳河縣北部和集安市、通化縣、靖宇縣的小部分地區。LH型主要分布在輝南縣、樺甸市、蛟河市、松原市、通榆縣、洮南縣、白城市和鎮賚縣的小部分區域。

2.3吉林省耕地質量評價結果合理性分析

通過對吉林省耕地質量自然等、利用等和經濟等的局域自相關分析結果進行分析，可以將耕地質量劃分為以下情況：耕地質量呈隨機

分布、異常分布、高值聚集分布和低值聚集分布。局域自相

關分析中，沒有顯著相關類型出現的區域即為耕地質量隨機分布區域。在自然等、利用等和經濟等中所占面積比例為49.03%、46.72%和42.89%。在這些區域中，耕地圖斑的耕地質量等別隨機出現且等別差異小（均在1等別差），說明臨近耕地質量圖斑的等別既具有一定的相似性，又有一定的空間差異，體現了耕地質量的區域一致性，符合耕地質量圖斑的空間分布規律。自然等、利用等和經濟等的高值聚集分布和低值聚集分布區域具有顯著正相關性，說明在這些區域內臨近耕地圖斑的質量等別基本一致，耕地質量評價的因素在該空間范圍內基本無變化。自然等、利用等和經濟等的正相關區域所占面積比例為48.98%、51.83%和55.30%，從這些區域的地理環境特征來看，這種分布特征也符合耕地質量的空間分布規律。自然等、利用等和經濟等的高值被低值包圍或低值被高值包圍的分布區域具有顯著負相關性。該類型分布特征表明，臨近圖斑的耕地質量變化劇烈，差異大，不符合區域一致性規律特征，是一種異常分布情況。自然等、利用等和經濟等的異常區域所占面積比例非常小，分別為1.99%、1.45%和1.80%。針對這些異常分布的情況，要從評價過程入手，檢查自然等、利用等和經濟等的評價因素獲取以及計算過程的正確性。

3結論

該研究對局域空間自相關方法在耕地質量中的應用進行了探索，將吉林省耕地質量自然等、利用等和經濟等作為屬性信息，分析了耕地質量的空間分布特征，并對耕地質量評價結果的合理性進行了分析，得到以下主要結論。

（1）吉林省耕地質量從整體來看，呈現出聚集性特征，并且東部長白山地區耕地質量低于中西部地區。呈正相關類型的區域以集中連片的方式出現，呈負相關類型的區域零星分布，不集中連片。

（2）從耕地質量評價結果角度來看，耕地質量在空間上呈隨機分布的區域，在耕地質量更新評價以及后續管理中要根據局地特征因地制宜地制定管理策略；耕地質量在空間呈現高值聚集的區域，要加強耕地質量的保護，嚴禁耕地占用；耕地質量在空間呈現低值聚集的區域，要結合區域限制因素進行綜合改良，提高耕地質量；耕地質量在空間呈現負相關的區域，在后續更新評價中要對評價因素等基礎信息進行認真搜集，審核其科學性和正確性，查缺補漏，使耕地質量評價結果更加科學、合理。

參考文獻

[1] 沈仁芳，陳美軍，孔祥斌，等.耕地質量的概念和評價與管理對策[J].土壤學報，2012，49（6）：1210-1217.

[2] 錢鳳魁，王秋兵，邊振興，等.凌源市耕地質量評價與立地條件分析[J].農業工程學報，2011，27（11）：325-329.

[3] 夏權，夏萍，馮東，等.安徽省含山縣土地整理耕地質量評價及其變化研究[J].國土資源遙感，2015，27（1）：182-186.

[4] 錢鳳魁，王秋兵，李娜.基于耕地質量與立地條件綜合評價的高標準基本農田劃定[J].農業工程學報，2015，31（18）：225-232.

[5] 沈仁芳，陳美軍，孔祥斌，等.耕地質量的概念和評價與管理對策[J].土壤學報，2012，49（6）：1210-1217.

[6] 李瑜，刁淑娟，邢宇.遙感和GIS技術在耕地質量調查中的應用[J].農業工程，2014，4（3）：86-88，99.

[7]

郭笑怡，高燕，張延玲，等.基于遙感與GIS的三江平原耕地質量評價與生產力分析[J].東北師大學報（自然科學版），2013，45（2）：139-145.

[8] 韓會慶，蔡廣鵬，張鳳太，等.基于GIS的喀斯特地區耕地質量評價：以貴州省綏陽縣為例[J].水土保持研究，2011，18（6）：129-131，135.

[9] 歐陽玲，王宗明，賈明明，等.基于遙感的吉林省中西部耕地數量和質量空間格局變化分析[J].農業工程學報，2016，32（13）：234-242.

[10] 魏曉峰，吳健平.基于ArcGIS的空間自相關分析模塊的開發與應用[J].測繪與空間地理信息，2005，28（6）：77-80.

[11] 郭敏，李淑杰.基于局部空間自相關的耕地質量空間集聚性和保護分區：以吉林省九臺市為例[J].江蘇農業科學，2017，45（3）：206-210.