detection/plot_state_2d.m

   1 function plot_state_2d(fig,red,s,tign,obs,time_now)
   2
   3
   4 % [m,n] = size(tign);
   5 % if m*n > 500^2
   6 %     tign = tign(1:10:end,1:10:end);
   7 % end
   8 figure(fig),clf
   9 for i=1:length(obs)
  10     x=obs(i);
  11     kk=find(x.data(:)>=7);
  12     if ~isempty(kk),
  13         x.data=double(x.data);
  14         x.data(x.data<7)=NaN;
  15         x.data=x.data+3*floor((time_now-x.time)/0.25);
  16         showmod14(x,0.5)
  17         hold on
  18     end
  19 end
  20 if ~iscell(tign),
  21     tign={tign};
  22 end
  23 if ~iscell(s),
  24     s={s};
  25 end
  26
  27 color={'k','b','-.r'};
  28
  29
  30 for i=1:length(tign)
  31     lon = red.fxlong;
  32     lat = red.fxlat;
  33     z = tign{i};
  34     [m,n] = size(z);
  35     if m*n > 500^2
  36         z = z(1:10:end,1:10:end);
  37         lon = lon(1:10:end,1:10:end);
  38         lat = lat(1:10:end,1:10:end);
  39     end
  40     [c,h]=contour(lon,lat,z,[time_now,time_now],color{i});
  41     set(h,'linewidth',2)
  42 end
  43 h_legend=legend(s{:});
  44 set(h_legend,'FontSize',16);
  45 a=[red.min_lon,red.max_lon,red.min_lat,red.max_lat];
  46 axis manual
  47 axis(a)
  48 grid on
  49 hold off
  50 end