detection/plot_state.m

   1     function plot_state(fig,red,s,T,obs,c)
   2         % fig   figure number
   3         % T     variable displayed
   4         % s     title string
   5         % obs   observations structure
   6         % red   reduced structure, for bounds
   7         % c     contour heights
   8         [m,n] = size(T);
   9         if m*n > 500^2
  10             d=7;
  11             T = T(1:d:end,1:d:end);
  12             red.fxlong = red.fxlong(1:d:end,1:d:end);
  13             red.fxlat = red.fxlat(1:d:end,1:d:end);
  14         end
  15         base_time=red.min_tign;
  16         fprintf('Figure %i %s',fig,s)
  17         figure(fig), clf
  18         if exist('T') && ~isempty(T),
  19             mint=min(T(:));
  20             maxt=max(T(:));
  21             fprintf(' tign min %g max %g',mint,maxt)
  22             num_hours=(maxt-mint)*24;
  23             step_hours=max(6,round(num_hours/100));
  24             hours = linspace(mint,maxt,30);
  25             %hours=[mint:step_hours/24:maxt];
  26             % T(T(:)>max(T(:))-tol)=NaN;
  27             h=contour3(red.fxlong,red.fxlat,T-base_time,hours-base_time);
  28             % set(h,'EdgeAlpha',0,'FaceAlpha',0.5); % show faces onl
  29             hold on
  30         end
  31         if exist('c') && ~isempty(c),
  32             fprintf(' contour at %s',num2str(c))
  33             contour3(red.fxlong,red.fxlat,T-base_time,c-base_time,'k');
  34             hold on
  35         end
  36         if exist('obs') && ~isempty(obs)
  37             fire_pixels3d(obs,base_time)
  38         end
  39         disp_bounds=[red.min_lon,red.max_lon,red.min_lat,red.max_lat];
  40         a=[disp_bounds,red.min_tign-base_time-1,red.max_tign-base_time];
  41         axis manual
  42         axis(a)
  43         xlabel('Longitude'),ylabel('Latitude'),zlabel('Days')
  44         title(s)
  45         grid on
  46         drawnow
  47         fprintf('\n')
  48 end
  49