a
    Re                     @   s   d dl Z d dlZd dlmZmZmZmZmZmZ d dl	Z	d dl	m
Z d dlmZ ddlmZ dZdddd	d
ddZG dd dZdd ZdS )    N)assert_assert_equalassert_array_equalassert_array_almost_equalassert_allclosesuppress_warnings)raises   )types-q=reflectZ	symmetricwrapZedgeconstant)mirrorr   grid-mirror	grid-wrapnearestgrid-constantc                   @   s  e Zd Zejddg dfdg dfdg dfdg d	fd
g dfdg dfdg dfgdd Zejddg dfdg dfdg dfdg dfd
g dfdg dfdg dfgdd Zejdg dejdeddd Z	ejdeddejd e
d!d" Zejdeddejd e
d#d$ Zejdeddejd e
d%d& Zejdeddejd e
d'd( Zejdeddejd e
d)d* Zejded+dd,d- Zejded+dd.d/ Zejded+dd0d1 Zejded+dd2d3 Zejded+dejd ejejgd4d5 Zejded+dd6d7 Zejded+dd8d9 Zejded+dd:d; Zejded+dd<d= Zejded+dd>d? Zejded+dd@dA Zejded+ddBdC Zejded+ddDdE Zejded+ddFdG Z ejded+ddHdI Z!ejded+ddJdK Z"ejded+ddLdM Z#ejded+ddNdO Z$ejded+ddPdQ Z%ejded+ddRdS Z&ejded+ddTdU Z'dVdW Z(ejdg dXejdeddYdZ Z)d[d\ Z*d]d^ Z+ejded+dejd ejejgd_d` Z,ejded+ddadb Z-dcdd Z.dedf Z/dgdh Z0ejj1die2j3v pe4d+j5djk dkdldmdn Z6ejded+ddodp Z7ejded+ddqdr Z8ejded+ddsdt Z9ejded+ddudv Z:ejded+dejd ejejgdwdx Z;ejded+ddydz Z<ejded+dd{d| Z=ejded+dd}d~ Z>ejded+ddd Z?ejded+ddd Z@ejded+ddd ZAejded+ddd ZBejded+ddd ZCejded+ddd ZDejded+ddd ZEejded+ddd ZFejded+ddd ZGejded+ddd ZHejded+ddd ZIejded+ddd ZJejded+ddd ZKejded+ddd ZLejded+ddd ZMejded+ddd ZNejded+ddd ZOejded+ddd ZPdd ZQdd ZRdd ZSdd ZTdd ZUejdg dejded+ddd ZVejded+ddd ZWejded+ddd ZXejded+ddd ZYejded+ddd ZZejded+ddd Z[ejded+dejd ejejgdd Z\ejded+dejdddgejd ejejgdd Z]ejded+ddd Z^ejded+ddd Z_ejded+dddĄ Z`ejded+dddƄ Zaejdg dejded+dddȄ Zbejdg dejded+dddʄ Zcdd̄ Zdejded+ddd΄ Zeejded+dejdddgddӄ Zfejdg dXejdedddՄ Zgejded+dddׄ Zhddل Ziddۄ Zjejded+dejd ejejgdd݄ Zkdd߄ Zldd Zmdd Znejdg dejdg ddd Zoejdddgejdg dejdg ddd Zpejdddgdd Zqejded+ddd Zrejded+ddd Zsejded+dejd ejejgdd Ztejded+ddd Zuejded+ddd Zvejded+ddd Zwejded+ddd Zxejded+ddd  Zydd Zzdd Z{dd Z|dS (  TestNdimageInterpolationzmode, expected_valuer   )      ?      @      @   r   r   r   r   )r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   )r   r   r   r   r   r   r   c              
   C   s6   dd }t g d}t|tj||d|ddd d S )Nc                 S   s   | d d fS )Nr         ? xr   r   Z/var/www/sunrise/env/lib/python3.9/site-packages/scipy/ndimage/tests/test_interpolation.pyshift&   s    z7TestNdimageInterpolation.test_boundaries.<locals>.shift)r	               @r   )   r	   )cvalmodeoutput_shapeordernumpyarrayr   ndimagegeometric_transformselfr%   Zexpected_valuer   datar   r   r   test_boundaries   s    z(TestNdimageInterpolation.test_boundaries)r	   r	   r    r!   )r!   r	   r    r!   )r   r	   r    r!   )r    r	   r    r!   )r   r	   r    r!   c              	   C   s4   dd }t g d}t|tj||d|dd d S )Nc                 S   s   | d d fS )Nr   g?r   r   r   r   r   r   :   s    z8TestNdimageInterpolation.test_boundaries2.<locals>.shiftr	   r    r!   r   r   r   )r$   r%   r&   r(   r-   r   r   r   test_boundaries2/   s    z)TestNdimageInterpolation.test_boundaries2r%   )r   r   r   r   r   r   r'      c                 C   s   t jddtd}t jdddd}tj||g||d}d	}t|}t j|||d
}tj||| g||d}	|dkrvdnd}
t	||	d|
d dS )z$Tests based on examples from gh-2640ir#   dtypei   i  )numr'   r%       r%   r   gh㈵>r   Hz>)rtolatolN)
r)   arangefloatZlinspacer+   map_coordinatesndimage_to_numpy_modegetpadr   )r.   r%   r'   r/   r   ynpadpad_modepaddedexpectedr>   r   r   r   test_boundary_spline_accuracyC   s    
z6TestNdimageInterpolation.test_boundary_spline_accuracyr    r6   c                 C   s(   t g |}tj||d}t|d d S )Nr'   r	   r)   onesr+   spline_filterr   r.   r6   r'   r/   outr   r   r   test_spline01W   s    z&TestNdimageInterpolation.test_spline01c                 C   s,   t dg|}tj||d}t|dg d S Nr	   rK   )r)   r*   r+   rN   r   rO   r   r   r   test_spline02^   s    z&TestNdimageInterpolation.test_spline02c                 C   s*   t g |}tj|||d}t|d d S )Noutputr	   rL   rO   r   r   r   test_spline03e   s    z&TestNdimageInterpolation.test_spline03c                 C   s,   t dg|}t||}t|g d d S )Nr   r	   r	   r	   r	   rL   rO   r   r   r   test_spline04l   s    z&TestNdimageInterpolation.test_spline04c                 C   sD   t ddg|}tj||d}t|g dg dg dg dg d S )Nr   rK   rW   rL   rO   r   r   r   test_spline05s   s    
z&TestNdimageInterpolation.test_spline05r   c                 C   s8   t dg}dd }tj|||j|d}t|dg d S )Nr	   c                 S   s   | S Nr   r   r   r   r   mapping   s    zDTestNdimageInterpolation.test_geometric_transform01.<locals>.mappingrK   r)   r*   r+   r,   shaper   r.   r'   r/   r[   rP   r   r   r   test_geometric_transform01}   s    z3TestNdimageInterpolation.test_geometric_transform01c                 C   s:   t dg}dd }tj|||j|d}t|g d d S )Nr   c                 S   s   | S rZ   r   r   r   r   r   r[      s    zDTestNdimageInterpolation.test_geometric_transform02.<locals>.mappingrK   rW   r)   rM   r+   r,   r]   r   r^   r   r   r   test_geometric_transform02   s    z3TestNdimageInterpolation.test_geometric_transform02c                 C   s:   t dg}dd }tj|||j|d}t|g d d S )Nr   c                 S   s   | d d fS Nr   r	   r   r   r   r   r   r[      s    zDTestNdimageInterpolation.test_geometric_transform03.<locals>.mappingrK   r   r	   r	   r	   r`   r^   r   r   r   test_geometric_transform03   s    z3TestNdimageInterpolation.test_geometric_transform03c                 C   s<   t g d}dd }tj|||j|d}t|g d d S )Nr   r	   r!   r    c                 S   s   | d d fS rb   r   r   r   r   r   r[      s    zDTestNdimageInterpolation.test_geometric_transform04.<locals>.mappingrK   r   r   r	   r!   r\   r^   r   r   r   test_geometric_transform04   s    z3TestNdimageInterpolation.test_geometric_transform04c                 C   s   t jg dg dg dg|d}t jg dg dg dg|d}|jjdkrd|d| 8 }|d| 8 }dd }tj|||j|d}t|| d S )	NrW   r5   rc   c              ?c                 S   s   | d | d d fS rb   r   r   r   r   r   r[      s    zDTestNdimageInterpolation.test_geometric_transform05.<locals>.mappingrK   )r)   r*   r6   kindr+   r,   r]   r   )r.   r'   r6   r/   rI   r[   rP   r   r   r   test_geometric_transform05   s(    

z3TestNdimageInterpolation.test_geometric_transform05c                 C   sX   t g dg dg dg}dd }tj|||j|d}t|g dg dg d	g d S )
Nre   r#   r4         r!   rn   r!   r4   c                 S   s   | d | d d fS rb   r   r   r   r   r   r[      s    zDTestNdimageInterpolation.test_geometric_transform06.<locals>.mappingrK   rf   r   r#   r4   rm   r   r!   rn   r!   r\   r^   r   r   r   test_geometric_transform06   s    

