Fragstats與Arcgis協(xié)同分析解析課件

1、應(yīng)用Fragstats軟件（本操作版本為Fragstats3.3）（在Google搜Fragstats 的第一個(gè)鏈接）Fragstats軟件主頁/landeco/research/fragstats/fragstats.html 景觀分析LandscapeAnalysisBy fhys第1頁，共54頁。在開始之前我們先再了解景觀分析（Landscape Analyst）的相關(guān)重要概念Fragstats軟件強(qiáng)大，可計(jì)算出50多個(gè)景觀指標(biāo)。這些指標(biāo)被分為三組級(jí)別：分別為Patch（斑塊/斑塊）反映單塊的結(jié)果特征Class（斑塊類型）反映屬于相同類型的塊的特征Landscape（整體景觀）反映整體的

2、結(jié)構(gòu)特征三種級(jí)別逐步擴(kuò)大尺度，高度相關(guān)這三個(gè)層次的指標(biāo)也基本是所有景觀分析中都會(huì)用到的指標(biāo)景觀分析LandscapeAnalysis舉例來說：每個(gè)斑塊的周長就屬于Patch級(jí)別；斑塊密度（某類型斑塊數(shù)量除以該類型所有斑塊面積之和）就屬于Class級(jí)別；整個(gè)區(qū)域內(nèi)的斑塊平均面積就屬于Landscape級(jí)別；其實(shí)在某種類型上也可以有斑塊平均面積，類似地計(jì)算就是了。第2頁，共54頁。下面列出一些重要的景觀分析指標(biāo)：景觀面積（TA）公式描述：CA等于某一斑塊類型中所有斑塊的面積之和（m2），除以10000后轉(zhuǎn)化為公頃（ha）；即某斑塊類型的總面積。生態(tài)意義：CA度量的是景觀的組分，也是計(jì)算其它指標(biāo)的基

3、礎(chǔ)。它有很重要的生態(tài)意義，其值的大小制約著以此類型斑塊作為聚居地（Habitation）的物種的豐度、數(shù)量、食物鏈及其次生種的繁殖等，如許多生物對(duì)其聚居地最小面積的需求是其生存的條件之一；不同類型面積的大小能夠反映出其間物種、能量和養(yǎng)分等信息流的差異，一般來說，一個(gè)斑塊中能量和礦物養(yǎng)分的總量與其面積成正比；為了理解和管理景觀，我們往往需要了解斑塊的面積大小，如所需要的斑塊最小面積和最佳面積是極其重要的兩個(gè)數(shù)據(jù)。屬性Fragstats中的單位范圍景觀面積（TA）ha-公頃（10000平方米）TA0拼斑塊所占景觀面積的比例（%LAND）百分比0 %LAND=1最大斑塊所占景觀面積的比例（LPI）百

4、分比0LPI0面積加權(quán)的平均形狀因子（AWMSI）。平均最近距離（MNN）m米MNN0香農(nóng)多樣性指數(shù)（SHDI）無SHDI=0景觀分析LandscapeAnalysis第3頁，共54頁。下面列出一些重要的景觀分析指標(biāo)：斑塊所占景觀面積的比例（%LAND）公式描述：%LAND等于某一斑塊類型的總面積占整個(gè)景觀面積的百分比。其值趨于0時(shí)，說明景觀中此斑塊類型變得十分稀少；其值等于100時(shí)，說明整個(gè)景觀只由一類斑塊組成。生態(tài)意義： %LAND度量的是景觀的組分，其在斑塊級(jí)別上與斑塊相似度指標(biāo)（LSIM）的意義相同。由于它計(jì)算的是某一斑塊類型占整個(gè)景觀的面積的相對(duì)比例，因而是幫助我們確定景觀中模地（M

