舒瑩，胡遠滿

摘要：取1986、1996和2001年3個時期覆蓋近代黃河三角洲的TM影像為數(shù)據(jù)源，結(jié)合1980年地形圖，經(jīng)監(jiān)督、非監(jiān)督分類及人工目視解譯，獲得研究區(qū)各時相土地覆被類型圖。運用地形圖和TM影像人工目視解譯勾畫3個時期居民點圖，通過3個時期的居民點分布圖向外等距離擴展形成輻射狀梯度圈，并以此剪裁相應(yīng)時相的土地覆被類型圖，形成不同距離輻射梯度土地覆被類型圖，計算各土地覆被類型圖的景觀指數(shù)，分析研究區(qū)15年間土地覆被動態(tài)過程、居民點動態(tài)過程以及居民點變化對景觀動態(tài)的影響。結(jié)果表明，隨著居民點建成區(qū)面積的增大，其對周圍地區(qū)影響的范圍亦不斷增大。

關(guān)鍵詞：近代黃河三角洲；居民點；土地覆被；緩沖區(qū)；梯度圈

中圖分類號：F301.2 文獻標識碼：A 文章編號：0439-8114（2013）05-1030-06

Temporal-Spatial Effects of Build-up Residential Area in Yellow River Delta on Land Use Pattern

SHU Ying1，HU Yuan-man2

（1.School of Environment and Energy Engineering， Anhui Institute of Architecture and Industry， Hefei 230023， China；

2.Institute of Applied Ecology， Chinese Academy of Sciences， Shenyang 110016，China）

Abstract： Using TM remote sensing images of yellow river delta in 1986， 1996 and 2001 and relief maps in 1980 as data sources， land cover landscape of the studied zone in different years were obtained by supervised classification， non-supervised classification and artificial visual interpretation. The map of residential area in the different years was pictured using TM remote sensing images and relief maps. The radiation transects were set around the residential area by buffer zone techniques， according to which the land cover landscapes of different years were cut out to form land cover landscapes of different radiation distance. The landscape indices were calculated based on the land cover landscapes cut out. Based on these data， the dynamics of land cover landscape and residential area， and the effects of residential area change on landscape in the 15-year period were analyzed. The results showed that the influencing zone of residential area was increased with the increasing of its area.

Key words： recent yellow river delta； residential area； land cover； buffer zone； grading diversification

隨著人類活動對土地利用方式日益深刻的影響與干擾，全球變化已成為當今世界一個重大的課題，土地利用與土地覆被變化是全球變化中的重要內(nèi)容之一[1]?；谶b感、GIS等對土地利用與土地覆被變化驅(qū)動力的分析大量涌現(xiàn)，結(jié)果表明土地利用與土地覆被變化的驅(qū)動力集中于自然、人口和社會經(jīng)濟條件幾個方面，自然因素是土地利用變化的內(nèi)在驅(qū)動力，人口增長、經(jīng)濟發(fā)展引起的城鎮(zhèn)建成區(qū)和農(nóng)村居民點用地的擴大是主要的外在驅(qū)動力[2-4]。當前對于居民點建成區(qū)的研究集中于建成區(qū)居民點自身，如建成區(qū)形態(tài)研究、擴展模式、動態(tài)模型研究和空間結(jié)構(gòu)研究等方面[5-7]；農(nóng)村居民點用地的變化、增長、驅(qū)動力分析及各個居民點的空間分布特征分析等方面[8，9]，但關(guān)于居民點建成區(qū)擴張對土地利用與土地覆被格局的影響研究甚少。居民點建成區(qū)從中心向周圍呈梯度擴張[10，11]，周圍土地類型和土地利用格局會隨之呈現(xiàn)景觀梯度分異，這就為利用緩沖區(qū)技術(shù)進行景觀梯度分異研究提供了依據(jù)。

