0-PhHeJdZddlZddlZddlmZddlmZm Z m Z ddl m Z ddl mZgd Zd Zd Zd Zd ZhdZe ddeefddddZe ddeefddddZe dddddZe dddddZe dddddZe dddddZdS)z$ Grayscale morphological operations N)ndimage)_footprint_is_sequencemirror_footprint pad_footprint)default_footprint) DEPRECATED)erosiondilationopeningclosing white_tophat black_tophatc .|d\}}||||||td|D]$}||||||%|ddD]9\}}t|D]$}||||||%:|S)zHelper to call `gray_func` for each footprint in a sequence. `gray_func` is a morphology function that accepts `footprint`, `output`, `mode` and `cval` keyword arguments (e.g. `scipy.ndimage.grey_erosion`). r) footprintoutputmodecvalrN)rangecopy) gray_funcimage footprintsoutrrfpnum_iter_s W/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/skimage/morphology/gray.py_iterate_gray_funcr s a=LB Ier#DtDDDD 1h  NN #((**3TMMMMM"122RR Hx R RA IchhjjBsD Q Q Q Q Q R Jctj|}|jdkr|S|j\}}|dzdkrPtjd|f|j}|rtj||f}ntj||f}|dz }|dzdkrKtj|df|j}|rtj||f}ntj||f}|S)aShift the binary image `footprint` in the left and/or up. This only affects 2D footprints with even number of rows or columns. Parameters ---------- footprint : 2D array, shape (M, N) The input footprint. shift_x, shift_y : bool or None Whether to move `footprint` along each axis. If ``None``, the array is not modified along that dimension. Returns ------- out : 2D array, shape (M + int(shift_x), N + int(shift_y)) The shifted footprint. r rr)npasarrayndimshapezerosdtypevstackhstack)rshift_xshift_ymn extra_row extra_cols r_shift_footprintr1#s& 9%%I~ ?DAq1uzzHaVY_55  : 9i"899II 9i"899I Q1uzzHaVY_55  : 9i"899II 9i"899I r!ctur tur|Sd}tj|tdt |rt fd|DSt |S)zShifts the footprints, whether it's a single array or a sequence. See `_shift_footprint`, which is called for each array in the sequence. zThe parameters `shift_x` and `shift_y` are deprecated since v0.23 and will be removed in v0.26. Use `pad_footprint` or modify the footprintmanually instead.) stacklevelc3DK|]\}}t||fVdS)N)r1).0rr.r+r,s r z$_shift_footprints..[s9XXUR&r7G<r?rAwrapignoremirrornearestreflectr@rLgrrcT|tj|}|tvrtd||dkrd}t |j||\}}t |||}t|d}t|s|dfg}ttj |||||}|S) a Return grayscale morphological erosion of an image. Morphological erosion sets a pixel at (i,j) to the minimum over all pixels in the neighborhood centered at (i,j). Erosion shrinks bright regions and enlarges dark regions. Parameters ---------- image : ndarray Image array. footprint : ndarray or tuple, optional The neighborhood expressed as a 2-D array of 1's and 0's. If None, use a cross-shaped footprint (connectivity=1). The footprint can also be provided as a sequence of smaller footprints as described in the notes below. out : ndarrays, optional The array to store the result of the morphology. If None is passed, a new array will be allocated. mode : str, optional The `mode` parameter determines how the array borders are handled. Valid modes are: 'reflect', 'constant', 'nearest', 'mirror', 'wrap', 'max', 'min', or 'ignore'. If 'max' or 'ignore', pixels outside the image domain are assumed to be the maximum for the image's dtype, which causes them to not influence the result. Default is 'reflect'. cval : scalar, optional Value to fill past edges of input if `mode` is 'constant'. Default is 0.0. .. versionadded:: 0.23 `mode` and `cval` were added in 0.23. Returns ------- eroded : array, same shape as `image` The result of the morphological erosion. Other Parameters ---------------- shift_x, shift_y : DEPRECATED .. deprecated:: 0.23 Notes ----- For ``uint8`` (and ``uint16`` up to a certain bit-depth) data, the lower algorithm complexity makes the :func:`skimage.filters.rank.minimum` function more efficient for larger images and footprints. The footprint can also be a provided as a sequence of 2-tuples where the first element of each 2-tuple is a footprint ndarray and the second element is an integer describing the number of times it should be iterated. For example ``footprint=[(np.ones((9, 1)), 1), (np.ones((1, 9)), 1)]`` would apply a 9x1 footprint followed by a 1x9 footprint resulting in a net effect that is the same as ``footprint=np.ones((9, 9))``, but with lower computational cost. Most of the builtin footprints such as :func:`skimage.morphology.disk` provide an option to automatically generate a footprint sequence of this type. For even-sized footprints, :func:`skimage.morphology.binary_erosion` and this function produce an output that differs: one is shifted by one pixel compared to the other. Examples -------- >>> # Erosion shrinks bright regions >>> import numpy as np >>> from skimage.morphology import footprint_rectangle >>> bright_square = np.array([[0, 0, 0, 0, 0], ... [0, 1, 1, 1, 0], ... [0, 1, 1, 1, 0], ... [0, 1, 1, 1, 0], ... [0, 0, 0, 0, 0]], dtype=np.uint8) >>> erosion(bright_square, footprint_rectangle((3, 3))) array([[0, 0, 0, 0, 0], [0, 0, 0, 0, 0], [0, 0, 1, 0, 0], [0, 0, 0, 0, 0], [0, 0, 0, 0, 0]], dtype=uint8) Nunsupported mode, got rIr?Fpad_endrrrrrrr) r# empty_like_SUPPORTED_MODES ValueErrorrGr(r=rrr ndi grey_erosionrrrr+r,rrs rr r sx {mE"" ###:$::;;; x*5;dCCJD$!)Wg>>Ii777I !) , ,%^$ "       C Jr!cr|tj|}|tvrtd||dkrd}t |j||\}}t |||}t|d}t|}t|s|dfg}ttj |||||}|S) af Return grayscale morphological dilation of an image. Morphological dilation sets the value of a pixel to the maximum over all pixel values within a local neighborhood centered about it. The values where the footprint is 1 define this neighborhood. Dilation enlarges bright regions and shrinks dark regions. Parameters ---------- image : ndarray Image array. footprint : ndarray or tuple, optional The neighborhood expressed as a 2-D array of 1's and 0's. If None, use a cross-shaped footprint (connectivity=1). The footprint can also be provided as a sequence of smaller footprints as described in the notes below. out : ndarray, optional The array to store the result of the morphology. If None is passed, a new array will be allocated. mode : str, optional The `mode` parameter determines how the array borders are handled. Valid modes are: 'reflect', 'constant', 'nearest', 'mirror', 'wrap', 'max', 'min', or 'ignore'. If 'min' or 'ignore', pixels outside the image domain are assumed to be the maximum for the image's dtype, which causes them to not influence the result. Default is 'reflect'. cval : scalar, optional Value to fill past edges of input if `mode` is 'constant'. Default is 0.0. .. versionadded:: 0.23 `mode` and `cval` were added in 0.23. Returns ------- dilated : uint8 array, same shape and type as `image` The result of the morphological dilation. Other Parameters ---------------- shift_x, shift_y : DEPRECATED .. deprecated:: 0.23 Notes ----- For ``uint8`` (and ``uint16`` up to a certain bit-depth) data, the lower algorithm complexity makes the :func:`skimage.filters.rank.maximum` function more efficient for larger images and footprints. The footprint can also be a provided as a sequence of 2-tuples where the first element of each 2-tuple is a footprint ndarray and the second element is an integer describing the number of times it should be iterated. For example ``footprint=[(np.ones((9, 1)), 1), (np.ones((1, 9)), 1)]`` would apply a 9x1 footprint followed by a 1x9 footprint resulting in a net effect that is the same as ``footprint=np.ones((9, 9))``, but with lower computational cost. Most of the builtin footprints such as :func:`skimage.morphology.disk` provide an option to automatically generate a footprint sequence of this type. For non-symmetric footprints, :func:`skimage.morphology.binary_dilation` and :func:`skimage.morphology.dilation` produce an output that differs: `binary_dilation` mirrors the footprint, whereas `dilation` does not. Examples -------- >>> # Dilation enlarges bright regions >>> import numpy as np >>> from skimage.morphology import footprint_rectangle >>> bright_pixel = np.array([[0, 0, 0, 0, 0], ... [0, 0, 0, 0, 0], ... [0, 0, 1, 0, 0], ... [0, 0, 0, 0, 0], ... [0, 0, 0, 0, 0]], dtype=np.uint8) >>> dilation(bright_pixel, footprint_rectangle((3, 3))) array([[0, 0, 0, 0, 0], [0, 1, 1, 1, 0], [0, 1, 1, 1, 0], [0, 1, 1, 1, 0], [0, 0, 0, 0, 0]], dtype=uint8) NrOrIrAFrPrrR) r#rSrTrUrGr(r=rrrr rV grey_dilationrXs rr r sz {mE"" ###:$::;;; x*5;dCCJD$!)Wg>>Ii777I!++I !) , ,%^$ #       C Jr!ct|d}t||||}t|t||||}|S)a| Return grayscale morphological opening of an image. The morphological opening of an image is defined as an erosion followed by a dilation. Opening can remove small bright spots (i.e. "salt") and connect small dark cracks. This tends to "open" up (dark) gaps between (bright) features. Parameters ---------- image : ndarray Image array. footprint : ndarray or tuple, optional The neighborhood expressed as a 2-D array of 1's and 0's. If None, use a cross-shaped footprint (connectivity=1). The footprint can also be provided as a sequence of smaller footprints as described in the notes below. out : ndarray, optional The array to store the result of the morphology. If None is passed, a new array will be allocated. mode : str, optional The `mode` parameter determines how the array borders are handled. Valid modes are: 'reflect', 'constant', 'nearest', 'mirror', 'wrap', 'max', 'min', or 'ignore'. If 'ignore', pixels outside the image domain are assumed to be the maximum for the image's dtype in the erosion, and minimum in the dilation, which causes them to not influence the result. Default is 'reflect'. cval : scalar, optional Value to fill past edges of input if `mode` is 'constant'. Default is 0.0. .. versionadded:: 0.23 `mode` and `cval` were added in 0.23. Returns ------- opening : array, same shape and type as `image` The result of the morphological opening. Notes ----- The footprint can also be a provided as a sequence of 2-tuples where the first element of each 2-tuple is a footprint ndarray and the second element is an integer describing the number of times it should be iterated. For example ``footprint=[(np.ones((9, 1)), 1), (np.ones((1, 9)), 1)]`` would apply a 9x1 footprint followed by a 1x9 footprint resulting in a net effect that is the same as ``footprint=np.ones((9, 9))``, but with lower computational cost. Most of the builtin footprints such as :func:`skimage.morphology.disk` provide an option to automatically generate a footprint sequence of this type. Examples -------- >>> # Open up gap between two bright regions (but also shrink regions) >>> import numpy as np >>> from skimage.morphology import footprint_rectangle >>> bad_connection = np.array([[1, 0, 0, 0, 1], ... [1, 1, 0, 1, 1], ... [1, 1, 1, 1, 1], ... [1, 1, 0, 1, 1], ... [1, 0, 0, 0, 1]], dtype=np.uint8) >>> opening(bad_connection, footprint_rectangle((3, 3))) array([[0, 0, 0, 0, 0], [1, 1, 0, 1, 1], [1, 1, 0, 1, 1], [1, 1, 0, 1, 1], [0, 0, 0, 0, 0]], dtype=uint8) FrPrMrrr)rr r r)rrrrrerodeds rr r nsTNi777I UIDt < < >> # Close a gap between two bright lines >>> import numpy as np >>> from skimage.morphology import footprint_rectangle >>> broken_line = np.array([[0, 0, 0, 0, 0], ... [0, 0, 0, 0, 0], ... [1, 1, 0, 1, 1], ... [0, 0, 0, 0, 0], ... [0, 0, 0, 0, 0]], dtype=np.uint8) >>> closing(broken_line, footprint_rectangle((3, 3))) array([[0, 0, 0, 0, 0], [0, 0, 0, 0, 0], [1, 1, 1, 1, 1], [0, 0, 0, 0, 0], [0, 0, 0, 0, 0]], dtype=uint8) FrPrMr\)rr r r)rrrrrdilateds rrrsTNi777Iuid>>>G '+I66CdQU V V VC Jr!ct||urQt||||}tj|jtrtj|||n||z}|St|||||}tj|jtrtj|||ntj||||S)a Return white top hat of an image. The white top hat of an image is defined as the image minus its morphological opening. This operation returns the bright spots of the image that are smaller than the footprint. Parameters ---------- image : ndarray Image array. footprint : ndarray or tuple, optional The neighborhood expressed as a 2-D array of 1's and 0's. If None, use a cross-shaped footprint (connectivity=1). The footprint can also be provided as a sequence of smaller footprints as described in the notes below. out : ndarray, optional The array to store the result of the morphology. If None is passed, a new array will be allocated. mode : str, optional The `mode` parameter determines how the array borders are handled. Valid modes are: 'reflect', 'constant', 'nearest', 'mirror', 'wrap', 'max', 'min', or 'ignore'. See :func:`skimage.morphology.opening`. Default is 'reflect'. cval : scalar, optional Value to fill past edges of input if `mode` is 'constant'. Default is 0.0. .. versionadded:: 0.23 `mode` and `cval` were added in 0.23. Returns ------- out : array, same shape and type as `image` The result of the morphological white top hat. Notes ----- The footprint can also be a provided as a sequence of 2-tuples where the first element of each 2-tuple is a footprint ndarray and the second element is an integer describing the number of times it should be iterated. For example ``footprint=[(np.ones((9, 1)), 1), (np.ones((1, 9)), 1)]`` would apply a 9x1 footprint followed by a 1x9 footprint resulting in a net effect that is the same as ``footprint=np.ones((9, 9))``, but with lower computational cost. Most of the builtin footprints such as :func:`skimage.morphology.disk` provide an option to automatically generate a footprint sequence of this type. See Also -------- black_tophat References ---------- .. [1] https://en.wikipedia.org/wiki/Top-hat_transform Examples -------- >>> # Subtract gray background from bright peak >>> import numpy as np >>> from skimage.morphology import footprint_rectangle >>> bright_on_gray = np.array([[2, 3, 3, 3, 2], ... [3, 4, 5, 4, 3], ... [3, 5, 9, 5, 3], ... [3, 4, 5, 4, 3], ... [2, 3, 3, 3, 2]], dtype=np.uint8) >>> white_tophat(bright_on_gray, footprint_rectangle((3, 3))) array([[0, 0, 0, 0, 0], [0, 0, 1, 0, 0], [0, 1, 5, 1, 0], [0, 0, 1, 0, 0], [0, 0, 0, 0, 0]], dtype=uint8) rMrr\)r r#rBr(rC logical_xorsubtract)rrrrropeneds rrrsV e|| 4@@@ =t , ,  N3C 0 0 0 0 0 6MC  %$T B B BC }SY%%) ucs+++++ E3C(((( Jr!c~||urct||||}tj|jtrtj|||ntj||||St|||||}tj|jtjrtj|||n||z}|S)aY Return black top hat of an image. The black top hat of an image is defined as its morphological closing minus the original image. This operation returns the dark spots of the image that are smaller than the footprint. Note that dark spots in the original image are bright spots after the black top hat. Parameters ---------- image : ndarray Image array. footprint : ndarray or tuple, optional The neighborhood expressed as a 2-D array of 1's and 0's. If None, use a cross-shaped footprint (connectivity=1). The footprint can also be provided as a sequence of smaller footprints as described in the notes below. out : ndarray, optional The array to store the result of the morphology. If None is passed, a new array will be allocated. mode : str, optional The `mode` parameter determines how the array borders are handled. Valid modes are: 'reflect', 'constant', 'nearest', 'mirror', 'wrap', 'max', 'min', or 'ignore'. See :func:`skimage.morphology.closing`. Default is 'reflect'. cval : scalar, optional Value to fill past edges of input if `mode` is 'constant'. Default is 0.0. .. versionadded:: 0.23 `mode` and `cval` were added in 0.23. Returns ------- out : array, same shape and type as `image` The result of the morphological black top hat. Notes ----- The footprint can also be a provided as a sequence of 2-tuples where the first element of each 2-tuple is a footprint ndarray and the second element is an integer describing the number of times it should be iterated. For example ``footprint=[(np.ones((9, 1)), 1), (np.ones((1, 9)), 1)]`` would apply a 9x1 footprint followed by a 1x9 footprint resulting in a net effect that is the same as ``footprint=np.ones((9, 9))``, but with lower computational cost. Most of the builtin footprints such as :func:`skimage.morphology.disk` provide an option to automatically generate a footprint sequence of this type. See Also -------- white_tophat References ---------- .. [1] https://en.wikipedia.org/wiki/Top-hat_transform Examples -------- >>> # Change dark peak to bright peak and subtract background >>> import numpy as np >>> from skimage.morphology import footprint_rectangle >>> dark_on_gray = np.array([[7, 6, 6, 6, 7], ... [6, 5, 4, 5, 6], ... [6, 4, 0, 4, 6], ... [6, 5, 4, 5, 6], ... [7, 6, 6, 6, 7]], dtype=np.uint8) >>> black_tophat(dark_on_gray, footprint_rectangle((3, 3))) array([[0, 0, 0, 0, 0], [0, 0, 1, 0, 0], [0, 1, 5, 1, 0], [0, 0, 1, 0, 0], [0, 0, 0, 0, 0]], dtype=uint8) rMrar\)rr#rBr(rCrbrcbool_)rrrrrcloseds rrresX e|| 4@@@ =t , , . N63C 0 0 0 0 0 K - - - - %$T B B BC }SY)) sEs+++++ u  Jr!)NN)__doc__r8numpyr#scipyrrVrrrrmiscr _shared.utilsr __all__r r1r=rGrTr r r rrrr!rrosJ OOOOOOOOOO######&&&&&& X W W"%%%P999(*      q  qqqqqh   u  uuuuupIYSIIIIIXIYSIIIIIXY)#YYYYYxY)#YYYYYYYr!