胡 丹，舒曉波，堯 波，曹慶安

(1.江西師范大學 a.地理與環(huán)境學院，b.科學技術學院，南昌 330022;2.江西省核工業(yè)地質(zhì)調(diào)查院，南昌 330038)

糧食安全關系社會穩(wěn)定和國家長治久安。作為人口眾多的發(fā)展中國家，中國的糧食生產(chǎn)和糧食安全問題已被國內(nèi)外所廣泛關注［1-4］。人均糧食占有量是衡量區(qū)域糧食供需平衡的重要指標，能夠反映一個國家或區(qū)域糧食安全程度。一般人均糧食占有量越高，糧食安全水平亦越高［5-6］。學術界已從不同的研究尺度、研究單元和研究方法，圍繞人均糧食與糧食安全、人均糧食的時空變化及其驅動機制等方面對區(qū)域人均糧食展開研究［7-13］。然而，當前研究大多針對全國尺度，省區(qū)內(nèi)人均糧食異質(zhì)性的研究亦相對集中于較發(fā)達地區(qū)，中部欠發(fā)達地區(qū)的研究則較為少見。

樣帶(transect)是沿著一個主導驅動因素在梯度上有規(guī)律變化的線狀研究區(qū)域或具有明顯差異特征的線狀地域類型［14］。樣帶分析有助于揭示區(qū)域內(nèi)部空間差異規(guī)律及其形成機理的地域差異性，已在全球變化與土地利用研究中得到了成功應用［15-17］。有學者指出，樣帶分析方法從自然地理學研究拓展至人文地理學研究領域，從深化格局動態(tài)與機理研究來看，值得進一步推廣和深入探索［14］。

江西省是傳統(tǒng)的農(nóng)業(yè)大省和糧食主產(chǎn)區(qū)。一方面，江西省經(jīng)濟相對落后，屬于典型的中部欠發(fā)達地區(qū)，其農(nóng)業(yè)產(chǎn)業(yè)化、現(xiàn)代化水平與發(fā)達地區(qū)有所差別，有必要豐富此類區(qū)域人均糧食的研究;另一方面，江西省行政邊界與鄱陽湖流域邊界的吻合度高達97.2%，是一個相對完整的自然地理單元［18-19］，這為基于樣帶分析人均糧食的空間異質(zhì)性提供了典型案例。本研究基于縣域單元，綜合運用樣帶、GIS空間分析、地理探測器等多種方法，分析近14年來江西省人均糧食占有量的時空格局演變及其形成機理，旨在為江西省糧食生產(chǎn)宏觀布局與優(yōu)化、促進區(qū)域農(nóng)業(yè)可持續(xù)發(fā)展提供科學依據(jù)，亦為樣帶分析在人文地理學研究或省區(qū)研究做進一步的嘗試和示范。

1 材料與方法

1.1 研究單元與數(shù)據(jù)來源

由于江西省11個設區(qū)市市區(qū)糧食生產(chǎn)比重較小，且數(shù)據(jù)難以獲取，所以，選取非市區(qū)類縣域作為研究單元，共包括江西省80個縣(市)，時間序列為1997—2010年。為保證研究單元的連貫性和一致性，2000年撤銷前的寧岡縣和2010年新設立的共青城市均未列入研究單元之內(nèi)。數(shù)據(jù)資料來自1998—2011年《江西省統(tǒng)計年鑒》，圖件數(shù)據(jù)來自1∶500 000中國基礎地理數(shù)據(jù)。

1.2 縣域人均糧食占有量模型

區(qū)域人均糧食占有量在較大程度上反映了區(qū)域的糧食供需狀況，人均糧食占有量模型為［11］:

式中:Cit和Pit分別表示i縣域在t年份的人均糧食占有量和總人口;Git表示i縣域相應年份近3年的糧食產(chǎn)量均值;Si表示i縣域人均糧食年均增長率。為便于不同年份縣域間的比較，參考劉玉等［10］的結果，基于人均糧食占有量將各縣域劃分為缺糧區(qū)(＜300 kg)、供需平衡區(qū)(300～400 kg)、一般余糧區(qū)(400～600 kg)和重要余糧區(qū)(＞600 kg)。