黃河三角洲是中國三大河口三角洲之一，地處海陸生態(tài)交錯帶，自然條件復(fù)雜，一直是學(xué)術(shù)研究的熱點，眾多研究結(jié)果表明黃河三角洲土地利用與土地覆被變化的主要外在驅(qū)動力是人為干擾引起人造景觀增加造成，關(guān)于人為干擾對區(qū)域土地利用格局影響的程度、范圍及對周圍土地斑塊的梯度影響的研究極少。近代黃河三角洲內(nèi)坐落著我國第二大油田——勝利油田和國家級自然保護區(qū)，區(qū)域內(nèi)土地利用類型豐富多樣，有利于研究建成區(qū)居民點對不同景觀類型的影響。本研究以1986、1996和2001年3個時期的TM遙感影像為數(shù)據(jù)源，以區(qū)域土地覆被為載體，研究建成區(qū)居民點緩沖區(qū)土地覆被景觀時空梯度分異。通過研究土地覆被景觀數(shù)據(jù)的幅度效應(yīng)，進一步探討不同面積幅度的梯度圈之間進行景觀特征比較的可行性； 最后，沿居民點建成區(qū)做緩沖區(qū)，形成從居民點建成區(qū)為中心向外圍擴散的輻射帶圈，分析帶圈梯度上的土地覆被景觀的梯度分異特征。通過土地覆被景觀的特征， 映射出居民點建成區(qū)發(fā)展對土地覆被變化的影響。

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

1.1 研究區(qū)概況

黃河三角洲根據(jù)其成陸時間可分為古代黃河三角洲、近代黃河三角洲和現(xiàn)代黃河三角洲。本研究的黃河三角洲區(qū)域主要指近代黃河三角洲（其中包括現(xiàn)代黃河三角洲），總面積60多萬公頃。近代黃河三角洲地處渤海灣畔，地理位置37°20′-38°12′ N，118°07′-119°10′E，以墾利寧海為頂點，北起套爾河口，南至淄脈溝口，境內(nèi)屬暖溫帶大陸性季風氣候。行政區(qū)劃上包括山東省東營市大部分和濱州市東北部地區(qū)，坐落有山東黃河三角洲國家級自然保護區(qū)，是我國最大的新生濕地生態(tài)系統(tǒng)，境內(nèi)有大小河流20余條。

1.2 分析數(shù)據(jù)的生成

原始數(shù)據(jù)以覆蓋整個黃河三角洲的Landsat TM遙感影像數(shù)據(jù)（1986年5月、1996年9月和2001年6月）為本研究基本數(shù)據(jù)源，輔助于該地區(qū)1980年1∶50 000地形圖。以地形圖為基準對遙感影像進行幾何校正，誤差控制在一個像元內(nèi)。再通過人工解譯和計算機自動解譯遙感影像，以及現(xiàn)場踏勘校正解譯結(jié)構(gòu)，形成研究區(qū)土地利用與土地覆被類型圖。其中的分類體系參照我國土地利用分類標準[12]，再考慮黃河三角洲實際情況以及研究內(nèi)容，將近代黃河三角洲分為15種土地類型：黃河、光板灘涂、湖泊水庫、林地、蝦蟹田、鹽田、草甸、蘆葦、農(nóng)田、檉柳蘆葦沼澤、灘涂、潮下帶、道路、油田和居民點。圖像處理使用軟件Erdas 8.6、Arcview 3.3和Arc/info 8.3。

不同尺度下求得的景觀指數(shù)會存在差異，尺度的不同可能直接影響到研究成果的可靠性，尺度一般用粒度和幅度進行分析。選取域值粒度以下的任一粒度大小都不會對研究結(jié)果產(chǎn)生質(zhì)的影響，由于分析數(shù)據(jù)來源于像元大小為30 m的TM影像，為最大限度地包含原始信息量，故數(shù)據(jù)分析粒度采用30 m為基本像元[13-20]。

1.3 研究方法

