o
    :ήc                     @   s   d dl mZ d dlmZmZ d dlZd dlmZ d dl	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d Zdd Zdd Zdd Zdd Zd d! ZdS )"    )testing)assert_array_equalassert_allcloseN)camera)random_noiseimg_as_floatc                  C   s,   d} t  }t|| d}t|t|| d d S )N*   )seed)r   r   r   )r	   camtest r   K/tmp/pip-target-vg8gfxp4/lib/python/skimage/util/tests/test_random_noise.pytest_set_seed	   s   r   c                  C   s~   d} t t }d}t|| d|d}||k}t|| t|  t| |jd |jd   }d}t	|| |ks=J d S )Nr   333333?saltr	   modeamountr      {Gz?)
r   r   r   r   nponessumfloatshapeabs)r	   r
   r   	cam_noisysaltmask
proportion	tolerancer   r   r   	test_salt   s   
 r    c                  C   s6   t jdd} t| ddd}t|g dg dg d S )N      r   r   )r   r   )r   r   r   )r   randomrandr   r   )imagenoisyr   r   r   test_salt_p1    s   r'   c                  C   sD   t jdd} t| dddd}d}tt |dkd |ks J dS )	zBEnsure images where size of a given dimension is 1 work correctly.r   i  r   皙?r   )r   r   r	   g?N)r   r#   r$   r   r   average)r%   r&   r   r   r   r   test_singleton_dim&   s   "r*   c            
      C   s   d} t t }|d d }d}t|| d|d}||k}t|| t|  t| |jd |jd   }d	}t	|| |ksCJ |d
k }t|| ddd}	t|	d
k | |jd |jd   }t	|| |ksoJ d S )Nr          @      ?r   pepperr   r   r   r   )
r   r   r   r   r   zerosr   r   r   r   )
r	   r
   data_signedr   r   
peppermaskr   r   
orig_zeroscam_noisy_signedr   r   r   test_pepper.   s$   
 r4   c                  C   s   d} t t }d}t|| d|dd}t||k|dk}t||k|dk}t|| t|  t|| t|  t	| |  |j
d |j
d	   }d
}t|| |ks^J d| |    k rodk srJ  J d S )Nr   r   zs&p      ?)r	   r   r   salt_vs_pepperr,           r   r   r   g
ףp=
?gffffff?)r   r   r   r   logical_andr   r   r   r/   r   r   r   )r	   r
   r   r   r   r1   r   r   r   r   r   test_salt_and_pepperI   s$   

,r9   c                  C   s   d} t dd }t|| dd}d|   k rdk s J  J t|| dd	d
}d| d   k r7dk s:J  J d|   k rGdk sJJ  J d S )Nr      r;         ?r   )r	   vargMb?g~jt?g333333?gQ?)r	   meanr=   gQ?g{Gz?g;On?)r   r/   r   r=   r>   )r	   datadata_gaussianr   r   r   test_gaussian`   s    $$rA   c                  C   s  d} t dd }t dd }d|d ddd f< d|dd d df< d|dd dd f< t|d	| |d
d}d|d dd df    k rLdk sOJ  J d|d ddd f    k rddk sgJ  J d|dd d df    k r|dk sJ  J d|dd dd f    k rdk sJ  J t |}tt t|d	| |d W d    n1 sw   Y  |d7 }d|d< tt t|d	| |d W d    d S 1 sw   Y  d S )Ni\  r:   r<   gMbP?r(   @   r5   g?localvarF)r   r	   
local_varsclipr7   gMb`?gRQ?gzG?g\(\?gRQ?g{Gz?gQ?)r   r	   rD   r.   )r   r   )r   r/   r   r=   
zeros_liker   raises
ValueError)r	   r?   rD   r@   bad_local_varsr   r   r   test_localvark   s4   0000
"rJ   c                  C   sb   d} t dd }t j| }|ddd}t |||  dd}t|d| ddd	}t|| d S )
Nr   r:   r(   g
cQ?r   r   speckleg{Gz?)r   r	   r>   r=   )r   r/   r#   default_rngnormalrE   r   r   )r	   r?   rngnoiseexpecteddata_speckler   r   r   test_speckle   s   
rR   c                  C   sl   d} t  }t|d| d}t|d| dd}tj| }|t|d d }t|t|dd	 t|| d S )
Nr   poisson)r   r	   Fr   r	   rE      g      p@r7   r,   )	r   r   r   r#   rL   rS   r   r   rE   )r	   r?   r   
cam_noisy2rN   rP   r   r   r   test_poisson   s   rW   c                  C   s   d} t  }t|d d }t|d| dd}t|d| dd}| dkr)| dks+J | dkr7| dks9J t|d| d	d}t|d| d	d}| d
krU| dksWJ | dkrc| dkseJ d S )Nr   r+   r,   rS   TrT   r7         Fgffffff?g?r   r   r   maxmin)r	   r?   r0   cam_poissoncam_poisson2r   r   r   test_clip_poisson       r^   c                  C   s   d} t  }t|d d }t|d| dd}t|d| dd}| dkr)| dks+J | dkr7| dks9J t|d| d	d}t|d| d	d}| d
krU| dk sWJ | dkrc| dk seJ d S )Nr   r+   r,   gaussianTrT   r7   rX   FgQ?gffffffֿ/$?/$rY   )r	   r?   r0   	cam_gauss
cam_gauss2r   r   r   test_clip_gaussian   r_   re   c                  C   s   d} t  }t|d d }t|d| dd}t|d| dd}| dkr)| dks+J | dkr7| dks9J t|d| d	d}t|d| d	d}| d
krU| dksWJ | d
krc| dk seJ d S )Nr   r+   r,   rK   TrT   r7   rX   Fra   rb   rY   )r	   r?   r0   cam_specklecam_speckle_sigr   r   r   test_clip_speckle   r_   rh   c                  C   sD   t d} tt t| d W d    d S 1 sw   Y  d S )N)rB   rB   perlin)r   r/   r   rG   KeyErrorr   )r?   r   r   r   test_bad_mode   s   
"rk   )skimage._sharedr   skimage._shared.testingr   r   numpyr   skimage.datar   skimage.utilr   r   r   r    r'   r*   r4   r9   rA   rJ   rR   rW   r^   re   rh   rk   r   r   r   r   <module>   s&    