o
    ^i:H                     @   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/RAG/RAG_venv/lib/python3.10/site-packages/skimage/segmentation/tests/test_slic.pytest_color_2d      
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   
r3   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_2d5   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 )Nr4   r5   r   z&channel_axis=-1 indicates multichannel)matchr;   )r   r   pytestraises
ValueErrorr	   r&   r(   r(   r)    test_gray2d_default_channel_axisP   s   

rD   皙?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_labelsX   s   rK   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_offsetrG   rH   seg3seg4seg5n_seg1n_seg4r(   r(   r)   ,test_slic_consistency_across_image_magnitude^   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_3dz   s$   

rm   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   ra   r   gzo ?g$I$I?MbP?r   r   r-   Fr9   )r   r   r   r   r!   rb   rc   rd   r   arangere   r    r	   r   r"   r#   rf   )
r%   r&   rg   rh   ri   shadesrk   shr'   rl   r(   r(   r)   test_gray_3d   s4   



rs   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   rE   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   rE   r   r   r7   )r   r   r;   r:   r     r   )r   r   spacingr:   r;   r   )
r   r   r   rw   rx   ry   r    r!   r	   r   )r%   r&   result_non_spacedresult_spacedseg_non_spaced
seg_spacedr(   r(   r)   test_spacing   s&   
	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   rw   rx   r@   rA   rB   r	   rC   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 )Nru   r   r   r   r   r   r   r   r   r   r   r   r   r   -C6?TFr   )r:   r   r<   r   r;   皙?)r:   r   r<   max_size_factorr   r;   )r   rw   rx   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  r+   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   r4   r5   r   r6   r7   r8   r   r   r   Fr9   )r   r   r   r   r    r!   r	   allravelrp   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   r4   r   r   r   r   r   r   r/   r   F)r   r   r   mask   r   r   r1   r%   mskr&   r'   r(   r(   r)   test_color_2d_mask0  (   

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   r1   r   r(   r(   r)   test_multichannel_2d_maskK  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   r4   r   r   r   r5   r   r6   r7   r8   r   r/   r   Fr   r   r:   r;   r<   r   r   r   r1   r   r(   r(   r)   test_gray_2d_maskg  s8   


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   rt   ru   rE   r   r   r   r   r   r   r   r   rv   )r   r   r;   r   )r   r   r   r   rw   rx   r    r!   ry   r	   r   )r%   r   r&   rz   r{   r(   r(   r)   test_list_sigma_mask  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   rE   r   r   r7   )r   r   r;   r:   r   rv   )r   r   r   r:   r;   r   )r   r   r   r   rw   rx   ry   r    r!   r	   r   )r%   r   r&   r   r   r   r   r(   r(   r)   test_spacing_mask  s*   

	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   ru   r   r   r   r   TF)r:   r   r<   r   r;   r   )r:   r   r<   r   r   r;   r   )r   r   rw   rx   r	   r   )r   r&   r   r   r   r   r   r(   r(   r)   test_enforce_connectivity_mask  sP   
	

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   r4   r   r   r   r   r   r   r   r/   r   T)r   r   r   r   r   r   r   r1   r   r(   r(   r)   test_slic_zero_mask  r   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   r4   r   r   r   r5   r   r6   r7   r8   r   r/   r   Fr   )r   r   r   r   r    r!   r2   r	   rp   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 )Nrn   r   r   r   r   r`   r   r   r   r/   r-   )r   r   r   	   )r   r   r   r   r!   rb   rc   rd   r   re   r    r2   r	   r   r"   r#   rf   )
r   r%   r&   rg   rh   ri   rj   rk   rl   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 )Nrn   r   r   r   r   r-   ro   r   r/   F)r   r   r;   r<   r   r   )r   r   r   r   r!   rb   rc   rd   r   linspacere   r    r2   r	   r   r"   r#   rf   )r   r%   r&   rg   rh   ri   rq   rk   rr   r'   rl   r(   r(   r)   test_gray_3d_mask*  s*   



&r   dtype)float16float32float64uint8ry   c                 C   s&   t jdd| }t|dd d d S )N   r   )r   r;   )r   r   randrO   r	   )r   r&   r(   r(   r)   test_dtype_supportD  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   r:   r   max_num_iter)r   r   r   standard_normalr   gaussianr	   min)r%   r&   r   superpr(   r(   r)   test_start_label_fixL  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   rx   nanr@   rA   rB   r	   isnan)r&   r   r(   r(   r)   'test_raises_ValueError_if_input_has_NaNf  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   rx   r@   rA   rB   r	   isfinite)r   r&   r   r(   r(   r)   'test_raises_ValueError_if_input_has_infp  s   
r   )rE   )/	itertoolsr   numpyr   r@   numpy.testingr   skimager   r   r   skimage._shared.testingr   r   skimage.segmentationr	   r*   r3   r=   rD   rK   r_   rm   rs   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>   sL    

0"2

