¿Cómo obtener una subimagen rectangular desde regionprops (Imagen, 'BoundingBox') en Matlab?

Question

Tengo algunas partículas que identifiqué en una imagen más grande, y necesito analizar imágenes más pequeñas para cada partícula. He utilizado la función 'BoundingBox' de regionprops, pero aún no he tenido éxito. ¿Cómo puedo hacer ahora una subimagen rectangular de la imagen I usando BoundingBox? Puedo usar BoundingBox para dibujar un rectángulo en la imagen original, pero los parámetros devueltos por BoundingBox parecen no ser de dimensión de píxel (x, y, ancho, alto), (x1, y1, x2, y2), etc. esperaría que volviera un cuadro delimitador. He escrito un código de muestra con coins.png para que sea más fácil de entender para cualquiera. ¿Me puede ayudar con esto? ¡Gracias! .

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

Los parámetros devueltos por regionprops son [y,x,width,height] en matriz coordenadas (véase también "unexpected Matlab"

Por lo tanto, para extraer el rectángulo, se escribe:

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))); 

Answer 2

De acuerdo con la documentación de REGIONPROPS:

BoundingBox es [ul_corner width], donde:

• ul_corner: está en forma [x y z ...] y especifica la esquina superior izquierda del cuadro delimitador

• anchura: está en forma [x_width y_width ...] y especifica el ancho del cuadro delimitador a lo largo de cada dimensión

Ahora puede utilizar IMCROP funciona como imcrop(I, rect) donde:

rect es un vector de posición de cuatro elementos [xmin ymin width height] que especifica el tamaño y la posición del rectángulo de recorte.

Por lo tanto:

s = regionprops(L, 'BoundingBox'); 

subImage = imcrop(I, s(1).BoundingBox); 
imshow(subImage)