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mZ d dlmZ ddlmZ G dd dZG dd	 d	Zd
d Zdd Zdd Zdd Zdd Zdd Zdd Zdd Zdd Zdd Zdd Zd d! Zd"d# Z d$d% Z!d&d' Z"d(d) Z#d*d+ Z$d,d- Z%d.d/ Z&d0d1 Z'd2d3 Z(d4d5 Z)d6d7 Z*d8d9 Z+d:d; Z,d<d= Z-d>d? Z.d@dA Z/dBdC Z0dDdE Z1dFdG Z2dHdI Z3dJdK Z4dLdM Z5dNdO Z6dPdQ Z7dRdS Z8dTdU Z9dVdW Z:dXdY Z;dZd[ Z<d\d] Z=d^d_ Z>d`da Z?dbdc Z@ddde ZAdfdg ZBdhdi ZCdjdk ZDdldm ZEdndo ZFdpdq ZGdrds ZHdtdu ZIdvdw ZJdxdy ZKdzd{ ZLd|d} ZMd~d ZNdd ZOdd ZPdd ZQdd ZRdd ZSdd ZTdd ZUdd ZVdd ZWdd ZXdd ZYdd ZZdd Z[dd Z\dd Z]dd Z^dd Z_dd Z`dd Zadd Zbdd Zcdd Zddd Zedd Zfdd Zgdd Zhdd Zidd Zjdd Zkdd Zldd Zmdd Zndd ZoddÄ Zpddń ZqddǄ ZrddɄ Zsdd˄ Ztdd̈́ Zuddτ Zvddф Zwddӄ ZxddՄ Zyddׄ Zzddل Z{ddۄ Z|dd݄ Z}dd߄ Z~dd ZG dd dZdS )    N)assert_assert_array_almost_equalassert_equalassert_almost_equalassert_array_equalsuppress_warnings)raises   )typesc                   @   s8   e Zd ZdZdd Zdd Zdd Zdd	 Zd
d ZdS )Test_measurements_statszDndimage._measurements._stats() is a utility used by other functions.c                 C   sx   g d}g d}ddg}dD ]V}t ||}t ||}tjj|||d\}}t|ddg t|dd	g qd S )
Nr   r	         r   r   r	   r	   r   r	   )   )r   r   labelsindexr         ?       @nparrayreshapendimage_measurementsZ_statsr   selfxr   r   shpcountssums r#   Y/var/www/sunrise/env/lib/python3.9/site-packages/scipy/ndimage/tests/test_measurements.pytest_a   s    
zTest_measurements_stats.test_ac                 C   sx   g d}g d}ddg}dD ]V}t ||}t ||}tjj|||d\}}t|ddg t|dd	g qd S )
Nr   r   r   	   r'   r   r'   r   r   r   r   r   r   r   r#   r#   r$   test_b   s    
zTest_measurements_stats.test_bc                 C   s   g d}g d}ddg}dD ]h}t ||}t ||}tjj|||dd\}}}t|ddg t|d	d
g t|dd
g qd S )Nr   r   r   r	   r   Tr   r   Zcenteredr   r   r         ?r   r   r   r   r   r    r!   r"   Zcentersr#   r#   r$   test_a_centered,   s    z'Test_measurements_stats.test_a_centeredc                 C   s   g d}g d}ddg}dD ]h}t ||}t ||}tjj|||dd\}}}t|ddg t|d	d
g t|dd
g qd S )Nr   r&   r   r'   r   Tr)   r   r   r   r*   r   r+   r#   r#   r$   test_b_centered9   s    z'Test_measurements_stats.test_b_centeredc                 C   s   g d}g d}ddg}dD ]h}t ||}t ||}tjj|||dd\}}}t|ddg t|d	d
g t|dd
g qd S )Nr   )        r.         "@r/   r.   r/   r   Tr)   r   r   r   r*   r   r+   r#   r#   r$   test_nonint_labelsF   s    z*Test_measurements_stats.test_nonint_labelsN)	__name__
__module____qualname____doc__r%   r(   r,   r-   r0   r#   r#   r#   r$   r      s   r   c                   @   s   e Zd ZdZdd ZdS )Test_measurements_selectzEndimage._measurements._select() is a utility used by other functions.c                 C   s  g d}g dddgfg dddgfg ddd	gfg}|D ]N\}}t jj|||d
}tt|dk t jj|||dd}tt|dk t|d ddg t jj|||dd}tt|dk t|d ddg t jj|||ddd}tt|dk t|d ddg t|d ddg t|d jjd t jj|||ddd}tt|dk t|d ddg t|d ddg t|d jjd q:d S )N)r   r	   r   r   r   r   r	   r&   r'   )r.   r.         @r6   r.   r6   r   T)r   r   find_maxr   )r   r   find_minr   )r   r   r8   Zfind_min_positions   i)r   r   r7   Zfind_max_positions)	r   r   _selectr   lenr   r   dtypekind)r   r   Zcasesr   r   resultr#   r#   r$   
test_basicW   sH    z#Test_measurements_select.test_basicN)r1   r2   r3   r4   r@   r#   r#   r#   r$   r5   T   s   r5   c                  C   s0   t g } t| \}}t|d t|d d S Nr	   r   onesr   labelr   r   dataoutnr#   r#   r$   test_label01z   s    

