o
    e                     @   s  d dl Z d dlZd dlZd dlmZ d dlmZm	Z	m
Z
mZmZ d dlmZmZmZ d dlmZ d dlmZmZmZmZmZmZmZmZmZmZmZmZ d dl m!Z! e"g dg dg d	g d	g d
g dg dg dg dg dg
Z#e#$ Z%de%dddf< e%$ &ej'd Z(e(dej)f g d Z*ej+dej,dZ-de-ddddf< e-$ d Z.ej/dej0dZ1de1ddddddf< de1d< e1$ Z2dd d!Z3dd#d$Z4dd%d&Z5d'd( Z6d)d* Z7d+d, Z8d-d. Z9d/d0 Z:d1d2 Z;d3d4 Z<d5d6 Z=ej>?d7e#dfe#dfe#dfe#d8fe1dfe1dfe1d9fgd:d; Z@d<d= ZAd>d? ZBd@dA ZCdBdC ZDdDdE ZEdFdG ZFdHdI ZGdJdK ZHej>?dLg dMg dNg dOgdPdQ ZIej>?dLdRdSgdTdU ZJej>?dLejKejLejL gdVdW ZMej>?dLdgddggd"fdXdY ZNej>?dLdZd[gddgddggd\d] ZOej>?dLd^gg d_gd`da ZPdbdc ZQddde ZRdfdg ZSdhdi ZTej>?dLg djdkdl ZUej>?dLg dmdndo ZVdpdq ZWdrds ZXdtdu ZYdvdw ZZdxdy Z[dzd{ Z\d|d} Z]d~d Z^dd Z_dd Z`dd Zadd Zbdd Zcdd Zddd Zedd Zfdd Zgdd Zhdd Zidd Zjdd Zkdd Zldd Zmdd Zndd Zodd Zpdd Zqdd Zrdd Zsdd Ztdd Zudd Zvdd Zwdd Zxdd Zydd Zzdd Z{dd Z|dd Z}dd Z~dd ZddÄ Zddń ZddǄ ZddɄ Zdd˄ Zdd̈́ Zddτ Zddф Zddӄ ZddՄ Zddׄ Zddل Zddۄ Zdd݄ Zdd߄ Zdd Zdd Zej>?dedd Zdd Zdd Zdd Zdd Zdd Zdd ZdS )    N)assert_allcloseassert_almost_equalassert_array_almost_equalassert_array_equalassert_equal)datadraw	transform)testing)
COL_DTYPESOBJECT_COLUMNSPROPS#_inertia_eigvals_to_axes_lengths_3D_parse_docs_props_to_dict_require_intensity_imageeuler_number	perimeterperimeter_croftonregionpropsregionprops_table)slic)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   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   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   	      g      $@.r   r      
   dtyper                   @   r&   r&   )r   r   r   r   r   c                    sZ   j \}}tjtjd||d  |ddtjd||d  |dddd\  fddS )	Nr   FZendpointr   ijZindexingc                    s   t |   |   S NnpsumpqXYimg VD:\Projects\ConvertPro\env\Lib\site-packages\skimage/measure/tests/test_regionprops.py<lambda>6       z%get_moment_function.<locals>.<lambda>shaper-   meshgridlinspace)r5   spacingrowscolsr6   r2   r7   get_moment_function2   s
   

rA   r   r   r   c              	      sx   j \}}}tjtjd||d  |ddtjd||d  |ddtjd||d  |dddd\  fdd	S )
Nr   Fr(   r   r   r)   r*   c                    s"   t |  |   |   S r+   r,   )r0   r1   rr3   r4   Zr5   r6   r7   r8   >      " z'get_moment3D_function.<locals>.<lambda>r:   )r5   r>   Zslicesr?   r@   r6   rD   r7   get_moment3D_function9   s   rG   c                    s   j \}}tjtjd||d  |ddtjd||d  |dddd\ t|d}|dd|dd |dd|dd  fdd	S )
Nr   Fr(   r   r)   r*   r>   c                    s"   t  |    |   S r+   r,   r/   r3   r4   cXcYr5   r6   r7   r8   J   rF   z-get_central_moment_function.<locals>.<lambda>)r;   r-   r<   r=   rA   )r5   r>   r?   r@   Mpqr6   rI   r7   get_central_moment_functionA   s   

rM   c               	   C   sp   t ttd } tD ],}z!t| | t| t|  | |kr+tt| |t| t|  W q	 ty5   Y q	w d S Nr   )r   SAMPLEINTENSITY_SAMPLEr   r   getattrlower	TypeErrorregionpropr6   r6   r7   test_all_propsM   s   
rW   c               
   C   st   t ttd } tD ].}z!t| | t| t|  | |kr+tt| |t| t|  W q	 ttfy7   Y q	w d S rN   )	r   	SAMPLE_3DINTENSITY_SAMPLE_3Dr   r   rQ   rR   NotImplementedErrorrS   rT   r6   r6   r7   test_all_props_3d_   s   
