基于J48決策樹(shù)算法的遙感土地利用變化分析

史澤鵬，馬中文，馬友華，王強(qiáng)，於忠祥

(1.安徽農(nóng)業(yè)大學(xué) 資源環(huán)境與信息技術(shù)研究所，合肥 230036；2.安徽農(nóng)業(yè)大學(xué)國(guó)土資源研究所，合肥 230036)

1 引 言

土地利用變化是人類(lèi)社會(huì)發(fā)展的一面鏡子，它反映了人類(lèi)與自然界相互影響、相互作用的一種特殊而又復(fù)雜的涉及土地資源自然屬性和人類(lèi)利用方式的社會(huì)自然現(xiàn)象[1]。早在1993年，“國(guó)際地圈與生物圈計(jì)劃(IGBP)”和“全球變化人類(lèi)影響和響應(yīng)計(jì)劃(HDP)”就共同將土地利用/土地覆蓋變化(LUCC)列為全球變化研究的核心計(jì)劃[2]。國(guó)際上利用遙感技術(shù)對(duì)土地利用的研究至少可以追溯至20世紀(jì)20年代。到了20世紀(jì)70年代，美國(guó)第一顆陸地球衛(wèi)星(Landsat)的發(fā)射，開(kāi)創(chuàng)了利用衛(wèi)星遙感技術(shù)進(jìn)行大范圍土地利用調(diào)查的新紀(jì)元[3]。到了80年代后期，隨著遙感技術(shù)與地理信息系統(tǒng)在土地研究中的廣泛應(yīng)用，土地利用動(dòng)態(tài)變化監(jiān)測(cè)得到了迅速的發(fā)展，并開(kāi)始了以此為基礎(chǔ)的土地利用優(yōu)化決策的定量分析研究[4]。到90年代中期，則開(kāi)始利用TM影像按州進(jìn)行土地利用調(diào)查[3]。而國(guó)內(nèi)這這方面的研究則起步較晚，直到20世紀(jì)70年代，才開(kāi)始采用航空遙感與衛(wèi)星相結(jié)合的方法進(jìn)行土地資源調(diào)查。其中“中科院780工程”的“騰沖縣大型綜合航空遙感試驗(yàn)研究項(xiàng)目”，取得了豐碩的成果。在“七五”、“八五”及“九五”期間，國(guó)家逐漸將RS、GIS和GPS技術(shù)綜合應(yīng)用研究列為國(guó)家科技攻關(guān)重點(diǎn)項(xiàng)目[5-7]。

隨著遙感和GIS技術(shù)的快速發(fā)展，很多學(xué)者開(kāi)始利用遙感影像和GIS技術(shù)獲取土地利用變化信息[8-9]。而選擇合適的遙感分類(lèi)方法則是遙感和GIS的土地利用變化分析領(lǐng)域的研究熱點(diǎn)。因此，本文以合肥市肥東縣為研究區(qū)域，在遙感和GIS的支持下，采用WEKA的J48決策樹(shù)算法和監(jiān)督分類(lèi)中的最大似然法兩種分類(lèi)方法提取肥東縣2000年～2011年的土地利用變化信息并進(jìn)行分析，對(duì)于認(rèn)識(shí)和了解肥東縣甚至合肥市區(qū)域土地利用變化和社會(huì)經(jīng)濟(jì)、生態(tài)環(huán)境的健康、協(xié)調(diào)、可持續(xù)發(fā)展具有十分重要的意義。

2 研究區(qū)概況

肥東縣地處江淮之間，居皖中腹地，位于省會(huì)合肥市東側(cè)，北緯31°35′～北緯32°16′，東經(jīng)117°19′～東經(jīng)117°52′。交通便捷，區(qū)位優(yōu)越，南臨巢湖，北接定遠(yuǎn)縣，東與全椒縣、巢湖市接壤，西與長(zhǎng)豐縣、合肥市相連，史稱(chēng)“吳楚要沖、包公故里”。地勢(shì)北高南低，縣域總面積2206km2，人口109.5萬(wàn)，轄18個(gè)鄉(xiāng)鎮(zhèn)、3個(gè)開(kāi)發(fā)園區(qū)、340個(gè)行政村。屬北亞熱帶季風(fēng)性氣候，年均降水量960ml，蒸發(fā)量1483ml，蒸發(fā)量過(guò)大使地下水位有所降低，光照充足，無(wú)霜期較長(zhǎng)，達(dá)235天，年平均氣溫在15℃。土壤類(lèi)型以水稻土和黃棕壤為主。受東南太平洋季風(fēng)氣候影響，水、熱、光資源均較為豐富，適宜多種植被生長(zhǎng)。