z3TestNdimageInterpolation.test_geometric_transform06c                 C   sX   t g dg dg dg}dd }tj|||j|d}t|g dg dg dg d S )Nre   rl   ro   c                 S   s   | d d | d fS rb   r   r   r   r   r   r[      s    zDTestNdimageInterpolation.test_geometric_transform07.<locals>.mappingrK   r   r   r   r   r\   r^   r   r   r   test_geometric_transform07   s    

z3TestNdimageInterpolation.test_geometric_transform07c                 C   sX   t g dg dg dg}dd }tj|||j|d}t|g dg dg d	g d S )
Nre   rl   ro   c                 S   s   | d d | d d fS rb   r   r   r   r   r   r[      s    zDTestNdimageInterpolation.test_geometric_transform08.<locals>.mappingrK   rs   rf   rp   r\   r^   r   r   r   test_geometric_transform08   s    

z3TestNdimageInterpolation.test_geometric_transform08c                 C   sv   t g dg dg dg}dd }|dkr<tj||d}n|}tj|||j|dd	}t|g d
g dg dg d S )Nre   rl   ro   c                 S   s   | d d | d d fS rb   r   r   r   r   r   r[      s    zDTestNdimageInterpolation.test_geometric_transform10.<locals>.mappingr	   rK   Fr'   	prefilterrs   rf   rp   )r)   r*   r+   rN   r,   r]   r   )r.   r'   r/   r[   filteredrP   r   r   r   test_geometric_transform10   s    

z3TestNdimageInterpolation.test_geometric_transform10c                 C   s>   t dgt j}dd }tj||dg|d}t|g d d S )Nr    c                 S   s   | d d fS Nr   r    r   r   r   r   r   r[     s    zDTestNdimageInterpolation.test_geometric_transform13.<locals>.mappingr   rK   rW   )r)   rM   float64r+   r,   r   r^   r   r   r   test_geometric_transform13   s    z3TestNdimageInterpolation.test_geometric_transform13c                 C   s6   g d}dd }t j||dg|d}t|g d d S )Nr	   rn   r    r4   r!   r#   r   r   c                 S   s   d| d  fS )Nr    r   r   r   r   r   r   r[     s    zDTestNdimageInterpolation.test_geometric_transform14.<locals>.mappingr   rK   r1   r+   r,   r   r^   r   r   r   test_geometric_transform14  s    z3TestNdimageInterpolation.test_geometric_transform14c                 C   s@   g d}dd }t j||dg|d}t|d d d g d d S )Nr1   c                 S   s   | d d fS rz   r   r   r   r   r   r[     s    zDTestNdimageInterpolation.test_geometric_transform15.<locals>.mappingrm   rK   r    r~   r^   r   r   r   test_geometric_transform15  s    z3TestNdimageInterpolation.test_geometric_transform15c                 C   sP   g dg dg dg}dd }t j||d|d}t|dd	gd
dgddgg d S )Nr1   rn   r4   r#   rm   g      "@
         c                 S   s   | d | d d fS Nr   r	   r    r   r   r   r   r   r[   !  s    zDTestNdimageInterpolation.test_geometric_transform16.<locals>.mappingr!   r    rK   r	   r!   rn   r#   	   r   r~   r^   r   r   r   test_geometric_transform16  s    
z3TestNdimageInterpolation.test_geometric_transform16c                 C   sD   g dg dg dg}dd }t j||d|d}t|g dg d S )Nr1   r   r   r   r   r   c                 S   s   | d d | d fS Nr   r    r	   r   r   r   r   r   r[   .  s    zDTestNdimageInterpolation.test_geometric_transform17.<locals>.mappingr	   r   rK   r~   r^   r   r   r   test_geometric_transform17(  s    
z3TestNdimageInterpolation.test_geometric_transform17c                 C   sD   g dg dg dg}dd }t j||d|d}t|dd	gg d S )
Nr1   r   r   c                 S   s   | d d | d d fS r   r   r   r   r   r   r[   ;  s    zDTestNdimageInterpolation.test_geometric_transform18.<locals>.mappingr	   r    rK   r	   r!   r~   r^   r   r   r   test_geometric_transform185  s    
z3TestNdimageInterpolation.test_geometric_transform18c                 C   sL   g dg dg dg}dd }t j||d|d}t|dd d d	f | d S )
Nr1   r   r   c                 S   s   | d | d d fS r   r   r   r   r   r   r[   H  s    zDTestNdimageInterpolation.test_geometric_transform19.<locals>.mappingr!   rm   rK   .r    r~   r^   r   r   r   test_geometric_transform19B  s    
z3TestNdimageInterpolation.test_geometric_transform19c                 C   sL   g dg dg dg}dd }t j||d|d}t|d d dd	f | d S )
Nr1   r   r   c                 S   s   | d d | d fS r   r   r   r   r   r   r[   U  s    zDTestNdimageInterpolation.test_geometric_transform20.<locals>.mappingr4   r   rK   r    .r~   r^   r   r   r   test_geometric_transform20O  s    
z3TestNdimageInterpolation.test_geometric_transform20c                 C   sR   g dg dg dg}dd }t j||d|d}t|d d dd d df | d S )	Nr1   r   r   c                 S   s   | d d | d d fS r   r   r   r   r   r   r[   b  s    zDTestNdimageInterpolation.test_geometric_transform21.<locals>.mappingr4   rm   rK   r    r~   r^   r   r   r   test_geometric_transform21\  s    
z3TestNdimageInterpolation.test_geometric_transform21c                 C   sb   t g dg dg dgt j}dd }dd }tj||d|d	}tj||d
|d	}t|| d S )Nr1   r   r   c                 S   s   | d d | d d fS r   r   r   r   r   r   mapping1o  s    zETestNdimageInterpolation.test_geometric_transform22.<locals>.mapping1c                 S   s   | d d | d d fS r   r   r   r   r   r   mapping2r  s    zETestNdimageInterpolation.test_geometric_transform22.<locals>.mapping2r   rK   r!   r   )r)   r*   r{   r+   r,   r   )r.   r'   r/   r   r   rP   r   r   r   test_geometric_transform22i  s    
z3TestNdimageInterpolation.test_geometric_transform22c                 C   sN   g dg dg dg}dd }t j||d|d}|tj}t|dd	g d S )
Nr1   r   r   c                 S   s   d| d d fS )Nr	   r   r    r   r   r   r   r   r[     s    zDTestNdimageInterpolation.test_geometric_transform23.<locals>.mappingr    rK   rn   r#   )r+   r,   astyper)   int32r   r^   r   r   r   test_geometric_transform23{  s    z3TestNdimageInterpolation.test_geometric_transform23c                 C   sJ   g dg dg dg}dd }t j||d|ddd	id
}t|ddg d S )Nr1   r   r   c                 S   s   || d | fS )Nr   r   )r   abr   r   r   r[     s    zDTestNdimageInterpolation.test_geometric_transform24.<locals>.mappingr   )r	   r   r    )r'   Zextra_argumentsZextra_keywordsrn   r#   r~   r^   r   r   r   test_geometric_transform24  s    
z3TestNdimageInterpolation.test_geometric_transform24c                 C   sT   t jg dg dgtd}dd }t g dg dg}ttj||dd	d
| d S )Nr	   r    r!   r   rn   r4   r5   c                 S   s   | d d | d d fS )Nr   r   r	   r   r   r   r   r   r[     s    zWTestNdimageInterpolation.test_geometric_transform_grid_constant_order1.<locals>.mappingg      ?g      ?      ?r   g      @r"   r   r	   r%   r'   )r)   r*   r@   r   r+   r,   )r.   r   r[   expected_resultr   r   r   -test_geometric_transform_grid_constant_order1  s    


zFTestNdimageInterpolation.test_geometric_transform_grid_constant_order1)r   r   r   r   r   c           
      C   s   t jdtddd}dd }d}t|}t j|||d}tt|| g|j	 }t
