makegrid_b.m

来自「麻省理工学院的人工智能工具箱,很珍贵,希望对大家有用!」· M 代码 · 共 31 行

function RP_Grid=makeGrid_b(data,loc,x,y)

%load ('..//AreaStations/TP_rellenados/Nacional/RMaesPre')
%load ('..//AreaStations//TP_rellenados/Nacional/RdataPre')
%data=data(1:5538,:);

%dum1=Maes(:,140:146);
%dum1(find(dum1==' '))='0';
%loc(:,1)=(str2num(dum1(:,1:2))+str2num(dum1(:,3:4))/60+...
%   str2num(dum1(:,5:6))/3600).*((-1).^(str2num(dum1(:,7))+1));
%dum1=Maes(:,147:152);
%dum1(find(dum1==' '))='0';
%loc(:,2)=(str2num(dum1(:,1:2))+str2num(dum1(:,3:4))/60+...
%   str2num(dum1(:,5:6))/3600);

%a cada punto del grid le asignamos los observatorios mas cercanos
RP_Grid=single(1);
RP_Grid(1:size(data,1),1:length(y)*length(x))=0;
k=0;
for j=1:length(y)
   for i=1:length(x)
      k=k+1;
      g=[x(i) y(j)];
      [lista,dd]=MLknn(g,loc,10,2);
      dd=1./dd.^2;dd=dd/sum(dd);
      b=repmat(dd,size(data(:,lista),1),1);
      a=nansum(data(:,lista).*b,2);
      RP_Grid(:,k)=a;  
   end
end