cycling/forecast_ellipse.m

   1 function e = forecast_ellipse(r,t)
   2
   3 %function e = detection_ellipse(r,t)
   4 % function returns eigenvectors for ellipse containg fire detections
   5 % inputs;
   6 %   r        struct, fire structure from subset_domain
   7 %   t        scalar, t = g(n).time, time of final sateliite data
   8 %output:
   9 %   e        struct with eigenvectors for ellipse contining forecast
  10 %               perimeter
  11
  12 f_mask = r.tign < t;
  13 %fprintf('Forecast are ~= %d \n',sum(f_mask(:)));
  14 f_lons = r.fxlong(f_mask);
  15 f_lats = r.fxlat(f_mask);
  16 forecast = [f_lons,f_lats];
  17
  18 ci = 1.96;
  19 e = ellipse_fit(forecast,ci);
  20 e.area = sum(f_mask(:));
  21 e.ellipse_area = det(e.d);
  22
  23
  24 end
  25
  26