3 研究材料與方法

3.1 研究數(shù)據(jù)和軟件

本文的研究數(shù)據(jù)主要有涵蓋研究區(qū)的2000年、2006年、2011年3期美國(guó)Landsat-7 ETM+影像和2011年德國(guó)RapidEye影像(表1)；地形圖(比例尺1∶50000)、肥東縣行政區(qū)劃圖(比例尺1∶50000)。

本文運(yùn)用的主要軟件有WEKA 3.6(用于J48數(shù)據(jù)挖掘)、ENVI 4.7(用于遙感圖像處理)、ArcGIS 10.0(用于土地利用變化分析)。

表1 遙感數(shù)據(jù)列表

3.2 研究方法

本文主要采用WEKA的J48決策樹(shù)算法和監(jiān)督分類(lèi)中的最大似然法兩種方法進(jìn)行遙感圖像解譯，提取肥東縣3期遙感影像的土地利用變化信息；并利用ArcGIS 10.0對(duì)遙感圖像解譯結(jié)果進(jìn)行空間疊加分析，為土地利用轉(zhuǎn)換矩陣的建立提供基礎(chǔ)；通過(guò)一些土地利用變化模型來(lái)進(jìn)行土地利用變化分析。本研究的技術(shù)流程圖如圖1所示。

4 遙感影像分類(lèi)處理

4.1 遙感數(shù)據(jù)預(yù)處理

因購(gòu)買(mǎi)的涵蓋研究區(qū)的RapidEye影像已經(jīng)經(jīng)過(guò)輻射校正和幾何校正等相關(guān)預(yù)處理，故本文主要對(duì)3期ETM+影像進(jìn)行大氣校正、幾何精校正、圖像裁剪等預(yù)處理。大氣校正選擇了ENVI中的Quick Atmospheric Correction功能進(jìn)行快速大氣校正；幾何精校正的檢驗(yàn)均方差RMS(Root Means Squared)小于0.7個(gè)像元，并選用了雙線性插值法(Bilinear)對(duì)影像進(jìn)行了重采樣；使用了包含研究區(qū)行政界矩形框的規(guī)則裁剪方式對(duì)研究區(qū)3期ETM+影像進(jìn)行切割。

圖1 技術(shù)流程圖

4.2 遙感圖像分類(lèi)

本文以全國(guó)遙感監(jiān)測(cè)土地利用/覆蓋分類(lèi)體系為基礎(chǔ)，并結(jié)合研究區(qū)實(shí)際和本文的研究目的建立了本研究的遙感分類(lèi)體系，共為6類(lèi)：耕地、林地、草地、水域、建設(shè)用地、未利用地。在分類(lèi)方法上，監(jiān)督分類(lèi)方面，選擇了監(jiān)督分類(lèi)中最常用、也是精度相對(duì)較高的最大似然法[10]；決策樹(shù)分類(lèi)方面，本文主要用WEKA的J48決策樹(shù)算法。WEKA的J48決策樹(shù)算法實(shí)際上是對(duì)Quinlan的C4.5決策樹(shù)算法的實(shí)現(xiàn)，并增加了很好的剪枝過(guò)程，其優(yōu)點(diǎn)在于訓(xùn)練出的決策樹(shù)規(guī)則易于理解且準(zhǔn)確度較高，但對(duì)前期訓(xùn)練樣本采集的準(zhǔn)確性要求也較高。

本文運(yùn)用WEKA的J48決策樹(shù)分類(lèi)算法的具體步驟是：通過(guò)在Google Earth中采集樣本(保證較高的樣本精確度)，用采樣樣本在ArcGIS 10.0中提取遙感影像柵格值，將提取的柵格值進(jìn)行標(biāo)準(zhǔn)化處理后，再利用WEKA軟件進(jìn)行數(shù)據(jù)挖掘，采用J48算法對(duì)樣本進(jìn)行訓(xùn)練，獲取決策樹(shù)規(guī)則，再將規(guī)則帶入ENVI4.7中進(jìn)行決策樹(shù)分類(lèi)，得到分類(lèi)結(jié)果，流程圖見(jiàn)圖2。