j||||d| }	tt
j||||d|	d	d
 d S )N   r5   r   c                 S   s   | d d | d d fS )Nr   皙?r	   gffffff@r   r   r   r   r   r[     s    zLTestNdimageInterpolation.test_geometric_transform_vs_padded.<locals>.mapping   r;   r   r<   r=   )r)   r?   r@   reshaperB   rC   rD   tupleslicendimr+   r,   r   )
r.   r'   r%   r   r[   rF   rG   xpcenter_slicer   r   r   r   "test_geometric_transform_vs_padded  s$    

z;TestNdimageInterpolation.test_geometric_transform_vs_paddedc                 C   sz   t dg}dd }|j|j t |t ||j fD ]4}tj|||j|d}|d u rd|n|}t	|dg q@d S )Nr	   c                 S   s   | S rZ   r   r   r   r   r   r[     s    zcTestNdimageInterpolation.test_geometric_transform_endianness_with_output_parameter.<locals>.mappingrT   )
r)   r*   r6   newbyteorder
empty_liker   r+   r,   r]   r   )r.   r/   r[   rP   returnedresultr   r   r   9test_geometric_transform_endianness_with_output_parameter  s    zRTestNdimageInterpolation.test_geometric_transform_endianness_with_output_parameterc                 C   sH   t dg}dd }tj||dd}t|jt du  t|dg d S )Nr	   c                 S   s   | S rZ   r   r   r   r   r   r[     s    zUTestNdimageInterpolation.test_geometric_transform_with_string_output.<locals>.mappingfrT   )r)   r*   r+   r,   r   r6   r   )r.   r/   r[   rP   r   r   r   +test_geometric_transform_with_string_output  s
    zDTestNdimageInterpolation.test_geometric_transform_with_string_outputc                 C   s   t g dg dg dg}t g dg dg dg}|jjdkr\|d|  }|d|  }t |j}|d	8 }tj|||d
}t|| d S )Nre   rl   ro   rs   rf   rp   rh   ri   r	   rK   )	r)   r*   r6   rj   indicesr]   r+   rA   r   )r.   r'   r6   r/   rI   idxrP   r   r   r   test_map_coordinates01  s    

