董秀春，李宗南，蔣 怡，劉忠友，郭 濤，李源洪

（四川省農業(yè)科學院遙感應用研究所，成都610066）

0 引言

隨著城市化水平提高、農村勞動力數(shù)量和結構變化、農業(yè)產(chǎn)業(yè)結構調整，都市郊區(qū)鄉(xiāng)村地區(qū)土地利用發(fā)生明顯變化［1-3］。使用土地利用現(xiàn)狀數(shù)據(jù)分析不同地區(qū)鄉(xiāng)村土地利用變化特征和驅動因素［4-5］，對區(qū)域生產(chǎn)布局調整、鄉(xiāng)村空間規(guī)劃和促進鄉(xiāng)村振興有重要意義［6-7］。已有研究結果顯示京津冀［8］、長三角［9］、珠三角［10］等地的郊區(qū)土地利用變化明顯，總體趨勢大致相同，主要表現(xiàn)為耕地面積減少，建設用地增加。受地理位置、資源稟賦、地方政策引導及產(chǎn)業(yè)發(fā)展差異，郊區(qū)土地利用類型變化存在區(qū)域性差異。如北京郊區(qū)昌平三鎮(zhèn)因城市綠化和兩道環(huán)城綠化隔離帶的建設，林地和草地面積大幅增加［11］；南京市郊區(qū)在社會經(jīng)濟和人口因素的驅動下，林地面積不斷減少，園地和草地面積持續(xù)增加［12］；洛陽市南郊因旅游產(chǎn)業(yè)發(fā)展，風景點游賞用地、商服用地、林草地等土地利用類型占比上升［13］。

成都市作為西南地區(qū)的中心城市，地處成都平原，具有優(yōu)越的地理位置、交通條件和經(jīng)濟發(fā)展等區(qū)位優(yōu)勢。在人口規(guī)模的增加、產(chǎn)業(yè)結構調整以及城市發(fā)展政策的變化等因素影響下，成都市土地利用發(fā)生明顯變化［14］。隨著現(xiàn)代農業(yè)轉型升級，成都市通過科學規(guī)劃帶動郊區(qū)農業(yè)為主導的鄉(xiāng)鎮(zhèn)積極發(fā)展現(xiàn)代農業(yè)，形成以優(yōu)質糧油、蔬菜、水果、花木、茶葉等特色產(chǎn)業(yè)為支撐的多個農業(yè)功能區(qū)［15-16］，農業(yè)產(chǎn)業(yè)結構呈現(xiàn)出產(chǎn)業(yè)空間聚集化、結構多元化的新特征，進一步改變郊區(qū)農業(yè)鄉(xiāng)鎮(zhèn)的土地利用格局。因此，對成都郊區(qū)農業(yè)鄉(xiāng)鎮(zhèn)土地利用變化特征進行研究，使用高分辨率遙感調查方法，獲取成都農業(yè)主導的鄉(xiāng)鎮(zhèn)土地利用空間信息，分析農業(yè)土地利用變化特征，為該地區(qū)農業(yè)產(chǎn)業(yè)綠色發(fā)展、結構優(yōu)化調整、鄉(xiāng)村振興等提供參考。

1 研究數(shù)據(jù)與方法

1.1 研究區(qū)概況

文章選擇原牟禮鎮(zhèn)為研究區(qū)，該區(qū)位于邛崍市東南面，距邛崍市城區(qū)25 km、成都市中心約60 km，交通便利，是典型的城市郊區(qū)鄉(xiāng)鎮(zhèn)，研究區(qū)位置見圖1。該區(qū)域屬亞熱帶濕潤氣候，四季分明，平均氣溫17.8℃，年均降雨量約1300 mm；區(qū)域內地勢平坦，平壩、淺丘相間，以平壩為主。該區(qū)耕地面積約23.33 km2，田間道路布局合理，排灌設施配套，隨著都市農業(yè)發(fā)展，逐步形成了糧油、水果、蔬菜以及水產(chǎn)養(yǎng)殖等多種農業(yè)產(chǎn)業(yè)，是成都郊區(qū)農業(yè)功能主導發(fā)展的鄉(xiāng)鎮(zhèn)。