5、atrix）或優(yōu)勢(shì)景觀元素的依據(jù)之一；也是決定景觀中的生物多樣性、優(yōu)勢(shì)種和數(shù)量等生態(tài)系統(tǒng)指標(biāo)的重要因素。最大斑塊所占景觀面積的比例（LPI）公式描述： LPI等于某一斑塊類型中的最大斑塊占據(jù)整個(gè)景觀面積的比例。生態(tài)意義：有助于確定景觀的模地或優(yōu)勢(shì)類型等。其值的大小決定著景觀中的優(yōu)勢(shì)種、內(nèi)部種的豐度等生態(tài)特征；其值的變化可以改變干擾的強(qiáng)度和頻率，反映人類活動(dòng)的方向和強(qiáng)弱。斑塊平均大小（MPS）公式描述： MPS在斑塊級(jí)別上等于某一斑塊類型的總面積除以該類型的斑塊數(shù)目；在景觀級(jí)別上等于景觀總面積除以各個(gè)類型的斑塊總數(shù)。生態(tài)意義： MPS代表一種平均狀況，在景觀結(jié)構(gòu)分析中反映兩方面的意義：景觀中MP

6、S值的分布區(qū)間對(duì)圖像或地圖的范圍以及對(duì)景觀中最小斑塊粒徑的選取有制約作用；另一方面MPS可以指征景觀的破碎程度，如我們認(rèn)為在景觀級(jí)別上一個(gè)具有較小MPS值的景觀比一個(gè)具有較大MPS值的景觀更破碎，同樣在斑塊級(jí)別上，一個(gè)具有較小MPS值的斑塊類型比一個(gè)具有較大MPS值的斑塊類型更破碎。研究發(fā)現(xiàn)MPS值的變化能反饋更豐富的景觀生態(tài)信息，它是反映景觀異質(zhì)性的關(guān)鍵。景觀分析LandscapeAnalysis第4頁，共54頁。下面列出一些重要的景觀分析指標(biāo)：面積加權(quán)的平均形狀因子（AWMSI）公式描述： AWMSI在斑塊級(jí)別上等于某斑塊類型中各個(gè)斑塊的周長與面積比乘以各自的面積權(quán)重之后的和；在景觀級(jí)別上

7、等于各斑塊類型的平均形狀因子乘以類型斑塊面積占景觀面積的權(quán)重之后的和。其中系數(shù)0.25是由柵格的基本形狀為正方形的定義確定的。公式表明面積大的斑塊比面積小的斑塊具有更大的權(quán)重。當(dāng)AWMSI=1時(shí)說明所有的斑塊形狀為最簡單的方形（采用矢量版本的公式時(shí)為圓形）；當(dāng)AWMSI值增大時(shí)說明斑塊形狀變得更復(fù)雜，更不規(guī)則。生態(tài)意義： AWMSI是度量景觀空間格局復(fù)雜性的重要指標(biāo)之一，并對(duì)許多生態(tài)過程都有影響。如斑塊的形狀影響動(dòng)物的遷移、覓食等活動(dòng)，影響植物的種植與生產(chǎn)效率；對(duì)于自然斑塊或自然景觀的形狀分析還有另一個(gè)很顯著的生態(tài)意義，即常說的邊緣效應(yīng)。平均最近距離（MNN）公式描述： MNN在斑塊級(jí)別上等于

8、從斑塊ij到同類型的斑塊的最近距離之和除以具有最近距離的斑塊總數(shù)；MNN在景觀級(jí)別上等于所有類型在斑塊級(jí)別上的MNN之和除以景觀中具有最近距離的斑塊總數(shù)。生態(tài)意義： MNN度量景觀的空間格局。一般來說MNN值大，反映出同類型斑塊間相隔距離遠(yuǎn)，分布較離散；反之，說明同類型斑塊間相距近，呈團(tuán)聚分布。另外，斑塊間距離的遠(yuǎn)近對(duì)干擾很有影響，如距離近，相互間容易發(fā)生干擾；而距離遠(yuǎn)，相互干擾就少。但景觀級(jí)別上的MNN在斑塊類型較少時(shí)應(yīng)慎用。景觀分析LandscapeAnalysis第5頁，共54頁。下面列出一些重要的景觀分析指標(biāo)：香農(nóng)多樣性指數(shù)（SHDI）公式描述： SHDI在景觀級(jí)別上等于各斑塊類型的面