rI   c                  C   s0   t g } t| \}}t|d t|d d S )Nr   r   zerosr   rD   r   r   rE   r#   r#   r$   test_label02   s    

rL   c                  C   s4   t dg} t| \}}t|dg t|d d S rA   rB   rE   r#   r#   r$   test_label03   s    rM   c                  C   s4   t dg} t| \}}t|dg t|d d S Nr	   r   rJ   rE   r#   r#   r$   test_label04   s    rO   c                  C   s6   t dg} t| \}}t|g d t|d d S )N   )r	   r	   r	   r	   r	   r	   rB   rE   r#   r#   r$   test_label05   s    rQ   c                  C   s8   t g d} t| \}}t|g d t|d d S )N)r	   r   r	   r	   r   r	   r	   r   r   r   r   r9   r9   r   r   r   rD   r   r   rE   r#   r#   r$   test_label06   s    rT   c               	   C   sx   t g dg dg dg dg dg dg} t| \}}t|g dg dg dg dg dg dg t|d d S )Nr   r   r   r   r   r   r   rS   rE   r#   r#   r$   test_label07   s     

rV   c               	   C   sx   t g dg dg dg dg dg dg} t| \}}t|g dg dg dg dg dg d	g t|d
 d S N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   r9   r9   r   r   r   r   r   r   r   r   r   r   r   rS   rE   r#   r#   r$   test_label08   s     

ra   c               	   C   s   t g dg dg dg dg dg dg} tdd}t| |\}}t|g dg dg dg d	g d	g d
g t|d d S )NrX   rY   rZ   r[   r\   r   r]   r^   )r   r   r   r   r   r   )r   r   r   r9   r9   r   r9   r   r   r   generate_binary_structurerD   r   r   rF   structrG   rH   r#   r#   r$   test_label09   s"    

rf   c                  C   sn   t g dg dg dg dg} tdd}t| |\}}t|g dg dg dg dg t|d d S )NrU   )r   r	   r	   r   r	   r   )r   r	   r	   r	   r	   r   r   r	   rb   rd   r#   r#   r$   test_label10   s    

rg   c               
   C   s   t D ]~} tg dg dg dg dg dg dg| }t|\}}g dg dg dg dg dg d	g}t|| t|d
 qd S rW   r
   r   r   r   rD   r   r   typerF   rG   rH   expectedr#   r#   r$   test_label11   s(    

rl   c               
   C   s   t D ]~} tg dg dg dg dg dg dg| }tj||d}g dg dg dg d	g d	g d
