土地利用景观转移矩阵方法汇总

最简单的土地利用/景观转移矩阵的制作（转）(2009-05-19 03:36:24)

写文章的过程中，同样也用到了土地利用/景观转移矩阵。转移矩阵的制作很多软件可以实现，有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”后选择“hel p”将

Dim Output as double

Dim pArea as Iarea

Set pArea = [shape]

Output = pArea.area

拷入到“field calculator”对话框下面的空白处

E、在对话框“field calculator”最下面的空白处填上“output”

F、在属性表点击“opti ons”按钮，→export（导出），属性以.dbf格式结果存储。

4 Excel进行转移矩阵制作

A、刚才存储的.dbf文件可以使用excel打开，打开的结果重新保存为.xls（excel）文件。

（dbf文件不能保存）

B、重新打开刚存储的excel文件。

C、选中所有数据，选择数据\数据透视表和数据透视图。将两期数据的id值分别拖入行列字段，然后以‘newarea’字段作为数据项拖入计算区域。得到的便是转移矩阵。如果发现矩阵的形式不美观，可以对单位进行修改。

方法二：作某一地物与其他不同年份地类

例如：

1、做水体和其它年份分类地物的交集

A. 分别打开刚刚产生的水域多边形和该地区其他年份的分类矢量文件

B. 选择gis里面的“AacToolbox”按钮下面的“Analysi s Tool s”工具下面的叠加分析模块（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.点击“poti ons”下面的“select by attribute”按钮，并选择“GRID_COD_1”

B.令GRID_COD_1＝1点击apply，此时会统计出该地区其他年份年仍然为水域的多边形

C.右键点击“newarea”选择统计功能（“stati stics”），就会得到从然为水域的面积

D.分别计算出GRID_COD_1＝2（林地）、3（草地）、4（城镇点）、5（耕地）、6（滩涂地）的面积就得到从由水域转入的面积。

以此类推

网上有人还采用的方法：

1、 Erdas：erdas--main--interpreter--gi s analysi s--matrix输入两个年份的分类图，然后进行重编码即可。

2、 Arcview：可以作个运算，其中一期的图*10或者100或者1000（土地利用类型取1位数时候×10，两位100),然后再加上另一期的土地利用图得出一个图的属性，肯定12，15，或者1221，或者1531之类的数据，应该就可以实现了。

3、 Matlab：听说，还没有亲眼见过，据说可以成功

4、ArcGIS：在一期数据的基础上，绘制出其变化的部分（很复杂，工作了量很大）。

根据你的数据类型选用不同的数据生成方法

若你的数据是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图层上。统计出面积，依进行，就可以得到土地利用类型转移矩阵。

刚才看了缘老大发的帖子，题目是：利用arcigis生成土地利用转移矩阵。方法流程讲的很详细。最近我也用了好多方法来做，其中个人认为最简便的方法是利用arcgis的spatial analyst tools---zonal---tabulate area

详细步骤请见------方法很简单。另提醒（矢量，栅格都可以哦）1.找到zonal--tebulate area→

2.导入两期的矢量或栅格图→,你就点OK,让arcgis运算即可。

3.结果显示，在arcgis里导出*.dbf→（这就是结果了）。导出来你用excel 打开，想插到那里就插到那里。

利用ARCGIS生成土地利用转移矩阵

上一篇/ 下一篇 2009-02-01 19:42:58

查看( 825 ) / 评论( 62 )

转移矩阵在土地利用变更调查时经常用到。arcgis里没有直接的工具来计算转移矩阵，但我们可以利用已有的工具间接来求。下面介绍一下利用arcgi s求转移矩阵的方法。

现有2005、2008两期土地利用现状数据：

1.jpg

2.jpg

3.jpg

4.jpg

在求转移矩阵前先做好下列工作：

1、备份数据。

2、分别将两期数据用di s sol v e工具按各自的地类编码合并要素，注意选中最下面的“create multipat features”选项

8.jpg

处理好数据后，用intersect工具求两期数据的交集

0.jpg

00.jpg

交集的结果：

5.jpg

6.jpg

将属性表增加面积字段，并计算面积

9.jpg

10.jpg

11.jpg 结果:

12.jpg

从上图中已经可以看出土地转移的情况：dlbm2005是2005年时该块地的地类代码，dlbm2008是2008年时该地块的地类代码，最后一列是相应地块的面积。