z/TestNdimageInterpolation.test_map_coordinates01c                 C   sb   t g dg dg dg}t |jt j}|d8 }tj|d|d}tj|||d}t|| d S )Nre   rl   ro   r   rK   )	r)   r*   r   r]   r{   r+   r   rA   r   )r.   r'   r/   r   out1out2r   r   r   test_map_coordinates02  s    
z/TestNdimageInterpolation.test_map_coordinates02c              	   C   sP  t jg dg dg dgdd}t |jd }t||}t|g dg dg d	g t|t|d
 t |d d d jd }t|d d d |}t|g dg dg t|t|d d d d
 t |d d d d df jd }t|d d d d df |}t|ddgddgddgg t|t|d d d d df d
 d S )Nre   rl   ro   FrK   r	   rs   rf   rp   r	   r	   r    r   r   r#   )r)   r*   r   r]   r+   rA   r   r   r.   r/   r   rP   r   r   r   test_map_coordinates03  s.    


"z/TestNdimageInterpolation.test_map_coordinates03c                 C   s   t ddgddgg}t ddgddgg}t |j}|d8 }|j|j t |t ||j fD ].}tj	|||d}|d u r|n|}t
|| qld S )Nr	   r    r#   r4   r   rT   )r)   r*   r   r]   r6   r   r   r   r+   rA   r   )r.   r/   rI   r   rP   r   r   r   r   r   5test_map_coordinates_endianness_with_output_parameter  s    zNTestNdimageInterpolation.test_map_coordinates_endianness_with_output_parameterc                 C   sP   t dgg}t |j}tj||dd}t|jt du  t|dgg d S Nr	   r   rT   )	r)   r*   r   r]   r+   rA   r   r6   r   r   r   r   r   'test_map_coordinates_with_string_output  s
    z@TestNdimageInterpolation.test_map_coordinates_with_string_outputwin32rm   z2do not run on 32 bit or windows (no sparse memory))reasonc              
   C   s   z`d}t j|d t jd||}d||d d |d d f< tj||d g|d ggdd W n0 ty } ztd	|W Y d }~n
d }~0 0 d S )
Ni0u  r    r5   r   r!   r   r	   rK   zNot enough memory available)	r)   emptyfloat32r   r+   rA   MemoryErrorpytestskip)r.   nr   er   r   r   test_map_coordinates_large_data$  s    $z8TestNdimageInterpolation.test_map_coordinates_large_datac                 C   s0   t dg}tj|dgg|d}t|dg d S rR   r)   r*   r+   affine_transformr   r.   r'   r/   rP   r   r   r   test_affine_transform012  s    z0TestNdimageInterpolation.test_affine_transform01c                 C   s2   t dg}tj|dgg|d}t|g d d S )Nr   r	   rK   rW   r)   rM   r+   r   r   r   r   r   r   test_affine_transform028  s    z0TestNdimageInterpolation.test_affine_transform02c                 C   s4   t dg}tj|dggd|d}t|g d d S )Nr   r	   r   rK   rc   r   r   r   r   r   test_affine_transform03>  s    z0TestNdimageInterpolation.test_affine_transform03c                 C   s6   t g d}tj|dggd|d}t|g d d S )Nre   r	   r   rK   rf   r   r   r   r   r   test_affine_transform04D  s    z0TestNdimageInterpolation.test_affine_transform04c                 C   s   t jg dg dg dg|d}t jg dg dg dg|d}|jjdkrd|d| 8 }|d| 8 }tj|ddgddggddg|d	}t|| d S )
NrW   r5   rc   rh   ri   r	   r   r   rK   )r)   r*   r6   rj   r+   r   r   r.   r'   r6   r/   rI   rP   r   r   r   test_affine_transform05J  s&    

z0TestNdimageInterpolation.test_affine_transform05c                 C   s^   t g dg dg dg}tj|ddgddggddg|d}t|g dg d	g d
g d S )Nre   rl   ro   r	   r   r   rK   rf   rp   rq   r   r   r   r   r   test_affine_transform06Z  s    

z0TestNdimageInterpolation.test_affine_transform06c                 C   s^   t g dg dg dg}tj|ddgddggddg|d}t|g dg dg dg d S )	Nre   rl   ro   r	   r   r   rK   rs   r   r   r   r   r   test_affine_transform07e  s    

z0TestNdimageInterpolation.test_affine_transform07c                 C   s^   t g dg dg dg}tj|ddgddggddg|d}t|g dg d	g d
g d S )Nre   rl   ro   r	   r   r   rK   rs   rf   rp   r   r   r   r   r   test_affine_transform08p  s    

z0TestNdimageInterpolation.test_affine_transform08c                 C   s|   t g dg dg dg}|dkr4tj||d}n|}tj|ddgddggddg|dd	}t|g d
g dg dg d S )Nre   rl   ro   r	   rK   r   r   Frv   rs   rf   rp   )r)   r*   r+   rN   r   r   r.   r'   r/   rx   rP   r   r   r   test_affine_transform09{  s    

z0TestNdimageInterpolation.test_affine_transform09c                 C   s8   t dgt j}tj|dggd|d}t|g d d S )Nr    r   r2   )r&   r'   r	   r	   r	   r   )r)   rM   r{   r+   r   r   r   r   r   r   test_affine_transform10  s
    z0TestNdimageInterpolation.test_affine_transform10c                 C   s2   g d}t j|dggdd|d}t|g d d S )Nr}   r    r   r2   rK   r1   r+   r   r   r   r   r   r   test_affine_transform11  s    z0TestNdimageInterpolation.test_affine_transform11c                 C   s<   g d}t j|dggdd|d}t|d d d g d d S )Nr1   r   r   rm   rK   r    r   r   r   r   r   test_affine_transform12  s    z0TestNdimageInterpolation.test_affine_transform12c                 C   sV   g dg dg dg}t j|ddgddggdd|d}t|dd	gd
dgddgg d S )Nr1   r   r   r	   r   r    r   rK   r!   rn   r#   r   r   r   r   r   r   r   test_affine_transform13  s    z0TestNdimageInterpolation.test_affine_transform13c                 C   sJ   g dg dg dg}t j|ddgddggdd|d}t|g dg d S )	Nr1   r   r   r    r   r	   r   rK   r   r   r   r   r   test_affine_transform14  s    z0TestNdimageInterpolation.test_affine_transform14c                 C   sJ   g dg dg dg}t j|ddgddggdd|d}t|dd	gg d S )
Nr1   r   r   r    r   r   rK   r	   r!   r   r   r   r   r   test_affine_transform15  s    z0TestNdimageInterpolation.test_affine_transform15c                 C   sR   g dg dg dg}t j|ddgddggdd|d	}t|d
d d df | d S )Nr1   r   r   r	           r   r   r   rK   .r    r   r   r   r   r   test_affine_transform16  s    z0TestNdimageInterpolation.test_affine_transform16c                 C   sR   g dg dg dg}t j|ddgddggdd|d}t|d d d	d
f | d S )Nr1   r   r   r   r   r	   r   rK   r    .r   r   r   r   r   test_affine_transform17  s    z0TestNdimageInterpolation.test_affine_transform17c                 C   sX   g dg dg dg}t j|ddgddggdd|d}t|d d dd d df | d S )	Nr1   r   r   r   r   r   rK   r    r   r   r   r   r   test_affine_transform18  s    z0TestNdimageInterpolation.test_affine_transform18c                 C   sn   t g dg dg dgt j}tj|ddgddggdd|d}tj|ddgddggdd	|d}t|| d S )
Nr1   r   r   r   r   r   rK          @r   )r)   r*   r{   r+   r   r   r   r   r   r   test_affine_transform19  s    
z0TestNdimageInterpolation.test_affine_transform19c                 C   sD   g dg dg dg}t j|dgdggdd|d}t|dd	g d S )
Nr1   r   r   r   r    r   rK   r	   r!   r   r   r   r   r   test_affine_transform20  s    z0TestNdimageInterpolation.test_affine_transform20c                 C   sD   g dg dg dg}t j|dgdggdd|d}t|dd	g d S )
Nr1   r   r   r    r   r   rK   r	   r   r   r   r   r   r   test_affine_transform21  s    z0TestNdimageInterpolation.test_affine_transform21c                 C   s:   t g d}tj|dggdgd|d}t|g d d S )Nre   r    r   )r!   rK   r   r	   r    r   r   r   r   r   test_affine_transform22  s
    z0TestNdimageInterpolation.test_affine_transform22c                 C   sD   t g d}tj|dggdgd|d}t|d d d g d d S )Nre   r   r   r   rK   r    rf   r   r   r   r   r   test_affine_transform23  s
    z0TestNdimageInterpolation.test_affine_transform23c                 C   sx   t g d}t 0}|td tj|dgd|d}W d    n1 sJ0    Y  tj|dggd|d}t|| d S )Nre   @The behavior of affine_transform with a 1-D array .* has changedr    r   rK   r)   r*   r   filterUserWarningr+   r   r   r.   r'   r/   supr   r   r   r   r   test_affine_transform24  s    2z0TestNdimageInterpolation.test_affine_transform24c                 C   sx   t g d}t 0}|td tj|dgd|d}W d    n1 sJ0    Y  tj|dggd|d}t|| d S )Nre   r   r   r   rK   r   r   r   r   r   test_affine_transform25  s    2z0TestNdimageInterpolation.test_affine_transform25c                 C   s   t g dg dg dg}|dkr4tj||d}n|}t d}t d }t ||f}t |g dgf}tj|||	 |d	d
}tj|||d	d
}	tj|||d	d
}
||	|
fD ] }t
|g dg dg dg qd S )Nre   rl   ro   r	   rK   r    r    r	   )r   r   r	   Frv   rs   rf   rp   )r)   r*   r+   rN   eyerM   hstackvstackr   Zravelr   )r.   r'   r/   rx   Ztform_originalZoffset_originaltform_h1tform_h2r   r   Zout3rP   r   r   r   test_affine_transform26  s4    




z0TestNdimageInterpolation.test_affine_transform26c                 C   s`   t g dg dg dg}t t dt d f}t |g dgf}tttj	|| d S )Nre   rl   ro   r    r   )rn   r    r	   )
r)   r*   r   r   rM   r   assert_raises
ValueErrorr+   r   )r.   r/   r  r  r   r   r   test_affine_transform273  s    
z0TestNdimageInterpolation.test_affine_transform27c              	   C   s   t d}t |t ||j |j|j fD ]p}t 0}|td t	j
|ddg|d}W d    n1 sv0    Y  |d u r|n|}t|ddgddgg q6d S )N)r    r    r   r	   rT   )r)   rM   r   r   r6   r   r   r   r   r+   r   r   )r.   r/   rP   r   r   r   r   r   r   9test_affine_transform_1d_endianness_with_output_parameter<  s    
2zRTestNdimageInterpolation.test_affine_transform_1d_endianness_with_output_parameterc                 C   sr   t dg}|j|j t |t ||j fD ]4}tj|dgg|d}|d u r\|n|}t|dg q8d S )Nr	   rT   )	r)   r*   r6   r   r   r   r+   r   r   )r.   r/   rP   r   r   r   r   r   >test_affine_transform_multi_d_endianness_with_output_parameterK  s    zWTestNdimageInterpolation.test_affine_transform_multi_d_endianness_with_output_parameterc                 C   s   t jdt jd}t d}tj|dgg|d t|d d | tt	& tj|dgg|dd W d    n1 sv0    Y  d S )Nrm   r5   )   r	   rT   )r   )rU   r&   )
r)   r?   r{   rM   r+   r   r   r   r   RuntimeErrorr.   r/   rP   r   r   r   "test_affine_transform_output_shapeV  s    
z;TestNdimageInterpolation.test_affine_transform_output_shapec                 C   sD   t dg}tj|dggdd}t|jt du  t|dg d S r   )r)   r*   r+   r   r   r6   r   r  r   r   r   (test_affine_transform_with_string_outputc  s    zATestNdimageInterpolation.test_affine_transform_with_string_outputr   ))r	   r   r   r	   )r   r	   )r!   )r    r#   c                 C   sr   t ddgddgg}t d}t d|d dd df< ||d d df< ttj||d|dt j||dd	 d S )
Nr   r	   r    r!   r    r!   r   r   r  Zaxis)r)   r*   zerosr   r   r+   r   roll)r.   r   r'   r   affiner   r   r   )test_affine_transform_shift_via_grid_wrapi  s    