由圖1 可知，研究區(qū)地形類型包括平壩區(qū)和淺丘，其中，西部為淺丘區(qū)，東部為河流沖擊平壩。調查顯示該區(qū)域在不同地形類型上發(fā)展不同的農業(yè)主導產(chǎn)業(yè)，因此在平壩、淺丘區(qū)域分別隨機建立2 個調查樣區(qū)，樣區(qū)范圍以道路進行劃分，總覆蓋面積約10 km2，具體樣區(qū)對應覆蓋區(qū)域見圖1。

圖1 研究區(qū)位置及調查樣方Fig.1 The location of study area and investigation samples area

1.2 數(shù)據(jù)來源

數(shù)據(jù)來源包括高分辨率遙感影像和地面調查數(shù)據(jù)。收集研究區(qū)分辨率優(yōu)于1 m 的高分辨率遙感影像，通過Google Earth 平臺獲取到2013 年1 月24 日RGB 合成影像，通過購買數(shù)據(jù)獲取2018 年2 月28 日WorldView-2 影像。使用2 期影像進行遙感調查，獲取調查區(qū)內各類土地利用類型的空間信息。2018 年2 月5 日，在研究區(qū)內開展外業(yè)調查，采集典型地物解譯標志，調查樣方內土地利用現(xiàn)狀。

1.3 研究方法

在ENVI 5.5 中對收集的WorldView-2 遙感影像進行數(shù)據(jù)融合、幾何配準等預處理，得到調查樣區(qū)所需的空間分辨率為0.5 m 的遙感影像。結合地面調查數(shù)據(jù)，調查樣區(qū)內土地利用類型包括農業(yè)用地和建設用地，其中，農業(yè)用地包括耕地（糧油作物、蔬菜和其他作物等）、林地、園地（以獼猴桃、葡萄、李子、柑橘等為主的果園用地）、設施大棚、坑塘水面（除天然河流外的人工坑塘、水產(chǎn)養(yǎng)殖以及灌溉蓄水池等），建設用地為建筑物和主要道路。分別對兩個時期的遙感影像進行目視解譯，獲取4 個調查樣區(qū)內各土地利用類型空間信息。

根據(jù)2013 年和2018 年兩期遙感解譯結果，通過地理空間信息技術和土地利用轉移矩陣方法［17］分析該地區(qū)近幾年各類土地利用類型之間的轉移變化情況；將兩期土地利用類型數(shù)據(jù)按淺丘、平壩進行統(tǒng)計，分析不同地形的土地利用變化特征。由此得出該農業(yè)鄉(xiāng)鎮(zhèn)的土地利用變化以及轉移情況。

土地利用轉移矩陣是對土地利用狀態(tài)以及各土地利用類型之間轉移情況的定量描述，研究采用百分比替代面積，直觀反映該地區(qū)期初、期末土地利用類型結構和研究期內轉移變化情況，見式（1）。矩陣的行表示研究期初土地利用類型，列表示研究期末土地利用類型；Sij表示研究期內第i類轉化為第j類的面積占調查土地總面積的百分比；對角線上Sii表示研究期內第i類保持不變的面積所占百分比。

式（1）中，n為土地利用類型的數(shù)量。根據(jù)Sij，按照式（2）計算各土地利用類型轉出（Oi）和轉入（Ij）的總面積所占百分比。

2 結果與分析

2.1 土地利用類型轉移及分析

