detect_ignition/new_likelihood/gauss_part.m

   1 function [g_part ] = gauss_part( fire_x, fire_y, pixel_x, pixel_y, sig )\r
   2 % compute gaussian part of the formula\r
   3 % inputs :\r
   4 %   fire_x, fire_y : [x y] location of pixel fire pixels (reported\r
   5 %   locations...)\r
   6 %   pixel_x, pixel_y : [x y] location of where satellite is really looking\r
   7 %   sig  :standard deviation\r
   8 %   \r
   9 % outputs:\r
  10 %   g_part: sum of the gaussians\r
  11 \r
  12 %r = 200.0;  %200 m pixel ->>> 300 m radius at 3 sigma uncertainty\r
  13 %sig = 1;\r
  14 \r
  15 dist2 = (fire_x-pixel_x)^2+(fire_y-pixel_y)^2;\r
  16 g_part = exp(-dist2/(2*sig^2));\r
  17 \r
  18 \r
  19 end\r