劉彥文，周 霞，韓冰華，陳 志，陳銳凱

(湖北科技學(xué)院資源環(huán)境科學(xué)與工程學(xué)院，湖北咸寧 437100)

農(nóng)田是糧食作物生產(chǎn)的主要陣地，是一種不可替代的不可再生的自然資源[1]，農(nóng)田利用的轉(zhuǎn)變引發(fā)了人們對(duì)糧食安全的擔(dān)憂[2]。國(guó)外很多學(xué)者對(duì)耕地和重要農(nóng)地保護(hù)必要性、政策變化、規(guī)劃制定、保護(hù)方法等方面進(jìn)行了諸多研究[3-7]。國(guó)內(nèi)農(nóng)地保護(hù)意識(shí)的覺醒始于1978年，但直到1986年耕地作為一項(xiàng)基本國(guó)策被第一次明確提出后才開始了真正意義上的立法保護(hù)[8]。1999年實(shí)施的《基本農(nóng)田保護(hù)條例》中明確要求，各級(jí)人民政府在土地利用總體規(guī)劃中必須明確基本農(nóng)田保護(hù)的布局安排、數(shù)量指標(biāo)和質(zhì)量要求。同年實(shí)施的《中華人民共和國(guó)土地管理法》第四章“耕地保護(hù)”第三十一條明確嚴(yán)格控制耕地轉(zhuǎn)為非耕地、非農(nóng)建設(shè)占用的按數(shù)量和質(zhì)量相當(dāng)?shù)脑瓌t進(jìn)行開墾補(bǔ)償。但根據(jù)中國(guó)國(guó)家統(tǒng)計(jì)局(National Bureau of Statistics of China，2010)和文獻(xiàn)[2]公布的數(shù)據(jù)統(tǒng)計(jì)：自1979年至2008年，耕地仍然減少2 282萬hm2，其中缺乏科學(xué)有效的保護(hù)方法是耕地持續(xù)減少的一個(gè)重要原因。為了有效保護(hù)耕地與農(nóng)田，國(guó)內(nèi)有關(guān)學(xué)者結(jié)合土地評(píng)價(jià)理論[9]、LESA綜合評(píng)價(jià)體系[10-11]、小波理論[12-13]、人工免疫系統(tǒng)[14]等，采用遙感與地理信息技術(shù)，從耕地的自然條件、連片性等角度出發(fā)，選擇質(zhì)量等別、坡度坡向、道路通達(dá)度、水利設(shè)施狀況以及與相鄰?fù)恋乩妙愋偷倪m宜性等因素作為評(píng)價(jià)指標(biāo)[15-18]，對(duì)農(nóng)田土地利用變化、遙感動(dòng)態(tài)監(jiān)測(cè)、基本農(nóng)田劃定及布局優(yōu)化、數(shù)據(jù)庫建設(shè)等進(jìn)行了相關(guān)研究[19-21]?；巨r(nóng)田作為耕地的重要組成部分，如何選擇合適的耕地作為基本農(nóng)田涉及到諸多因素，已有文獻(xiàn)多注重劃定方法的理論研究，筆者從工程化應(yīng)用角度出發(fā)，探討了RS與GIS技術(shù)在基本農(nóng)田布局優(yōu)化中的應(yīng)用。

1 材料與方法

1.1研究區(qū)概況嘉魚縣位于湖北省東南部、長(zhǎng)江中游南岸，地理坐標(biāo)為113°39′～114°22′E，29°48′～30°19′N，地處江漢沖積平原，屬副熱帶濕潤(rùn)季風(fēng)氣候，日照充足，雨量充沛，土地肥沃，全縣轄8個(gè)鄉(xiāng)鎮(zhèn)1個(gè)農(nóng)場(chǎng)，2013年總?cè)丝?7.07萬，是長(zhǎng)江經(jīng)濟(jì)帶重要開發(fā)區(qū)之一，也是國(guó)家食品安全示范縣和全國(guó)首屆國(guó)土資源節(jié)約集約模范縣。地勢(shì)自西南丘崗向東北平原交迭過渡，總體海拔較低，大部分在19～50 m，全縣大體形成“一山三水四分田，兩分道路與莊園”的地貌格局。近年來，嘉魚縣嘉魚農(nóng)業(yè)種植發(fā)生了很大的改變，農(nóng)民由單一種糧邁向糧菜瓜果交織、多季多種種植，其中蔬菜和瓜果的種植與經(jīng)營(yíng)形成規(guī)模，成為特色，享譽(yù)全國(guó)。

