ImageIO find intensity

 

# Create left ventricle mask

labels, nlabels = ndi.label(mask)

lv_val = labels[128, 128]

lv_mask = np.where(labels == lv_val, 1, 0)

# Find bounding box of left ventricle

bboxes = ndi.find_objects(lv_mask)

print('Number of objects:', len(bboxes))

print('Indices for first box:', bboxes[0])

# Crop to the left ventricle (index 0)

im_lv = im[bboxes[0]]

# Plot the cropped image

plt.imshow(im_lv, cmap='rainbow')

format_and_render_plot()