9、積比乘以其值的自然對(duì)數(shù)之后的和的負(fù)值。SHDI=0表明整個(gè)景觀僅由一個(gè)斑塊組成；SHDI增大，說明斑塊類型增加或各斑塊類型在景觀中呈均衡化趨勢(shì)分布。生態(tài)意義： SHDI是一種基于信息理論的測量指數(shù)，在生態(tài)學(xué)中應(yīng)用很廣泛。該指標(biāo)能反映景觀異質(zhì)性，特別對(duì)景觀中各斑塊類型非均衡分布狀況較為敏感，即強(qiáng)調(diào)稀有斑塊類型對(duì)信息的貢獻(xiàn)，這也是與其它多樣性指數(shù)不同之處。在比較和分析不同景觀或同一景觀不同時(shí)期的多樣性與異質(zhì)性變化時(shí)，SHDI也是一個(gè)敏感指標(biāo)。如在一個(gè)景觀系統(tǒng)中，土地利用越豐富，破碎化程度越高，其不定性的信息含量也越大，計(jì)算出的SHDI值也就越高。景觀生態(tài)學(xué)中的多樣性與生態(tài)學(xué)中的物種多樣性有緊密的聯(lián)

10、系，但并不是簡單的正比關(guān)系，研究發(fā)現(xiàn)在一景觀中二者的關(guān)系一般呈正態(tài)分布。香農(nóng)均度指數(shù)（SHEI）、散布與并列指數(shù)（IJI）、蔓延度指數(shù)（CONTAG）。都是Patch、Class、Landscape這三個(gè)層次上的計(jì)算得到的景觀分析LandscapeAnalysis第6頁，共54頁。通過上述關(guān)于層次與指標(biāo)概念介紹我們給景觀分析下一個(gè)定義景觀生態(tài)學(xué)是研究景觀單元的類型組成、空間配置及其與生態(tài)學(xué)過程相互作用的綜合性學(xué)科。強(qiáng)調(diào)空間格局、生態(tài)學(xué)過程與尺度之間的相互作用是景觀生態(tài)學(xué)研究的核心所在。 Wikipedia景觀分析LandscapeAnalysis第7頁，共54頁。由于計(jì)算量比較大，所以只截取部

11、分區(qū)域操作景觀分析LandscapeAnalysis行列數(shù)柵格大小圖像色深GRID格式第8頁，共54頁。設(shè)置：柵格大小、行列數(shù)、背景等如果是GRID就自動(dòng)檢測無需輸入選擇：ArcGrid圖像數(shù)據(jù)（或其它）（ASCII和img格式使用的區(qū)別）Standard分析方式（也可以用移動(dòng)窗口的方式）自動(dòng)保存結(jié)果文件景觀分析LandscapeAnalysis選擇GRID格式文件選擇輸出文件自動(dòng)保存輸出文件先選擇處理單個(gè)文件（Landscape），批處理（batch）后面會(huì)講到選擇Standard，移動(dòng)窗口的方法詳細(xì)見軟件的幫助文檔第9頁，共54頁。選擇：ArcGrid圖像數(shù)據(jù)（或其它）（ASCII和img

12、格式使用的區(qū)別）Standard分析方式（也可以用移動(dòng)窗口的方式）輸出數(shù)據(jù)選上所有級(jí)別（單個(gè)斑塊、斑塊類型、整體景觀）自動(dòng)保存結(jié)果文件為了對(duì)應(yīng)ArcGIS中的柵格轉(zhuǎn)矢量的方法，設(shè)置為4cell rule（附：如果設(shè)為8，則若在arcgis中把原圖轉(zhuǎn)成多邊形來與Fragstats生成的結(jié)果來匹配會(huì)出現(xiàn)斑塊數(shù)不同的狀況）景觀分析LandscapeAnalysis（4cell rule表示有共鄰邊才歸為一個(gè)斑塊，而8cell表示rule 即使對(duì)角相鄰也歸為一個(gè)斑塊）輸出與原圖位置相同，像元值為該像元所在斑塊的斑塊序號(hào) 的柵格圖設(shè)置你想要得到的景觀參數(shù)的級(jí)別設(shè)置分類文件，暫不用，后面會(huì)講最后，OK第1

