o
    h                     @   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   x   g d}g d}ddg}dD ]+}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_measurements_statsr   selfxr   r   shpcountssums r%   `/var/www/html/scripts/venv/lib/python3.10/site-packages/scipy/ndimage/tests/test_measurements.pytest_a   s   
zTest_measurements_stats.test_ac                 C   r   )
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      g d}g d}ddg}dD ]4}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   centeredr   r   r         ?r   r    r!   r   r   r"   r#   r$   centersr%   r%   r&   test_a_centered,      z'Test_measurements_stats.test_a_centeredc                 C   r+   )Nr   r(   r   r)   r   Tr,   r   r   r   r.   r   r/   r%   r%   r&   test_b_centered9   r2   z'Test_measurements_stats.test_b_centeredc                 C   r+   )Nr   )        r4         "@r5   r4   r5   r   Tr,   r   r   r   r.   r   r/   r%   r%   r&   test_nonint_labelsF   r2   z*Test_measurements_stats.test_nonint_labelsN)	__name__
__module____qualname____doc__r'   r*   r1   r3   r6   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 ]\}}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)   )r4   r4         @r<   r4   r<   r   T)r   r   find_maxr   )r   r   find_minr   )r   r   r>   find_min_positions   i)r   r   r=   find_max_positions)	r   r   _selectr   lenr   r   dtypekind)r    r!   casesr   r   resultr%   r%   r&   
test_basicW   sJ   z#Test_measurements_select.test_basicN)r7   r8   r9   r:   rI   r%   r%   r%   r&   r;   T   s    r;   c                  C   0   t g } t| \}}t|d t|d d S Nr	   r   onesr   labelr   r   dataoutnr%   r%   r&   test_label01z      

rS   c                  C   rJ   )Nr   r   zerosr   rN   r   r   rO   r%   r%   r&   test_label02   rT   rW   c                  C   s4   t dg} t| \}}t|dg t|d d S rK   rL   rO   r%   r%   r&   test_label03      rX   c                  C   s4   t dg} t| \}}t|dg t|d d S Nr	   r   rU   rO   r%   r%   r&   test_label04   rY   r[   c                  C   s6   t dg} t| \}}t|g d t|d d S )N   )r	   r	   r	   r	   r	   r	   rL   rO   r%   r%   r&   test_label05   s   r]   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   r@   r@   r   r   r   rN   r   r   rO   r%   r%   r&   test_label06   s   r`   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_   rO   r%   r%   r&   test_label07       

rb   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   r@   r@   r   r   r   r   r   r   r   r   r   r   r   r_   rO   r%   r%   r&   test_label08   rc   rn   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 )Nre   rf   rg   rh   ri   r   rj   rk   )r   r   r   r   r   r   )r   r   r   r@   r@   r   r@   r   r   r   generate_binary_structurerN   r   r   rP   structrQ   rR   r%   r%   r&   test_label09   s"   

rs   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 )Nra   )r   r	   r	   r   r	   r   )r   r	   r	   r	   r	   r   r   r	   ro   rq   r%   r%   r&   test_label10   s   

rt   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 rd   r
   r   r   r   rN   r   r   typerP   rQ   rR   expectedr%   r%   r&   test_label11   *   

ry   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 )Nre   rf   rg   rh   ri   outputrj   rk   rl   rm   r   ru   )rw   rP   rR   rx   r%   r%   r&   test_label11_inplace   rz   r}   c               	   C   s|   t D ]9} 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	   ri   r	   ru   rv   r%   r%   r&   test_label12  s&   

r~   c                  C   sp   t D ]3} 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	   ru   rv   r%   r%   r&   test_label13  s"   

r   c                  C   sL   t dg} tD ]}t jdg|d}tj| |d}t|d t|d qd S )Nr\   rE   r{   r	   )r   rM   r
   rV   r   rN   r   r   rP   tr|   rR   r%   r%   r&   test_label_output_typed-  s   
r   c                  C   sD   t dg} tD ]}tj| |d\}}t|d |j|ksJ qd S )Nr\   r{   r	   )r   rM   r
   r   rN   r   rE   r   r%   r%   r&   test_label_output_dtype6  s   
r   c                  C   s>   t dg} tD ]}t dg|}tttftj| |d qd S )Nr\   
   r{   )	r   rM   r
   rV   assert_raisesRuntimeError
