ó  ‰\c@s>ddlZddlZddlmZd„Zd„ZdS(iÿÿÿÿNi(twarncCs6|}x)t|jƒD]}|jd|ƒ}qW|S(s0Integral 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. taxis(trangetndimtcumsum(timagetSti((s9lib/python2.7/site-packages/skimage/transform/integral.pytintegral_imagesc Cstjtj|ƒƒ}tjtj|ƒƒ}|jd}|j}tj||dgƒ}|dk}|dk}|||||}|||||}tj||dkƒrÌtdƒ‚ntj|ƒ}d|j}t t |dƒdƒ} xút |ƒD]ì} t | ƒdj | ƒ} g| D]} | dk^q5} dt | ƒ}gt |ƒD]'}tj||d| dkƒ^qj}|tj| ƒ|d| }|gt |ƒD].}||sî||t||ƒnd^qÆ7}qW|S(sµUse 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.]) iis1end coordinates must be greater or equal to startit1iÿÿÿÿ(tnpt atleast_2dtarraytshapettiletanyt IndexErrortzerosRtlentbinRtzfilltsumtinvertttuple(tiitstarttendtrowst total_shapetstart_negativest end_negativesRtbit_permtwidthRtbinarytbitt bool_masktsigntrtbadt corner_points((s9lib/python2.7/site-packages/skimage/transform/integral.pyt integrate's6"          7C(tnumpyR t collectionst _shared.utilsRRR((((s9lib/python2.7/site-packages/skimage/transform/integral.pyts