r[   c                  C   s<   t td j} | dksJ t tddd j} | dksJ d S )Nr   H   r   r   rH   )r   rO   
num_pixels)r^   r6   r6   r7   test_num_pixelso   s   r_   c                   C   s   t tjdtd t tjdtjd tt t tjdtd W d    n1 s,w   Y  tt t tjdtj	d W d    n1 sKw   Y  tt t tjdt
d W d    d S 1 sjw   Y  d S )Nr   r   r   )r   r-   zerosintZuintpytestraisesrS   floatfloat64boolr6   r6   r6   r7   
test_dtypew   s   "rh   c                   C   s   t tjdtd t tjdtd t tjdtd t tjdtd t tjdtd tt t tjdtd W d    d S 1 sGw   Y  d S )Nr`   r   )r   r   r   )r   r   r   r'   rB   )r   r   r   r   )r   r-   ra   rb   rc   rd   rS   r6   r6   r6   r7   	test_ndim   s   "ri   c                  C   s   d} t td j}t||  dk sJ d}t tddgdd j}t|| dk s,J tjdtjd}d|d	d
d	d
f< t |d j}t|dtd	  dk sUJ t|td dk scJ d S )N   r   r   r   皙?rH      rm   r   r      i  ư>r   rO   feret_diameter_maxr-   absra   uint8sqrt)comparator_resulttest_resultZcomparator_result_spacingZtest_result_spacingr5   rr   r6   r6   r7   test_feret_diameter_max   s     rx   c               
   C   s   d} t td j}t||  dk sJ d}tjdtjd}d|ddddf< t ||d	d j}t|t|d d
 |d |d k d |d d
 |d |d k  d   dk s]J d S )Nrj   r   r   r]   rl   r   r   rn   rH   ro   rp   rq   )rv   rw   r>   r5   rr   r6   r6   r7   test_feret_diameter_max_spacing   s   
ry   c                  C   sP  t jdt jd} d| ddddf< t | fd }t|d j}t |t d d	k s/J d
}t||dd j}t |t |d d d |d d d  |d d d   d	k s^J t |t |d d d |d d d  |d d d   d	ksJ t |t |d d d |d d d  |d d d   d	ksJ d S )Nrl   r   r   r   rn   r   r   i  rp   r   rH      ro   )r-   ra   rt   Zdstackr   rr   rs   ru   )r5   Zimg_3drr   r>   r6   r6   r7   test_feret_diameter_max_3d   sD   


r{   zsample,spacingr   r   )r   r   r   c                 C   s>   t | |dd j}t| |rt|nd }||ksJ d S )NrH   r   r   )r   arear-   r.   prod)sampler>   r}   Zdesiredr6   r6   r7   	test_area   s   r   c                  C   s  t td j} t| ddtjd tjd f t tddd j} t| ddtjd tjd f t }d|d d df< t |d j} t| ddtjd tjd d f t |ddd j} t| ddtjd tjd d f t td j} t| d t tddd j} t| d d S )	Nr   r   r|   rH   r   r   )r   r   r      r   r   )      ?r   r"   )r   rO   bboxr   r;   copyrX   )r   
SAMPLE_modr6   r6   r7   	test_bbox   s   ""
r   c                  C   s.   t jtddd} t| d j}t|tj d S )Nr!   constantmoder   )r-   padrO   r   	area_bboxr   size)padded	bbox_arear6   r6   r7   test_area_bbox   s   r   c                  C   s@   d} t jtddd}t|| dd j}t|tjt |   d S )N)r   r   r!   r   r   rH   r   )r-   r   rO   r   r   r   r   r~   )r>   r   r   r6   r6   r7   test_area_bbox_spacing   s   r   c                  C   s  t td j} t| d d t| d d t| d d t| d d	 t| d
 d t| d d t| d d ttdd}t|dd| d  t|dd| d  t|dd| d  t|dd| d  t|dd| d
  t|dd| d  t|dd| d  d S )Nr   r   r   g    @{@r   r   g
r'   g@UUUUUr]   g@8_r   r   gpq@r|   gKh/A@r   r   ga2rH   r   r   r   )r   rO   moments_centralr   rM   )mu
centralMpqr6   r6   r7   test_moments_central   s    r   c                  C   s   d} t t| d}tt| dd j}t|d |dd t|d |dd t|d |d	d	 t|d
 |dd	 t|d |dd t|d |d	d t|d |dd d S )Ng?g?rH   r   r   r   r   r   r'   r   r]   r   r|   r   )rM   rO   r   r   r   )r>   r   r   r6   r6   r7   test_moments_central_spacing
  s   r   c                  C   s^   t td j} t| d ttdd}|dd|dd }|dd|dd }t||f|  d S )Nr   )g@g98"@r'   rH   r   )r   rO   centroidr   rA   )r   rL   rK   rJ   r6   r6   r7   test_centroid  s   