g}t|| t|d qd S )NrX   rY   rZ   r[   r\   outputr]   r^   r_   r`   r   rh   )rj   rF   rH   rk   r#   r#   r$   test_label11_inplace   s(    

ro   c               	   C   s|   t D ]r} tg dg dg dg dg dg| }t|\}}g dg dg dg dg dg}t|| t|d qd S )N)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	   rh   ri   r#   r#   r$   test_label12  s$    

rp   c                  C   sp   t D ]f} tg dg dg dg dg| }t|\}}g dg dg dg dg}t|| t|d qd S )N)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	   rh   ri   r#   r#   r$   test_label13  s     

rq   c                  C   sL   t dg} tD ]6}t jdg|d}tj| |d}t|d t|d qd S )NrP   r=   rm   r	   )r   rC   r
   rK   r   rD   r   r   rF   trn   rH   r#   r#   r$   test_label_output_typed-  s    
ru   c                  C   sD   t dg} tD ].}tj| |d\}}t|d |j|ksJ qd S )NrP   rm   r	   )r   rC   r
   r   rD   r   r=   rs   r#   r#   r$   test_label_output_dtype6  s
    
rv   c                  C   s>   t dg} tD ](}t dg|}tttftj| |d qd S )NrP   
   rm   )	r   rC   r
   rK   assert_raisesRuntimeError
ValueErrorr   rD   )rF   rt   rn   r#   r#   r$   test_label_output_wrong_size>  s    r{   c               	   C   s  t tjtjtdd} t tjtjtdd}t tjtjtdd}| d} |d}|d}d}t| j	d D ]t}| |d d d d f }t|j	d D ]J}||d d d d f }t
t||d ||d d d d f  |d7 }qqd S )	NrF   zlabel_inputs.txtzlabel_strels.txtzlabel_results.txt)   r}   )r|   r9   r9   r   r	   )r   Zloadtxtospathjoindirname__file__r   rangeshaper   r   rD   )rF   Zstrelsresultsrr:   djsr#   r#   r$   test_label_structuring_elementsF  s*    




(r   c                  C   sH   ddd} t t jt dkrDd}t jt ||}| | d S )Nffffff?r   c                 S   s8   | |k}t |j}t|t||\}}t|}d S )N)r<   r   r   rD   rc   find_objects)imgZthreshsizemaskZrankZlaco_r#   r#   r$   SEZ  s    

ztest_ticket_742.<locals>.SEr:   )r9     r   )r   r   )r   r=   Zintprandomrandprodr   )r   r   ar#   r#   r$   test_ticket_742Y  s
    
r   c                  C   sl   t d} d| ddddf< d| ddddf< d| d< d| d< d| d< t| t d	d dkshJ dS )
z/Github issue #3025 - improper merging of labels)<   i@  r	   N  i  )$   r   )#   i  )r   i  )r9   r9   )r   rK   r   rD   rC   )r   r#   r#   r$   test_gh_issue_3025h  s    
r   c                  C   sD   t jdd} t| dk\}}t|jt jt jfv  t	| d S )Nrw   r*   )
r   r   r   r   rD   r   r=   int32int64r   )Z
test_arrayrD   Zno_featuresr#   r#   r$   test_label_default_dtypes  s    r   c                  C   s*   t jg td} t| }t|dgk d S )Nrr   r#   )r   rC   intr   r   r   rF   rG   r#   r#   r$   test_find_objects01{  s    
r   c                  C   s(   t jg td} t| }t|g k d S )Nrr   )r   rK   r   r   r   r   r   r#   r#   r$   test_find_objects02  s    
r   c                  C   s4   t jdgtd} t| }t|tddd fg d S )Nr	   rr   r   r   rC   r   r   r   r   slicer   r#   r#   r$   test_find_objects03  s    
r   c                  C   s(   t jdgtd} t| }t|g  d S )Nr	   rr   )r   rK   r   r   r   r   r   r#   r#   r$   test_find_objects04  s    
r   c                  C   s4   t jdgtd} t| }t|tddd fg d S )NrP   rr   r   r   r   r#   r#   r$   test_find_objects05  s    
r   c                  C   sJ   t g d} t| }t|tddd ftddd ftddd fg d S )NrR   r   r	   r   r   rP   r   r   r   r   r   r   r   r   r#   r#   r$   test_find_objects06  s    
r   c               	   C   sF   t g dg dg dg dg dg dg} t| }t|g  d S )NrU   )r   r   r   r   r   r   r#   r#   r$   test_find_objects07  s    

r   c               
   C   s   t g dg dg dg dg dg dg} t| }t|tddd tddd ftddd td	d
d ftdd
d tdd	d ftd
dd tdd
d fg d S )NrX   r]   r^   r_   r`   r   r	   r9   r   rP   r   r   r   r#   r#   r$   test_find_objects08  s    

r   c               
   C   s   t g dg dg dg dg dg dg} t| }t|tddd tddd ftddd td	d
d fd td
dd tdd
d fg d S )NrX   r]   r^   rU   r`   r   r	   r9   r   rP   r   r   r   r#   r#   r$   test_find_objects09  s    

r   c               	   C   s   t g dg dg dg dg dg dg} tj| dd}g d}tt| | i }|D ]}t | |k||< q^tj| dd}t|| d	S )
z Test dictionary keys and entriesrX   r]   r^   rU   r`   r   )Zignore_value)r	   r   r   N)r   r   r   value_indicesr   listkeyswhere)rF   viZ	true_keysZtruevikr#   r#   r$   test_value_indices01  s    
r   c                  C   sN   t jdt jd} d}tt|d t|  W d   n1 s@0    Y  dS )zTest input checking)rP   r   rr   z(Parameter 'arr' must be an integer array)matchN)r   rK   float32rx   rz   r   r   )rF   msgr#   r#   r$   test_value_indices02  s    r   c                  C   s   dD ]} t jddg ddg  ddg  t jd| }t |}t|}tt|	 t| |D ] }t 
||k}t|| | qfqdS )z2Test different input array shapes, from 1-D to 4-D))r   )   r   )r9   r9   r   )r9   r9   r   r      r	   r   r9   rr   N)r   r   r   r   uniquer   r   r   r   r   r   )r   r   ZtrueKeysr   r   ZtrueNdxr#   r#   r$   test_value_indices03  s    0