1.3.1 不同幅度景觀圖的生成 本研究的幅度是指景觀分析區(qū)域面積的大小。相應(yīng)地，幅度效應(yīng)指景觀指數(shù)對不同景觀分析面積大小的響應(yīng)。 根據(jù)研究區(qū)形狀特點，在近代黃河三角洲內(nèi)部截取一個最大的正方形，此為分析的最大區(qū)域（總共包括1 600×1 600個像元）。再以正方形4個頂點為出發(fā)點，分別向內(nèi)裁減50行和50列的帶寬，形成第二個區(qū)域。依次用同樣的方法在第二個區(qū)域的基礎(chǔ)上形成第三個區(qū)域，在第三個區(qū)域上形成第四個區(qū)域，以此類推，形成一系列不同幅度的景觀圖，共包括16個區(qū)域，記為區(qū)域1、2…16，最小的區(qū)域包括100×100個像元（圖1）。

1.3.2 居民點建成區(qū)緩沖區(qū)的生成 居民點建成區(qū)圖層形成過程：根據(jù)1980年1∶50 000地形圖對1986年TM影像進行人工目視勾繪，1996年和2001年在1986年基礎(chǔ)上通過判斷影像紋理和顏色特征進行勾繪。

緩沖區(qū)形成過程：以全部居民點建成區(qū)為對象，向外圍形成1 km寬度的緩沖區(qū)，剪裁建成區(qū)及其緩沖區(qū)，生成第一個梯度圈；再以第一梯度圈為基礎(chǔ)，向外圍形成1 km寬度的緩沖區(qū)，剪裁第一梯度圈及其外圍緩沖區(qū)，生成第二個梯度圈，以此類推，直至覆蓋整個研究區(qū)域。研究區(qū)域西南部靠近內(nèi)陸，居民點建成區(qū)密集，第一個梯度圈即可延伸覆蓋研究區(qū)西南部；北部和東部與渤海毗鄰，其中坐落著黃河三角洲國家級自然保護區(qū)，居民點建成區(qū)分布少，尤其是東部黃河入海處，故在北部和東部需生成較多梯度圈，1986年生成29個梯度圈，1996年生成34個梯度圈，2001年生成30個梯度圈（黃河1996年8月經(jīng)人工改道由原來東部入海改為東北方向入海，所以2001年近代黃河三角洲東部相對1996年有一定蝕退）。基于這些梯度圈剪裁各時相土地覆被景觀圖，形成從居民點建成區(qū)向外圍的一系列輻射的梯度圈景觀，如圖2所示。

1.3.3 景觀指數(shù)選取與計算 景觀格局指數(shù)（Landscape indices or metrics）用簡單數(shù)字描述復(fù)雜的景觀格局。在景觀指數(shù)中，大部分指數(shù)之間景觀格局信息重復(fù)[15，21，22]，指數(shù)之間存在一定的相關(guān)性，故可采用少數(shù)的指數(shù)代替所有指數(shù)進行分析。針對景觀水平常用的景觀指數(shù)，利用景觀格局分析軟件APACK 2.2計算不同幅度景觀圖的指數(shù)值。根據(jù)指數(shù)計算結(jié)果選取其中受幅度影響較小的指數(shù)，并考慮各指數(shù)的生態(tài)學(xué)意義及重要程度，最終選定進行土地覆被景觀的梯度分異指數(shù)。以選取的指數(shù)為介質(zhì)探討居民點建成區(qū)緩沖區(qū)土地覆被景觀時空梯度分異。各指數(shù)的計算公式見該軟件的幫助文件。

2 結(jié)果與分析

黃河三角洲在1986、1996、2001年3個時期總面積先少量增加后大幅減少，1986～1996年研究區(qū)增加土地面積4 006.9 hm2（其中新生灘涂面積5 015.9 hm2，沿海蝕退面積1 009.0 hm2），1996年到2001年土地凈減少16 960.5 hm2（其中新生灘涂面積2 623.0 hm2，沿海蝕退面積19 583.5 hm2）。這是由于近年來黃河流量銳減，上游帶來泥沙日益減少，淤積速度大大小于蝕退速度所造成。

