cycling/plot_paths2d.m

   1 function plot_paths2d(path_struct)
   2 %plots paths in 2d setting
   3 % input: path_struct - structure with graph, paths, red, etc...
   4 %        path_struct = graph_dets(w,cull)
   5 r = path_struct.red;
   6
   7 % draw filled contour of fuels
   8 fuel_map = r.nfuel_cat;
   9 %fuel_map = simp_fuels(fuel_map);
  10 figure,contourf(r.fxlong,r.fxlat,fuel_map,'LineWidth',0.5);
  11 colormap summer;
  12
  13 %draw paths on figure
  14 paths = path_struct.paths;
  15 pts = path_struct.points;
  16 hold on
  17 for i = 1:length(paths)
  18     p = paths(i).p;
  19     plot(pts(p,2),pts(p,1),'r','LineWidth',2);
  20 end %for i
  21 hold off
  22
  23 %plot elevation data and paths
  24
  25 % draw filled contour of fuels
  26 elev = r.fhgt;
  27 %fuel_map = simp_fuels(fuel_map);
  28 figure,contour(r.fxlong,r.fxlat,elev,'k','LineWidth',0.5,'ShowText','on');
  29
  30 %draw paths on figure
  31 paths = path_struct.paths;
  32 pts = path_struct.points;
  33 hold on
  34 for i = 1:length(paths)
  35     p = paths(i).p;
  36     plot(pts(p,2),pts(p,1),'r','LineWidth',2);
  37 end %for i
  38 hold off
  39
  40 %plot paths with gradient of surface
  41 %smooth surface first
  42 r.tign = imgaussfilt(r.tign,1/2)
  43 [dx,dy]= fire_gradients(r.fxlong,r.fxlat,r.tign,1);
  44
  45 figure,contour(r.fxlong,r.fxlat,r.tign)
  46 hold on
  47 quiver(r.fxlong,r.fxlat,dy,dx)
  48 for i = 1:length(paths)
  49     p = paths(i).p;
  50     plot(pts(p,2),pts(p,1),'r','LineWidth',1);
  51 end %for i
  52 hold off
  53
  54
  55 end %function