1.  
  2. % load background image
  3.  
  4. directory = '/group/project/VISION/web/AV2/DATA2/frame';
  5.  
  6.     data = [ directory, 'bg' , '.jpg'];
  7.     bgim = imread(data, 'jpeg'); %import image
  8.  
  9.  
  10. %magic numbers
  11.  
  12. threshold = 13;
  13. opening = 5;
  14. connopening = 200;
  15.  
  16. for i = 1:146
  17.  
  18.     im = latest_im(directory, i);
  19.  
  20.     subim = abs(bgim -im);
  21.     binim = subim > threshold;
  22.     binim = binim(:,:,1) | binim(:,:,2) | binim(:,:,3);
  23.  
  24.     erodedim = bwmorph(binim,'open', opening);
  25.     erodedim = bwmorph(binim,'close', opening);
  26.  
  27.  
  28.     erodedim = bwareaopen(erodedim,connopening);
  29.     labeledim = bwlabel(erodedim, 8);
  30.  
  31. %     bwim = im2bw(subim, threshold);
  32. figure(1);    
  33. imagesc(labeledim);
  34.  
  35. pause(0.1);
  36.  
  37. end