zBTestNdimageInterpolation.test_affine_transform_shift_via_grid_wrapc                 C   s|   t g dg dg}t d}t d|d dd df< |j|d d df< ttj||d|d|d d dd d df  d S )Nr   r!   r   rn   r  r    r   r   r   )r)   r*   r  r   r]   r   r+   r   )r.   r'   r   r  r   r   r   #test_affine_transform_shift_reflectx  s    

z<TestNdimageInterpolation.test_affine_transform_shift_reflectc                 C   s.   t dg}tj|dg|d}t|dg d S )Nr	   rK   r   r)   r*   r+   r   r   r   r   r   r   test_shift01  s    z%TestNdimageInterpolation.test_shift01c                 C   s0   t dg}tj|dg|d}t|g d d S )Nr   r	   rK   rc   r)   rM   r+   r   r   r   r   r   r   test_shift02  s    z%TestNdimageInterpolation.test_shift02c                 C   s.   t dg}tj|d|d}t|g d d S )Nr   r   rK   r   r  r   r   r   r   test_shift03  s    z%TestNdimageInterpolation.test_shift03c                 C   s0   t g d}tj|d|d}t|g d d S )Nre   r	   rK   rf   r  r   r   r   r   test_shift04  s    z%TestNdimageInterpolation.test_shift04c                 C   s   t jg dg dg dg|d}t jg dg dg dg|d}|jjdkrd|d| 8 }|d| 8 }tj|ddg|d}t|| d S )	NrW   r5   rc   rh   ri   r   r	   rK   r)   r*   r6   rj   r+   r   r   r   r   r   r   test_shift05  s"    

