土地利用轉(zhuǎn)移矩陣生成的幾種方法

1、土地利用轉(zhuǎn)移矩陣生成的幾種方法土地利用類(lèi)型轉(zhuǎn)換矩陣為相同區(qū)域在不同時(shí)間段內(nèi)的土地利用類(lèi)型的相互轉(zhuǎn)換關(guān)系，一般用二維表來(lái)表達(dá)，從二維表中可以快速查看各個(gè)地類(lèi)間相互轉(zhuǎn)化的具體情況。比如某一類(lèi)別的土地有百分之多少（或者面積）分別轉(zhuǎn)化成了其他的土地類(lèi)型，現(xiàn)在某類(lèi)型的土地分別是由過(guò)去的哪些類(lèi)別轉(zhuǎn)化而來(lái)的等等。還可以生成變化統(tǒng)計(jì)柵格圖（掩膜圖像），它描述了前后兩幅土地分類(lèi)圖之間的地類(lèi)發(fā)生轉(zhuǎn)變的位置和類(lèi)別。土地利用類(lèi)型轉(zhuǎn)換矩陣可以從兩幅柵格圖中計(jì)算得到，也可以從兩個(gè)矢量文件中計(jì)算獲得。下面介紹土地利用轉(zhuǎn)移矩陣生成的幾種方法：ENVI下從兩幅分類(lèi)結(jié)果的柵格圖中計(jì)算土地利用類(lèi)型轉(zhuǎn)換矩陣1、準(zhǔn)備數(shù)據(jù)兩個(gè)時(shí)相的土地

2、利用分類(lèi)結(jié)果，它是單波段、專(zhuān)題類(lèi)型的偽彩色圖像（ENVIClassification）2、計(jì)算轉(zhuǎn)換矩陣打開(kāi)兩個(gè)土地利用分類(lèi)結(jié)果。（1）在主菜單中，選擇BasicToolsChangeDetectionChangeDetectionStatistics（2）分別在InitialState對(duì)話(huà)框和finalstate對(duì)話(huà)框中選擇前一時(shí)相和后一時(shí)相的土地利用結(jié)果。（3）在DefineEquivalentClasses對(duì)話(huà)框中（圖1）,如果兩個(gè)土地利用分類(lèi)名稱(chēng)一致，系統(tǒng)自動(dòng)將InitialStateClass和FinalStateClass對(duì)應(yīng)，否則手動(dòng)選擇，單擊AddPair按鈕選擇。圖1Defin

3、eEquivalentClasses對(duì)話(huà)框（4）選擇對(duì)應(yīng)的地物類(lèi)型之后，單擊0K按鈕，出現(xiàn)圖2對(duì)話(huà)框。選擇生成圖表表示單位（ReportType）：像素（Pixels）、百分比（Percent）和面積（Area）。選擇OutputClassificationMaskImages?為YES,輸出掩膜圖像，選擇輸入路徑及文件名。圖2選擇數(shù)據(jù)參數(shù)(5)單擊OK,執(zhí)行土地利用類(lèi)型轉(zhuǎn)換矩陣計(jì)算過(guò)程。3、查看結(jié)果(1)如圖3為得到的土地利用類(lèi)型轉(zhuǎn)換矩陣結(jié)果。橫字段表示前一時(shí)間段(InitialState)的土地利用類(lèi)別，縱字段為后一時(shí)間段(FinalState)的土地利用類(lèi)別。橫字段和縱字段交叉處表示變化

