def upperGeodesicDilate3D(imIn, imMask, imOut, n=1, se=m3D.CUBOCTAHEDRON1):
"""
Performs a upper geodesic dilation of 3D image 'imIn' above 'imMask'.
The result is put inside 'imOut', 'n' controls the size of the dilation.
'se' specifies the type of structuring element used to perform the
computation (CUBOCTAHEDRON by default).
Warning! 'imMask' and 'imOut' must be different.
"""
m3D.logic3D(imIn, imMask, imOut, "sup")
if imIn.getDepth() == 1:
for i in range(n):
m3D.diff3D(imOut, imMask, imOut)
m3D.dilate3D(imOut, imOut, se=se)
m3D.logic3D(imMask, imOut, imOut, "sup")
else:
imWrk1 = m3D.image3DMb(imIn)
imWrk2 = m3D.image3DMb(imIn, 1)
for i in range(n):
m3D.generateSupMask3D(imOut, imMask, imWrk2, True)
m3D.convertByMask3D(imWrk2, imWrk1, 0, m3D.computeMaxRange3D(imWrk1)[1])
m3D.logic3D(imOut, imWrk1, imOut, "inf")
m3D.dilate3D(imOut, imOut, se=se)
m3D.logic3D(imOut, imMask, imOut, "sup")
def lowerGeodesicErode3D(imIn, imMask, imOut, n=1, se=m3D.CUBOCTAHEDRON1):
"""
Performs a lower geodesic erosion of 3D image 'imIn' under 'imMask'.
The result is put inside 'imOut', 'n' controls the size of the erosion.
'se' specifies the type of structuring element used to perform the
computation (CUBOCTAHEDRON by default).
The binary lower geodesic erosion is realised using the fact that the
dilation is the dual operation of the erosion.
Warning! 'imMask' and 'imOut' must be different.
"""
if imIn.getDepth() == 1:
m3D.diff3D(imMask, imIn, imOut)
lowerGeodesicDilate3D(imOut, imMask, imOut, n, se=se)
m3D.diff3D(imMask, imOut, imOut)
else:
imWrk1 = m3D.image3DMb(imIn)
imWrk2 = m3D.image3DMb(imIn, 1)
m3D.logic3D(imIn, imMask, imOut, "inf")
for i in range(n):
m3D.generateSupMask3D(imOut, imMask, imWrk2, False)
m3D.convertByMask3D(imWrk2, imWrk1, 0, m3D.computeMaxRange3D(imWrk1)[1])
m3D.logic3D(imOut, imWrk1, imOut, "sup")
m3D.erode3D(imOut, imOut, se=se)
m3D.logic3D(imOut, imMask, imOut, "inf")
def upperGeodesicDilate3D(imIn, imMask, imOut, n=1, se=m3D.CUBOCTAHEDRON1):
"""
Performs a upper geodesic dilation of 3D image 'imIn' above 'imMask'.
The result is put inside 'imOut', 'n' controls the size of the dilation.
'se' specifies the type of structuring element used to perform the
computation (CUBOCTAHEDRON by default).
Warning! 'imMask' and 'imOut' must be different.
"""
m3D.logic3D(imIn, imMask, imOut, "sup")
if imIn.getDepth() == 1:
for i in range(n):
m3D.diff3D(imOut, imMask, imOut)
m3D.dilate3D(imOut, imOut, se=se)
m3D.logic3D(imMask, imOut, imOut, "sup")
else:
imWrk1 = m3D.image3DMb(imIn)
imWrk2 = m3D.image3DMb(imIn, 1)
for i in range(n):
m3D.generateSupMask3D(imOut, imMask, imWrk2, True)
m3D.convertByMask3D(imWrk2, imWrk1, 0,
m3D.computeMaxRange3D(imWrk1)[1])
m3D.logic3D(imOut, imWrk1, imOut, "inf")
m3D.dilate3D(imOut, imOut, se=se)
m3D.logic3D(imOut, imMask, imOut, "sup")
def lowerGeodesicErode3D(imIn, imMask, imOut, n=1, se=m3D.CUBOCTAHEDRON1):
"""
Performs a lower geodesic erosion of 3D image 'imIn' under 'imMask'.
The result is put inside 'imOut', 'n' controls the size of the erosion.
'se' specifies the type of structuring element used to perform the
computation (CUBOCTAHEDRON by default).
The binary lower geodesic erosion is realised using the fact that the
dilation is the dual operation of the erosion.
Warning! 'imMask' and 'imOut' must be different.
"""
if imIn.getDepth() == 1:
m3D.diff3D(imMask, imIn, imOut)
lowerGeodesicDilate3D(imOut, imMask, imOut, n, se=se)
m3D.diff3D(imMask, imOut, imOut)
else:
imWrk1 = m3D.image3DMb(imIn)
imWrk2 = m3D.image3DMb(imIn, 1)
m3D.logic3D(imIn, imMask, imOut, "inf")
for i in range(n):
m3D.generateSupMask3D(imOut, imMask, imWrk2, False)
m3D.convertByMask3D(imWrk2, imWrk1, 0,
m3D.computeMaxRange3D(imWrk1)[1])
m3D.logic3D(imOut, imWrk1, imOut, "sup")
m3D.erode3D(imOut, imOut, se=se)
m3D.logic3D(imOut, imMask, imOut, "inf")
def simpleLevelling3D(imIn, imMask, imOut, grid=m3D.DEFAULT_GRID3D):
"""
Performs a simple levelling of 3D image 'imIn' controlled by image 'imMask'
and puts the result in 'imOut'. This operation is composed of two
geodesic reconstructions. This filter tends to level regions in the
image of homogeneous grey values.
"""
imWrk1 = m3D.image3DMb(imIn)
imWrk2 = m3D.image3DMb(imIn)
mask_im = m3D.image3DMb(imIn, 1)
m3D.logic3D(imIn, imMask, imWrk1, "inf")
m3D.build3D(imIn, imWrk1, grid=grid)
m3D.logic3D(imIn, imMask, imWrk2, "sup")
m3D.dualBuild3D(imIn, imWrk2, grid=grid)
m3D.generateSupMask3D(imIn, imMask, mask_im, False)
m3D.convertByMask3D(mask_im, imOut, 0, m3D.computeMaxRange3D(imIn)[1])
m3D.logic3D(imOut, imWrk1, imWrk1, "inf")
m3D.negate3D(imOut, imOut)
m3D.logic3D(imOut, imWrk2, imOut, "inf")
m3D.logic3D(imWrk1, imOut, imOut, "sup")
def simpleLevelling3D(imIn, imMask, imOut, grid=m3D.DEFAULT_GRID3D):
"""
Performs a simple levelling of 3D image 'imIn' controlled by image 'imMask'
and puts the result in 'imOut'. This operation is composed of two
geodesic reconstructions. This filter tends to level regions in the
image of homogeneous grey values.
"""
imWrk1 = m3D.image3DMb(imIn)
imWrk2 = m3D.image3DMb(imIn)
mask_im = m3D.image3DMb(imIn, 1)
m3D.logic3D(imIn, imMask, imWrk1, "inf")
m3D.build3D(imIn, imWrk1, grid=grid)
m3D.logic3D(imIn, imMask, imWrk2, "sup")
m3D.dualBuild3D(imIn, imWrk2, grid=grid)
m3D.generateSupMask3D(imIn, imMask, mask_im, False)
m3D.convertByMask3D(mask_im, imOut, 0, m3D.computeMaxRange3D(imIn)[1])
m3D.logic3D(imOut, imWrk1, imWrk1, "inf")
m3D.negate3D(imOut, imOut)
m3D.logic3D(imOut, imWrk2, imOut, "inf")
m3D.logic3D(imWrk1, imOut, imOut, "sup")
