转移矩阵计算

1、写文章的过程中，同样也用到了土地利用/景观转移矩阵。转移矩阵的制作很多软件可以实现，有Arcview、Arcmap、还有人用matlab。但是网上都写的都不详细，在这里，本文才用arcmap的叠加操作+excel的数据统计功能，可以轻松完成这一过程。相信看完了本文，初学者或者对GIS了解不多的人，也能轻松的制作出转移矩阵。方法一：1数据准备准备好两期的coverage或shape文件。(注意：拓扑关系要建好，不然无法计算)，采用Arcmap打开上述两期文件。2叠加操作选才AArcmap里面的"ArcToolbox”按钮下面的"AnalysisTools”工具下面的叠加分析模

2、块(Overlay)下面的交叉分析功能(Intersect)。选择已经打开的两期数据进行叠加分析。叠加的结果存在一个可以找到的地方。同时，把叠加后的结果添加在Arcmap里。，划ArcTOilS+31Analyst二i,国Tg工,+电Extract上C1VL电工1vIdentitjr?IntersectSymmetricalDl££«户VniQH产Upda.l«3输出叠加文件的属性数据A、右键打开intersect产生的矢量文件的数据属性表(openattributetable)。B、点击"options”按钮，选择"addfield

3、”,然后给出一个新name"newarea”,数据类型为doubleC、右键点击刚刚产生的"newarea",并选择"calculatevalues”D、然后点击“fieldcalculator"对话框里面的"Advanced"后选择"help"将DimOutputasdoubleDimpAreaasIareaSetpArea=shapeOutput=pArea.area拷入到“Reldcalculator"对话框下面的空白处E、在对话框“fieldcalculator"最下面的空白处

4、填上"output"Ficrlci4包工匚1u工曰七ormStk.l.p«DUPFIDj(itftlbBUD_I防看Fr*nLc£kcVtlKSrrr&tCo-d«Di»Outp-utmedubla口in>Ar««MItf*«5*tpAr«*s.hfeptOkklput-V*F、在属性表点击“options”按钮，一export（导出），属性以.dbf格式结果存储。4Excel进行转移矩阵制作A、刚才存储的.dbf文件可以使用excel打开，打开的结果重新保存为.xls（excel

5、）文件。（dbf文件不能保存）B、重新打开刚存储的excel文件。C、选中所有数据，选择数据数据透视表和数据透视图。将两期数据的id值分别拖入行列字段，然后以newarea'字段作为数据项拖入计算区域。得到的便是转移矩阵。如果发现矩阵的形式不美观，可以对单位进行修改。方法二：作某一地物与其他不同年份地类例如：1、做水体和其它年份分类地物的交集A.分别打开刚刚产生的水域多边形和该地区其他年份的分类矢量文件B.选才?gis里面的"AacToolbox”按钮下面的"AnalysisTools”工具下面的叠加分析模块（Overlay）下面的交叉分析功能（Intersect）

6、C.分别输入上面抽出的水域多边形和该地区其他年份区的分类后的多边形进行intersect分析D.右键打开intersect产生的矢量文件的数据属性表（openattributetable）E.点击“options”按钮，选择"addfield”,然后给出一个新name"newarea”,数据类型为doubleF.右键点击刚刚产生的"newarea",并选择"calculatevalues”G.然后点击“fieldcalculator"对话框里面的"Advanced"后选择"help"将DimOu

7、tputasdoubleDimpAreaasIareaSetpArea=shapeOutput=pArea.area拷入到“fieldcalculator"对话框下面的空白处H.在对话框“fieldcalculator"最下面填上"output"2、然后分别计算水体在该地区其他年份年转入和转出情况A.点击“potions”下面的"selectbyattribute"按钮，并选择"GRID_COD_1B.令GRID_COD_1=1点击apply,此时会统计出该地区其他年份年仍然为水域的多边形C.右键点击“newarea”选择统

8、计功能（“statistics"）,就会得到从然为水域的面积D.分别计算出GRID_COD_1=2（林地）、3（草地）、4（城镇点）、5（耕地）、6（滩涂地）的面积就得到从由水域转入的面积。以此类推网上有人还采用的方法：1、 Erdas：erdas-main-interpreter-gisanalysis-matrix输入两个年份的分类图，然后进行重编码即可。2、 Arcview:可以作个运算，其中一期的图*10或者100或者1000（土地利用类型取1位数时彳炭X10,两位100），然后再加上另一期的土地利用图得出一个图的属性，肯定12,15,或者1221,或者1531之类的数据，应

9、该就可以实现了。3、 Matlab:听说，还没有亲眼见过，据说可以成功4、ArcGIS:在一期数据的基础上，绘制出其变化的部分（很复杂，工作了量很大）。看了这些，我想本文所采用的方法，应该是最为简单的一种。如果有什么问题，可以给我留言一起探讨。同时借此文感谢小曾曾(82444233)、嘉许(330323039)转移矩阵在土地利用变更调查时经常用到。arcgis里没有直接的工具来计算转移矩阵，但我们可以利用已有的工具间接来求。下面介绍一下利用arcgis求转移矩阵的方法。现有2005、2008两期土地利用现状数据(polygonshp)：在求转移矩阵前先做好下列工作：1、备份数据。2、分别将两期数据用dissoke工具按各自的地类编码合并要素，注意选中最下面的“createmultipatfeatures选项处理好数据后，用intersect工具求两期数据的交集交集的结果：将属性表增加面积字段，并计算面积结果：从上图中已经可以看出土地转移的情况：dlbm2005是2005年时该块地的地类代码，dlbm2008是2008年时该地块的地类代码，最后一列是相应地块的面积。接下来将属性表导出到dbf文件，然后在excel中打开，执行数据菜单中的数据透视表和数据透视图”命令点完成，然后将""dlbm2005作为行字段，dlbm2008作为列