13、0頁，共54頁。景觀分析LandscapeAnalysis最后，OK會(huì)發(fā)現(xiàn)提示“cannot load class properties file”，就是沒有選擇設(shè)置分類信息的原因，但也不影響GRID格式的分析，因?yàn)镕ragstats會(huì)默認(rèn)自動(dòng)識(shí)別分類信息第11頁，共54頁。對(duì)三個(gè)級(jí)別的想要輸出的指標(biāo)量進(jìn)行選擇按需選擇Patch（斑塊）Class（斑塊類型）Fragstats軟件強(qiáng)大，可計(jì)算出50多個(gè)景觀指標(biāo)。這些指標(biāo)被分為三組級(jí)別：分別為Patch（斑塊/斑塊）反映單塊的結(jié)果特征Class（斑塊類型）反映屬于相同的塊的特征Landscape（整體景觀）反映整體的結(jié)構(gòu)特征三種級(jí)別逐步擴(kuò)大尺度，

14、高度相關(guān)景觀分析LandscapeAnalysisLandscape（整體景觀）第12頁，共54頁。執(zhí)行成功景觀分析LandscapeAnalysis設(shè)置patch景觀參數(shù)設(shè)置patch景觀參數(shù)設(shè)置patch景觀參數(shù)設(shè)置總體參數(shù)運(yùn)行相同效果ID 文件默認(rèn)保存到了與源文件同一個(gè)文件夾下第13頁，共54頁。如果沒有設(shè)置autosave的話，點(diǎn)Save run as 保存結(jié)果文件執(zhí)行成功，查看結(jié)果景觀分析LandscapeAnalysis第14頁，共54頁。得到的文件：四個(gè)結(jié)果文件、一個(gè)柵格文件（GRID格式）景觀分析LandscapeAnalysis用記事本打開生成的文件注意查看幫助文檔說明的單位

15、Adj文件是各類型統(tǒng)計(jì)的矩陣class文件是每種類型的對(duì)應(yīng)參數(shù)結(jié)果，這里為class級(jí)別上的斑塊密度（PD：patch density）第15頁，共54頁。Land文件，得到形狀指數(shù)LSI和景觀聚集度指數(shù)con，而且都是MN（平均值）Patch文件，因?yàn)榘邏K數(shù)數(shù)量多，所以文件也稍大，這里算出面積，周長和回旋半徑景觀分析LandscapeAnalysis用記事本打開生成的文件注意查看幫助文檔說明的單位第16頁，共54頁。此外還會(huì)生成一個(gè)柵格圖，其像元值為其斑塊在Fragstats中的編號(hào)順序ID Image圖像ID Image圖像的屬性表這個(gè)文件名字是原圖名字后面加上id4，如果是8cell r

16、ule的話會(huì)是id8（Fragstats自動(dòng)加的）景觀分析LandscapeAnalysis注意：必須前面有選上這個(gè)才有這個(gè)圖像第17頁，共54頁。導(dǎo)入到ArcGIS中用Excel打開（按逗號(hào)”,”分割）ArcCatalog中新建個(gè)人數(shù)據(jù)庫并導(dǎo)入single tableArcGIS中在對(duì)應(yīng)layer中join導(dǎo)入的屬性表景觀分析LandscapeAnalysis一定要是“打開”而不能是直接拖進(jìn)去選擇patch文件第18頁，共54頁。導(dǎo)入到ArcGIS中景觀分析LandscapeAnalysis用Excel打開（按逗號(hào)”,”分割）ArcCatalog中新建個(gè)人數(shù)據(jù)庫并導(dǎo)入single table

17、ArcGIS中在對(duì)應(yīng)layer中join導(dǎo)入的屬性表選 分隔符然后下一步第19頁，共54頁。導(dǎo)入到ArcGIS中景觀分析LandscapeAnalysis用Excel打開（按逗號(hào)”,”分割）ArcCatalog中新建個(gè)人數(shù)據(jù)庫并導(dǎo)入single tableArcGIS中在對(duì)應(yīng)layer中join導(dǎo)入的屬性表選擇 逗號(hào)下一步第20頁，共54頁。導(dǎo)入到ArcGIS中景觀分析LandscapeAnalysis用Excel打開（按逗號(hào)”,”分割）ArcCatalog中新建個(gè)人數(shù)據(jù)庫并導(dǎo)入single tableArcGIS中在對(duì)應(yīng)layer中join導(dǎo)入的屬性表默認(rèn)，點(diǎn)完成第21頁，共54頁。導(dǎo)入到A