r   c                  C   s.   t D ]$} tg | }t|}t|d qd S )Nr.   )r
   r   r   r   sumr   rj   inputrn   r#   r#   r$   
test_sum01  s    
r   c                  C   s2   t D ](} tddg| }t|}t|d qd S )Nr   r   r.   )r
   r   rK   r   r   r   r   r#   r#   r$   
test_sum02  s    
r   c                  C   s.   t D ]$} tg | }t|}t|d qd S )Nr   )r
   r   rC   r   r   r   r   r#   r#   r$   
test_sum03  s    
r   c                  C   s2   t D ](} tddg| }t|}t|d qd S )Nr	   r         @r
   r   r   r   r   r   r   r#   r#   r$   
test_sum04  s    
r   c                  C   s:   t D ]0} tddgddgg| }t|}t|d qd S )Nr	   r   r9   r   g      $@r   r   r#   r#   r$   
test_sum05  s    
r   c                  C   s>   t g t} tD ](}t g |}tj|| d}t|d qd S )Nr   r.   r   r   boolr
   r   r   r   r   rj   r   rn   r#   r#   r$   
test_sum06  s
    r   c                  C   sF   t ddgt} tD ],}t ddg|}tj|| d}t|d qd S )Nr   r   r   r.   )r   rC   r   r
   rK   r   r   r   r   r#   r#   r$   
test_sum07  s
    r   c                  C   sF   t ddgt} tD ],}t ddg|}tj|| d}t|d qd S )Nr	   r   r   r   r   r   r   r#   r#   r$   
test_sum08#  s
    r   c                  C   sN   t ddgt} tD ]4}t ddgddgg|}tj|| d}t|d qd S )Nr	   r   r   r9   r   r         @)r   r   r   r
   r   r   r   r   r#   r#   r$   
test_sum09+  s
    r   c                  C   sD   t ddgt} t ddgddggt}tj|| d}t|d d S Nr	   r   r   r9   r   r          @)r   r   r   r   r   r   r   r   rn   r#   r#   r$   
test_sum103  s    r   c                  C   sR   t ddgt j} tD ]6}t ddgddgg|}tj|| dd}t|d qd S )Nr	   r   r9   r   r         @)r   r   int8r
   r   r   r   r   r#   r#   r$   
test_sum11:  s    r   c                  C   sb   t ddgddggt j} tD ]>}t ddgddgg|}tj|| g dd}t|g d qd S Nr	   r   r   r9   r      r   r   )r   r.   g      @)r   r   r   r
   r   r   r   r   r#   r#   r$   
test_sum12C  s
    r   c                  C   s   t ddgddggt j} tD ]b}t ddgddgg|}tj|| g dd}tj|| g dd}||k srJ t|g d qd S r   )	r   r   r   r
   r   r   Z
sum_labelsallr   )r   rj   r   Z
output_sumZoutput_labelsr#   r#   r$   test_sum_labelsK  s    
r   c                  C   sN   t ddgt} tD ]4}t ddgddgg|}tj|| d}t|d qd S r   )r   r   r   r
   r   meanr   r   r#   r#   r$   test_mean01W  s
    r   c                  C   sD   t ddgt} t ddgddggt}tj|| d}t|d d S Nr	   r   r   r9   r   r   r   )r   r   r   r   r   r   r   r#   r#   r$   test_mean02_  s    r   c                  C   sN   t ddg} tD ]6}t ddgddgg|}tj|| dd}t|d qd S )Nr	   r   r9   r   r   r   )r   r   r
   r   r   r   r   r#   r#   r$   test_mean03f  s    r   c                  C   s   t ddgddggt j} t jddn tD ]X}t ddgddgg|}tj|| g dd}t|d	dg d
