a J5dg@s*dZddlmZddlZddlmZddlmZmZddl m Z m Z m Z m Z mZmZmZddlmZmZddlmZdd lmZgd ZGd d d ZeZejZejZejZejZej Z d&ddZ!d'ddZ"ddZ#ddZ$ddZ%ddZ&ddZ'GdddZ(e(j)Z)d(dd Z*d!d"Z+d)d$d%ZdS)*z-Support for various GEOS geometry operations )warnN) polylabel)GeometryTypeErrorShapelyDeprecationWarning)GeometryCollection LineStringMultiLineString MultiPointPointPolygonshape) BaseGeometryBaseMultipartGeometryorient)prep)cascaded_union linemergeoperator polygonizepolygonize_full transform unary_union triangulatevoronoi_diagramsplitnearest_pointsvalidatesnap shared_paths clip_by_rectr substringc@s>eZdZddZddZddZddd Zd d Zd d ZdS)CollectionOperatorc Cs<t|tr|Sz t|WSttfy6t|YS0dSN) isinstancer r ValueErrorAttributeErrorr)selfobr)G/var/www/html/django/DPS/env/lib/python3.9/site-packages/shapely/ops.pyshapeup,s   zCollectionOperator.shapeupc srt|ddp|}zAz1CollectionOperator.polygonize..)getattriter TypeErrorshapelyrr,)r'linessourceobs collectionr)r1r*r5s  ( zCollectionOperator.polygonizec slt|ddp|}z.)r4r5r6r7r)r'r8r9r:r)r1r*rEs  (z"CollectionOperator.polygonize_fullFcCsd}t|dddkr|}n\t|dr:tdd|jD}n.__iter__cSsg|] }|jqSr)r=r?r)r)r*r2lr3zCannot linemerge )directed)r4hasattrrr,r&r%r7Z line_merge)r'r8rAr9r)r)r*r]s   zCollectionOperator.linemergecCstdtddtj|ddS)zReturns the union of a sequence of geometries .. deprecated:: 1.8 This function was superseded by :meth:`unary_union`. zKThe 'cascaded_union()' function is deprecated. Use 'unary_union()' instead.) stacklevelNZaxis)rrr7 union_allr'r,r)r)r*rss z!CollectionOperator.cascaded_unioncCstj|ddS)zReturns the union of a sequence of geometries Usually used to convert a collection into the smallest set of polygons that cover the same area. NrE)r7rFrGr)r)r*rszCollectionOperator.unary_unionN)F) __name__ __module__ __qualname__r+rrrrrr)r)r)r*r"+s   r"FcCs tj|||d}dd|jDS)aCreates the Delaunay triangulation and returns 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_edgescSsg|]}|qSr)r))r/gr)r)r*r2r3ztriangulate..)r7Zdelaunay_trianglesr,)geomrLedgesr;r)r)r*rsrc Csztj||||d}WnPtjyf}z6d}|d|d7}|rH|d7}t||WYd}~n d}~00|jdkr|t|gS|S)aT Constructs 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) )rLZ extend_torMz;Could not create Voronoi Diagram with the specified inputs (z).z0 Try running again with default tolerance value.Nr)r7Zvoronoi_polygonsZ GEOSExceptionr%r<r)rOenveloperLrPresulterrZerrstrr)r)r*rs*    rcCs t|Sr#)r7Zis_valid_reasonrOr)r)r*rsrcst|jr |S|jdvr(zz|jdvr>t|tt|jWS|jdkrt|jtt|jj}tfdd|jD}t|||WSWnty$|jdvrt|fdd|jDYS|jdkr t|jfdd|jjD}tfd d|jD}t|||YSYn0nH|j d sB|jd kr^t|fd d|j DSt d |jddS)aqApplies `func` to all coordinates of `geom` and 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 rrVr c3s(|] }t|tt|jVqdSr#)typezipr>r/ringfuncr)r* sztransform..csg|] }|qSr)r)r/cr[r)r*r2r3ztransform..csg|] }|qSr)r)r^r[r)r*r2!r3c3s*|]"}t|fdd|jDVqdS)csg|] }|qSr)r)r^r[r)r*r2#r3z'transform...N)rWr>rYr[r)r*r]"sZMultircsg|]}t|qSr))r)r/partr[r)r*r2)r3zType z not recognizedN) is_emptyr<rWrXr>ZexteriorlistZ interiorsr6 startswithr,r)r\rOshellZholesr)r[r*rs0%         rcCsLt||}|dur,|jr$tdntdt|d}t|d}||fS)zReturns the calculated nearest points in the input geometries The points are returned in the same order as the input geometries. Nz!The first input geometry is emptyz"The second input geometry is emptyr)r7Z shortest_linerar%Z get_point)g1g2seqp1p2r)r)r*r.s    rcCst|||S)aSnap one geometry to another with a given tolerance Vertices of the first geometry are snapped to vertices of the second geometry. The resulting snapped geometry is returned. The input geometries are not modified. Parameters ---------- g1 : geometry The first geometry g2 : geometry The second geometry tolerance : float The snapping tolerance Example ------- >>> square = Polygon([(1,1), (2, 1), (2, 2), (1, 2), (1, 1)]) >>> line = LineString([(0,0), (0.8, 0.8), (1.8, 0.95), (2.6, 0.5)]) >>> result = snap(line, square, 0.5) >>> result.wkt 'LINESTRING (0 0, 1 1, 2 1, 2.6 0.5)' )r7r)rfrgrLr)r)r*r?srcCs0t|tstdt|ts$tdt||S)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$rrr7r)rfrgr)r)r*rZs   rc@sHeZdZeddZeddZeddZeddZed d Zd S) SplitOpcsNttstdt|ts$tdj|}tfddt|DS)z!Split a Polygon with a LineStringz First argument must be a Polygonz$Second argument must be a LineStringcsg|]}|r|qSr))containsZrepresentative_point)r/Zpgpolyr)r*r2sz4SplitOp._split_polygon_with_line..)r$r rrboundaryunionrr)rnsplitterrpr)rmr*_split_polygon_with_liners    z SplitOp._split_polygon_with_linecCs|jdvr|j}t|ts"tdt|ts>t|ts>td||}|ddkr^tdn(|ddksv|ddkr||S|gSd S) zCSplit a LineString with another (Multi)LineString or (Multi)Polygon)r MultiPolygon#First argument must be a LineStringz@Second argument must be either a LineString or a MultiLineStringr1z2Input geometry segment overlaps with the splitter.0N) r<ror$rrrZrelater% difference)r0rqZrelationr)r)r*_split_line_with_lines        zSplitOp._split_line_with_linec CsVt|tstdt|ts$td||ds6|gS|jd|jdkrP|gS||}t|j}d}tt |dD]}||}||d}|d|d}|d|d} |d| dd} || 7}||krt|d |dt||dd gS||krxt|d |d|jdgt|jdg||dd gSqx|gS) zSplit a LineString with a PointrtzSecond argument must be a Pointz 0********rrKrerC?N) r$rrr Zrelate_patternr>projectrbrangelen) r0rqZdistance_on_liner>current_positioniZpoint1Zpoint2ZdxZdysegment_lengthr)r)r*_split_line_with_points2       , zSplitOp._split_line_with_pointcCsht|tstdt|ts$td|g}|jD]2}g}tdd|D]}|t||qF|}q0|S)z$Split a LineString with a MultiPointrtz$Second argument must be a MultiPointcSs|j Sr#)ra)xr)r)r*r3z5SplitOp._split_line_with_multipoint..) r$rrr r,filterextendrkr)r0rqchunksptZ new_chunkschunkr)r)r*_split_line_with_multipoints   z#SplitOp._split_line_with_multipointcs|jdvr"tfdd|jDS|jdkrvjdvr>tj}qjdkrPtj}qjdkrbtj}qtdjd nB|jd krjdkrtj}qtd jd ntd |jd t||S)aS Splits a geometry by another geometry and returns 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 Example ------- >>> pt = Point((1, 1)) >>> line = LineString([(0,0), (2,2)]) >>> result = split(line, pt) >>> result.wkt 'GEOMETRYCOLLECTION (LINESTRING (0 0, 1 1), LINESTRING (1 1, 2 2))' )rrscs$g|]}t|jD]}|qqSr))rkrr,)r/r`rrqr)r*r2r3z!SplitOp.split..r)rrr rsr r zSplitting a LineString with a z is not supportedr zSplitting a Polygon with a z Splitting z geometry is not supported) r<rr,rkryrrrrr)rOrqZ split_funcr)rr*rs0          z SplitOp.splitN) rHrIrJ staticmethodrrryrrrr)r)r)r*rkqs   ) rkc Cst|tstd|jd||kr0|||S|sV||jkrV||jkrV||j|S|s~| |jkr~| |jkr~|d|S|r|dkr|dkr|d|S|r| dkr| dkr|d|S|r||j9}||j9}|dkr|krnnt|||jkr||S|dkr*|krLnnt|||jkrL||S||}||}|dkrt|j|}|dkr|j|}||k}|r||}}|dkrd}|r|j|jfg}n|j|jfg}t |j }d} t ||ddD]r\} } || kr|kr"nn | | n| |kr2qf| | d| dd| d| ddd7} q|r| |j|jft |}n| |j|jft|S)a_Return 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.rreNrCrz)r$rrr<Z interpolatelengthabsryrbr>rXappendreversed) rOZ start_distZend_dist normalizedZ start_pointZ end_pointreverseZ vertex_listr>Zcurrent_distancerirjr)r)r*r!.sd;       . 0            4 r!cCs|jr |St|||||S)a]Returns 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 ----- Requires GEOS >= 3.5.0 New in 1.7. )rar7r )rOZxminZyminZxmaxZymaxr)r)r*r sr ?csBt|tr(|ttfdd|jSt|tfr>t|S|S)aA 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. 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 cs t|Sr#rrUsignr)r*rr3zorient..)r$r __class__rbmapr,r orient_)rOrr)rr*rs    r)rKF)NrKF)F)r),__doc__warningsrr7Zshapely.algorithms.polylabelrZshapely.errorsrrZshapely.geometryrrrr r r r Zshapely.geometry.baser rZshapely.geometry.polygonrrZshapely.preparedr__all__r"rrrrrrrrrrrrrrkrr!r r)r)r)r*s:  $   _  :K;