ValueErrorr   rN   )rP   r   r|   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 ]:}| |d d d d f }t|j	d D ]%}||d d d d f }t
t||d ||d d d d f  |d7 }qYqEd S )	NrP   zlabel_inputs.txtzlabel_strels.txtzlabel_results.txt)   r   )r   r@   r@   r   r	   )r   loadtxtospathjoindirname__file__r   rangeshaper   r   rN   )rP   strelsresultsrrA   djsr%   r%   r&   test_label_structuring_elementsF  s.   




(
r   c                  C   sL   ddd} t t jt dkr$d}t jt ||}| | d S d S )Nffffff?r   c                 S   s8   | |k}t |j}t|t||\}}t|}d S )N)rD   r   r   rN   rp   find_objects)imgthreshsizemaskranklaco_r%   r%   r&   SEZ  s   

ztest_ticket_742.<locals>.SErA   )r@     r   )r   r   )r   rE   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s4J dS )
z/Github issue #3025 - improper merging of labels)<   i@  r	   N  i  )$   r   )#   i  )r   i  )r@   r@   )r   rV   r   rN   rM   )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 )Nr   r.   )
r   r   r   r   rN   r   rE   int32int64r   )
test_arrayrN   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 )Nr   r%   )r   rM   intr   r   r   rP   rQ   r%   r%   r&   test_find_objects01{  s   
r   c                  C   s(   t jg td} t| }t|g k d S )Nr   )r   rV   r   r   r   r   r   r%   r%   r&   test_find_objects02  s   
r   c                  C   4   t jdgtd} t| }t|tddd fg d S )Nr	   r   r   r   rM   r   r   r   r   slicer   r%   r%   r&   test_find_objects03     
r   c                  C   s(   t jdgtd} t| }t|g  d S )Nr	   r   )r   rV   r   r   r   r   r   r%   r%   r&   test_find_objects04     
r   c                  C   r   )Nr\   r   r   r   r   r%   r%   r&   test_find_objects05  r   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 )Nr^   r   r	   r   r   r\   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 )Nra   )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 )Nre   rj   rk   rl   rm   r   r	   r@   r   r\   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 )Nre   rj   rk   ra   rm   r   r	   r@   r   r\   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 entriesre   rj   rk   ra   rm   r   )ignore_value)r	   r   r   N)r   r   r   value_indicesr   listkeyswhere)rP   vi	true_keystruevikr%   r%   r&   test_value_indices01  s   
r   c                  C   sP   t jdt jd} d}tt|d t|  W d   dS 1 s!w   Y  dS )zTest input checking)r\   r   r   z(Parameter 'arr' must be an integer array)matchN)r   rV   float32r   r   r   r   )rP   msgr%   r%   r&   test_value_indices02  s
   "r   c                  C   s   dD ]B} t jddg ddg  ddg  t jd| }t |}t|}tt|	 t| |D ]}t 
||k}t|| | q3qdS )z2Test different input array shapes, from 1-D to 4-D))r   )   r   )r@   r@   r   )r@   r@   r   r      r	   r   r@   r   N)r   r   r   r   uniquer   r   r   r   r   r   )r   r   trueKeysr   r   trueNdxr%   r%   r&   test_value_indices03  s   0

r   c                  C   .   t D ]} tg | }t|}t|d qd S )Nr4   )r
   r   r   r   sumr   rw   inputr|   r%   r%   r&   
test_sum01  
   
r   c                  C   2   t D ]} tddg| }t|}t|d qd S )Nr   r   r4   )r
   r   rV   r   r   r   r   r%   r%   r&   
test_sum02  
   
r   c                  C   r   )Nr   )r
   r   rM   r   r   r   r   r%   r%   r&   
test_sum03  r   r   c                  C   r   )Nr	   r         @r
   r   r   r   r   r   r   r%   r%   r&   
test_sum04  r   r   c                  C   s:   t D ]} tddgddgg| }t|}t|d qd S )Nr	   r   r@   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   r4   r   r   boolr
   r   r   r   r   rw   r   r|   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   r4   )r   rM   r   r
   rV   r   r   r   r   r%   r%   r&   