1.2數(shù)據(jù)來源數(shù)據(jù)來源主要包括遙感數(shù)據(jù)平臺(tái)下載、國(guó)土部門提供和圖書館查閱等，收集的資料主要包括研究區(qū)社會(huì)經(jīng)濟(jì)數(shù)據(jù)、自然環(huán)境數(shù)據(jù)、土地利用數(shù)據(jù)、衛(wèi)星遙感數(shù)據(jù)以及Google Earth影像合成加工數(shù)據(jù)。獲取方式包括向有關(guān)部門搜集和購(gòu)買、通過文獻(xiàn)綜合研究等手段收集和整理，以及通過實(shí)地考察與調(diào)研進(jìn)行資料分析判斷與綜合取舍。通過與已有相關(guān)數(shù)據(jù)、互聯(lián)網(wǎng)收集資料等進(jìn)行交叉驗(yàn)證，以檢驗(yàn)數(shù)據(jù)來源的真實(shí)可靠性。

1.3研究方法收集原有的相關(guān)遙感影像資料和現(xiàn)成地類資料成果，以遙感技術(shù)為主要手段，利用多時(shí)相的衛(wèi)星遙感資料，根據(jù)地類的可解譯程度，對(duì)工作區(qū)的土地利用類型進(jìn)行監(jiān)督與非監(jiān)督分類解譯。在DEM上提取坡度、坡向、流域水系等相關(guān)數(shù)據(jù)，對(duì)解譯出的地類與上述信息進(jìn)行疊置分析，初步判斷最適宜的農(nóng)田布局。對(duì)土地利用總體規(guī)劃庫、土地利用現(xiàn)狀數(shù)據(jù)庫、基本農(nóng)田基期數(shù)據(jù)庫綜合分析提取的基本農(nóng)田邊界，將此邊界與最適宜耕地布局進(jìn)行空間疊置分析，人工目視劃定最終基本農(nóng)田分布(圖1)。

圖1 技術(shù)路線Fig.1 Technical route

2 結(jié)果與分析

2.1遙感數(shù)據(jù)處理結(jié)果對(duì)Landsat OLI多光譜數(shù)據(jù)進(jìn)行圖像增強(qiáng)、幾何校正、矩形裁剪，采用5×5濾波處理消除微地形、地物陰影、地貌差異以及成像過程等在遙感影像上造成的噪聲，再進(jìn)行輻射校正與大氣校正，大氣校正參數(shù)設(shè)置為single scale factor：1，研究區(qū)平均高程為0.028 km，從ENVI自帶OLI光譜庫中制作多光譜對(duì)應(yīng)波譜曲線庫LC8MultilSli.sli。依據(jù)《土地利用現(xiàn)狀分類》標(biāo)準(zhǔn)，將土地利用類型分為草地、林地、耕地、建設(shè)用地、水域和其他6種類型，分類結(jié)果見圖2A，表1是對(duì)分類信息的統(tǒng)計(jì)，由地理空間云獲取研究區(qū)DEM數(shù)據(jù)，將研究區(qū)坡度分為≤2°、>2°～6°、>6°～15°、>15°～25°、>25° 5個(gè)坡度級(jí)，生成研究區(qū)坡度見圖2B，對(duì)應(yīng)坡度統(tǒng)計(jì)見表2。

