Matlab 画散点热力图 城市表层土壤重金属污染分析

Matlab 画散点热力图 城市表层土壤重金属污染分析

利用建模题目 A题 城市表层土壤重金属污染分析 的数据大致来画金属污染的热力图

数据从第一列开始分别是：

数据已经按功能区分类排序好

https://pan.baidu.com/s/1tBmrLiqRAF6tTLRenUxgUQ 数据链接
数据默认放到matlab的bin文件夹下

ww=load('金属的位置以及区域.csv');
%开始画热力图 以Ni离子为例
a=ww;
x=a(:,11);y=a(:,12);nongdu=a(:,7);   %第十一十二列为xy值，第七列为Ni浓度值
xij=min(x(:,:)):(max(x(:,:))-min(x(:,:)))/70:max(x(:,:));
yij=min(y(:,:)):(max(y(:,:))-min(y(:,:)))/70:max(y(:,:));
[xi,yj]=meshgrid(xij,yij);
z=griddata(x,y,nongdu,xi,yj,'nearest');  %这里有三种方法 nearest  linear cubic  可以自行选择
surf(xi,yj,z);
colorbar;
hold on;
%以下是散点图
x1=a(1:44,11);y1=a(1:44,12);z1=a(1:44,7);
scatter3(x1,y1,z1,'k^')
hold on
x2=a(45:80,11);y2=a(45:80,12);z2=a(45:80,7);
scatter3(x2,y2,z2,'m^')
hold on
x3=a(45:80,11);y3=a(45:80,12);z3=a(45:80,7);
scatter3(x3,y3,z3,'r^')
hold on
x4=a(81:146,11);y4=a(81:146,12);z4=a(81:146,7);
scatter3(x4,y4,z4,'c^')
hold on
x5=a(147:284,11);y5=a(147:284,12);z5=a(147:284,7);
scatter3(x5,y5,z5,'g^')
hold on
shading interp
title('Ni 离子浓度分布');
zlabel('Ni (μg/g)')
legend('地形图','生活区','工业区','山区','交通区','公园路绿地区','Location','Best');  %图例

运行效果：热力图和散点图在一起，也可以自行分开