景观转移矩阵的制作

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

2、erlay）下面的交叉分析功能（Intersect）。选择已经打开的两期数据进行叠加分析。叠加的结果存在一个可以找到的地方。同时，把叠加后的结果添加在Arcmap里。 3 输出叠加文件的属性数据 A、右键打开intersect产生的矢量文件的数据属性表（open attribute table）。 B、点击“options”按钮，选择“add field”，然后给出一个新name“newarea”，数据类型为double C、右键点击刚刚产生的“newarea”，并选择“calculate values” D、然后点击“field calculator”对话框里面的“Advanced”后选择“

3、help”将 Dim Output as double Dim pArea as Iarea Set pArea = shape Output = pArea.area 拷入到“field calculator”对话框下面的空白处 E、在对话框“field calculator”最下面的空白处填上“output” F、在属性表点击“options”按钮，export（导出），属性以.dbf格式结果存储。 4Excel进行转移矩阵制作 A、刚才存储的.dbf文件可以使用excel打开，打开的结果重新保存为.xls（excel）文件。 （dbf文件不能保存） B、重新打开刚存储的excel文件。

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

5、ct分析 D.右键打开intersect产生的矢量文件的数据属性表（open attribute table） E.点击“options”按钮，选择“add field”，然后给出一个新name“newarea”，数据类型为double F.右键点击刚刚产生的“newarea”，并选择“calculate values” G.然后点击“field calculator”对话框里面的“Advanced”后选择“help”将Dim Output as double Dim pArea as Iarea Set pArea = shape Output = pArea.area拷入到“field c

6、alculator”对话框下面的空白处 H.在对话框“field calculator”最下面填上“output” 2、然后分别计算水体在该地区其他年份年转入和转出情况 A.点击“potions”下面的“select by attribute”按钮，并选择“GRID_COD_1” B.令GRID_COD_11 点击apply，此时会统计出该地区其他年份年仍然为水域的多边形 C.右键点击“newarea”选择统计功能（“statistics”），就会得到从然为水域的面积 D.分别计算出GRID_COD_12（林地）、3（草地）、4（城镇点）、5（耕地）、6（滩涂地）的面积就得到从由水域转入的面积。 以此类推 网上有人还采用的方法： 1、Erdas：erdas-main-interpreter-gis analysis-matrix输入两个年份的分类图，然后进行重编码即可。 2、Arcview：可以作个运算，其中一期的图*10或者100或者1000（土地利用类型取1位数时候10，两位100),然后再加上另一期的土地利用图得出一个图的属性，肯定