注：A.OLI數(shù)據(jù)解譯結(jié)果；B.坡度分布Note：A.OLI data interpretation results；B.Gradient distribution圖2 遙感數(shù)據(jù)處理結(jié)果Fig.2 Remote sensing data processing results

土地類型Soil type編碼Code像元數(shù)The number of pixel像元累積The cumulative pixel百分比Percent∥%累積百分比Cumulative percent∥%草地Grass191 28091 2808.048.04耕地Arable land2249 054340 33421.9429.98建設(shè)用地Construction land3104 605444 9399.2239.20林地Woodland4373 632818 57132.9272.12水域Water area5237 6471 056 21820.9493.06其他Others678 7921 135 0106.94100

表2 坡度統(tǒng)計(jì)結(jié)果

由表1可知，在地類構(gòu)成中，林地占比為32.92%，比重最大，由圖2解譯結(jié)果可知主要分布在縣域西南角，縣域西南角在坡度圖上起伏明顯，該區(qū)是縣域內(nèi)主要的丘陵分布區(qū)，常年有大量植被覆蓋，耕地與水域分別排在第二位和第三位，占比分別為21.94%和20.94%，耕地成連片狀，自縣域中部延伸分布至東北角，該區(qū)地勢(shì)低平，平均坡度在6°以下，是全縣主要的糧食作物和蔬菜主產(chǎn)區(qū)。長(zhǎng)江水域自縣域西側(cè)由北向南貫穿始終，斧頭湖、西涼湖緊鄰其右，蜀茶湖、三湖連江水庫、蜜泉湖、大巖湖等主要湖泊水庫星羅分布期間，為便捷灌溉提供了根本保障。

2.2基本農(nóng)田布局優(yōu)化

2.2.1基本農(nóng)田布局優(yōu)化流程。基本農(nóng)田劃定以二調(diào)基本農(nóng)田上圖成果為基礎(chǔ)，套合新一輪土地利用總體規(guī)劃基本農(nóng)田保護(hù)區(qū)范圍，先確定基本農(nóng)田基期邊界，再將基本農(nóng)田基期邊界與遙感影像上解譯坡度信息、土地利用分類信息、Google Earth人工解譯的路網(wǎng)信息等進(jìn)行疊加分析，評(píng)價(jià)基期基本農(nóng)田劃定的合理性。運(yùn)用解譯信息等對(duì)基期劃定結(jié)果進(jìn)行布局優(yōu)化，劃定基本農(nóng)田保護(hù)布局初步方案，通過數(shù)據(jù)庫導(dǎo)出全縣及分鄉(xiāng)鎮(zhèn)匯總面積，將該套數(shù)據(jù)與全縣及各鄉(xiāng)鎮(zhèn)基本農(nóng)田保護(hù)數(shù)量指標(biāo)進(jìn)行對(duì)比，確定基本農(nóng)田補(bǔ)劃數(shù)量；最后根據(jù)外業(yè)底圖調(diào)查結(jié)果，在ArcGIS中確定最終調(diào)入地塊、調(diào)出地塊和繼續(xù)保留為基本農(nóng)田地塊的空間位置、面積、地類、質(zhì)量等級(jí)等現(xiàn)狀信息，優(yōu)化生成最終成果。

