hist_af.m

   1 function pixels=hist_af(af)
   2 for i=1:10
   3     count(i)=sum(af(:)==i-1);
   4 end
   5 pixels.unknown= count(1)+count(2)+count(3)+count(7);
   6 pixels.water  = count(4);
   7 pixels.cloud  = count(5);
   8 pixels.land   = count(6);
   9 pixels.fire_low=count(8);
  10 pixels.fire_med=count(9);
  11 pixels.fire_high=count(10);
  12 pixels.total = prod(size(af));
  13
  14