matlab赌徒破产模型转移矩阵,基于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.

更多需求，请查看个人介绍