0001 function cells = labelcell(bw, fr, minSize) 0002 0003 region = regionprops(bw,'Centroid','Area','Orientation','MajorAxisLength','MinorAxisLength'); 0004 0005 cells = horzcat(cat(1, region.Centroid), ... 0006 cat(1, region.Area), ... 0007 cat(1, region.Orientation),... 0008 cat(1, region.MajorAxisLength),... 0009 cat(1, region.MinorAxisLength)); 0010 0011 if ~isempty(cells) 0012 filterIndex = (cells(:,3) <= minSize); 0013 cells(filterIndex,:) = []; 0014 ids = 1:1:size(cells,1); 0015 frNum = fr * ones(size(cells,1),1); 0016 cells = horzcat(frNum, ids', cells); 0017 end 0018 0019 0020 0021 0022 end