4、值，如有2520900平方米林地用地變化為草地。圖3土地利用類(lèi)型轉(zhuǎn)換矩陣(2)還可以為每一個(gè)地類(lèi)生成一個(gè)變換掩膜圖像，圖4所示為其中一個(gè)地類(lèi)的掩膜圖像。掩膜圖像的灰度值表示變化類(lèi)型，如這里的2草地表示林地變化為草地的像元。圖4變化掩膜圖像來(lái)源若你的數(shù)據(jù)是Raster格式：則有如下方法1ErdasImagine-lnterpreter-GisAnalysis-Matrix,輸入兩個(gè)時(shí)相的Raster數(shù)據(jù)即可做這一步之前記得先對(duì)兩時(shí)相的數(shù)據(jù)進(jìn)行重編碼(nterpreter-GisAnalysis-Recode)一般運(yùn)行如果出現(xiàn)錯(cuò)誤肯定是重編碼沒(méi)做好,請(qǐng)繼續(xù)查證。2先在Erdas中利用Modeler

5、計(jì)算如下公式NC(I,J)=NC(I)*10+NC(J),(JI)其中：NC(I,J)表示i,j兩年份的土地利用變化圖；NC(i)表示i年份遙感分類(lèi)影像；NC(j)表示j年份的遙感分類(lèi)影像。在此計(jì)算的基礎(chǔ)上，將以上變化影像圖轉(zhuǎn)化為BIL格式，再利用ARC/INFOGRID模塊將影像轉(zhuǎn)為GRID格式，然后利用GRID模塊中的屬性表(vat)查看命令對(duì)影像灰度值進(jìn)行統(tǒng)計(jì)，最后得出土地利用轉(zhuǎn)化舉證。(注:此方法本人尚未實(shí)現(xiàn)過(guò),不知可行否)。若數(shù)據(jù)是Vector格式ErdasImagine-lnterpreter-GisAnalysisMatrix,輸入兩個(gè)時(shí)相的Vector數(shù)據(jù)即可此時(shí)注意輸出柵格大

6、小不應(yīng)設(shè)的太小要不一運(yùn)行就會(huì)提示你的空間不足做這一步之前，請(qǐng)做好前期的地理編碼。ArcView3.3加載spatialanalysis模塊把兩時(shí)相的Vector圖轉(zhuǎn)成grid格式(當(dāng)然中間有一些單位的設(shè)置根據(jù)你做的圖的分辨率來(lái)設(shè)置即可)analysis-mapcaculate直接計(jì)算即可。3把兩期解譯完的Vector文件在arctoolboxoverlayunion中疊加，注意：兩個(gè)文件不能用同一個(gè)字段名，比如一個(gè)用93Type,另一個(gè)時(shí)相則用OOType疊加后的文件在Arcmap中打開(kāi)，選中文件，然后點(diǎn)右Property空間查詢(xún)，輸入條件語(yǔ)句，比如：93Type=1And00Type=2；查

7、詢(xún)結(jié)果即為第一種類(lèi)型轉(zhuǎn)化為第二種類(lèi)型的圖形，可以另建一圖層比如：12，把查詢(xún)結(jié)果復(fù)制到12圖層上。統(tǒng)計(jì)出面積，依進(jìn)行，就可以得到土地利用類(lèi)型轉(zhuǎn)移矩陣。一、數(shù)據(jù)準(zhǔn)備(圖1)準(zhǔn)備兩幅不同時(shí)相的土地利用現(xiàn)狀圖(shp格式)，每幅圖的屬性表都要有一個(gè)表示土地利用類(lèi)型的字段，并且要使用不同的名稱(chēng)加以區(qū)分，女Type1995Type2000。土地利用類(lèi)型名稱(chēng)必須統(tǒng)一，并且完整，如都使用“城鎮(zhèn)用地”、“有林地”等。、數(shù)據(jù)融合（圖2）在ArcMap里分別打開(kāi)兩個(gè)時(shí)相的圖層，打開(kāi)ArcToolbox，選擇DataManagementTools|Generalization|Dissolve工具。InputFea

