Mahotas – Getting Border of label at given region
Last Updated :
24 Feb, 2023
In this article we will see how we can get the border of the label in the labelled image within the given region in mahotas. For this we are going to use the fluorescent microscopy image from a nuclear segmentation benchmark. We can get the image with the help of command given below
mahotas.demos.nuclear_image()
Below is the nuclear_image
In order to do this we will use mahotas.labelled.border method
Syntax : mahotas.labelled.border(labelled_image, i, j)
Argument : It takes numpy.ndarray object i.e labelled image and two integer as argument
Note : A pixel is on the border if it has value i (or j)
Return : It returns numpy.ndarray object i.e labelled image with border label
Note : The input of the this should be the filtered image object which is labelled
In order to filter the image we will take the image object which is numpy.ndarray and filter it with the help of indexing, below is the command to do this
image = image[:, :, 0]
Example 1 :
Python3
import mahotas
import numpy as np
from pylab import imshow, show
import os
f = mahotas.demos.load( 'nuclear' )
f = f[:, :, 0 ]
f = mahotas.gaussian_filter(f, 4 )
f = (f> f.mean())
labelled, n_nucleus = mahotas.label(f)
print ( "Labelled Image" )
imshow(labelled)
show()
relabeled = mahotas.labelled.border(labelled, 0 , 20 )
print ( "Labels With borders at given point" )
imshow(relabelled)
show()
|
Output :
Example 2 :
Python3
import numpy as np
import mahotas
from pylab import imshow, show
img = mahotas.imread( 'dog_image.png' )
img = img[:, :, 0 ]
gaussian = mahotas.gaussian_filter(img, 15 )
gaussian = (gaussian > gaussian.mean())
labelled, n_nucleus = mahotas.label(gaussian)
print ( "Labelled Image" )
imshow(labelled)
show()
relabeled = mahotas.labelled.border(labelled, 1 , 0 )
print ( "Labels With borders at given point" )
imshow(relabelled)
show()
|
Output :
Share your thoughts in the comments
Please Login to comment...