11年国赛A题土壤中重金属的二维与三维分布图_matlab插值
更多内容请关注公众号:狗头的成长日记
三维作图
A=xlsread('二维','sheet1'); %读取excel文件中代表该类元素的列向量
x1=0:180:28654;
y1=0:120:18449;
x=A(:,1); %读取x坐标值列向量
y=A(:,2); %读取y坐标值列向量
z=A(:,4); %读取文件中代表该类元素的浓度列向量
[x1,y1]=meshgrid(x1,y1); %生成坐标矩阵
z1=griddata(x,y,z,x1,y1,'v4'); %插值
meshc(x1,y1,z1) %画带有等高线的三维分布图
title('图3 城市土壤Cd的空间分布特征 ')
二维作图
A=xlsread('二维','sheet1'); %读取excel文件
x=A(:,1);
y=A(:,2);
[x,y]=meshgrid(0:100:29000,100:30:19000); %建立坐标矩阵
z3=griddata(A(:,1),A(:,2),A(:,3),x,y,'v4'); %做出二维地势分布图
%subplot(2,2,1) %将图片放在2*2图片窗口上的早上角
figure
[c,h] =contour(x,y,z3); %绘制Z3等高线
clabel(c,h)
ylabel('y');
title('图2 城市土壤As的空间分布特征');
z4=griddata(A(:,1),A(:,2),A(:,4),x,y,'v4');
%subplot(2,2,2)
figure
[c,h] =contour(x,y,z4);
clabel(c,h)
ylabel('y');
title('图4 城市土壤Cd的空间分布特征')
z5=griddata(A(:,1),A(:,2),A(:,5),x,y,'v4');
%subplot(2,2,3)
figure
[c,h] =contour(x,y,z5);
clabel(c,h)
ylabel('y');
title('图6 城市土壤Cr的空间分布特征')
z6=griddata(A(:,1),A(:,2),A(:,6),x,y,'v4');
%subplot(2,2,4)
figure
[c,h] =contour(x,y,z6);
clabel(c,h)
ylabel('y');
title('图8 城市土壤Cu的空间分布特征')
figure
z7=griddata(A(:,1),A(:,2),A(:,7),x,y,'v4');
%subplot(2,2,1)
[c,h] =contour(x,y,z7);
clabel(c,h)
ylabel('y');
title('图10 城市土壤Hg的空间分布特征')
z8=griddata(A(:,1),A(:,2),A(:,8),x,y,'v4');
%subplot(2,2,2)
figure
[c,h] =contour(x,y,z8);
clabel(c,h)
ylabel('y');
title('图12 城市土壤Ni的空间分布特征')
z9=griddata(A(:,1),A(:,2),A(:,9),x,y,'v4');
%subplot(2,2,3)
figure
[c,h] =contour(x,y,z9);
clabel(c,h)
ylabel('y');
title('图14 城市土壤Pb的空间分布特征')
z10=griddata(A(:,1),A(:,2),A(:,10),x,y,'v4');
%subplot(2,2,4)
figure
[c,h] =contour(x,y,z10);
clabel(c,h)
ylabel('y');
title('图16 城市土壤Zn的空间分布特征')