圖2 J48決策樹(shù)分類(lèi)流程圖

4.3 分類(lèi)后處理及分類(lèi)精度驗(yàn)證

本文選用了ENVI中的Majority分析以及結(jié)合人機(jī)交互方式分別對(duì)分類(lèi)結(jié)果進(jìn)行了處理，并利用ENVI的精度評(píng)價(jià)功能，選取50個(gè)地面數(shù)據(jù)并結(jié)合Google Earth目視解譯等方法，組成驗(yàn)證樣本，構(gòu)建混淆矩陣，對(duì)兩種分類(lèi)結(jié)果進(jìn)行了精度驗(yàn)證，精度比較見(jiàn)表2。從表2中可以看出決策樹(shù)分類(lèi)的總精確度和Kappa系數(shù)均高于監(jiān)督分類(lèi)，且決策樹(shù)分類(lèi)的總精確度均在85%以上。故本文采用J48決策樹(shù)分類(lèi)的結(jié)果進(jìn)行土地利用變化信息分析。

表2 監(jiān)督分類(lèi)與決策樹(shù)分類(lèi)精度比較

5 土地利用變化分析

土地利用變化主要體現(xiàn)在土地利用類(lèi)型的變化、數(shù)量和質(zhì)量的變化、土地利用結(jié)構(gòu)和程度的變化以及土地利用區(qū)域差異的變化等方面[11]。本文主要從土地利用類(lèi)型的數(shù)量、結(jié)構(gòu)和程度3個(gè)方面對(duì)肥東縣土地利用變化進(jìn)行分析。

5.1 土地利用類(lèi)型轉(zhuǎn)換矩陣分析

按照本文建立的分類(lèi)體系和遙感影像分類(lèi)步驟，通過(guò)利用決策樹(shù)分類(lèi)的結(jié)果，提取了肥東縣2000年、2006年、2011年的土地利用信息，并利用ArcGIS的制圖功能得到肥東縣2000年、2006年、2011年的土地利用分類(lèi)圖(圖3)。

圖3 2000年～2011年肥東縣土地利用分類(lèi)圖

從3期土地利用分類(lèi)圖中可以明顯看出，肥東縣建設(shè)用地不斷增加，耕地不斷減少。其中，建設(shè)用地主要由肥東縣城向四周擴(kuò)散；2000年～2006年間，主要向西南方向擴(kuò)張；2006年～2011年間，主要向西北、東南方向擴(kuò)張；林地、草地、水域、未利用地均呈減少趨勢(shì)。

為了量化各土地利用類(lèi)型的變化情況，將3期分類(lèi)結(jié)果在ArcGIS中進(jìn)行空間疊加并進(jìn)行處理，最后構(gòu)建土地利用類(lèi)型轉(zhuǎn)換矩陣(見(jiàn)表3、表4)。

表3 2000年～2006年肥東縣土地利用類(lèi)型轉(zhuǎn)換矩陣 單位：hm2

表4 2006年～2011年肥東縣土地利用類(lèi)型轉(zhuǎn)換矩陣 單位：hm2

從表3、表4中可以看出，2000年～2006年和2006年～2011年間，肥東縣土地利用類(lèi)型轉(zhuǎn)換面積最大的為均為耕地轉(zhuǎn)換為建設(shè)用地，轉(zhuǎn)化面積分別為26080.41hm2和31607.48hm2。同時(shí)段耕地轉(zhuǎn)換為其他土地利用類(lèi)型空間分布圖見(jiàn)圖4。

圖4 肥東縣耕地轉(zhuǎn)換為其他用地類(lèi)型圖

5.2 土地利用類(lèi)型的數(shù)量變化分析

本文主要通過(guò)分析土地利用變化幅度、單一土地利用類(lèi)型變化率和綜合土地利用年變化率，來(lái)了解土地利用類(lèi)型的總量變化和總體態(tài)勢(shì)。

5.2.1 土地利用變化幅度

土地利用變化幅度反映的是某研究區(qū)一定時(shí)間范圍內(nèi)土地利用類(lèi)型的數(shù)量變化情況[12]。土地利用變化率更好地反映了土地利用變化幅度。結(jié)合肥東縣土地利用類(lèi)型轉(zhuǎn)換矩陣的相關(guān)數(shù)據(jù)以及土地利用變化幅度、變化率公式(式1、式2)，計(jì)算得到肥東縣2個(gè)研究時(shí)段的土地利用變化幅度及變化率(圖5、圖6)。

