在计算两幅图像像元水平上的相关性时,由于图像的分辨率高,导致直接画散点图根本就体现不出来大致的趋势,因此需要用密度散点图来进行刻画。散点密度图的核心是在x和y的范围内找到在这个范围内的点的数量,需要合理选择x和y的区间。
首先本文这里以两幅ET的数据作为样例,分辨率为1000m,值的区间再8到1100,两者的最大差距不超过800,选择以5作为两幅图像的间隔,提取如下代码
% author yinlichang3064@163.com
filemean='I:\蒸散发\各ET的均值和趋势\yr_SSEB_2003-2013年ET的均值.tif';
sseb=importdata(filemean);
filemean='I:\蒸散发\各ET的均值和趋势\yr_SSEBopYRB_2003-2013年ET的均值.tif';
ssebop=importdata(filemean);
sy1=find(sseb>0 & ssebop>0); %找到两者的有效值
sseb_1=double(sseb(sy1));
ssebop_1=double(ssebop(sy1));
valuesum=zeros(47961,3); %以47961是因为8:5:1100的长度为219,219的平方为47961
k=1;
for i=8:5:1100 %sseb
sy=find(sseb_1>=i & sseb_1
j1=i+800;
j2=i-800;
if j1>1100
j1=1100;