makegrid.m

麻省理工学院的人工智能工具箱,很珍贵,希望对大家有用!

function RP_Grid=makeGrid(data,Maes)

%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);

%nacional10
x=[-10:.5:5];
y=[44:-.5:35];
%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