接下来将属性表导出到dbf文件，然后在excel中打开，执行数据菜单中的“数据透视表和数据透视图”命令

13.jpg

14.jpg

15.jpg

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

生成土地利用变化转移矩阵的方法 这里是网上搜到的生成土地利用变化转移矩阵的几种方法，以飨来者： A 栅格数据做转移矩阵 B 矢量数据做转移矩阵 1 来源： https://www.docsj.com/doc/aa1395114.html,/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图层上。统计出面积，依进行，就可以得到土地利用类型转移矩阵。 最后输出土地利用变化图，如下图所示：

土地利用转移矩阵

土地利用转移矩阵 现有ABCD...期等土地利用数据： 1）矢量操作：使用union或intersect，将A期和B期的矢量图层进行叠加，然后数据库操作-> 字段计算器-> “新字段= A期.土地利用类型字段*100＋B期.土地利用类型字段”-> 然后根据“新字段提供的信息”进行面积统计等后期操作。 注：*100 是为了区分AB两期的土地利用类型，一位的土地利用类型分类可以*10、100...， 两位的分类可以*100、*1000...，但不要超过10000，因为太大了就变成连续型栅格了； union 相当于集合操作的并集，两期图层取并集； intersect 相当于集合操作的交集，两期图层取交集； 当各期的图层的范围相同时两者都一样，图层范围不相同时，请仔细选择使用； 面积统计：在字段计算器中完成。 2）栅格操作：在空间分析模块中-> 栅格计算器-> 输入：[A期] * 100 + [B期] -> 生成新的栅格图层-> 然后再根据“新生成的栅格图层”进行面积统计等后期操作。 注：预处理：若是“连续栅格”，通过->空间分析模块的重分类-> 转成“唯一值的栅格”-> 再进行计算； *100 （同上）； [A期]、[B期]...是栅格图层名； 某种土地利用类型的面积＝某种土地利用类型的栅格数量* (栅格分辨率* 栅格分辨率)。 3）A -> B的土地利用变化： 101 A期类型1 -> B期类型1（没变化） 102 A期类型1 -> B期类型2（变化） 103 A期类型1 -> B期类型3（变化） ... 201 A期类型2 -> B期类型1（变化） 202 A期类型2 -> B期类型2（没变化） 203 A期类型2 -> B期类型3（变化） ... n0n A期类型n -> B期类型n（？？？） 或 n00n A期类型n -> B期类型n

ArcGIS土地利用转移矩阵

一、数据准备（图1） 准备两幅不同时相的土地利用现状图（shp格式），每幅图的属性表都要有一个表示土地利用类型的字段，并且要使用不同的名称加以区分，如Type1995，Type2000。土地利用类型名称必须统一，并且完整，如都使用“城镇用地”、“有林地”等。 二、数据融合（图2） 在ArcMap里分别打开两个时相的图层，打开ArcToolbox，选择Data Management Tools | Generalization | Dissolve工具。Input Feature选择要融合的图层，Output Feature Class选择输出结果存储的位置及名称，Dissolve Field(s)选择土地利用类型字段（如Type1995），然后勾选Creat multipart features选项，点击OK完成。重复此过程，对另一时相数据进行融合。此步骤使相同利用类型的记录融合为一个记录，以提高后面步骤的计算速度。 三、叠置分析（图3） 在ArcMap中打开两个时相融合后的数据，在ArcToolbox中选择Analysis Tools | Overlay | Intersect工具，Input Features选择两个时相的图层，Output Feature Class 选择叠加结果存储的位置及名称，其余选项可以忽略，单击【OK】完成。 四、计算面积并导出属性表（图4-6） 在ArcMap中打开叠加后的图层数据，在该图层上右键打开属性表，选择Option | Add field… 新建一个字段，命名为NewArea。 在Editer工具条中选择Editer | Start Editing，然后在属性表中NewArea字段上单击右键选择Calculate Geometry… ，在打开的Calculate Geometry对话框中，Property选择Area，Units选择要使用的面积单位，单击【OK】完成图斑面积计算。依次选择Editer | Save Edits / End Editing保存和退出编辑状态。 在属性表中选择Option | Export… 将属性表保存为dbf文件。 五、制作转移矩阵（图7-10）（以Excel2007为例） 在Excel中打开上一步保存的dbf，另存为Excel格式并打开。在Excel中选中所有数据（不要点左上角，只选择有效数据），点击【插入】选项卡，选择【数据透视表】|【数据透视表】，点击【确定】。 在打开的数据透视表中按图示将字段拖入相应区域。 Excel自动计算矩阵，将该表稍事整饰就得到美观的土地利用转移矩阵。矩阵中 r(I, j)就表示i类型向j类型转移的土地面积，空值表示i类型向j类型没有转移。 阅读全文(131) / 评论 / 扔小纸条 / 给曾经善良过留言 收藏: QQ书签 https://www.docsj.com/doc/aa1395114.html, / 订阅: Google 抓虾