2.2.2基本農(nóng)田調(diào)出分析。嘉魚縣調(diào)出的基本農(nóng)田主要優(yōu)先考慮由于自然原因(如旱、澇)不能繼續(xù)作為基本農(nóng)田的地塊、由于地類改變不適宜作為基本農(nóng)田的地塊等多個(gè)方面。例如，現(xiàn)狀用途為坑塘、林地、村莊等的地塊，共調(diào)出地類圖斑51個(gè)。丘陵地區(qū)基本農(nóng)田調(diào)出主要是由于地類改變所致，以陸溪鎮(zhèn)、高鐵嶺鎮(zhèn)最為典型。濱湖地勢(shì)較低的地區(qū)，湖區(qū)水位的季節(jié)性變化使得土地不適宜作為基本農(nóng)田，進(jìn)而將土地利用類型規(guī)劃為水塘；丘陵坡地地勢(shì)較高的地區(qū)，交通水利設(shè)施不健全，且分布零星，進(jìn)而將土地利用類型規(guī)劃為林地。嘉魚縣平原地區(qū)基本農(nóng)田調(diào)出主要是由于建設(shè)用地占用、二調(diào)地類調(diào)繪錯(cuò)誤，以新街鎮(zhèn)、高鐵嶺鎮(zhèn)最為典型。為了保證城鎮(zhèn)后備用地的儲(chǔ)備，促進(jìn)當(dāng)?shù)亟?jīng)濟(jì)的長(zhǎng)遠(yuǎn)發(fā)展，對(duì)部分城鎮(zhèn)周邊的土地做出了置換，保證等質(zhì)等量對(duì)換、增減達(dá)到平衡。嘉魚縣在“各鄉(xiāng)(鎮(zhèn))2010—2020年土地利用總體規(guī)劃”劃定的基本農(nóng)田的基礎(chǔ)上，共調(diào)出基本農(nóng)田23.62 hm2，其中村莊為1.42 hm2，公路用地為9.13 hm2，果園為1.47 hm2，旱地為0.01 hm2，建制鎮(zhèn)為1.47 hm2，坑塘水面為5.95 hm2，其他林地為3.98 hm2，有林地為0.21 hm2。

2.2.3基本農(nóng)田調(diào)入分析。為了保證基本農(nóng)田“數(shù)量不減少、質(zhì)量不降低”，就要調(diào)入一部分質(zhì)量較高的耕地。依據(jù)集中連片、質(zhì)量由高到低基本原則，調(diào)入的基本農(nóng)田全部為質(zhì)量等級(jí)較好的耕地，如現(xiàn)狀為良好的水田、旱地等。基本農(nóng)田調(diào)入?yún)^(qū)優(yōu)先安排在已實(shí)施完成的土地開發(fā)整理項(xiàng)目區(qū)域內(nèi)，數(shù)量不夠時(shí)則安排在交通便利、灌溉設(shè)施較齊備的城鎮(zhèn)擴(kuò)展區(qū)外圍，以保證嘉魚縣基本農(nóng)田從分布區(qū)域到耕作質(zhì)量都更趨于合理和提高。數(shù)量方面在保證調(diào)出調(diào)入平衡的基礎(chǔ)上，部分村鎮(zhèn)調(diào)整后的基本農(nóng)田面積有所增加。全縣在“各鄉(xiāng)(鎮(zhèn))2010—2020年土地利用總體規(guī)劃”劃定的基本農(nóng)田的基礎(chǔ)上，共調(diào)入地類圖斑47個(gè)，調(diào)入基本農(nóng)田139.18 hm2，其中水田為41.56 hm2，水澆地為89.69 hm2，旱地為7.92 hm2。

2.3基本農(nóng)田優(yōu)化面積對(duì)比分析嘉魚縣基本農(nóng)田基期面積為25 696.52 hm2，新一輪規(guī)劃修編下達(dá)的保護(hù)指標(biāo)為25 435.03 hm2，優(yōu)化前基本農(nóng)田劃定面積為25 478.35 hm2。調(diào)出基本農(nóng)田面積為23.62 hm2，全部從非耕地中調(diào)出；調(diào)入基本農(nóng)田面積為139.18 hm2，一般耕地占調(diào)入面積的比例為100%，因此，基本農(nóng)田上圖面積最終為25 593.90 hm2，略高于本輪規(guī)劃確定的基本農(nóng)田保護(hù)目標(biāo)。

