Path='D:\GlobalFVC\paper\密度散点图\rrr.xls'
ast_mod=xlsread(Path);
y=ast_mod(:,1);
x=ast_mod(:,2);
[M N]=size(x);
maxx=1;
minx=0;
maxy=1;
miny=0;
tj=zeros(101,101); %经各种算我的数据分成400比较正好。x,y都将近9万个数。图像得到的。
xfd=(maxx-minx)/100;
yfd=(maxy-miny)/100;
for i=1:M
i_tmp=int32((x(i)-minx)/xfd)+1;
j_tmp=int32((y(i)-miny)/yfd)+1;
tj(i_tmp,j_tmp)=tj(i_tmp,j_tmp)+1;
end
tjmax=max(max(tj)); %我程序需要找出最大值
[i,j]=find(tj==tjmax);
i=mean(i) ;%因为最大值不只一个,取个平均。因为我不需要那个确切的点。
j=mean(j);
colormap(jet);
image(tj);
colorbar
hold on;
Path2='D:\GlobalFVC\paper\密度散点图\rrr2.xls';
ast_mod2=xlsread(Path2);
x2=ast_mod2(:,1).*100;
y2=ast_mod2(:,2).*100;
plot(x2,y2,'m*');
评论11