ArcGIS土地利用转移矩阵

一、数据准备(图1) 准备两幅不同时相得土地利用现状图(shp格式),每幅图得属性表都要有一个表示土地利用类型得字段,并且要使用不同得名称加以区分,如Type1995,Type2000。土地利用类型名称必须统一,并且完整,如都使用“城镇用地”、“有林地”等。 二、数据融合(图2) 在ArcMap里分别打开两个时相得图层,打开ArcToolbox,选择Data Management Tools | Generalization | Dissolve工具。Input Feature选择要融合得图层,Output Feature Class选择输出结果存储得位置及名称,Dissolve Field(s)选择土地利用类型字段(如Type1995),然后勾选Creat multipart features选项,点击OK完成。重复此过程,对另一时相数据进行融合。此步骤使相同利用类型得记录融合为一个记录,以提高后面步骤得计算速度。

三、叠置分析(图3) 在ArcMap中打开两个时相融合后得数据,在ArcToolbox中选择Analysis Tools | Overlay | Intersect工具,Input Features选择两个时相得图层,Output Feature Class 选择叠加结果存储得位置及名称,其余选项可以忽略,单击【OK】完成。 四、计算面积并导出属性表(图4-6) 在ArcMap中打开叠加后得图层数据,在该图层上右键打开属性表,选择Option | Add

field… 新建一个字段,命名为NewArea。 在Editer工具条中选择Editer | Start Editing,然后在属性表中NewArea字段上单击右键选择Calculate Geometry… ,在打开得Calculate Geometry对话框中,Property 选择Area,Units选择要使用得面积单位,单击【OK】完成图斑面积计算。依次选择Editer | Save Edits / End Editing保存与退出编辑状态。

土地利用转移矩阵生成的几种方法

土地利用转移矩阵生成的几种方法 土地利用转移矩阵生成的几种方法 根据你的数据类型选用不同的数据生成方法 若你的数据是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图层上。统计出面积，依进行，就可以得到土地利用类型转移矩阵。

利用栅格计算器进行土地利用类型转移矩阵计算

利用栅格计算器进行土地利用类型转移矩阵计算 1.数据准备 准备好两期的土地利用分类数据裁剪_86和裁剪_95，用arcmap 打开这两期影像文件：

2.转移矩阵计算 打开Arctoolbox window 窗口，在Arctoolbox window窗口中选择Spatial Analyst Tools>>Map Algebra>>Raster Calculator. 双击Raster Calculator，出现如下显示窗口。

表达式输入框 在表达式输入框中输入表达式:"裁剪_86" * 10 + "裁剪_95"，在Raster Calculator中双击map algebra expression中的文件即为 选中此文件，Output raster中设置输出文件目录及文件名fangfa_1， 选择OK。 如图得到栅格计算器结果文件fangfa_1，右击layers中文件fangfa_1, 选择open attribute table,查看fangfa_1的属性表。

在fangfa_1属性表中选择table options>>export,选择路径和文件名，输出一个.dbf文件表格。 3.结果分析 按照刚刚的输出路径，找到该表，用excel打开该表格。每个像元大小为30*30，所以需要将表中count字段乘以900即为面积（单位为米）。

需要注意的是，“23”表示由裁剪_86年的第2类型转化为裁剪 _95年的第3类型，“56”表示由裁剪_86年的第5类型转化为裁剪_95 年的第6类型，等等。 调整excel表中的数值为下表所示： 调整后的转移矩阵结果表中右边即为86年到95年土地利用类型转移矩阵结果。

ArcGIS生成景观转移矩阵