2.1 景觀的幅度效應(yīng)及指數(shù)的選取

幅度是空間數(shù)據(jù)對空間地理現(xiàn)象或過程進行度量和評價的一個重要影響因素，研究表明幅度變化直接影響景觀格局[10，16，22]。選取常用的13個指數(shù)進行幅度效應(yīng)分析，分別為歸一化平均斑塊面積（Average normalized patch area，ANPA）、平均斑塊周長面積比（Average patch perimeter-area ratio，APA）、校正的斑塊周長面積比（Average corrected patch perimeter-area ratio，ACPA）、邊界密度（Edge density，ED）、雙對數(shù)分維數(shù)（Fractal dimension，F(xiàn)D）、多樣性（Shannon-Weaver diversity，SWD）、均一度（Shannon-Weaver evenness，SWE）、優(yōu)勢度（Dominance，DO）、蔓延度（Contagion，CI）、邊界分布均勻度（Edge density evenness，EDE）、角秒矩（Angular second moment，ASM）、反差矩（Inverse difference moment，IDM）和聚集度（Aggregation index，AI）。根據(jù)景觀格局指數(shù)的生態(tài)學(xué)意義（Fragstats的分類標準）[23，24]，以上指數(shù)可分為4個類別，其中歸一化平均斑塊面積、平均斑塊周長面積比、校正的斑塊周長面積比和邊界密度為基本指數(shù)，雙對數(shù)分維數(shù)為形狀指數(shù)，聚集度、蔓延度、角秒矩、反差矩和邊界分布均勻度為蔓延度指數(shù)，優(yōu)勢度、多樣性和均一性屬于多樣性指數(shù)。

2.1.1 景觀指數(shù)的幅度效應(yīng) 幅度效應(yīng)決定著梯度分析結(jié)果的可靠性，在研究區(qū)域景觀特征之前，先以2001年影像為例分析土地覆被景觀幅度效應(yīng)，以保證研究的完整性以及后續(xù)研究成果的可行性與可靠性。對2001年度分析區(qū)域圖進行指數(shù)計算，結(jié)果如圖3。從圖3可以看出，2001年平均斑塊周長面積比、邊界密度、反差矩、雙對數(shù)分維數(shù)和聚集度指數(shù)的幅度效應(yīng)不明顯，其他指數(shù)幅度效應(yīng)較為明顯；從區(qū)域16到區(qū)域12（即包括像元從100×100個到500×500個）各指數(shù)受幅度影響較敏感，區(qū)域12之前指數(shù)值變化趨緩。說明當景觀面積大于區(qū)域12（包括像元500×500個，面積22 500 hm2）并且景觀之間面積相差不大時，運用上述景觀指數(shù)法進行景觀特征分析可行有效。

2.1.2 景觀指數(shù)的選取 景觀指數(shù)之間存在信息重復(fù)，指數(shù)之間存在顯著相關(guān)性?；局笖?shù)中校正斑塊周長面積比不是一個很好的格局指數(shù)[25]，而基本指數(shù)中的歸一化平均斑塊面積與其他指數(shù)不存在顯著相關(guān)性，代表了景觀的獨特性，在景觀對比研究中應(yīng)重視此指數(shù)，解釋其景觀生態(tài)學(xué)意義[23]，故在基本指數(shù)中選取歸一化平均斑塊面積指數(shù)；形狀指數(shù)表示構(gòu)成景觀斑塊形狀的復(fù)雜性，雙對數(shù)分維數(shù)能預(yù)示大多形狀指數(shù)的變化趨勢，是景觀水平上常用的形狀指數(shù)；蔓延度指數(shù)表明景觀類型在空間上的聚集程度或類型間鑲嵌程度，聚集度和反差矩是常用的兩個蔓延度指數(shù)，從圖3中可看出這兩個指數(shù)受幅度影響小，研究結(jié)果可以剔除其幅度的影響；多樣性指數(shù)關(guān)系到景觀類型數(shù)以及類型的面積百分比，包括的大多數(shù)指數(shù)存在明顯的相關(guān)性，其中Shannon-Wiener多樣性指數(shù)多見于景觀格局分析。由此共選擇歸一化平均斑塊面積、雙對數(shù)分維數(shù)、聚集度、反差矩和多樣性指數(shù)這5個指數(shù)分析居民點緩沖區(qū)對區(qū)域土地利用格局的影響。

