基于matlab的土地利用转移矩阵及土地利用转移空间分布

%土地利用转移矩阵
假设目前有两个不同年份的土地利用数据，分别是1990年和2000年的，要求不同土地利用类型间的相互转换关系，本文以matlab为例进行求解，非常的快速方面，代码如下

[a,R]=geotiffread('D:\qixiang\lucc\8km土地利用类型\1990年土地利用类型.tif');%先导入纬度数据
info=geotiffinfo('D:\qixiang\lucc\8km土地利用类型\1990年土地利用类型.tif');
data1990=importdata('D:\qixiang\lucc\8km土地利用类型\1990年土地利用类型.tif');%土地类型分别是123456。
data2000=importdata('D:\qixiang\lucc\8km土地利用类型\2000年土地利用类型.tif');
zy=zeros(size(data1990,1),size(data1990,2))；
counsum_value=[];
for i=1:6
    sy=find(data1990==i);
    countsum=[];
    for j=1:6
        sy1=find(data2000==j);
        sy_inter=intersect(sy,sy1);
        countsy=length(sy_inter);
        countsum=[countsum,countsy];
        zy(sy_inter)=i*10+j;%得到转移前后的量，比如i=1；j=2，则结果是12，表示由1转移到2
    end
    counsum_value=[counsum_value;countsum];
end
xlswrite('D:\qixiang\lucc\8km土地利用类型\990-2000年的土地利用转移矩阵.xlsx',counsum_value)
filename='1990-2000年土地利用空间转移分布图.tif';
geotiffwrite(filename,zy,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);

结果如下：

image.png

表示的含义见下图：

image.png
空间转移图如下所示：
QQ截图20181101105616.jpg

图中第一个数表示1990年的土地利用类型，第二个数是2000年的土地利用类型。如21，则表示1990年土地利用为2,2000年土地利用为1，表示该像元在1990-2000年间土地利用类型由2变成了1.