z%TestNdimageInterpolation.test_shift05c                 C   s   t jg dg dg dg|d}t jg dg dg dg|d}|jjdkrd|d| 8 }|d| 8 }d}||d d df< tj|ddg|||d	}t|| d S )
NrW   r5   rc   rh   ri   g      @r   r	   )r'   r%   r$   r  )r.   r'   r%   r6   r/   rI   r$   rP   r   r   r   test_shift_with_nonzero_cval  s&    

z5TestNdimageInterpolation.test_shift_with_nonzero_cvalc                 C   sP   t g dg dg dg}tj|ddg|d}t|g dg dg d	g d S )
Nre   rl   ro   r   r	   rK   rf   rp   rq   r  r   r   r   r   test_shift06  s    

z%TestNdimageInterpolation.test_shift06c                 C   sP   t g dg dg dg}tj|ddg|d}t|g dg dg dg d S )Nre   rl   ro   r	   r   rK   rs   r  r   r   r   r   test_shift07  s    

z%TestNdimageInterpolation.test_shift07c                 C   sP   t g dg dg dg}tj|ddg|d}t|g dg dg dg d S )	Nre   rl   ro   r	   rK   rs   rf   rp   r  r   r   r   r   test_shift08  s    

z%TestNdimageInterpolation.test_shift08c                 C   sn   t g dg dg dg}|dkr4tj||d}n|}tj|ddg|dd}t|g dg d	g d
g d S )Nre   rl   ro   r	   rK   Frv   rs   rf   rp   )r)   r*   r+   rN   r   r   r   r   r   r   test_shift09  s    