ARCGIS转移矩阵 2012-03-08 20:28:13| 分类：ARCGIS | 标签：|字号大中小订阅 土地利用转移矩阵生成的几种方法 查阅相关的资料，也没有得到土地利用类型转换矩阵确切的定义，我理解为不同时间段内同一区域内土地利用类型的相互转换关系，一般用二维表来表达，从二维表中可以快速查看各个地类间相互转化的具体情况。比如某一类别的土地有百分之多少（或者面积）分别转化成了其他的土地类型，现在某类型的土地分别是由过去的哪些类别转化而来的等等。还可以生成变化统计栅格图（掩膜图像），它描述了前后两幅土地分类图之间的地类发生转变的位置和类别。 土地利用类型转换矩阵可以从两幅栅格图中计算得到，也可以从两个矢量文件中计算获得。下面介绍在ENVI下从两幅分类结果的栅格图中计算土地利用类型转换矩阵。 1、准备数据 两个时相的土地利用分类结果，它是单波段、专题类型的伪彩色图像（ENVI Classification）。 2、计算转换矩阵 打开两个土地利用分类结果。 （1）在主菜单中，选择Basic Tools → Change Detection → Change Detection Statistics。 （2）分别在Initial State对话框和final state对话框中选择前一时相和后一时相的土地利用结果。 （3）在Define Equivalent Classes对话框中（图1），如果两个土地利用分类名称一致，系统自动将Initial State Class和Final State Class对应，否则手动选择，单击Add Pair按钮选择。 （4）选择对应的地物类型之后，单击OK按钮，出现图2对话框。选择生成图表表示单位（Report Type）：像素（Pixels）、百分比（Percent）和面积（Area）。选择Output Classification Mask Images?为YES，输出掩膜图像，选择输入路径及文件名。 （5）单击OK，执行土地利用类型转换矩阵计算过程。

利用栅格计算器进行土地利用类型转移矩阵计算

利用栅格计算器进行土地利用类型转移矩阵计算 1. 数据准备 准备好两期的土地利用分类数据裁剪 _86和裁剪_95,用arcmap 耳护 吕dil 无― ■誉U jQ.啊ixifsj C"f 04H 4 Er|p □ * ?吕 * r* *- ■已* ■ ■ 衣、e ii a'laj - Sk Uja -吊“/■直* 136糾■咖 T… j* m u 2： +■ h < . jj-；-■: 一飞既耳 ? w ?也* ? f A ?號.￡ -丄-$ ? 1 D 应 1 “I ￡1 Si 世 她1 ?? M E 叮 車. T,如 4M D WK . 打开这两期影像文件: hte tidta ZiriCLWIawHVjtf]? Lii^ Bt^hFMilbA 已

2. 转移矩阵计算 打开Arctoolbox window 审窗口，在Arctoolbox window 窗口中选择Spatial Analyst Tools>>Map Algebra>>Raster Calculator. ArcTczi'fcicw CcnvwGtan Tciali -(XsM l Urr : 口仃晶商孑$ 篡? Ora T ook S/ ? (dhiing -?灯z匚口dMpi T?li ?-Ana ^8 Tools T Li/jar商粘曲疋：怕T OO I E MutidiELhsiori T DD I K --电I Netwc rk I oo 去 ??旳p ml衍b*址5取起 O 务heEHit* Tgk Q号时细r Tfl-eh -Ar ■tly- ■ Tcjoh 7 輪Cixdtianal ■-电El”弓阿 -bizlance “ 4(^ E*Hu?pr> 岭< I3*suritf^wcr *令11河「口1?卵 r i|^ lrrr^rpdaiicn ~ 輪Loc^l 岭M-zp A get r s --\ Cjikulatx?r Malh MmHiwiine N?ghte

ENVI下土地利用类型转换(变化)矩阵(转)

ENVI下土地利用类型转换（变化）矩阵（转） 土地利用类型转换矩阵可以从两幅栅格图中计算得到，也可以从两个矢量文件中计算获得。下面介绍在ENVI下从两幅分类结果的栅格图中计算土地利用类型转换矩阵。 1、准备数据 两个时相的土地利用分类结果，它是单波段、专题类型的伪彩色图像（ENVI Classification）。 2、计算转换矩阵 打开两个土地利用分类结果。 （1）在主菜单中，选择Basic Tools →Change Detection →Change Detection Statistics。（2）分别在Initial State对话框和final state对话框中选择前一时相和后一时相的土地利用结果。 （3）在Define Equivalent Classes对话框中（图1），如果两个土地利用分类名称一致，系统自动将Initial State Class和Final State Class对应，否则手动选择，单击Add Pair按钮选择。 （4）选择对应的地物类型之后，单击OK按钮，出现图2对话框。选择生成图表表示单位（Report Type）：像素（Pixels）、百分比（Percent）和面积（Area）。选择Output Classification Mask Images?为YES，输出掩膜图像，选择输入路径及文件名。 （5）单击OK，执行土地利用类型转换矩阵计算过程。 图1 Define Equivalent Classes对话框