18、rcGIS中景觀分析LandscapeAnalysis用Excel打開（按逗號(hào)”,”分割）ArcCatalog中新建個(gè)人數(shù)據(jù)庫并導(dǎo)入single tableArcGIS中在對(duì)應(yīng)layer中join導(dǎo)入的屬性表顯示有些奇怪，需要去掉表格中的空格第22頁，共54頁。導(dǎo)入到ArcGIS中景觀分析LandscapeAnalysis用Excel打開（按逗號(hào)”,”分割）ArcCatalog中新建個(gè)人數(shù)據(jù)庫并導(dǎo)入single tableArcGIS中在對(duì)應(yīng)layer中join導(dǎo)入的屬性表界面顯示清楚了替換空格的前提是你的文件目錄路徑不包含空格第23頁，共54頁。導(dǎo)入到ArcGIS中景觀分析Landscape

19、Analysis用Excel打開（按逗號(hào)”,”分割）ArcCatalog中新建個(gè)人數(shù)據(jù)庫并導(dǎo)入single tableArcGIS中在對(duì)應(yīng)layer中join導(dǎo)入的屬性表因?yàn)檫@次都是同一個(gè)文件的，可刪除第一列的路徑信息。若是多個(gè)文件則需按第一行來分類第24頁，共54頁。導(dǎo)入到ArcGIS中景觀分析LandscapeAnalysis用Excel打開（按逗號(hào)”,”分割）ArcCatalog中新建個(gè)人數(shù)據(jù)庫并導(dǎo)入single tableArcGIS中在對(duì)應(yīng)layer中join導(dǎo)入的屬性表另存為xls表格文件第25頁，共54頁。導(dǎo)入到ArcGIS中景觀分析LandscapeAnalysis之所以要新建

20、一個(gè)，主要是為了到時(shí)能導(dǎo)出一個(gè)數(shù)據(jù)庫文件，如果直接join excel表格的話，當(dāng)excel位置變化后會(huì)丟失join用Excel打開（按逗號(hào)”,”分割）ArcCatalog中新建個(gè)人數(shù)據(jù)庫并導(dǎo)入single tableArcGIS中在對(duì)應(yīng)layer中join導(dǎo)入的屬性表第26頁，共54頁。導(dǎo)入到ArcGIS中景觀分析LandscapeAnalysis用Excel打開（按逗號(hào)”,”分割）ArcCatalog中新建個(gè)人數(shù)據(jù)庫并導(dǎo)入single tableArcGIS中在對(duì)應(yīng)layer中join導(dǎo)入的屬性表為新建的個(gè)人數(shù)據(jù)庫導(dǎo)入剛才保存的xls表格文件第27頁，共54頁。導(dǎo)入到ArcGIS中景觀分析

21、LandscapeAnalysis選擇那個(gè)表格，雙擊進(jìn)去，選中$結(jié)尾的文件名（其實(shí)是excel中的一個(gè)表格頁）用Excel打開（按逗號(hào)”,”分割）ArcCatalog中新建個(gè)人數(shù)據(jù)庫并導(dǎo)入single tableArcGIS中在對(duì)應(yīng)layer中join導(dǎo)入的屬性表第28頁，共54頁。導(dǎo)入到ArcGIS中景觀分析LandscapeAnalysis用Excel打開（按逗號(hào)”,”分割）ArcCatalog中新建個(gè)人數(shù)據(jù)庫并導(dǎo)入single tableArcGIS中在對(duì)應(yīng)layer中join導(dǎo)入的屬性表再設(shè)置一下輸出名字，就OK第29頁，共54頁。導(dǎo)入到ArcGIS中景觀分析LandscapeAnal

22、ysis用Excel打開（按逗號(hào)”,”分割）ArcCatalog中新建個(gè)人數(shù)據(jù)庫并導(dǎo)入single tableArcGIS中在對(duì)應(yīng)layer中join導(dǎo)入的屬性表把剛才的ID柵格文件和new personal database中的表格拖到ArcGIS中第30頁，共54頁。導(dǎo)入到ArcGIS中景觀分析LandscapeAnalysis用Excel打開（按逗號(hào)”,”分割）ArcCatalog中新建個(gè)人數(shù)據(jù)庫并導(dǎo)入single tableArcGIS中在對(duì)應(yīng)layer中join導(dǎo)入的屬性表Join表格柵格的value對(duì)應(yīng)表格的ObjectID第31頁，共54頁。導(dǎo)入到ArcGIS中景觀分析Lands

