detect_ignition/new_likelihood/test3.m

   1 %test 3   -   put reported detection on the fire line, move the actual\r
   2 %location away\r
   3 \r
   4 % second test script looking for tails???\r
   5 \r
   6 \r
   7 domain_size = 1000;\r
   8 \r
   9 \r
  10 %create fire mask\r
  11 mask = zeros(domain_size,domain_size);\r
  12 strip = ones(domain_size,1);\r
  13 mask(:,500) = strip;\r
  14 \r
  15 %make heat map\r
  16 heat = zeros(domain_size,domain_size);\r
  17 heat_strip = zeros(1,domain_size);\r
  18 a = 0.02;\r
  19 for i = 500:domain_size\r
  20     heat_strip(1,i) = exp(-(i-500)*a);\r
  21 end\r
  22 for i = 1:domain_size\r
  23     heat(i,:) = heat_strip;\r
  24 end\r
  25 %mesh(heat)\r
  26 contour(heat)\r
  27 %make probability matrix\r
  28 detection_probs = zeros(domain_size,domain_size);\r
  29 \r
  30 for i = 1:domain_size\r
  31     for j = 1:domain_size\r
  32         detection_probs(i,j) = detection_probability(heat(i,j));\r
  33     end\r
  34 end\r
  35 \r
  36 %move actual pixel across domain, passing through location of reported\r
  37 %location\r
  38 \r
  39 radius = 50;\r
  40 like = [];\r
  41 for i = 1:domain_size\r
  42     dist2 = (500-i)^2;\r
  43     l = exp(-dist2/(3*radius)^2)*detection_probs(500,i);\r
  44     like(i) = log(l);\r
  45 end\r
  46 \r
  47 close all\r
  48 plot(like)\r