图2 选择数据参数 3、查看结果 （1）如图3为得到的土地利用类型转换矩阵结果。横字段表示前一时间段（Initial State）的土地利用类别，纵字段为后一时间段（Final State）的土地利用类别。横字段和纵字段交叉处表示变化值，如有2520900平方米林地用地变化为草地。 图3 土地利用类型转换矩阵 （2）还可以为每一个地类生成一个变换掩膜图像，图4所示为其中一个地类的掩膜图像。掩膜图像的灰度值表示变化类型，如这里的2{草地}表示林地变化为草地的像元。

转移矩阵计算

写文章的过程中，同样也用到了土地利用/景观转移矩阵。转移矩阵的制作很多软件可以实现，有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”

土地利用矢量栅格转移矩阵计算

最简单的土地利用/景观转移矩阵的制作(转)(2009-05-19 03:36:24) 标签:杂谈分类:相关资料 写文章的过程中,同样也用到了土地利用/景观转移矩阵。转移矩阵的制作很多软件可以实现,有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”

土地利用景观转移矩阵方法汇总

最简单的土地利用/景观转移矩阵的制作（转）(2009-05-19 03:36:24) 写文章的过程中，同样也用到了土地利用/景观转移矩阵。转移矩阵的制作很多软件可以实现，有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”后选择“hel p”将 Dim Output as double Dim pArea as Iarea Set pArea = [shape] Output = pArea.area 拷入到“field calculator”对话框下面的空白处 E、在对话框“field calculator”最下面的空白处填上“output”

土地利用转移矩阵生成的几种方法

土地利用转移矩阵生成的几种方法 土地利用类型转换矩阵为相同区域在不同时间段内的土地利用类型的相互转换关系，一般用二维表来表达，从二维表中可以快速查看各个地类间相互转化的具体情况。比如某一类别的土地有百分之多少（或者面积）分别转化成了其他的土地类型，现在某类型的土地分别是由过去的哪些类别转化而来的等等。还可以生成变化统计栅格图（掩膜图像），它描述了前后两幅土地分类图之间的地类发生转变的位置和类别。 土地利用类型转换矩阵可以从两幅栅格图中计算得到，也可以从两个矢量文件中计算获得。下面介绍土地利用转移矩阵生成的几种方法： ENVI下从两幅分类结果的栅格图中计算土地利用类型转换矩阵 1、准备数据 两个时相的土地利用分类结果，它是单波段、专题类型的伪彩色图像（ENVI Classification）。 2、计算转换矩阵 打开两个土地利用分类结果。 （1）在主菜单中，选择Basic Tools → Change Detection → Change Detection Statistics。 （2）分别在Initial State对话框和final state对话框中选择前一时相和后一时相的土地利用结果。 （3）在Define Equivalent Classes对话框中（图1），如果两个土地利用分类名称一致，系统自动将Initial State Class和Final State Class对应，否则手动选择，单击Add Pair按钮选择。 图1 Define Equivalent Classes对话框 （4）选择对应的地物类型之后，单击OK按钮，出现图2对话框。选择生成图表表示单位（Report Type）：像素（Pixels）、百分比（Percent）和面积（Area）。选择Output Classification Mask Images?为YES，输出掩膜图像，选择输入路径及文件名。

ArcGIS土地利用转移矩阵

一、数据准备(图1) 准备两幅不同时相的土地利用现状图(shp格式),每幅图的属性表都要有一个表示土地利用类型的字段,并且要使用不同的名称加以区分,如Type1995,Type2000。土地利用类型名称必须统一,并且完整,如都使用“城镇用地”、“有林地”等。 二、数据融合(图2) 在ArcMap里分别打开两个时相的图层,打开ArcToolbox,选择Data Management Tools | Generalization | Dissolve工具。Input Feature选择要融合的图层,Output Feature Class选择输出结果存储的位置及名称,Dissolve Field(s)选择土地利用类型字段(如Type1995),然后勾选Creat multipart features选项,点击OK完成。重复此过程,对另一时相数据进行融合。此步骤使相同利用类型的记录融合为一个记录,以提高后面步骤的计算速度。