test_sum07     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#  r   r   c                  C   N   t ddgt} tD ]}t ddgddgg|}tj|| d}t|d q
d S )Nr	   r   r   r@   r   r         @)r   r   r   r
   r   r   r   r   r%   r%   r&   
test_sum09+     r   c                  C   D   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%   r%   r&   
test_sum103     r   c                  C   sR   t ddgt j} tD ]}t ddgddgg|}tj|| dd}t|d qd S )Nr	   r   r@   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   r@   r      r   r   )r   r4   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 ]1}t ddgddgg|}tj|| g dd}tj|| g dd}||k s9J t|g d qd S r   )	r   r   r   r
   r   r   
sum_labelsallr   )r   rw   r   
output_sumoutput_labelsr%   r%   r&   test_sum_labelsK  s   
r  c                  C   r   r   )r   r   r   r
   r   meanr   r   r%   r%   r&   test_mean01W  r   r	  c                  C   r   Nr	   r   r   r@   r   r   r   )r   r   r   r   r  r   r   r%   r%   r&   test_mean02_  r   r  c                  C   N   t ddg} tD ]}t ddgddgg|}tj|| dd}t|d q	d S )Nr	   r   r@   r   r   r   )r   r   r
   r   r  r   r   r%   r%   r&   test_mean03f     r  c                  C   s   t ddgddggt j} t jdd8 tD ],}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    d S 1 sNw   Y  d S )Nr	   r   r   ignorer  r@   r   r   r   r         @)
r   r   r   errstater
   r   r  r   r   isnanr   r%   r%   r&   test_mean04o  s   "r  c                  C   r   r
  )r   r   r   r
   r   minimumr   r   r%   r%   r&   test_minimum01z  r   r  c                  C   D   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  r   r  c                  C   r  )Nr	   r   r@   r   r   r   )r   r   r
   r   r  r   r   r%   r%   r&   test_minimum03  r  r  c                  C   ^   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   r@   r   r   r@   r  r   )r   r   r4   )r   r   r
   r   r  r   r   r%   r%   r&   test_minimum04     r  c                  C   r   )Nr	   r   r   r@   r   r   r   )r   r   r   r
   r   maximumr   r   r%   r%   r&   test_maximum01  r   r   c                  C   r  r  )r   r   r   r   r  r   r   r%   r%   r&   test_maximum02  r   r!  c                  C   r  )Nr	   r   r@   r   r   r   )r   r   r
   r   r  r   r   r%   r%   r&   test_maximum03  r  r"  c                  C   r  )Nr	   r   r@   r   r  r   )r   r   r4   )r   r   r
   r   r  r   r   r%   r%   r&   test_maximum04  r  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	   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   medianr   r   r   r|   r%   r%   r&   test_median01  s   

r0  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   r|   r%   r%   r&   test_median02  s   

r3  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   r1  r/  r%   r%   r&   test_median03  s   

r4  c                  C   s:   t jddgtd} tj| t ddgd}t|dg d S )Nr	   r   r   r   r   )r   asarrayr   r   r.  rM   r   r2  r%   r%   r&   test_median_gh12836_bool  s   r7  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   r   r5  r	   r   g     P@)r   r6  r   r   r.  rM   r   r2  r%   r%   r&   test_median_no_int_overflow  s   r:  c               
   C      t jdd9 tD ]-} t g | }t }|td t|}W d    n1 s*w   Y  t	t 
| q	W d    d S 1 sBw   Y  d S Nr  r  zMean of empty slice)r   r  r
   r   r   filterRuntimeWarningr   variancer   r  rw   r   supr|   r%   r%   r&   test_variance01     "rB  c                  C   0   t D ]} tdg| }t|}t|d qd S Nr	   r4   r
   r   r   r   r?  r   r   r%   r%   r&   test_variance02  
   
