def conjugateDirectionalErode(imIn, imOut, d, size, grid=mamba.DEFAULT_GRID, edge=mamba.FILLED):
"""
Performs the erosion of image 'imIn' in the conjugate direction 'd' of the
grid and puts the result in image 'imOut'. The images can be binary, greyscale
or 32-bit images. 'size' is a multiple of the distance between two adjacent
points in the conjugate directions. 'edge' is set to FILLED by default (can
be changed).
Note that this operator is not equivalent to successive erosions by a doublet
of points.
"""
imWrk1 = mamba.imageMb(imIn)
imWrk2 = mamba.imageMb(imIn)
mamba.copy(imIn, imOut)
j3 = mamba.transposeDirection(d, grid=grid)
j2 = mamba.rotateDirection(d, grid=grid)
j1 = mamba.transposeDirection(j2, grid=grid)
val = mamba.computeMaxRange(imIn)[1] * int(edge == mamba.FILLED)
for i in range(size):
mamba.copy(imOut, imWrk1)
mamba.copy(imOut, imWrk2)
mamba.linearErode(imWrk1, imWrk1, d, n=1, grid=grid, edge=edge)
mamba.shift(imWrk1, imWrk1, j1, 1, val, grid=grid)
mamba.logic(imWrk1, imOut, imWrk1, "inf")
mamba.linearErode(imWrk2, imWrk2, j2, n=1, grid=grid, edge=edge)
mamba.shift(imWrk2, imWrk2, j3, 1, val, grid=grid)
mamba.logic(imWrk2, imOut, imWrk2, "inf")
mamba.logic(imWrk1, imWrk2, imOut, "sup")
def conjugateDirectionalDilate(imIn, imOut, d, size, grid=mamba.DEFAULT_GRID, edge=mamba.EMPTY):
"""
Performs the dilation of image 'imIn' in the conjugate direction 'd' of the
grid and puts the result in image 'imOut'. The images can be binary, greyscale
or 32-bit images. 'size' is a multiple of the distance between two adjacent
points in the conjugate directions. 'edge' is set to EMPTY by default (can
be changed).
Note that the linear structuring element is not connected. Points connecting
adjacent points in the conjugate directions are not present. This is normal
if we want to insure that the result is identical when we iterate n size 1
operations to get a size n one (the same remark applies to the erosion).
"""
imWrk1 = mamba.imageMb(imIn)
imWrk2 = mamba.imageMb(imIn)
mamba.copy(imIn, imOut)
j3 = mamba.transposeDirection(d, grid=grid)
j2 = mamba.rotateDirection(d, grid=grid)
j1 = mamba.transposeDirection(j2, grid=grid)
val = mamba.computeMaxRange(imIn)[1] * int(edge == mamba.FILLED)
for i in range(size):
mamba.copy(imOut, imWrk1)
mamba.copy(imOut, imWrk2)
mamba.linearDilate(imWrk1, imWrk1, d, 1, grid=grid, edge=edge)
mamba.shift(imWrk1, imWrk1, j1, 1, val, grid=grid)
mamba.logic(imWrk1, imOut, imWrk1, "sup")
mamba.linearDilate(imWrk2, imWrk2, j2, 1, grid=grid, edge=edge)
mamba.shift(imWrk2, imWrk2, j3, 1, val, grid=grid)
mamba.logic(imWrk2, imOut, imWrk2, "sup")
mamba.logic(imWrk1, imWrk2, imOut, "inf")
def conjugateDirectionalErode(imIn,
imOut,
d,
size,
grid=mamba.DEFAULT_GRID,
edge=mamba.FILLED):
"""
Performs the erosion of image 'imIn' in the conjugate direction 'd' of the
grid and puts the result in image 'imOut'. The images can be binary, greyscale
or 32-bit images. 'size' is a multiple of the distance between two adjacent
points in the conjugate directions. 'edge' is set to FILLED by default (can
be changed).
Note that this operator is not equivalent to successive erosions by a doublet
of points.
"""
imWrk1 = mamba.imageMb(imIn)
imWrk2 = mamba.imageMb(imIn)
mamba.copy(imIn, imOut)
j3 = mamba.transposeDirection(d, grid=grid)
j2 = mamba.rotateDirection(d, grid=grid)
j1 = mamba.transposeDirection(j2, grid=grid)
val = mamba.computeMaxRange(imIn)[1] * int(edge == mamba.FILLED)
for i in range(size):
mamba.copy(imOut, imWrk1)
mamba.copy(imOut, imWrk2)
mamba.linearErode(imWrk1, imWrk1, d, n=1, grid=grid, edge=edge)
mamba.shift(imWrk1, imWrk1, j1, 1, val, grid=grid)
mamba.logic(imWrk1, imOut, imWrk1, "inf")
mamba.linearErode(imWrk2, imWrk2, j2, n=1, grid=grid, edge=edge)
mamba.shift(imWrk2, imWrk2, j3, 1, val, grid=grid)
mamba.logic(imWrk2, imOut, imWrk2, "inf")
mamba.logic(imWrk1, imWrk2, imOut, "sup")
def conjugateDirectionalDilate(imIn, imOut, d, size, grid=mamba.DEFAULT_GRID, edge=mamba.EMPTY):
"""
Performs the dilation of image 'imIn' in the conjugate direction 'd' of the
grid and puts the result in image 'imOut'. The images can be binary, greyscale
or 32-bit images. 'size' is a multiple of the distance between two adjacent
points in the conjugate directions. 'edge' is set to EMPTY by default (can
be changed).
Note that the linear structuring element is not connected. Points connecting
adjacent points in the conjugate directions are not present. This is normal
if we want to insure that the result is identical when we iterate n size 1
operations to get a size n one (the same remark applies to the erosion).
"""
imWrk1 = mamba.imageMb(imIn)
imWrk2 = mamba.imageMb(imIn)
mamba.copy(imIn, imOut)
j3 = mamba.transposeDirection(d, grid=grid)
j2 = mamba.rotateDirection(d, grid=grid)
j1 = mamba.transposeDirection(j2, grid=grid)
val = mamba.computeMaxRange(imIn)[1] * int(edge == mamba.FILLED)
for i in range(size):
mamba.copy(imOut, imWrk1)
mamba.copy(imOut, imWrk2)
mamba.linearDilate(imWrk1, imWrk1, d, 1, grid=grid, edge=edge)
mamba.shift(imWrk1, imWrk1, j1, 1, val, grid=grid)
mamba.logic(imWrk1, imOut, imWrk1, "sup")
mamba.linearDilate(imWrk2, imWrk2, j2, 1, grid=grid, edge=edge)
mamba.shift(imWrk2, imWrk2, j3, 1, val, grid=grid)
mamba.logic(imWrk2, imOut, imWrk2, "sup")
mamba.logic(imWrk1, imWrk2, imOut, "inf")
def linearOpen(imIn, imOut, dir, n, grid=mamba.DEFAULT_GRID, edge=mamba.FILLED):
"""
Performs an opening by a segment of size 'n' in direction 'dir'.
'edge' is set to 'FILLED' by default.
"""
mamba.linearErode(imIn, imOut, dir, n, edge=edge, grid=grid)
mamba.linearDilate(imOut, imOut, mamba.transposeDirection(dir, grid=grid), n, grid=grid)
def linearClose(imIn, imOut, dir, n, grid=mamba.DEFAULT_GRID, edge=mamba.FILLED):
"""
Performs a closing by a segment of size 'n' in direction 'dir'.
If 'edge' is set to 'EMPTY', the operation must be modified to remain extensive.
"""
imWrk = mamba.imageMb(imIn)
if edge==mamba.EMPTY:
mamba.copy(imIn, imWrk)
mamba.linearDilate(imIn, imOut, dir, n, grid=grid)
mamba.linearErode(imOut, imOut, mamba.transposeDirection(dir, grid=grid), n, edge=edge, grid=grid)
if edge==mamba.EMPTY:
mamba.logic(imOut, imWrk, imOut, "sup")
def linearOpen(imIn,
imOut,
dir,
n,
grid=mamba.DEFAULT_GRID,
edge=mamba.FILLED):
"""
Performs an opening by a segment of size 'n' in direction 'dir'.
'edge' is set to 'FILLED' by default.
"""
mamba.linearErode(imIn, imOut, dir, n, edge=edge, grid=grid)
mamba.linearDilate(imOut,
imOut,
mamba.transposeDirection(dir, grid=grid),
n,
grid=grid)
def infFarNeighbor3D(imIn, imInOut, nb, amp, grid=m3D.DEFAULT_GRID3D, edge=mamba.FILLED):
"""
Performs an infimum operation between the 'imInOut' 3D image pixels and their neighbor 'nb'
at distance 'amp' according to 'grid' in 3D image 'imIn'. The result is put in 'imInOut'.
"grid' value can be CUBIC, CENTER_CUBIC or FACE_CENTER_CUBIC. 'edge' value can be EMPTY
or FILLED.
"""
(width,height,length) = imIn.getSize()
if length!=len(imInOut):
mamba.raiseExceptionOnError(core.MB_ERR_BAD_SIZE)
imWrk = mamba.imageMb(imIn[0])
# Computing limits according to the scanning direction.
scan = grid.convertFromDir(nb,0)[0]
if scan == 0:
startPlane, endPlane, scanDir = 0, length, 1
startFill, endFill = 0, 0
elif scan == 1:
startPlane, endPlane, scanDir = 0, length - amp, 1
startFill, endFill = max(length - amp, 0), length
else:
startPlane, endPlane, scanDir = length - 1, amp - 1, -1
startFill, endFill = 0, min(amp, length)
# Performing the shift operations given by the getShiftDirList method.
if edge == mamba.EMPTY:
fillValue = 0
else:
fillValue = mamba.computeMaxRange(imIn[0])[1]
for i in range(startPlane, endPlane, scanDir):
j = i + amp * scan
dirList = grid.getShiftDirsList(nb, amp, i)
if len(dirList) == 1:
mamba.infFarNeighbor(imIn[j], imInOut[i], dirList[0][0] , dirList[0][1], grid=dirList[0][2], edge=edge)
else:
d = mamba.transposeDirection(dirList[1][0], dirList[1][2])
mamba.shift(imIn[j], imWrk, d, dirList[1][1], fillValue, grid=dirList[1][2])
mamba.infFarNeighbor(imWrk, imInOut[i], dirList[0][0] , dirList[0][1], grid=dirList[0][2], edge=edge)
# Filling the necessary planes.
if edge == mamba.EMPTY:
for i in range(startFill, endFill):
imInOut[i].fill(fillValue)
def linearClose(imIn,
imOut,
dir,
n,
grid=mamba.DEFAULT_GRID,
edge=mamba.FILLED):
"""
Performs a closing by a segment of size 'n' in direction 'dir'.
If 'edge' is set to 'EMPTY', the operation must be modified to remain extensive.
"""
imWrk = mamba.imageMb(imIn)
if edge == mamba.EMPTY:
mamba.copy(imIn, imWrk)
mamba.linearDilate(imIn, imOut, dir, n, grid=grid)
mamba.linearErode(imOut,
imOut,
mamba.transposeDirection(dir, grid=grid),
n,
edge=edge,
grid=grid)
if edge == mamba.EMPTY:
mamba.logic(imOut, imWrk, imOut, "sup")