三、叠置分析(图3) 在ArcMap中打开两个时相融合后的数据,在ArcToolbox中选择Analysis Tools | Overlay | Intersect工具,Input Features选择两个时相的图层,Output Feature Class 选择叠加结果存储的位置及名称,其余选项可以忽略,单击【OK】完成。 四、计算面积并导出属性表(图4-6) 在ArcMap中打开叠加后的图层数据,在该图层上右键打开属性表,选择Option |

Add field… 新建一个字段,命名为NewArea。 在Editer工具条中选择Editer | Start Editing,然后在属性表中NewArea字段上单击右键选择Calculate Geometry… ,在打开的Calculate Geometry对话框 中,Property选择Area,Units选择要使用的面积单位,单击【OK】完成图斑面积计算。依次选择Editer | Save Edits / End Editing保存与退出编辑状态。

土地利用矩阵的转移与分析

宿迁学院 2011 学年— 2012 学年第 1 学期 土地利用矩阵的转移与分析实验作业报告书 专业：资源环境班级： 姓名：学号： 实验地点：实验室D203 任课教师： 实验题目：土地利用矩阵的转移与分析 实验环境： ArcMap 实验目的、实验内容、实验方法和步骤、实验结果及结论分析等： 一、实验目的： 通过对宿迁1980到2000土地利用转移矩阵计算，分析宿迁这20年来土地利用的变化 二、实验内容： 利用Arcmap软件，计算宿迁1980年到2000年，土地利用转移的面积，以此分析宿迁土地利用。

三、实验方法和步骤： 1、双击打开ArcMap，选择“一个新的空地图”，单击“确定”。 2、加载实验文件所在文件夹下“lu19802000_Unew.shp”文件。

3、选中文件，如图所示： 4、右击选择“打开属性表”选项，打开如图所示表：

5、在右下脚“选项”中选择“通过属性选择”跳出如图所示表： 6、在左边“字段”与右边“唯一值”以及中间运算式中依次选中如下字段，形成代码："useid3" = 11 AND "id_new2000" =21

7、选择“应用”后，单击“选中的”，得出如下图所示表： 8、右击选中“newarea”，在下拉菜单中选择“统计”，得出如下图所示： 9、在表格中记录“总计”一栏数据（四舍五入）。

四、实验结果及结论分析 从计算的LUCC表格可以看出，对角线的面积最大，可以看出宿迁从1980年到2000年土地类型变化不大。土地类型为耕地，林地转化城镇用地的面积较大，说明宿迁这20年中，随着城市的不断发展，宿迁城镇不断扩大，耕地林地渐渐减少。 通过本次的实验，我初步掌握了Arcmap软件，并学会了运用Arcmap软件计算宿迁1980年到2000年土地利用转移的面积，初步分析了宿迁市的土地利用的转移用途等。我学会了遇到难题需要及时发现及时解决，要懂得谦逊请教其他人，不能盲目的做题，需要从多角度来思考问题。

ArcGIS土地利用转移矩阵

土地利用转移矩阵生成的几种方法 查阅相关的资料，也没有得到土地利用类型转换矩阵确切的定义，我理解为不同时间段内同一区域内土地利用类型的相互转换关系，一般用二维表来表达，从二维表中可以快速查看各个地类间相互转化的具体情况。比如某一类别的土地有百分之多少（或者面积）分别转化成了其他的土地类型，现在某类型的土地分别是由过去的哪些类别转化而来的等等。还可以生成变化统计栅格图（掩膜图像），它描述了前后两幅土地分类图之间的地类发生转变的位置和类别。 土地利用类型转换矩阵可以从两幅栅格图中计算得到，也可以从两个矢量文件中计算获得。下面介绍在ENVI下从两幅分类结果的栅格图中计算土地利用类型转换矩阵。 1、准备数据 两个时相的土地利用分类结果，它是单波段、专题类型的伪彩色图像（ENVI Classification）。 2、计算转换矩阵 打开两个土地利用分类结果。 （1）在主菜单中，选择Basic Tools →Change Detection →Change Detection Statistics。 （2）分别在Initial State对话框和final state对话框中选择前一时相和后一时相的土地利用结果。 （3）在Define Equivalent Classes对话框中（图1），如果两个土地利用分类名称一致，系统自动将Initial State Class和Final State Class对应，否则手动选择，单击Add Pair按钮选择。 （4）选择对应的地物类型之后，单击OK按钮，出现图2对话框。选择生成图表表示单位（Report Type）：像素（Pixels）、百分比（Percent）和面积（Area）。选择Output Classification Mask Images?为YES，输出掩膜图像，选择输入路径及文件名。 （5）单击OK，执行土地利用类型转换矩阵计算过程。 图1 Define Equivalent Classes对话框 图2 选择数据参数 3、查看结果

