detect_ignition/new_likelihood/make_heats.m

   1 function [ heat_map ] = make_heats( tign , t_now )\r
   2 % [ heat_map ] = make_heats( tign )\r
   3 % function returns an array of heat output in pixels\r
   4 %inputs\r
   5 %   tign : matrix with fire arrival times for a simulation\r
   6 %   t_now : present time, (time of satellite imaging)\r
   7 %output\r
   8 %   heap_map : matrix with heat in each pixel, based on exponential decay\r
   9 %       with time since fire arrival\r
  10 \r
  11 [m n] = size(tign);\r
  12 heat_map = zeros(m,n);\r
  13 \r
  14 %const\r
  15 T_0 = 1;\r
  16 %decay = -0.3 --> 1.1% of heat output at 15 hours\r
  17 decay = -0.3; % ?? -->  %set to10 for patch  ,0.02 for hill2 test\r
  18 for ii = 1:m\r
  19     for jj = 1:n\r
  20         if tign(ii,jj) < t_now\r
  21             heat_map(ii,jj) = 0;\r
  22         else\r
  23             heat_map(ii,jj) = T_0*exp((tign(ii,jj)-t_now)*decay);\r
  24         end\r
  25     end\r
  26 end\r
  27 \r
  28 %figure\r
  29 %contour(heat_map)\r
  30 \r
  31 end\r
  32 \r