生成土地利用变化转移矩阵的方法

生成土地利用变化转移矩阵的方法
这里是网上搜到的生成土地利用变化转移矩阵的几种方法，以飨来者：
A 栅格数据做转移矩阵
B 矢量数据做转移矩阵
1 来源：
/benben-sky/blog/static/24530388200811100256763/作者DAHONGME
根据你的数据类型选用不同的数据生成方法
若你的数据是Raster格式：则有如下方法
1 Erdas Imagine----Interpreter---Gis Analysis---Matrix，输入两个时相的Raster 数据即可
做这一步之前记得先对两时相的数据进行重编码（nterpreter---Gis Analysis---Recode）
一般运行如果出现错误肯定是重编码没做好，请继续查证。

2 先在Erdas中利用Modeler 计算如下公式
NC(I,J)=NC(I)*10+NC(J),(J>I)
其中：NC(I,J)表示i，j 两年份的土地利用变化图；NC（i）表示i年份遥感分类影像；NC（j）表示j年份的遥感分类影像。

在此计算的基础上，将以上变化影像图转化为BIL格式，再利用ARC/INFO GRID模块将影像转为GRID格式，然后利用GRID模块中的属性表（vat）查看命令对影像灰度值进行统计，最后得出土地利用转化举证。

（注:此方法本人尚未实现过，不知可行否）。

若数据是Vector格式
1 Erdas Imagine----Interpreter---Gis Analysis---Matrix，输入两个时相的Vector 数据即可
此时注意输出栅格大小不应设的太小要不一运行就会提示你的空间不足
做这一步之前，请做好前期的地理编码。

2 ArcView3.3加载spatial analysis模块
把两时相的Vector图转成grid格式（当然中间有一些单位的设置根据你做的图的分辨率来设置即可）analysis---mapcaculate 直接计算即可。

3 把两期解译完的Vector文件在arctoolbox——overlay——union中叠加，注意：两个文件不能用同一个字段名，比如一个用93Type，另一个时相则用
00Type
叠加后的文件在Arcmap中打开，选中文件，然后点右键——Property——空间查询，输入条件语句，比如：93Type=…1‟And 00Type=…2‟；查询结果即为第一种类型转化为第二种类型的图形，可以另建一图层比如：12，把查询结果复制到12图层上。

统计出面积，依进行，就可以得到土地利用类型转移矩阵。

最后输出土地利用变化图，如下图所示：
2 来源：/blog/static/4537466320094271121474/
作者：许科研
转移矩阵在土地利用变更调查时经常用到。

arcgis里没有直接的工具来计算转移矩阵，但我们可以利用已有的工具间接来求。

下面介绍一下利用arcgis求转移矩阵的方法。

现有2005、2008两期土地利用现状数据（polygon.shp）
分别将两期数据用Data Management toolbox > Generalization > Dissolve工具，按各自的地类编码合并要素，注意选中最下面的“create multipat features”选项处理
好数据后，用intersect工具求两期数据的交集交集的结果：将属性表增加面积字段，并计算面积结果：
从上图中已经可以看出土地转移的情况：dlbm2005是2005年时该块地的地类代码，dlbm2008是2008年时该地块的地类代码，最后一列是相应地块的面积。

接下来将属性表导出到dbf文件，然后在excel中打开，执行数据菜单中的“数据透视表和数据透视图”命令点完成，然后将""dlbm2005作为行字段，dlbm2008作为列字段""，area作为数据拖到相应的区域生成的数据透视表将数据复制到新的工作表中，适当的修整后即生成了土地利用转移矩阵。

3 来源：/124793909
作者：西海
写文章的过程中，同样也用到了土地利用/景观转移矩阵。

转移矩阵的制作很多软件可以实现，有Arcview、Arcmap、还有人用matlab。

但是网上都写的都不详细，在这里，本文才用arcmap的叠加操作＋excel的数据统计功能，可以轻松完成这一过程。

相信看完了本文，初学者或者对GIS了解不多的人，也能轻松的制作出转移矩阵。

方法一：
1数据准备
准备好两期的coverage或shape文件。

（注意：拓扑关系要建好，不然无法计算），采用Arcmap打开上述两期文件。

2叠加操作
选择Arcmap里面的“ArcToolbox”按钮下面的“Analysis Tools”工具下面的叠加分析模块（Overlay）下面的交叉分析功能（Intersect）。

选择已经打开的两期数据进行叠加分析。

叠加的结果存在一个可以找到的地方。

同时，把叠加后的结果添加在Arcmap里。

3输出叠加文件的属性数据
A、右键打开intersect产生的矢量文件的数据属性表（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：在一期数据的基础上，绘制出其变化的部分（很复杂，工作了量很大）。