遥感专题讲座——影像信息提取(六、土地利用类型转换混淆矩阵)

土地利用类型转换（变化）混淆矩阵 查阅相关的资料，也没有得到土地利用类型转换矩阵确切的定义，我理解为不同时间段内同一区域内土地利用类型的相互转换关系，一般用二维表来表达，从二维表中可以快速查看各个地类间相互转化的具体情况。比如某一类别的土地有百分之多少（或者面积）分别转化成了其他的土地类型，现在某类型的土地分别是由过去的哪些类别转化而来的等等。还可以生成变化统计栅格图（掩膜图像），它描述了前后两幅土地分类图之间的地类发生转变的位置和类别。 土地利用类型转换矩阵可以从两幅栅格图中计算得到，也可以从两个矢量文件中计算获得。下面介绍在ENVI 下从两幅分类结果的栅格图中计算土地利用类型转换矩阵。 1、准备数据 两个时相的土地利用分类结果，它是单波段、专题类型的伪彩色图像（ENVI Classification）。 2、计算转换矩阵 打开两个土地利用分类结果。 （1）在主菜单中，选择Basic Tools → Change Detection → Change Detection Statistics。 （2）分别在Initial State对话框和final state对话框中选择前一时相和后一时相的土地利用结果。 （3）在Define Equivalent Classes对话框中（图1），如果两个土地利用分类名称一致，系统自动将Initial State Class和Final State Class对应，否则手动选择，单击Add Pair按钮选择。 （4）选择对应的地物类型之后，单击OK按钮，出现图2对话框。选择生成图表表示单位（Report Type）：像素（Pixels）、百分比（Percent）和面积（Area）。选择Output Classification Mask Images?为YES，输出掩膜图像，选择输入路径及文件名。 （5）单击OK，执行土地利用类型转换矩阵计算过程。

土地利用转移矩阵生成的几种方法

土地利用转移矩阵生成的几种方法 遥感2010-11-25 17:35:36 阅读308 评论0 字号：大中小订阅 本文引用自liangshzh0816《ARCGIS转移矩阵》 查阅相关的资料，也没有得到土地利用类型转换矩阵确切的定义，我理解为不同时间段内同一区域内土地利用类型的相互转换关系，一般用二维表来表达，从二维表中可以快速查看各个地类间相互转化的具体情况。比如某一类别的土地有百分之多少（或者面积）分别转化成了其他的土地类型，现在某类型的土地分别是由过去的哪些类别转化而来的等等。还可以生成变化统计栅格图（掩膜图像），它描述了前后两幅土地分类图之间的地类发生转变的位置和类别。 土地利用类型转换矩阵可以从两幅栅格图中计算得到，也可以从两个矢量文件中计算获得。下面介绍在ENVI下从两幅分类结果的栅格图中计算土地利用类型转换矩阵。 1、准备数据 两个时相的土地利用分类结果，它是单波段、专题类型的伪彩色图像（ENVI Classification）。 2、计算转换矩阵 打开两个土地利用分类结果。 （1）在主菜单中，选择Basic Tools → Change Detection → Change Detection Statistics。 （2）分别在Initial State对话框和final state对话框中选择前一时相和后一时相的土地利用结果。 （3）在Define Equivalent Classes对话框中（图1），如果两个土地利用分类名称一致，系统自动将Initial State Class和Final State Class对应，否则手动选择，单击Add Pair按钮选择。 （4）选择对应的地物类型之后，单击OK按钮，出现图2对话框。选择生成图表表示单位（Report Type）：像素（Pixels）、百分比（Percent）和面积（Area）。选择Output Classification Mask Images?为YES，输出掩膜图像，选择输入路径及文件名。 （5）单击OK，执行土地利用类型转换矩阵计算过程。