\iOdfSrSSKrSSKJr SSKJr SSKJrJrJ r J r J r J r J r SSKJr SSKJr /SQr"S S 5r\"5r\R*r\R,r\R.r\R0rSS jrSS jrS rSrSrSrSr"SS5r \ RBr!SSjr"Sr#SSjr$g)z-Support for various GEOS geometry operations.N) polylabel)GeometryTypeError)GeometryCollection LineStringMultiLineString MultiPointPointPolygonshape) BaseGeometry)prep) clip_by_rect linemergenearest_pointsoperatororient polygonizepolygonize_full shared_pathssnapsplit substring transform triangulate unary_unionvalidatevoronoi_diagramc6\rSrSrSrSrSrS SjrSrSr g) CollectionOperator&c[U[5(aU$[U5$![[4a [ U5s$f=fN) isinstancer r ValueErrorAttributeErrorr)selfobs N/var/www/html/kml_chatgpt/mouzaenv/lib/python3.13/site-packages/shapely/ops.pyshapeupCollectionOperator.shapeup'sA b, ' 'I &Ry / &!"~% &s $AAcN[USS5=(d Un[U5nUVs/sHo0RU5PM nn[R "U5nUR $![a U/nNQf=fs snf!UVs/sHo0RU5PM nnfs snf=f)zCreate polygons from a source of lines. The source may be a MultiLineString, a sequence of LineString objects, or a sequence of objects than can be adapted to LineStrings. geomsN)getattriter TypeErrorr)shapelyrr,)r&linessourcelineobs collections r(rCollectionOperator.polygonize0s .7% :&\F399&$<<%&C9'',   XF :&9&$<<%&C9C9s4 A%A8% A52A=4A55A==B$BB$c6[USS5=(d Un[U5nUVs/sHo0RU5PM nn[R "U5$![a U/nNEf=fs snf!UVs/sHo0RU5PM nnfs snf=f)a@Create polygons from a source of lines. The polygons and leftover geometries are returned as well. The source may be a MultiLineString, a sequence of LineString objects, or a sequence of objects than can be adapted to LineStrings. Returns a tuple of objects: (polygons, cut edges, dangles, invalid ring lines). Each are a geometry collection. Dangles are edges which have one or both ends which are not incident on another edge endpoint. Cut edges are connected at both ends but do not form part of polygon. Invalid ring lines form rings which are invalid (bowties, etc). r,N)r-r.r/r)r0r)r&r1r2r3r4s r(r"CollectionOperator.polygonize_full@s .7% :&\F399&$<<%&C9&&s++  XF :&9&$<<%&C9C9s4 AA, A)&A1(A))A11B6BBcSn[USS5S:XaUnOx[US5(a0[URVs/sHoDRPM sn5nO7[US5(a&[UVs/sHoDRPM sn5nUc[ SU35e[R"X2S9$s snfs snf![ a [U5nNFf=f)zMerge all connected lines from a source. The source may be a MultiLineString, a sequence of LineString objects, or a sequence of objects than can be adapted to LineStrings. Returns a LineString or MultiLineString when lines are not contiguous. N geom_typerr,__iter__zCannot linemerge )directed) r-hasattrrr,coordsr%r$r0 line_merge)r&r1r<r2lss r(rCollectionOperator.linemergeYs 5+t ,0A AF UG $ $$%++%F+Bii+%FGF UJ ' ' 0(e)De))e)DE >089 9!!&<<&G*E! 0(/ 0s)B4* B>3B9B>9B>>CCc,[R"USS9$)zReturn the union of a sequence of geometries. Usually used to convert a collection into the smallest set of polygons that cover the same area. N)axis)r0 union_all)r&r,s r(rCollectionOperator.unary_unionos   T22NF) __name__ __module__ __qualname____firstlineno__r)rrrr__static_attributes__rGrFr(rr&s& ,2=,3rFrcV[R"XUS9n[UR5$)aCreate the Delaunay triangulation and return a list of geometries. The source may be any geometry type. All vertices of the geometry will be used as the points of the triangulation. From the GEOS documentation: tolerance is the snapping tolerance used to improve the robustness of the triangulation computation. A tolerance of 0.0 specifies that no snapping will take place. If edges is False, a list of Polygons (triangles) will be returned. Otherwise the list of LineString edges is returned. ) tolerance only_edges)r0delaunay_triangleslistr,)geomrOedgesr5s r(rrs(++DRWXJ    !!rFc[R"XXS9nURS:wa [ U/5$U$![Ra)nSnUSU<S3- nU(aUS- n[U5UeSnAff=f)aConstruct a Voronoi Diagram [1] from the given geometry. Returns a list of geometries. Parameters ---------- geom: geometry the input geometry whose vertices will be used to calculate the final diagram. envelope: geometry, None clipping envelope for the returned diagram, automatically determined if None. The diagram will be clipped to the larger of this envelope or an envelope surrounding the sites. tolerance: float, 0.0 sets the snapping tolerance used to improve the robustness of the computation. A tolerance of 0.0 specifies that no snapping will take place. edges: bool, False If False, return regions as polygons. Else, return only edges e.g. LineStrings. GEOS documentation can be found at [2] Returns ------- GeometryCollection geometries representing the Voronoi regions. Notes ----- The tolerance `argument` can be finicky and is known to cause the algorithm to fail in several cases. If you're using `tolerance` and getting a failure, try removing it. The test cases in tests/test_voronoi_diagram.py show more details. References ---------- [1] https://en.wikipedia.org/wiki/Voronoi_diagram [2] https://geos.osgeo.org/doxygen/geos__c_8h_source.html (line 730) )rO extend_torPz;Could not create Voronoi Diagram with the specified inputs (z).z0 Try running again with default tolerance value.Nr)r0voronoi_polygons GEOSExceptionr$r:r)rSenveloperOrTresulterrerrstrs r(rrsV *))  //!6(++ M  *NAcWB-  H HF c) *s5A2 $A--A2c.[R"U5$)z%Return True if the geometry is valid.)r0is_valid_reason)rSs r(rrs  " "4 ((rFc JUR(aU$URS;aURS;a*[U5"[U"[UR6665$URS:Xa[UR 5"[U"[UR R6665nUR Vs/sH-n[U5"[U"[UR6665PM/ nn[U5"X$5$gURRS5(dURS:Xa5[U5"URVs/sHn[X5PM sn5$[SUR<S35es snf![Ga URS;a6[U5"URVs/sHoP"U6PM Os snfsn5s$URS:Xa[UR 5"UR RVs/sHoP"U6PM Os snfsn5nUR VVs/sH7n[U5"URVs/sHoP"U6PM Os snfsn5PM9 Os snnfnnn[U5"X$5s$gf=fs snf) aApply `func` to all coordinates of `geom`. Returns a new geometry of the same type from the transformed coordinates. `func` maps x, y, and optionally z to output xp, yp, zp. The input parameters may iterable types like lists or arrays or single values. The output shall be of the same type. Scalars in, scalars out. Lists in, lists out. For example, here is an identity function applicable to both types of input. def id_func(x, y, z=None): return tuple(filter(None, [x, y, z])) g2 = transform(id_func, g1) Using pyproj >= 2.1, this example will accurately project Shapely geometries: import pyproj wgs84 = pyproj.CRS('EPSG:4326') utm = pyproj.CRS('EPSG:32618') project = pyproj.Transformer.from_crs(wgs84, utm, always_xy=True).transform g2 = transform(project, g1) Note that the always_xy kwarg is required here as Shapely geometries only support X,Y coordinate ordering. Lambda expressions such as the one in g2 = transform(lambda x, y, z=None: (x+1.0, y+1.0), g1) also satisfy the requirements for `func`. )r r LinearRingr )r rrar MultirzType z not recognizedN) is_emptyr:typezipr>exterior interiorsr/ startswithr,rr)funcrSshellringholescparts r(rrs/L }}  ~~II  0~~!FFDz#tS$++->'?"@AA9,T]]+CsDMM .Dz%// -*  " "7 + +t~~AU/UDzTZZHZT9T0ZHII%'9 IJJ- 0~~!FFDzT[["A[48["ABB9,T]]+t}}?S?S,T?S!T1X?S,TU!% .J$++>+Qa+>? .Dz%// - 0Isf9FAF;4F/F J F4J;G J=JH$ #J?I>I-,I>=JJc[R"X5nUc'UR(a [S5e[S5e[R"US5n[R"US5nX44$)zReturn the calculated nearest points in the input geometries. The points are returned in the same order as the input geometries. z!The first input geometry is emptyz"The second input geometry is emptyr)r0 shortest_linercr$ get_point)g1g2seqp1p2s r(rrsc    'C { ;;@A AAB B   3 "B   3 "B 8OrFc0[R"XU5$)a Snaps an input geometry (g1) to reference (g2) geometry's vertices. Parameters ---------- g1 : geometry The first geometry g2 : geometry The second geometry tolerance : float The snapping tolerance Refer to :func:`shapely.snap` for full documentation. )r0r)rsrtrOs r(rr.s << **rFc[U[5(d [S5e[U[5(d [S5e[R"X5$)aFind paths shared between the two given lineal geometries. Returns a GeometryCollection with two elements: - First element is a MultiLineString containing shared paths with the same direction for both inputs. - Second element is a MultiLineString containing shared paths with the opposite direction for the two inputs. Parameters ---------- g1 : geometry The first geometry g2 : geometry The second geometry z#First geometry must be a LineStringz$Second geometry must be a LineString)r#rrr0r)rsrts r(rr@sG" b* % % EFF b* % % FGG    ''rFcd\rSrSr\S5r\S5r\S5r\S5r\S5r Sr g) SplitOpiXcf[U[5(d [S5e[U[[45(d [S5eUR R U5n[U5n[U5Vs/sH*o0RUR55(dM(UPM, sn$s snf)z"Split a Polygon with a LineString.z First argument must be a Polygonz+Second argument must be a (Multi)LineString) r#r rrrboundaryunionr rcontainsrepresentative_point)polysplitterr~pgs r(_split_polygon_with_line SplitOp._split_polygon_with_lineYs$((#$FG G(Z$ABB#$QR R ##H- Dz $E* *2mmBprojectrRrangelen) r3rdistance_on_liner>current_positionipoint1point2dxdysegment_lengths r(_split_line_with_pointSplitOp._split_line_with_points$ ++#$IJ J(E**#$EF F""8[99 6M [[^xq1 16M <<1dkk"s6{Q'AYFE]FVAY&BVAY&B !eb!em3N  . 3"6'AE?3Z1uw5PQQ!4vgA(//!2D1EEF 23fUWoEF( v rFc*[U[5(d [S5e[U[5(d [S5eU/nURH?n/n[ SU5H'nUR [RXS55 M) UnMA U$)z%Split a LineString with a MultiPoint.rz$Second argument must be a MultiPointc$UR(+$r")rc)xs r(5SplitOp._split_line_with_multipoint..s ajj.rF) r#rrrr,filterextendr{r)r3rchunkspt new_chunkschunks r(_split_line_with_multipoint#SplitOp._split_line_with_multipoints$ ++#$IJ J(J//#$JK K..BJ 8&A!!'"@"@"KLB F ! rFc URS;aN[URVVs/sH+n[R X!5RHo3PM M- snn5$URS:Xa|URS;a[R nOURS:Xa[R nOURS:Xa[RnO|[SURS35eURS:Xa:URS ;a[RnO2[S URS35e[S URS 35e[U"X55$s snnf) aSplit a geometry by another geometry and return a collection of geometries. This function is the theoretical opposite of the union of the split geometry parts. If the splitter does not split the geometry, a collection with a single geometry equal to the input geometry is returned. The function supports: - Splitting a (Multi)LineString by a (Multi)Point or (Multi)LineString or (Multi)Polygon - Splitting a (Multi)Polygon by a LineString It may be convenient to snap the splitter with low tolerance to the geometry. For example in the case of splitting a line by a point, the point must be exactly on the line, for the line to be correctly split. When splitting a line by a polygon, the boundary of the polygon is used for the operation. When splitting a line by another line, a ValueError is raised if the two overlap at some segment. Parameters ---------- geom : geometry The geometry to be split splitter : geometry The geometry that will split the input geom Examples -------- >>> import shapely.ops >>> from shapely import Point, LineString >>> pt = Point((1, 1)) >>> line = LineString([(0,0), (2,2)]) >>> result = shapely.ops.split(line, pt) >>> result.wkt 'GEOMETRYCOLLECTION (LINESTRING (0 0, 1 1), LINESTRING (1 1, 2 2))' )rrr)rrr rr rzSplitting a LineString with a z is not supportedr )rrzSplitting a Polygon with a z Splitting z geometry is not supported) r:rr,r{rrrrrr)rSrrnr split_funcs r(r SplitOp.splitsVN >>@ @%#zzWzt7==3P3V3Va3VzW ^^| +!!& %:: ##w.$;; ##|3$@@ '4X5G5G4HI$$ ^^y (!!%FF$== '1(2D2D1EEVW $T^^,,FG "*T"<==EXs2E rGN) rIrJrKrL staticmethodrrrrrrMrGrFr(r{r{Xsh  *B''R"J>J>rFr{c[U[5(d[SURS35eX:XaUR XS9$U(d8XR :a)X R :aUR UR US9$U(d2U*UR :a!U*UR :aUR SUS9$U(aUS:aUS:aUR SUS9$U(aU*S:aU*S:aUR SUS9$U(aXR -nX R -nUSs=:aU:a0O O-[ U5U-UR :XaUR U5$USs=:aU:a0O O-[ U5U-UR :XaUR U5$UR U5nUR U5nUS:aUR U-nUS:aUR U-nX:nU(aX!p![US5nU(a[UR6/nO[UR6/n[UR5nSn [XSS5HOupXs=:aU:aO OURU 5 OX:a O&XSU S- S-U SU S- S--S-- n MQ U(a-UR[UR65 [U5nO!UR[UR65 [U5$) aReturn a line segment between specified distances along a LineString. Negative distance values are taken as measured in the reverse direction from the end of the geometry. Out-of-range index values are handled by clamping them to the valid range of values. If the start distance equals the end distance, a Point is returned. If the start distance is actually beyond the end distance, then the reversed substring is returned such that the start distance is at the first coordinate. Parameters ---------- geom : LineString The geometry to get a substring of. start_dist : float The distance along `geom` of the start of the substring. end_dist : float The distance along `geom` of the end of the substring. normalized : bool, False Whether the distance parameters are interpreted as a fraction of the geometry's length. Returns ------- Union[Point, LineString] The substring between `start_dist` and `end_dist` or a Point if they are at the same location. Raises ------ TypeError If `geom` is not a LineString. Examples -------- >>> from shapely.geometry import LineString >>> from shapely.ops import substring >>> ls = LineString((i, 0) for i in range(6)) >>> ls.wkt 'LINESTRING (0 0, 1 0, 2 0, 3 0, 4 0, 5 0)' >>> substring(ls, start_dist=1, end_dist=3).wkt 'LINESTRING (1 0, 2 0, 3 0)' >>> substring(ls, start_dist=3, end_dist=1).wkt 'LINESTRING (3 0, 2 0, 1 0)' >>> substring(ls, start_dist=1, end_dist=-3).wkt 'LINESTRING (1 0, 2 0)' >>> substring(ls, start_dist=0.2, end_dist=-0.6, normalized=True).wkt 'LINESTRING (1 0, 2 0)' Returning a `Point` when `start_dist` and `end_dist` are at the same location. >>> substring(ls, 2.5, -2.5).wkt 'POINT (2.5 0)' z;Can only calculate a substring of LineString geometries. A z was provided.) normalizedrrpNrr)r#rrr: interpolatelengthabsmaxtupler>rRreappendreversed) rS start_distend_distr start_point end_pointreverse vertex_listr>current_distancervrws r(rrsv dJ ' '   0   BB J++5(kk:Q CC ZK4;;6H9 ;Sj99  aHMj99  q(hY!^j99kk! KKA  S_x%?4;;%N)) A " "s8}z'AT[['P ++"":.K  *IA~[[:- !|;;)#G'HZ#Ji../0 k0012 $++ FfQRj)  38 3   r "  ) ebem1RURU]q4HHSPP *5+"4"456{+ 5)"2"234 k ""rFcXUR(aU$[R"XX#U5$)aReturn the portion of a geometry within a rectangle. The geometry is clipped in a fast but possibly dirty way. The output is not guaranteed to be valid. No exceptions will be raised for topological errors. Parameters ---------- geom : geometry The geometry to be clipped xmin : float Minimum x value of the rectangle ymin : float Minimum y value of the rectangle xmax : float Maximum x value of the rectangle ymax : float Maximum y value of the rectangle Notes ----- New in 1.7. )rcr0r)rSxminyminxmaxymaxs r(rrs%2 }}   D ==rFc0[R"XS:S9$)aReturn a properly oriented copy of the given geometry. The signed area of the result will have the given sign. A sign of 1.0 means that the coordinates of the product's exterior rings will be oriented counter-clockwise. It is recommended to use :func:`shapely.orient_polygons` instead. Parameters ---------- geom : Geometry The original geometry. May be a Polygon, MultiPolygon, or GeometryCollection. sign : float, optional. The sign of the result's signed area. Returns ------- Geometry r) exterior_cw)r0orient_polygons)rSsigns r(rrs,  " "4AX >>rF)rF)NrFrH)g?)%__doc__r0shapely.algorithms.polylabelrshapely.errorsrshapely.geometryrrrrr r r shapely.geometry.baser shapely.preparedr __all__rrrrrrrrrrrrrr{rrrrrGrFr(rs32,/! (O3O3d    **    "" "&8v) HKV"+$(0@>@>F  {#|><?rF