cycling/combine_rans.m

   1 function c = combine_rans(n,w)
   2
   3 c = ran_cone(w);
   4 for i = 1:n
   5     c1 = ran_cone(w);
   6     a1 = sum(sum(c1<max(c1(:))));
   7     a = sum(sum(c<max(c(:))));
   8     c = a/(a+a1)*c+a1/(a+a1)*c1;
   9     %c = max(c,c1);
  10     %c = 1/2*(c+c1);
  11 end
  12
  13 %c = imgaussfilt(c,1);
  14 c = smooth_up(c);
  15 % figure(123),mesh(c);
  16 %figure(124),contour(c,20,'k')
  17
  18 t_min = min(w.tign_g(:));
  19 t_max = max(w.tign_g(:));
  20 c_min = min(c(:));
  21 c_max = max(c(:));
  22
  23 % days = round((t_max-t_min)/24/3600);
  24 % stepper = (c_max-c_min)/days;
  25 % for i = 1:days
  26 %    c_min = min(c(:));
  27 %    c_max = max(c(:));
  28 %    msk =  c>c_min+i*stepper;
  29 %    figure(217),scatter(w.fxlong(msk),w.fxlat(msk))
  30 %    pause(3)
  31 %    sum(msk(:));
  32 %    m_rand = 1+(-1)^i/2
  33 %    %m_rand = 1/10*(1+sin((i+1)*pi*rand/2))^2;
  34 %    c(~msk)= m_rand/10*c(~msk);
  35 % end
  36
  37 c_min = min(c(:));
  38 c_max = max(c(:));
  39
  40 c = t_min + (t_max-t_min)/(c_max-c_min)*(c-c_min);
  41 %figure(123),mesh(w.fxlong,w.fxlat,c);
  42 figure(124),contour(w.fxlong,w.fxlat,c,20,'k')
  43
  44 %quick_mesh(c)
  45
  46
  47 end%function