23、capeAnalysis用Excel打開（按逗號(hào)”,”分割）ArcCatalog中新建個(gè)人數(shù)據(jù)庫并導(dǎo)入single tableArcGIS中在對(duì)應(yīng)layer中join導(dǎo)入的屬性表Join表格后第32頁，共54頁。ArcGIS中顯示舉例原土地利用圖回旋半徑圖景觀分析LandscapeAnalysis第33頁，共54頁。ArcGIS中顯示舉例原土地利用圖回旋半徑圖景觀分析LandscapeAnalysisFragstats中幫助文檔的描述第34頁，共54頁。ArcGIS中顯示舉例類似的，可以對(duì)class文件與源圖像按照類型進(jìn)行匹配來join到圖像的屬性表中景觀分析LandscapeAnalysis

24、第35頁，共54頁。ASCII類型ArcGIS中將柵格轉(zhuǎn)換為ASCII后是不能直接在fragstats分析的，需要打開txt文件，把前面的幾行描述信息刪掉，只保留柵格信息，才能被Fragstats分析。須刪掉前面那6行信息ASCII類型景觀分析LandscapeAnalysis第36頁，共54頁。ERDAS和8bit類型貌似Fragstats不支持ERDAS9以后的img格式，所以會(huì)報(bào)錯(cuò)ERDAS類型其幫助文檔是這樣描述的：8bit類型雖然不可以直接用ERDAS類型，但可以選擇8bit類型來對(duì)img圖像進(jìn)行分析（實(shí)際就可以是ERDAS9的文件）但是這時(shí)候一定要選擇分類文件，否則即使能夠執(zhí)行，但

25、結(jié)果也是不正確的景觀分析LandscapeAnalysis第37頁，共54頁。指定Class類型文件這里可以指定你自己寫的txt或文件或fdc文件，其格式很簡單ClassID, ClassName, Status, isBackgroundStatus 代表是否參與計(jì)算， isBackground代表是否作為背景景觀分析LandscapeAnalysis第38頁，共54頁。批處理如果要對(duì)多個(gè)鎮(zhèn)區(qū)分別計(jì)算其景觀指數(shù)的話，需要在Fragstats用到批處理，格式如下：景觀分析LandscapeAnalysisGRID文件的批處理格式Fragstats幫助文檔的關(guān)于其他文件格式的批處理格式，基本類似

26、第39頁，共54頁。批處理會(huì)遇到的問題Fragstats不能對(duì)一整幅圖像自動(dòng)分各鎮(zhèn)區(qū)來計(jì)算對(duì)應(yīng)的景觀指數(shù)，所以先要把各鎮(zhèn)區(qū)的土地利用圖裁剪出來景觀分析LandscapeAnalysis批量裁剪柵格圖的方法：1、split工具分割矢量圖（這里要注意有些鎮(zhèn)區(qū)包含多個(gè)polygon時(shí)，因?yàn)閮蓚€(gè)polygon名字相同，會(huì)導(dǎo)致分割出9999的錯(cuò)誤，要修改其名字為為唯一值才可以，如在鎮(zhèn)區(qū)名前加上FID，因?yàn)镕ID肯定是唯一的）2、右鍵點(diǎn)擊extract by mask工具來批量裁剪柵格圖（這個(gè)需時(shí)較長）3、在excel中批量生成Fragstats中需要的文件名注意：文件路徑都不能有空格，否則Fragsta

27、ts分析時(shí)會(huì)報(bào)錯(cuò)第40頁，共54頁。批處理會(huì)遇到的問題Fragstats中的問題景觀分析LandscapeAnalysisFragstats批處理的方法和問題：1、在excel中批量生成Fragstats中需要的文件名（手動(dòng)輸入太麻煩）2、Fragstats對(duì)批量的柵格圖生成了結(jié)果文件后，以城市用地的斑塊密度為例，有些鎮(zhèn)區(qū)包含城市用地和非城市用地，而有些鎮(zhèn)區(qū)卻只有城市用地或只有非城市用地，這時(shí)下圖的情況只有城市用地這一行有些polygon的面積比柵格圖的一個(gè)像元還小，造成裁剪出來也是無實(shí)際值，F(xiàn)ragstats就無法識(shí)別，變成沒有這個(gè)鎮(zhèn)區(qū)一樣通常我們需要的只是城市用地的斑塊密度，所以要把這兩個(gè)