r   c                  C   s\   d} t t| d}|dd|dd }|dd|dd }tt| dd j}t|||f d S )Nr   rH   r   r   )rA   rO   r   r   r   r>   rL   rK   rJ   r   r6   r6   r7   test_centroid_spacing&  s   r   c                  C   s   t td j} t| d ttdd}|ddd|ddd }|ddd|ddd }|ddd|ddd }t|||f|  d S )Nr   )gI?:?r   rB   rH   r   )r   rX   r   r   rG   )r   MpqrcZrK   rJ   r6   r6   r7   test_centroid_3d1  s   
r   r>   ) @皙@ffffff@)r$   r$   r$   )r   r   r   c                 C   sz   t t| d}|ddd|ddd }|ddd|ddd }|ddd|ddd }tt| dd j}t||||f dS )!Test the _normalize_spacing code.rH   r   r   N)rG   rX   r   r   r   )r>   r   r   rK   rJ   r   r6   r6   r7   test_spacing_parameter_3d?  s   r   )r   y              ?y      ?        c                 C   sF   t jtdd tt| dd j W d   dS 1 sw   Y  dS )r   z3Element of spacing isn't float or integer type, gotmatchrH   r   N)rc   rd   rS   r   rO   r   rH   r6   r6   r7   $test_spacing_parameter_complex_inputO  s   "r   c                 C   sB   t t tt| dd j W d   dS 1 sw   Y  dS )r   rH   r   Nrc   rd   
ValueErrorr   rO   r   rH   r6   r6   r7   test_spacing_parameter_nan_infX  s   "r   c                 C   sF   t jtdd tt| dd j W d    d S 1 sw   Y  d S )Nz%spacing isn't a scalar nor a sequencer   rH   r   r   rH   r6   r6   r7   test_spacing_mismtaching_shapea  s   "r   r   r   c                 C   sZ   t t| d}|dd|dd }|dd|dd }ttt| dd j}t|||f dS )r   rH   r   r   intensity_imager>   N)rA   rP   r   rO   centroid_weightedr   r   r6   r6   r7   test_spacing_parameter_2di  s
   r   	bad input)r   r   r   c                 C   sD   t t ttt| dd j W d   dS 1 sw   Y  dS )r   r   r   N)rc   rd   r   r   rO   rP   r   rH   r6   r6   r7   #test_spacing_parameter_2d_bad_inputw  s   "r   c                  C   s   t td j} | dksJ d S )Nr   }   )r   rO   area_convexr}   r6   r6   r7   test_area_convex  s   r   c                  C   s0   d} t t| dd j}|dt|  ksJ d S )Nr   r   rH   r   r   )r   rO   r   r-   r~   )r>   r}   r6   r6   r7   test_area_convex_spacing  s   r   c                  C   sb   t td j} tg dg dg dg dg dg dg dg d	g d	g d
g
}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   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   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   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   rO   Zimage_convexr-   arrayr   )r5   refr6   r6   r7   test_image_convex  s   r   c                  C   s   t jdt jd} t ddgddgddgg}d| |d d df |d d df f< t| d j}t|| t| dd	d j}t|| d S )
Nr`   r   r   r   r   r   r   )r   333333?rH   )r-   ra   int8r   r   coordsr   )r   r   prop_coordsr6   r6   r7   test_coordinates  s   $
r   )Nr   r   r'   )r   r   c                 C   s   t jdt jd}t ddgddgddgg}d||d d df |d d df f< t|| dd j}| d u r8|}n|t |  }t|| d S )	Nr`   r   r   r   r   r   r   rH   r-   ra   r   r   r   Zcoords_scaledr   r>   r   r   r   Zdesired_coordsr6   r6   r7   test_coordinates_scaled  s   $r   )Nr   r   )皙?r   r   c                 C   s   t jdt jd}t g dg dg dg}d||d d df |d d df |d d df f< t|| d	d j}| d u r?|}n|t |  }t|| d S )
Nr%   r   rB   )r   r   r   )r   r   r   r   r   r   rH   r   r   r6   r6   r7   test_coordinates_scaled_3d  s   2r   c                  C   sR   t jtddd} tj\}}t| d j}tdd| tdd| f}t|| d S )N)r   r   )r!   r   r   r   r   r   r!   r-   r   rO   r;   r   slicer   )r   nrowncolresultexpectedr6   r6   r7   
test_slice  s
   
r   c                  C   sh   t jtddd} tj\}}t| d j}tdd| tdd| f}d}t| |dd j}t|| d S )	Nr   r   r   r   r   r!   )r   r   rH   r   )r   r   r   r   r   r>   r6   r6   r7   test_slice_spacing  s   
r   c                  C   s   t td j} t| d t tddd j} t| d tjdtd}d|d< t |d j} t| d t |d	dd j} t| d d S )
Nr   gPq~q?      ?r   rH   )r!   r!   r   r   r   r   r   r   )r   rO   Zeccentricityr   r-   ra   rb   )Zepsr5   r6   r6   r7   test_eccentricity  s   


