我试图找到图像中最大的对象,并删除图像中比它小的任何其他对象。
这是我所拥有的,但我无法让它工作。
l=bwlabel(BW);
%the area of all objects in the image is calculated
stat = regionprops(l,'Area','PixelIdxList');
[maxValue,index] = max([stat.Area]);
%remove any connected areas smaller than the biggest object
BW2=bwareaopen(BW,[maxValue,index],8);
subplot(5, 5, 4);
imshow(BW2, []);
我正在使用数字乳房X光检查,例如这些。我正在尝试从图像中删除除乳房区域之外的所有对象。
bwconncomp
代替,因为它返回单独单元格中区域的坐标索引,其中每个单元格的大小很容易辨别:
>> BW = [1 0 0; 0 0 0; 0 1 1]; % two regions
>> CC = bwconncomp(BW)
CC =
Connectivity: 8
ImageSize: [3 3]
NumObjects: 2
PixelIdxList: {[1] [2x1 double]}
PixelIdxList
字段是一个元胞数组,其中包含每个区域的坐标索引。 每个数组的长度就是每个区域的大小:
>> numPixels = cellfun(@numel,CC.PixelIdxList)
numPixels =
1 2
>> [biggestSize,idx] = max(numPixels)
biggestSize =
2
idx =
2
然后你就可以轻松地用这个组件制作一个新图像:
BW2 = false(size(BW));
BW2(CC.PixelIdxList{idx}) = true;
编辑:根据评论,需要裁剪输出图像以使区域到达边缘可以使用“BoundingBox”选项通过
regionprops
来解决:
s = regionprops(BW2, 'BoundingBox');
它会给你一个矩形
s.BoundingBox
,你可以用它来裁剪BW3 = imcrop(BW2,s.BoundingBox);
。
如果您想继续使用
bwlabel
方法,您可以使用此 -
代码
BW = im2bw(imread('coins.png')); %%// Coins photo from MATLAB Library
[L, num] = bwlabel(BW, 8);
count_pixels_per_obj = sum(bsxfun(@eq,L(:),1:num));
[~,ind] = max(count_pixels_per_obj);
biggest_blob = (L==ind);
%%// Display the images
figure,
subplot(211),imshow(BW)
subplot(212),imshow(biggest_blob)
输出