z%TestNdimageInterpolation.test_shift09c                 C   s>   t ddgddgg}ttj||d|dt j||dd d S )	Nr   r	   r    r!   r   r   r  r  )r)   r*   r   r+   r   r  r.   r   r'   r   r   r   r   test_shift_grid_wrap  s    
z-TestNdimageInterpolation.test_shift_grid_wrapc                 C   s:   t dd}ttj||d|dtj||d|d d S )N   )rn   r   r   r   r   )r)   r?   r   r   r+   r   r%  r   r   r   test_shift_grid_constant1  s
    z2TestNdimageInterpolation.test_shift_grid_constant1c                 C   sL   t jg dg dgtd}t g dg dg}ttj|dddd	| d S )
Nr   r   r5   r   r   )r   r   r   r	   r   )r)   r*   r@   r   r+   r   )r.   r   r   r   r   r   test_shift_grid_constant_order1  s    

z8TestNdimageInterpolation.test_shift_grid_constant_order1c                 C   sH   t g dg dg}ttj||jd|d|d d dd d df  d S )Nr   r  r   r   r   )r)   r*   r   r+   r   r]   )r.   r'   r   r   r   r   test_shift_reflect  s    
z+TestNdimageInterpolation.test_shift_reflectrw   FTc                 C   sr   t d}td||d}ttj||d d fi |d |d  ttj|| d d fi |d |d  d S )Nr	  r   )r%   r'   rw   r    r	   r   r   )r)   r?   dictr   r+   r   )r.   r'   rw   r   kwargsr   r   r   test_shift_nearest_boundary  s    
$&z4TestNdimageInterpolation.test_shift_nearest_boundaryc           
      C   s   t jdtddd}d}d}t|}t j|||d}tt|| g|j	 }t
j||||d| }	tt
j||||d|	dd	 d S )
Nr   r5   r   )r   gffffffr:   r;   r   r<   r   )r)   r?   r@   r   rB   rC   rD   r   r   r   r+   r   r   )
r.   r'   r%   r   r   rF   rG   r   r   r   r   r   r   test_shift_vs_padded%  s     
z-TestNdimageInterpolation.test_shift_vs_paddedc                 C   s   dddgfD ]}t ttddt}tj|||d}t	|j
d tt |dd d f dk tt |dd d f dt k tt |dd d f d	t k tt |dt k tt |d
t k qd S )Nr       )rn   rn   rK   )r   r   r   r   r'  rn   r   )r)   r*   listranger   r   r@   r+   zoomr   r]   r   alleps)r.   r'   zarrr   r   r   
test_zoom1:  s    ""z#TestNdimageInterpolation.test_zoom1c                 C   s2   t dd}tt|dd}t|| d S )Nr   r   r    r   )r)   r?   r   r+   r2  r   )r.   r6  rP   r   r   r   
test_zoom2F  s    z#TestNdimageInterpolation.test_zoom2c                 C   s^   t ddgg}t|d}t|d}t|t ddgddgg t|t g dg d S )Nr	   r    r   r   )r	   r	   r    r    )r)   r*   r+   r2  r   )r.   r6  r   r   r   r   r   
test_zoom3K  s
    z#TestNdimageInterpolation.test_zoom3c                 C   s   t jg dg dg dg|d}|jjdkr8|d| 8 }t 4}|td tj|ddgd	d
|d}W d    n1 sx0    Y  t	|d d dd d df | d S )Nr1   r   r   r5   rh   ri   r   r   r   r   rK   r    )
r)   Zasarrayr6   rj   r   r   r   r+   r   r   )r.   r'   r6   r/   r   rP   r   r   r   test_zoom_affine01S  s     
$z+TestNdimageInterpolation.test_zoom_affine01c                 C   s&   d}t jt||fd| dd d S )Nrm         ?r   r;   )r+   r2  r)   r  )r.   Zdimr   r   r   test_zoom_infinityc  s    z+TestNdimageInterpolation.test_zoom_infinityc                 C   s6   t d}d}tj||dd}t d}t|| d S )N)r	   rn   rn   )r;  r   r   r#   )r$   )r	   r   r   )r)   r  r+   r2  r   )r.   r6  r2  rP   refr   r   r   test_zoom_zoomfactor_oneh  s
    

z1TestNdimageInterpolation.test_zoom_zoomfactor_onec                 C   s*   t d}d}t||}t|jd d S )N)r!   r   r/  )gUUUUUU?gE]t?g(\?)r   r7      )r)   r  r+   r2  r   r]   )r.   r6  r2  rP   r   r   r   test_zoom_output_shape_roundoffq  s    
z8TestNdimageInterpolation.test_zoom_output_shape_roundoffr2  )r   )r!   rn   )rm   r    )rm   rm   )r   r   r   r   r   r   r   r   c                 C   sD   t jddgddggtd}ttj||d|dt |t | d S )Nr   r	   r    r!   r5   r9   )r)   r*   r@   r   r+   r2  kronrM   )r.   r2  r%   r   r   r   r   test_zoom_by_int_order0w  s    	
z0TestNdimageInterpolation.test_zoom_by_int_order0r]   r  )r   r   )r   r   r   r   r   c                 C   sF   t jt |td|}ttj||d|ddt |t 	| d S )Nr5   r   T)r'   r%   	grid_mode)
r)   r?   prodr@   r   r   r+   r2  rA  rM   )r.   r]   r2  r%   r   r   r   r   test_zoom_grid_by_int_order0  s
    z5TestNdimageInterpolation.test_zoom_grid_by_int_order0c                 C   sZ   t jdtdd}tjtdd$ tj|d|ddf W d    n1 sL0    Y  d S )	Nr   r5   )r!   r!   zIt is recommended to use mode)matchr    T)r%   rC  )	r)   r?   r@   r   r   Zwarnsr   r+   r2  )r.   r%   r   r   r   r   test_zoom_grid_mode_warnings  s
    z5TestNdimageInterpolation.test_zoom_grid_mode_warningsc                 C   s@   t jg dg dg dgt jd}tj|d|d}t|| d S )Nrs   r   r	   r	   r   r5   r   rK   r)   r*   r{   r+   rotater   r   r   r   r   test_rotate01  s    
z&TestNdimageInterpolation.test_rotate01c                 C   sh   t jg dg dg dgt jd}t jg dg dg dg dgt jd}tj|d|d}t|| d S )Nrs   )r   r	   r   r   r5   r   r   r   r   r	   r   Z   rK   rI  r.   r'   r/   rI   rP   r   r   r   test_rotate02  s    