2.4基本農(nóng)田優(yōu)化質(zhì)量對(duì)比分析變動(dòng)地塊質(zhì)量首先要滿足現(xiàn)狀庫質(zhì)量等別統(tǒng)計(jì)要求，但已有相關(guān)成果中劃劣不劃優(yōu)等現(xiàn)象排除不易，該研究運(yùn)用變化地塊同地類面積差值進(jìn)行衡量，減少說明質(zhì)量下降，增加說明質(zhì)量上升。嘉魚縣基本農(nóng)田優(yōu)化過程中，調(diào)出的基本農(nóng)田中無耕地，而調(diào)入的基本農(nóng)田全部為耕地，其中水田占40.45%，水澆地占51.92%，旱地占7.57%，該次基本農(nóng)田優(yōu)化凈調(diào)入耕地139.18 hm2，參照嘉魚縣農(nóng)用地分等技術(shù)報(bào)告，因補(bǔ)劃的基本農(nóng)田都是優(yōu)質(zhì)耕地，優(yōu)化后的基本農(nóng)田平均質(zhì)量等別高于優(yōu)化前的平均質(zhì)量等別(表3)。

表3 基本農(nóng)田優(yōu)化質(zhì)量對(duì)比情況

2.5基本農(nóng)田優(yōu)化地類構(gòu)成分析從基本農(nóng)田構(gòu)成上看，優(yōu)化前基本農(nóng)田面積為25 478.35 hm2，其中耕地面積為25 429.05 hm2，占基本農(nóng)田面積的99.81%(表4)，其他農(nóng)用地面積為47.61 hm2，占基本農(nóng)田面積的0.19%。在耕地中，水田占48.32%、水澆地占39.81%，旱地占11.87%。調(diào)整后基本農(nóng)田面積為25 593.90 hm2，其中耕地面積為25 561.36 hm2，占基本農(nóng)田面積的99.87%，其他農(nóng)用地面積32.54 hm2，占基本農(nóng)田面積的0.13%。在耕地中，水田占48.32%、水澆地占39.81%，旱地占11.87%。調(diào)整后的基本農(nóng)田中耕地所占的比重有所提高，該次基本農(nóng)田調(diào)整中，無基本農(nóng)田調(diào)出，除了官橋鎮(zhèn)、魚岳鎮(zhèn)、簰洲灣鎮(zhèn)無基本農(nóng)田調(diào)入外，嘉魚縣其他5個(gè)鎮(zhèn)均有不同規(guī)模的耕地調(diào)入，尤其以渡普鎮(zhèn)和潘家灣鎮(zhèn)調(diào)入的基本農(nóng)田面積最多，分別調(diào)入為20.87和92.93 hm2。

表4基本農(nóng)田布局優(yōu)化地類構(gòu)成變化分析

Table4Analysisonchangeofcompositionofbasicfarmlandlayoutoptimization

地類Land type占比Proportion∥%優(yōu)化前Before optimization優(yōu)化后After optimization耕地Arable land99.8199.87園地Garden00其他農(nóng)用地Other agricultural land0.190.13

2.6基本農(nóng)田調(diào)布局優(yōu)化集中連片性分析優(yōu)化前基本農(nóng)田布局比較分散，主要分布在嘉魚縣東南部地勢(shì)較平緩、水利設(shè)施完善的平原地區(qū)和西北部地區(qū)，在全縣西南部地區(qū)和中南部地區(qū)，因自然條件和灌溉條件差、土地利用率低，基本農(nóng)田分布較少，集中連片程度總體不高。優(yōu)化后刪除了面積較小的零星像元213個(gè)(多分布在丘陵區(qū))，基本農(nóng)田布局得到改善，補(bǔ)充的像元主要分布在魚岳鎮(zhèn)以北的地區(qū)，基本農(nóng)田總體集中連片程度有所提高(圖3)。

圖3 基本農(nóng)田分布Fig.3 Basic farmland distribution

3 結(jié)論與討論