rG  c                  C   r   Nr	   r@   r   rF  r   r%   r%   r&   test_variance03  r   rJ  c                  C   (   t ddgt} t| }t|d d S )Nr	   r   g      ?)r   r   r   r   r?  r   r   r|   r%   r%   r&   test_variance04  r   rM  c                  C   >   g d} t D ]}tg d|}t|| d}t|d qd S N)r   r   r@   )r	   r@   r  r   r   rF  r   r%   r%   r&   test_variance05     rP  c                  C   t   g d} t jdd& tD ]}t g d|}t|| g d}t|g d qW d    d S 1 s3w   Y  d S N)r   r   r@   r@   r   r  r  )r	   r@   r  r   r  )r   r@   r   )r   r   r4   )r   r  r
   r   r   r?  r   r   r%   r%   r&   test_variance06     "rT  c               
   C   r;  r<  )r   r  r
   r   r   r=  r>  r   standard_deviationr   r  r@  r%   r%   r&   test_standard_deviation01"  rC  rW  c                  C   rD  rE  r
   r   r   r   rV  r   r   r%   r%   r&   test_standard_deviation02,  rH  rY  c                  C   s8   t D ]} tddg| }t|}t|td qd S rI  )r
   r   r   r   rV  r   sqrtr   r%   r%   r&   test_standard_deviation033  s
   
r[  c                  C   rK  )Nr	   r   r.   )r   r   r   r   rV  r   rL  r%   r%   r&   test_standard_deviation04:  r   r\  c                  C   rN  rO  rX  r   r%   r%   r&   test_standard_deviation05@  rQ  r]  c                  C   rR  rS  r   r  r
   r   r   rV  r   r   r%   r%   r&   test_standard_deviation06H  rU  r_  c                  C   sl   dg} t jdd# tD ]}t dg|}t|| dg}t|dg qW d    d S 1 s/w   Y  d S )Nr	   r  r  gl `yr   r^  r   r%   r%   r&   test_standard_deviation07Q  s   "r`  c                  C   r   )Nr	   r   r   r@   r   r   r   r   )r   r   r   r
   r   minimum_positionr   r   r%   r%   r&   test_minimum_position01Z  r   rc  c                  C   @   t D ]} t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   rb  r   r   r%   r%   r&   test_minimum_position02b     

rk  c                  C   6   t g dg dg dgt} t| }t|d d S re  r   r   r   r   rb  r   rL  r%   r%   r&   test_minimum_position03k     

ro  c                  C   rm  )Nrf  )r@   r   r	   r   rh  ra  rn  rL  r%   r%   r&   test_minimum_position04s  rp  rq  c                  C   J   g d} t D ]}tg dg dg dg|}t|| }t|d qd S )Nr	   r   r   r   rf  rg  )r	   r\   r   r@   )r   r   rj  r   r%   r%   r&   test_minimum_position05{     
rt  c                  C   L   g d} t D ]}tg dg dg dg|}t|| d}t|d qd S )Nr	   r   r@   r   rf  rg  rh  r   r   r	   rj  r   r%   r%   r&   test_minimum_position06     
ry  c                  C   sb   g d} t D ](}tg dg dg dg|}t|| ddg}t|d d t|d	 d
 qd S )Nrw  rf  rg  rh  r   r@   r   rx  r	   ri  rj  r   r%   r%   r&   test_minimum_position07     
r{  c                  C   r   )Nr	   r   r   r@   r   r   )r	   r   )r   r   r   r
   r   maximum_positionr   r   r%   r%   r&   test_maximum_position01  s   r~  c                  C   rd  )Nrf  r@   r   r  r   rh  ri  r
   r   r   r   r}  r   r   r%   r%   r&   test_maximum_position02  rl  r  c                  C   rm  )Nrf  r  rh  ra  )r   r   r   r   r}  r   rL  r%   r%   r&   test_maximum_position03  rp  r  c                  C   rr  )Nrs  rf  r  rh  r	   r	   r  r   r%   r%   r&   test_maximum_position04  ru  r  c                  C   rv  )Nrs  rf  r  rh  r	   ra  r  r   r%   r%   r&   test_maximum_position05  rz  r  c                  C   sb   g d} t D ](}tg dg dg dg|}t|| ddg}t|d d t|d d	 qd S )
Nrs  rf  r  rh  r	   r   r   ra  r  r  r   r%   r%   r&   test_maximum_position06  r|  r  c                  C   sh   t g d} tD ](}t g dg dg dg|}t|| ddg}t|d d t|d	 d
 q	d S )N)r   r  r4         @rf  r  rh  r   r  r   ra  r	   )r   r@   )r   r   r
   r   r}  r   r   r%   r%   r&   test_maximum_position07  s   
