11年国赛A题土壤中重金属分布的二维与三维分布图_matlab插值

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的空间分布特征')