在RS與GIS支持下，利用土地利用現(xiàn)狀數(shù)據(jù)庫、土地利用總體規(guī)劃數(shù)據(jù)庫和基本農(nóng)田劃定專題數(shù)據(jù)庫中提取的基本農(nóng)田保護(hù)邊界進(jìn)行空間疊加分析，以工程化應(yīng)用為基礎(chǔ)，結(jié)合自然生態(tài)特點(diǎn)與人文社會(huì)發(fā)展，進(jìn)行了嘉魚縣基本農(nóng)田劃定與布局優(yōu)化，優(yōu)化后研究區(qū)共劃定基本農(nóng)田保護(hù)區(qū)1 748個(gè)，基本農(nóng)田保護(hù)片(塊)2 281個(gè)，保護(hù)圖斑11 617個(gè)，面積25 593.90 hm2，經(jīng)現(xiàn)場(chǎng)調(diào)查核實(shí)，所劃成果與實(shí)地現(xiàn)狀一致吻合。

我國(guó)目前對(duì)基本農(nóng)田的劃定，主要是上級(jí)對(duì)下級(jí)按照一定的比例將基本農(nóng)田指標(biāo)分解到各個(gè)省、市、縣、鄉(xiāng)鎮(zhèn)，由于經(jīng)費(fèi)有限、時(shí)間緊迫等主客觀方面的原因，為了完成上級(jí)下達(dá)的目標(biāo)并按時(shí)完成劃定任務(wù)，基本農(nóng)田劃定過程中的實(shí)地調(diào)研等環(huán)節(jié)往往強(qiáng)度不足，具體劃定時(shí)未經(jīng)過全面細(xì)致的調(diào)查思考就將現(xiàn)有耕地中的絕大部分劃為基本農(nóng)田。這可能

存在確實(shí)把較為肥沃作物產(chǎn)量較高的土地劃為基本農(nóng)田，但也存在將類似丘陵、易旱易澇等不適合作為基本農(nóng)田的地塊劃入其中的情況，用統(tǒng)一的基本農(nóng)田劃定數(shù)量來要求全部的地區(qū)是不科學(xué)的，在一定程度上也致使基本農(nóng)田的保護(hù)往往有名無實(shí)。具體劃定過程中人為因素影響較大，基本農(nóng)田首先是耕地，但耕地不可能全部是基本農(nóng)田，決策層一般給出發(fā)展戰(zhàn)略導(dǎo)向，具體劃定時(shí)由于籠統(tǒng)的戰(zhàn)略導(dǎo)向缺乏工程化具體步驟支持，在保證數(shù)量的前提下，空間定位選擇隨意性較大。草率劃定使區(qū)域發(fā)展短期內(nèi)的部分利益可以有效滿足，但對(duì)區(qū)域的可持續(xù)發(fā)展不利。深入挖掘遙感數(shù)據(jù)的潛在信息，利用時(shí)序遙感數(shù)據(jù)蘊(yùn)藏的物候信息、NDVI及其他植被指數(shù)與土地利用的關(guān)聯(lián)信息等，制定土地質(zhì)量客觀性衡量指標(biāo)與工程性步驟化優(yōu)化方法，為后續(xù)提升基本農(nóng)田保護(hù)質(zhì)量、建設(shè)高層級(jí)的高標(biāo)準(zhǔn)農(nóng)田區(qū)提供參考。

[1] CHANG K，YING Y H.External benefits of preserving agricultural land：Taiwan's rice fields[J].The social science journal，2005，42(2)：285-293.

[2] LICHTENBERG E，DING C G.Assessing farmland protection policy in China[J].Land use policy，2008，25(1)：59-68.

[3] BRYANT C R，RUSSWURM L H.North American farmland protection strategies in retrospect[J].Geojournal，1982，6(6)：501-511.

[4] PIERCE J T，F(xiàn)URUSETH O J.Farmland protection planning in British Columbia[J].Geojournal，1982，6(6)：555-560.

