1-PhddlZddlZddlmZddlZddlmZddlm Z m Z ddl m Z ddl mZmZmZddlmZdd lmZd d lmZd d lmZmZmZd d lmZd dlmZmZd.dZd/dZ d0dZ! d1dZ"d2dZ#dZ$dZ%dZ&dZ'd3dZ(d.dZ)d4d!Z*d5d"Z+d5d#Z,d6d&Z-d7d)Z.d ddddej/ddfej/ej/d*d+Z0d5d,Z1d-Z2dS)8N)combinations_with_replacement)ndimage)spatialstatsgaussian)_supported_float_type safe_as_intwarn)integral_image) img_as_float)_hessian_matrix_det) _corner_fast_corner_moravec_corner_orientations)peak_local_max)_prepare_grayscale_input_2D_prepare_grayscale_input_nDconstantcPfdtjD}|S)aCompute derivatives in axis directions using the Sobel operator. Parameters ---------- image : ndarray Input image. mode : {'constant', 'reflect', 'wrap', 'nearest', 'mirror'}, optional How to handle values outside the image borders. cval : float, optional Used in conjunction with mode 'constant', the value outside the image boundaries. Returns ------- derivatives : list of ndarray Derivatives in each axis direction. c@g|]}tj|S))axismodecval)ndisobel).0irimagers V/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/skimage/feature/corner.py z(_compute_derivatives..'s9;< %ad666)rangendim)r!rr derivativess``` r"_compute_derivativesr(sI(@Eej@Q@QK r$rcc|dkr|jdkrtd|dvrtd|dtjs6t t |jkrtdt |}t|}|dkrt|}fd t|dD}|S) aCompute structure tensor using sum of squared differences. The (2-dimensional) structure tensor A is defined as:: A = [Arr Arc] [Arc Acc] which is approximated by the weighted sum of squared differences in a local window around each pixel in the image. This formula can be extended to a larger number of dimensions (see [1]_). Parameters ---------- image : ndarray Input image. sigma : float or array-like of float, optional Standard deviation used for the Gaussian kernel, which is used as a weighting function for the local summation of squared differences. If sigma is an iterable, its length must be equal to `image.ndim` and each element is used for the Gaussian kernel applied along its respective axis. mode : {'constant', 'reflect', 'wrap', 'nearest', 'mirror'}, optional How to handle values outside the image borders. cval : float, optional Used in conjunction with mode 'constant', the value outside the image boundaries. order : {'rc', 'xy'}, optional NOTE: 'xy' is only an option for 2D images, higher dimensions must always use 'rc' order. This parameter allows for the use of reverse or forward order of the image axes in gradient computation. 'rc' indicates the use of the first axis initially (Arr, Arc, Acc), whilst 'xy' indicates the usage of the last axis initially (Axx, Axy, Ayy). Returns ------- A_elems : list of ndarray Upper-diagonal elements of the structure tensor for each pixel in the input image. Examples -------- >>> from skimage.feature import structure_tensor >>> square = np.zeros((5, 5)) >>> square[2, 2] = 1 >>> Arr, Arc, Acc = structure_tensor(square, sigma=0.1, order='rc') >>> Acc array([[0., 0., 0., 0., 0.], [0., 1., 0., 1., 0.], [0., 4., 0., 4., 0.], [0., 1., 0., 1., 0.], [0., 0., 0., 0., 0.]]) See also -------- structure_tensor_eigenvalues References ---------- .. [1] https://en.wikipedia.org/wiki/Structure_tensor xyrz)Only "rc" order is supported for dim > 2.r)r+zorder z( is invalid. Must be either "rc" or "xy"z2sigma must have as many elements as image has axesrrcBg|]\}}t||zS)sigmarrr)rder0der1rrr0s r"r#z$structure_tensor..~sB D$ E4@@@r$) r& ValueErrornpisscalartuplelenrr(reversedr)r!r0rrorderr'A_elemss ``` r"structure_tensorr;.s z }}aDEEE L  Q%QQQRRR ;u  Ve  u:: # #TUU U ' . .E&u4dCCCK }}{++ 7 QGGG Nr$reflectcH ttj}|djdkr|dkrt d|dvrt d|t j|r |fjz}td|Drd nd }d tj dz tfd |D}t|||| }tjtjfi| j t fdt# D fdt# D t# } |dkrt%| } fdt'| dD} | S)a,Compute the Hessian via convolutions with Gaussian derivatives. In 2D, the Hessian matrix is defined as: H = [Hrr Hrc] [Hrc Hcc] which is computed by convolving the image with the second derivatives of the Gaussian kernel in the respective r- and c-directions. The implementation here also supports n-dimensional data. Parameters ---------- image : ndarray Input image. sigma : float or sequence of float, optional Standard deviation used for the Gaussian kernel, which sets the amount of smoothing in terms of pixel-distances. It is advised to not choose a sigma much less than 1.0, otherwise aliasing artifacts may occur. mode : {'constant', 'reflect', 'wrap', 'nearest', 'mirror'}, optional How to handle values outside the image borders. cval : float, optional Used in conjunction with mode 'constant', the value outside the image boundaries. order : {'rc', 'xy'}, optional This parameter allows for the use of reverse or forward order of the image axes in gradient computation. 'rc' indicates the use of the first axis initially (Hrr, Hrc, Hcc), whilst 'xy' indicates the usage of the last axis initially (Hxx, Hxy, Hyy) Returns ------- H_elems : list of ndarray Upper-diagonal elements of the hessian matrix for each pixel in the input image. In 2D, this will be a three element list containing [Hrr, Hrc, Hcc]. In nD, the list will contain ``(n**2 + n) / 2`` arrays. Fcopyrr++order='xy' is only supported for 2D images.r,unrecognized order: c3"K|] }|dkV dS)rN)rss r" z0_hessian_matrix_with_gaussian..s&--!A------r$drc3"K|] }|zV dS)NrC)rrDsq1_2s r"rEz0_hessian_matrix_with_gaussian..s'22q222222r$)r0rrtruncatec3FK|]}dg|zdgzdg|z dz zzVdS)rrNrC)rdr&s r"rEz0_hessian_matrix_with_gaussian..sCMMAA37aS=A3$(Q,#77MMMMMMr$c6g|]}|Sr9rC)rrL gaussian_r!orderss r"r#z1_hessian_matrix_with_gaussian..s,HHHq5q 222HHHr$cHg|]\}}||SrNrC)rax0ax1rP gradientsrQs r"r#z1_hessian_matrix_with_gaussian..sD C  )C.s 444r$)rr dtypeastyper&r3r4r5allmathsqrtr6dict functoolspartialrgaussian_filterr%r8r)r!r0rrr9 float_dtyperJ sigma_scaled common_kwargsaxesH_elemsrPrUr&rQrIs` @@@@@r"_hessian_matrix_with_gaussianrdsP   E' 44K LL5L 1 1E zA~~%4--FGGG L  777888 {5&5:%--u-----6qq3H  ! E2222E22222L|$THUUUM!#"5GGGGI :D MMMMtMMM M MFHHHHHHE$KKHHHI ;;D }}~~5dA>>G Nr$c" t|}t|j}||d}|jdkr|dkrt d|dvrt d||d}t d td |rt||||| St|||| }tj | t|j}|dkrt|} fd t|dD} | S)a Compute the Hessian matrix. In 2D, the Hessian matrix is defined as:: H = [Hrr Hrc] [Hrc Hcc] which is computed by convolving the image with the second derivatives of the Gaussian kernel in the respective r- and c-directions. The implementation here also supports n-dimensional data. Parameters ---------- image : ndarray Input image. sigma : float Standard deviation used for the Gaussian kernel, which is used as weighting function for the auto-correlation matrix. mode : {'constant', 'reflect', 'wrap', 'nearest', 'mirror'}, optional How to handle values outside the image borders. cval : float, optional Used in conjunction with mode 'constant', the value outside the image boundaries. order : {'rc', 'xy'}, optional For 2D images, this parameter allows for the use of reverse or forward order of the image axes in gradient computation. 'rc' indicates the use of the first axis initially (Hrr, Hrc, Hcc), whilst 'xy' indicates the usage of the last axis initially (Hxx, Hxy, Hyy). Images with higher dimension must always use 'rc' order. use_gaussian_derivatives : boolean, optional Indicates whether the Hessian is computed by convolving with Gaussian derivatives, or by a simple finite-difference operation. Returns ------- H_elems : list of ndarray Upper-diagonal elements of the hessian matrix for each pixel in the input image. In 2D, this will be a three element list containing [Hrr, Hrc, Hcc]. In nD, the list will contain ``(n**2 + n) / 2`` arrays. Notes ----- The distributive property of derivatives and convolutions allows us to restate the derivative of an image, I, smoothed with a Gaussian kernel, G, as the convolution of the image with the derivative of G. .. math:: \frac{\partial }{\partial x_i}(I * G) = I * \left( \frac{\partial }{\partial x_i} G \right) When ``use_gaussian_derivatives`` is ``True``, this property is used to compute the second order derivatives that make up the Hessian matrix. When ``use_gaussian_derivatives`` is ``False``, simple finite differences on a Gaussian-smoothed image are used instead. Examples -------- >>> from skimage.feature import hessian_matrix >>> square = np.zeros((5, 5)) >>> square[2, 2] = 4 >>> Hrr, Hrc, Hcc = hessian_matrix(square, sigma=0.1, order='rc', ... use_gaussian_derivatives=False) >>> Hrc array([[ 0., 0., 0., 0., 0.], [ 0., 1., 0., -1., 0.], [ 0., 0., 0., 0., 0.], [ 0., -1., 0., 1., 0.], [ 0., 0., 0., 0., 0.]]) Fr>rr+r@r,rANzuse_gaussian_derivatives currently defaults to False, but will change to True in a future version. Please specify this argument explicitly to maintain the current behavior)category stacklevel)r0rrr9r/cNg|]!\}}tj||"S)r)r4gradient)rrSrTrUs r"r#z"hessian_matrix..Ns@ C IcN---r$)rr rVrWr&r3r FutureWarningrdr r4rjr%r8r) r!r0rrr9use_gaussian_derivativesr_gaussian_filteredrbrcrUs @r"hessian_matrixrns]\   E' 44K LL5L 1 1E zA~~%4--FGGG L  777888'#(   C#       , TE    !e$TJJJ -..I   D }}~~5dA>>G Nr$Tcnt|}t|j}||d}|jdkr3|r1t |}t jt||Stt||d}t j |S)aCompute the approximate Hessian Determinant over an image. The 2D approximate method uses box filters over integral images to compute the approximate Hessian Determinant. Parameters ---------- image : ndarray The image over which to compute the Hessian Determinant. sigma : float, optional Standard deviation of the Gaussian kernel used for the Hessian matrix. approximate : bool, optional If ``True`` and the image is 2D, use a much faster approximate computation. This argument has no effect on 3D and higher images. Returns ------- out : array The array of the Determinant of Hessians. References ---------- .. [1] Herbert Bay, Andreas Ess, Tinne Tuytelaars, Luc Van Gool, "SURF: Speeded Up Robust Features" ftp://ftp.vision.ee.ethz.ch/publications/articles/eth_biwi_00517.pdf Notes ----- For 2D images when ``approximate=True``, the running time of this method only depends on size of the image. It is independent of `sigma` as one would expect. The downside is that the result for `sigma` less than `3` is not accurate, i.e., not similar to the result obtained if someone computed the Hessian and took its determinant. Fr>r)rl) rr rVrWr&r r4arrayr_symmetric_imagernlinalgdet)r!r0 approximater_integralhessian_mat_arrays r"hessian_matrix_detrwUsH   E' 44K LL5L 1 1E zQ;!%((x+He<<===, 5%% H H H  y}}.///r$ct|dkr{|\}}}tjdg|jR|j}||zdz |dd<tj|dz||z dz dzz}|dxx|z cc<|dxx|zcc<|St |}tj|ddddf}tt|j dz }tj ||j dz f|zS)aCompute eigenvalues from the upper-diagonal entries of a symmetric matrix. Parameters ---------- S_elems : list of ndarray The upper-diagonal elements of the matrix, as returned by `hessian_matrix` or `structure_tensor`. Returns ------- eigs : ndarray The eigenvalues of the matrix, in decreasing order. The eigenvalues are the leading dimension. That is, ``eigs[i, j, k]`` contains the ith-largest eigenvalue at position (j, k). rNrr.) r7r4emptyshaperVrZrqrreigvalshr6r%r& transpose)S_elemsM00M01M11eigshsqrtdetmatrices leading_axess r"_symmetric_compute_eigenvaluesrs$ 7||q S#xSY339/QQQ736cCi1_$::;; Q8 Q8 #G,,y!!(++C2I6U49q=1122 |D49q="2\"ABBBr$c(|d}tj|j|j|jfz|dj}t t t|jdD]$\}\}}|||d||f<|||d||f<%|S)aConvert the upper-diagonal elements of a matrix to the full symmetric matrix. Parameters ---------- S_elems : list of array The upper-diagonal elements of the matrix, as returned by `hessian_matrix` or `structure_tensor`. Returns ------- image : array An array of shape ``(M, N[, ...], image.ndim, image.ndim)``, containing the matrix corresponding to each coordinate. rrVr.)r4zerosr|r&rV enumeraterr%)rr!symmetric_imageidxrowcols r"rqrqs AJEh uz5:..gaj6FO%%eEJ&7&7;;66Zc3*1S# &)0S# && r$c t|S)aCompute eigenvalues of structure tensor. Parameters ---------- A_elems : list of ndarray The upper-diagonal elements of the structure tensor, as returned by `structure_tensor`. Returns ------- ndarray The eigenvalues of the structure tensor, in decreasing order. The eigenvalues are the leading dimension. That is, the coordinate [i, j, k] corresponds to the ith-largest eigenvalue at position (j, k). Examples -------- >>> from skimage.feature import structure_tensor >>> from skimage.feature import structure_tensor_eigenvalues >>> square = np.zeros((5, 5)) >>> square[2, 2] = 1 >>> A_elems = structure_tensor(square, sigma=0.1, order='rc') >>> structure_tensor_eigenvalues(A_elems)[0] array([[0., 0., 0., 0., 0.], [0., 2., 4., 2., 0.], [0., 4., 0., 4., 0.], [0., 2., 4., 2., 0.], [0., 0., 0., 0., 0.]]) See also -------- structure_tensor r)r:s r"structure_tensor_eigenvaluesrsD *' 2 22r$c t|S)aCompute eigenvalues of Hessian matrix. Parameters ---------- H_elems : list of ndarray The upper-diagonal elements of the Hessian matrix, as returned by `hessian_matrix`. Returns ------- eigs : ndarray The eigenvalues of the Hessian matrix, in decreasing order. The eigenvalues are the leading dimension. That is, ``eigs[i, j, k]`` contains the ith-largest eigenvalue at position (j, k). Examples -------- >>> from skimage.feature import hessian_matrix, hessian_matrix_eigvals >>> square = np.zeros((5, 5)) >>> square[2, 2] = 4 >>> H_elems = hessian_matrix(square, sigma=0.1, order='rc', ... use_gaussian_derivatives=False) >>> hessian_matrix_eigvals(H_elems)[0] array([[ 0., 0., 2., 0., 0.], [ 0., 1., 0., 1., 0.], [ 2., 0., -2., 0., 2.], [ 0., 1., 0., 1., 0.], [ 0., 0., 2., 0., 0.]]) r)rcs r"hessian_matrix_eigvalsrs< *' 2 22r$c t||||dd}t|\}}tjdd5dtjz tj||z||z z zcdddS#1swxYwYdS)a~Compute the shape index. The shape index, as defined by Koenderink & van Doorn [1]_, is a single valued measure of local curvature, assuming the image as a 3D plane with intensities representing heights. It is derived from the eigenvalues of the Hessian, and its value ranges from -1 to 1 (and is undefined (=NaN) in *flat* regions), with following ranges representing following shapes: .. table:: Ranges of the shape index and corresponding shapes. =================== ============= Interval (s in ...) Shape =================== ============= [ -1, -7/8) Spherical cup [-7/8, -5/8) Through [-5/8, -3/8) Rut [-3/8, -1/8) Saddle rut [-1/8, +1/8) Saddle [+1/8, +3/8) Saddle ridge [+3/8, +5/8) Ridge [+5/8, +7/8) Dome [+7/8, +1] Spherical cap =================== ============= Parameters ---------- image : (M, N) ndarray Input image. sigma : float, optional Standard deviation used for the Gaussian kernel, which is used for smoothing the input data before Hessian eigen value calculation. mode : {'constant', 'reflect', 'wrap', 'nearest', 'mirror'}, optional How to handle values outside the image borders cval : float, optional Used in conjunction with mode 'constant', the value outside the image boundaries. Returns ------- s : ndarray Shape index References ---------- .. [1] Koenderink, J. J. & van Doorn, A. J., "Surface shape and curvature scales", Image and Vision Computing, 1992, 10, 557-564. :DOI:`10.1016/0262-8856(92)90076-F` Examples -------- >>> from skimage.feature import shape_index >>> square = np.zeros((5, 5)) >>> square[2, 2] = 4 >>> s = shape_index(square, sigma=0.1) >>> s array([[ nan, nan, -0.5, nan, nan], [ nan, -0. , nan, -0. , nan], [-0.5, nan, -1. , nan, -0.5], [ nan, -0. , nan, -0. , nan], [ nan, nan, -0.5, nan, nan]]) r)F)r0rrr9rlignore)divideinvalidg@N)rnrr4errstatepiarctan)r!r0rrHl1l2s r" shape_indexr sD    !&    A$A & &FB Hh 7 7 7@@be ry"r'b2g)>???@@@@@@@@@@@@@@@@@@s-A88A<?A<ct|j}||d}t|||\}}t|||\}}t|||\}} ||dzz||dzzzd|z|z|zz } |dz|dzz} t j||} | dk} | | | | z | | <| S)aCompute Kitchen and Rosenfeld corner measure response image. The corner measure is calculated as follows:: (imxx * imy**2 + imyy * imx**2 - 2 * imxy * imx * imy) / (imx**2 + imy**2) Where imx and imy are the first and imxx, imxy, imyy the second derivatives. Parameters ---------- image : (M, N) ndarray Input image. mode : {'constant', 'reflect', 'wrap', 'nearest', 'mirror'}, optional How to handle values outside the image borders. cval : float, optional Used in conjunction with mode 'constant', the value outside the image boundaries. Returns ------- response : ndarray Kitchen and Rosenfeld response image. References ---------- .. [1] Kitchen, L., & Rosenfeld, A. (1982). Gray-level corner detection. Pattern recognition letters, 1(2), 95-102. :DOI:`10.1016/0167-8655(82)90020-4` Fr>r-rrr)r rVrWr(r4 zeros_like)r!rrr_imyimximxyimxximyyimyx numerator denominatorresponsemasks r"corner_kitchen_rosenfeldr[sB( 44K LL5L 1 1E#E4@@@HC%c4@@@JD$%c4@@@JD$sAv sAv -D30DDIq&36/K}U+666H ! Dt_{4'88HTN Or$k皙?ư>ct||d\}}}||z|dzz }||z} |dkr ||| dzzz } n d|z| |zz } | S)aCompute Harris corner measure response image. This corner detector uses information from the auto-correlation matrix A:: A = [(imx**2) (imx*imy)] = [Axx Axy] [(imx*imy) (imy**2)] [Axy Ayy] Where imx and imy are first derivatives, averaged with a gaussian filter. The corner measure is then defined as:: det(A) - k * trace(A)**2 or:: 2 * det(A) / (trace(A) + eps) Parameters ---------- image : (M, N) ndarray Input image. method : {'k', 'eps'}, optional Method to compute the response image from the auto-correlation matrix. k : float, optional Sensitivity factor to separate corners from edges, typically in range `[0, 0.2]`. Small values of k result in detection of sharp corners. eps : float, optional Normalisation factor (Noble's corner measure). sigma : float, optional Standard deviation used for the Gaussian kernel, which is used as weighting function for the auto-correlation matrix. Returns ------- response : ndarray Harris response image. References ---------- .. [1] https://en.wikipedia.org/wiki/Corner_detection Examples -------- >>> from skimage.feature import corner_harris, corner_peaks >>> square = np.zeros([10, 10]) >>> square[2:8, 2:8] = 1 >>> square.astype(int) array([[0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 1, 1, 1, 1, 1, 1, 0, 0], [0, 0, 1, 1, 1, 1, 1, 1, 0, 0], [0, 0, 1, 1, 1, 1, 1, 1, 0, 0], [0, 0, 1, 1, 1, 1, 1, 1, 0, 0], [0, 0, 1, 1, 1, 1, 1, 1, 0, 0], [0, 0, 1, 1, 1, 1, 1, 1, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0]]) >>> corner_peaks(corner_harris(square), min_distance=1) array([[2, 2], [2, 7], [7, 2], [7, 7]]) r)rOrr)r;) r!methodrepsr0ArrArcAccdetAtraceArs r" corner_harrisrspB%UE>>>MCc 9sAv D 3YF }}!fai-'t8v|, Or$ct||d\}}}||ztj||z dzd|dzzzz dz }|S)a`Compute Shi-Tomasi (Kanade-Tomasi) corner measure response image. This corner detector uses information from the auto-correlation matrix A:: A = [(imx**2) (imx*imy)] = [Axx Axy] [(imx*imy) (imy**2)] [Axy Ayy] Where imx and imy are first derivatives, averaged with a gaussian filter. The corner measure is then defined as the smaller eigenvalue of A:: ((Axx + Ayy) - sqrt((Axx - Ayy)**2 + 4 * Axy**2)) / 2 Parameters ---------- image : (M, N) ndarray Input image. sigma : float, optional Standard deviation used for the Gaussian kernel, which is used as weighting function for the auto-correlation matrix. Returns ------- response : ndarray Shi-Tomasi response image. References ---------- .. [1] https://en.wikipedia.org/wiki/Corner_detection Examples -------- >>> from skimage.feature import corner_shi_tomasi, corner_peaks >>> square = np.zeros([10, 10]) >>> square[2:8, 2:8] = 1 >>> square.astype(int) array([[0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 1, 1, 1, 1, 1, 1, 0, 0], [0, 0, 1, 1, 1, 1, 1, 1, 0, 0], [0, 0, 1, 1, 1, 1, 1, 1, 0, 0], [0, 0, 1, 1, 1, 1, 1, 1, 0, 0], [0, 0, 1, 1, 1, 1, 1, 1, 0, 0], [0, 0, 1, 1, 1, 1, 1, 1, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0]]) >>> corner_peaks(corner_shi_tomasi(square), min_distance=1) array([[2, 2], [2, 7], [7, 2], [7, 7]]) r)rOr)r;r4rZ)r!r0rrrrs r"corner_shi_tomasirsYl%UE>>>MCcsbgsSyQ&6S!V&CDDDIH Or$ct||d\}}}||z|dzz }||z}tj||j}tj|}|dk} || || z || <d|| z|| dzz || <||fS)uCompute Foerstner corner measure response image. This corner detector uses information from the auto-correlation matrix A:: A = [(imx**2) (imx*imy)] = [Axx Axy] [(imx*imy) (imy**2)] [Axy Ayy] Where imx and imy are first derivatives, averaged with a gaussian filter. The corner measure is then defined as:: w = det(A) / trace(A) (size of error ellipse) q = 4 * det(A) / trace(A)**2 (roundness of error ellipse) Parameters ---------- image : (M, N) ndarray Input image. sigma : float, optional Standard deviation used for the Gaussian kernel, which is used as weighting function for the auto-correlation matrix. Returns ------- w : ndarray Error ellipse sizes. q : ndarray Roundness of error ellipse. References ---------- .. [1] Förstner, W., & Gülch, E. (1987, June). A fast operator for detection and precise location of distinct points, corners and centres of circular features. In Proc. ISPRS intercommission conference on fast processing of photogrammetric data (pp. 281-305). https://cseweb.ucsd.edu/classes/sp02/cse252/foerstner/foerstner.pdf .. [2] https://en.wikipedia.org/wiki/Corner_detection Examples -------- >>> from skimage.feature import corner_foerstner, corner_peaks >>> square = np.zeros([10, 10]) >>> square[2:8, 2:8] = 1 >>> square.astype(int) array([[0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 1, 1, 1, 1, 1, 1, 0, 0], [0, 0, 1, 1, 1, 1, 1, 1, 0, 0], [0, 0, 1, 1, 1, 1, 1, 1, 0, 0], [0, 0, 1, 1, 1, 1, 1, 1, 0, 0], [0, 0, 1, 1, 1, 1, 1, 1, 0, 0], [0, 0, 1, 1, 1, 1, 1, 1, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0]]) >>> w, q = corner_foerstner(square) >>> accuracy_thresh = 0.5 >>> roundness_thresh = 0.3 >>> foerstner = (q > roundness_thresh) * (w > accuracy_thresh) * w >>> corner_peaks(foerstner, min_distance=1) array([[2, 2], [2, 7], [7, 2], [7, 7]]) r)rOrrrr)r;r4rrV) r!r0rrrrrwqrs r"corner_foerstnerrsD%UE>>>MCc 9sAv D 3YF e4:...A aA Q;D4j6$<'AdG$t*nvd|q00AdG a4Kr$ 333333?cnt|}tj|}t|||}|S)aExtract FAST corners for a given image. Parameters ---------- image : (M, N) ndarray Input image. n : int, optional Minimum number of consecutive pixels out of 16 pixels on the circle that should all be either brighter or darker w.r.t testpixel. A point c on the circle is darker w.r.t test pixel p if `Ic < Ip - threshold` and brighter if `Ic > Ip + threshold`. Also stands for the n in `FAST-n` corner detector. threshold : float, optional Threshold used in deciding whether the pixels on the circle are brighter, darker or similar w.r.t. the test pixel. Decrease the threshold when more corners are desired and vice-versa. Returns ------- response : ndarray FAST corner response image. References ---------- .. [1] Rosten, E., & Drummond, T. (2006, May). Machine learning for high-speed corner detection. In European conference on computer vision (pp. 430-443). Springer, Berlin, Heidelberg. :DOI:`10.1007/11744023_34` http://www.edwardrosten.com/work/rosten_2006_machine.pdf .. [2] Wikipedia, "Features from accelerated segment test", https://en.wikipedia.org/wiki/Features_from_accelerated_segment_test Examples -------- >>> from skimage.feature import corner_fast, corner_peaks >>> square = np.zeros((12, 12)) >>> square[3:9, 3:9] = 1 >>> square.astype(int) array([[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, 0, 0, 0, 0, 0, 0], [0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0], [0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0], [0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0], [0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0], [0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0], [0, 0, 0, 1, 1, 1, 1, 1, 1, 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, 0, 0, 0, 0, 0, 0, 0, 0, 0]]) >>> corner_peaks(corner_fast(square, 9), min_distance=1) array([[3, 3], [3, 8], [8, 3], [8, 8]]) )rr4ascontiguousarrayr)r!n thresholdrs r" corner_fastrps8t ( . .E   ' 'EE1i00H Or$ Gz?c |dz dz}t|j}||d}tj||dd}t ||z}tjd| }tjd| }tjd | }tjd | } |dzdz } tj d|z | | } tj || | } tj | |dz| |dzf\} }tj || }t|D]\}\}}||z dz }||zdz}||z dz }||zdz}|||||f}t|dd \}}||zdd dd f}||zdd dd f}||zdd dd f}tj|}tj|}tj|}tj||z} tj|| z}!tj||z}"tj|| z}#tj||z}$tj|| z}%||d <| x|d<|d<||d<||d <|x|d<|d<||d<|!|"z |$|#z f|dd<|%|"z| |#zf| dd< tj||}&tj|| }'n:#tjj$r#tjtjf||ddf<YwxYw| |&dz }(||&dz })| |'dz }*||'dz }+|)|)z},|)|(z}-|(|(z}.|+|+z}/|+|*z}0|*|*z}1tj||.zd|z|-zz ||,zz}2tj||1zd|z|0zz||/zz}3|2tjdkr%|3tjdkr tj}4n|2dkr tj}4n|3|2z }4t+|4| kt+|4| kz }5|5d kr||&dz||&dzf||ddf<?|5dkr!tjtjf||ddf<f|5dkr||'dz||'dzf||ddf<||z}|S)uDetermine subpixel position of corners. A statistical test decides whether the corner is defined as the intersection of two edges or a single peak. Depending on the classification result, the subpixel corner location is determined based on the local covariance of the grey-values. If the significance level for either statistical test is not sufficient, the corner cannot be classified, and the output subpixel position is set to NaN. Parameters ---------- image : (M, N) ndarray Input image. corners : (K, 2) ndarray Corner coordinates `(row, col)`. window_size : int, optional Search window size for subpixel estimation. alpha : float, optional Significance level for corner classification. Returns ------- positions : (K, 2) ndarray Subpixel corner positions. NaN for "not classified" corners. References ---------- .. [1] Förstner, W., & Gülch, E. (1987, June). A fast operator for detection and precise location of distinct points, corners and centres of circular features. In Proc. ISPRS intercommission conference on fast processing of photogrammetric data (pp. 281-305). https://cseweb.ucsd.edu/classes/sp02/cse252/foerstner/foerstner.pdf .. [2] https://en.wikipedia.org/wiki/Corner_detection Examples -------- >>> from skimage.feature import corner_harris, corner_peaks, corner_subpix >>> img = np.zeros((10, 10)) >>> img[:5, :5] = 1 >>> img[5:, 5:] = 1 >>> img.astype(int) array([[1, 1, 1, 1, 1, 0, 0, 0, 0, 0], [1, 1, 1, 1, 1, 0, 0, 0, 0, 0], [1, 1, 1, 1, 1, 0, 0, 0, 0, 0], [1, 1, 1, 1, 1, 0, 0, 0, 0, 0], [1, 1, 1, 1, 1, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 1, 1, 1, 1, 1], [0, 0, 0, 0, 0, 1, 1, 1, 1, 1], [0, 0, 0, 0, 0, 1, 1, 1, 1, 1], [0, 0, 0, 0, 0, 1, 1, 1, 1, 1], [0, 0, 0, 0, 0, 1, 1, 1, 1, 1]]) >>> coords = corner_peaks(corner_harris(img), min_distance=2) >>> coords_subpix = corner_subpix(img, coords, window_size=7) >>> coords_subpix array([[4.5, 4.5]]) rrFr>rr) pad_widthrconstant_values)rrr)rr-rz)rr)rr)rr)rrN)r rVrWr4padr rrfisfmgridrrr(sumrrsolve LinAlgErrornanspacinginfint)6r!corners window_sizealphawextr_N_dotN_edgeb_dotb_edge redundancy t_crit_dot t_crit_edgeyxcorners_subpixr y0x0minymaxyminxmaxxwindowwinywinx winx_winx winx_winy winy_winyAxxAxyAyybxx_xbxx_ybxy_xbxy_ybyy_xbyy_yest_dotest_edgery_dotrx_dotry_edgerx_edgerxx_dotrxy_dotryy_dotrxx_edgerxy_edgeryy_edgevar_dotvar_edget corner_classs6 r" corner_subpixr sx !O !D' 44K LL5L 1 1E F5Dz1 M M ME'D.))G HV; / / /E XfK 0 0 0F HT - - -E Xd+ . . .Fa!#JQY J??J'++eZ<>>N ))XFXF 8BDy1}Dy1}Dy1}Dy1}tDy$t)+,)&zJJJ dD[!B$"*- D[!B$"*- D[!B$"*- fYfYfYy1}%%y1}%%y1}%%y1}%%y1}%%y1}%%d %(D(d eDkd t &))t vd|t 5=%%-/aaaEM55=0qqq  iooeU33Gyvv66HHy$   #%626>N1aaa4 H  WQZWQZhqk/hqk/6/6/6/W$W$W$&  !i-'"9 9IR R   RZ]] " "x"*Q--'?'?AA \\AA7"A1{?++c!j..A.AA 2  #% ?BO#CN1aaa4 Q  #%626>N1aaa4 Q  #% #3R(1+5E#EN1aaa4 dN sAL3MM)num_peaks_per_labelp_normc  Rtj|rd}t|||||tj|||  } t | rt j| } t} t| D]K\}}|| vrB| ||| }| || |Ltj | t| dd|} |r| Stj|t}d|t| j<|S)aFind peaks in corner measure response image. This differs from `skimage.feature.peak_local_max` in that it suppresses multiple connected peaks with the same accumulator value. Parameters ---------- image : (M, N) ndarray Input image. min_distance : int, optional The minimal allowed distance separating peaks. * : * See :py:meth:`skimage.feature.peak_local_max`. p_norm : float Which Minkowski p-norm to use. Should be in the range [1, inf]. A finite large p may cause a ValueError if overflow can occur. ``inf`` corresponds to the Chebyshev distance and 2 to the Euclidean distance. Returns ------- output : ndarray or ndarray of bools * If `indices = True` : (row, column, ...) coordinates of peaks. * If `indices = False` : Boolean array shaped like `image`, with peaks represented by True values. See also -------- skimage.feature.peak_local_max Notes ----- .. versionchanged:: 0.18 The default value of `threshold_rel` has changed to None, which corresponds to letting `skimage.feature.peak_local_max` decide on the default. This is equivalent to `threshold_rel=0`. The `num_peaks` limit is applied before suppression of connected peaks. To limit the number of peaks after suppression, set `num_peaks=np.inf` and post-process the output of this function. Examples -------- >>> from skimage.feature import peak_local_max >>> response = np.zeros((5, 5)) >>> response[2:4, 2:4] = 1 >>> response array([[0., 0., 0., 0., 0.], [0., 0., 0., 0., 0.], [0., 0., 1., 1., 0.], [0., 0., 1., 1., 0.], [0., 0., 0., 0., 0.]]) >>> peak_local_max(response) array([[2, 2], [2, 3], [3, 2], [3, 3]]) >>> corner_peaks(response) array([[2, 2]]) N) min_distance threshold_abs threshold_relexclude_border num_peaks footprintlabelsr )rprrirT)r4isinfrrr7rcKDTreesetrquery_ball_pointremoveupdatedeleter6rboolT)r!rrrrindicesrrrr r coordstreerejected_peaks_indicesrpoint candidatespeakss r" corner_peaksr'fsEX x   !##%&/   F 6{{ Vv&&!$#F++ : :JC000!225LF2SS !!#&&&&--j99965)?#@#@qIII*9*U M%t , , ,E!E%// Lr$ct|}t|j}||d}t t j||S)aCompute Moravec corner measure response image. This is one of the simplest corner detectors and is comparatively fast but has several limitations (e.g. not rotation invariant). Parameters ---------- image : (M, N) ndarray Input image. window_size : int, optional Window size. Returns ------- response : ndarray Moravec response image. References ---------- .. [1] https://en.wikipedia.org/wiki/Corner_detection Examples -------- >>> from skimage.feature import corner_moravec >>> square = np.zeros([7, 7]) >>> square[3, 3] = 1 >>> square.astype(int) array([[0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 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, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0]]) >>> corner_moravec(square).astype(int) array([[0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0], [0, 0, 1, 1, 1, 0, 0], [0, 0, 1, 2, 1, 0, 0], [0, 0, 1, 1, 1, 0, 0], [0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0]]) Fr>)rr rVrWrr4r)r!rr_s r"corner_moravecr)sQX   E' 44K LL5L 1 1E 2/66 D DDr$cBt|}t|||S)a Compute the orientation of corners. The orientation of corners is computed using the first order central moment i.e. the center of mass approach. The corner orientation is the angle of the vector from the corner coordinate to the intensity centroid in the local neighborhood around the corner calculated using first order central moment. Parameters ---------- image : (M, N) array Input grayscale image. corners : (K, 2) array Corner coordinates as ``(row, col)``. mask : 2D array Mask defining the local neighborhood of the corner used for the calculation of the central moment. Returns ------- orientations : (K, 1) array Orientations of corners in the range [-pi, pi]. References ---------- .. [1] Ethan Rublee, Vincent Rabaud, Kurt Konolige and Gary Bradski "ORB : An efficient alternative to SIFT and SURF" http://www.vision.cs.chubu.ac.jp/CV-R/pdf/Rublee_iccv2011.pdf .. [2] Paul L. Rosin, "Measuring Corner Properties" http://users.cs.cf.ac.uk/Paul.Rosin/corner2.pdf Examples -------- >>> from skimage.morphology import octagon >>> from skimage.feature import (corner_fast, corner_peaks, ... corner_orientations) >>> square = np.zeros((12, 12)) >>> square[3:9, 3:9] = 1 >>> square.astype(int) array([[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, 0, 0, 0, 0, 0, 0], [0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0], [0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0], [0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0], [0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0], [0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0], [0, 0, 0, 1, 1, 1, 1, 1, 1, 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, 0, 0, 0, 0, 0, 0, 0, 0, 0]]) >>> corners = corner_peaks(corner_fast(square, 9), min_distance=1) >>> corners array([[3, 3], [3, 8], [8, 3], [8, 8]]) >>> orientations = corner_orientations(square, corners, octagon(3, 2)) >>> np.rad2deg(orientations) array([ 45., 135., -45., -135.]) )rr)r!rrs r"corner_orientationsr+ s$~ ( . .E w 5 55r$)rr)rrrr))rr<rr))rrrr)N)rT)rrr)rrrr)r)rr)rr)3r\rY itertoolsrnumpyr4scipyrrrr_shared.filtersr _shared.utilsr r r transformr utilr_hessian_det_appxr corner_cyrrrpeakrrrr(r;rdrnrwrrqrrrrrrrrr rr'r)r+rCr$r"r6s 333333 &&&&&&DDDDDDDDDD&&&&&&222222JJJJJJJJJJ JJJJJJJJ6UUUUpXXXXxSWqqqqh.0.0.0.0bCCCD8"3"3"3J333BN@N@N@N@b0000fMMMM`;;;;|QQQQh>>>>Brrrrn f p 6pppppf/E/E/E/Ed@6@6@6@6@6r$