生成土地利用变化转移矩阵的方法
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
生成土地利用变化转移矩阵的方法
这里是网上搜到的生成土地利用变化转移矩阵的几种方法,以飨来者:
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进行转移矩阵制作