使用gridata插值拟合
clear;
data=xlsread('F:\data1');
x=data(:,1);
y=data(:,2);
z=data(:,3);
a=data(:,4);
As=data(:,5);
Cd=data(:,6);
Cr=data(:,7);
Cu=data(:,8);
Hg=data(:,9);
Ni=data(:,10);
Pb=data(:,11);
Zn=data(:,12);
[X,Y,Z]=griddata(x,y,z,linspace(0,28654),linspace(0,28654)','v4');%插值
[X,Y,A]=griddata(x,y,a,linspace(0,28654),linspace(0,28654)','nearest');%插值
[X,Y,As_1]=griddata(x,y,As,linspace(0,28654),linspace(0,28654)','v4');%插值
% [X,Y,Cd_1]=griddata(x,y,Cd,linspace(0,28654),linspace(0,28654)','v4');%插值
% [X,Y,Cr_1]=griddata(x,y,Cr,linspace(0,28654),linspace(0,28654)','v4');%插值
% [X,Y,Cu_1]=griddata(x,y,Cu,linspace(0,28654),linspace(0,28654)','v4');%插值
% [X,Y,Hg_1]=griddata(x,y,Hg,linspace(0,28654),linspace(0,28654)','v4');%插值
% [X,Y,Ni_1]=griddata(x,y,Ni,linspace(0,28654),linspace(0,28654)','v4');%插值
% [X,Y,Pb_1]=griddata(x,y,Pb,linspace(0,28654),linspace(0,28654)','v4');%插值
% [X,Y,Zn_1]=griddata(x,y,Zn,linspace(0,28654),linspace(0,28654)','v4');%插值
%
%scatter(x,y,5,z);
%title('散点图');
%figure
%pcolor(X,Y,A);shading interp%伪彩色图
%title('伪彩色图');
contourf(X,Y,A)
grid on
title('区域图');
colorbar();
xlabel('1-生活区 2-工业区 3-山区 4-交通区 5-公园绿地');
figure,contourf(X,Y,Z)%等高线图
colorbar();
title('地形等高线图')
grid on
figure,surf(X,Y,Z)
view(-20,45);
colorbar(); |