8、ture選擇要融合的圖層，OutputFeatureClass選擇輸出結(jié)果存儲(chǔ)的位置及名稱(chēng)，DissolveField（s）選擇土地利用類(lèi)型字段（如Type1995）,然后勾選Creatmultipartfeatures選項(xiàng)，點(diǎn)擊OK完成。重復(fù)此過(guò)程，對(duì)另一時(shí)相數(shù)據(jù)進(jìn)行融合。此步驟使相同利用類(lèi)型的記錄融合為一個(gè)記錄，以提高后面步驟的計(jì)算速度。4*SL扎+厲I14*13EiHI口：r百4冷伸“町說(shuō)UE3斬國(guó);侮*tu-*口1i*L-rtf-Mi-iid*v母i”證口訕Ei.4-i.i24-.J1“乜Ifim三、疊置分析（圖3）在ArcMap中打開(kāi)兩個(gè)時(shí)相融合后的數(shù)據(jù)，在ArcToolbox中選擇

9、AnalysisTools|Overlay|Intersect工具，InputFeatures選擇兩個(gè)時(shí)相的圖層，OutputFeatureClass選擇疊加結(jié)果存儲(chǔ)的位置及名稱(chēng)，其余選項(xiàng)可以忽略，單擊【0K】完成。E.Lb7*乂:1H.14.ue-LLH卄IMl-liT即口ii廠(chǎng)訶siZi叭！T|-ii|mJi*!1MtwihflT-rfSi3-4四、計(jì)算面積并導(dǎo)出屬性表（圖4-6）在ArcMap中打開(kāi)疊加后的圖層數(shù)據(jù)，在該圖層上右鍵打開(kāi)屬性表，選擇Option|Addfield新建一個(gè)字段，命名為NewArea。44Q如仏liAiiuQXs曲L.Pri卜4站Mai(4bru.i*Irfal

10、HImJhb腫|JDcfQ命):FK嚇、if鼻騎-rpc!S|S|5*孕廠(chǎng)jd-*Atri1:一在Editer工具條中選擇Editer|StartEditing然后在屬性表中NewArea字段上單擊右鍵選擇CalculateGeometry.，在打開(kāi)的CalculateGeometry對(duì)話(huà)框中，Property選擇Area，Units選擇要使用的面積單位，單擊【0K】完成圖斑面積計(jì)算。依次選擇Editer|SaveEdits/EndEditing保存和退出編輯狀態(tài)。$LiflUilinj：+tilt1JDer念：iMT,t*.:段?？趇Vr+lillLlll-rkliaiLLjrLl*SJIt

11、11kl;hl護(hù)I.L第一步書(shū):H.匚時(shí)F&lreSeeL牛THHlhbcEvil.fl在屬性表中選擇Option|Export將屬性表保存為dbf文件。oajkftWi一-F5A,*!ASaliriLI3R-l屮Ml(4:hi.03jf51tJiltdBllHSfabMldHCj.訴礙JkiU-iM.l*La4rifXk|Q缶ZEfii*0g.R.名hl-1+41*4AE土二*五、制作轉(zhuǎn)移矩陣（圖7-10）（以Excel2007為例）在Excel中打開(kāi)上一步保存的dbf,另存為Excel格式并打開(kāi)。在Excel中選中所有數(shù)據(jù)（不要點(diǎn)左上角，只選擇有效數(shù)據(jù)），點(diǎn)擊【插入】選項(xiàng)卡，選擇【數(shù)據(jù)透視表】|【數(shù)據(jù)透視表】，點(diǎn)擊【確定】。A-二*L-pcdtGatfl-Wk!OKfl臼疋召*汐幵葷1?.貝辭斥筑蹲1昭廻肖塩2125270&辛工亦秦CHE3匸MT燒iKS出田TflE)朮S3炸和iKffi*a)OftfflwrwwrtijiiJVVRbVVO-0WO62.OCOQl?ZJ03932廻片SF4*誠(chéng)如在打開(kāi)的數(shù)據(jù)透視表中按圖示將字段拖入相應(yīng)區(qū)域。Excel自動(dòng)計(jì)算矩陣，將該表稍事整飾就得到美觀的土地利用轉(zhuǎn)移矩陣。矩陣中r(l,j)就表示i類(lèi)型向j類(lèi)型轉(zhuǎn)移的土地面積，空