dg tt 	|d  q,W d    n1 s0    Y  d S )Nr	   r   r   ignorer   r9   r   r   r   r         @)
r   r   r   errstater
   r   r   r   r   isnanr   r#   r#   r$   test_mean04o  s    r   c                  C   sN   t ddgt} tD ]4}t ddgddgg|}tj|| d}t|d qd S r   )r   r   r   r
   r   minimumr   r   r#   r#   r$   test_minimum01z  s
    r   c                  C   sD   t ddgt} t ddgddggt}tj|| d}t|d d S Nr	   r   r   r   r   r   )r   r   r   r   r   r   r   r#   r#   r$   test_minimum02  s    r   c                  C   sN   t ddg} tD ]6}t ddgddgg|}tj|| dd}t|d qd S )Nr	   r   r9   r   r   r   )r   r   r
   r   r   r   r   r#   r#   r$   test_minimum03  s    r   c                  C   s^   t ddgddgg} tD ]>}t ddgddgg|}tj|| g dd}t|g d qd S )Nr	   r   r9   r   r   r9   r   r   )r   r   r.   )r   r   r
   r   r   r   r   r#   r#   r$   test_minimum04  s    r   c                  C   sN   t ddgt} tD ]4}t ddgddgg|}tj|| d}t|d qd S )Nr	   r   r   r9   r   r   r   )r   r   r   r
   r   maximumr   r   r#   r#   r$   test_maximum01  s
    r   c                  C   sD   t ddgt} t ddgddggt}tj|| d}t|d d S r   )r   r   r   r   r   r   r   r#   r#   r$   test_maximum02  s    r   c                  C   sN   t ddg} tD ]6}t ddgddgg|}tj|| dd}t|d qd S )Nr	   r   r9   r   r   r   )r   r   r
   r   r   r   r   r#   r#   r$   test_maximum03  s    r   c                  C   s^   t ddgddgg} tD ]>}t ddgddgg|}tj|| g dd}t|g d qd S )Nr	   r   r9   r   r   r   )r   r   r.   )r   r   r
   r   r   r   r   r#   r#   r$   test_maximum04  s    r   c                  C   s"   t g d} tt| d d S )N)r|   r|   )r   r   r   r   r   )r   r#   r#   r$   test_maximum05  s    r   c                  C   sj   t g dg dg dg dg} t g dg dg dg dg}tj| |g dd	}t|g d
 d S )Nr	   r   r   r	   rP   r9   r   r   r   r   r   r}   r'   r9   r   r   r	   r	   r   r   r   r   r   r   r9   r9   r   r   )r	   r   r9   r   )r   r   r   )r   r   r   medianr   r   r   rn   r#   r#   r$   test_median01  s    

r  c                  C   s:   t g dg dg dg dg} t| }t|d d S )Nr   r   r   r   r   r   r   r   r   r   r   rn   r#   r#   r$   test_median02  s    

r  c                  C   s`   t g dg dg dg dg} t g dg dg dg dg}tj| |d}t|d	 d S )
Nr   r   r   r   r   r   r   r   r   r  r   r#   r#   r$   test_median03  s    

r  c                  C   s:   t jddgtd} tj| t ddgd}t|dg d S )Nr	   rr   r   r   r   )r   asarrayr   r   r   rC   r   r  r#   r#   r$   test_median_gh12836_bool  s    r  c                  C   s<   t jddgt jd} tj| t ddgd}t|dg d S )NA   F   rr   r  r	   r   g     P@)r   r  r   r   r   rC   r   r  r#   r#   r$   test_median_no_int_overflow  s    r  c               
   C   s   t jddp tD ]Z} t g | }t &}|td t|}W d    n1 sT0    Y  t	t 
| qW d    n1 s0    Y  d S Nr   r   zMean of empty slice)r   r   r
   r   r   filterRuntimeWarningr   variancer   r   rj   r   suprn   r#   r#   r$   test_variance01  s    (r  c                  C   s0   t D ]&} tdg| }t|}t|d qd S Nr	   r.   r
   r   r   r   r  r   r   r#   r#   r$   test_variance02  s    
r  c                  C   s2   t D ](} tddg| }t|}t|d qd S Nr	   r9   r   r  r   r#   r#   r$   test_variance03  s    
r  c                  C   s(   t ddgt} t| }t|d d S )Nr	   r   g      ?)r   r   r   r   r  r   r   rn   r#   r#   r$   test_variance04  s    
r  c                  C   s>   g d} t D ],}tg d|}t|| d}t|d qd S N)r   r   r9   )r	   r9   r   r   r   r  r   r#   r#   r$   test_variance05  s
    r  c                  C   sr   g d} t jddJ tD ]4}t g d|}t|| g d}t|g d qW d    n1 sd0    Y  d S N)r   r   r9   r9   r   r   r   )r	   r9   r   rw   r   )r   r9   r   )r   r   r.   )r   r   r
   r   r   r  r   r   r#   r#   r$   test_variance06  s    r  c               
   C   s   t jddp tD ]Z} t g | }t &}|td t|}W d    n1 sT0    Y  t	t 
| qW d    n1 s0    Y  d S r  )r   r   r
   r   r   r  r  r   standard_deviationr   r   r  r#   r#   r$   test_standard_deviation01"  s    (r  c                  C   s0   t D ]&} tdg| }t|}t|d qd S r  r
   r   r   r   r  r   r   r#   r#   r$   test_standard_deviation02,  s    
r!  c                  C   s8   t D ].} tddg| }t|}t|td qd S r  )r
   r   r   r   r  r   sqrtr   r#   r#   r$   test_standard_deviation033  s    
r#  c                  C   s(   t ddgt} t| }t|d d S )Nr	   r   r*   )r   r   r   r   r  r   r  r#   r#   r$   test_standard_deviation04:  s    
r$  c                  C   s>   g d} t D ],}tg d|}t|| d}t|d qd S r  r   r   r#   r#   r$   test_standard_deviation05@  s
    r%  c                  C   sr   g d} t jddJ tD ]4}t g d|}t|| g d}t|g d qW d    n1 sd0    Y  d S r  r   r   r
   r   r   r  r   r   r#   r#   r$   test_standard_deviation06H  s    r'  c                  C   sj   dg} t jddD tD ].}t dg|}t|| dg}t|dg qW d    n1 s\0    Y  d S )Nr	   r   r   gl `yr   r&  r   r#   r#   r$   test_standard_deviation07Q  s    r(  c                  C   sN   t ddgt} tD ]4}t ddgddgg|}tj|| d}t|d qd S )Nr	   r   r   r9   r   r   r   r   )r   r   r   r
   r   minimum_positionr   r   r#   r#   r$   test_minimum_position01Z  s
    r+  c                  C   s@   t D ]6} tg dg dg dg| }t|}t|d qd S NrP   r   r   rP   r9   r}   r   r   r	   rP   r	   r	   r	   r   r
   r   r   r   r*  r   r   r#   r#   r$   test_minimum_position02b  s    