根據(jù)圖2 兩個時期遙感解譯結果可知，2013—2018 年研究區(qū)土地利用發(fā)生明顯變化，表1 顯示了2013—2018 年土地利用類型轉移情況，樣區(qū)內土地利用總體結構以農業(yè)用地為主，面積比例均超過83%；農業(yè)用地結構發(fā)生明顯變化。2013 年耕地、園地、設施大棚、坑塘水面和林地在農業(yè)用地的占比分別為70.2%、22.0%、0.5%、1.8%和5.6%，到2018 年變?yōu)?5.4%、43.1%、4.7%、3.6%和3.0%；園地、設施大棚和坑塘水面用地占比上升，耕地和林地用地占比下降。該地區(qū)農業(yè)產(chǎn)業(yè)結構調整，糧油、水果、蔬菜產(chǎn)業(yè)為當?shù)剞r業(yè)發(fā)展的主導產(chǎn)業(yè)，農業(yè)用地結構相應發(fā)生變化。

圖2 不同時期的土地利用類型遙感解譯結果Fig.2 Land use types of remote sensing interpreted results in 2013 and 2018

2013—2018 年樣區(qū)內各土地利用類型均發(fā)生了轉入和轉出，其中最為突出的是耕地轉為園地，占調查總面積的15.9%。從土地利用類型轉出來看，耕地轉出24.8%，主要轉為園地、設施大棚和建設用地，分別占耕地轉出面積的64.1%、15.3%和12.5%；林地主要轉為園地，占轉出面積的64.5%；建設用地主要轉為耕地，占轉出面積的47.2%；其余地類轉出較少。從土地利用類型轉入來看，園地轉入19.1%，主要來源于耕地和林地兩大類，分別占轉入面積的83.2%和11.5%；建設用地轉入4.1%，主要來源于耕地，占轉入面積的75.5%；設施大棚轉入3.9%，主要來源于耕地，占轉入面積的97.4%；耕地轉入3.5%，主要來源于園地和建設用地，分別占轉入面積的37.1%和40.0%。根據(jù)土地利用轉移情況，該地區(qū)土地利用類型之間發(fā)生明顯的相互轉換，主要表現(xiàn)為部分耕地被園地、設施大棚、建設用地占用，部分林地變?yōu)閳@地，部分居民點、建設用地和園地又改為耕地；對比建設用地轉移和變化情況，該地區(qū)農業(yè)用地的類型結構變化較大，表明影響農業(yè)鄉(xiāng)鎮(zhèn)地區(qū)土地利用變化的主要因素是農業(yè)產(chǎn)業(yè)結構變化。

根據(jù)走訪調查，除了城市化水平提高、農村勞動力數(shù)量和結構變化、農業(yè)產(chǎn)業(yè)結構調整等因素之外，認為影響該地區(qū)土地利用發(fā)生變化的原因主要包括：（1）政府引導，該地區(qū)充分發(fā)揮資源優(yōu)勢，積極發(fā)展都市現(xiàn)代農業(yè)、優(yōu)化產(chǎn)業(yè)布局、推進土地流轉和適度規(guī)模經(jīng)營，農業(yè)產(chǎn)業(yè)結構由傳統(tǒng)單一化逐漸演變多元化，產(chǎn)業(yè)規(guī)?；厔菁訌?，改變了該地區(qū)農業(yè)土地利用格局。（2）郊區(qū)農業(yè)鄉(xiāng)鎮(zhèn)涉農投入增加，隨著城市化發(fā)展，居民生活水平提高、農產(chǎn)品消費結構改變［18］，為滿足市場需求，服務城市市民休閑，該地區(qū)設施農業(yè)、觀光休閑農業(yè)投入增加，進一步影響該地區(qū)土地利用。

表1 2013—2018 年各土地利用類型轉移Table 1 Transition matrix of the land use types in study area from 2013 to 2018 %

2.2 不同地形類型土地利用變化分析

