关于土地分类格式互转、土地利用转移矩阵计算详解
一、ENVI、ArcGIS格式互转
1.1 ENVI分类格式转TIF
1.2 ArcGIS Coverage格式转TIF
1.3 其他分类读取
二、关于拼接、重分类
下面我们以从中国科学院地理科学与资源研究所下载各省土地利用数据为例进行操作演示:
三、转移矩阵制作
3.1 ArcGIS里操作:

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

拖动code1到列标签,拖动code2到行标签,area到数值求和:

数据稍加整饰即可得到转移矩阵,但是EXCEL读取数据的行数有限,写入最大行数是 1048576,本次数据已经远超这个数值,因此这样算出来是不准的,但是这个方法适用于数据量比较小的分类结果。
3.2 ENVI里操作:
ENVI里操作需要满足三个条件:1.ENVI标准分类格式 2.数据坐标投影一致 3.数据行列数一致。


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

然后对1980年的数据(layerstack导出的数据)进行应用掩膜:

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


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

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



设置忽略背景值为0:

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

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




导出TXT:

EXCEL导入TXT:


稍加整饰即可:

四、变化图谱
发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/229875.html原文链接:https://javaforall.net
