Question

I have some particles that I’ve identified in a larger image, and need to parse into smaller images for each particle. I’ve used the regionprops ‘BoundingBox’ function, but haven’t been successful yet. How can I now make a rectangular subimage of image I using BoundingBox? I can use BoundingBox to draw a rectangle on the original image, but the parameters returned by BoundingBox seem not to be of pixel dimension (x,y, width, height), (x1, y1, x2, y2), etc, which I would expect a bounding box to return. I’ve written some sample code using coins.png to make it easier for anyone to understand. Can you help me with this? Thanks!

figure(1);
I = imread('coins.png');
bw = im2bw(I, graythresh(I));
bw2 = imfill(bw,'holes');
imshow(bw2);


figure(2);
L = bwlabel(bw2);
imshow(label2rgb(L, @jet, [.7 .7 .7]))

figure(3);
imshow(I);
s = regionprops(L, 'BoundingBox');
rectangle('Position', s(1).BoundingBox);

Answer 1

The parameters returned by regionprops are [y,x,width,height] in matrix coordinates (see also “unexpected Matlab”.

Thus, to extract the rectangle, you write:

subImage = I(round(s(1).BoundingBox(2):s(1).BoundingBox(2)+s(1).BoundingBox(4)),...
       round(s(1).BoundingBox(1):s(1).BoundingBox(1)+s(1).BoundingBox(3)));