o j @s ddlZddlZddlZddlmZmZmZmZm Z ddl m Z ddl m Z ddl mZddl mZddlmZddlmZdd lmZdd lmZdd lmZejd ejejgd dZddZddZ ddZ!ddZ"ddZ#ddZ$ddZ%ddZ&ejd ej'ejej(gdd Z)ejd ej'ejej(gd!d"Z*d#d$Z+ejd%d d&gejd ej,ej-ej(gejd'gd(d)d*Z.ej/0de 1Z2e3e2Z4e5e4d+d,Z4d-d.Z6d/d0Z7ejd ej'ejej(gd1d2Z8d3d4Z9d5d6Z:ejd7d&dd8gfd dd9gfd:d;dZ;ejd7d&d?d@gfd dAd8gfdBdCdDgfgdEdFZdJZ?dKZ@dLZAdMdNZBdOdPZCejd ej'ejej(gdQdRZDdSdTZEejd ejejFej'ejej(gdUdVZGdWdXZHdYdZZId[d\ZJd]d^ZKd_d`ZLdadbZMdcddZNejjOe ejPe dekdfdgejdhdidjgdkdlZQejdmdnej,fdoejRfdpejRfdqejRfdreSfgdsdtZTdudvZUdwdxZVdydzZWejd ej'ejej(gd{d|ZXd}d~ZYddZZddZ[ddZ\ddZ]ddZ^ddZ_ddZ`ddZaddZbddZcejd ej'ejej(gddZdddZeddZfddZgddZhejd ej'ejej(gddZiddZjddZkejd ej'ejej(gddZlddZmddZnddZoddZpddZqddZrejdejsejtejugddZvddZwddZxdS)N)assert_allcloseassert_almost_equalassert_array_almost_equalassert_array_equal assert_equal)Version)data)exposure)util)rgb2gray)intensity_range) dtype_range)expected_warnings)_supported_float_typedtypecCsXtjddg|d}tjtddtj|dd\}}WddS1s%wYdS)Ndrz+Incorrect value for `source_range` argumentmatchZfoobar source_range)nparraypytestraises ValueErrorr histogramrim frequencies bin_centersr"e/var/www/html/Deteccion_Ine/venv/lib/python3.10/site-packages/skimage/exposure/tests/test_exposure.pytest_wrong_source_ranges"r$cCsjtjddgtjd}t|\}}t|tdd|ddks"J|ddks*Jt|ddddS)Nrrrerrrint8r rrarangerr r!r"r"r#test_negative_overflow# r+cCsjtjddgtjd}t|\}}t|tdd|ddks"J|ddks*Jt|ddddS)Nirrrr&r'r*r"r"r#test_all_negative_image,r,r-cCsTtjddgtjd}t|\}}tt|t|t|ddt|dddS)N rrrr)rrr(r rrlenr*r"r"r#test_int_range_image5s r0cCsxtjddgddggtjd}tj|dd\}}t|jdD] }tt||t|qt|ddt|dddS) Nr.rfrr) channel_axisr) rrr(r rrangeshaperr/)rr r!chr"r"r#!test_multichannel_int_range_image=s r7cCsrtjddgtjd}tj|dd\}}t|tddt|ddt|ddt|d dt|jd dS) Nr.rrrrrr&r%r8) rruint8r rrr)rr5r*r"r"r#test_peak_uint_range_dtypeFr;cCsrtjddgtjd}tj|dd\}}t|tddt|dd t|d d t|d d t|jd dS)Nr.rrrrr&rr9) rrr(r rrr)rr5r*r"r"r#test_peak_int_range_dtypePr<rBcCsHtjdddtjd}tj|dd\}}t|tddt|jddS)Nrr8rrrr9) rlinspacer:r rrr)rr5r*r"r"r#test_flat_uint_range_dtypeZrEcCsHtjdddtjd}tj|dd\}}t|tddt|jddS)Nr=r>r8rrrr9) rrDr(r rrr)rr5r*r"r"r#test_flat_int_range_dtypearFrGcCsLtjddg|d}tj|dd\}}|j|ksJt|tddddS)Nr.rrZnbins?)rrr rrrr)rr"r"r#"test_peak_float_out_of_range_imagehsrLcCsntjddg|d}d}tj||dd\}}|j|ksJtt|ddtt|ddtt |ddS) Nr.rrr)rJrg?) rrr rrrminmaxrr/)rrrJr r!r"r"r#"test_peak_float_out_of_range_dtypeqsrQcCsvtjgdtjd}tj|ddd\}}td}d|d<d |d <t||tj|dd d\}}|d }t||dS) N)rrCrCrrF)r normalizer8r&rrTg@)rrr:r rzerosr)rr r!expectedr"r"r#test_normalize|s   rVrimager3)rr&rc Csd}|d|d}t|\}}tjtj|||dtj|||df|d}tj|||d\}} t|tjrAt| t ||d|dd|ksKJ|dddksUJ|dddks_J|dd|ksiJdS) z-Check that all channels use the same binning.)r1r1rr&r)Zaxis)rr3rN) r rstackfullr rZ issubdtypeintegerrr)) rrr3r5Z channel_sizeZiminZimaxrr r!r"r"r#(test_multichannel_hist_common_bins_uint8s$  r[g@rcCs&tt}tjtdd}t||dS)z)Check integer bins used for uint8 images.rMrIN)r equalize_hist test_img_intr)Zimg_eq0Zimg_eq1r"r"r#test_equalize_uint8_approxs r^cCs.tt}t|}t|\}}t|dS)N)r Z img_as_ubytetest_imgr r\cumulative_distributioncheck_cdf_slope)imgimg_eqcdf bin_edgesr"r"r#test_equalize_ubytes   rfcCs\ttj|dd}t|}|jt|ksJt|\}}t ||jt|ks,JdS)NFcopy) r img_as_floatr_astyper r\rrr`ra)rrbrcrdrer"r"r#test_equalize_floats  rkcCsltt}ttj}d|ddddf<tj||d}t|}t|\}}t |||k r4JdS)Nr&ri)mask) r rir_rrTr5r r\r`raall)rbrlZ img_mask_eqrcrdrer"r"r#test_equalize_maskeds   rncCsDtddt|}t||d\}}d|krdks JJdS)z=Slope of cdf which should equal 1 for an equalized histogram.rr&rN皙?N)rrDr/Zpolyfit)rdZnorm_intensityZslopeZ interceptr"r"r#ras raztest_input,expectedr&rC)r.r.rpcC.tjddgtjd}t||d}t||dS)Nrr&r range_values)rrr:r rZ test_inputrUrWoutr"r"r#test_intensity_range_uint8 rv皙?皙?r)333333?皙?rzr{cCrq)Nrxryrrrrrfloat64r rrtr"r"r#test_intensity_range_floatrwr~cCs0tjddgtjd}t|ddd}t|ddS)NrxryrrT)rsZ clip_negative)rr&r|rWrur"r"r#"test_intensity_range_clipped_floatsriii?icCs@tjgdtjd}t|}|jtjksJt|gddS)N3r2rrrC)rrr:r rescale_intensityrrrr"r"r#test_rescale_stretch s rcCs*tgd}t|}t|gddS)NgI@gY@g c@)rrKr&rrr rrrr"r"r#test_rescale_shrinks rcCsFtjgd|d}tj|dd}t|gddd|j|jks!JdS)NrrrrCin_range)ryr{g333333?)decimal)rrr rrrrrWrur"r"r#test_rescale_in_rangesrcCs.tgd}tj|dd}t|gddS)Nr)rr2r)rKr&r&rrr"r"r#test_rescale_in_range_clip"srcCsFtjgd|d}tj|dd}|jt|jksJt|gddS)a?Check that output range is correct. .. versionchanged:: 0.17 This function used to return dtype matching the input dtype. It now matches the output. .. versionchanged:: 0.19 float16 and float32 inputs now result in float32 output. Formerly they would give float64 outputs. )rr.rrr out_range)rgO@rN)rrr rrrrrr"r"r#test_rescale_out_range(srcCs<tjdttdgtjd}tj|dd}t|dttgdS)Nrrruint10r)rr uint10_maxuint16r rr uint16_maxrr"r"r#test_rescale_named_in_range<srcC4tjdtgtjd}tj|dd}t|dtgdS)Nrrrr)rrrrr rrrrr"r"r#test_rescale_named_out_rangeBrcCr)Nrruint12r)rrrrr rr uint12_maxrr"r"r#test_rescale_uint12_limitsHrrcCr)NrrZuint14r)rrrrr rr uint14_maxrr"r"r#test_rescale_uint14_limitsNrrcCs<tjdtjd}t|}t|sJt||dS)NrSrSr)rrTr:r risnanrmrrr"r"r#test_rescale_all_zerosTs rcCs4tjddgtjd}tj|dd}t|ddgdS)Nrrrr)rrrr rrrr"r"r#test_rescale_constant[rrcCs6td}t|}t|sJt||dS)Nr)ronesr rrrmrrr"r"r#test_rescale_same_valuesas  rz1.25z,Older NumPy throws a few extra warnings herereasonzin_range,out_range)rWr)rrWcCs^tjdtddd}tj|d<tdgt|||WddS1s(wYdS)N rrMr)r&r&zWOne or more intensity levels are NaN\. Rescaling will broadcast NaN to the full image\.)rr)floatreshapenanrr r)rrrWr"r"r#test_rescale_nan_warninghs "rzout_range, out_dtyper:rrrrcCs4tjgdtjd}tj||d}|j|ksJdS)Nr=rrrr)rrr(r rr)rZ out_dtyperW output_imager"r"r#test_rescale_output_dtypezs rcCsFtjgdtjd}tj|tjd}t|gd|jtjks!JdS)Nrrrrr>rC)rrr(r rr:rrrWrr"r"r#test_rescale_no_overflowsrcCsBtjgdtjd}tj|dd}t|gd|jtksJdS)Nrrrrr)rrr(r rrrrrr"r"r#test_rescale_float_outputsrcCsRtjgdtjd}tttj|dd}WddS1s"wYdS)NrrZflatr)rrr(rrrr r)rW_r"r"r#*test_rescale_raises_on_incorrect_out_ranges "rcCsttj|dd}t|}t|||f}tj |dddd}|j |j ks)J|j t |ks2J|tj kr9dnd}tt||d |tt||d dd S) zTest a grayscale float imageFrg)9r{Gz?r>) kernel_size clip_limitrJrMrSg)\Y@ghs?N)r rir astronautrjr rdstackr equalize_adapthistr5rrfloat16rpeak_snrnorm_brightness_err)rrbadaptedZ snr_decimalr"r"r#test_adapthist_grayscalesrcCstt}tjdd}tdt|\}}t |dks"JWdn1s,wYtj |dd}| dks@J| dksHJ|j |j ksPJt|}tt||d d tt||d d dS) zTest an RGB color uint16 imageT)recordalwaysrNrr?x&Y[@r&{Gz?rS)r Z img_as_uintrrwarningscatch_warnings simplefilterr rr/rrOrPr5rrrr)rbwhistr!r full_scaler"r"r#test_adapthist_colors  rcCstt}tj|jd|jdftd}t||f}t |}|j|jks*J|ddddddf}t |}|j|jksDJt t ||ddt t||dddS) zTest an RGBA color imagerr&rNrMrrSga+e?)r rirrrrr5rrr rrrrr)rbalpharrr"r"r#test_adapthist_alphas  rcCstt}t|}d}t|dd|dd|f}t|g|jd|}tj |ddd}tj |ddd}|j|jksAJ|j|jksIJt t |||jdddks^Jd S) z Test for n-dimensional consistency with float images Note: Currently if img.ndim == 3, img.shape[2] > 4 must hold for the image not to be interpreted as a color image by @adapt_rgb rrr1g?)rrrSrN) r rirrr rrr5r rmeanabs)rbaZimg2dZimg3dZ adapted2dZ adapted3dr"r"r#test_adapthist_grayscale_Nds.rcCstd}|d7}|tj}t|d}t|t|ks!Jtd}|d7}|tj}t|d}t|t|ksBJdS)z#Test constant image, float and uintrrSrMrxN) rrTrjrr rrOrPr})rbrr"r"r#test_adapthist_constants      rcCsttt}|d}d||jdd|jddf<d}tdddD]4}tj||d d }t |d d |f|d d |fd ksDJt ||d d f||d d fd ksYJq%d S) zTest border processinggY@rrrSr&rrGrKrNrx) r r rirrr5r4r rr)rbZ border_indexrrr"r"r#test_adapthist_borderss *,rcCsbt}t|}tj|dd}tj|dd}t||tj|dd}tj|dd}t||dS)Nrrr&)rmoonr rir rr)Zimg_uZimg_fZ img_clahe0Z img_clahe1r"r"r#test_adapthist_clip_limit s  rcCsv|jdkrt|t|}}t|}t|}d|jt||}t |j j \}}dt ||S)zPeak signal to noise ratio of two images Parameters ---------- img1 : array-like img2 : array-like Returns ------- peak_snr : float Peak signal to noise ratio rMrrp) ndimr rhr risizerZsquaresumr rtypelog)img1img2ZmserZmax_r"r"r#rs   rcCsJ|jdkrt|t|}}t||}|t|jjd}|S)zNormalized Absolute Mean Brightness Error between two images Parameters ---------- img1 : array-like img2 : array-like Returns ------- norm_brightness_error : float Normalized absolute mean brightness error rMr&)rr rrrr rr)rrZambeZnber"r"r#r1s rcCsNtjdtd}tjtddt|dWddS1s wYdS)Nrrz Incorrect value of `kernel_size`r)rMrMrM)rrrrrrr rrbr"r"r#$test_adapthist_incorrect_kernel_sizeEs"rcCs.tddg}t|d}|j|jksJdS)"Check that the shape is maintainedr&g?N)rrr adjust_gammar5)rbresultr"r"r#test_adjust_gamma_1x1_shapeOs rcCs2tjddtjdd}t|d}t||dS)z4Same image should be returned for gamma equal to onerr8rrr&N)rr)r:rr rr)rWrr"r"r#test_adjust_gamma_oneVs rcCsTtjdddj|dd}t|d}|jj}t|t |d|j|jks(JdS)z6White image should be returned for gamma equal to zerorrCrFrgr&N) rrandomuniformrjr rrrrr )rrWrr"r"r#test_adjust_gamma_zero]s  rcCLtjddtjdd}tjgdtjdd}t|d}t||dS)z\Verifying the output with expected results for gamma correction with gamma equal to halfrr8rr(rr $'*-0257:<>@BDFrIKMNPQSTVWYrH\]^`abrr%r2gijklmopqrstuvwyz{|}~rr>rr4r?r=rCrrGrKrOrRrUrXr[r^r`rbrergrirkrmrorqrrrtrvrwryr{r|r~rrrrrrrrrrrr@r@rArArrrrrrrrrrrrrrrrrrrrrrrrCrKNrr)r:rrr rrrWrUrr"r"r#test_adjust_gamma_less_oneg& rcCr)z[Verifying the output with expected results for gamma correction with gamma equal to tworr8rr(rrrrrrrrrrrrr&r&r&r&r&r&r&r&rSrSrSrSrSrSrMrMrMrMrrrrr1r1r1r1rrrrrrr rrr.r. rrrr rrrrrrrrrrprprrrrrrrrrr!"#rr%&rr()rr+,r./rr1rr4r6r8rrr;r=r?rArCrErrr Jr Lr r Or rRrUrrXr[rrr_rrcrr2rhrrrrr r!r"r$r%r&r(xr)r+r,r.rr>rr0r2r3r5r6r7r8r:r;r=r>r@rArCrDrErFrHrIrKrMrNrPrQrSrUrVrXrZr[r]r_r`rbrdrergrirkrlrnrprrrsrurwryrzr|r~rrrrrrrrr@rrrrrrrrrrrrrrCrSNrrr"r"r#test_adjust_gamma_greater_onevrrcCsTtdddtjd}ttt|dWddS1s#wYdS)NrrCrrr) rr)r:rrrrr rrWr"r"r#test_adjust_gamma_negatives "rcCs6dtjdtjd}ttj|ddddksJdS)NrCrrr&ro)gammaZgain)rrr:rmr rrr"r"r#test_adjust_gamma_u8_overflows"rcCs@tjddg|d}t|d}|j|jksJ|j|ksJdS)rr&rN)rrr adjust_logr5rrrbrr"r"r#test_adjust_log_1x1_shapes rc Csrtdddtjd}tjgdgdgdgdgd gd gd gd gtjd }t|d}t||dS)z|Verifying the output with expected results for logarithmic correction with multiplier constant multiplier equal to unityrrCrr)rr1rrrrrr)rrrrrrr r )rrrrrrrr")r&r)r-r/r2r6r9r=)rArrHrKrOrSrWrZ)r^rarerhrlrorrrv)ryr|rrrrrr@)rrrrrrrrrr&Nrr)r:rrr rrrr"r"r#test_adjust_logs  rc Cttdddtjd}tjgdgdgdgdgd gd gd gd gtjd }t|dd}t||dS)zVerifying the output with expected results for inverse logarithmic correction with multiplier constant multiplier equal to unityrrCrr)rrSr1rrrrrp)rrrrrrrr)rrrrrrrH)rrrrrHrrr2)rnr#r'r*r.rr3)r?r;r?rCrGrKrPrU)rYr^rcrgrlrqrvr{)rrrrArrrrrr&TNrrr"r"r#test_adjust_inv_log rcCsBtjddg|d}t|dd}|j|jksJ|j|ksJdS)rr&rr1N)rrr adjust_sigmoidr5rrr"r"r#test_adjust_sigmoid_1x1_shapesrc Cr)zlVerifying the output with expected results for sigmoid correction with cutoff equal to one and gain of 5rrCrr)r&r&r&rSrSrSrSrS)rMrMrMrrrr1r1)r1rrrrrrr.)r.rrrrrrr)rrprrrrrr)rrrrrrrr)rrrrrrr r)rrrrrr"r'r+rr&r1Nrr)r:rrr rrrr"r"r#test_adjust_sigmoid_cutoff_onerrc Csttdddtjd}tjgdgdgdgdgd gd gd gd gtjd }t|dd }t||d S)znVerifying the output with expected results for sigmoid correction with cutoff equal to zero and gain of 10rrCrr)rr6r?rGrQrZrbrj)rqrxr~rrrrAr)rrrrrrrrrrrrrrrr)rrrrrrrr)rrrrrrrrrr.Nrrr"r"r#test_adjust_sigmoid_cutoff_zerorrc Cr)znVerifying the output with expected results for sigmoid correction with cutoff equal to half and gain of 10rrCrr)r&r&rSrSrMrMrr1)r1rrrr.rrr)rrrrrrrr)rrrr rrrr')r>r?r@rIrRr[rcrk)rrrxr~rrrrAr)rrrrrrrrrrrKr.Nrrr"r"r#test_adjust_sigmoid_cutoff_halfrrc Csvtdddtjd}tjgdgdgdgdgd gd gd gd gtjd }t|ddd}t||dS)zvVerifying the output with expected results for inverse sigmoid correction with cutoff equal to half and gain of 10rrCrr)rrrrrrrr)rrrrrrrr)rrrArrrr}rw)rprirarYrPrFr>r5)r.r%rrrr rr)rrrrrrrr)rrrrr.rrr)r1rrrMrMrSrSr&rrKr.TNrrr"r"r##test_adjust_inv_sigmoid_cutoff_halfs rcCstddd}t|sJd|d<t|sJtj|ddr"J|dtj}t|s1Jtj|ddr:J|tjd}t|sIJtj|ddrRJdS) Nrg{Gz?rr&r)Zupper_percentilerCr8)rrDr is_low_contrastrjr:rrr"r"r#test_is_low_contrast/srcCs:tjdtd}t|sJd|dd<t|rJdS)Nrrr&r1)rrTboolr rrr"r"r#test_is_low_contrast_boolean?s r exposure_funccCsTtddddtj}tt ||WddS1s#wYdS)Nrrrr)rr)rrjr}rrr)rrWr"r"r#test_negative_inputKs  "rcCstjdddddlm}|jtddgddggdd}t|\}}gd }gd }t||s4Jt||ss