ΔU=Uj-Ui

(1)

(2)

其中，ΔU表示研究時(shí)段內(nèi)某一土地利用類(lèi)型的變化幅度；S表示研究時(shí)段內(nèi)某一土地利用類(lèi)型的變化率；Ui表示研究初期某一土地利用類(lèi)型的面積；Uj表示研究末期某一土地利用類(lèi)型的面積。

圖5 2000年～2011年肥東縣土地利用類(lèi)型數(shù)量變化幅度

從圖5、圖6中可以看出，2000年～2006年和2006年～2011年期間，建設(shè)用地面積均增加，除2000年～2006年期間水域面積略有增加外，其他土地利用類(lèi)型均呈減少趨勢(shì)。兩個(gè)研究時(shí)段的土地利用類(lèi)型變化幅度和變化率最大的均是建設(shè)用地，分別為25332.06hm2、64.00%和31632.18hm2、48.73%。

圖6 2000年～2011年肥東縣土地利用變化率

5.2.2 單一土地利用類(lèi)型年變化率

單一土地利用類(lèi)型年變化率，又叫單一土地利用類(lèi)型動(dòng)態(tài)度，它是以研究時(shí)段內(nèi)單一土地利用類(lèi)型的變化面積為基礎(chǔ)，可以直觀地反映單一土地利用類(lèi)型的動(dòng)態(tài)度。單一土地利用類(lèi)型年變化量和年變化率計(jì)算模型(式3、式4)：

(3)

(4)

式3、式4中，L表示研究區(qū)某一土地利用類(lèi)型年變化量；K表示研究區(qū)某一土地利用類(lèi)型年變化率；T表示研究時(shí)段年數(shù)；ΔU同式(1)；S同式(2)。

通過(guò)公式(3)、公式(4)計(jì)算得到肥東縣2個(gè)研究時(shí)段各土地利用類(lèi)型的年變化量和變化率，如表5所示。

從表5中可以看出，2000年～2006年和2006年～2011年期間，建設(shè)用地的年變化量和年變化率最大的均為建設(shè)用地，分別為4222.01hm2、10.67%和6326.44hm2、9.75%。

表5 2000年～2011年肥東縣單一土地利用類(lèi)型年變化量及變化率

5.2.3 綜合土地利用年變化率

綜合土地利用年變化率，又叫綜合土地利用動(dòng)態(tài)度，它能反映研究區(qū)某一時(shí)間段土地利用總體趨勢(shì)[5]。其計(jì)算模型(式5)：

(5)

式5中，Z表示研究區(qū)某一時(shí)間段內(nèi)綜合土地利用年變化率；ΔZx-y表示研究時(shí)段內(nèi)第x類(lèi)土地利用類(lèi)型轉(zhuǎn)換為非x類(lèi)土地利用類(lèi)型面積的絕對(duì)值；Zx表示研究初期第x類(lèi)土地利用類(lèi)型面積；T表示研究時(shí)段年數(shù)。

按照土地利用年變化率可以將中國(guó)土地利用的時(shí)間動(dòng)態(tài)特征劃分為4種類(lèi)型[13]，見(jiàn)表6。

表6 中國(guó)土地利用時(shí)間動(dòng)態(tài)特征類(lèi)型

通過(guò)公式5計(jì)算得到2000年～2006年、2006年～2011年兩個(gè)研究時(shí)段的綜合土地利用年變化率分別為7.68和8.48，再結(jié)合表6，可以發(fā)現(xiàn)研究時(shí)段內(nèi)，肥東縣土地利用變化均處于慢速變化類(lèi)型，但很明顯變化速度越來(lái)越快。由此可見(jiàn)，在不久的將來(lái)，肥東縣土地利用變化很快就會(huì)進(jìn)入快速變化時(shí)期。

5.3 土地利用結(jié)構(gòu)變化分析

運(yùn)用遙感分類(lèi)結(jié)果和ArcGIS的空間數(shù)據(jù)庫(kù)等功能，得到肥東縣2000年、2006年、2011年土地利用結(jié)構(gòu)信息(表7)。

表7 2000年～2011年肥東縣土地利用類(lèi)型面積及比例

