(&hOd@dZddlZddlmZddlmZddlmZmZm Z m Z m Z m Z m Z ddlmZddlmZgdZGd d ZeZejZejZejZejZdd ZddZdZdZdZdZdZGddZ e j!Z!ddZ"dZ#ddZ$dS)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_diagramc.eZdZdZdZdZddZdZdS) CollectionOperatorct|tr|S t|S#ttf$rt |cYSwxYwN) isinstancer r ValueErrorAttributeErrorr)selfobs F/var/www/html/reinick/venv/lib/python3.11/site-packages/shapely/ops.pyshapeupzCollectionOperator.shapeup's\ b, ' ' &I &Ryy / & & &!"~~%%% &s( A  A ct|ddp|} t|}n#t$r|g}YnwxYwfd|D}n#fd|D}wxYwtj|}|jS)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. geomsNc:g|]}|Sr(.0liner%s r' z1CollectionOperator.polygonize..<%999$4<<%%999)getattriter TypeErrorshapelyrr*)r%linessourceobs collections` r'rzCollectionOperator.polygonize0s ..7% :&\\FF   XFFF :999&999CC9999&999C 9 9 9 9',, !&A 6A6AAct|ddp|} t|}n#t$r|g}YnwxYwfd|D}n#fd|D}wxYwtj|S)aCreate 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*Nc:g|]}|Sr,r-r.s r'r1z6CollectionOperator.polygonize_full..Vr2r3)r4r5r6r7r)r%r8r9r:s` r'rz"CollectionOperator.polygonize_full@s ..7% :&\\FF   XFFF :999&999CC9999&999C 9 9 9 9&s+++r<Fc|d}t|dddkr|}nyt|drtd|jD}nJt|dr: td|D}n#t$rt|}YnwxYw|t d|t j|| S) aMerge 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*cg|] }|j Sr,coordsr/lss r'r1z0CollectionOperator.linemerge..es%F%F%FBbi%F%F%Fr3__iter__cg|] }|j Sr,rBrDs r'r1z0CollectionOperator.linemerge..hs)D)D)D"))D)D)Dr3zCannot linemerge )directed)r4hasattrrr*r$r#r7 line_merge)r%r8rHr9s r'rzCollectionOperator.linemergeYs 5+t , ,0A A AFF UG $ $ 0$%F%F%+%F%F%FGGFF UJ ' ' 0 0()D)De)D)D)DEE! 0 0 0(// 0 >88899 9!&8<<<= 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 rrfr c lg|]0}t|tt|j1Sr,)typeziprCr/ringfuncs r'r1ztransform..sLDJJsDD#t{*;$<=>>r3cg|]}|Sr,r,r/crls r'r1ztransform..s"A"A"A448"A"A"Ar3cg|]}|Sr,r,rns r'r1ztransform..s,T,T,T!TT1X,T,T,Tr3cbg|]+}t|fd|jD,S)cg|]}|Sr,r,rns r'r1z(transform...s>>>Qa>>>r3)rhrCrjs r'r1ztransform..sPDJJ>>>>$+>>>??r3Multirc0g|]}t|Sr,)r)r/partrls r'r1ztransform..s#HHHT9T400HHHr3zType z not recognizedN) is_emptyr@rhrirCexterior interiorsr6 startswithr*r)rlrXshellholess` r'rrsKL }  ~III  0~!FFF!tDzz#ttS$+->'?"@AAA9,,+T]++CsDM A.B>>.ksF   dmmBz5SplitOp._split_line_with_multipoint..s aj.r3) r"rrrr*filterextendrr)r0rchunkspt new_chunkschunks r'_split_line_with_multipointz#SplitOp._split_line_with_multipoints$ ++ K#$IJJ J(J// L#$JKK K.  BJ 8 8&AA M M!!'"@"@"K"KLLLLFF r3c|jdvr tfd|jDS|jdkr^jdvr tj}njdkr tj}njdkr tj}nitdjd|jd kr.jd vr tj}n0td jdtd |jd t||S)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))' )rrcZg|]'}t|jD]}|(Sr,)rrr*)r/rurrs r'r1z!SplitOp.split..s7WWWt7==x3P3P3VWWaWWWWr3r)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*rrrrrr)rXr split_funcs ` r'rz SplitOp.splitsUN >@ @ @%WWWWtzWWW ^| + +!& %: #w..$; #|33$@ '$X5G$$$ ^y ( (!%FFF$= 'W(2DWWW $GT^GGG "**T8"<"<===r3N) rOrPrQ staticmethodrrrrrr,r3r'rrXs  \ *\B''\'R\"J>J>\J>J>J>r3rct|tstd|jd||kr|||S|s2||jkr'||jkr||j|S|s/| |jkr#| |jkr|d|S|r#|dkr|dkr|d|S|r%| dkr| dkr|d|S|r||jz}||jz}|dcxkr|kr3nn0t ||z|jkr||S|dcxkr|kr3nn0t ||z|jkr||S||}||}|dkr |j|z}|dkr |j|z}||k}|r||}}t|d}|rt|j g}nt|j g}t|j }d} t||ddD]b\} } || cxkr|krnn| | n| |krn0| | d| dz dz| d| dz dzzdzz } c|r1| t|j t|}n!| t|j t|S) aaReturn 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.) normalizedrr}Nrr)r"rrr@ interpolatelengthabsmaxtuplerCrWriappendreversed) rX start_distend_distr start_point end_pointreverse vertex_listrCcurrent_distancerrs r'rrsv dJ ' '  0 0 0 0   X zBBB :J$+55(dk:Q:Q CCC :ZK4;66H9 ;S;Sj999 : aHMMj999 : q((hY!^^j999 dk! DKA         S__x%?4;%N%N))) A " " " " " " " " "s8}}z'AT['P'P +++"":..K  **IA~~[:- !||;)8#G4'H Z##J3i./0 k012 $+  FffQRRj))QQB ( 3 3 3 38 3 3 3 3 3   r " " " "  ) ) Ebebem1RURU]q4HHSPP55+"45666{++ 5)"23444 k " ""r3cD|jr|Stj|||||S)aCReturn 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. )rvr7r)rXxminyminxmaxymaxs r'rrs,2 }  dD$ = ==r3?c6tj||dkS)a6Return 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)r7orient_polygons)rXsigns r'rrs,  "4TAX > > >>r3)rRF)NrRFrN)r)%__doc__r7shapely.algorithms.polylabelrshapely.errorsrshapely.geometryrrrrr r r shapely.geometry.baser shapely.preparedr __all__rrrrrrrrrrrrrrrrrrr,r3r'rs%33222222,,,,,,/.....!!!!!!   (O3O3O3O3O3O3O3O3d      *   " """"&8888v))) HKHKHKV"+++$(((0@>@>@>@>@>@>@>@>F  {#{#{#{#|>>><??????r3