1.3 空間自相關分析

空間自相關最為常用的指標是Moran’s I，全局Moran’s I可以用來驗證整個研究區(qū)域的空間模式和度量屬性值在整個區(qū)域空間的分布態(tài)勢;局部Moran’s I是用來揭示空間地域單元與其臨近空間單元屬性特征值之間的相似性或相關性，可以用于識別“熱點區(qū)域”以及數(shù)據(jù)的異質(zhì)檢驗［20-21］。利用GeoDa軟件進行空間自相關分析及Z值顯著性檢驗，空間權重矩陣統(tǒng)一采用鄰接標準定義，即縣域單元i和j有公共邊界，空間權重wij取1，否則取0。

1.4 樣帶選擇及地統(tǒng)計趨勢分析

為更好地揭示江西省縣域人均糧食占有量的時空格局及其內(nèi)在機理，結合江西省自然、人文景觀要素，選取了贛北邊境線、浙贛線、京九線和撫河沿線1 km內(nèi)的縣域，組成了2條東西向和2條南北向樣帶(圖1)。地統(tǒng)計趨勢分析圖中的每一根豎棒代表了一個數(shù)據(jù)點的值(高度)和位置，這些點被投影到一個東西向與南北向的正交平面上，通過投影點繪制出一條最佳擬合線，用以模擬特定方向上的變化趨勢［22-23］。樣帶縣域人均糧食占有量的地統(tǒng)計趨勢分析在ArcGIS 9.2中完成。

圖1 樣帶選取Fig.1 The selection of transects

1.5 地理探測器機理分析方法

地理事物的空間分布深受社會經(jīng)濟或自然要素影響，而地理探測器能夠有效診斷出地理事物空間分布的影響因子的決定力大小，已在地方病形成原因的探測中得到應用［14，24］。地理探測器模型為:

式中:PD，U為人均糧食占有量影響因素探測力指標為次一級區(qū)域樣本數(shù);n為整個區(qū)域樣本數(shù);m為次級區(qū)域個數(shù);為整個區(qū)域人均糧食占有量的方差為次一級區(qū)域的方差。假設，模型成立的取值區(qū)間為［0，1］。PD，U=0，表明人均糧食分布呈隨機分布;PD，U值越大，分區(qū)因素對人均糧食格局的影響越大。

2 結果與分析

2.1 江西省縣域人均糧食占有量的時空變化

江西省縣域人均糧食占有量空間集聚特征明顯，主要年份(1997，2005，2010年)全局 Moran’s I分別為0.53，0.51和0.53，均通過1%的顯著性檢驗，顯示出較強的正自相關性。

全省縣域人均糧食占有量在主要年份中均表現(xiàn)出較為明顯的圈層分布特征，且呈加強趨勢(圖2)。其中，重點余糧區(qū)位于贛中平原地區(qū)，一般余糧區(qū)位于次外圍地區(qū)，供需平衡區(qū)與缺糧區(qū)則鑲嵌分布于最外圍(贛南、贛北地區(qū))。數(shù)量上，1997年全省重要余糧區(qū)20個、一般余糧區(qū)29個、供需平衡區(qū)20個、缺糧區(qū)11個，2005年依次為23，21，21，15 個，2010 年依次為 27，20，17，16個，全省大部分縣域為余糧區(qū)。

圖2 江西省縣域人均糧食占有量的時空變化Fig.2 The spatial and temporal variations of county per captain grain possession in Jiangxi Province