[5] FAZAL S.The need for preserving farmland：A case study from a predominantly agrarian economy(India)[J].Landscape and urban planning，2001，55(1)：1-13.

[6] STOBBE T E，EAGLE A J，COTTELEER G，et al.Farmland preservation verdicts-rezoning agricultural land in British Columbia[J].Canadian journal of agricultural economics，2011，59(4)：555-572.

[7] MARTELLOZZO F，RAMANKUTTY N，HALL R J，et al.Urbanization and the loss of prime farmland：A case study in the Calgary-Edmonton corridor of Alberta[J].Regional environmental change，2015，15(5)：881-893.

[8] CHENG L，JIANG P H，CHEN W，et al.Farmland protection policies and rapid urbanization in China：A case study for Changzhou City[J].Land use policy，2015，48：552-566.

[9] 聶艷，吳西子，于婧，等.基于土地評(píng)價(jià)和空間聚類的基本農(nóng)田劃定方法研究：以湖北省鶴峰縣為例[J].中國(guó)土地科學(xué)，2013，27(12)：39-45.

[10] 錢鳳魁，王秋兵.基于農(nóng)用地分等和LESA方法的基本農(nóng)田劃定[J].水土保持研究，2011，18(3)：251-255.

[11] 邊振興，楊子?jì)桑X鳳魁，等.基于LESA體系的高標(biāo)準(zhǔn)基本農(nóng)田建設(shè)時(shí)序研究[J].自然資源學(xué)報(bào)，2016(3)：436-446.

[12] 沈明霞，姬長(zhǎng)英，張瑞合.基于小波變換的農(nóng)田景物邊緣檢測(cè)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào)，2001，32(2)：27-29，41.

[13] 羅輝，張鐵中，楊麗.小波變換在農(nóng)田圖像邊緣檢測(cè)中的應(yīng)用[J].農(nóng)機(jī)化研究，2006(2)：157-160.

[14] LIU X P，LI X，TAN Z Z，et al.Zoning farmland protection under spatial constraints by integrating remote sensing，GIS and artificial immune systems[J].International journal of geographical information science，2011，25(11)：1829-1848.

[15] 喬家君，毛磊.基于RS、GIS村域農(nóng)田數(shù)據(jù)庫設(shè)計(jì)研究：以河南省吳溝村為例[J].農(nóng)業(yè)系統(tǒng)科學(xué)與綜合研究，2009，25(3)：312-316.

[16] 盧德彬，涂建軍，華娟，等.GIS技術(shù)在永久性基本農(nóng)田劃定中的應(yīng)用研究[J].農(nóng)機(jī)化研究，2012(4)：65-68.

[17] 涂建軍，盧德彬.基于GIS與耕地質(zhì)量組合評(píng)價(jià)模型劃定基本農(nóng)田整備區(qū)[J].農(nóng)業(yè)工程學(xué)報(bào)，2012，28(2)：234-238.

[18] 董秀茹，尤明英，王秋兵.基于土地評(píng)價(jià)的基本農(nóng)田劃定方法[J].農(nóng)業(yè)工程學(xué)報(bào)，2011，27(4)：336-339.

[19] 唐永航，孫世宏，孫健.基于GIS的縣級(jí)基本農(nóng)田劃區(qū)定界數(shù)據(jù)庫建設(shè)：以鄞州區(qū)為例[J].中國(guó)土地科學(xué)，2013，27(10)：83-87.

[20] 王曉輝，吳多朋，王錚.基于GPS和GIS技術(shù)的基本農(nóng)田動(dòng)態(tài)監(jiān)測(cè)研究[J].城市勘測(cè)，2013(5)：13-15.

[21] 敖為赳，陳一帆，關(guān)濤，等.GF-1衛(wèi)星數(shù)據(jù)在永久基本農(nóng)田非糧化監(jiān)測(cè)中的應(yīng)用[J].安徽農(nóng)業(yè)科學(xué)，2016，44(18)：250-255.