9i)ASrSSKrSSKrSSKJr SSKJr SSK J r SSK J r S r \ Vs0sHoS U-_M snrS r \R\ Vs0sHoU_M sn5 S rS rSSS.SjjrSSS.SjjrSSSSS.Sjrgs snfs snf)z#Miscellaneous morphology functions.N)ndimage)cKDTree)warn)_remove_objects_by_distance)erosiondilationopeningclosinggrey_)binary_erosionbinary_dilationbinary_openingbinary_closing black_tophat white_tophatcJ^[R"T5SU4Sjj5nU$)a_Decorator to add a default footprint to morphology functions. Parameters ---------- func : function A morphology function such as erosion, dilation, opening, closing, white_tophat, or black_tophat. Returns ------- func_out : function The function, using a default footprint of same dimension as the input image with connectivity 1. ch>Uc![R"URS5nT"U/UQ7SU0UD6$)Nr footprint)ndigenerate_binary_structurendim)imagerargskwargsfuncs W/var/www/html/land-doc-ocr/venv/lib/python3.13/site-packages/skimage/morphology/misc.pyfunc_out#default_footprint..func_out.s9  55ejj!DIE@@Y@@@N) functoolswraps)rrs` rdefault_footprintr%s)"__TAA Or!cUR[:XdN[R"UR[R5(d[ SURS35egg)Nz4Only bool or integer image types are supported. Got .)dtypeboolnp issubdtypeinteger TypeError)ars r_check_dtype_supportedr/7sL HH  bhh C C EbhhZq Q  !D r!outc,[U5 UcUR5nOXSS&US:XaU$UR[:XaZ[R "UR U5n[R"U[RS9n[R"XUS9 OUn[R"UR55n[U5S:XaUR[:wa [S5 Xa:nXunSX8'U$![a [S5ef=f)aRemove objects smaller than the specified size. Expects ar to be an array with labeled objects, and removes objects smaller than min_size. If `ar` is bool, the image is first labeled. This leads to potentially different behavior for bool and 0-and-1 arrays. Parameters ---------- ar : ndarray (arbitrary shape, int or bool type) The array containing the objects of interest. If the array type is int, the ints must be non-negative. min_size : int, optional (default: 64) The smallest allowable object size. connectivity : int, {1, 2, ..., ar.ndim}, optional (default: 1) The connectivity defining the neighborhood of a pixel. Used during labelling if `ar` is bool. out : ndarray Array of the same shape as `ar`, into which the output is placed. By default, a new array is created. Raises ------ TypeError If the input array is of an invalid type, such as float or string. ValueError If the input array contains negative values. Returns ------- out : ndarray, same shape and type as input `ar` The input array with small connected components removed. See Also -------- skimage.morphology.remove_objects_by_distance Examples -------- >>> from skimage import morphology >>> a = np.array([[0, 0, 0, 1, 0], ... [1, 1, 1, 0, 0], ... [1, 1, 1, 0, 1]], bool) >>> b = morphology.remove_small_objects(a, 6) >>> b array([[False, False, False, False, False], [ True, True, True, False, False], [ True, True, True, False, False]]) >>> c = morphology.remove_small_objects(a, 7, connectivity=2) >>> c array([[False, False, False, True, False], [ True, True, True, False, False], [ True, True, True, False, False]]) >>> d = morphology.remove_small_objects(a, 6, out=a) >>> d is a True Nrr()outputz{Negative value labels are not supported. Try relabeling the input with `scipy.ndimage.label` or `skimage.morphology.label`.rz[Only one label was provided to `remove_small_objects`. Did you mean to use a boolean array?)r/copyr(r)rrrr* zeros_likeint32labelbincountravel ValueErrorlenr) r.min_size connectivityr1rccscomponent_sizes too_smalltoo_small_masks rremove_small_objectsrC?sx2 {ggiA1}  yyD11"''<H mmBbhh/ ", ++ciik2 ?q SYY$%6  3  *I^NC J#   *   s $C==Dc6[U5 UR[:wa[S[5 Ub UR[:wa [ S5eOUR [SS9n[R"XS9 [X1X#S9n[R"X3S9 U$)aRemove contiguous holes smaller than the specified size. Parameters ---------- ar : ndarray (arbitrary shape, int or bool type) The array containing the connected components of interest. area_threshold : int, optional (default: 64) The maximum area, in pixels, of a contiguous hole that will be filled. Replaces `min_size`. connectivity : int, {1, 2, ..., ar.ndim}, optional (default: 1) The connectivity defining the neighborhood of a pixel. out : ndarray Array of the same shape as `ar` and bool dtype, into which the output is placed. By default, a new array is created. Raises ------ TypeError If the input array is of an invalid type, such as float or string. ValueError If the input array contains negative values. Returns ------- out : ndarray, same shape and type as input `ar` The input array with small holes within connected components removed. Examples -------- >>> from skimage import morphology >>> a = np.array([[1, 1, 1, 1, 1, 0], ... [1, 1, 1, 0, 1, 0], ... [1, 0, 0, 1, 1, 0], ... [1, 1, 1, 1, 1, 0]], bool) >>> b = morphology.remove_small_holes(a, 2) >>> b array([[ True, True, True, True, True, False], [ True, True, True, True, True, False], [ True, False, False, True, True, False], [ True, True, True, True, True, False]]) >>> c = morphology.remove_small_holes(a, 2, connectivity=2) >>> c array([[ True, True, True, True, True, False], [ True, True, True, False, True, False], [ True, False, False, True, True, False], [ True, True, True, True, True, False]]) >>> d = morphology.remove_small_holes(a, 2, out=a) >>> d is a True Notes ----- If the array type is int, it is assumed that it contains already-labeled objects. The labels are not kept in the output image (this function always outputs a bool image). It is suggested that labeling is completed after using this function. z\Any labeled images will be returned as a boolean array. Did you mean to use a boolean array?zout dtype must be boolTr5r0) r/r(r)r UserWarningr-astyper* logical_notrC)r.area_thresholdr>r1s rremove_small_holesrJsv2 xx4  3    99 45 5 ii4i(NN2 sL JCNN3 Jr!)priorityp_normspacingr1c ^^US:a[SU35e[R"UR[R5(d[SUR35eUcUR SS9nOXPLaXSS&UR SS9nTbP[R"T5mTRUR4:wdTR5S::a [S5e[R"U5nTb![R"TST:H5(a(US::a"[R"URS 5nO+[R"URUR5n[R"XXS 9[R "XXS 9:gR 5Un Xyn Xy)n U R"S:XaU$U [R$"Xj5n U [R$"Xk5n Uc|[R&"UR[R(S S 9(d3[R*"UR-[R(S S95nO[R*"U5nXjSn U S:a[SU <35eU [R$"X&U SS9SSS2n [R2"XR5mTb)[5UU4Sj[7UR555m[9[R:"T[R<S9R>S9n[AUU U UUUURS9 URBULaUREUR5USS&U$![.aMn [R0"U5S -4nURU:wa[SUSURS35U eeSn A ff=f)aRemove objects, in specified order, until remaining are a minimum distance apart. Remove labeled objects from an image until the remaining ones are spaced more than a given distance from one another. By default, smaller objects are removed first. Parameters ---------- label_image : ndarray of integers An n-dimensional array containing object labels, e.g. as returned by :func:`~.label`. A value of zero is considered background, all other object IDs must be positive integers. min_distance : int or float Remove objects whose distance to other objects is not greater than this positive value. Objects with a lower `priority` are removed first. priority : ndarray, optional Defines the priority with which objects are removed. Expects a 1-dimensional array of length :func:`np.amax(label_image) + 1 ` that contains the priority for each object's label at the respective index. Objects with a lower value are removed first until all remaining objects fulfill the distance requirement. If not given, priority is given to objects with a higher number of samples and their label value second. p_norm : int or float, optional The Minkowski distance of order p, used to calculate the distance between objects. The default ``2`` corresponds to the Euclidean distance, ``1`` to the "Manhattan" distance, and ``np.inf`` to the Chebyshev distance. spacing : sequence of float, optional The pixel spacing along each axis of `label_image`. If not specified, a grid spacing of unity (1) is implied. out : ndarray, optional Array of the same shape and dtype as `image`, into which the output is placed. By default, a new array is created. Returns ------- out : ndarray Array of the same shape as `label_image`, for which objects that violate the `min_distance` condition were removed. See Also -------- skimage.morphology.remove_small_objects Remove objects smaller than the specified size. Notes ----- The basic steps of this algorithm work as follows: 1. Find the indices for of all given objects and separate them depending on if they point to an object's border or not. 2. Sort indices by their label value, ensuring that indices which point to the same object are next to each other. This optimization allows finding all parts of an object, simply by stepping to the neighboring indices. 3. Sort boundary indices by `priority`. Use a stable-sort to preserve the ordering from the previous sorting step. If `priority` is not given, use :func:`numpy.bincount` as a fallback. 4. Construct a :class:`scipy.spatial.cKDTree` from the boundary indices. 5. Iterate across boundary indices in priority-sorted order, and query the kd-tree for objects that are too close. Remove ones that are and don't take them into account when evaluating other objects later on. The performance of this algorithm depends on the number of samples in `label_image` that belong to an object's border. Examples -------- >>> import skimage as ski >>> ski.morphology.remove_objects_by_distance(np.array([2, 0, 1, 1]), 2) array([0, 0, 1, 1]) >>> ski.morphology.remove_objects_by_distance( ... np.array([2, 0, 1, 1]), 2, priority=np.array([0, 1, 9]) ... ) array([2, 0, 0, 0]) >>> label_image = np.array( ... [[8, 0, 0, 0, 0, 0, 0, 0, 0, 9, 9], ... [8, 8, 8, 0, 0, 0, 0, 0, 0, 9, 9], ... [0, 0, 0, 0, 0, 0, 0, 0, 9, 0, 0], ... [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], ... [0, 0, 3, 0, 0, 0, 1, 0, 0, 0, 0], ... [2, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0], ... [0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0], ... [0, 0, 0, 0, 0, 0, 0, 0, 0, 7, 7]] ... ) >>> ski.morphology.remove_objects_by_distance( ... label_image, min_distance=3 ... ) array([[8, 0, 0, 0, 0, 0, 0, 0, 0, 9, 9], [8, 8, 8, 0, 0, 0, 0, 0, 0, 9, 9], [0, 0, 0, 0, 0, 0, 0, 0, 9, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 7, 7]]) rzmin_distance must be >= 0, was z,`label_image` must be of integer dtype, got NC)orderzV`spacing` must contain exactly one positive factor for each dimension of `label_image`rr)rsafe)castingFrEzfound object with negative ID stable)kindzBshape of `priority` must be (np.amax(label_image) + 1,), expected z, got z insteadc3:># UHnTUTU-v M g7fr").0dimrMunraveled_indicess r -remove_objects_by_distance..s#" =Lc c "WS\ 1_sr3)data)r1border_indices inner_indiceskdtree min_distancerLshape)#r;r*r+r(r,r5r:arrayrbrmin flatnonzeroallrrmaximum_filterminimum_filtersizeargsortcan_castintpr9rG IndexErroramax unravel_indextuplerangerasarrayfloat64Trbasereshape) label_imagerarKrLrMr1 out_raveledindicesrborderr^r_ smallest_iderrorexpected_shaper`rZs ` @rremove_objects_by_distancer~sTa:<.IJJ ==**BJJ 7 7:;;L;L:M N   {S) A))#)&K((7# ==SXXK '7;;=A+=6  nn[)G 266'!*"788fk11#((A> 11#((CHHE  34   c 7 8 eggF_NG$Ma $BJJ{/J$KLN!"**[-G"HIM{{399bggv>{{;#5#5bggE#5#JKH{{;/HQ/0KQ9+IJJ( JJxN ;<8 LTrT R ((C!" =B388_"  "**%6bjjIKK LF %#!s"$$SYY/A JA  ''+.24 >>^ +*+6(..1AK   s#M?? O AOO)@r)__doc__numpyr*r#scipyrr scipy.spatialr _shared.utilsr_misc_cyrfuncsskimage2ndimageupdater%r/rCrJr~)xs0rrs) ! 1 6+015agk>51 e,e1e,-4 ``FSTSt  NS2-s B B