28、得到的值分離開來第41頁，共54頁。批處理會(huì)遇到的問題Fragstats中的問題景觀分析LandscapeAnalysis有些鎮(zhèn)區(qū)有一行，有些有兩行，有些連一行都沒有，這樣不規(guī)則的數(shù)據(jù)分離，可以用到excel中的宏命令來處理Excel宏命令步驟1、office2007為例子2、視圖3、查看宏4、輸入自定義的宏名字，點(diǎn)編輯第42頁，共54頁。批處理會(huì)遇到的問題景觀分析LandscapeAnalysis關(guān)鍵是找先到定位單位格位置的函數(shù)，見下圖Sub forclass()需要先清除excel內(nèi)所有空格Dim t As IntegerDim s As StringDim count As Intege

29、rDim nomalCity_FD As String 對(duì)應(yīng)LID的位置顯示PD值Dim orderCity_FD As String 去掉TYPE為0的重復(fù)城市Dim order As String 序號(hào)t = 1 累計(jì)城市數(shù)量s = D:split1995outputextract_b 指定文件名和位置count = 350 設(shè)定得到的文件數(shù)量設(shè)定輸出列-nomalCity_FD = forderCity_FD = horder = gRange(nomalCity_FD + 1).Value = PD(city)Range(orderCity_FD + 1).Value = PD(cit

30、y)Range(order + 1).Value = FID序號(hào)Range(order + 2).Value = 1For n = 2 To count -If Range(A + CStr(n).Value Like D:split1995outputextract_b + CStr(t) Then If Range(b + CStr(n).Value = 1 Then Range(nomalCity_FD + CStr(n).Value = Range(c + CStr(n).Value Range(orderCity_FD + CStr(t + 1).Value = Range(c +

31、CStr(n).Value Else Range(nomalCity_FD + CStr(n).Value = - Range(orderCity_FD + CStr(t + 1).Value = - End IfElse t = t + 1 累加一項(xiàng) Range(order + CStr(t + 1).Value = t If Range(b + CStr(n).Value = 1 Then Range(nomalCity_FD + CStr(n).Value = Range(c + CStr(n).Value Range(orderCity_FD + CStr(t + 1).Value =

32、 Range(c + CStr(n).Value Else Range(nomalCity_FD + CStr(n).Value = - Range(orderCity_FD + CStr(t + 1).Value = 0 End IfEnd IfNext n -End Sub第43頁，共54頁。批處理會(huì)遇到的問題景觀分析LandscapeAnalysis得到的分離的結(jié)果第44頁，共54頁。Moving window選擇Moving window景觀分析LandscapeAnalysis按需要選擇圓形（Round）或方形（Square）范圍概念：首先要指明的是，選擇moving window后

33、只能分析landscape級(jí)別的指數(shù)，因?yàn)閙oving window就是源圖像中每個(gè)像元都按指定范圍只作為一個(gè)整體來分析，并把分析的值保存在對(duì)應(yīng)點(diǎn)坐標(biāo)位置，形成新的圖像。每選一個(gè)landscape級(jí)別的指數(shù)就會(huì)多生成一個(gè)grid圖像第45頁，共54頁。Moving window景觀分析LandscapeAnalysis選擇兩個(gè)參數(shù)：shapeIndex（形狀指數(shù)）和ContiguityIndex（聚集度指數(shù)）第46頁，共54頁。Moving window景觀分析LandscapeAnalysisMoving window計(jì)算得到的文件與之前standard模式下的不同，它生成的是GRID格式文件，每個(gè)GRID圖像對(duì)應(yīng)一個(gè)參數(shù)：ContiguityIndex（聚集度指數(shù)）和shapeIndex（形狀指數(shù)）第47頁，共54頁。Moving window景觀分析LandscapeAnalysisMoving window得到的圖像的范圍會(huì)比源圖像范圍較小，因?yàn)楫?dāng)像元點(diǎn)離邊界的距離小于指定的范圍（radius）時(shí)，就會(huì)使得該像元點(diǎn)的指定范圍不能構(gòu)成完整的區(qū)域，見Fragstats的幫助文檔第48頁，共54頁。Moving window景觀分