為更清晰地反映人均糧食占有量格局的時空演變，借助LISA圖來表示各主要年份縣域屬性的相互關系(圖3)。2000年人均糧食占有量高高集聚與低低集聚縣域總數(shù)為28個、高低集聚與低高集聚縣域總數(shù)為3個，2005 年分別有25，5 個，2010 年分別有20，5 個，說明研究時段內(nèi)人均糧食占有量的空間均質(zhì)性有小幅下降，而異質(zhì)性變化不大。從“熱區(qū)”與“冷區(qū)”的時空變化來看，研究期內(nèi)“熱區(qū)”格局相對穩(wěn)定于贛中地區(qū)，且呈縮小態(tài)勢(2000，2005和2010年人均糧食占有量“熱區(qū)”分別有18，12和9個);“冷區(qū)”格局則相對離散，1997年與2005年集中于贛南、贛東北、贛西北的山地丘陵地區(qū)，而2010年僅集中于贛南、贛東北地區(qū)，研究期內(nèi)“冷區(qū)”范圍變化不大(2000，2005，2010年“冷區(qū)”分別有10，13，11個)。主要年份中均未發(fā)生變化的縣域總數(shù)占60%以上，說明大部分區(qū)域仍然保持原有格局。

圖3 1997，2005，2010年江西省縣域人均糧食占有量的LISA聚集圖Fig.3 The LISA cluster map of county per captain grain possession in Jiangxi Province(1997，2005，2010)

從人均糧食年均增長來看(圖4)，1997—2010年，43個縣域人均糧食占有量有所提高，主要集中于贛中平原和贛北鄱陽湖周邊等農(nóng)業(yè)生產(chǎn)條件較好的地區(qū);37個縣域有所減少，主要集中于贛南、贛西北和贛東北等條件較差的丘陵、山地區(qū)(圖4a)。2005—2010年是江西省糧食生產(chǎn)大發(fā)展時期，人均糧食增長的縣域覆蓋范圍(57個)較1997—2005年(36個)明顯擴大，幾乎覆蓋贛中與贛北地區(qū)(圖4b，4c)。

圖4 江西省縣域人均糧食占有量的年均增速Fig.4 The average annual growth of county per captain grain possession in Jiangxi Province

2.2 各樣帶縣域人均糧食占有量的變化趨勢

為深入揭示縣域人均糧食占有量的時空分異規(guī)律，在選取贛北邊境線、浙贛線、京九線和撫河沿線4個樣帶的基礎上，采取地統(tǒng)計趨勢線方法，繪制出1997，2005和2010年各樣帶人均糧食占有量變化趨勢線(圖5)。結果顯示:(1)贛北邊境樣帶是江西省縣域人均糧食占有量較低的樣帶。1997年分縣人均糧食呈東西高、中部低的“凹”型趨勢，2005年呈西部低、中部高、東部較低的趨勢，2010年則變?yōu)闁|西低、中部高的“凸”型趨勢，主要是修水縣的降低與浮梁、鄱陽縣的增長所致。(2)浙贛線樣帶，橫貫上饒、鷹潭、撫州、南昌、新余、宜春、萍鄉(xiāng)7市(區(qū))，是刻畫江西省東西向地理分異的典型樣帶。1997，2005和2010年樣帶呈東西低、中部高的“凸”型趨勢，且凸出態(tài)勢有所趨緩，主要是由東鄉(xiāng)、余江縣的增長所造成。(3)京九線樣帶，貫穿江西省贛北、贛中、贛南，經(jīng)九江、南昌、吉安和贛州4個市(區(qū))，跨平原、丘陵、山地等多種地貌，與贛江干流高度吻合，是全省地理要素南北差異顯著的典型樣帶。樣帶均呈北部較低、中部高、南部低的趨勢，且各主要年份樣帶趨勢差異不明顯。(4)撫河沿線樣帶，南北貫穿撫州市全境，經(jīng)撫州、南昌2市(區(qū))，對全省樣帶分析在空間上進行補充。與京九線樣帶趨勢相似，撫河樣帶亦均呈北部較低、中部高、南部低的趨勢，但中部凸出態(tài)勢有所減緩。

