other/Matlab/detection/fire_pixels.m

   1 function fire_pixels(obs,base_time)
   2 % fire_pixels(obs,base_time)
   3 % fire_pixels(obs,base_time,dim)
   4 % in
   5 %   obs         structure array of observations
   6 %   base_time   the start of simulations
   7 for i=1:length(obs)
   8     x=obs(i);
   9     kk=find(x.data(:)>=7);
  10     if ~isempty(kk),
  11         rlon=0.5*abs(x.lon(end)-x.lon(1))/(length(x.lon)-1);
  12         rlat=0.5*abs(x.lat(end)-x.lat(1))/(length(x.lat)-1);
  13         lon1=x.xlon(kk)-rlon;
  14         lon2=x.xlon(kk)+rlon;
  15         lat1=x.xlat(kk)-rlat;
  16         lat2=x.xlat(kk)+rlat;
  17         X=[lon1,lon2,lon2,lon1]';
  18         Y=[lat1,lat1,lat2,lat2]';
  19         Z=ones(size(X))*(x.time-base_time);
  20         cmap=cmapmod14;
  21         C=cmap(x.data(kk)'+1,:);
  22         C=reshape(C,length(kk),1,3);
  23         patch(X,Y,Z,C);
  24         hold on
  25     end
  26 end
  27 hold off
  28 end