從表7中可以看出，肥東縣的耕地面積占土地總面積的比例由2000年的55.41%下降達(dá)到2006年的45.75%，到2011年，耕地面積比例則下降到39.73%；建設(shè)用地占土地總面積中的比例由2000年的17.94%上升到2006年的29.43%，2011年時(shí)，建設(shè)用地則上升到43.77%；草地、林地和未利用地在2000年～2011年時(shí)間段均呈減少趨勢(shì)。

5.4 土地利用程度變化分析

一個(gè)地區(qū)的土地利用程度的變化是該區(qū)域中多種土地利用類(lèi)型綜合變化的結(jié)果，它可以定量地反映該區(qū)域土地利用的綜合水平和變化趨勢(shì)[14]。在中國(guó)資源環(huán)境數(shù)據(jù)庫(kù)中，劉紀(jì)遠(yuǎn)等從生態(tài)學(xué)的角度提出了土地利用程度分級(jí)標(biāo)準(zhǔn)(表8)[15]。

表8 土地利用類(lèi)型分級(jí)指數(shù)表

土地利用程度及其變化可分別用土地利用程度綜合指數(shù)(式6)[ 12、16]、土地利用程度變化量(式7)和變化率(式8)[12]表示。

(6)

(7)

(8)

上式中，Ma和Mb分別表示a時(shí)段和b時(shí)段土地利用程度綜合指數(shù)，土地利用程度指數(shù)介于100～400之間；Ax表示第x級(jí)土地利用程度綜合指數(shù)；Cx表示第x級(jí)土地利用類(lèi)型面積百分比；Cxa和Cxb分別表示a時(shí)段和b時(shí)段第x級(jí)土地利用類(lèi)型面積百分比；ΔMb-a表示土地利用程度變化量；R表示土地利用程度變化率。

若ΔMb-a>0，表示該區(qū)域土地利用處于發(fā)展期；若ΔMb-a<0，表示該區(qū)域土地利用處于衰退期或調(diào)整期。根據(jù)上式計(jì)算得到表9。

表9 2000年～2011年肥東縣土地利用程度綜合指數(shù)、變化量及變化率

從表9中可以看出，肥東縣2000年、2006年、2011年的土地利用程度綜合指數(shù)分別為289.64、303.57、327.53，這與土地利用程度綜合指數(shù)的最大值400相比，說(shuō)明研究時(shí)段內(nèi)肥東縣土地利用程度綜合指數(shù)值都比較高且呈增加趨勢(shì)，與肥東縣地形地貌和經(jīng)濟(jì)發(fā)展水平相符合；2000年～2006年和2006年～2011年肥東縣土地利用程度變化量、變化率分別為13.12、4.81%和23.96、7.89%，研究時(shí)段內(nèi)變化量均大于0，說(shuō)明肥東縣在研究時(shí)段內(nèi)土地利用處于發(fā)展期，且后一階段變化量和變化率均大于前一階段，說(shuō)明肥東縣土地利用變化越來(lái)越快。

6 結(jié)束語(yǔ)

在遙感圖像解譯方法上，采用WEKA的J48決策樹(shù)算法的分類(lèi)精度要高于監(jiān)督分類(lèi)中的最大似然法。J48決策樹(shù)算法增加了很好的剪枝過(guò)程，其優(yōu)點(diǎn)在于訓(xùn)練出的的決策樹(shù)規(guī)則比較容易理解并且準(zhǔn)確度較高，因此，利用其訓(xùn)練出的規(guī)則進(jìn)行土地利用分類(lèi)，可以得到較高的分類(lèi)精度。但是運(yùn)用該方法的前提是一定要保證前期采集的訓(xùn)練樣本的準(zhǔn)確性，否則就會(huì)降低最終的分類(lèi)精度。

