(&hqddlZddlmZddlZddlmZmZmZddl m Z ddl m Z m Z mZgdZGddeZe d Ze d Ze d Ze d Ze d Ze dZe dZe dZe dZe eddZe dZe dZe dZe dZe dZe dZ e dgd'dZ!e dgd'dZ"e dZ#e dZ$Gd d!e Z%e d(d#Z&e d$Z'e d)d&Z(dS)*N)IntEnum)_geometry_helpers geos_versionlib) ParamEnum)deprecate_positionalmultithreading_enabled requires_geos) GeometryTypeforce_2dforce_3dget_coordinate_dimensionget_dimensionsget_exterior_ring get_geometryget_interior_ringget_mget_num_coordinatesget_num_geometriesget_num_interior_ringsget_num_points get_parts get_point get_precision get_ringsget_srid get_type_idget_xget_yget_z set_precisionset_sridc6eZdZdZdZdZdZdZdZdZ dZ d Z d Z d S) r z'The enumeration of GEOS geometry types.rN) __name__ __module__ __qualname____doc__MISSINGPOINT LINESTRING LINEARRINGPOLYGON MULTIPOINTMULTILINESTRING MULTIPOLYGONGEOMETRYCOLLECTIONL/var/www/html/reinick/venv/lib/python3.11/site-packages/shapely/_geometry.pyr r *sD11G EJJGJOLr:r c &tj|fi|S)aReturn the type ID of a geometry. Possible values are: - None (missing) is -1 - POINT is 0 - LINESTRING is 1 - LINEARRING is 2 - POLYGON is 3 - MULTIPOINT is 4 - MULTILINESTRING is 5 - MULTIPOLYGON is 6 - GEOMETRYCOLLECTION is 7 Parameters ---------- geometry : Geometry or array_like Geometry or geometries to get the type ID of. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. See Also -------- GeometryType Examples -------- >>> import shapely >>> from shapely import LineString, Point >>> shapely.get_type_id(LineString([(0, 0), (1, 1), (2, 2), (3, 3)])) 1 >>> shapely.get_type_id([Point(1, 2), Point(2, 3)]).tolist() [0, 0] )rrgeometrykwargss r;rr;sJ ?8 . .v . ..r:c &tj|fi|S)aReturn the inherent dimensionality of a geometry. The inherent dimension is 0 for points, 1 for linestrings and linearrings, and 2 for polygons. For geometrycollections it is the max of the containing elements. Empty collections and None values return -1. Parameters ---------- geometry : Geometry or array_like Geometry or geometries to get the dimensionality of. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. Examples -------- >>> import shapely >>> from shapely import GeometryCollection, Point, Polygon >>> point = Point(0, 0) >>> shapely.get_dimensions(point) 0 >>> polygon = Polygon([(0, 0), (0, 10), (10, 10), (10, 0), (0, 0)]) >>> shapely.get_dimensions(polygon) 2 >>> shapely.get_dimensions(GeometryCollection([point, polygon])) 2 >>> shapely.get_dimensions(GeometryCollection([])) -1 >>> shapely.get_dimensions(None) -1 )rrr=s r;rrcsB  h 1 1& 1 11r:c &tj|fi|S)aReturn the dimensionality of the coordinates in a geometry (2, 3 or 4). The return value can be one of the following: * Return 2 for geometries with XY coordinate types, * Return 3 for XYZ or XYM coordinate types (distinguished by :meth:`has_z` or :meth:`has_m`), * Return 4 for XYZM coordinate types, * Return -1 for missing geometries (``None`` values). Note that with GEOS < 3.12, if the first Z coordinate equals ``nan``, this function will return ``2``. Geometries with M coordinates are supported with GEOS >= 3.12. Parameters ---------- geometry : Geometry or array_like Geometry or geometries to get the coordinate dimension of. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. Examples -------- >>> import shapely >>> from shapely import Point >>> shapely.get_coordinate_dimension(Point(0, 0)) 2 >>> shapely.get_coordinate_dimension(Point(0, 0, 1)) 3 >>> shapely.get_coordinate_dimension(None) -1 )rrr=s r;rrsD  ' ; ;F ; ;;r:c &tj|fi|S)a Return the total number of coordinates in a geometry. Returns 0 for not-a-geometry values. Parameters ---------- geometry : Geometry or array_like Geometry or geometries to get the number of coordinates of. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. Examples -------- >>> import shapely >>> from shapely import GeometryCollection, LineString, Point >>> point = Point(0, 0) >>> shapely.get_num_coordinates(point) 1 >>> shapely.get_num_coordinates(Point(0, 0, 0)) 1 >>> line = LineString([(0, 0), (1, 1)]) >>> shapely.get_num_coordinates(line) 2 >>> shapely.get_num_coordinates(GeometryCollection([point, line])) 3 >>> shapely.get_num_coordinates(None) 0 )rrr=s r;rrs>  "8 6 6v 6 66r:c &tj|fi|S)aBReturn the SRID of a geometry. Returns -1 for not-a-geometry values. Parameters ---------- geometry : Geometry or array_like Geometry or geometries to get the SRID of. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. See Also -------- set_srid Examples -------- >>> import shapely >>> from shapely import Point >>> point = Point(0, 0) >>> shapely.get_srid(point) 0 >>> with_srid = shapely.set_srid(point, 4326) >>> shapely.get_srid(with_srid) 4326 )rrr=s r;rrs: < + +F + ++r:c Ltj|tj|fi|S)aUReturn a geometry with its SRID set. Parameters ---------- geometry : Geometry or array_like Geometry or geometries to set the SRID of. srid : int The SRID to set on the geometry. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. See Also -------- get_srid Examples -------- >>> import shapely >>> from shapely import Point >>> point = Point(0, 0) >>> shapely.get_srid(point) 0 >>> with_srid = shapely.set_srid(point, 4326) >>> shapely.get_srid(with_srid) 4326 )rr"npintc)r>sridr?s r;r"r"s&: <"'$-- : :6 : ::r:c &tj|fi|S)aReturn the x-coordinate of a point. Parameters ---------- point : Geometry or array_like Non-point geometries will result in NaN being returned. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. See Also -------- get_y, get_z, get_m Examples -------- >>> import shapely >>> from shapely import MultiPoint, Point >>> shapely.get_x(Point(1, 2)) 1.0 >>> shapely.get_x(MultiPoint([(1, 1), (1, 2)])) nan )rrpointr?s r;rr2 9U % %f % %%r:c &tj|fi|S)aReturn the y-coordinate of a point. Parameters ---------- point : Geometry or array_like Non-point geometries will result in NaN being returned. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. See Also -------- get_x, get_z, get_m Examples -------- >>> import shapely >>> from shapely import MultiPoint, Point >>> shapely.get_y(Point(1, 2)) 2.0 >>> shapely.get_y(MultiPoint([(1, 1), (1, 2)])) nan )rrrIs r;rr-rKr:c &tj|fi|S)aaReturn the z-coordinate of a point. Parameters ---------- point : Geometry or array_like Non-point geometries or geometries without Z dimension will result in NaN being returned. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. See Also -------- get_x, get_y, get_m Examples -------- >>> import shapely >>> from shapely import MultiPoint, Point >>> shapely.get_z(Point(1, 2, 3)) 3.0 >>> shapely.get_z(Point(1, 2)) nan >>> shapely.get_z(MultiPoint([(1, 1, 1), (2, 2, 2)])) nan )rr rIs r;r r Is8 9U % %f % %%r:z3.12.0c &tj|fi|S)aReturn the m-coordinate of a point. .. versionadded:: 2.1.0 Parameters ---------- point : Geometry or array_like Non-point geometries or geometries without M dimension will result in NaN being returned. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. See Also -------- get_x, get_y, get_z Examples -------- >>> import shapely >>> from shapely import Point, from_wkt >>> shapely.get_m(from_wkt("POINT ZM (1 2 3 4)")) 4.0 >>> shapely.get_m(from_wkt("POINT M (1 2 4)")) 4.0 >>> shapely.get_m(Point(1, 2, 3)) nan >>> shapely.get_m(from_wkt("MULTIPOINT M ((1 1 1), (2 2 2))")) nan )rrrIs r;rrhsB 9U % %f % %%r:c Ltj|tj|fi|S)atReturn the nth point of a linestring or linearring. Parameters ---------- geometry : Geometry or array_like Geometry or geometries to get the point of. index : int or array_like Negative values count from the end of the linestring backwards. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. See Also -------- get_num_points Examples -------- >>> import shapely >>> from shapely import LinearRing, LineString, MultiPoint, Point >>> line = LineString([(0, 0), (1, 1), (2, 2), (3, 3)]) >>> shapely.get_point(line, 1) >>> shapely.get_point(line, -2) >>> shapely.get_point(line, [0, 3]).tolist() [, ] The function works the same for LinearRing input: >>> shapely.get_point(LinearRing([(0, 0), (1, 1), (2, 2), (0, 0)]), 1) For non-linear geometries it returns None: >>> shapely.get_point(MultiPoint([(0, 0), (1, 1), (2, 2), (3, 3)]), 1) is None True >>> shapely.get_point(Point(1, 1), 0) is None True )rrrErFr>indexr?s r;rrs'T =275>> < ` for other keyword arguments. See Also -------- get_point get_num_geometries Examples -------- >>> import shapely >>> from shapely import LineString, MultiPoint >>> shapely.get_num_points(LineString([(0, 0), (1, 1), (2, 2), (3, 3)])) 4 >>> shapely.get_num_points(MultiPoint([(0, 0), (1, 1), (2, 2), (3, 3)])) 0 >>> shapely.get_num_points(None) 0 )rrr=s r;rrs>  h 1 1& 1 11r:c &tj|fi|S)aaReturn the exterior ring of a polygon. Parameters ---------- geometry : Geometry or array_like Geometry or geometries to get the exterior ring of. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. See Also -------- get_interior_ring Examples -------- >>> import shapely >>> from shapely import Point, Polygon >>> shapely.get_exterior_ring(Polygon([(0, 0), (0, 10), (10, 10), (10, 0), (0, 0)])) >>> shapely.get_exterior_ring(Point(1, 1)) is None True )rrr=s r;rrs2   4 4V 4 44r:c Ltj|tj|fi|S)a}Return the nth interior ring of a polygon. The number of interior rings in non-polygons equals zero. Parameters ---------- geometry : Geometry or array_like Geometry or geometries to get the interior ring of. index : int or array_like Negative values count from the end of the interior rings backwards. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. See Also -------- get_exterior_ring get_num_interior_rings Examples -------- >>> import shapely >>> from shapely import Point, Polygon >>> polygon_with_hole = Polygon( ... [(0, 0), (0, 10), (10, 10), (10, 0), (0, 0)], ... holes=[[(2, 2), (2, 4), (4, 4), (4, 2), (2, 2)]] ... ) >>> shapely.get_interior_ring(polygon_with_hole, 0) >>> shapely.get_interior_ring(polygon_with_hole, 1) is None True >>> polygon = Polygon([(0, 0), (0, 10), (10, 10), (10, 0), (0, 0)]) >>> shapely.get_interior_ring(polygon, 0) is None True >>> shapely.get_interior_ring(Point(0, 0), 0) is None True )rrrErFrPs r;rrs(N  275>> D DV D DDr:c &tj|fi|S)aReturn number of internal rings in a polygon. Returns 0 for not-a-geometry values. Parameters ---------- geometry : Geometry or array_like Geometry or geometries to get the number of interior rings of. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. See Also -------- get_exterior_ring get_interior_ring Examples -------- >>> import shapely >>> from shapely import Point, Polygon >>> polygon = Polygon([(0, 0), (0, 10), (10, 10), (10, 0), (0, 0)]) >>> shapely.get_num_interior_rings(polygon) 0 >>> polygon_with_hole = Polygon( ... [(0, 0), (0, 10), (10, 10), (10, 0), (0, 0)], ... holes=[[(2, 2), (2, 4), (4, 4), (4, 2), (2, 2)]] ... ) >>> shapely.get_num_interior_rings(polygon_with_hole) 1 >>> shapely.get_num_interior_rings(Point(0, 0)) 0 >>> shapely.get_num_interior_rings(None) 0 )rrr=s r;rr'sJ  %h 9 9& 9 99r:c Ltj|tj|fi|S)aReturn the nth geometry from a collection of geometries. Parameters ---------- geometry : Geometry or array_like Geometry or geometries to get the nth geometry of. index : int or array_like Negative values count from the end of the collection backwards. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. Notes ----- - simple geometries act as length-1 collections - out-of-range values return None See Also -------- get_num_geometries, get_parts Examples -------- >>> import shapely >>> from shapely import Point, MultiPoint >>> multipoint = MultiPoint([(0, 0), (1, 1), (2, 2), (3, 3)]) >>> shapely.get_geometry(multipoint, 1) >>> shapely.get_geometry(multipoint, -1) >>> shapely.get_geometry(multipoint, 5) is None True >>> shapely.get_geometry(Point(1, 1), 0) >>> shapely.get_geometry(Point(1, 1), 1) is None True )rrrErFrPs r;rrRs(N  Hbgenn ? ? ? ??r: return_indexFctj|tj}tj|}|jdkrt d|rt j|St j|dS)aGet parts of each GeometryCollection or Multi* geometry object. A copy of each geometry in the GeometryCollection or Multi* geometry object is returned. Note: This does not return the individual parts of Multi* geometry objects in a GeometryCollection. You may need to call this function multiple times to return individual parts of Multi* geometry objects in a GeometryCollection. Parameters ---------- geometry : Geometry or array_like Geometry or geometries to get the parts of. return_index : bool, default False If True, will return a tuple of ndarrays of (parts, indexes), where indexes are the indexes of the original geometries in the source array. Notes ----- .. deprecated:: 2.1.0 A deprecation warning is shown if ``return_index`` is specified as a positional argument. This will need to be specified as a keyword argument in a future release. Returns ------- ndarray of parts or tuple of (parts, indexes) See Also -------- get_geometry, get_rings Examples -------- >>> import shapely >>> from shapely import MultiPoint >>> shapely.get_parts(MultiPoint([(0, 1), (2, 3)])).tolist() [, ] >>> parts, index = shapely.get_parts([MultiPoint([(0, 1)]), MultiPoint([(4, 5), (6, 7)])], return_index=True) >>> parts.tolist() [, , ] >>> index.tolist() [0, 1, 1] dtyper%Array should be one dimensionalrrEasarrayobject_ atleast_1dndim ValueErrorrrr>rWs r;rrsudz("*555H}X&&H}:;;;5 *8444  &x 0 0 33r:ctj|tj}tj|}|jdkrt d|rt j|dSt j|ddS)a>Get rings of Polygon geometry object. For each Polygon, the first returned ring is always the exterior ring and potential subsequent rings are interior rings. If the geometry is not a Polygon, nothing is returned (empty array for scalar geometry input or no element in output array for array input). Parameters ---------- geometry : Geometry or array_like Geometry or geometries to get the rings of. return_index : bool, default False If True, will return a tuple of ndarrays of (rings, indexes), where indexes are the indexes of the original geometries in the source array. Notes ----- .. deprecated:: 2.1.0 A deprecation warning is shown if ``return_index`` is specified as a positional argument. This will need to be specified as a keyword argument in a future release. Returns ------- ndarray of rings or tuple of (rings, indexes) See Also -------- get_exterior_ring, get_interior_ring, get_parts Examples -------- >>> import shapely >>> from shapely import Polygon >>> polygon_with_hole = Polygon( ... [(0, 0), (0, 10), (10, 10), (10, 0), (0, 0)], ... holes=[[(2, 2), (2, 4), (4, 4), (4, 2), (2, 2)]] ... ) >>> shapely.get_rings(polygon_with_hole).tolist() [, ] With ``return_index=True``: >>> polygon = Polygon([(0, 0), (2, 0), (2, 2), (0, 2), (0, 0)]) >>> rings, index = shapely.get_rings( ... [polygon, polygon_with_hole], ... return_index=True ... ) >>> rings.tolist() [, , ] >>> index.tolist() [0, 1, 1] rYr%r[T) extract_ringsrr\rbs r;rrszz("*555H}X&&H}:;;;I *84HHHH  &xt D D DQ GGr:c &tj|fi|S)aReturn number of geometries in a collection. Returns 0 for not-a-geometry values. The number of geometries in points, linestrings, linearrings and polygons equals one. Parameters ---------- geometry : Geometry or array_like Geometry or geometries to get the number of geometries of. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. See Also -------- get_num_points get_geometry Examples -------- >>> import shapely >>> from shapely import MultiPoint, Point >>> shapely.get_num_geometries(MultiPoint([(0, 0), (1, 1), (2, 2), (3, 3)])) 4 >>> shapely.get_num_geometries(Point(1, 1)) 1 >>> shapely.get_num_geometries(None) 0 )rrr=s r;rrs>  !( 5 5f 5 55r:c &tj|fi|S)a0Get the precision of a geometry. If a precision has not been previously set, it will be 0 (double precision). Otherwise, it will return the precision grid size that was set on a geometry. Returns NaN for not-a-geometry values. Parameters ---------- geometry : Geometry or array_like Geometry or geometries to get the precision of. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. See Also -------- set_precision Examples -------- >>> import shapely >>> from shapely import Point >>> point = Point(1, 1) >>> shapely.get_precision(point) 0.0 >>> geometry = shapely.set_precision(point, 1.0) >>> shapely.get_precision(geometry) 1.0 >>> shapely.get_precision(None) nan )rrr=s r;rr7sF  X 0 0 0 00r:ceZdZdZdZdZdS)SetPrecisionModerr%r&N)r,r-r. valid_output pointwisekeep_collapsedr9r:r;rhrh]sLINNNr:rhric bt|trt|}n#t j|st d|tjkr'tdkrtj dtdtj ||t j|fi|S)adReturn geometry with the precision set to a precision grid size. By default, geometries use double precision coordinates (grid_size = 0). Coordinates will be rounded if the precision grid specified is less precise than the input geometry. Duplicated vertices will be dropped from lines and polygons for grid sizes greater than 0. Line and polygon geometries may collapse to empty geometries if all vertices are closer together than ``grid_size`` or if a polygon becomes significantly narrower than ``grid_size``. Spikes or sections in polygons narrower than ``grid_size`` after rounding the vertices will be removed, which can lead to multipolygons or empty geometries. Z values, if present, will not be modified. Notes ----- * subsequent operations will always be performed in the precision of the geometry with higher precision (smaller "grid_size"). That same precision will be attached to the operation outputs. * input geometries should be geometrically valid; unexpected results may occur if input geometries are not. * the geometry returned will be in :ref:`mild canonical form `, and the order of vertices can change and should not be relied upon. * returns None if geometry is None. Parameters ---------- geometry : Geometry or array_like Geometry or geometries to set the precision of. grid_size : float Precision grid size. If 0, will use double precision (will not modify geometry if precision grid size was not previously set). If this value is more precise than input geometry, the input geometry will not be modified. mode : {'valid_output', 'pointwise', 'keep_collapsed'}, default 'valid_output' This parameter determines the way a precision reduction is applied on the geometry. There are three modes: 1. `'valid_output'` (default): The output is always valid. Collapsed geometry elements (including both polygons and lines) are removed. Duplicate vertices are removed. 2. `'pointwise'`: Precision reduction is performed pointwise. Output geometry may be invalid due to collapse or self-intersection. Duplicate vertices are not removed. In GEOS this option is called NO_TOPO. .. note:: 'pointwise' mode requires at least GEOS 3.10. It is accepted in earlier versions, but the results may be unexpected. 3. `'keep_collapsed'`: Like the default mode, except that collapsed linear geometry elements are preserved. Collapsed polygonal input elements are removed. Duplicate vertices are removed. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. See Also -------- get_precision Examples -------- >>> import shapely >>> from shapely import LineString, Point >>> shapely.set_precision(Point(0.9, 0.9), 1.0) >>> shapely.set_precision(Point(0.9, 0.9, 0.9), 1.0) >>> shapely.set_precision(LineString([(0, 0), (0, 0.1), (0, 1), (1, 1)]), 1.0) >>> shapely.set_precision(LineString([(0, 0), (0, 0.1), (0.1, 0.1)]), 1.0, mode="valid_output") >>> shapely.set_precision(LineString([(0, 0), (0, 0.1), (0.1, 0.1)]), 1.0, mode="pointwise") >>> shapely.set_precision(LineString([(0, 0), (0, 0.1), (0.1, 0.1)]), 1.0, mode="keep_collapsed") >>> shapely.set_precision(None, 1.0) is None True zmode only accepts scalar values)r' rz+'pointwise' is only supported for GEOS 3.10r&) stacklevel) isinstancestrrh get_valuerEisscalar TypeErrorrjrwarningswarn UserWarningrr!rF)r> grid_sizemoder?s r;r!r!csd$;))$// [  ;9::: )))lZ.G.G 9      Xy"'$-- J J6 J JJr:c &tj|fi|S)aForce the dimensionality of a geometry to 2D. Parameters ---------- geometry : Geometry or array_like Geometry or geometries to force to 2D. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. Examples -------- >>> import shapely >>> from shapely import LineString, Point, Polygon, from_wkt >>> shapely.force_2d(Point(0, 0, 1)) >>> shapely.force_2d(Point(0, 0)) >>> shapely.force_2d(LineString([(0, 0, 0), (0, 1, 1), (1, 1, 2)])) >>> shapely.force_2d(from_wkt("POLYGON Z EMPTY")) >>> shapely.force_2d(None) is None True )rr r=s r;r r s6 < + +F + ++r:c tj|rtdt j||fi|S)aForce the dimensionality of a geometry to 3D. 2D geometries will get the provided Z coordinate; Z coordinates of 3D geometries are unchanged (unless they are nan). Note that for empty geometries, 3D is only supported since GEOS 3.9 and then still only for simple geometries (non-collections). Parameters ---------- geometry : Geometry or array_like Geometry or geometries to force to 3D. z : float or array_like, default 0.0 The Z coordinate value to set on the geometry. **kwargs See :ref:`NumPy ufunc docs ` for other keyword arguments. Examples -------- >>> import shapely >>> from shapely import LineString, Point >>> shapely.force_3d(Point(0, 0), z=3) >>> shapely.force_3d(Point(0, 0, 0), z=3) >>> shapely.force_3d(LineString([(0, 0), (0, 1), (1, 1)])) >>> shapely.force_3d(None) is None True z1It is not allowed to set the Z coordinate to NaN.)rEisnananyrarr )r>zr?s r;r r sIB x{{NLMMM <! . .v . ..r:)F)ri)rz))rtenumrnumpyrEshapelyrrr shapely._enumrshapely.decoratorsrr r __all__r rrrrrr"rrr rrrrrrrrrrrrhr!r r r9r:r;rs8888888888######    8     7   "$/$/$/N 2 2 2F!<!<!;;;D&&&6&&&6&&&<x&&&J)=)=)=X222H5556&E&E&ER$:$:$:T&@&@&@d~&'':4:4:4(':4L~&''EHEHEH('EHP666B"1"1"1Jy [K[K[K[K|,,,:"/"/"/"/"/"/r: