6Dh6X (SSKJr SSKJr SSKJr SSKJr \R"S5r \R"S5r SSK J r SS S /S Q4S jrSS jrSrSrSrSrSSjrSrS SjrSrSrSrS!Sjr/S QSSSSSSSSSS4 Sjrg)"N) graph_objs) exceptions)optional_importsnumpyzscipy.interpolate)measurezTernary contour ploti&i )abcc[UUU[S[[USS9SSSS9[[USS9SSSS9[[USS9SSSS9S9S S 9$) a) Layout of ternary contour plot, to be passed to ``go.FigureWidget`` object. Parameters ========== title : str or None Title of ternary plot width : int Figure width. height : int Figure height. pole_labels : str, default ['a', 'b', 'c'] Names of the three poles of the triangle. r)text{Gz?outside)titlemin linewidthticks)sumaaxisbaxiscaxisF)rwidthheightternary showlegend)dict)rrr pole_labelss X/var/www/html/env/lib/python3.13/site-packages/plotly/figure_factory/_ternary_contour.py_ternary_layoutr s|$  A/TQi A/TQi A/TQi  ! cUS:Hn[RUSS9nSX-- n[RU5S:a [S5e[RX!X@-5nU$)a Replaces zero ternary coordinates with delta and normalize the new triplets (a, b, c). Parameters ---------- ternary_data : ndarray of shape (N, 3) delta : float Small float to regularize logarithm. Notes ----- Implements a method by J. A. Martin-Fernandez, C. Barcelo-Vidal, V. Pawlowsky-Glahn, Dealing with zeros and missing values in compositional data sets using nonparametric imputation, Mathematical Geology 35 (2003), pp 253-278. raxisr zMThe provided value of delta led to negativeternary coords.Set a smaller delta)npany ValueErrorwhere) ternary_datadelta zero_maskis_any_coord_zerounity_complements r_replace_zero_coordsr.7so*!Iyq1544 vv!# 1  88I.>.MNL r!cD[RU5n[RUSUS- 5[RS5- nS[RS5- [RUSUS-USS-- 5-n[R X45nU$)a Perform Isometric Log-Ratio on barycentric (compositional) data. Parameters ---------- barycentric: ndarray of shape (3, N) Barycentric coordinates. References ---------- "An algebraic method to compute isometric logratio transformation and back transformation of compositional data", Jarauta-Bragulat, E., Buenestado, P.; Hervada-Sala, C., in Proc. of the Annual Conf. of the Intl Assoc for Math Geology, 2003, pp 31-30. rr r?)r%asarraylogsqrtstack) barycentricx_0x_1 ilr_tdatas r_ilr_transformr:Ys **[)K &&Q+a.0 1BGGAJ >C bggaj266+a.;q>"AKPQNVWDW"WXX#$I r!c 8[RU5n[R/SQ/SQ/SQ/5n[RS5S- n[RS5n[RS[RX#/U*U//5U5n[R S[R U5R SS 9-5*n[RSU[RSU-US-X0S-U455n[R U5nU$) a Perform inverse Isometric Log-Ratio (ILR) transform to retrieve barycentric (compositional) data. Parameters ---------- x : array of shape (2, N) Coordinates in ILR space. References ---------- "An algebraic method to compute isometric logratio transformation and back transformation of compositional data", Jarauta-Bragulat, E., Buenestado, P.; Hervada-Sala, C., in Proc. of the Annual Conf. of the Intl Assoc for Math Geology, 2003, pp 31-30. )?r r0)gr r0)r=r0rg? ik, kj -> ijr rr#)r%arrayr4einsumr3exprr5)xmatrixstSkZlog_barycentric iilr_tdatas r _ilr_inverserJrs"  A XX}noF GF  QA A >288aVqb!W,=#> BB BFF2JNNN** ++Aii!a%!A$,aD!)D EO(J r!c[RS[RS5S- /SS/SS//5n[RUSS2S4USS2S4[RS5/5nU[RR U54$)z] Returns the transformation matrix from barycentric to Cartesian coordinates and conversely. r<rrr N)r%r?r4oneslinalginv) tri_vertsMs r _transform_barycentric_cartesianrRsw 3 Q/!Q!Q@AI )AqD/9QT?BGGAJ?@A biimmA r!c[U[[R45(d [ S5e[R U5nUR SS;a [ S5e[U5S:Xa[[RURSS9SSS 9(d4[RURSS9S SS 9(d S n[ U5e[U5S :Xa Uup#SX#-- nOXRSS9- up#n[R[RX#U45S:5(a [ S 5e[RX#U45$)zI Check ternary coordinates and return the right barycentric coordinates. zOData should be either an array of shape (n,m),or a list of n m-lists, m=2 or 3r)rrLzCA point should have 2 (a, b) or 3 (a, b, c)barycentric coordinatesrLr#r r)rtoldz=The sum of coordinates should be 1 or 100 for all data pointsrz+Barycentric coordinates should be positive.) isinstancelistr%ndarrayr'r2shapelenallcloserr&r5)b_coordsmsgABCs r_prepare_barycentric_coordras) hrzz 2 3 3 /  zz(#H~~a& T   X!  HLLaL0!$ ? HLLaL0#D AMo 8} K\\q\11a vvbhhay!A%&&FGG 88Q1I r!ilrc US:Xa#[5up4[RSX05nO(US:Xa[U5n[ U5nO [ S5eUSSupgUR 5UR5pUR 5UR5p[S[[R[U5555n [RXU 5n [RXU 5n[RX5unn[RUSSRXU4SS 9n[RUSSRXU4S S 9nUX4$) aw Transform data points with Cartesian or ILR mapping, then Compute interpolation on a regular grid. Parameters ========== coordinates : array-like Barycentric coordinates of data points. values : 1-d array-like Data points, field to be represented as contours. interp_mode : 'ilr' (default) or 'cartesian' Defines how data are interpolated to compute contours. cartesianr>rbz&interp_mode should be cartesian or ilrNrcubic)methodnearest)rRr%r@r.r:r'rmaxintr4rZlinspacemeshgrid scipy_interpgriddataT) coordinatesvalues interp_moderQinvM coord_pointsxxyyx_minx_maxy_miny_maxn_interpgr_xgr_ygrid_xgrid_ygrid_z grid_z_others r _compute_gridrs?k!24yy@  *;7 %k2 ABB "1 FB668RVVX5668RVVX53BGGCK012H ;;uX .D ;;uX .D[[,NFF " "RaFV$4W#F ((RaFV$4Y)L 4 r!c S[R[RU[RUS55[RU[RUS55- 5-$)Nr<r )r%absdotroll)rBys r _polygon_arearsE q"''!Q-0266!RWWQ]3KKL LLr!cVU[RR5;a[RUnOA[R"SR [RR555e[ RSSU5n[ RUVs/sHoDSPM sn5n[ RUVs/sHoDSPM sn5n[ RXS5nSUS;a:UVs/sH-n[R"[R"U55PM/ nnUS/n [USSUSS5HGupXZS- XZpX- X- - n[R"XjS- XjUSS9nU RU5 MI U $s snfs snfs snf)zI Return a list of ``ncontours`` colors from the ``colormap`` colorscale. zRColorscale must be a valid Plotly Colorscale.The available colorscale names are {}rr #Nrgb) colortype)clrs PLOTLY_SCALESkeysr PlotlyErrorformatr%rkr? searchsorted label_rgb hex_to_rgbzipfind_intermediate_colorappend) ncontourscolormapcmaprqpair vals_cmapcolsindscolcolorsindvalval1val2interms r_colorsrsy4%%**,,!!(+$$ 44:F4;M;M;R;R;T4U  [[Ay )Fd3ddq'd34I 88.!W. /D ??9 -D d1g~@DEts34E1gYFQR&*-Qw'd*-** qM49f   c . M4.Fs"F F!4F&c[R[RXS- 5S:5=(a0 [R[RXS- 5S:5nU$)zm Utility function for _contour_trace Contours with an area of the order as 1 pixel are considered spurious. rr)r%allr)rBr too_smalls r_is_invalid_contourr sN rvvaA$h'!+,MqQ4x8H18L1MI r!c [RU5nU[RU5R5U[RU5R 5pT[R U5nSU-Xc'[R U5nSU-Xs'////4upp////4upp[ U5GHaunn[R"UU5n[R"UU5nURU5 U RU5 U RU/[U5-5 U RU/[U5-5 U RUVs/sH)n[URSURS5PM+ sn5 URUVs/sH)n[URSURS5PM+ sn5 U RUU/[U5-5 URUU/[U5-5 GMd [U5[U 5::aXX4$XX4$s snfs snf)aC Utility function for _contour_trace. In ``im`` only one part of the domain has valid values (corresponding to a subdomain where barycentric coordinates are well defined). When computing contours, we need to assign values outside of this domain. We can choose a value either smaller than all the values inside the valid domain, or larger. This value must be chose with caution so that no spurious contours are added. For example, if the boundary of the valid domain has large values and the outer value is set to a small one, all intermediate contours will be added at the boundary. Therefore, we compute the two sets of contours (with an outer value smaller of larger than all values in the valid domain), and choose the value resulting in a smaller total number of contours. There might be a faster way to do this, but it works... rr r) r%isnan logical_notrricopy enumerater find_contoursextendrZrro)imrqrmask_nanim_minim_maxzz_minzz_max all_contours1 all_values1 all_areas1 all_colors1 all_contours2 all_values2 all_areas2 all_colors2ircontour_level1contour_level2contours r_extract_contoursrs $xx|H 2>>( #$((* 2>>( #$((* WWR[F6zF WWR[F6zF: ?F JJv* +A .q 1E s6{Q    -:-J-:-JiI&,DfU0C5fKD{#L KK1lQ. /EKk*5E!Hk!n56E!a%L+a.9:+l+[9I-JKL YK GGEs  F&Fc ([US-U5n [RXUS-5n U SU SpU SSn U SSn UcSnOU/U-n [X+U 5upnn[R U5SSS2n[ UUUUU UUSUU U U U 5 upnnn[R S/US-45n/n[5unnUR5UR5- UR- nUR5UR5- UR- nUGHpnUURunnUUnUS:XaP[RU[RUU-UU-[RUR5455nONUS:XaH[![RUU-UR5-UU-UR5-455nUS:XaF[R#/S Q5n[R#/S Q5n[R#/S Q5n OWunnn [%UU5(aGM%US :XaUUOUn!['S UUU S ['U!SSS9SUUSSSU-S9 n"US :XaSU"S'UR)U"5 GMs UU4$)a Contour trace in Cartesian coordinates. Parameters ========== x, y : array-like Cartesian coordinates z : array-like Field to be represented as contours. ncontours : int or None Number of contours to display (determined automatically if None). colorscale : None or str (Plotly colormap) colorscale of the contours. linecolor : rgb color Color used for lines. If ``colorscale`` is not None, line colors are determined from ``colorscale`` instead. interp_mode : 'ilr' (default) or 'cartesian' Defines how data are interpolated to compute contours. If 'irl', ILR (Isometric Log-Ratio) of compositional data is performed. If 'cartesian', contours are determined in Cartesian space. coloring : None or 'lines' How to display contour. Filled contours if None, lines if ``lines``. vmin, vmax : float Bounds of interval of values used for the colorspace Notes ===== rrrr Nzrgb(150, 150, 150)rdrb)r rr)rr rrlinesscatterternaryspline)rrYrtoselfTskipz%.3f) typerr r modelinefill fillcolorr hoverinfonamer)rr%rkrargsortrrrRrirsizerorr5rMrYrJr?rrr)#rBrzr colorscale linecolorrrcoloringrrrrqrrrrrrorderrtracesrQrsdxdyr y_contour x_contourr bar_coordsrr r _coltraces# r_contour_tracersXY]J /F[[y1} 5F!!9fRjy Ab\F Ab\F( y(7H 673Li JJy !$B$ 'EDV58   D@Li[ NNQC+ ,EF.0GAt %%'AEEG qvv %B %%'AEEG qvv %B+E244 9 + %"y."y."'')//:RSTJE !%"y.15572BNQUUW4LMNJ A:#A#A#A GAq! y) 4 4 $,$7z% Y!D: '#   w  E&M eKN ; r!iFBlueredc [c [S5e[R"S5n U c [S5eUcSnUcSn[ U5nUR 5UR 5p[XU S9unnn[X#XJS9n[UUUUUU U UUUS 9 unn[R"UUS 9nU (aS OS nUunnnUS S -US -S-US-S-nURUUUSUUS[US:g5S.S.UUS9 U(aPU (dUn[SS/S/S/UR 5UR 5USS.SS.5nURU5 U$)a Ternary contour plot. Parameters ---------- coordinates : list or ndarray Barycentric coordinates of shape (2, N) or (3, N) where N is the number of data points. The sum of the 3 coordinates is expected to be 1 for all data points. values : array-like Data points of field to be represented as contours. pole_labels : str, default ['a', 'b', 'c'] Names of the three poles of the triangle. width : int Figure width. height : int Figure height. ncontours : int or None Number of contours to display (determined automatically if None). showscale : bool, default False If True, a colorbar showing the color scale is displayed. coloring : None or 'lines' How to display contour. Filled contours if None, lines if ``lines``. colorscale : None or str (Plotly colormap) colorscale of the contours. linecolor : None or rgb color Color used for lines. ``colorscale`` has to be set to None, otherwise line colors are determined from ``colorscale``. title : str or None Title of ternary plot interp_mode : 'ilr' (default) or 'cartesian' Defines how data are interpolated to compute contours. If 'irl', ILR (Isometric Log-Ratio) of compositional data is performed. If 'cartesian', contours are determined in Cartesian space. showmarkers : bool, default False If True, markers corresponding to input compositional points are superimposed on contours, using the same colorscale. Examples ======== Example 1: ternary contour plot with filled contours >>> import plotly.figure_factory as ff >>> import numpy as np >>> # Define coordinates >>> a, b = np.mgrid[0:1:20j, 0:1:20j] >>> mask = a + b <= 1 >>> a = a[mask].ravel() >>> b = b[mask].ravel() >>> c = 1 - a - b >>> # Values to be displayed as contours >>> z = a * b * c >>> fig = ff.create_ternary_contour(np.stack((a, b, c)), z) >>> fig.show() It is also possible to give only two barycentric coordinates for each point, since the sum of the three coordinates is one: >>> fig = ff.create_ternary_contour(np.stack((a, b)), z) Example 2: ternary contour plot with line contours >>> fig = ff.create_ternary_contour(np.stack((a, b, c)), z, coloring='lines') Example 3: customize number of contours >>> fig = ff.create_ternary_contour(np.stack((a, b, c)), z, ncontours=8) Example 4: superimpose contour plot and original data as markers >>> fig = ff.create_ternary_contour(np.stack((a, b, c)), z, coloring='lines', ... showmarkers=True) Example 5: customize title and pole labels >>> fig = ff.create_ternary_contour(np.stack((a, b, c)), z, ... title='Ternary plot', ... pole_labels=['clay', 'quartz', 'fledspar']) NzH The create_ternary_contour figure factory requires the scipy packageskimagezS The create_ternary_contour figure factory requires the scikit-image packageF)rr)rrrr)rrrrrrrr)datalayoutr rz: %{a:.3f}
z: %{b:.3f}
rz3: %{c:.3f}
z: %{marker.color:.3f}markerszrgb(120, 120, 120)r)rr)rrr)rr r rmarkeropacity hovertemplaterT)cmincmaxr showscale)rrr r rr)rm ImportErrorr get_modulerarrirr rgoFigureadd_scatterternaryrjr add_trace)rprqrrrrrrrrrrr showmarkers sk_measurerrrr|r}r contour_tracerfigrrr r rcolorbars rcreate_ternary_contourrsB L  ",,Y7J     ,[9K::<5&{ TFD$ VF"0   "M; ))v 6CaAGGAq!A   a.    a.    0  0 $2SWAT=UV # $J(VVV"JJL"JJL",!%  "   h Jr!)gMb@?)rb)N)NElectriczrgb(150,150,150)llrNrr ) plotly.colorsrrplotly.graph_objsrr plotlyrrr r%rmrrr r.r:rJrRrarrrrrrrrr!rrs.#  )**+>?  !C_#RD2<@(\M:/Cd6Hz   P   ur!