转移矩阵计算
写文章的过程中,同样也用到了土地利用/景观转移矩阵。转移矩阵的制作很多软件可以实现,有Arcview、Arcmap、还有人用matlab。但是网上都写的都不详细,在这里,本文才用arcmap的叠加操作,excel的数据统计功能,可以轻松完成这一过程。相信看完了本文,初学者或者对GIS了解不多的人,也能轻松的制作出转移矩阵。
方法一:
1 数据准备
准备好两期的coverage或shape文件。(注意:拓扑关系要建好,不然无法计算),采用Arcmap打开上述两期文件。
2 叠加操作
选择Arcmap里面的“ArcToolbox”按钮下面的“Analysis Tools”工具下面的叠加
分析
定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析
模块(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 O
utput 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:在一期数据的基础上,绘制出其变化的部分(很复杂,工作了量很大)。
看了这些,我想本文所采用的方法,应该是最为简单的一种。如果有什么问题,可以给我留言一起探讨。
同时借此文感谢小曾曾(82444233)、嘉许(330323039)
转移矩阵在土地利用变更调查时经常用到。arcgis里没有直接的工具来计算转移矩阵,但我们可以利用已有的工具间接来求。下面介绍一下利用arcgis求转移矩阵的方法。现有2005、2008两期土地利用现状数据(polygon shp): 在求转移矩阵前先做好下列工作: 1、备份数据。 2、分别将两期数据用dissolve工具按各自的地类编码合并要素,注意选中最下面的“create multipat features”选项 处理好数据后,用intersect工具求两期数据的交集 交集的结果: 将属性表增加面积字段,并计算面积 结果: 从上图中已经可以看出土地转移的情况:dlbm2005是2005年时该块地的地类代码,dlbm2008是2008年时该地块的地类代码,最后一列是相应地块的面积。接下来将属性表导出到dbf文件,然后在excel中打开,执行数据菜单中的“数据透视表和数据透视图”命令 点完成,然后将""dlbm2005作为行字段,dlbm2008作为列字段"",area作为数据拖到相应的区域生成的数据透视表将数据复制到新的工作表中,适当的修整后即生成了土地利用转移矩阵。
在arctoolbox
表格
关于规范使用各类表格的通知入职表格免费下载关于主播时间做一个表格详细英语字母大小写表格下载简历表格模板下载
就是转移矩阵的结果。