o
    hJ                     @   sn  d dl mZ d dlZd dlZd dlmZ d dlmZm	Z	m
Z
 d dlmZmZ d dlmZ e dd Zd	d
 Zdd Zdd ZdCd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'ej()d8g d9d:d; Z*d<d= Z+d>d? Z,ej()d@ej- ej-gdAdB Z.dS )D    )productN)assert_equal)datafiltersimg_as_float)run_in_parallelexpected_warnings)slicc                  C   .  t jd} t d}d|d dd ddf< d|dd d ddf< d|dd dd df< |d| j|jd 7 }d||dk< d||dk < t|ddd	dd
}ttt 	|d t|j|jd d  t|d dd df d t|dd d df d t|d ddd f d t|dd dd f d d S )Nr               
      {Gz?size   F)
n_segmentssigmaenforce_connectivitystart_labelr   
nprandomdefault_rngzerosnormalshaper	   r   lenuniquerngimgseg r(   _/var/www/html/scripts/venv/lib/python3.10/site-packages/skimage/segmentation/tests/test_slic.pytest_color_2d   s"   

r*   c                  C   sL  t jd} t d}d|d dd dddf< d|d ddd ddf< d|dd d dddf< d|dd dd ddf< |d	| j|jd
 7 }t j|dd|d}t|dddd}tt	t 
|d t|j|jd d  t|d dd df d t|dd d df d t|d ddd f d t|dd dd f d d S )Nr   r   r      r   r   r   r      r,   r   r   outF)r   r   r   r   r   r   r   r   r   r    r!   clipr	   r   r"   r#   r$   r(   r(   r)   test_multichannel_2d"   s   
r2   c               	   C   s$  t jd} t d}d|d dd df< d|dd d df< d|dd dd f< |d| j|jd 7 }d	||d	k< d||dk < t|dd
d	d ddd}ttt 	|d
 t|j|j t|d dd df d t|dd d df d t|d ddd f d	 t|dd dd f d d S )Nr   r   r   Q?r   q=
ףp?      ?F%uk?r   r   r   Fr   r   compactnesschannel_axisconvert2labr   r   r   r   r$   r(   r(   r)   test_gray_2d6   s"   

r<   c                  C   sd   t d} d| d dd df< tjtdd t|  W d    n1 s%w   Y  t| d d d S )Nr3   r4   r   z&channel_axis=-1 indicates multichannel)matchr:   )r   r   pytestraises
ValueErrorr	   r&   r(   r(   r)    test_gray2d_default_channel_axisJ   s   

rC   皙?c                 C   s(   | j }t| |k}|| |k sJ d S )N)r   r   sum)seg1seg2allowed_mismatch_ratior   ndiffr(   r(   r)   _check_segment_labelsT   s   rJ   c                  C   s   t  d dd df } d| tj }t| }||  }|d }t| }t|}t|}t|}t|}	tj	|| tj	|| tj	||	 | }
| }t
|
| |
 dk s_J d S )N      i  g      ?)r   catastyper   uint16r   maxr	   testingassert_array_equalabs)	img_uint8
img_uint16img_float32img_float32_normimg_float32_offsetrF   rG   seg3seg4seg5n_seg1n_seg4r(   r(   r)   ,test_slic_consistency_across_image_magnitudeZ   s    r^   c            	      C   s
  t jd} t d}g }|jd d D ]}|d }|td |t|d f qtt| }ttddgfd  }t	||D ]\}}|||< q>|d| j
|jd 7 }d||dk< d||dk < t|dd	dd
}ttt |d	 t	|td	D ]\}}t|| | qwd S )Nr   r   r      r   r   r   r   r   r   r   r,   )r   r   r   )r   r   r   r   r!   appendslicelistr   zipr    r	   r   r"   r#   range)	r%   r&   slicesdim_sizemidpointcolorsscr'   r(   r(   r)   test_color_3dv   s$   

rl   c            
   	   C   s   t jd} t d}g }|jD ]}|d }|td |t|d f qtt| }t 	ddd}t