表2 顯示了不同地形條件下土地利用及變化情況，研究區(qū)各土地利用類型在不同時期、不同地形下所占比例不同，其變化趨勢在空間上表現(xiàn)出明顯差異。從空間分布來看，2013—2018 年耕地和園地在平壩區(qū)和淺丘區(qū)均有不同程度變化，設施大棚、坑塘水面、林地和建設用地僅在一個區(qū)域表現(xiàn)出明顯的變化；從變化量來看，耕地、園地和林地是淺丘區(qū)大于平壩區(qū)，設施大棚則相反。最明顯土地利用轉換在不同地形所占比例結果（表3），進一步顯示了主要發(fā)生變化的土地利用類型之間轉移方向在空間上的差異。淺丘區(qū)以耕地和林地轉為園地為主，分別占調查區(qū)總面積的12.8%和2.2%；平壩區(qū)以耕地轉為設施大棚和果園為主，分別占3.4%和3.1%。該空間變化特征表明：平壩區(qū)果園分布較為零散、有增長趨勢，但高度集約化的蔬菜種植更具規(guī)模優(yōu)勢，同時，該區(qū)域設施農業(yè)投入增加，蔬菜產(chǎn)業(yè)不斷發(fā)展，形成優(yōu)勢產(chǎn)業(yè)；淺丘區(qū)果園擴張明顯，趨向于規(guī)?；?、產(chǎn)業(yè)化發(fā)展，水果種植逐漸成為該區(qū)域的農業(yè)主導產(chǎn)業(yè)。

表2 2013—2018 年不同地形類型土地利用及變化Table 2 Land use changes in different terrain during from 2013 to 2018%

表3 不同地形類型中最明顯土地利用類型轉換Table 3 Results of the most obvious conversion between land use types in different terrain

3 結論與討論

以2013 年1 月24 日的Google Earth 和2018 年2 月28 日WorldView-2 高分辨率遙感影像為數(shù)據(jù)源，結合地面調查，在成都郊區(qū)農業(yè)鄉(xiāng)鎮(zhèn)以遙感抽樣調查方式隨機建立平壩、淺丘共4 個樣區(qū)，通過人工目視解譯，獲取調查樣區(qū)兩個時期的土地利用類型空間信息，采用GIS 技術和轉移矩陣分析其變化特征，得出以下結論：（1）該地區(qū)農業(yè)用地結構發(fā)生明顯變化，農業(yè)產(chǎn)業(yè)類型發(fā)展是影響土地利用變化的主要因素；（2）該地區(qū)農業(yè)土地利用變化特征明顯，具體表現(xiàn)耕地轉為果園、設施大棚、坑塘水面和建設用地，少量的居民點、建設用地和園地轉為耕地；（3）淺丘、平壩不同地形農區(qū)因發(fā)展的主導產(chǎn)業(yè)不同，農業(yè)土地利用變化特征存在差異。其中，淺丘區(qū)以耕地和林地轉為園地為主；平壩區(qū)以耕地變?yōu)樵O施大棚和果園為主。

成都郊區(qū)農業(yè)產(chǎn)業(yè)結構發(fā)生變化，由早期單一的糧油功能主導，演變?yōu)榧Z油、蔬菜、水果、特色作物、養(yǎng)殖等多元結構，在一定程度上保證成都城市瓜果、蔬菜、肉類等農產(chǎn)品供應。近年來，成都市郊區(qū)農業(yè)在土地流轉、適度規(guī)模經(jīng)營的推動下，農業(yè)產(chǎn)業(yè)規(guī)模化、聚集化趨勢加強，形成新的、以蔬菜、水果、花木、茶葉等為主的勞動密集型的郊區(qū)現(xiàn)代農業(yè)產(chǎn)業(yè)。密集型農業(yè)因高度集約耕種和大量施用化肥農藥對土地利用及干預上升，從而引發(fā)土壤退化、環(huán)境污染等一系列農業(yè)生態(tài)問題［19-20］。為保障城市郊區(qū)農業(yè)可持續(xù)發(fā)展，有待進一步通過遙感、大數(shù)據(jù)及GIS 技術進一步分析郊區(qū)農業(yè)土地利用時空格局、監(jiān)測作物種植類型變化，評價農業(yè)土地利用綜合效應。