关于土地分类格式互转、土地利用转移矩阵计算详解

一、ENVI、ArcGIS格式互转

1.1 ENVI分类格式转TIF

ENVI分类格式使用File下面的Save as另存为TIF格式：

对于部分5.3版本可能会出现报错：

那我们可以使用Classic模式：

然后使用ArcGIS打开：但是可能会发现分类图像的值域范围可能和ENVI不一样：
然后我们符号系统里选唯一值：
这样值域范围就对应了：

还可以设置标签名称：

1.2 ArcGIS Coverage格式转TIF

对于Coverage这种格式，ArcMAP打开：

右键导出：

ENVI打开后：右键新建色带

注意需要将默认色带删除，自己重新建立色带并设置为唯一值对应：

导出为ENVI标准分类格式：

1.3 其他分类读取

对应清华大学GCL10或者30，我们加载数据后同样的操作即可：

二、关于拼接、重分类

下面我们以从中国科学院地理科学与资源研究所下载各省土地利用数据为例进行操作演示：

下图是下载的1970年代末期（1980年）中国土地利用现状遥感监测数据：

全部加载后：

因为各省数据的名称都是一样的，直接拼接容易报错或者不出结果，我们先把数据图层进行重命名：

镶嵌至新栅格：

拼接完成后选唯一值：

拼接后的结果：

这是这套数据的空间参数：

这是这套数据的分类编码：

我们下面按照一级编码对这个数据进行重分类，注意新旧值的对应关系：
重分类后的结果：

三、转移矩阵制作

3.1 ArcGIS里操作：

上面我们把1980年的数据进行了重分类，接下来我们按照同样的方法来处理2018年的数据：

两期数据都处理好了，然后进行栅格转矢量，建议不勾选简化面：

两期数据的矢量进行相交分析：

相交分析的结果添加面积字段：

计算几何面积：

然后使用EXCEL打开相交分析数据的.dbf文件，创建数据透视表：

拖动code1到列标签，拖动code2到行标签，area到数值求和：

数据稍加整饰即可得到转移矩阵，但是EXCEL读取数据的行数有限，写入最大行数是 1048576，本次数据已经远超这个数值，因此这样算出来是不准的，但是这个方法适用于数据量比较小的分类结果。

对于面特别多的数据，可以先进行融合操作，根据分类字段将面合并，然后进行相交分析，：

但是此方法步骤繁琐，下面我们介绍在ENVI下的处理方法。

3.2 ENVI里操作：

ENVI里操作需要满足三个条件：1.ENVI标准分类格式 2.数据坐标投影一致 3.数据行列数一致。

首先我们打开数据后查看元数据后的行列数，发现两者的行列数并不一致：

接下来我们使用layerstack工具，将这个两个数据组合在一起，但是因为这个数据的坐标系统ENVI并不识别：

为了方便后续操作我们将这个坐标系统新建入库：

接着使用layerstack工具，导入这两个数据，坐标系选择我们刚才新建的，注意重采样方法，选最邻近，因为这是土地里分类数据，
最邻近可以不改变原有分类数据的像元值，假如选双线性，那么，数据中某两个像元的值一个是3一个是5，可能会插值出一个4，这样每个类别的个数就改变了。

注意layerstack的时候，选第二个可以只保留两个数据的公共相交区域：

通过layerstack工具，两个数据的行列数一致，现在将每个数据导出，可以仍然layerstack，选择spectral subset，选中其中一个数据先导出，同样的方法导出另一个数据：

导出的数据背景值可能不是0，这种情况需要做一下掩膜，首先构建掩膜文件，我们先构建1980的MASK：

我们的数据有效值域是1-7：

然后对1980年的数据（layerstack导出的数据）进行应用掩膜：

右键新建色带后导出Class格式：

可以编辑头文件对导出的Class格式编辑分类名称：

对于某些5.3版本可能会报错，我们可以用classic模式进行编辑：

设置忽略背景值为0：

同样的操作对2018年的数据进行编辑，然后两者数据的都统一了：

接下来我们进行变化检测，计算转移矩阵，由于ENVI计算出的转移矩阵是反向的，因此我们前时相选2018，后时相选1980，这样计算出来的转移矩阵才符合平常的习惯：

导出TXT：

EXCEL导入TXT：

稍加整饰即可：

四、变化图谱