r   c                  C   s\   t td j} t| d d}t t|dd j} tj| d d  }t|t t|  d S )Nr   g3&#@r   r   rH   r$   r   )r   rO   Zequivalent_diameter_arear   r-   pir.   r~   )Zdiameterr>   Zequivalent_arear6   r6   r7   test_equivalent_diameter_area  s   
r   c                  C   s  dD ]>} t t| dd j}|dksJ t }d|d< t || dd j}|dks*J ttd}|dks5J t|d}|dks@J qttd}|dksLJ ttd}|dksWJ td	}d|d
dd
dd
df< t|d}|dkstJ d|ddddddf< t|d}|dksJ d S )N)r'   )r   g?rH   r   r"   r   r   r   r   )d   r   r   (   <   -   7   )r   rO   r   r   rX   r-   ra   )r>   enr   ZSAMPLE_3D_2r6   r6   r7   test_euler_number  s,   






r   c                  C   s8   t td j} t| d t tddd j} t| d d S )Nr   皙?)r!   r   rH   )r   rO   extentr   )r   r6   r6   r7   test_extent     
r   c                  C   sh   t td j} tg d}t| | tt t tddd j W d    d S 1 s-w   Y  d S )Nr   )g?z~?gu/,?g'?gHU-T?g^L>gw-g?b1۾r]   rH   )	r   rO   Z
moments_hur-   r   r   r
   rd   rZ   )hur   r6   r6   r7   test_moments_hu  s   

"r   c                  C   sJ   t td j} t| t t td j} t| tddddddf  d S )Nr   r   r   r   )r   rO   imager   rX   r5   r6   r6   r7   
test_image+  s   
$r   c                  C   s4   t td j} t| d t td j} t| d d S )Nr   r   )r   rO   labelr   rX   r   r6   r6   r7   
test_label3  s   
r   c                  C   s$   t td j} | ttksJ d S rN   )r   rO   area_filledr-   r.   r   r6   r6   r7   test_area_filled;  s   r   c                  C   s4   t  } d| d< t| d j}|tt ksJ d S )Nr   r   )rO   r   r   r   r-   r.   )r   r}   r6   r6   r7   test_area_filled_zero@  s   r   c                  C   st   t  } d| d< d}tt |dd j}|tt t| ks!J t| |dd j}|tt t| ks8J d S )Nr   r   )r   r   rH   )rO   r   r   r   r-   r.   r~   )r   r>   r}   r6   r6   r7   test_area_filled_spacingG  s    r   c                  C   s8   t td j} t| t t tddd j} t| t d S )Nr   r   rH   )r   rO   image_filledr   r   r6   r6   r7   test_image_filledS  r   r   c                  C   s   t td j} d}t| | t tddd j} t| d|  ddlm} tjdtjd}|d	d	d
dt	dd\}}d|||f< t |ddd j}t |d d d ddd j
}t|| dksaJ t |d d d d df ddd j} t| |dd d S )Nr   g8D0@r   rH   r   ellipse)rm      r   r   r"   r   r   rotationr   r'   rk   r|   decimal)r   rO   axis_major_lengthr   skimage.drawr   r-   ra   rt   deg2radaxis_minor_lengthrs   lengthZtarget_lengthr   r5   rrccZlength_wo_spacingr6   r6   r7   test_axis_major_lengthZ  s$   
$r  c                  C       t ttdd j} t| d d S )Nr   r   r   )r   rO   rP   Zintensity_maxr   Z	intensityr6   r6   r7   test_intensity_maxq     
r	  c                  C   r  )Nr  r   gqq?)r   rO   rP   Zintensity_meanr   r  r6   r6   r7   test_intensity_meanw  r
  r  c                  C   r  )Nr  r   r   )r   rO   rP   Zintensity_minr   r  r6   r6   r7   test_intensity_min}  r
  r  c                  C   s   t td j} d}t| | t tddd j} t| d|  ddlm} tjdtjd}|d	d
dd	t	dd\}}d|||f< t |ddd j}t |d d d ddd j}t
|| dksaJ t |d d d ddd j} t| |dd d S )Nr   g7z#@r   rH   r   r   )r      r   r!   r&   r      r   r   r'   r   rk   r]   r   )r   rO   r   r   r   r   r-   ra   rt   r   rs   r  r6   r6   r7   test_axis_minor_length  s$   
r  c                  C   sr  t td j} t| d d t| d d t| d d t| d d	 t| d
 d t| d d t| d d t| d d t| d d t| d d ttdd}t|dd| d  t|dd| d  t|dd| d  t|dd| d  t|dd| d
  t|dd| d  t|dd| d  t|dd| d  t|dd| d  t|dd| d  d S )Nr   r   r   g      R@r   r   g     @@r   g     @r   g    @V@r   r   g     y@r'   g     l@r|   g    @m@r   g     x@r]   g     A@r   g     P@rH   r   r   r   )r   rO   momentsr   rA   )mrL   r6   r6   r7   test_moments  s,   r  c                  C   s   d} t t| dd j}tt| d}t|d |dd t|d |dd t|d |dd t|d	 |dd
 t|d |dd t|d |dd t|d |dd t|d |dd t|d |dd t|d |d