2.2 研究區(qū)景觀特征時間變化

研究區(qū)景觀特征時間變化如圖4所示，可以看出研究區(qū)15年中，土地覆被景觀變化較為平緩。其中均一化平均斑塊面積指數(shù)經(jīng)歷了先上升后下降趨勢，但變化幅度不大，在0.5左右波動，說明3個年份在斑塊形狀上相似。聚集度和反差矩指數(shù)在研究年份間均出現(xiàn)下降趨勢，其中聚集度指數(shù)值1986年為0.976、1996年為0.973、2001年為0.972，反差矩值1986年為0.950、1996年為0.945、2001年為0.941，說明隨著年份的增加斑塊破碎化程度增加、研究區(qū)域結(jié)構(gòu)趨于簡單。雙對數(shù)分維數(shù)呈現(xiàn)上升趨勢，但上升速度極為緩慢，說明研究年段斑塊邊緣幾何形狀僅略趨于復(fù)雜。進行研究的五個指數(shù)中變化幅度最大的為多樣性指數(shù)，呈明顯上升趨勢，說明研究年段整個研究區(qū)多樣性不斷上升，各種景觀覆被類型所占比例越來越均衡。

2.3 居民點建成區(qū)緩沖區(qū)土地覆被景觀變化的梯度分異

由于黃河三角洲地理位置的特殊性，部分區(qū)域淤積（黃河入?？跒槭澜缟铣赏磷羁斓暮涌谌侵蓿?、部分區(qū)域蝕退（海水沖刷），造成3個年份研究區(qū)面積不同；各年度居民點數(shù)量面積亦不同，故3個年度居民點緩沖區(qū)數(shù)量有所差別。圖5-圖9中1代表居民點第一個緩沖區(qū)形成的第一梯度圈，依次類推，直至覆蓋整個研究區(qū)，即2001年居民點外圍第30緩沖區(qū)形成的第30梯度圈覆蓋整個研究區(qū)，得出的為整個研究區(qū)指數(shù)值，1996年第34梯度圈為整個研究區(qū)指數(shù)值，1986年第29梯度圈為整個研究區(qū)指數(shù)值。根據(jù)景觀指數(shù)幅度分析結(jié)果，當景觀面積大于區(qū)域12時，運用上述景觀指數(shù)法進行景觀特征分析可行有效。測量所得3個年份居民點建成區(qū)所有緩沖區(qū)面積中，面積最小的為1986年第一個緩沖區(qū)形成的第一梯度圈，面積為212 175.9 hm2，大于面積22 500 hm2，即指數(shù)的計算結(jié)果基本不受幅度影響。

2.3.1 歸一化平均斑塊面積 從居民點緩沖區(qū)梯度分異來看，1986年從梯度圈1到梯度圈2有明顯增加趨勢，往后趨于平緩；1996年從梯度圈1到梯度圈5有明顯增加趨勢，往后趨于平緩；2001年從梯度圈1到梯度圈5有明顯增加趨勢，往后趨于平緩（圖5）。說明對于歸一化平均斑塊面積，1986年居民點建成區(qū)影響范圍為周圍2 km，1996年和2001年隨著居民點面積擴大，影響范圍為周圍5 km。