我们以刚才两期数据为例，来计算一下建设用地的变化图谱，比如计算下1980-2018年间，其他地类转为建设用地的分布：
使用地图代数工具，输入下列公式，该公式是将1980所有像元乘以10，这样比如耕地从1变成10，再加上2018，那么耕地保持不变的像元会变成11，耕地转成林地的像元变成12，…

计算完成后在符号系统里只显示与转入建设用地相关的值，并附上属性标签：

可以切换到布局视图插入多个数据框：

调整数据框位置;

复制数据到每个数据框：

接下来插入图例，指北针、比例尺这些要素就可以出图了：

好了今天的介绍就到这里了，更多文字点击博客首页：

https://blog.csdn.net/qq_46071146

关于土地分类格式互转、土地利用转移矩阵、变化图谱计算详解相关推荐

python语言的格式框架_django框架模板语言使用方法详解
本文实例讲述了django框架模板语言使用方法.分享给大家供大家参考,具体如下: 模板功能作用:生成html界面内容,模版致力于界面如何显示,而不是程序逻辑.模板不仅仅是一个html文件,还包括了页 ...
请写出sfr和sbit的语句格式_单片机sfr和sbit关键字，详解
sfr用于将一个单片机的特殊功能寄存器(special function register)赋值给一个变量,这样在后面的程序中就可以中这个变量指引(refer to)该寄存器. sbit与sfr用法类 ...
snort + barnyard2如何正确读取snort.unified2格式的数据集并且入库MySQL（图文详解）...
不多说,直接上干货! 为什么,要写这篇论文? 是因为,目前科研的我,正值研三,致力于网络安全.大数据.机器学习研究领域! 论文方向的需要,同时不局限于真实物理环境机器实验室的攻防环境.也不局限于真实物 ...
arcgis使用arcpy计算土地利用LUCC变化转移矩阵、马尔科夫模型变化表，最好用的方法TabulateArea实现两期数据变化统计
arcgis使用arcpy计算土地利用LUCC变化转移矩阵.马尔科夫模型变化表,最好用的方法TabulateArea实现两期数据变化统计: 下载链接:LUCC变化转移矩阵Arcpy脚本 coding ...
土地利用转移矩阵、土壤侵蚀模数变化等此类用地类型变化的Arcgis操作步骤
土地利用转移矩阵.土壤侵蚀模数变化等此类用地类型变化的Arcgis操作步骤总的步骤:两期栅格转面-融合-相交 1.两期土地利用栅格或矢量图 2.栅格转面有时会出错,是因为栅格数据不是整型或字符串 ...
土地利用转移矩阵图怎么做_土地利用转移矩阵生成的几种方法
土地利用转移矩阵生成的几种方法土地利用类型转换矩阵为相同区域在不同时间段内的土地利用类型的相互转换关系,一般用二维表来表达,从二维表中可以快速查看各个地类间相互转化的具体情况.比如某一类别的土地有 ...
ArcGis空间分析学习：土地利用动态变化分析
说明:本实验主要参考GIS空间分析实验教程-田永中,适合初学者参阅目录一.实验理论 (1)实验目的 (2)实验内容 (3)实验原理 (4)实验方法 (5)实验流程 (6)实验材料二.实验步骤 ( ...
ArcGIS中实现土地利用转移矩阵
土地利用转移矩阵,就是根据同一地区不同时相的土地覆盖现状的变化关系,求得一个二维矩阵.通过对得到的转移矩阵进行分析,能够得到2个时相,不同的地类之间相互转化的情况,它描述了不同的土地利用的类型在不同年 ...
波段合成，去除黑边并制作土地利用转移矩阵（附练习数据下载）
这次呢,我们来做一做第七届全国大学生GIS技能大赛试题A(上午),因为这次的大赛其实有蛮多套试题,所以我们可能会多写一点类似的文章哦~~试题A告诉我们:土地利用每年都在发生变化,在统计一些大区域的土地 ...
ArcGIS中的土地利用变化分析详解
本篇主要是针对矢量数据的分析. 一.不同时期的土地利用矢量数据,如何分析其图形及属性变化? 土地利用图(左图为1993年,右图为2003年) 思路如下: 可以先对2个图层进行Union操作,然后在结果 ...

关于土地分类格式互转、土地利用转移矩阵、变化图谱计算详解