B \i@s0ddlZddlZddlmZddZddZdS)N)warncCs(|}xt|jD]}|j|d}qW|S)a0Integral image / summed area table. The integral image contains the sum of all elements above and to the left of it, i.e.: .. math:: S[m, n] = \sum_{i \leq m} \sum_{j \leq n} X[i, j] Parameters ---------- image : ndarray Input image. Returns ------- S : ndarray Integral image/summed area table of same shape as input image. References ---------- .. [1] F.C. Crow, "Summed-area tables for texture mapping," ACM SIGGRAPH Computer Graphics, vol. 18, 1984, pp. 207-212. )Zaxis)rangendimZcumsum)ZimageSir9lib/python3.7/site-packages/skimage/transform/integral.pyintegral_imagesr c sZtttt|}jd}j}t||dg}dk}|dk}||||||||}t|dkrtdt|}dj}t t |ddd} xt |D]} t | dd | } dd| Ddt fd dt |D|td|fd dt |D7}qW|S) aUse an integral image to integrate over a given window. Parameters ---------- ii : ndarray Integral image. start : List of tuples, each tuple of length equal to dimension of `ii` Coordinates of top left corner of window(s). Each tuple in the list contains the starting row, col, ... index i.e `[(row_win1, col_win1, ...), (row_win2, col_win2,...), ...]`. end : List of tuples, each tuple of length equal to dimension of `ii` Coordinates of bottom right corner of window(s). Each tuple in the list containing the end row, col, ... index i.e `[(row_win1, col_win1, ...), (row_win2, col_win2, ...), ...]`. Returns ------- S : scalar or ndarray Integral (sum) over the given window(s). Examples -------- >>> arr = np.ones((5, 6), dtype=np.float) >>> ii = integral_image(arr) >>> integrate(ii, (1, 0), (1, 2)) # sum from (1, 0) to (1, 2) array([ 3.]) >>> integrate(ii, [(3, 3)], [(4, 5)]) # sum from (3, 3) to (4, 5) array([ 6.]) >>> # sum from (1, 0) to (1, 2) and from (3, 3) to (4, 5) >>> integrate(ii, [(1, 0), (3, 3)], [(1, 2), (4, 5)]) array([ 3., 6.]) rz1end coordinates must be greater or equal to startrNcSsg|] }|dkqS)1r).0bitrrr wszintegrate..cs&g|]}t|ddkqS)r r)npany)r r) bool_maskstartrr r{scs,g|]$}|s$t|ndqS)r)tuple)r r)bad corner_pointsiisignrr rs)rZ atleast_2dZarrayshapeZtiler IndexErrorZzerosrlenbinrzfillsuminvert) rrendZrowsZ total_shapeZstart_negativesZ end_negativesrZbit_permwidthrZbinaryr)rrrrrrr integrate's6"           r$)Znumpyr collectionsZ _shared.utilsrr r$rrrr s