進一步分析發(fā)現(xiàn)，贛北邊境樣帶的趨勢陡變，表明贛北地區(qū)人均糧食增長范圍較大，已至贛北邊境地區(qū);贛北邊境與浙贛線樣帶的變化趨勢要比京九線與撫河樣帶更為明顯，表明14年來江西省人均糧食占有量的東西向變化較南北向要大;主要年份中，多數(shù)樣帶均呈中間高、兩端相對較低的趨勢，表明贛中宜春、撫州等傳統(tǒng)農(nóng)業(yè)優(yōu)勢地區(qū)人均糧食居高不下以及贛南、贛東北、贛西北等山地農(nóng)業(yè)生產(chǎn)受限，同時驗證了全省人均糧食占有量呈圈層狀分布的空間格局。

2.3 縣域人均糧食占有量空間格局形成機理分析

糧食生產(chǎn)過程是自然再生產(chǎn)和經(jīng)濟再生產(chǎn)的復合統(tǒng)一體，自然環(huán)境的演替、要素投入的調(diào)整、糧食生產(chǎn)經(jīng)濟效益的變化以及區(qū)域人口數(shù)量的變化等均會導致區(qū)域人均糧食地域格局的改變［10，25］。農(nóng)業(yè)自然地理條件是糧食生產(chǎn)布局的基礎，耕地資源豐富、水熱條件優(yōu)越的地區(qū)，人均糧食通常較高;農(nóng)業(yè)機械、化肥和農(nóng)藥等要素投入水平與效率的增加，有助于降低糧食生產(chǎn)成本，提高生產(chǎn)效率;經(jīng)濟發(fā)展一方面為農(nóng)民非農(nóng)就業(yè)增加機會，一方面通過比較經(jīng)濟效益，大大降低農(nóng)戶種糧積極性，耕地的非農(nóng)化與非糧化現(xiàn)象嚴重［26-29］。部分學者利用不同模型綜合多種要素，對區(qū)域人均糧食變化的驅動機制進行了研究［10，25，30］。在借鑒劉玉等［10］研究的基礎上，選取了1997年人均糧食占有量(糧食生產(chǎn)基礎)、人均耕地面積、種植結構(糧食播面比)、復種指數(shù)、灌溉指數(shù)、地均化肥投入、農(nóng)機總動力投入、農(nóng)林牧漁勞動力規(guī)模、人均GDP、農(nóng)民人均收入10項要素，對全省縣域人均糧食空間格局形成機理進行探測研究。

圖5 1997，2005和2010年樣帶人均糧食占有量變化趨勢線Fig.5 The change trend of per captain grain possession in transect(1997，2005，2010)

通過地理探測器模型計算得出江西省及樣帶內(nèi)各要素對縣域人均糧食空間分布影響的PD，U值(表1)。結果表明，江西省縣域人均糧食的空間格局主要受1997年人均糧食占有量與人均耕地面積的影響，其次為農(nóng)機總動力投入、農(nóng)村居民人均純收入、種植結構、復種指數(shù)，而人均GDP、灌溉指數(shù)、農(nóng)林牧漁勞動力規(guī)模、地均化肥投入等要素的影響不大。研究發(fā)現(xiàn)，全省人均糧食對自然地理條件表現(xiàn)出較大的依賴，人均糧食較高的區(qū)域仍主要集中于人均耕地面積量大的傳統(tǒng)農(nóng)區(qū)，而在此基礎上，適當加大農(nóng)機投入、提高復種指數(shù)與糧食播種面積等有助于提高區(qū)域人均糧食占有量。對發(fā)達地區(qū)的研究顯示，人均GDP對糧食生產(chǎn)具有明顯的副作用，而本研究中人均GDP的作用幾乎為零，這主要是江西省社會經(jīng)濟水平相對落后、農(nóng)業(yè)生產(chǎn)還未能擺脫自然地理格局對其發(fā)展的制約所致［10，28］。4條樣帶縣域人均糧食地理探測器結果與全省大體一致，均主要受1997年人均糧食占有量與人均耕地面積的影響，但相互間亦存在部分差異。贛北邊境樣帶與撫河樣帶縣域具有一定的農(nóng)業(yè)生產(chǎn)自然條件，隨著農(nóng)用物資投入、復種指數(shù)與糧食播面比的提高，糧食生產(chǎn)會有較大幅度的提升，兩條樣帶內(nèi)種植結構、復種指數(shù)、地均化肥投入和農(nóng)林牧漁勞動力規(guī)模的PD，U值均明顯高于全省。京九線樣帶貫穿江西南北，最能反映全省自然與人文景觀的差異，其1997年人均糧食占有量與人均耕地面積的PD，U值是全省的近1.5倍，更加突出了江西省人均糧食較大依賴于區(qū)域自然地理條件的特點。浙贛線樣帶經(jīng)過全省7個市(區(qū))，是經(jīng)濟發(fā)展水平較高的一條樣帶，人均GDP的PD，U值顯著高于全省及其他樣帶。

