cycling/smooth_up.m

   1 function sm_up = smooth_up(tign,a,b)
   2 %applies gradual blurring at different contours of tign
   3 % a,b are optional parameters
   4
   5 t0 = min(tign(:));
   6 t1 = max(tign(:));
   7
   8 steps = 30;
   9 t = linspace(t0,t1,steps);
  10 [m,n]= size(tign);
  11
  12 %compute different std. for imgausfilt
  13 %mx = 2;%
  14 if ~exist('a','var')
  15     a = 300;
  16     b = 1/3;
  17 end
  18 mx = max(1.0,max(m,n)/a);
  19 %mn = 1/2;
  20 mn = b*mx;
  21 fprintf('Smooth paramaters : %f   %f   \n',mn,mx);
  22 m = (mn-mx)/(steps-2);
  23 y = @(t) mx + m*(t-2);
  24
  25 %keep track of the change in ignition time
  26 t_diff(1) = 0;
  27
  28 t_temp = tign;
  29 % figure(136),mesh(lon,lat,tign),title('Original')
  30 % figure(137),mesh(lon,lat,tign),title('Smoothing')
  31 for i = 2:steps
  32     m1 = tign < t(i-1);
  33     m2 = tign < t(i);
  34     msk = logical(m2-m1);
  35     t_blur = imgaussfilt(tign,y(i));
  36     t_temp(msk) = t_blur(msk);
  37     %figure(138),contour(t_temp,steps,'k')
  38     %pause(5/steps)
  39     t_diff(i) = min(t_temp(:))-t0;
  40     t_temp(m2) = t_temp(m2)-i/steps*t_diff(i);
  41     %t_temp = t_temp-i/steps*t_diff(i);
  42 end
  43     %shift down
  44 %     pixel_time_diff = t_temp-tign;
  45 %     time_slope = pixel_time_diff./t_diff(i);
  46 %     t_temp3 = t_temp - time_slope.*pixel_time_diff;
  47 %     t_shift = max(t_temp3,t_temp);
  48 %figure,plot(t_diff),title('Change in tign')
  49 sm_up = t_temp;