2.3.2 雙對數(shù)分維數(shù) 隨著居民點緩沖區(qū)范圍的擴大，指數(shù)值先增大后趨緩。1986年從梯度圈1到梯度圈4變化幅度明顯，往后趨緩；1996年從梯度圈1到梯度圈5變化幅度明顯，往后趨緩；2001年從梯度圈1到梯度圈6變化幅度明顯，往后趨緩（圖6）。說明對于雙對數(shù)分維數(shù)，1986年居民點建成區(qū)影響范圍為周圍4 km，隨著居民點建成區(qū)面積的增大，1996年和2001年分別擴大為5 km和6 km。

2.3.3 聚集度 3個年度聚集度變化趨勢基本相同，先逐漸下降后基本保持平穩(wěn)。1986年從梯度圈1到梯度圈5指數(shù)值不斷下降，往后趨于平緩；1996年從梯度圈1到梯度圈4指數(shù)值不斷下降，往后趨于平緩；2001年從梯度圈1到梯度圈6指數(shù)值不斷下降，往后趨于平緩（圖7）。說明對于聚集度指數(shù)，1986年居民點建成區(qū)影響范圍為周圍5 km，1996年居民點建成區(qū)面積增加，影響范圍卻縮小為4 km，2001年居民點建成區(qū)面積繼續(xù)增加，影響范圍擴大為6 km。

2.3.4 反差矩 3個年份變化趨勢相似，隨著緩沖區(qū)帶擴大，數(shù)值先減少后增加再趨于平緩。1986年從梯度圈1到梯度圈5數(shù)值顯著下降，往后趨緩；1996年從梯度圈1到梯度圈5數(shù)值顯著下降，往后趨緩；2001年從梯度圈1到梯度圈6數(shù)值顯著下降，變化幅度明顯，往后趨緩（圖8）。說明對于反差矩指數(shù)，1986年居民點建成區(qū)影響范圍為5 km，1996年居民點建成區(qū)面積增加，影響范圍保持不變，為周圍5 km，2001年居民點建成區(qū)影響范圍擴大為6 km。

2.3.5 多樣性 與前4個指數(shù)相同，3個年度多樣性指數(shù)值呈現(xiàn)相似的變化趨勢，隨著緩沖區(qū)范圍的擴大，多樣性指數(shù)值不斷擴大直至相對穩(wěn)定。1986年從梯度圈1到梯度圈5數(shù)值顯著上升，往后趨緩；1996年從梯度圈1到梯度圈5數(shù)值顯著上升，往后趨緩；2001年從梯度圈1到梯度圈6數(shù)值顯著上升，變化幅度明顯，往后趨緩（圖9）。說明對于多樣性指數(shù)，1986年居民點建成區(qū)影響范圍為周圍5 km，1996年居民點建成區(qū)面積增加，影響范圍保持為5 km，2001年居民點建成區(qū)面積繼續(xù)增加，影響范圍擴大為6 km。

3 討論

土地利用覆被變化能直接反映人類活動這一當今全球環(huán)境變化的主導(dǎo)因子。影響黃河三角洲的土地利用覆被變化的人類活動集中在道路、居民點、油井、農(nóng)田、葦田和蝦蟹田開發(fā)等，關(guān)于居民點建成區(qū)對土地利用覆被變化的研究為本次研究的重點。

景觀研究受到粒度和幅度影響，數(shù)據(jù)主要來源于分辨率為30 m的TM影像，考慮原數(shù)據(jù)粒度特點采用分辨率30 m；通過分析，當研究面積大于22 500 hm2時，景觀研究受幅度影響明顯減小，研究區(qū)1986年、1996年和2001年居民點面積分別為23 102.4、33 110.5和39 397.8 hm2，均大于幅度影響范圍，以它們?yōu)橹行男纬傻木用顸c緩沖區(qū)面積更是大于幅度影響范圍。居民點建成區(qū)對周圍地區(qū)具有輻射影響特性，它的變化會帶動周圍景觀的變化，如居民點建成區(qū)周圍道路、農(nóng)田等將會有所增加。在此基礎(chǔ)上進一步研究15年以來居民點建成區(qū)變化對土地利用格局影響，通過分析，表現(xiàn)出以下規(guī)律。

