o
    eL                     @   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(   TD:\Projects\ConvertPro\env\Lib\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seg2Zallowed_mismatch_ratior   Zndiffr(   r(   r)   _check_segment_labelsT   s   rH   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   Zuint16r   maxr	   testingZassert_array_equalabs)Z	img_uint8Z
img_uint16Zimg_float32Zimg_float32_normZimg_float32_offsetrF   rG   Zseg3Zseg4Zseg5Zn_seg1Zn_seg4r(   r(   r)   ,test_slic_consistency_across_image_magnitudeZ   s    rP   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$   

r^   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   rR   r   gzo ?g$I$I?MbP?r   r   r,   Fr8   )r   r   r   r   r!   rS   rT   rU   r   arangerV   r    r	   r   r"   r#   rW   )
r%   r&   rX   rY   rZ   shadesr\   shr'   r]   r(   r(   r)   test_gray_3d   s(   



rd   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$   

rm   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   rh   ri   rj   r    r!   r	   r   )r%   r&   result_non_spacedresult_spacedseg_non_spaced
seg_spacedr(   r(   r)   test_spacing   s2   




rv   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 )Nrn   ro   r   r   Tr   )r:   r;   r   )r   rh   ri   r?   r@   rA   r	   rB   r(   r(   r)   test_invalid_lab_conversion   s   
"rw   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 )Nrf   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   rh   ri   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   rp   Fr8   )r   r   r   r   r    r!   r	   allravelra   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   re   rf   rD   r   r   r   r   r   r   r   r   rg   )r   r   r:   r   )r   r   r   r   rh   ri   r    r!   rj   r	   r   )r%   r   r&   rk   rl   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   rn   ro   )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   rg   )r   r   rq   r9   r:   r   )r   r   r   r   rh   ri   rj   r    r!   r	   r   )r%   r   r&   rr   rs   rt   ru   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   rf   rx   ry   r   rz   TF)r9   r   r;   r   r:   r{   )r9   r   r;   r|   r   r:   r   )r   r   rh   ri   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.   rp   Fr   )r   r   r   r   r    r!   r1   r	   ra   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 )Nr_   r   r   r   r   rQ   r   r   r   r.   r,   )r   r   r   	   )r   r   r   r   r!   rS   rT   rU   r   rV   r    r1   r	   r   r"   r#   rW   )
r   r%   r&   rX   rY   rZ   r[   r\   r]   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 )Nr_   r   r   r   r   r,   r`   r   r.   F)r   r   r:   r;   r   r   )r   r   r   r   r!   rS   rT   rU   r   ZlinspacerV   r    r1   r	   r   r"   r#   rW   )r   r%   r&   rX   rY   rZ   rb   r\   rc   r'   r]   r(   r(   r)   test_gray_3d_mask  s*   




&r   dtype)Zfloat16Zfloat32Zfloat64Zuint8rj   c                 C   s&   t jdd| }t|dd d d S )N   r   )r   r:   )r   r   ZrandrL   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   Zmax_num_iter)r   r   r   Zstandard_normalr   Zgaussianr	   min)r%   r&   r   Z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   ri   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   ri   r?   r@   rA   r	   isfinite)r   r&   r   r(   r(   r)   'test_raises_ValueError_if_input_has_inf0  s   
r   )rD   )/	itertoolsr   numpyr   r?   Znumpy.testingr   Zskimager   r   r   Zskimage._shared.testingr   r   Zskimage.segmentationr	   r*   r2   r<   rC   rH   rP   r^   rd   rm   rv   rw   r   r   r   r   r   r   r   r   r   r   r   r   r   markZparametrizer   r   r   r   r   r(   r(   r(   r)   <module>   sP    


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