在土地利用變化分析方面，肥東縣2000年～2011年的11年間土地利用類(lèi)型在數(shù)量、結(jié)構(gòu)和程度方面發(fā)生了較大的變化。在數(shù)量方面，建設(shè)用地在兩個(gè)研究時(shí)段的土地利用變化幅度和變化率均最大；2000年～2006年和2006年～2011年肥東縣的綜合土地利用年變化率分別為7.68%和8.48%，盡管兩個(gè)研究時(shí)段均處于土地慢速變化型的范圍內(nèi)，但可以預(yù)見(jiàn)，隨著經(jīng)濟(jì)的快速發(fā)展，肥東縣的土地利用變化速度正在逐漸加快。在結(jié)構(gòu)方面，建設(shè)用地在肥東縣土地總面積中所占比例從2000年的17.94%上升到2011年的43.77%，耕地在肥東縣土地總面積中所占比例從2000年的55.41%下降到2011年的39.73%，草地、水域、林地和未利用地所占比例在研究時(shí)段內(nèi)均呈減少趨勢(shì)。在程度方面，2000年、2006年和2011年肥東縣的土地利用程度綜合指數(shù)分別為289.64、303.57和327.53，說(shuō)明肥東縣的土地利用綜合水平隨著經(jīng)濟(jì)發(fā)展在逐步提高且速度越來(lái)越快。

參考文獻(xiàn)：

[1] 蔡運(yùn)龍.土地利用/土地覆被變化研究：尋求新的綜合途徑[J].地理研究，2001，20(6)：645-652.

[2] TURNER II B L，MASS R H，SKOLE D L.Relating land use and global land-cover chrange[R].IGBP Report No.24 and HDP Report No.5，1993.

[3] 王瑞國(guó).遙感土地利用動(dòng)態(tài)監(jiān)測(cè)變化信息自動(dòng)提取與實(shí)證研究[D].長(zhǎng)春：東北農(nóng)業(yè)大學(xué)，2007.

[4] 李學(xué)梅，李忠峰.土地利用/覆蓋變化研究進(jìn)展及其意義[J].安徽農(nóng)業(yè)科學(xué)，2008，36(6)：2462-2464.

[5] 劉紀(jì)遠(yuǎn)，布和敖斯?fàn)?中國(guó)土地利用變化現(xiàn)代過(guò)程時(shí)空特征的研究：基于衛(wèi)星遙感數(shù)據(jù)[J].第四紀(jì)研究，2000，20(3)：229-239.

[6] 莊大方，劉紀(jì)遠(yuǎn).中國(guó)土地利用程度的區(qū)域分異模型研究[J].自然資源學(xué)報(bào)，1997，12(2)：105-111.

[7] 鄭國(guó)強(qiáng)，江南，于興修.長(zhǎng)江下游沿線土地利用區(qū)域結(jié)構(gòu)演化分析[J].自然資源學(xué)報(bào)，2003.18(5)：568-572.

[8] QUARMBY N A，CUSHNLE J L.Monitoring urban land cover changes at the urban hinge from SPOT HRV iraagery south-east England[J].International Journal of Remote Sensing，1989(10)：953-963.

[9] 戴聲佩，張勃.利用Landsat影像構(gòu)建河西綠洲土地利用信息圖譜—以張掖市甘州區(qū)為例[J].遙感信息，2012，27(5)：107-114.

[10] 史澤鵬，王玉佳，馬友華，等.遙感影像土地利用/覆蓋分類(lèi)方法研究進(jìn)展[J].中國(guó)農(nóng)學(xué)通報(bào)，2012，28(12)：273-278.

[11] 劉紀(jì)遠(yuǎn).中國(guó)資源環(huán)境遙感宏觀調(diào)查與動(dòng)態(tài)研究[M].北京，中國(guó)科學(xué)技術(shù)出版社，1996：158-188.

[12] 王秀蘭，包玉梅.土地利用動(dòng)態(tài)變化研究方法探討[J].地理科學(xué)進(jìn)展，1999，18(1)：81-86.

[13] 于思遠(yuǎn)，劉紀(jì)遠(yuǎn)，張?jiān)鱿?，?中國(guó)土地利用時(shí)空特征分析[J].地理學(xué)報(bào)，2001，56(6)：631-639.

[14] 劉靜玉，錢(qián)樂(lè)祥，苗長(zhǎng)虹，等.基于RS與GIS技術(shù)的豫西山地典型區(qū)域LUCC動(dòng)態(tài)監(jiān)測(cè)研究[J].河南大學(xué)學(xué)報(bào)(自然科學(xué)版)，2005，35(2)：47-51.

[15] LIU J Y.The macro investigation and dynamic research of the resource and environment[M]，Beijing：Science and Technology Press，1996：158-188.

[16] 莊大方，劉紀(jì)遠(yuǎn).中國(guó)土地利用程度的區(qū)域分異模型研究[J].自然資源學(xué)報(bào)，1997，12(2)：105-111.