1-PhddlZddlmZddlmZmZmZddl m Z ddl m Z ddl mZmZmZmZmZmZmZeeefZd Z d"d Ze d"dddZ d#dZd$dZdZ d%dZdZejfdZdZ d&dZ!dZ"dZ#e d'ddddddZ$d Z% d(ddd!Z&dS))N)ndimage)SimilarityTransformAffineTransformProjectiveTransform) _warp_fast) block_reduce)get_bound_method_class safe_as_intwarnconvert_to_float_to_ndimage_mode_validate_interpolation_orderchannel_as_last_axisc*t|}t|}|j}||jkr&|d||jz zz }t j||}n:||jdz kr||jdfz}n||jkrt d||fS)a'Validate resize output shape according to input image. Parameters ---------- image: ndarray Image to be resized. output_shape: iterable Size of the generated output image `(rows, cols[, ...][, dim])`. If `dim` is not provided, the number of channels is preserved. Returns ------- image: ndarray The input image, but with additional singleton dimensions appended in the case where ``len(output_shape) > input.ndim``. output_shape: tuple The output image converted to tuple. Raises ------ ValueError: If output_shape length is smaller than the image number of dimensions Notes ----- The input image is reshaped if its number of dimensions is not equal to output_shape_length. rrzIoutput_shape length cannot be smaller than the image number of dimensions)tuplelenshapendimnpreshape ValueError)image output_shape output_ndim input_shapes X/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/skimage/transform/_warps.py_preprocess_resize_output_shaper!s>&&Ll##K+KUZt{UZ788  5+..  Q & &#u{2&88 uz ! ! )   , reflectTFc  t||\}}|j} |j} | tjkr|tj}|Y| tk oLtj| tj o|dk o&tdt|| D}| tkr|rtdtj | |} t| |}|dkrt||}t!|} |r|tjd| dz dz }n~tj|tj| z}tj |dkrtdtj |dk| dkzrt)dt+j|||| } n|} d | D}t+j| ||| |d }t1||||||S) a Resize image to match a certain size. Performs interpolation to up-size or down-size N-dimensional images. Note that anti-aliasing should be enabled when down-sizing images to avoid aliasing artifacts. For downsampling with an integer factor also see `skimage.transform.downscale_local_mean`. Parameters ---------- image : ndarray Input image. output_shape : iterable Size of the generated output image `(rows, cols[, ...][, dim])`. If `dim` is not provided, the number of channels is preserved. In case the number of input channels does not equal the number of output channels a n-dimensional interpolation is applied. Returns ------- resized : ndarray Resized version of the input. Other parameters ---------------- order : int, 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', 'edge', 'symmetric', 'reflect', 'wrap'}, optional Points outside the boundaries of the input are filled according to the given mode. Modes match the behaviour of `numpy.pad`. cval : float, optional Used in conjunction with mode 'constant', the value outside the image boundaries. clip : bool, optional Whether to clip the output to the range of values of the input image. This is enabled by default, since higher order interpolation may produce values outside the given input range. preserve_range : bool, optional Whether to keep the original range of values. Otherwise, the input image is converted according to the conventions of `img_as_float`. Also see https://scikit-image.org/docs/dev/user_guide/data_types.html anti_aliasing : bool, optional Whether to apply a Gaussian filter to smooth the image prior to downsampling. It is crucial to filter when downsampling the image to avoid aliasing artifacts. If not specified, it is set to True when downsampling an image whose data type is not bool. It is also set to False when using nearest neighbor interpolation (``order`` == 0) with integer input data type. anti_aliasing_sigma : {float, tuple of floats}, optional Standard deviation for Gaussian filtering used when anti-aliasing. By default, this value is chosen as (s - 1) / 2 where s is the downsampling factor, where s > 1. For the up-size case, s < 1, no anti-aliasing is performed prior to rescaling. Notes ----- Modes 'reflect' and 'symmetric' are similar, but differ in whether the edge pixels are duplicated during the reflection. As an example, if an array has values [0, 1, 2] and was padded to the right by four values using symmetric, the result would be [0, 1, 2, 2, 1, 0, 0], while for reflect it would be [0, 1, 2, 1, 0, 1, 2]. Examples -------- >>> from skimage import data >>> from skimage.transform import resize >>> image = data.camera() >>> resize(image, (100, 100)).shape (100, 100) Nrc3(K|] \}}||kVdSN).0xys r zresize..s*EEdaAEEEEEEEr"z.anti_aliasing must be False for boolean imagesrr zFAnti-aliasing standard deviation must be greater than or equal to zerozWAnti-aliasing standard deviation greater than zero but not down-sampling along all axes)cvalmodecg|]}d|z Srr')r(fs r zresize..s+++aAE+++r"T)orderr-r, grid_mode)r!rdtyperfloat16astypefloat32bool issubdtypeintegeranyziprdividerrrmaximum atleast_1d ones_liker ndigaussian_filterzoom_clip_warp_output)rrr1r-r,clippreserve_range anti_aliasinganti_aliasing_sigmar input_typefactorsndi_modefiltered zoom_factorsouts r resizerNGsRh:%NNE<+KJRZ RZ((d" " F]:rz::IuzJ FEEc, &D&DEEEEE  TmIJJJi \22G )*e < >> from skimage import data >>> from skimage.transform import rescale >>> image = data.camera() >>> rescale(image, 0.1).shape (51, 51) >>> rescale(image, 0.5).shape (256, 256) Nrz4Supply a single scale, or one value per spatial axisr)r1r-r,rDrErFrG) rr>rrr concatenateasarrayrr=roundrN) rscaler1r-r,rDrErFrGrO multichannel orig_shapers r rescalerWs v M% Et+L 5zzA~~ XUuz!9!9 ": ZZ5:>99VWW W  1NEA3<00EEK((J:bhuz'9::A>>L*%b> R     %#/    r"constantc |jd|jd} } |jtjkr|tj}|tj| | fdz dz }ntj|}t|} ttj |} t| } | | z| z}d}|rtjddgd| dz g| dz | dz g| dz dgg}| |}|dddf }|dddf }|dddf }|dddf }||z dz}||z dz}tj ||f}||f}t|}||z}d|jd <t|||||||| S) a Rotate image by a certain angle around its center. Parameters ---------- image : ndarray Input image. angle : float Rotation angle in degrees in counter-clockwise direction. resize : bool, optional Determine whether the shape of the output image will be automatically calculated, so the complete rotated image exactly fits. Default is False. center : iterable of length 2 The rotation center. If ``center=None``, the image is rotated around its center, i.e. ``center=(cols / 2 - 0.5, rows / 2 - 0.5)``. Please note that this parameter is (cols, rows), contrary to normal skimage ordering. Returns ------- rotated : ndarray Rotated version of the input. Other parameters ---------------- order : int, 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', 'edge', 'symmetric', 'reflect', 'wrap'}, optional Points outside the boundaries of the input are filled according to the given mode. Modes match the behaviour of `numpy.pad`. cval : float, optional Used in conjunction with mode 'constant', the value outside the image boundaries. clip : bool, optional Whether to clip the output to the range of values of the input image. This is enabled by default, since higher order interpolation may produce values outside the given input range. preserve_range : bool, optional Whether to keep the original range of values. Otherwise, the input image is converted according to the conventions of `img_as_float`. Also see https://scikit-image.org/docs/dev/user_guide/data_types.html Notes ----- Modes 'reflect' and 'symmetric' are similar, but differ in whether the edge pixels are duplicated during the reflection. As an example, if an array has values [0, 1, 2] and was padded to the right by four values using symmetric, the result would be [0, 1, 2, 2, 1, 0, 0], while for reflect it would be [0, 1, 2, 1, 0, 1, 2]. Examples -------- >>> from skimage import data >>> from skimage.transform import rotate >>> image = data.camera() >>> rotate(image, 2).shape (512, 512) >>> rotate(image, 2, resize=True).shape (530, 530) >>> rotate(image, 90, resize=True).shape (512, 512) rrNg@?) translation)rotation)rrrr )rr1r-r,rDrE)rr3rr4r5r6arrayrRrdeg2radinverseminmaxaroundparamswarp)ranglerNcenterr1r-r,rDrErowscolstform1tform2tform3tformrcornersmincminrmaxcmaxrout_rowsout_colsr[tform4s r rotateruJs\QQ$D {bj   RZ((~4,''#-3F## V 4 4 4F "*U*;*; < < >> a = np.arange(15).reshape(3, 5) >>> a array([[ 0, 1, 2, 3, 4], [ 5, 6, 7, 8, 9], [10, 11, 12, 13, 14]]) >>> downscale_local_mean(a, (2, 3)) array([[3.5, 4. ], [5.5, 4.5]]) )r rmean)rrIr,rDs r downscale_local_meanrxsX w 6 66r"c|j\}}|\}}tj||z dz||z dzz} |dz tjdz}||tj| |z zztj||z ||z z} || tj| zz|d<|| tj| zz|d<|S)Nr ).r).r)Trsqrtlogexparctan2cossin) xyrfr\strengthradiusr)r*x0y0rhothetas r _swirl_mappingrs 4DAq FB '1r6a-1r6a-/ 0 0CaZ"&)) #F x"&#"7"77 7"*QVQQSV:T:T TEcBF5MM))BvJcBF5MM))BvJ Ir"dc |-tj|jddddddz }||||d} t|t| ||||| |  S)aZPerform a swirl transformation. Parameters ---------- image : ndarray Input image. center : (column, row) tuple or (2,) ndarray, optional Center coordinate of transformation. strength : float, optional The amount of swirling applied. radius : float, optional The extent of the swirl in pixels. The effect dies out rapidly beyond `radius`. rotation : float, optional Additional rotation applied to the image. Returns ------- swirled : ndarray Swirled version of the input. Other parameters ---------------- output_shape : tuple (rows, cols), optional Shape of the output image generated. By default the shape of the input image is preserved. order : int, 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', 'edge', 'symmetric', 'reflect', 'wrap'}, optional Points outside the boundaries of the input are filled according to the given mode, with 'reflect' used as the default. Modes match the behaviour of `numpy.pad`. cval : float, optional Used in conjunction with mode 'constant', the value outside the image boundaries. clip : bool, optional Whether to clip the output to the range of values of the input image. This is enabled by default, since higher order interpolation may produce values outside the given input range. preserve_range : bool, optional Whether to keep the original range of values. Otherwise, the input image is converted according to the conventions of `img_as_float`. Also see https://scikit-image.org/docs/dev/user_guide/data_types.html Nr r)rfr\rr)map_argsrr1r-r,rDrE)rr]rrdr) rrfrrr\rr1r-r,rDrE warp_argss r swirlrsz~%+&&rr*44R4014 I  !   %    r"cj|jdkr |ddddtjf|dd<dS||dd<dS)a5Copy b into each color layer of a, such that:: a[:,:,0] = a[:,:,1] = ... = b Parameters ---------- a : (M, N) or (M, N, P) ndarray Target array. b : (M, N) Source array. Notes ----- Color images are stored as an ``(M, N, 3)`` or ``(M, N, 4)`` arrays. N)rrnewaxis)abs r _stackcopyr\sG" v{{AAArz!"!!!!!!r"ct|}|d|d}}t|||g}t|dkr||dtj||}tj||f|ddj}||}|jd||fdd}t|d|dt|d|dt|dkrt|d|d <|S) aBuild the source coordinates for the output of a 2-D image warp. Parameters ---------- coord_map : callable like GeometricTransform.inverse Return input coordinates for given output coordinates. Coordinates are in the shape (P, 2), where P is the number of coordinates and each element is a ``(row, col)`` pair. shape : tuple Shape of output image ``(rows, cols[, bands])``. dtype : np.dtype or string dtype for return value (sane choices: float32 or float64). Returns ------- coords : (ndim, rows, cols[, bands]) array of dtype `dtype` Coordinates for `scipy.ndimage.map_coordinates`, that will yield an image of shape (orows, ocols, bands) by drawing from source points according to the `coord_transform_fn`. Notes ----- This is a lower-level routine that produces the source coordinates for 2-D images used by `warp()`. It is provided separately from `warp` to give additional flexibility to users who would like, for example, to re-use a particular coordinate mapping, to use specific dtypes at various points along the the image-warping process, or to implement different post-processing logic than `warp` performs after the call to `ndi.map_coordinates`. Examples -------- Produce a coordinate map that shifts an image up and to the right: >>> from skimage import data >>> from scipy.ndimage import map_coordinates >>> >>> def shift_up10_left20(xy): ... return xy - np.array([-20, 10])[None, :] >>> >>> image = data.astronaut().astype(np.float32) >>> coords = warp_coords(shift_up10_left20, image.shape) >>> warped_image = map_coordinates(image, coords) rrrr r3r)r.)r.)r .) r rappendremptyindicesrr{swapaxesrrange) coord_maprr3rgrh coords_shapecoords tf_coordss r warp_coordsrss9b   Eq58$DJJd+L 5zzQE!H%%% Xl% 0 0 0F D$>q!DDIvf~y0111vf~y0111 5zzQuQxv Mr"cX|r&tj|}tj|r$tj}tj}||}ntj}tj}||}|dko2||cxko|knc o!|||cxko||knc} | r:|j|}t||}t ||}tj|tj|}}tj ||||dSdS)aClip output image to range of values of input image. Note that this function modifies the values of `output_image` in-place and it is only modified if ``clip=True``. Parameters ---------- input_image : ndarray Input image. output_image : ndarray Output image, which is modified in-place. Other parameters ---------------- mode : {'constant', 'edge', 'symmetric', 'reflect', 'wrap'} Points outside the boundaries of the input are filled according to the given mode. Modes match the behaviour of `numpy.pad`. cval : float Used in conjunction with mode 'constant', the value outside the image boundaries. clip : bool Whether to clip the output to the range of values of the input image. This is enabled by default, since higher order interpolation may produce values outside the given input range. rX)rMN) rr`isnannanminnanmaxrar3typerRrD) input_image output_imager-r,rDmin_valmin_funcmax_funcmax_val preserve_cvals r rCrCsb6 B&%% 8G   yHyHh{++GGvHvH(;'' J  It....w..... I&&$HHHH((<2H2HHHHH   )$))$//D'4((G'4((G:g.. 70C0C  gwLAAAAAA;BBr"c i|jdkrtd|jt|j|}|dkrDt ||}|jt jkr|t j }t j |j} || }nt|}d} |dkrtd|dvrtsqd} tt jrjdkr} nwttrj} nZt#drJjd kr?t'tvr)t jjj} | | |j} |jt j krd nd } |jdkrt1| || |||| } nt|jd krig} t3|jdD]8}| t1| |d|f| |||| 9t j| } | tt jrjdkrt9tt jr}nx|jdks |jd krtdfd}t;| d kr*t;|dkr|d|d| df}t=||}|dk}t?|}tAj!||||||} tE|| |||| S)aWarp an image according to a given coordinate transformation. Parameters ---------- image : ndarray Input image. inverse_map : transformation object, callable ``cr = f(cr, **kwargs)``, or ndarray Inverse coordinate map, which transforms coordinates in the output images into their corresponding coordinates in the input image. There are a number of different options to define this map, depending on the dimensionality of the input image. A 2-D image can have 2 dimensions for gray-scale images, or 3 dimensions with color information. - For 2-D images, you can directly pass a transformation object, e.g. `skimage.transform.SimilarityTransform`, or its inverse. - For 2-D images, you can pass a ``(3, 3)`` homogeneous transformation matrix, e.g. `skimage.transform.SimilarityTransform.params`. - For 2-D images, a function that transforms a ``(M, 2)`` array of ``(col, row)`` coordinates in the output image to their corresponding coordinates in the input image. Extra parameters to the function can be specified through `map_args`. - For N-D images, you can directly pass an array of coordinates. The first dimension specifies the coordinates in the input image, while the subsequent dimensions determine the position in the output image. E.g. in case of 2-D images, you need to pass an array of shape ``(2, rows, cols)``, where `rows` and `cols` determine the shape of the output image, and the first dimension contains the ``(row, col)`` coordinate in the input image. See `scipy.ndimage.map_coordinates` for further documentation. Note, that a ``(3, 3)`` matrix is interpreted as a homogeneous transformation matrix, so you cannot interpolate values from a 3-D input, if the output is of shape ``(3,)``. See example section for usage. map_args : dict, optional Keyword arguments passed to `inverse_map`. output_shape : tuple (rows, cols), optional Shape of the output image generated. By default the shape of the input image is preserved. Note that, even for multi-band images, only rows and columns need to be specified. order : int, optional The order of interpolation. The order has to be in the range 0-5: - 0: Nearest-neighbor - 1: Bi-linear (default) - 2: Bi-quadratic - 3: Bi-cubic - 4: Bi-quartic - 5: Bi-quintic Default is 0 if image.dtype is bool and 1 otherwise. mode : {'constant', 'edge', 'symmetric', 'reflect', 'wrap'}, optional Points outside the boundaries of the input are filled according to the given mode. Modes match the behaviour of `numpy.pad`. cval : float, optional Used in conjunction with mode 'constant', the value outside the image boundaries. clip : bool, optional Whether to clip the output to the range of values of the input image. This is enabled by default, since higher order interpolation may produce values outside the given input range. preserve_range : bool, optional Whether to keep the original range of values. Otherwise, the input image is converted according to the conventions of `img_as_float`. Also see https://scikit-image.org/docs/dev/user_guide/data_types.html Returns ------- warped : double ndarray The warped input image. Notes ----- - The input image is converted to a `double` image. - In case of a `SimilarityTransform`, `AffineTransform` and `ProjectiveTransform` and `order` in [0, 3] this function uses the underlying transformation matrix to warp the image with a much faster routine. Examples -------- >>> from skimage.transform import warp >>> from skimage import data >>> image = data.camera() The following image warps are all equal but differ substantially in execution time. The image is shifted to the bottom. Use a geometric transform to warp an image (fast): >>> from skimage.transform import SimilarityTransform >>> tform = SimilarityTransform(translation=(0, -10)) >>> warped = warp(image, tform) Use a callable (slow): >>> def shift_down(xy): ... xy[:, 1] -= 10 ... return xy >>> warped = warp(image, shift_down) Use a transformation matrix to warp an image (fast): >>> matrix = np.array([[1, 0, 0], [0, 1, -10], [0, 0, 1]]) >>> warped = warp(image, matrix) >>> from skimage.transform import ProjectiveTransform >>> warped = warp(image, ProjectiveTransform(matrix=matrix)) You can also use the inverse of a geometric transformation (fast): >>> warped = warp(image, tform.inverse) For N-D images you can pass a coordinate array, that specifies the coordinates in the input image for every element in the output image. E.g. if you want to rescale a 3-D cube, you can do: >>> cube_shape = np.array([30, 30, 30]) >>> rng = np.random.default_rng() >>> cube = rng.random(cube_shape) Setup the coordinate array, that defines the scaling: >>> scale = 0.1 >>> output_shape = (scale * cube_shape).astype(int) >>> coords0, coords1, coords2 = np.mgrid[:output_shape[0], ... :output_shape[1], :output_shape[2]] >>> coords = np.array([coords0, coords1, coords2]) Assume that the cube contains spatial data, where the first array element center is at coordinate (0.5, 0.5, 0.5) in real space, i.e. we have to account for this extra offset when scaling the image: >>> coords = (coords + 0.5) / scale - 0.5 >>> warped = warp(cube, coords) Nrz'Cannot warp empty image with dimensionsr aZBi-quadratic interpolation behavior has changed due to a bug in the implementation of scikit-image. The new version now serves as a wrapper around SciPy's interpolation functions, which itself is not verified to be a correct implementation. Until skimage's implementation is fixed, we recommend to use bi-linear or bi-cubic interpolation instead.)rr)rr__name__r_ float32_t float64_t)rr1r-r,r.)matrixz\Only 2-D images (grayscale or color) are supported, when providing a callable `inverse_map`.c|iSr&r')args inverse_maprs r rzwarp..coord_maps"{D5H555r"r) prefilterr-r1r,)#sizerrrr3rrr4r5r6r]r r isinstancendarrayHOMOGRAPHY_TRANSFORMSrchasattrrr linalginv__self__rrrrdstackrrrrr@map_coordinatesrC)rrrrr1r-r,rDrErwarpedrctypedimsdimrrrrJs `` r rdrdsn zQBEKPPP )%+u = =E qyy 77 ;"* $ $LL,,E(5;''K" "<00 F zz  B    x k2: . . @;3D3N3N FF  %: ; ; @ 'FF K , , @$ 11&{337LLLY]];#7#>??F  ]]5;//F#(;"*#<#>> from skimage import data >>> from skimage.transform import warp_polar >>> image = data.checkerboard() >>> warped = warp_polar(image) Perform a log-polar warp on a grayscale image: >>> warped = warp_polar(image, scaling='log') Perform a log-polar warp on a grayscale image while specifying center, radius, and output shape: >>> warped = warp_polar(image, (100,100), radius=100, ... output_shape=image.shape, scaling='log') Perform a log-polar warp on a color image: >>> image = data.astronaut() >>> warped = warp_polar(image, scaling='log', channel_axis=-1) Nr z@Input array must be 2-dimensional when `channel_axis=None`, got rzHInput array must be 3-dimensional when `channel_axis` is specified, got rZihrrrr}z*Scaling value must be in {'linear', 'log'})rrrf)rr)rrrr]rr|intceilr rr}rpird)rrfrrrrOkwargsrUwhheightwidthrmap_funcrrrs r warp_polarrSsF t+L zQ| 5(-  5 5    zQ< =05  = =   ~(5;''+a/36 ~x $$RaR(1,1A1%%BGFOO$$ "<00 aQ(6>( E  26&>>)%EFFFBE "G#VLLI  x"+,JPF Mr"c|r7tjd||dz|}tjd||dz|}n|dz |dz }}tjtjd|dz ||ddd|f}|dkr tj}nd|z|z }tjtj|||z ||ddd|f}tj|ddtjf|tjddf} tj|dd tjf|tjdd f} tj| | z d} | | dd z} | S) aCreate a 2D weight matrix for resizing with the local mean. Parameters ---------- old_size: int Old size. new_size: int New size. grid_mode : bool Whether to use grid data model of pixel/voxel model for average weights computation. dtype: dtype Output array data type. Returns ------- weights: (new_size, old_size) array Rows sum to 1. rr)numr3rZrrX)constant_valuesNrT)axiskeepdims)rlinspacepadinfminimumrr=sum) old_sizenew_sizer2r3 old_breaks new_breaksoldnewvalupperlowerweightss r _local_mean_weightsrs* [H(Q,eLLL [H(Q,eLLL a<ASV KS3Y5 9 9 9 H    !88&CC)c/CV KS3Y5 9 9 9 H     Jz!""bj.1:bj!""n3M N NE Jz#2#rz/2Jrz3B34O P PEj**G w{{D{111G Nr"c|||j ks ||jkrtdtj||d}|jd}t |}||jdz kr||fz}nY||jkr?|||krtd||jz}|d|||dz|fz}ntd|}nt ||\}}t||}|j}tt|j|D]Q\} \} } | | krt| | ||} tj || | gdgg} tj| d| }R|tj|d|}|S)a Resize an array with the local mean / bilinear scaling. Parameters ---------- image : ndarray Input image. If this is a multichannel image, the axis corresponding to channels should be specified using `channel_axis`. output_shape : iterable Size of the generated output image. When `channel_axis` is not None, the `channel_axis` should either be omitted from `output_shape` or the ``output_shape[channel_axis]`` must match ``image.shape[channel_axis]``. If the length of `output_shape` exceeds image.ndim, additional singleton dimensions will be appended to the input ``image`` as needed. grid_mode : bool, optional Defines ``image`` pixels position: if True, pixels are assumed to be at grid intersections, otherwise at cell centers. As a consequence, for example, a 1d signal of length 5 is considered to have length 4 when `grid_mode` is False, but length 5 when `grid_mode` is True. See the following visual illustration: .. code-block:: text | pixel 1 | pixel 2 | pixel 3 | pixel 4 | pixel 5 | |<-------------------------------------->| vs. |<----------------------------------------------->| The starting point of the arrow in the diagram above corresponds to coordinate location 0 in each mode. preserve_range : bool, optional Whether to keep the original range of values. Otherwise, the input image is converted according to the conventions of `img_as_float`. Also see https://scikit-image.org/docs/dev/user_guide/data_types.html Returns ------- resized : ndarray Resized version of the input. See Also -------- resize, downscale_local_mean Notes ----- This method is sometimes referred to as "area-based" interpolation or "pixel mixing" interpolation [1]_. When `grid_mode` is True, it is equivalent to using OpenCV's resize with `INTER_AREA` interpolation mode. It is commonly used for image downsizing. If the downsizing factors are integers, then `downscale_local_mean` should be preferred instead. References ---------- .. [1] http://entropymine.com/imageworsener/pixelmixing/ Examples -------- >>> from skimage import data >>> from skimage.transform import resize_local_mean >>> image = data.camera() >>> resize_local_mean(image, (100, 100)).shape (100, 100) Nzinvalid channel_axisrrzWCannot reshape along the channel_axis. Use channel_axis=None to reshape along all axes.zZlen(output_shape) must be image.ndim or (image.ndim - 1) when a channel_axis is specified.) rrrmoveaxisrrr!rr3 enumerater;r tensordot)rrr2rErOncrresizedr3rrrrproducts r resize_local_meanrsJ 5:+ % %)C)C344 4 E<44 [_,'' %*q. ( ('2%/LL EJ & &L)R// C (%*4L]l]+l<==.IIREQ L4  ?| T Tw77G ME&/EK0N0N&O&O11""x x   %h)UKK,w4&2$@@+gr400+gr<88 Nr")Nr#rTFNN)FNNrXrTF)rT) NrrrNNr#rTF)NNNrXrTFr&)TF)'numpyrscipyrr@ _geometricrrr _warps_cyrmeasurer _shared.utilsr r r rrrrrr!rNrWrurxrrrfloat64rrCrdrrrrrr'r"r rs QQQQQQQQQQ!!!!!!""""""-o?RS///j   JJJJZ   rrrrrrp     {{{{|,7,7,7,7^&       QQQQh.)+ KKKK\8B8B8B|   VVVVr>> o  oooood111j9>rOSrrrrrrrr"