||D ]\}}|||< q6|d| j|jd 7 }d||dk< d||dk < t|dd	dd d
dd}ttt |d	 t
|td	D ]\}}	t|| |	 qrd S )Nr   r   r   r`   r   gzo ?g$I$I?MbP?r   r   r,   Fr8   )r   r   r   r   r!   ra   rb   rc   r   arangerd   r    r	   r   r"   r#   re   )
r%   r&   rf   rg   rh   shadesrj   shr'   rk   r(   r(   r)   test_gray_3d   s(   



rr   c                  C   s   t jd} t g dg dgt}|d| j|jd 7 }t g dg dgt}tdg t	|dg dd dd	}W d    n1 sDw   Y  t
|| d S )
Nr   r   r   r   r   r   r   r   r   r   r   r   r   rD   r   z5Input image is 2D: sigma number of elements must be 2r   )r   2   r   )r   r   r:   r   )r   r   r   arrayfloatr    r!   intr   r	   r   )r%   r&   result_sigma	seg_sigmar(   r(   r)   test_list_sigma   s$   

r{   c               	   C   s   t jd} t g dg dgt}t g dg dgt}t g dg dgt}|d| j|jd	 7 }t|d
dd ddd}t|d
dddgdd dd}t	|| 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   rD   r   r   r6   )r   r   r:   r9   r     r   )r   r   spacingr9   r:   r   )
r   r   r   rv   rw   rx   r    r!   r	   r   )r%   r&   result_non_spacedresult_spacedseg_non_spaced
seg_spacedr(   r(   r)   test_spacing   s2   




r   c                  C   s\   t g dg dgtd } tt t| dddd W d    d S 1 s'w   Y  d S )Nr|   r}   r   r   Tr   )r:   r;   r   )r   rv   rw   r?   r@   rA   r	   rB   r(   r(   r)   test_invalid_lab_conversion   s   
"r   c               
   C   s   t g dg dg dgt} t| dddddd d	}t| dddddd d	}t| ddddd
dd d}t g dg dg dgt}t g dg dg dgt}t|| t|| t|| d S )Nrt   r   r   r   r   r   r   r   r   r   r   r   r   r   -C6?TFr   )r9   r   r;   r   r:   皙?)r9   r   r;   max_size_factorr   r:   )r   rv   rw   r	   r   )r&   segments_connectedsegments_disconnectedsegments_connected_low_maxresult_connectedresult_disconnectedr(   r(   r)   test_enforce_connectivity   sL   




r   c                  C   r
   )Nr   r   r   r   r   r   r   r   T)r   r   	slic_zeror   r   r   r   r$   r(   r(   r)   test_slic_zero   s   
r   c               	   C   s   t jd} t d}d|d dd df< d|dd d df< d|dd dd f< |d| j|jd 7 }d	||d	k< d||dk < t|dd
d	d ddd}t | t 	|j
ksZJ d S )Nr   r3   r4   r   r5   r6   r7   r   r   r~   Fr8   )r   r   r   r   r    r!   r	   allravelro   r   r$   r(   r(   r)   test_more_segments_than_pixels   s   

"r   c                  C   s  t jd} t d}d|ddddf< t d}d|d dd ddf< d|dd d ddf< d|dd dd df< |d| j|jd	 7 }t j|dd|d
 t|ddd|d}tt	t 
|d t|j|jd d  t|ddddf d t|ddddf d t|ddddf d t|ddddf d t|d dd d f d t|dd d d f d t|d d d df d t|d d dd f d d S )Nr   r3   r   r   r   r   r   r   r.   r   F)r   r   r   mask   r   r   r0   r%   mskr&   r'   r(   r(   r)   test_color_2d_mask  s,   


r   c                  C   s  t jd} t d}d|ddddf< t d}d|d dd dddf< d|d ddd ddf< d|dd d ddd	f< d|dd dd d	d
f< |d| j|jd 7 }t j|dd|d t|dd|d}tt	t 
|d t|j|jd d  t|ddddf d t|ddddf d t|ddddf d t|ddddf d t|d dd d f d t|dd d d f d t|d d d df d t|d d dd f d d S )Nr   )r   r   r   r   r   r+   r   r   r-   r,   r   r   r.   F)r   r   r   r   r   r   r0   r   r(   r(   r)   test_multichannel_2d_mask*  s.   

r   c               	   C   s  t jd} t d}d|ddddf< t d}d|d dd df< d|dd d df< d	|dd dd f< |d
| j|jd 7 }t j|dd|d t|dddd d|d}tt	t 
|d t|j|j t|ddddf d t|ddddf d t|ddddf d t|ddddf d t|d dd d f d t|dd d d f d t|d d d df d t|d d dd f d d S )Nr   r3   r   r   r   r4   r   r5   r6   r7   r   r.   r   Fr   r   r9   r:   r;   r   r   r   r0   r   r(   r(   r)   test_gray_2d_maskG  s,   


r   c                  C   s   t jd} t d}d|d d ddf< t g dg dgt}|d| j|jd 7 }t g d	g d	gt}t	|d
ddgd |d}t
|| d S )Nr   )r   r-   r   r   rs   rt   rD   r   r   r   r   r   r   r   r   ru   )r   r   r:   r   )r   r   r   r   rv   rw   r    r!   rx   r	   r   )r%   r   r&   ry   rz   r(   r(   r)   test_list_sigma_maskb  s$   


r   c               	   C   s   t jd} t d}d|d d ddf< t g dg dgt}t g dg dgt}t g d	g d
gt}|d| j|jd 7 }t	|ddd d|d}t	|ddddgdd |d}t
|| 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   rD   r   r   r6   )r   r   r:   r9   r   ru   )r   r   r   r9   r:   r   )r   r   r   r   rv   rw   rx   r    r!   r	   r   )r%   r   r&   r   r   r   r   r(   r(   r)   test_spacing_maskp  s6   





r   c               
   C   s   t d} d| d d ddf< t g dg dg dgt}t|ddd	d
| d d}t|ddd
d
| d d}t|ddd	d
d| d d}t g dg dg dgt}t g dg dg dgt}t|| t|| t|| d S )N)r   r-   r   r   rt   r   r   r   r   TF)r9   r   r;   r   r:   r   )r9   r   r;   r   r   r:   r   )r   r   rv   rw   r	   r   )r   r&   r   r   r   r   r   r(   r(   r)   test_enforce_connectivity_mask  sJ   





r   c                  C   s  t jd} t d}d|ddddf< t d}d|d dd ddf< d|dd d ddf< d|dd dd df< |d| j|jd	 7 }t j|dd|d
 t|ddd|d}tt	t 
|d t|j|jd d  t|ddddf d t|ddddf d t|ddddf d t|ddddf d t|d dd d f d t|dd d d f d t|d d d df d t|d d dd f d d S )Nr   r3   r   r   r   r   r   r   r   r.   r   T)r   r   r   r   r   r   r   r0   r   r(   r(   r)   test_slic_zero_mask  s,   


r   c               	   C   s   t jd} t d}d|ddddf< t d}d|d dd df< d|dd d df< d	|dd dd f< |d
| j|jd 7 }t j|dd|d t|dddd d|d}t |ddddf j	d }t 
|ddddf  |kszJ d S )Nr   r3   r   r   r   r4   r   r5   r6   r7   r   r.   r~   Fr   )r   r   r   r   r    r!   r1   r	   ro   r   r   r   )r%   r   r&   r'   expectedr(   r(   r)   #test_more_segments_than_pixels_mask  s   


 *r   c            
      C   s8  t d} d| ddddddf< t jd}t d}g }| jD ]}|d }|td |t|d f q"tt| }ttddgfd  }t	||D ]\}}|||< qL|d|j
|jd	 7 }t j|dd|d
 t|dd| d}	ttt |	d t	|tddD ]\}}t|	| ddddddf | qd S )Nrm   r   r   r   r   r_   r   r   r   r.   r,   )r   r   r   	   )r   r   r   r   r!   ra   rb   rc   r   rd   r    r1   r	   r   r"   r#   re   )
r   r%   r&   rf   rg   rh   ri   rj   rk   r'   r(   r(   r)   test_color_3d_mask  s&   



&r   c                  C   s4  t d} d| ddddddf< t jd}t d}g }|jD ]}|d }|td |t|d f q"tt| }t 	ddd}t
||D ]\}}|||< qH|d|j|jd 7 }t j|dd|d	 t|ddd d
| d}	ttt |	d t
|tddD ]\}}
t|	| ddddddf |
 qd S )Nrm   r   r   r   r   r,   rn   r   r.   F)r   r   r:   r;   r   r   )r   r   r   r   r!   ra   rb   rc   r   linspacerd   r    r1   r	   r   r"   r#   re   )r   r%   r&   rf   rg   rh   rp   rj   rq   r'   rk   r(   r(   r)   test_gray_3d_mask  s*   




&r   dtype)float16float32float64uint8rx   c                 C   s&   t jdd| }t|dd d d S )N   r   )r   r:   )r   r   randrN   r	   )r   r&   r(   r(   r)   test_dtype_support  s   r   c               	   C   sV   t jd} | ddk}tj|dd}d}t||dddd	d
d}| |ks)J dS )a  Tests the fix for a bug producing a label < start_label (gh-6240).

    For the v0.19.1 release, the `img` and `slic` call as below result in two
    non-contiguous regions with value 0 despite `start_label=1`. We verify that
    the minimum label is now `start_label` as expected.
    r   )r,      r   r   )r   Nr-   r   Tr   )r   r:   r   r9   r   max_num_iter)r   r   r   standard_normalr   gaussianr	   min)r%   r&   r   superpr(   r(   r)   test_start_label_fix  s   	r   c                  C   sl   t jdtd} t j| d< tt t| d d W d    n1 s"w   Y  t |  }t| |d d d S N)r   r   )r   )r   r   r>   )r   r:   )	r   r   rw   NaNr?   r@   rA   r	   isnan)r&   r   r(   r(   r)   'test_raises_ValueError_if_input_has_NaN&  s   
r   infc                 C   sh   t jdtd}| |d< tt t|d d W d    n1 s!w   Y  t |}t||d d d S r   )r   r   rw   r?   r@   rA   r	   isfinite)r   r&   r   r(   r(   r)   'test_raises_ValueError_if_input_has_inf0  s   
r   )rD   )/	itertoolsr   numpyr   r?   numpy.testingr   skimager   r   r   skimage._shared.testingr   r   skimage.segmentationr	   r*   r2   r<   rC   rJ   r^   rl   rr   r{   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   markparametrizer   r   r   r   r   r(   r(   r(   r)   <module>   sP    


$#

