detect_ignition/new_likelihood/test6.m

   1 function [] = test6()
   2
   3
   4 close all
   5
   6 %make splines
   7 fprintf('Making splines \n');
   8
   9 [p_like_spline,p_deriv_spline,n_deriv_spline] = make_spline(72,2000);
  10 save splines.mat p_like_spline p_deriv_spline n_deriv_spline
  11
  12 %make fire data
  13 cone_slope = 20;
  14 fire_cone = @(x,y) cone_slope*sqrt(( x.^2 + y.^2));
  15 g = 100;
  16 grid_size =2*g+1;
  17 position = linspace(-10,10,grid_size);
  18 [x,y]= meshgrid(position,position);
  19
  20 z = fire_cone(x,y);
  21 fire_top = z > 100;
  22 z(fire_top) = 100;
  23 mesh(x,y,z);
  24 hold on
  25 xlabel('x');
  26
  27
  28 slice_time = [47 49 51];
  29 fires = -1*ones(size(x));
  30
  31
  32 contour3(x,y,z,[49 49],'k')
  33 fires = 5*ones(size(x));
  34 num_pts = 1000;
  35 x_coords = 1+round(2*g*rand(1,num_pts));
  36 y_coords = 1+round(2*g*rand(1,num_pts));
  37 %figure,scatter(x_coords,y_coords);
  38 %make fire mask
  39 hold on
  40 radius = slice_time(2)/10;
  41 angle = linspace(0,2*pi);
  42 perim = zeros(length(angle),2);
  43 perim(:,1) = radius*cos(angle)';
  44 perim(:,2) = radius*sin(angle)';
  45 plot(perim(:,1),perim(:,2),'k')
  46
  47 for i = 1:num_pts
  48     u = x(x_coords(i),y_coords(i));
  49     v = y(x_coords(i),y_coords(i));
  50     zt = norm([u v]);
  51     if abs(zt - radius) < 2  && zt < radius %49
  52         fires(x_coords(i),y_coords(i)) = 1;
  53         scatter(u,v,'r*');
  54     else
  55         if rand < 0.98
  56             fires(x_coords(i),y_coords(i)) = -1;
  57             scatter(u,v,'b');
  58         else
  59             fires(x_coords(i),y_coords(i)) = 1;
  60             scatter(u,v,'r*');
  61         end
  62     end
  63 end
  64
  65 %make all fire detections / ground
  66 fires = -ones(grid_size,grid_size);
  67
  68 % evaluate and plot likelihoods
  69 t = slice_time(2) - z;
  70
  71 new_like = input_num('Use new like? yes = 1',1);
  72 if new_like == 1
  73     [like,deriv]= evaluate_likes(fires,t,p_like_spline,p_deriv_spline,n_deriv_spline);
  74     %using old likelihood
  75 else
  76     stretch = [0.5,10,5,10];
  77     [like,deriv] = like2(fires,t,stretch);
  78 end
  79 fprintf('paused for plotting, etc... \n');
  80 figure,mesh(like),title('like')
  81 figure,mesh(deriv),title('deriv')
  82 figure,plot(t(100,:),like(100,:)),title('like slice')
  83 figure,plot(t(100,:),deriv(100,:)),title('deriv slice')
  84
  85
  86
  87 end % function
  88