0-Ph}RddlZddlmZddlmZmZ d ejdd Z d d Z dS) N)ndimage)_validate_interpolation_order_fix_ndimage_modereflect) reduce_funcc tjtt||| jdkrR fdt jdD}tjtj |d}ntj dk}tj |d}||} n6 ||d} n'#t$rtj||d } YnwxYw| S) a# Return the intensity profile of an image measured along a scan line. Parameters ---------- image : ndarray, shape (M, N[, C]) The image, either grayscale (2D array) or multichannel (3D array, where the final axis contains the channel information). src : array_like, shape (2,) The coordinates of the start point of the scan line. dst : array_like, shape (2,) The coordinates of the end point of the scan line. The destination point is *included* in the profile, in contrast to standard numpy indexing. linewidth : int, optional Width of the scan, perpendicular to the line order : int in {0, 1, 2, 3, 4, 5}, optional The order of the spline interpolation, default is 0 if image.dtype is bool and 1 otherwise. The order has to be in the range 0-5. See `skimage.transform.warp` for detail. mode : {'constant', 'nearest', 'reflect', 'mirror', 'wrap'}, optional How to compute any values falling outside of the image. cval : float, optional If `mode` is 'constant', what constant value to use outside the image. reduce_func : callable, optional Function used to calculate the aggregation of pixel values perpendicular to the profile_line direction when `linewidth` > 1. If set to None the unreduced array will be returned. Returns ------- return_value : array The intensity profile along the scan line. The length of the profile is the ceil of the computed length of the scan line. Examples -------- >>> x = np.array([[1, 1, 1, 2, 2, 2]]) >>> img = np.vstack([np.zeros_like(x), x, x, x, np.zeros_like(x)]) >>> img array([[0, 0, 0, 0, 0, 0], [1, 1, 1, 2, 2, 2], [1, 1, 1, 2, 2, 2], [1, 1, 1, 2, 2, 2], [0, 0, 0, 0, 0, 0]]) >>> profile_line(img, (2, 1), (2, 4)) array([1., 1., 2., 2.]) >>> profile_line(img, (1, 0), (1, 6), cval=4) array([1., 1., 1., 2., 2., 2., 2.]) The destination point is included in the profile, in contrast to standard numpy indexing. For example: >>> profile_line(img, (1, 0), (1, 6)) # The final point is out of bounds array([1., 1., 1., 2., 2., 2., 2.]) >>> profile_line(img, (1, 0), (1, 5)) # This accesses the full first row array([1., 1., 1., 2., 2., 2.]) For different reduce_func inputs: >>> profile_line(img, (1, 0), (1, 3), linewidth=3, reduce_func=np.mean) array([0.66666667, 0.66666667, 0.66666667, 1.33333333]) >>> profile_line(img, (1, 0), (1, 3), linewidth=3, reduce_func=np.max) array([1, 1, 1, 2]) >>> profile_line(img, (1, 0), (1, 3), linewidth=3, reduce_func=np.sum) array([2, 2, 2, 4]) The unreduced array will be returned when `reduce_func` is None or when `reduce_func` acts on each pixel value individually. >>> profile_line(img, (1, 2), (4, 2), linewidth=3, order=0, ... reduce_func=None) array([[1, 1, 2], [1, 1, 2], [1, 1, 2], [0, 0, 0]]) >>> profile_line(img, (1, 0), (1, 3), linewidth=3, reduce_func=np.sqrt) array([[1. , 1. , 0. ], [1. , 1. , 0. ], [1. , 1. , 0. ], [1.41421356, 1.41421356, 0. ]]) ) linewidthc \g|](}tjd|fdk)S).r prefilterordermodecval)ndimap_coordinates).0irimagerr perp_liness W/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/skimage/measure/profile.py z profile_line..ks[      c1f !)        r)rrrrr)axisN)arrr)rdtyper_line_profile_coordinatesndimrangeshapenp transposeasarrayrrflip TypeErrorapply_along_axis) rsrcdstr rrrr pixels intensitiesrs ` ``` @r profile_liner.s\~ *%+u = =E T " "D*3yIIIJ zQ         5;q>**    bj00)<<$ :%dQU    WV! $ $ $F  O%+f1555KK O O O-kvANNNKKK O s C!!!DDctj|tx\}}}tj|tx\}}}||z \}}tj||} t tjtj||dz} tj||| } tj||| } dz tj| zdz dz tj | zdz tj fd| D} tj fd| D}tj | |gS)aReturn the coordinates of the profile of an image along a scan line. Parameters ---------- src : 2-tuple of numeric scalar (float or int) The start point of the scan line. dst : 2-tuple of numeric scalar (float or int) The end point of the scan line. linewidth : int, optional Width of the scan, perpendicular to the line Returns ------- coords : array, shape (2, N, C), float The coordinates of the profile along the scan line. The length of the profile is the ceil of the computed length of the scan line. Notes ----- This is a utility method meant to be used internally by skimage functions. The destination point is included in the profile, in contrast to standard numpy indexing. )rrrcHg|]}tj|z |zSr$linspace)rrow_ir row_widths rrz-_line_profile_coordinates..A    K )59+.r6r) r$r&floatarctan2intceilhypotr3sincosstack)r*r+r src_rowsrc_coldst_rowdst_cold_rowd_colthetalengthline_colline_row perp_rows perp_colsr9r5s ` @@rr r s|0 Z59999GWsZ59999GWs9LE5 Jue $ $E %//!344 5 5F{7GV44H{7GV44H Q"&%..014IQ"&--/!3I     !   I      !   I 8Y * + ++r)rNrr )r) numpyr$scipyrr _shared.utilsrrmeanr.r r1rrrRs LLLLLLLL   @@@@@@F4,4,4,4,4,4,r