r2  c                  C   s6   t g dg dg dgt} t| }t|d d S r,  r   r   r   r   r*  r   r  r#   r#   r$   test_minimum_position03k  s    

r4  c                  C   s6   t g dg dg dgt} t| }t|d d S )Nr-  )r9   r}   r	   r   r/  r)  r3  r  r#   r#   r$   test_minimum_position04s  s    

r5  c                  C   sJ   g d} t D ]8}tg dg dg dg|}t|| }t|d qd S )Nr	   r   r   r   r-  r.  )r	   rP   r   r9   )r   r   r1  r   r#   r#   r$   test_minimum_position05{  s    
r7  c                  C   sL   g d} t D ]:}tg dg dg dg|}t|| d}t|d qd S )Nr	   r   r9   r   r-  r.  r/  r   r   r	   r1  r   r#   r#   r$   test_minimum_position06  s    
r:  c                  C   sb   g d} t D ]P}tg dg dg dg|}t|| ddg}t|d d t|d	 d
 qd S )Nr8  r-  r.  r/  r   r9   r   r9  r	   r0  r1  r   r#   r#   r$   test_minimum_position07  s    
r;  c                  C   sN   t ddgt} tD ]4}t ddgddgg|}tj|| d}t|d qd S )Nr	   r   r   r9   r   r   )r	   r   )r   r   r   r
   r   maximum_positionr   r   r#   r#   r$   test_maximum_position01  s    r=  c                  C   s@   t D ]6} tg dg dg dg| }t|}t|d qd S )Nr-  r9   r}   r   r   r/  r0  r
   r   r   r   r<  r   r   r#   r#   r$   test_maximum_position02  s    

r@  c                  C   s6   t g dg dg dgt} t| }t|d d S )Nr-  r>  r/  r)  )r   r   r   r   r<  r   r  r#   r#   r$   test_maximum_position03  s    

rA  c                  C   sJ   g d} t D ]8}tg dg dg dg|}t|| }t|d qd S )Nr6  r-  r>  r/  r	   r	   r?  r   r#   r#   r$   test_maximum_position04  s    
rC  c                  C   sL   g d} t D ]:}tg dg dg dg|}t|| d}t|d qd S )Nr6  r-  r>  r/  r	   r)  r?  r   r#   r#   r$   test_maximum_position05  s    
rD  c                  C   sb   g d} t D ]P}tg dg dg dg|}t|| ddg}t|d d t|d d	 qd S )
Nr6  r-  r>  r/  r	   r   r   r)  rB  r?  r   r#   r#   r$   test_maximum_position06  s    
rE  c                  C   sh   t g d} tD ]P}t g dg dg dg|}t|| ddg}t|d d t|d	 d
 qd S )N)r   r   r.         @r-  r>  r/  r   rF  r   r)  r	   )r   r9   )r   r   r
   r   r<  r   r   r#   r#   r$   test_maximum_position07  s    
rG  c                  C   s   t ddgt} tD ]t}t ddgddgg|}tj|| d}tj|| d}tj|| d}tj|| d}tj	|| d}t
|||||f qd S )Nr	   r   r   r9   r   r   )r   r   r   r
   r   extremar   r   r*  r<  r   r   rj   r   Zoutput1Zoutput2Zoutput3Zoutput4Zoutput5r#   r#   r$   test_extrema01  s    rJ  c                  C   s   t ddg} tD ]~}t ddgddgg|}tj|| dd}tj|| dd}tj|| dd}tj|| dd}tj|| dd}t	|||||f qd S )Nr	   r   r9   r   r   )
