散乱插值及MATLAB绘图函数

转载自：https://blog.csdn.net/helloworld_wj/article/details/77511007

常见的散乱插值方法有Kriging（克里金）插值，Shepard插值等

这里使用了matlab中的工具直接进行插值并绘图，示例代码如下：

clc,clear
A=xlsread('C:\Users\Administrator\Desktop\集训\8.23\3. 2011~A资料\2011A~题目\data.xls','附件1','A4:E322');
C=xlsread('C:\Users\Administrator\Desktop\集训\8.23\3. 2011~A资料\2011A~题目\data.xls','附件2','A4:I322');
x=A(:,2);y=A(:,3);z=A(:,4);qy=A(:,5);
As=C(:,2);Cd=C(:,3);Cr=C(:,4);Cu=C(:,5);
Hg=C(:,6);Ni=C(:,7);Pb=C(:,8);Zn=C(:,9);

[X,Y,Z]=griddata(x,y,z,linspace(0,28654)’,linspace(0,18449),‘v4’);
figure,mesh(X,Y,Z)
title(‘采样地形图’);

[X,Y,AS]=griddata(x,y,As,linspace(0,28654)’,linspace(0,18449),‘v4’);
figure,mesh(X,Y,AS)
title(‘As浓度随平面地形分布图’);

[X,Y,CD]=griddata(x,y,Cd,linspace(0,28654)’,linspace(0,18449),‘v4’);
figure,mesh(X,Y,CD)
title(‘Cd浓度随平面地形分布图’);

[X,Y,CR]=griddata(x,y,Cr,linspace(0,28654)’,linspace(0,18449),‘v4’);
figure,mesh(X,Y,CR)
title(‘Cr浓度随平面地形分布图’);

[X,Y,CU]=griddata(x,y,Cu,linspace(0,28654)’,linspace(0,18449),‘v4’);
figure,mesh(X,Y,CU)
title(‘Cu浓度随平面地形分布图’);

[X,Y,NI]=griddata(x,y,Ni,linspace(0,28654)’,linspace(0,18449),‘v4’);
figure,mesh(X,Y,NI)
title(‘Ni浓度随平面地形分布图’);

[X,Y,HG]=griddata(x,y,Hg,linspace(0,28654)’,linspace(0,18449),‘v4’);
figure,mesh(X,Y,HG)
title(‘Hg浓度随平面地形分布图’);

[X,Y,PB]=griddata(x,y,Pb,linspace(0,28654)’,linspace(0,18449),‘v4’);
figure,mesh(X,Y,PB)
title(‘Pb浓度随平面地形分布图’);

[X,Y,ZN]=griddata(x,y,Zn,linspace(0,28654)’,linspace(0,18449),‘v4’);
figure,mesh(X,Y,ZN)
title(‘Zn浓度随平面地形分布图’);

[X,Y,QY]=griddata(x,y,qy,linspace(0,28654)’,linspace(0,18449),‘v4’);
figure,surf(X,Y,QY)
title(‘区域分布图’);

• 1
• 2
• 3
• 4
• 5
• 6
• 7
• 8
• 9
• 10
• 11
• 12
• 13
• 14
• 15
• 16
• 17
• 18
• 19
• 20
• 21
• 22
• 23
• 24
• 25
• 26
• 27
• 28
• 29
• 30
• 31
• 32
• 33
• 34
• 35
• 36
• 37
• 38
• 39
• 40
• 41
• 42
• 43
• 44
• 45
• 46

目标为已知坐标、海拔、地区类型以及各重金属元素含量的319个调查点。

各二维三维绘图函数在《精通MATLABR2014b》中有简略介绍，也可以自行上网搜索。