detection/showmod14.m

   1 function showmod14(v,varargin)
   2 % showmod14(v,alpha,string)
   3 % display MODIS14 data
   4 % from loading file converted by geotiff2mat.py
   5 % input
   6 %   v.data     MODIS 14 data
   7 %   v.geo      geotransform from geotiff
   8 %   v.file    optional title string
   9 %   v.axis    min, max longgitude. latitude
  10
  11 % newplot
  12
  13 if length(varargin)>0,
  14     alpha=varargin{1};
  15 else
  16     alpha=1;
  17 end
  18 pixels_fire=[sum(v.data(:)==7),sum(v.data(:)==8),sum(v.data(:)==9)];
  19 if length(varargin)>1,
  20     t=varargin{2};
  21 else
  22     t = sprintf('%s %s fire pixels %i %i %i',v.file,...
  23     datestr(v.time,'yyyy-mm-dd HH:MM:SS'),pixels_fire);
  24 end
  25 cmap=cmapmod14;
  26 alphadata=zeros(size(v.data));
  27 alphamap=any(cmap,2)*alpha;
  28 % alphamap=ones(size(alphamap));
  29 d=v.data;
  30 d(isnan(d))=0;
  31 a=alphamap(d+1);
  32 image('Xdata',[v.lon(1),v.lon(end)],'Ydata',[v.lat(1),v.lat(end)],...
  33     'Cdata',v.data,'Alphadata',a);
  34 a=gca;set(a,'Ydir','normal')
  35 colormap(cmap);
  36 xlabel('Longitude (deg)')
  37 ylabel('Latitude (deg)')
  38 title(t,'Interpreter','none')
  39 if isfield(v,'axis'),
  40     axis(v.axis);
  41 end
  42 daspect([1,cos(mean(v.lat(:))*pi/180),1]);
  43
  44 % grid on
  45 % show fire pixels
  46 % hold on, [i,j,c]=find(data > 6); plot(lat(j),lon(i),'+r'), hold off
  47 end
  48