o
    hk                     @   sx   d dl Zd dlZd dlmZ d dlmZmZmZm	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 )    N)assert_array_equal)disk_level_setinverse_gaussian_gradientmorphological_chan_vese%morphological_geodesic_active_contourc                  C   s4   t jddddf } | d d}t | d S )N      r   
   )npmgridsumexp)coordssqrdistances r   f/var/www/html/scripts/venv/lib/python3.10/site-packages/skimage/segmentation/tests/test_morphsnakes.pygaussian_blob   s   r   c                  C   s   t d} t d}tt t| d|d W d    n1 s!w   Y  tt t| d|d W d    d S 1 s>w   Y  d S )N)r
   r
      )r
   	      num_iterinit_level_setr   zerospytestraises
ValueErrorr   r   imglsr   r   r   &test_morphsnakes_incorrect_image_shape      

"r"   c                  C   s   t d} t d}tt t| d|d W d    n1 s!w   Y  tt t| d|d W d    d S 1 s>w   Y  d S )N)   r$   r$   r$   r   r   r   r   r   r   r   test_morphsnakes_incorrect_ndim   r#   r%   c                  C   s   t d} t| jddd}t j| jt jd}t j| jt jd}t| d|d}t|| t| d|d}t|| t| d|dd	d
d}t|| |j	|j	  krZ|j	  krZt jks]J  J d S )N   r'      r)   r   centerradiusdtyper   r   r   r   )r   r   balloon	threshold	smoothing)
r   r   r   shapeint8onesr   r   r   r.   )r    r!   	ref_zerosref_onesacwe_lsgac_lsgac_ls2r   r   r   test_morphsnakes_black'   s"   



2r;   c                  C   sr   t  } t| jddd}t| jddd}t| d|d}t| d|d}t|| |j|j  kr4tjks7J  J d S )Nr(   r   r*   r   r
   r   )r   r   r3   r   r   r.   r   r4   )r    ls1ls2acwe_ls1acwe_ls2r   r   r   'test_morphsnakes_simple_shape_chan_vese>   s   
&r@   c                  C   s   t ddddt} t| ddd}t | jddd}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jd}t|d|dd}t	|| |j
tjksWJ d S )Nr&   r(   g      @r*   g      $@g      ?)alphasigmar   )r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   )r   r   r   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   )r   astypefloatr   r3   r   arrayr4   r   r   r.   )r    gimgr!   refr9   r   r   r   5test_morphsnakes_simple_shape_geodesic_active_contourK   s.   

rH   c                  C   s   t d} t| dd}t jg dg dg dg dg dg dgt jd}t| dd}t jg dg d	g d
g d
g d
g d	gt jd}t|| t|| d S )N)r   r   r   checkerboard)r   r   r   r   r   r   )r   r   r   r   r   r   r-   disk)r   r   r   r   r   r   )r   r   r   r   r   r   )r   r   r   r   r   r   )r   r   r   rE   r4   r   r   )imagecheckerboard_lscheckerboard_refdisk_lsdisk_refr   r   r   test_init_level_setsd   s.   



rP   c                     s~   t d} g   fdd}t| dd|d} d dksJ | dks%J t d d	  d
d  D ]
\}}||ks<J q2d S )N)   rQ   rQ   c                    s     |   d S )N)appendr   )x	evolutionr   r   callback   s   z%test_morphsnakes_3d.<locals>.callbackr)   rJ   )iter_callbackr   Q   r/   r   )r   r   r   r   zip)rK   rV   r!   v1v2r   rT   r   test_morphsnakes_3dz   s   
"r\   )numpyr   r   numpy.testingr   skimage.segmentationr   r   r   r   r   r"   r%   r;   r@   rH   rP   r\   r   r   r   r   <module>   s    