r   r   r
   r   rH  r   r   r*  r<  r   rI  r#   r#   r$   test_extrema02  s&    rK  c                  C   s   t ddgddgg} tD ]}t ddgddgg|}tj|| g dd}tj|| g dd}tj|| g dd}tj|| g dd}tj|| g dd}t	|d | t	|d | t	|d | t	|d | qd S )Nr	   r   r9   r   r   r   r   )
r   r   r
   r   rH  r   r   r*  r<  r   rI  r#   r#   r$   test_extrema03  s,    rL  c                  C   s   g d} t D ]}tg dg dg dg|}t|| ddg}t|| ddg}t|| ddg}t|| ddg}t|| ddg}t	|d | t	|d | t	|d | t	|d | qd S )	Nr6  r-  r>  r/  r	   r   r   r9   )
r
   r   r   r   rH  r   r   r*  r<  r   rI  r#   r#   r$   test_extrema04  s*    
rM  c                  C   sB   ddg} t D ]0}tddgddgg|}t|}t||  qd S )Nr.   r	   r   r
   r   r   r   center_of_massr   rk   rj   r   rn   r#   r#   r$   test_center_of_mass01(  s
    
rQ  c                  C   sB   ddg} t D ]0}tddgddgg|}t|}t||  qd S rN   rN  rP  r#   r#   r$   test_center_of_mass020  s
    
rR  c                  C   sB   ddg} t D ]0}tddgddgg|}t|}t||  qd S )Nr   r	   rN  rP  r#   r#   r$   test_center_of_mass038  s
    
rS  c                  C   sB   ddg} t D ]0}tddgddgg|}t|}t||  qd S rN   rN  rP  r#   r#   r$   test_center_of_mass04@  s
    
rT  c                  C   sB   ddg} t D ]0}tddgddgg|}t|}t||  qd S )Nr*   r	   rN  rP  r#   r#   r$   test_center_of_mass05H  s
    
rU  c                  C   s8   ddg} t ddgddggt}t|}t||  d S )Nr*   r	   r   r9   r   r   r   r   rO  r   rk   r   rn   r#   r#   r$   test_center_of_mass06P  s    
rX  c                  C   sB   ddg} ddg}t ddgddggt}t|| }t|| d S )Nr	   r   r*   r.   r   r9   rV  r   rk   r   rn   r#   r#   r$   test_center_of_mass07W  s
    rZ  c                  C   sD   ddg} ddg}t ddgddggt}t|| d}t|| d S )Nr	   r   r*   r   rP   r9   rV  rY  r#   r#   r$   test_center_of_mass08_  s
    r[  c                  C   sH   ddg} ddg}t ddgddggt}t|| ddg}t|| d S )Nr	   r   )r*   r.   )r*   r   rV  rY  r#   r#   r$   test_center_of_mass09g  s
    r\  c                  C   s2   t d} t d}t|ddd}t||  d S )Nrw   r   )r   rC   Zaranger   	histogramr   rW  r#   r#   r$   test_histogram01o  s    