基于緩沖區(qū)技術(shù)的居民點建成區(qū)輻射梯度圈景觀分析結(jié)果顯示，居民點建成區(qū)存在明顯的緩沖區(qū)效應(yīng)。研究年段內(nèi)，隨著居民點建成區(qū)面積不斷增加，各年份緩沖區(qū)梯度圈景觀指數(shù)變化區(qū)間也不斷變化，數(shù)據(jù)變化說明隨著居民點建成區(qū)范圍的增大，各指數(shù)的響應(yīng)范圍均有所擴大；1986年歸一化平均斑塊面積指數(shù)受影響范圍為居民點建成區(qū)周圍2 km距離，雙對數(shù)分維數(shù)為4 km，聚集度指數(shù)為5 km，反差矩指數(shù)為5 km，多樣性指數(shù)為5 km；1996年歸一化平均斑塊面積指數(shù)受影響范圍為居民點建成區(qū)周圍5 km距離，雙對數(shù)分維數(shù)為5 km，聚集度指數(shù)為4 km，反差矩指數(shù)為5 km，多樣性指數(shù)為5 km；2001年歸一化平均斑塊面積指數(shù)受影響范圍為居民點建成區(qū)周圍5 km距離，雙對數(shù)分維數(shù)為6 km，聚集度指數(shù)為6 km，反差矩指數(shù)為6 km，多樣性指數(shù)為6 km。1996年雖然居民點建成區(qū)面積增大，但是部分指數(shù)受影響的距離和1986年持平，聚集度指數(shù)影響距離甚至還出現(xiàn)比1986年下降的趨勢，但從影響面積來看，1986年居民點建成區(qū)緩沖區(qū)5形成的第五梯度圈包含的范圍為435 008 hm2，1996年居民點建成區(qū)緩沖區(qū)5形成的第五梯度圈包含的范圍為451 238 hm2，影響面積1996年較大；隨著居民點建成區(qū)面積的擴大，2001年各指數(shù)相應(yīng)距離明顯大于1986年和1996年，故相應(yīng)的范圍也將大于1986年和1996年。綜上得出，隨著居民點建成區(qū)面積的增大，各指數(shù)相應(yīng)的范圍呈現(xiàn)擴大的趨勢，即隨著居民點建成區(qū)面積的增大，其對周圍地區(qū)影響的范圍亦不斷增大。

3S技術(shù)在土地利用覆被變化中的運用為獲得區(qū)域發(fā)展時間和空間的動態(tài)過程信息提供了更加便捷的方法，在此基礎(chǔ)上運用梯度和景觀指數(shù)相結(jié)合的方法，分析居民點建成區(qū)變化對區(qū)域生態(tài)環(huán)境的影響范圍，可以量化居民點建成區(qū)對景觀生態(tài)影響的梯度特征，為區(qū)域規(guī)劃、區(qū)域生態(tài)環(huán)境保護和管理提供決策依據(jù)。

參考文獻：

[1] 李秀彬.中國近20年來耕地面積的變化及其政策啟示[J].自然資源學(xué)報，1999，14（4）：329-333.

[2] 張 明.區(qū)域土地利用結(jié)構(gòu)及其驅(qū)動因子的統(tǒng)計分析[J].自然資源學(xué)報，1999，14（4）：381-384.

[3] 張惠遠，趙昕奕，蔡運龍，等.喀斯特山區(qū)土地利用變化的人類驅(qū)動機制研究——以貴州省為例[J].地理研究，1999，18（2）：136-142.

[4] 史培軍，陳 晉，潘耀忠.深圳市土地利用變化機制分析[J]. 地理學(xué)報，2000，55（2）：151-160.

[5] 高 峻，宋永昌.上海西南城市干道兩側(cè)地帶景觀動態(tài)研究[J].應(yīng)用生態(tài)學(xué)報，2001，12（4）：605-609.

