1、%function Z=IDW(data,station,lon,lat,e,num)%data1*n%station data n4*2 %lon lat ()1*n %e 2% num 412 %e=2num=12arcgis% DATA=IDW2015(P_m,latlon,lon,lat,2,12);%n1 n2=size(lon);n3 n4=size(data);% if station(1,2)station(1,1) %x x=station(:,2); % y=station(:,1); % else% y=station(:,2); %wei% x=station(:,1)
2、; %jing% endz=data; %for i1=1:n1 X=lon(i1,:); % Y=lat(i1,:); % t1=length(x); t2=length(X); if(t1=0&t2=0) %r for j=1:t2 r=sqrt(X(j)-x).2+(Y(j)-y).2).e; I J=sort(r); m=length(I); if(mnum) 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);% if isnan(Z(:,j)% Z(:,j)=z(J(1);% end end endend