r^  c                  C   s@   g d} g d}t g d}t|ddd| d}t|| d S )N)r	   r	   r	   r	   r   r   r   r   )r   r   r   r	   r	   )r	   r	   r9   r   r9   r9   r9   r9   r   r   rP   r	   r   r   r   r]  r   rY  r#   r#   r$   test_histogram02v  s
    r`  c                  C   sZ   g d} g d}g d}t g d}t|ddd| d}t|d | t|d	 | d S )
N)r	   r   r	   r	   r   r   r   r   )r   r	   r   r	   r	   )r   r   r   r9   r   )r	   r	   r9   r   r9   rP   r9   r9   r   r   rP   r0  r	   r_  )r   Z	expected1Z	expected2r   rn   r#   r#   r$   test_histogram03~  s    ra  c                  C   s  t g dg dg dg} t g dg dg dg}tj| |ddgd}t|d	d
g tj| |ddgd}t|ddg tj| |ddgd}t|t ddg tj	| |ddgd}t|d	d
g tj
| |ddgd}t|ddg tj| |ddgd}t|ddg d S )N)rP   r   r   r   r   )r   r'   r   r   r   )r   r   r   r9   rP   )r	   r	   r   r   r   )r   r   r   r   r   r	   r   r   r6   r   r   r   rP   r9   r'   )r   r   r   r   r   r  r  r   r"  r   r   r   )r   Zlblr   varZstdZmedminmaxr#   r#   r$   test_stat_funcs_2d  s    re  c                   @   sL   e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd Zdd Z	dd Z
dS )TestWatershedIftc                 C   s   t g dg dg dg dg dg dg dg dgt j}t g dg dg dg dg dg dg dg dgt j}tj||g dg dg dgd}g dg d	g d	g d	g d	g d	g dg dg}t|| d S )
N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	   	structurer|   r|   r|   r|   r|   r|   r|   r|   r	   r	   r	   r	   r	   r|   r   r   uint8r   r   watershed_iftr   r   rF   markersrG   rk   r#   r#   r$   test_watershed_ift01  sH    

z%TestWatershedIft.test_watershed_ift01c                 C   s   t g dg dg dg dg dg dg dg dgt j}t g dg dg dg dg dg dg dg dgt j}t||}g dg dg dg dg dg dg dg dg}t|| d S )	Nrg  rh  ri  rj  rk  ro  )r|   r|   r	   r	   r	   r|   r|   rp  rq  rt  r#   r#   r$   test_watershed_ift02  sB    

z%TestWatershedIft.test_watershed_ift02c              
   C   s   t g dg dg dg dg dg dg dgt j}t g dg dg dg dg dg dg dgt j}t||}g dg dg dg dg dg dg dg}t|| d S )	Nrg  rh  r   r	   r   r	   r   r	   r   r   r   r   r   r9   r   r   r   r   r   r   r   r   r|   ro  )r|   r|   r   r|   r9   r|   r|   r|   r   r   r9   r9   r9   r|   rq  rt  r#   r#   r$   test_watershed_ift03  s<    

z%TestWatershedIft.test_watershed_ift03c              
   C   s   t g dg dg dg dg dg dg dgt j}t g dg dg dg dg dg dg dgt j}tj||g dg dg dgd}g dg d	g d	g d	g d	g d	g dg}t|| d S )
Nrg  rh  rx  ry  rz  rl  rm  ro  r{  rq  rt  r#   r#   r$   test_watershed_ift04  sF    

z%TestWatershedIft.test_watershed_ift04c              
   C   s   t g dg dg dg dg dg dg dgt j}t g dg dg dg dg dg dg dgt j}tj||g dg dg dgd}g dg d	g d	g d	g d	g d	g dg}t|| d S )
Nrg  rh  rx  )r   r   r9   r   r   r   r   rz  rl  rm  ro  )r|   r9   r9   r   r   r   r|   rq  rt  r#   r#   r$   test_watershed_ift05  sF    

z%TestWatershedIft.test_watershed_ift05c              	   C   s   t g dg dg dg dg dg dgt j}t g dg dg dg dg dg dgt j}tj||g dg dg dgd}g dg dg dg dg d	g d	g}t|| d S )
Nri  rh  rg  rj  rk  rl  rm  rp  ro  rq  rt  r#   r#   r$   test_watershed_ift06/  s@    

z%TestWatershedIft.test_watershed_ift06c              	   C   s   d}t j|t jd}| }t g dg dg dg dg dg dgt j|d< t g dg dg dg dg dg dgt j}t j|t jd}| }tj||g d	g d	g d	g|d
 g dg dg dg dg dg dg}t	|| d S )N)r}   r   rr   ri  rh  rg  .rj  rk  rl  )rn  rn   rp  ro  )
r   rK   rr  Z	transposer   r   int16r   rs  r   )r   r   rF   ru  rG   rk   r#   r#   r$   test_watershed_ift07H  sL    

z%TestWatershedIft.test_watershed_ift07c                 C   s^   t ddgddggt j}t ddgddggt j}t||}ddgddgg}t|| d S )N   r   r	   )r   r   Zuint16r   r   rs  r   rt  r#   r#   r$   test_watershed_ift08g  s    

z%TestWatershedIft.test_watershed_ift08N)r1   r2   r3   rv  rw  r|  r}  r~  r  r  r  r#   r#   r#   r$   rf    s   rf  )os.pathr~   numpyr   Znumpy.testingr   r   r   r   r   r   Zpytestr   rx   Zscipy.ndimager    r
   r   r5   rI   rL   rM   rO   rQ   rT   rV   ra   rf   rg   rl   ro   rp   rq   ru   rv   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  r  r  r  r  r  r!  r#  r$  r%  r'  r(  r+  r2  r4  r5  r7  r:  r;  r=  r@  rA  rC  rD  rE  rG  rJ  rK  rL  rM  rQ  rR  rS  rT  rU  rX  rZ  r[  r\  r^  r`  ra  re  rf  r#   r#   r#   r$   <module>   s    E&								
	
			

		

