detection/cmapmod14.m

   1 function [cmap,varargout]=cmapmod14
   2 % cmap=cmapmod14
   3 % colormap for active fires detection level 2 MOD14 files and derived
   4 % black = notprocessed
   5 % blue  = water
   6 % gray  = cloud
   7 % green = land
   8 % red   = fire (saturation is confidence level)
   9
  10 imax=20;
  11 cpap=zeros(7+3*(imax+1),3);
  12 cmap(1:7,:)= [ ...
  13     0   0   0      %0 not processed (missing input data), black
  14     0   0   0      %1 not used, black
  15     0   0   0      %2 not processed (other reason)
  16     0   0   0.2    %3 water, dark blue
  17     0.2 0.4 0.4    %4 cloud, purple gray
  18     0   0.3   0    %5 non-fire clear land, green
  19     0   0   0      %6 unknown
  20 ];
  21
  22  %7 low-confidence fire
  23  %8 nominal-confidence fire
  24  %9 high-confidence fire
  25  % then repeat transition to old
  26 red=[1 0 0];
  27 yellow=[1  1  0];
  28
  29 atten=5;
  30 for i=0:imax
  31     w=exp(-atten*i/imax)-exp(-atten);
  32     col=(1-w)*yellow+w*red;
  33     cmap(8+3*i  ,:)=0.4*col;
  34     cmap(8+3*i+1,:)=0.7*col;
  35     cmap(8+3*i+2,:)=1*col;
  36 end
  37 if size(cmap,1)>256,
  38     error('imax too large')
  39 end
  40 if nargout>1,
  41     varargout(1)={imax};
  42 end