此外，糧食政策、市場價格等對區(qū)域糧食生產(chǎn)亦具有一定影響。2004年以來，國家不斷強調(diào)農(nóng)業(yè)生產(chǎn)的戰(zhàn)略問題，陸續(xù)出臺糧食直補、增加農(nóng)業(yè)財政支出、取消農(nóng)業(yè)稅等各項惠農(nóng)政策，放開糧食購銷市場，加強和改善糧食宏觀調(diào)控。陳飛等［31］研究發(fā)現(xiàn)各項惠農(nóng)政策對糧食生產(chǎn)均具有顯著正向影響，其中，農(nóng)業(yè)支出政策和農(nóng)村固定資產(chǎn)投資是拉動我國糧食產(chǎn)量增長的最重要因素。劉克春［32］基于江西省農(nóng)戶調(diào)查，認為政府現(xiàn)行的糧食直接補貼、最低收購價政策，調(diào)節(jié)了以糧食生產(chǎn)為主要收入來源的農(nóng)戶的糧食種植決策行為，提高了其糧食生產(chǎn)積極性，一定程度上促進了糧食種植面積的擴大。市場價格直接影響糧農(nóng)收入預期以決策其糧食種植行為，糧價低迷或糧價相對于其他農(nóng)產(chǎn)品和工業(yè)品價格低，均會嚴重降低糧農(nóng)積極性，而農(nóng)業(yè)生產(chǎn)資料價格走高亦會負向調(diào)節(jié)農(nóng)戶的糧食種植決策行為?？梢灶A見，在政府資金、政策傾斜及日益完善的市場經(jīng)濟等條件下，區(qū)域糧食生產(chǎn)地域功能將得到強化，區(qū)域人均糧食占有量將進一步提高。

表1 江西省及不同樣帶各要素地理探測器分析的PD，U值Tab.1 PD，Uof geographic detector analysis between per captain grain possession and different elements in Jiangxi Province

3 結論與討論

(1)江西省縣域人均糧食占有量空間集聚特征明顯，顯示出較強的正自相關性。全省縣域人均糧食占有量在主要年份中均表現(xiàn)出較為明顯的圈層分布特征，且呈加強趨勢。研究期內(nèi)“熱區(qū)”格局相對穩(wěn)定于贛中地區(qū)，且呈縮小態(tài)勢;“冷區(qū)”格局則相對離散，“冷區(qū)”范圍變化不大。全省大部分縣域為余糧區(qū)，過半的縣域人均糧食占有量有所提高，主要集中于贛中平原和贛北鄱陽湖周邊等農(nóng)業(yè)生產(chǎn)條件較好的地區(qū)，2005年來，人均糧食增長的縣域覆蓋范圍明顯擴大。

