最简单的土地利用景观转移矩阵的制作
方法一:
1数据准备
准备好两期的coverage或shape文件。(注意:拓扑关系要建好，不然无法计算)，采用Arcmap打开上述两期文件。
2叠加操作
选择Arcmap里面的“ArcToolbox”按钮下面的“Analysis Tools”工具下面的叠加分析模块(Overlay)下面的交叉分析功能(Intersect)。选择已经打开的两期数据进行叠加分析。叠加的结果存在一个可以找到的地方。同时，把叠加后的结果添加在Arcmap里。
3输出叠加文件的属性数据
A、右键打开intersect产生的矢量文件的数据属性
表
关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf
(open attribute table)。 B、点击“options”按钮，选择“add field”，然后给出一个新name“newarea”，数据类型为double C、右键点击刚刚产生的“newarea”，并选择“calculate values”
D、然后点击“field calculator”对话框里面的“Advanced”后选择“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格式结果存储。 4 Excel进行转移矩阵制作
A、刚才存储的.dbf文件可以使用excel打开，打开的结果重新保存为.xls(excel)文件。 (dbf文件不能保存)
B、重新打开刚存储的excel文件。
C、选中所有数据，选择数据\数据透视表和数据透视图。将两期数据的id值分别拖入行列字段，然后以‘newarea’字段作为数据项拖入计算区域。得到的便是转移矩阵。如果发现矩阵的形式不美观，可以对单位进行修改。
方法二:作某一地物与其他不同年份地类
例如:
1、做水体和其它年份分类地物的交集
A. 分别打开刚刚产生的水域多边形和该地区其他年份的分类矢量文件 B. 选择gis里面的“AacToolbox”按钮下面的“Analysis Tools”工具下面的叠加分析模块(Overlay)下面的交叉分析功能(Intersect)
C. 分别输入上面抽出的水域多边形和该地区其他年份区的分类后的多边形进行intersect分析
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 calculator”对话框下面的空白处
H.在对话框“field calculator”最下面填上“output”
2、然后分别计算水体在该地区其他年份年转入和转出情况
A.点击“potions”下面的“select by attribute”按钮，并选择“GRID_COD_1”
B.令GRID_COD_1,1点击apply，此时会统计出该地区其他年份年仍然为水域的多边形
C.右键点击“newarea”选择统计功能(“statistics”)，就会得到从然为水域的面积 D.分别计算出GRID_COD_1,2(林地)、3(草地)、4(城镇点)、5(耕地)、6(滩涂地)的面积就得到从由水域转入的面积。
以此类推
网上有人还采用的方法:
1、 Erdas:erdas--main--interpreter--gis analysis--matrix输入两个年份的分类图，然后进行重编码即可。
2、 Arcview:可以作个运算，其中一期的图*10或者100或者1000(土地利用类型取1位数时候×10，两位100),然后再加上另一期的土地利用图得出一个图的属性，肯定12，15，或者1221，或者1531之类的数据，应该就可以实现了。
3、 Matlab:听说，还没有亲眼见过，据说可以成功
4、ArcGIS:在一期数据的基础上，绘制出其变化的部分(很复杂，工作了量很大)。
浙公网安备 33021202002483