z&TestNdimageInterpolation.test_rotate02c                 C   s   t jg dg dg dg|d}t jg dg dg dg dg dg|d}|jjdkrp|d| 8 }|d| 8 }tj|d|d	}t|| d S )
Nr   r   r   r   r   r   r	   r	   r   r   r5   rL  rM  rh   ri   rN  rK   )r)   r*   r6   rj   r+   rJ  r   r   r   r   r   test_rotate03  s&    

z&TestNdimageInterpolation.test_rotate03c                 C   sd   t jg dg dg dgt jd}t jg dg dg dgt jd}tj|dd|d}t|| d S )NrQ  rR  r5   r   r   r	   r   r   rN  F)r   r'   rI  rO  r   r   r   test_rotate04  s    

z&TestNdimageInterpolation.test_rotate04c                 C   s   t d}tdD ]>}t jg dg dg dg dgt jd|d d d d |f< qt jg dg dg dgt jd}tj|d|d	}tdD ] }t|d d d d |f | qd S )
N)r   r!   r!   r!   rL  rM  r5   rs   rH  rN  rK   r)   r   r1  r*   r{   r+   rJ  r   r.   r'   r/   irI   rP   r   r   r   test_rotate05  s$    


z&TestNdimageInterpolation.test_rotate05c                 C   s   t d}tdD ]8}t jg dg dg dgt jd|d d d d |f< qt jg dg dg dg dgt jd}tj|d|d	}tdD ] }t|d d d d |f | qd S )
N)r!   r   r!   r!   rs   rH  r5   rL  rM  rN  rK   rV  rW  r   r   r   test_rotate06  s$    


z&TestNdimageInterpolation.test_rotate06c                 C   s   t jg dg dg dggd t jd}| }t jg dg dg dg dg dggd t jd}|g d}tj|dd	|d
}t|| d S )NrQ  rR  r    r5   rL  rM  )r    r	   r   rN  r  )axesr'   r)   r*   r{   Z	transposer+   rJ  r   rO  r   r   r   test_rotate07  s,    

z&TestNdimageInterpolation.test_rotate07c                 C   s   t jg dg dg dggd t jd}| }t jg dg dg dggd t jd}| }tj|ddd|d	}t|| d S )
NrQ  rR  r    r5   rT  rN  r  F)r[  r   r'   r\  rO  r   r   r   test_rotate08  s(    

z&TestNdimageInterpolation.test_rotate08c                 C   sh   t jg dg dg dgd t jd}tt& tj|dd|jfd W d    n1 sZ0    Y  d S )NrQ  rR  r    r5   rN  r   )r[  )r)   r*   r{   r  r  r+   rJ  r   )r.   r/   r   r   r   test_rotate09  s    

z&TestNdimageInterpolation.test_rotate09c              
   C   s   t jdt jdd}t g dg dg dg dg dgg dg dg d	g d
g dgg dg dg dg dg dgg}tj|ddd}t|| d S )N-   r5   )r!   rn   r!   )r   r   r   )g|S2@g|S2@g>)!@)g9HD%@g9HD'@g9HD)@)gc~@gc~@gzLb @)g[ic6*@g[ic6,@g[ic6.@)g      5@g      6@g      7@)gRK<@gRK=@gRK>@)g,PA@g,PgB@g,PB@)gݠ'?@grn@@grnГ@@)gGA@gG9B@gGB@r   F)Zangler   )r)   r?   r{   r   r*   r+   rJ  r   )r.   r/   rI   rP   r   r   r   test_rotate10  s,    
z&TestNdimageInterpolation.test_rotate10c                 C   s4   t t dd}tt|dd}t|| d S )Nrn   )rn   r	      iL)r)   Ztiler?   r+   rJ  r   )r.   r   r   r   r   r   test_rotate_exact_180&  s    z.TestNdimageInterpolation.test_rotate_exact_180N)}__name__
__module____qualname__r   markZparametrizer0   r3   r1  rJ   r
   rQ   rS   rV   rX   rY   r_   ra   rd   rg   r)   r{   Z
complex128rk   rr   rt   ru   ry   r|   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   ZskipifsysplatformZintpitemsizer   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r  r  r  r  r  r  r  r  r  r  r  r  r  r   r!  r"  r#  r$  r&  r(  r)  r*  r-  r.  r7  r8  r9  r:  r<  r>  r@  rB  rE  rG  rK  rP  rS  rU  rY  rZ  r]  r^  r_  ra  rc  r   r   r   r   r      s  
































	
	
	










 





























	





	
	
	
	

	
	








r   c                  C   s*   t dd} tj| dt dd dS )zTicket #643r   r   r    r   rT   N)r)   r?   r   r+   r2  r  r   r   r   r   test_zoom_output_shape,  s    rk  )rh  r)   Znumpy.testingr   r   r   r   r   r   r   r   r  Zscipy.ndimager+    r
   r4  rB   r   rk  r   r   r   r   <module>   s4    
          