[6] 吳啟焰，任東明.改革開放以來我國城市地域結(jié)構(gòu)演變與持續(xù)發(fā)展研究——以南京都市區(qū)為例[J].地理科學(xué)，1994，19（2）：13-18.

[7] 陳 力.舊城更新中城市形態(tài)的延續(xù)與創(chuàng)新[J]. 華僑大學(xué)學(xué)報（自然科學(xué)版），1997（1）：58-61.

[8] 陳 泳.近現(xiàn)代蘇州城市形態(tài)演化研究[J].城市規(guī)劃匯刊，2003（6）：62-71.

[9] 胡賢輝，楊鋼橋，張 霞，等.農(nóng)村居民點用地數(shù)量變化及驅(qū)動機制研究——基于湖北仙桃市的實證[J].資源科學(xué)，2007，29（3）：191-197.

[10] 孫華生，黃敬峰，金 艷，等.基于GIS技術(shù)的縣域居民點空間分布特征分析及其優(yōu)化布局[J].浙江大學(xué)學(xué)報（農(nóng)業(yè)與生命科學(xué)版），2007，33（3）：348-354.

[11] MATTHEW L， WU J G. A gradient analysis of urban landscape pattern： A case study from the Phoenix metropolitan region，Arizona，USA[J]. Landscape Ecology，2002，17（4）：327-339.

[12] 孫 娟，夏漢平，藍崇鈺，等.基于緩沖帶的貴港市城市景觀格局梯度分析[J].生態(tài)學(xué)報，2006，26（3）：654-662.

[13] 劉黎明.土地資源學(xué)[M].北京：中國農(nóng)業(yè)大學(xué)出版社，2002.84-86.

[14] HAY G J，MARCEAU D J，DUBE P， et al. A multiscale framework for landscape analysis： Object-specific analysis and upscaling[J]. Landscape Ecology，2001，16（6）：471-490.

[15] 布仁倉，李秀珍，胡遠滿，等.尺度分析對景觀格局指標的影響[J].應(yīng)用生態(tài)學(xué)報，2003，14（2）：2181-2186.

[16] WU J G. Effects of changing scale on landscape pattern analysis： Scaling relations[J]. Landscape Ecology，2004，19（2）：125-138.

[17] SAURA S. Effects of remote sensor spatial resolution and data aggregation on selected fragmentation indices[J]. Landscape Ecology，2004，19（2）：197-209.

[18] WIMBERLY M C， OHMANN J L. A multi-scale assessment of human and environmental constraints on forest land cover change on the Oregon（USA） coast range[J]. Landscape Ecology， 2004，19（6）：631-646.

[19] 申衛(wèi)軍，鄔建國，任 海，等.空間幅度變化對景觀格局分析的影響[J].生態(tài)學(xué)報，2003，23（11）：2219-2231.

[20] 張 娜.生態(tài)學(xué)中的尺度問題：內(nèi)涵與分析方法[J].生態(tài)學(xué)報，2006，26（7）：2340-2355.

[21] RIITTERS K H， ONEILL R V， HUNSAKER C T， et al. A factor analysis of landscape pattern and structure metrics[J]. Landscape Ecology，1995，10（1）：23-39.

[22] TRAUB B， KLEINN C. Measuring fragmentation and structural diversity[J]. Forstw Centralblatt，1999，118（1）：39-50.

[23] 布仁倉，胡遠滿，常 禹，等.景觀指數(shù)之間的相關(guān)分析[J].生態(tài)學(xué)報，2005，25（10）：2764-2775.

[24] MCGARIGAL K， CUSHMAN S A， NEEL M C M， et al. Fragstats： Spatial pattern analysis program for categorical maps computer software program produced by the authors at the University of Massachusetts，Amherst[EB/OL]. http：//www.umassedu/landeco/research/fragstats/fragstats.html，2002.

[25] 王緒高，李秀珍，賀紅士，等.中性景觀模型與真實景觀的一致性[J].應(yīng)用生態(tài)學(xué)報，2004，16（6）：973-978.