d d S )N)r   333333?rH   r   r  r  r   r   r   r   r   r  r'   r|   r   r]   r   )r   rO   r  rA   r   )r>   r  rL   r6   r6   r7   test_moments_spacing  s   r  c                  C   sf   t td j} t| d d t| d d t| d d t| d d	 t| d
 d t| d d d S )Nr   r   e/
?r   cZr'   ,lC@r|   ^ ^H?r   /O懵?r]   ⪹Rgr   rO   Zmoments_normalizedr   )nur6   r6   r7   test_moments_normalized  s   r   c                  C   sn   d} t t| dd j}t|d d t|d d t|d d	 t|d
 d t|d d t|d d d S )Nr   rH   r   r   r  r   r  r'   r  r|   r  r   r  r]   r  r  )r>   r  r6   r6   r7   test_moments_normalized_spacing  s   r!  c                  C   s  t td j} d}t| | t tddd j} t| | tjdtd}t |d j}t|tjd  t |ddd j}t|t	d	t
d
  t t|d j}t|tj d  t t|ddd j}t|t	d	t
d
   t t|d j}t|tj d  t t|ddd j}t|t	d	t
d
   t tt|d j}t|tjd  t tt|ddd j}t|t	d	t
d
  d S )Nr   g/9vr   rH   r   r   r   r|   r   g      ?)r   rO   orientationr   r-   eyerb   mathr   Zarccosru   ZflipudZfliplr)ZorientZtarget_orientZdiagZorient_diagr6   r6   r7   test_orientation  s,   

