cycling/filled_contours.m

   1 function filled_countours(ps,an)
   2 % plots filled contours for forecast and analysis times in seconds
   3
   4
   5
   6 tmax = max(ps.points(:,3));
   7 tmax = ps.red.end_datenum;
   8 tmin = min(min(ps.red.tign_g(:)),min(an(:)));
   9 clines = linspace(tmin,tmax,20);
  10
  11 if isfield(ps.red,'red') && length(ps.red.fxlong) >2000
  12     figure(190),contourf(ps.red.red.fxlong(1:10:end,1:10:end),ps.red.red.fxlat(1:10:end,1:10:end) ...
  13         ,ps.red.red.tign(1:10:end,1:10:end),'k')
  14 else
  15     %figure,contourf(ps.red.fxlong,ps.red.fxlat,ps.red.tign_g,clines,'k')
  16     figure(190),contourf(ps.red.fxlong,ps.red.fxlat,ps.red.tign_g,'k')
  17 end
  18 hold on
  19
  20 points_set = ps.points(:,3)<=tmax;
  21 scatter(ps.points(points_set,2),ps.points(points_set,1),'*r')
  22 xlabel('Lon [degree]'),ylabel('Lat [degre]')
  23 title('Forecast and active fire detections')
  24 %title('"Ground Truth" and artificial fire detections')
  25 hold off
  26
  27
  28 %figure,contourf(ps.red.fxlong,ps.red.fxlat,an,clines,'k')
  29 figure(191),contourf(ps.red.fxlong,ps.red.fxlat,an,'k')
  30 hold on
  31 scatter(ps.points(:,2),ps.points(:,1),'*r')
  32 xlabel('Lon [degree]'),ylabel('Lat [degre]')
  33 title('Analysis and active fire detections')
  34 %title('Estimate and artificial fire detections')
  35 hold off
  36
  37 end