o
    Eb                     @   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 pue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   H/usr/lib/python3/dist-packages/scipy/ndimage/tests/test_interpolation.pyshift&      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      | d d fS )Nr   g?r   r   r   r   r   r   :   r    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r;dnd}
t	||	d|
d dS )z$Tests based on examples from gh-2640ir$   dtypei   i  )Z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(   r0   r   ynpadpad_modeZpaddedexpectedr?   r   r   r   test_boundary_spline_accuracyC   s   
z6TestNdimageInterpolation.test_boundary_spline_accuracyr!   r8   c                 C   s(   t g |}tj||d}t|d d S )Nr(   r	   r*   onesr,   spline_filterr   r/   r8   r(   r0   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 Nr   r   r   r   r   mapping      zDTestNdimageInterpolation.test_geometric_transform01.<locals>.mappingrK   r*   r+   r,   r-   shaper   r/   r(   r0   r\   rP   r   r   r   test_geometric_transform01}   s   z3TestNdimageInterpolation.test_geometric_transform01c                 C   :   t dg}dd }tj|||j|d}t|g d d S )Nr   c                 S   rZ   r[   r   r   r   r   r   r\      r]   zDTestNdimageInterpolation.test_geometric_transform02.<locals>.mappingrK   rW   r*   rM   r,   r-   r_   r   r`   r   r   r   test_geometric_transform02      z3TestNdimageInterpolation.test_geometric_transform02c                 C   rb   )Nr   c                 S   r2   Nr   r	   r   r   r   r   r   r\      r    zDTestNdimageInterpolation.test_geometric_transform03.<locals>.mappingrK   r   r	   r	   r	   rc   r`   r   r   r   test_geometric_transform03   re   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   r2   rf   r   r   r   r   r   r\      r    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r2|d| 8 }|d| 8 }dd }tj|||j|d}t|| d S )	NrW   r7   rg   c              ?c                 S      | d | d d fS rf   r   r   r   r   r   r\         zDTestNdimageInterpolation.test_geometric_transform05.<locals>.mappingrK   )r*   r+   r8   kindr,   r-   r_   r   )r/   r(   r8   r0   rI   r\   rP   r   r   r   test_geometric_transform05   s(   

z3TestNdimageInterpolation.test_geometric_transform05c                 C   X   t g dg dg dg}dd }tj|||j|d}t|g dg dg d	g d S )
Nri   r$   r6         r"   ru   r"   r6   c                 S   rn   rf   r   r   r   r   r   r\      ro   zDTestNdimageInterpolation.test_geometric_transform06.<locals>.mappingrK   rj   r   r$   r6   rt   r   r"   ru   r"   r^   r`   r   r   r   test_geometric_transform06      


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 )Nri   rs   rv   c                 S   s   | d d | d fS rf   r   r   r   r   r   r\      ro   zDTestNdimageInterpolation.test_geometric_transform07.<locals>.mappingrK   r   r   r   r   r^   r`   r   r   r   test_geometric_transform07   rz   z3TestNdimageInterpolation.test_geometric_transform07c                 C   rr   )
Nri   rs   rv   c                 S      | d d | d d fS rf   r   r   r   r   r   r\         zDTestNdimageInterpolation.test_geometric_transform08.<locals>.mappingrK   r{   rj   rw   r^   r`   r   r   r   test_geometric_transform08   rz   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 )Nri   rs   rv   c                 S   r}   rf   r   r   r   r   r   r\      r~   zDTestNdimageInterpolation.test_geometric_transform10.<locals>.mappingr	   rK   Fr(   	prefilterr{   rj   rw   )r*   r+   r,   rN   r-   r_   r   )r/   r(   r0   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\     r    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	   ru   r!   r6   r"   r$   r   r   c                 S   s   d| d  fS )Nr!   r   r   r   r   r   r   r\     r    zDTestNdimageInterpolation.test_geometric_transform14.<locals>.mappingr   rK   r3   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 )Nr3   c                 S   s   | d d fS r   r   r   r   r   r   r\     r    zDTestNdimageInterpolation.test_geometric_transform15.<locals>.mappingrt   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 )Nr3   ru   r6   r$   rt   g      "@
         c                 S   s   | d | d d fS Nr   r	   r!   r   r   r   r   r   r\   !  ro   zDTestNdimageInterpolation.test_geometric_transform16.<locals>.mappingr"   r!   rK   r	   r"   ru   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 )Nr3   r   r   r   r   r   c                 S   s   | d d | d fS Nr   r!   r	   r   r   r   r   r   r\   .  ro   zDTestNdimageInterpolation.test_geometric_transform17.<locals>.mappingr	   r   rK   r   r`   r   r   r   test_geometric_transform17(     
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 )
Nr3   r   r   c                 S      | d d | d d fS r   r   r   r   r   r   r\   ;  r~   zDTestNdimageInterpolation.test_geometric_transform18.<locals>.mappingr	   r!   rK   r	   r"   r   r`   r   r   r   test_geometric_transform185  r   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 )
Nr3   r   r   c                 S   s   | d | d d fS r   r   r   r   r   r   r\   H  ro   zDTestNdimageInterpolation.test_geometric_transform19.<locals>.mappingr"   rt   rK   .r!   r   r`   r   r   r   test_geometric_transform19B     
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 )
Nr3   r   r   c                 S   s   | d d | d fS r   r   r   r   r   r   r\   U  ro   zDTestNdimageInterpolation.test_geometric_transform20.<locals>.mappingr6   r   rK   r!   .r   r`   r   r   r   test_geometric_transform20O  r   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 )	Nr3   r   r   c                 S      | d d | d d fS r   r   r   r   r   r   r\   b  r~   zDTestNdimageInterpolation.test_geometric_transform21.<locals>.mappingr6   rt   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 )Nr3   r   r   c                 S   r   r   r   r   r   r   r   mapping1o  r~   zETestNdimageInterpolation.test_geometric_transform22.<locals>.mapping1c                 S   r   r   r   r   r   r   r   mapping2r  r~   zETestNdimageInterpolation.test_geometric_transform22.<locals>.mapping2r   rK   r"   r   )r*   r+   r   r,   r-   r   )r/   r(   r0   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 )
Nr3   r   r   c                 S   s   d| d d fS )Nr	   r   r!   r   r   r   r   r   r\        zDTestNdimageInterpolation.test_geometric_transform23.<locals>.mappingr!   rK   ru   r$   )r,   r-   astyper*   Z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 )Nr3   r   r   c                 S   s   || d | fS )Nr   r   )r   abr   r   r   r\     r   zDTestNdimageInterpolation.test_geometric_transform24.<locals>.mappingr   )r	   r   r!   )r(   Zextra_argumentsZextra_keywordsru   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   ru   r6   r7   c                 S   s   | d d | d d fS )Nr   r   r	   r   r   r   r   r   r\     r~   zWTestNdimageInterpolation.test_geometric_transform_grid_constant_order1.<locals>.mappingg      ?g      ?      ?r   g      @r#   r   r	   r&   r(   )r*   r+   rA   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   r7   r   c                 S   s   | d d | d d fS )Nr   皙?r	   gffffff@r   r   r   r   r   r\     r~   zLTestNdimageInterpolation.test_geometric_transform_vs_padded.<locals>.mapping   r<   r   r=   r>   )r*   r@   rA   reshaperC   rD   rE   tupleslicendimr,   r-   r   )
r/   r(   r&   r   r\   rG   rH   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 ]}tj|||j|d}|d u r2|n|}t	|dg q d S )Nr	   c                 S   rZ   r[   r   r   r   r   r   r\     r]   zcTestNdimageInterpolation.test_geometric_transform_endianness_with_output_parameter.<locals>.mappingrT   )
r*   r+   r8   newbyteorder
empty_liker   r,   r-   r_   r   )r/   r0   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   rZ   r[   r   r   r   r   r   r\     r]   zUTestNdimageInterpolation.test_geometric_transform_with_string_output.<locals>.mappingfrT   )r*   r+   r,   r-   r   r8   r   )r/   r0   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 )Nri   rs   rv   r{   rj   rw   rl   rm   r	   rK   )	r*   r+   r8   rp   indicesr_   r,   rB   r   )r/   r(   r8   r0   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 )Nri   rs   rv   r   rK   )	r*   r+   r   r_   r   r,   r   rB   r   )r/   r(   r0   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 )Nri   rs   rv   FrK   r	   r{   rj   rw   r	   r	   r!   r   r   r$   )r*   r+   r   r_   r,   rB   r   r   r/   r0   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 rF|n|}t
|| q6d S )Nr	   r!   r$   r6   r   rT   )r*   r+   r   r_   r8   r   r   r   r,   rB   r   )r/   r0   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,   rB   r   r8   r   r   r   r   r   'test_map_coordinates_with_string_output  s
   z@TestNdimageInterpolation.test_map_coordinates_with_string_outputwin32rt   z2do not run on 32 bit or windows (no sparse memory))reasonc              
   C   s   z1d}t j|d t jd||}d||d d |d d f< tj||d g|d ggdd W d S  tyC } ztd	|d }~ww )
Ni0u  r!   r7   r   r"   r   r	   rK   zNot enough memory available)	r*   emptyZfloat32r   r,   rB   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(   r0   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   rg   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 )Nri   r	   r   rK   rj   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r2|d| 8 }|d| 8 }tj|ddgddggddg|d	}t|| d S )
NrW   r7   rg   rl   rm   r	   r   r   rK   )r*   r+   r8   rp   r,   r   r   r/   r(   r8   r0   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 )Nri   rs   rv   r	   r   r   rK   rj   rw   rx   r   r   r   r   r   test_affine_transform06Z     


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 )	Nri   rs   rv   r	   r   r   rK   r{   r   r   r   r   r   test_affine_transform07e  r   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 )Nri   rs   rv   r	   r   r   rK   r{   rj   rw   r   r   r   r   r   test_affine_transform08p  r   z0TestNdimageInterpolation.test_affine_transform08c                 C   s|   t g dg dg dg}|dkrt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 )Nri   rs   rv   r	   rK   r   r   Fr   r{   rj   rw   )r*   r+   r,   rN   r   r   r/   r(   r0   r   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   r4   )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   r4   rK   r3   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 )Nr3   r   r   rt   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 )Nr3   r   r   r	   r   r!   r   rK   r"   ru   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 )	Nr3   r   r   r!   r   r	   r   rK   r   r   r   r   r   test_affine_transform14     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 )
Nr3   r   r   r!   r   r   rK   r	   r"   r   r   r   r   r   test_affine_transform15  r   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 )Nr3   r   r   r	           r   r   r   rK   .r!   r   r   r   r   r   test_affine_transform16     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 )Nr3   r   r   r   r   r	   r   rK   r!   .r   r   r   r   r   test_affine_transform17  r   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 )	Nr3   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 )
Nr3   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 )
Nr3   r   r   r   r!   r   rK   r	   r"   r   r   r   r   r   test_affine_transform20     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 )
Nr3   r   r   r!   r   r   rK   r	   r   r   r   r   r   r   test_affine_transform21  r  z0TestNdimageInterpolation.test_affine_transform21c                 C   s:   t g d}tj|dggdgd|d}t|g d d S )Nri   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 )Nri   r   r   r   rK   r!   rj   r   r   r   r   r   test_affine_transform23  s
   z0TestNdimageInterpolation.test_affine_transform23c                 C   x   t g d}t }|td tj|dgd|d}W d    n1 s%w   Y  tj|dggd|d}t|| d S )Nri   @The behavior of affine_transform with a 1-D array .* has changedr!   r   rK   r*   r+   r   filterUserWarningr,   r   r   r/   r(   r0   supr   r   r   r   r   test_affine_transform24     z0TestNdimageInterpolation.test_affine_transform24c                 C   r  )Nri   r  r   r   rK   r  r  r   r   r   test_affine_transform25  r  z0TestNdimageInterpolation.test_affine_transform25c                 C   s   t g dg dg dg}|dkrt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 q[d S )Nri   rs   rv   r	   rK   r!   r!   r	   )r   r   r	   Fr   r{   rj   rw   )r*   r+   r,   rN   eyerM   hstackvstackr   Zravelr   )r/   r(   r0   r   Ztform_originalZoffset_originaltform_h1tform_h2r   r   Zout3rP   r   r   r   test_affine_transform26  s6   




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 )Nri   rs   rv   r!   r  )ru   r!   r	   )
r*   r+   r  r  rM   r  assert_raises
ValueErrorr,   r   )r/   r0   r  r  r   r   r   test_affine_transform273  s   
z0TestNdimageInterpolation.test_affine_transform27c              	   C   s   t d}t |t ||j |j|j fD ]8}t }|td t	j
|ddg|d}W d    n1 s;w   Y  |d u rF|n|}t|ddgddgg qd S )N)r!   r!   r  r	   rT   )r*   rM   r   r   r8   r   r   r	  r
  r,   r   r   )r/   r0   rP   r  r   r   r   r   r   9test_affine_transform_1d_endianness_with_output_parameter<  s   
zRTestNdimageInterpolation.test_affine_transform_1d_endianness_with_output_parameterc                 C   sr   t dg}|j|j t |t ||j fD ]}tj|dgg|d}|d u r.|n|}t|dg qd S )Nr	   rT   )	r*   r+   r8   r   r   r   r,   r   r   )r/   r0   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}|j}tj|dgg|d t|d d | t	t
 tj|dgg|dd W d    d S 1 s?w   Y  d S )Nrt   r7   )   r	   rT   )r   )rU   r'   )r*   r@   r   rM   r_   r,   r   r   r   r   RuntimeError)r/   r0   rP   Zoshaper   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   r8   r   )r/   r0   rP   r   r   r   (test_affine_transform_with_string_outputd  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_wrapj  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   ru   r"  r!   r   r   r   )r*   r+   r$  r  r_   r   r,   r   )r/   r(   r   r&  r   r   r   #test_affine_transform_shift_reflecty  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   rg   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 )Nri   r	   rK   rj   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r2|d| 8 }|d| 8 }tj|ddg|d}t|| d S )	NrW   r7   rg   rl   rm   r   r	   rK   r*   r+   r8   rp   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r2|d| 8 }|d| 8 }d}||d d df< tj|ddg|||d	}t|| d S )
NrW   r7   rg   rl   rm   g      @r   r	   )r(   r&   r%   r0  )r/   r(   r&   r8   r0   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 )
Nri   rs   rv   r   r	   rK   rj   rw   rx   r*  r   r   r   r   test_shift06     


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 )Nri   rs   rv   r	   r   rK   r{   r*  r   r   r   r   test_shift07  r4  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 )	Nri   rs   rv   r	   rK   r{   rj   rw   r*  r   r   r   r   test_shift08  r4  z%TestNdimageInterpolation.test_shift08c                 C   sn   t g dg dg dg}|dkrtj||d}n|}tj|ddg|dd}t|g dg d	g d
g d S )Nri   rs   rv   r	   rK   Fr   r{   rj   rw   )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   )ru   r   r   r   r   )r*   r@   r   r   r,   r   r8  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   r7   r   r   )r   r   r   r	   r   )r*   r+   rA   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_reflectr   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(   r   r!   r	   r   r   )r*   r@   dictr   r,   r   )r/   r(   r   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   r7   r   )r   gffffffr;   r<   r   r=   r   )r*   r@   rA   r   rC   rD   rE   r   r   r   r,   r   r   )
r/   r(   r&   r   r   rG   rH   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 ]f}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!      )ru   ru   rK   )r   r   r   r   r:  ru   r   )r*   r+   listranger   r   rA   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,   rE  r   )r/   rI  rP   r   r   r   
test_zoom2G  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,   rE  r   )r/   rI  r   r   r   r   r   
test_zoom3L  s
   z#TestNdimageInterpolation.test_zoom3c                 C   s   t jg dg dg dg|d}|jjdkr|d| 8 }t }|td tj|ddgd	d
|d}W d    n1 s<w   Y  t	|d d dd d df | d S )Nr3   r   r   r7   rl   rm   r  r   r   r   rK   r!   )
r*   Zasarrayr8   rp   r   r	  r
  r,   r   r   )r/   r(   r8   r0   r  rP   r   r   r   test_zoom_affine01T  s"   
"z+TestNdimageInterpolation.test_zoom_affine01c                 C   s&   d}t jt||fd| dd d S )Nrt         ?r   r<   )r,   rE  r*   r$  )r/   Zdimr   r   r   test_zoom_infinityd  s   "z+TestNdimageInterpolation.test_zoom_infinityc                 C   s6   t d}d}tj||dd}t d}t|| d S )N)r	   ru   ru   )rN  r   r   r$   )r%   )r	   r   r   )r*   r$  r,   rE  r   )r/   rI  rE  rP   refr   r   r   test_zoom_zoomfactor_onei  s
   

z1TestNdimageInterpolation.test_zoom_zoomfactor_onec                 C   s*   t d}d}t||}t|jd d S )N)r"   r   rB  )gUUUUUU?gE]t?g(\?)r   r9      )r*   r$  r,   rE  r   r_   )r/   rI  rE  rP   r   r   r   test_zoom_output_shape_roundoffr  s   
z8TestNdimageInterpolation.test_zoom_output_shape_roundoffrE  )r   )r"   ru   )rt   r!   )rt   rt   )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"   r7   r:   )r*   r+   rA   r   r,   rE  kronrM   )r/   rE  r&   r   r   r   r   test_zoom_by_int_order0x  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 )Nr7   r   T)r(   r&   	grid_mode)
r*   r@   ZprodrA   r   r   r,   rE  rT  rM   )r/   r_   rE  r&   r   r   r   r   test_zoom_grid_by_int_order0  s
   z5TestNdimageInterpolation.test_zoom_grid_by_int_order0c                 C   s\   t jdtdd}tjtdd tj|d|ddf W d    d S 1 s'w   Y  d S )	Nr   r7   )r"   r"   zIt is recommended to use mode)matchr!   T)r&   rV  )	r*   r@   rA   r   r   Zwarnsr
  r,   rE  )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 )Nr{   r   r	   r	   r   r7   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 )Nr{   )r   r	   r   r   r7   r   r   r   r   r	   r   Z   rK   r[  r/   r(   r0   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r8|d| 8 }|d| 8 }tj|d|d	}t|| d S )
Nr   r   r   r   r   r   r	   r	   r   r   r7   r^  r_  rl   rm   r`  rK   )r*   r+   r8   rp   r,   r\  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 )Nrc  rd  r7   r   r   r	   r   r   r`  F)r   r(   r[  ra  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 | qFd S )
N)r   r"   r"   r"   r^  r_  r7   r{   rZ  r`  rK   r*   r   rD  r+   r   r,   r\  r   r/   r(   r0   irI   rP   r   r   r   test_rotate05  s&   


z&TestNdimageInterpolation.test_rotate05c                 C   s   t d}tdD ]}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 | qFd S )
N)r"   r   r"   r"   r{   rZ  r7   r^  r_  r`  rK   rh  ri  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 )Nrc  rd  r!   r7   r^  r_  )r!   r	   r   r`  r   )axesr(   r*   r+   r   Z	transposer,   r\  r   ra  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 )
Nrc  rd  r!   r7   rf  r`  r   F)rm  r   r(   rn  ra  r   r   r   test_rotate08  s(   

z&TestNdimageInterpolation.test_rotate08c                 C   sj   t jg dg dg dgd t jd}tt tj|dd|jfd W d    d S 1 s.w   Y  d S )Nrc  rd  r!   r7   r`  r   )rm  )r*   r+   r   r  r  r,   r\  r   )r/   r0   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-   r7   )r"   ru   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,   r\  r   )r/   r0   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 )Nru   )ru   r	      iL)r*   Ztiler@   r,   r\  r   )r/   r   r   r   r   r   test_rotate_exact_180'  s   z.TestNdimageInterpolation.test_rotate_exact_180N)}__name__
__module____qualname__r   ZmarkZparametrizer1   r5   rD  rJ   r
   rQ   rS   rV   rX   rY   ra   rd   rh   rk   r*   r   Z
complex128rq   ry   r|   r   r   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/  r1  r2  r3  r5  r6  r7  r9  r;  r<  r=  r@  rA  rJ  rK  rL  rM  rO  rQ  rS  rU  rW  rY  r]  rb  re  rg  rk  rl  ro  rp  rq  rs  ru  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,   rE  r$  r   r   r   r   test_zoom_output_shape-  s   r|  )ry  r*   Znumpy.testingr   r   r   r   r   r   r   r   r  Zscipy.ndimager,    r
   rG  rC   r   r|  r   r   r   r   <module>   s6     
          