Wählen Sie das größte Objekt in einem Bild

Question

Ich versuche, das größte Objekt in einem Bild zu finden und alle anderen Objekte im Bild zu entfernen, die kleiner als das Bild sind.

Dies ist, was ich habe, aber ich kann es nicht zur Arbeit bringen.

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, []); 

Ich arbeite mit digitalen Mammogrammen wie these. Ich versuche alle Objekte aus dem Bild zu entfernen, mit Ausnahme der Brustregion.

Answer 1

Verwenden bwconncomp statt, da es die Koordinatenindizes für Region in einer separaten Zelle zurückkehrt, wo die Größe von jedem leicht zu erkennen ist:

>> 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]} 

Das PixelIdxList Feld ist ein Zellenfeld mit dem Indizes der Koordinaten für jeden Region. Die Länge jeder Anordnung ist die Größe der jeweiligen Region:

>> numPixels = cellfun(@numel,CC.PixelIdxList) 
numPixels = 
    1  2 
>> [biggestSize,idx] = max(numPixels) 
biggestSize = 
    2 
idx = 
    2 

Dann können Sie einfach ein neues Bild machen mit nur diese Komponente:

BW2 = false(size(BW)); 
BW2(CC.PixelIdxList{idx}) = true; 

---

EDIT: Aus den Kommentaren, die Notwendigkeit, das Ausgabebild zu beschneiden, so dass die Region an die Ränder kommt, kann mit regionprops unter Verwendung der Option "BoundingBox" adressiert werden:

s = regionprops(BW2, 'BoundingBox'); 

das gibt Ihnen ein Rechteck s.BoundingBox, die Sie verwenden können, um mit BW3 = imcrop(BW2,s.BoundingBox); zu beschneiden.

Answer 2

Wenn Sie mit dem bwlabel Ansatz fortsetzen möchten, können Sie diese verwenden -

-Code

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) 

Ausgabe