(2)樣帶分析能夠細膩刻畫全省人均糧食的時空分異規(guī)律。主要年份中，多數(shù)樣帶均呈中間高、兩端相對較低的“凸”型趨勢，表明贛中宜春、撫州等傳統(tǒng)農(nóng)業(yè)優(yōu)勢地區(qū)人均糧食居高不下以及贛南、贛東北、贛西北等山地農(nóng)業(yè)生產(chǎn)受限，同時，驗證了全省人均糧食占有量呈圈層狀分布的空間格局。贛北邊境與浙贛線樣帶的變化趨勢要比京九線與撫河樣帶更為明顯，表明14年來江西省人均糧食占有量的東西向變化較南北向要大。贛北邊境樣帶的趨勢陡變，則表明全省人均糧食增長范圍已擴至贛北邊境地區(qū)。

(3)江西省縣域人均糧食的空間格局主要受1997年人均糧食占有量與人均耕地面積的影響，其次為農(nóng)機總動力投入、農(nóng)村居民人均純收入、種植結構、復種指數(shù);較低的經(jīng)濟水平使得全省人均糧食對自然地理條件表現(xiàn)出較大的依賴。4條樣帶縣域人均糧食地理探測器結果與全省大體一致，均主要受1997年人均糧食占有量與人均耕地面積的影響。不同樣帶地理探測分析的結果有助于深入認識不同農(nóng)業(yè)生產(chǎn)自然條件、不同經(jīng)濟發(fā)展水平等縣域人均糧食驅動機制的差異。此外，糧食政策、市場價格等對區(qū)域糧食生產(chǎn)亦有一定影響。

(4)以樣帶為基礎，通過地統(tǒng)計趨勢研究與地理探測器等方法，揭示了全省及樣帶縣域人均糧食占有量的時空格局及形成機理。樣帶分析能夠較好地應用于區(qū)域人文地理研究，可以進一步推廣至其他區(qū)域。例如，京九線樣帶能很大程度反映出江西全省地理事物的空間分布，而浙贛線樣帶則可以反映出經(jīng)濟發(fā)展較好地區(qū)地理事物的分布特征。地理探測器方法主要基于各要素與變量之間的匹配程度，以量化各要素對變量的空間格局所產(chǎn)生的影響，相比傳統(tǒng)的回歸分析、因子分析等，具有更好的地理空間解釋能力。另外，以縣域為單元，在一定程度上忽略了更小尺度上的相互作用，指標選取上亦受數(shù)據(jù)限制，對自然環(huán)境因素考慮不足。由此，基于多源數(shù)據(jù)(如遙感、實地調(diào)查等)的區(qū)域人均糧食研究將是進一步深化的研究方向。

［1］Brown L.Who Will Feed China?Wake-up Call for A Small Planet［M］.New York:W.W.Norton ＆ Company，1995.

［2］黨安榮，閻守邕，吳宏岐，等.基于GIS的中國土地生產(chǎn)潛力研究［J］.生態(tài)學報，2000，20(6):910-915.

［3］劉彥隨，翟榮新.中國糧食生產(chǎn)時空格局動態(tài)及其優(yōu)化策略探析［J］.地域研究與開發(fā)，2009，28(1):1-5.

［4］封志明，楊艷昭，張晶，等.從柵格到縣域:中國糧食生產(chǎn)的資源潛力區(qū)域差異分析［J］.自然資源學報，2007，22(5):747-755.

［5］殷培紅，方修琦，馬玉玲，等.21世紀初我國糧食供需的新空間格局［J］.自然資源學報，2006，21(4):625-631.

［6］張利國.我國區(qū)域糧食安全演變:1949—2008［J］.經(jīng)濟地理，2011，31(5):833-838.

［7］中國農(nóng)業(yè)科學院.人均400公斤糧食必不可少［J］.中國農(nóng)業(yè)科學，1986，19(5):1-7.

［8］杜忠潮，車自力.陜西省人均糧食變化及糧食安全研究［J］.經(jīng)濟地理，2002，22(6):744-749.

［9］郭柏林.中國人均糧食變化的地理研究［J］.地理科學，1996，16(3):238-244.

［10］劉玉，劉彥隨，郭麗英.環(huán)渤海地區(qū)縣域人均糧食占有量空間格局演化［J］.地理科學，2011，31(1):102-109.