r%  c            	      C   s6  t g dg dg dg dg} t g dg dg dg dg}t g d	g d
g dg dg}t | ||f}t|}dd |D }t jj||d ddd t|d d t|d d t|d d d}ttt|dd j	}t
t|d}t|d |dd t|d |dd t|d |dd t|d |dd t|d |dd t|d  |dd t|d! |dd t|d" |dd t|d# |dd t|d$ |dd t|d% |dd t|d& |dd t|d' |dd t|d( |dd t|d) |dd t|d* |dd d S )+N)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   c                 S   s   g | ]}|j qS r6   )r"  ).0rV   r6   r6   r7   
<listcomp>  s    z/test_orientation_continuity.<locals>.<listcomp>r   r   g{Gz?ZrtolZatolge{g-DT!r   g4,1g뿩g	@r   r   rH   r  r  r   r   r   r  r'   r|   r   r   r]   r   r   r   r   r   r   r   r   )r-   r   Zhstackr   r
   r   r   rO   rP   moments_weighted_centralrM   )	Zarr1Zarr2Zarr3r   propsZorientationsr>   wmur   r6   r6   r7   test_orientation_continuity  sF   """r/  c                  C   s   t td j} d}t| | t tddd j} t| d|  ttddd} t| d	 tt t td
dd j} W d    d S 1 sFw   Y  d S )Nr   g֟K@r   rH   r   doubler#   )ZneighborhoodgXg	jG@r]   )r   rO   r   r   astyper
   rd   rZ   )perZ
target_perr6   r6   r7   test_perimeter"     

"r3  c                  C   s   t td j} d}t| | t tddd j} t| d|  ttddd} t| d tt t td	dd j} W d    d S 1 sFw   Y  d S )
Nr   gχ?N@r   rH   r   r0  )Z
directionsgP@r]   )r   rO   r   r   r1  r
   rd   rZ   )r2  Ztarget_per_crofr6   r6   r7   test_perimeter_crofton0  r4  r5  c                  C   s<   t td j} d}t| | t tddd j} t| | d S )Nr   g;On?)r   r   rH   )r   rO   solidityr   )r6  Ztarget_solidityr6   r6   r7   test_solidity>  s
   
r7  c                  C   s   t } tt| d dd }tt| dddf dd }tt| dd }tj|jt|j tj|jt|jdddf  t	|jtj
fksLJ t	|jtj
dfksYJ t	|jtj
| j	d fksiJ tt| dd	}t|t|jks|J dS )
zBTest for https://github.com/scikit-image/scikit-image/issues/6860.).r   r  r   .r   Nr   )r   )r   
properties)#INTENSITY_FLOAT_SAMPLE_MULTICHANNELr   rO   r-   r
   r   r   Zsqueezer   r;   ndimr   lenr   )r   Zrp0Zrp1Zrpmtabler6   r6   r7   )test_multichannel_centroid_weighted_tableG  s&   



r=  c                  C   s  t ttdd j} tg dg dg dg dg}tjdd t| | ttd	d
}t	|dd|d  t	|dd|d  t	|dd|d  t	|dd|d  t	|dd|d  t	|dd|d	  t	|dd|d  t	|dd|d  t	|dd|d  t	|dd|d  t	|dd|d  t	|dd|d  t	|dd|d  t	|dd|d  t	|dd|d  t	|dd|d  d S )Nr  r   )     R@g($=gp8"@g)g0   Ngc"UgiC*Nڠ@g]}ᘎ)g|}@g(kbgz*@gqo)g}gF%uݓg~2	@g"0V0r   )	precisionr'   rH   r  r   r  r   r   r   r   r  r|   r   r   r]   r   r*  r   r+  r   r   )
r   rO   rP   r,  r-   r   Zset_printoptionsr   rM   r   )r.  r   r   r6   r6   r7   test_moments_weighted_central\  s>   


r@  c               	   C   s  t g dg dg dg dg dg dg} d}tttdd j}d}t|| ttdd}|dd	|dd }|d	d|dd }t||f| d
}tt|d}|dd	|dd }|d	d|dd }ttt|dd j}t|||f t|dt |  t| | |dd j}t|dt |  d}tt|d}|dd	|dd }|d	d|dd }ttt|dd j}t|||f t| | |dd j}t||t |  d S )N)r   r   r   r   r   r   )r   r   r   r   r   r   )      @rA  r  r   )g8 u)@g1<S"@r'   rH   r   r   r   r   )g?gffffff?)	r-   r   r   rO   rP   r   r   rA   r   )Zsample_for_spacingZtarget_centroid_wspacingr   Ztarget_centroidrL   rK   rJ   r>   r6   r6   r7   test_centroid_weighted  sL   

rB  c                  C   sl   t ttdd j} tg d}t| | tt	 t tddd j W d    d S 1 s/w   Y  d S )Nr  r   )ga'R?gg	?y?g"
-?g!mT?gODګ>gpcgʝ~ܾr]   rH   )
r   rO   rP   Zmoments_weighted_hur-   r   r   r
   rd   rZ   )Zwhur   r6   r6   r7   test_moments_weighted_hu  s   

	"rC  c                  C   s  t ttdd j} tg dg dg dg dg}t| | ttdd}t|dd|d	  t|dd
|d  t|dd|d  t|dd|d  t|d
d|d  t|d
d
|d  t|d
d|d  t|d
d|d  t|dd|d  t|dd
|d  t|dd|d  t|dd|d  t|dd|d  t|dd
|d  t|dd|d  t|dd|d  d S )Nr  r   r>  g     ؅@g     @g    P@g     y@g     @g    @g    .y!Ag     |@g    E@g    Ag   MAg    P@g    \Ag    ?Ag   p1zAr'   rH   r  r   r  r   r   r   r   r  r|   r   r   r]   r   r*  r   r+  r   r   	r   rO   rP   Zmoments_weightedr-   r   r   rA   r   )wmr   rL   r6   r6   r7   test_moments_weighted  s<   

rJ  c                  C   s  t ttdd j} tg dg dg dg dg}t| | d}t tt|dd j}tt|d	}t|d
 |dd t|d |dd t|d |dd t|d |dd t|d |dd t|d |dd t|d |dd t|d |dd t|d |dd t|d |dd t|d |dd t|d |dd t|d |dd t|d |dd t|d |dd t|d |dddd d S )Nr  r   rD  rE  rF  rG  r)  r   rH   r  r  r   r   r   r   r   r  r'   r|   r   r   r]   r   r*  r   r+  r   r   r&   r   rH  )rI  r   r>   r.  rL   r6   r6   r7   test_moments_weighted_spacing  sL   

rK  c                  C   sP   t ttdd j} ttjtjddgtjdddgg dg d	g}t| | d S )
Nr  r   !%su?Smvlc6r?,>AkS*]?NCyicd?\}$bg7"h ;i\q?A׸R)r   rO   rP   moments_weighted_normalizedr-   r   nanr   )wnur   r6   r6   r7    test_moments_weighted_normalized  s   
r^  c                  C   s  d} t tt| dd j}ttjtjddgtjdddgg d	g d
g t|d d t|d d t|d d t|d d t|d d t|d d t|d d t|d d t|d d t|d d t|d d t|d d t|d d d S )Nr   r   r   rL  rM  rN  rO  rP  rQ  rV  r   r   r'   r|   r   r   rR  r]   rS  r   rT  r*  rU  r   rW  r+  rX  r   rY  rZ  )r   rO   rP   r[  r-   r   r\  r   )r>   r]  r6   r6   r7   (test_moments_weighted_normalized_spacing  s,   r_  c                  C   s@   t td } tddg}t t|dd }t| j|j|  d S )Nr   i   i   )offset)r   rO   r-   r   r   r   )r-  r`  Zprops_offsetr6   r6   r7   test_offset_features#  s   ra  c                  C   sP   t jdtd} d| d d d d f< t| }t|dksJ |d jdks&J d S )Nr   r   r   r   r   )r-   emptyrb   r   r;  r   apsr6   r6   r7   test_label_sequence+  s
   rf  c                  C   s*   t jdtd} t| }t|dksJ d S )Nr   r   r   )r-   ra   rb   r   r;  rc  r6   r6   r7   test_pure_background3  s   rg  c                     sJ   t t  fdd} tt |   W d    d S 1 sw   Y  d S )Nc                      s    d j  d S rN   )image_intensityr6   re  r6   r7   get_intensity_image<     z)test_invalid.<locals>.get_intensity_image)r   rO   rc   rd   AttributeError)rj  r6   ri  r7   test_invalid9  s
   "rm  c                  C   sL   t dgdgg} tt tt|  W d    d S 1 sw   Y  d S Nr   )r-   r   rc   rd   r   r   rO   )Zwrong_intensity_sampler6   r6   r7   test_invalid_sizeC  s   "ro  c                  C   s   t jdtd} d| ddddf< d| ddddf< t| }|d }t| }|d }|d }t||kd	d
 t||kd	d d S )N)r   r   r   r   r      r   2   c   TzSame regionprops are not equalzDifferent regionprops are equal)r-   ra   rb   r   r   )Zarrregionsr1r2Zr3r6   r6   r7   test_equalsI  s   rv  c                     sX   t td   fdd D } t ttdd   fdd D }t| t|k s*J d S )Nr   c                       i | ]}| | qS r6   r6   r&  r0   rU   r6   r7   
<dictcomp>[      z*test_iterate_all_props.<locals>.<dictcomp>r  c                    rw  r6   r6   rx  ry  r6   r7   rz  ^  r{  )r   rO   rP   r;  )p0p1r6   ry  r7   test_iterate_all_propsY  s
   r~  c                  C   sZ   t  } t| d }|j}d|jd d< |j}t|| d|_|j}t||ks+J d S )Nr   r   r   F)	rO   r   r   r   Z_label_imager   Z_cache_activer-   any)r   rU   Zf0f1r6   r6   r7   
test_cachec  s   
r  c                  C   s   dd } t | j}ttd }t }dd t|D }t|}t|}|rAt|| |d }d|vs4J t|dd	ks?J d S t|d d S )
Nc                   S      dS )fooNr6   r6   r6   r6   r7   r  u  s    z&test_docstrings_and_props.<locals>.foor   c                 S   s   g | ]	}| d s|qS )_)
startswith)r&  r  r6   r6   r7   r'  }  r9   z-test_docstrings_and_props.<locals>.<listcomp>r[  	iteration
r   )	rg   __doc__r   rO   r   dirr;  r   split)r  Zhas_docstringsrU   docsr-  Znr_docs_parsedZnr_propsZdsr6   r6   r7   test_docstrings_and_propst  s   

r  c               	   C   s   t t} t| }|tdgtdgtdgtdgtdgdks(J t t} t| ddd}|tdgtd	gtdgtdgtdgtdgd
ksXJ d S Nr   r   r   rj   )r   zbbox-0zbbox-1zbbox-2zbbox-3r   r}   r   +r8  	separatorr\   )r   r}   bbox+0bbox+1bbox+2bbox+3)r   rO   r   r-   r   )rs  outr6   r6   r7   test_props_to_dict  s   r  c               	   C   s   t t} | tdgtdgtdgtdgtdgdks$J t tddd} | tdgtd	gtdgtdgtdgtdgd
ksPJ d S r  )r   rO   r-   r   r  r6   r6   r7   test_regionprops_table  s   r  c                  C   s,   t tdd} |  D ]	}|dsJ q
d S )N)local_centroidr8  r  )r   rO   keysr  )r  keyr6   r6   r7   1test_regionprops_table_deprecated_vector_property  s   r  c                  C   s&   t tdd} t|  dgksJ d S )N)r   r  r   )r   rO   listr  r  r6   r6   r7   1test_regionprops_table_deprecated_scalar_property  s   r  c               	   C   s   t tt} tttt d}t D ]a\}}t| D ]X\}}|| }t	|s0|t
v s0|tju r:t||| |  qt|tjrC|jnt|f}t|D ]$}dtt|f| }	t|dkrb|n|d }
t||
 ||	 |  qMqqd S )Nr  -r   r   )r   rO   INTENSITY_FLOAT_SAMPLEr   r   r  items	enumerater-   Zisscalarr   Zobject_r   
isinstanceZndarrayr;   r;  Zndindexjoinmapstrr   )rs  Z	out_tablerV   r    iregrpr;   indZmodified_proplocr6   r6   r7   (test_regionprops_table_equal_to_original  s&   


r  c                  C   s   t tjdtdddd} t| dksJ t| d dksJ t| d	 dks(J t| d
 dks2J t| d dks<J t| d dksFJ t| d dksPJ d S )Nr   r   r  r  r  r&   r   r   r}   r  r  r  r  )r   r-   ra   rb   r;  r  r6   r6   r7   !test_regionprops_table_no_regions  s   r  c                   C   s&   t t tt t ksJ d S r+   )setr   r  unionr   r   valuesr6   r6   r6   r7   test_column_dtypes_complete  s   &r  c               
   C   s   d} t ttdd }tD ]j}|| }|tv rt| tksJ qtt|d }t	|tj
rGt| tksFJ | d| d|  dt|  qt	|tjrft| tkseJ | d| d|  dt|  qJ | d| d|  dt|  d S )Nzmismatch with expected type,r  r   z dtype  )r   rO   rP   r   r   objecttyper-   ZravelZ
issubdtypeZfloatingre   integerrb   )msgrU   colrC   tr6   r6   r7   test_column_dtypes_correct  s(   r  c                 C   s
   t | S )z"a short test for an extra propertyr,   )
regionmaskr6   r6   r7   
pixelcount  s   
r  c                 C   s   t ||  S r+   )r-   median)r  rh  r6   r6   r7   intensity_median  rk  r  c                 C   r  rn  r6   )r  rh  Zsuperfluousr6   r6   r7   too_many_args     r  c                   C   r  rn  r6   r6   r6   r6   r7   too_few_args  r  r  c                  C   s.   t ttfdd } | jttdkksJ d S )Nextra_propertiesr   r   )r   rO   r  r-   r.   ry  r6   r6   r7   test_extra_properties  s   r  c                  C   s4   t tttfdd } | jtttdk ksJ d S N)r   r  r   r   )r   rO   rP   r  r-   r  ry  r6   r6   r7   test_extra_properties_intensity  s    r  intensity_propc                 C   sb   t td }tt}t||  W d    n1 sw   Y  d|  d}|t|jks/J d S )Nr   zAttribute 'z<' unavailable when `intensity_image` has not been specified.)r   rO   rc   rd   rl  rQ   r  value)r  rU   eZexpected_errorr6   r6   r7   test_intensity_image_required  s   
r  c                  C   sH   t t tttfdd } | j}W d    d S 1 sw   Y  d S Nr  r   )rc   rd   rl  r   rO   r  rU   r  r6   r6   r7   +test_extra_properties_no_intensity_provided  s   "r  c                  C   s   t t tttfdd } | j}W d    n1 sw   Y  t t tttfdd } | j}W d    d S 1 s>w   Y  d S r  )rc   rd   rl  r   rO   r  r  r  r6   r6   r7   test_extra_properties_nr_args   s   "r  c                  C   sN   t ttttfdd } | jtttdk ksJ | jttdkks%J d S r  )r   rO   rP   r  r  r-   r  r.   ry  r6   r6   r7   test_extra_properties_mixed)  s   r  c                  C   sH   t ttdttfd} t| d tddg t| d tddg d S )	Nr   )r   r8  r  r  r$   rA  r  r   r   )	r   SAMPLE_MULTIPLEINTENSITY_SAMPLE_MULTIPLEr  r  r   r-   r   r   r  r6   r6   r7   test_extra_properties_table2  s   r  c            	      C   s   t  ddddddf } | d }t| tdd}t|d }t||tgd| }t|| tgd| }t	t
 dg D ]'}|| }|| }t|t|krYt|| q?t|t|d d	d	d
 q?dS )z2Test that computing multichannel properties works.Nr   ).r   r   )Zstart_labelr   r  r  g-q=r(  )r   Z	astronautr   r1  re   r-   maxr   r  r  r   r  r;   r   r   Zasarray)	ZastroZastro_greenlabelsZsegment_idxrU   Zregion_multirV   r0   Zp_multir6   r6   r7   test_multichannel<  s2   r  c            	      C   s   d} t j|  t}tj|g ddd}tjg ddd}t	||j
}t|d }|j}t|}td	d
 | D dd}t||D ]\}}t|| d| k sQJ qAt|j|d  dk s_J t|j|d  dk slJ dS )zuVerify that estimated axis lengths are correct.

    Uses an ellipsoid at an arbitrary position and orientation.
    )rm   r   rq  ))r  rj   )r  r  )r   rm   r   )Z	pad_widthr   )r   r  r   r   )r   Zdimensionalityr   c                 S   s   g | ]}d | qS )r   r6   )r&  hr6   r6   r7   r'  q  s    z2test_3d_ellipsoid_axis_lengths.<locals>.<listcomp>T)reverseg{Gz?gHz>r   N)r   Z	ellipsoidr1  rb   r-   r   r	   ZEuclideanTransformndiZaffine_transformparamsr   Zinertia_tensor_eigvalsr   sortedziprs   r   r   )	Zhalf_lengthsr  Rr  ZevsZaxis_lengthsZexpected_lengthsZax_len_expectedZax_lenr6   r6   r7   test_3d_ellipsoid_axis_lengthsY  s   
r  )r'   )rB   )r$  numpyr-   rc   Zscipy.ndimageZndimager  Znumpy.testingr   r   r   r   r   Zskimager   r   r	   Zskimage._sharedr
   Zskimage.measure._regionpropsr   r   r   r   r   r   r   r   r   r   r   r   Zskimage.segmentationr   r   rO   r   rP   r1  rf   r  Znewaxisr9  r#  Zint32r  r  ra   rt   rX   rY   rA   rG   rM   rW   r[   r_   rh   ri   rx   ry   r{   markZparametrizer   r   r   r   r   r   r   r   r   r   r   r\  infr   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/  r3  r5  r7  r=  r@  rB  rC  rJ  rK  r^  r_  ra  rf  rg  rm  ro  rv  r~  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r6   r6   r6   r7   <module>   sN   8











	


!&	$)"


		
