反距离加权插值matlab_精品文档.txt
- 文档编号:1520207
- 上传时间:2022-10-22
- 格式:TXT
- 页数:2
- 大小:1.41KB
反距离加权插值matlab_精品文档.txt
《反距离加权插值matlab_精品文档.txt》由会员分享,可在线阅读,更多相关《反距离加权插值matlab_精品文档.txt(2页珍藏版)》请在冰豆网上搜索。
%�������Ȩ��ֵ�����Զ�����������
functionZ=IDW(data,station,lon,lat,e,num)
%ע��dataΪ1*n������
%stationdata��Ӧ����ľ�γ��n4*2��һ��γ���ڶ��о���
%lonlat��Ҫ��ֵ���ľ�γ�����ݸ�ʽ���ŷ�(���ŷţ����ŷ�~)��1*n
%e��������Ϊ2
%������ѡ��numΪ���ֵ����ĵ�����Ϊ4~12
%e=2��num=12Ϊarcgis��Ĭ��ֵ
%���ӣ��������ֵ
%DATA=IDW2015(P_m,latlon,lon',lat',2,12);
%
[n1n2]=size(lon);
[n3n4]=size(data);
%ifstation(1,2)>station(1,1)%%%%%%%%ֻ�������ض�����Ŀ�ģ���֤����Ϊx
x=station(:
2)';%��֪���꾭��
y=station(:
1)';%��֪����γ��
%else
%y=station(:
2)';%��֪����wei��
%x=station(:
1)';%��֪����jing��
%end
z=data;%��֪���꺯��ֵ
fori1=1:
n1
X=lon(i1,:
);%Ҫ��ֵ�ľ���
Y=lat(i1,:
);%Ҫ��ֵ��γ��
t1=length(x);
t2=length(X);
if(t1~=0&t2~=0)
%�������r��������ֵ
forj=1:
t2
r=sqrt((X(j)-x).^2+(Y(j)-y).^2).^e;
[IJ]=sort(r);
m=length(I);
if(m>num)
R=I(1:
num);
z11=z(:
J);
z1=z11(:
1:
num);
else
display('����������!
');
stop;
end
a=sum(z1./repmat(R,n3,1),2);%%%gai
b=sum(1./R);
Z(:
j)=a./repmat(b,n3,1);
%ifisnan(Z(:
j))
%Z(:
j)=z(J
(1));
%end
end
end
end
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 距离 加权 matlab_ 精品 文档
![提示](https://static.bdocx.com/images/bang_tan.gif)