［11］劉玉，高秉博，潘瑜春.中國縣域人均糧食占有量空間變異特征研究［J］.中國農(nóng)業(yè)大學學報，2012，17(3):154-162.

［12］李裕瑞，卞新民.江蘇省糧食生產(chǎn)地域格局變化研究［J］.地域研究與開發(fā)，2008，27(2):113-117.

［13］熊友云，張明軍，劉園園，等.中國糧食產(chǎn)量省區(qū)差距變化及其成因分析［J］.自然資源學報，2009，24(6):965-974.

［14］劉彥隨，楊忍.中國縣域城鎮(zhèn)化的空間特征與形成機理［J］.地理學報，2012，67(8):1011-1020.

［15］Koch G W，Scholes R J，Steffen W L，et al.The IGBP Terrestrial Transects:Science Plan［R］.Stockholm:IGBP，1995.

［16］Long H L，Heilig G K，Li X B，et al.Socio-economic Development and Land-use Change:Analysis of Rural Housing Land Transition in the Transect of the Yangtse River，China［J］.Land Use Policy，2007，24(1):141-153.

［17］龍花樓.中國農(nóng)村宅基地轉型的理論與證實［J］.地理學報，2006，61(10):1093-1100.

［18］安和平，王瑞璞.我國欠發(fā)達地區(qū)農(nóng)業(yè)增長的因素分析［J］.中國農(nóng)村經(jīng)濟，1999(7):62-65.

［19］俞勇軍，陸玉麒.江西省經(jīng)濟區(qū)位特征及其區(qū)域發(fā)展空間結構研究［J］.經(jīng)濟地理，2003，23(4):462-471.

［20］Anselin L.The Local Indicators of Spatial Association-LISA［J］.Geographical Analysis，1995，27(2):93-115.

［21］徐建華.計量地理學［M］.北京:高等教育出版社，2006.

［22］Egenhofer M J，Golledge R G.Spatial and Temporal Reasoning in Geographic Information Systems［M］.New York:Oxford University Press，1998.

［23］Burrough P A，McDonnell R A.Principles of Geographical Information Systems［M］.New York:Oxford University Press，1998.

［24］Wang J F，Li X H，Christakos G，et al.Geographical Detectors-based Health Risk Assessment and Its Application in the Neural Tube Defects Study of the He Shun Region，China［J］.International Journal of Geographical Information Science，2010，24(1):107-127.

［25］郭淑敏，馬帥，陳印軍.我國糧食主產(chǎn)區(qū)糧食生產(chǎn)影響因素研究［J］.農(nóng)業(yè)現(xiàn)代化研究，2007，28(1):83-87.

［26］胡冰川，吳強，周曙東.糧食生產(chǎn)的投入產(chǎn)出影響因素分析——基于江蘇省糧食生產(chǎn)的實證研究［J］.長江流域資源與環(huán)境，2006，15(1):71-75.

［27］劉彥隨.中國東部沿海地區(qū)鄉(xiāng)村轉型發(fā)展與新農(nóng)村建設［J］.地理學報，2007，62(6):563-570.

［28］李裕瑞，呂愛清，卞新民.江蘇省人均糧食地域格局變化特征及驅動機制［J］.資源科學，2008，30(3):423-430.

［29］李裕瑞，劉彥隨，龍花樓.江蘇省糧食生產(chǎn)時空變化的影響機制［J］.地理科學進展，2009，28(1):125-131.

［30］許月卿，李秀彬.河北省糧食生產(chǎn)灰色關聯(lián)動態(tài)分析［J］.地理研究，2002，21(3):339-346.

［31］陳飛，范慶泉，高鐵梅.農(nóng)業(yè)政策、糧食產(chǎn)量與糧食生產(chǎn)調(diào)整能力［J］.經(jīng)濟研究，2010(11):101-114.

［32］劉克春.糧食生產(chǎn)補貼政策對農(nóng)戶糧食種植決策行為的影響與作用機理分析——以江西省為例［J］.中國農(nóng)村經(jīng)濟，2010(2):12-21.