r  c                  C   s   t ddgt} tD ]:}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   r@   r   r   )r   r   r   r
   r   extremar  r  rb  r}  r   r   rw   r   output1output2output3output4output5r%   r%   r&   test_extrema01  s   r  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   r@   r   r   )
r   r   r
   r   r  r  r  rb  r}  r   r  r%   r%   r&   test_extrema02  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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   r@   r   r  r   r   )
r   r   r
   r   r  r  r  rb  r}  r   r  r%   r%   r&   test_extrema03  s.   r  c                  C   s   g d} t D ]Z}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 )	Nrs  rf  r  rh  r	   r   r   r@   )
r
   r   r   r   r  r  r  rb  r}  r   r  r%   r%   r&   test_extrema04  s,   
r  c                  C   sB   ddg} t D ]}tddgddgg|}t|}t||  qd S )Nr4   r	   r   r
   r   r   r   center_of_massr   rx   rw   r   r|   r%   r%   r&   test_center_of_mass01(     
r  c                  C   sB   ddg} t D ]}tddgddgg|}t|}t||  qd S rZ   r  r  r%   r%   r&   test_center_of_mass020  r  r  c                  C   sB   ddg} t D ]}tddgddgg|}t|}t||  qd S )Nr   r	   r  r  r%   r%   r&   test_center_of_mass038  r  r  c                  C   sB   ddg} t D ]}tddgddgg|}t|}t||  qd S rZ   r  r  r%   r%   r&   test_center_of_mass04@  r  r  c                  C   sB   ddg} t D ]}tddgddgg|}t|}t||  qd S )Nr.   r	   r  r  r%   r%   r&   test_center_of_mass05H  r  r  c                  C   s8   ddg} t ddgddggt}t|}t||  d S )Nr.   r	   r   r@   r   r   r   r   r  r   rx   r   r|   r%   r%   r&   test_center_of_mass06P  s   
r  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.   r4   r   r@   r  r   rx   r   r|   r%   r%   r&   test_center_of_mass07W  s
   r  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   r\   r@   r  r  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.   r4   )r.   r   r  r  r%   r%   r&   test_center_of_mass09g  s
   r  c                  C   s2   t d} t d}t|ddd}t||  d S )Nr   r   )r   rM   aranger   	histogramr   r  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	   r@   r   r@   r@   r@   r@   r   r   r\   r	   r   r   r   r  r   r  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   r@   r   )r	   r	   r@   r   r@   r\   r@   r@   r   r   r\   ri  r	   r  )r   	expected1	expected2r   r|   r%   r%   r&   test_histogram03~  s   r  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)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?  rV  r   rZ  r.  r  r  )r   lblr  varstdmedminmaxr%   r%   r&   test_stat_funcs_2d  s   r  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    rP   markersrQ   rx   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 )	Nr  r  r  r  r  r  )r   r   r	   r	   r	   r   r   r  r  r  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 )	Nr  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&   test_watershed_ift03  s<   

z%TestWatershedIft.test_watershed_ift03c              
   C      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 )
Nr  r  r  r  r  r  r  r  r  r  r  r%   r%   r&   test_watershed_ift04  F   

z%TestWatershedIft.test_watershed_ift04c              
   C   r  )
Nr  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&   test_watershed_ift05  r  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 )
Nr  r  r  r  r  r  r  r  r  r  r  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   r   r  r  r  .r  r  r  )r  r|   r  r  )
r   rV   r  	transposer   r   int16r   r  r   )r    r   rP   r  rQ   rx   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   uint16r   r   r  r   r  r%   r%   r&   test_watershed_ift08g  s   

z%TestWatershedIft.test_watershed_ift08N)r7   r8   r9   r  r  r  r  r  r  r  r  r%   r%   r%   r&   r    s    r  )os.pathr   numpyr   numpy.testingr   r   r   r   r   r   pytestr   r   scipy.ndimager    r
   r   r;   rS   rW   rX   r[   r]   r`   rb   rn   rs   rt   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   r   r   r   r   r   r   r   r   r   r  r  r	  r  r  r  r  r  r  r  r   r!  r"  r#  r&  r0  r3  r4  r7  r:  rB  rG  rJ  rM  rP  rT  rW  rY  r[  r\  r]  r_  r`  rc  rk  ro  rq  rt  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  r  r  r%   r%   r%   r&   <module>   s     E&								
	
			

		

