a b/functions/clean_segs.m 

  1
  
function t_yLV_clean=clean_segs(t_yLVh)





  
  

  2
  
%%





  
  

  3
  
for k=1:size(t_yLVh,3)





  
  

  4
  
    





  
  

  5
  
    BW=t_yLVh(:,:,k);





  
  

  6
  
    BW2= imfill(BW,'holes');





  
  

  7
  
    cc = bwconncomp(BW2,4);





  
  

  8
  
    stats = regionprops(cc, 'Area');





  
  

  9
  
    [area_max,idx] = max([stats.Area]);





  
  

  10
  
    temp= ismember(labelmatrix(cc), idx);





  
  

  11
  
    t_yLV_clean(:,:,k) = imfill(temp,'holes');





  
  

  12
  
    t_yLV_clean(:,:,k)=bwconvhull(t_yLV_clean);





  
  

  13
  
end





  
  

  14
  





  
  

  15