
    \iq                     :   S SK r S SKJr  S SKrS SKJrJrJr  S SK	J
r
  S SKJrJrJr  / SQr " S S\5      r\S	 5       r\S
 5       r\S 5       r\S 5       r\S 5       r\S 5       r\S 5       r\S 5       r\S 5       r\\" S5      S 5       5       r\S 5       r\S 5       r\S 5       r\S 5       r\S 5       r\S 5       r \" S/5      S$S j5       r!\" S/5      S$S j5       r"\S 5       r#\S 5       r$ " S S \
5      r%\S%S! j5       r&\S" 5       r'\S&S# j5       r(g)'    N)IntEnum)_geometry_helpersgeos_versionlib)	ParamEnum)deprecate_positionalmultithreading_enabledrequires_geos)GeometryTypeforce_2dforce_3dget_coordinate_dimensionget_dimensionsget_exterior_ringget_geometryget_interior_ringget_mget_num_coordinatesget_num_geometriesget_num_interior_ringsget_num_points	get_parts	get_pointget_precision	get_ringsget_sridget_type_idget_xget_yget_zset_precisionset_sridc                   <    \ rS rSrSrSrSrSrSrSr	Sr
S	rS
rSrSrg)r   *   z'The enumeration of GEOS geometry types.r                         N)__name__
__module____qualname____firstlineno____doc__MISSINGPOINT
LINESTRING
LINEARRINGPOLYGON
MULTIPOINTMULTILINESTRINGMULTIPOLYGONGEOMETRYCOLLECTION__static_attributes__r-       T/var/www/html/kml_chatgpt/mouzaenv/lib/python3.13/site-packages/shapely/_geometry.pyr   r   *   s3    1GEJJGJOLr=   r   c                 0    [         R                  " U 40 UD6$ )a  Return 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 <ufuncs.kwargs>` 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]

)r   r   geometrykwargss     r>   r   r   ;   s    J ??8.v..r=   c                 0    [         R                  " U 40 UD6$ )an  Return 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 <ufuncs.kwargs>` 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

)r   r   r@   s     r>   r   r   c   s    B h1&11r=   c                 0    [         R                  " U 40 UD6$ )a  Return 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 <ufuncs.kwargs>` 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

)r   r   r@   s     r>   r   r      s    D '';F;;r=   c                 0    [         R                  " U 40 UD6$ )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 <ufuncs.kwargs>` 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

)r   r   r@   s     r>   r   r      s    > ""86v66r=   c                 0    [         R                  " U 40 UD6$ )a  Return 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 <ufuncs.kwargs>` 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

)r   r   r@   s     r>   r   r      s    : <<+F++r=   c                 Z    [         R                  " U [        R                  " U5      40 UD6$ )a  Return 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 <ufuncs.kwargs>` 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

)r   r"   npintc)rA   sridrB   s      r>   r"   r"      s"    : <<"''$-:6::r=   c                 0    [         R                  " U 40 UD6$ )a  Return 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 <ufuncs.kwargs>` 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

)r   r   pointrB   s     r>   r   r         2 99U%f%%r=   c                 0    [         R                  " U 40 UD6$ )a  Return 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 <ufuncs.kwargs>` 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

)r   r   rL   s     r>   r   r   -  rN   r=   c                 0    [         R                  " U 40 UD6$ )a  Return 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 <ufuncs.kwargs>` 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

)r   r    rL   s     r>   r    r    I  s    8 99U%f%%r=   z3.12.0c                 0    [         R                  " U 40 UD6$ )as  Return 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 <ufuncs.kwargs>` 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

)r   r   rL   s     r>   r   r   h  s    B 99U%f%%r=   c                 Z    [         R                  " U [        R                  " U5      40 UD6$ )a  Return 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 <ufuncs.kwargs>` 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)
<POINT (1 1)>
>>> shapely.get_point(line, -2)
<POINT (2 2)>
>>> shapely.get_point(line, [0, 3]).tolist()
[<POINT (0 0)>, <POINT (3 3)>]

The function works the same for LinearRing input:

>>> shapely.get_point(LinearRing([(0, 0), (1, 1), (2, 2), (0, 0)]), 1)
<POINT (1 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

)r   r   rH   rI   rA   indexrB   s      r>   r   r     s#    T ==2775><V<<r=   c                 0    [         R                  " U 40 UD6$ )a  Return the number of points in a linestring or linearring.

Returns 0 for not-a-geometry values. The number of points in geometries
other than linestring or linearring equals zero.

Parameters
----------
geometry : Geometry or array_like
    Geometry or geometries to get the number of points of.
**kwargs
    See :ref:`NumPy ufunc docs <ufuncs.kwargs>` 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

)r   r   r@   s     r>   r   r     s    > h1&11r=   c                 0    [         R                  " U 40 UD6$ )a  Return 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 <ufuncs.kwargs>` 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)]))
<LINEARRING (0 0, 0 10, 10 10, 10 0, 0 0)>
>>> shapely.get_exterior_ring(Point(1, 1)) is None
True

)r   r   r@   s     r>   r   r     s    2   4V44r=   c                 Z    [         R                  " U [        R                  " U5      40 UD6$ )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 <ufuncs.kwargs>` 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)
<LINEARRING (2 2, 2 4, 4 4, 4 2, 2 2)>
>>> 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

)r   r   rH   rI   rS   s      r>   r   r     s%    N   2775>DVDDr=   c                 0    [         R                  " U 40 UD6$ )aH  Return 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 <ufuncs.kwargs>` 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

)r   r   r@   s     r>   r   r   '  s    J %%h9&99r=   c                 Z    [         R                  " U [        R                  " U5      40 UD6$ )a  Return 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 <ufuncs.kwargs>` 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)
<POINT (1 1)>
>>> shapely.get_geometry(multipoint, -1)
<POINT (3 3)>
>>> shapely.get_geometry(multipoint, 5) is None
True
>>> shapely.get_geometry(Point(1, 1), 0)
<POINT (1 1)>
>>> shapely.get_geometry(Point(1, 1), 1) is None
True

)r   r   rH   rI   rS   s      r>   r   r   R  s%    N Hbggen???r=   return_indexc                    [         R                  " U [         R                  S9n [         R                  " U 5      n U R                  S:w  a  [        S5      eU(       a  [        R                  " U 5      $ [        R                  " U 5      S   $ )aP  Get 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()
[<POINT (0 1)>, <POINT (2 3)>]
>>> parts, index = shapely.get_parts([MultiPoint([(0, 1)]), MultiPoint([(4, 5), (6, 7)])], return_index=True)
>>> parts.tolist()
[<POINT (0 1)>, <POINT (4 5)>, <POINT (6 7)>]
>>> index.tolist()
[0, 1, 1]

dtyper&   Array should be one dimensionalr   rH   asarrayobject_
atleast_1dndim
ValueErrorr   r   rA   rZ   s     r>   r   r     si    d zz("**5H}}X&H}}:;; **844&&x033r=   c                    [         R                  " U [         R                  S9n [         R                  " U 5      n U R                  S:w  a  [        S5      eU(       a  [        R                  " U SS9$ [        R                  " U SS9S   $ )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()
[<LINEARRING (0 0, 0 10, 10 10, 10 0, 0 0)>,
 <LINEARRING (2 2, 2 4, 4 4, 4 2, 2 2)>]

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()
[<LINEARRING (0 0, 2 0, 2 2, 0 2, 0 0)>,
 <LINEARRING (0 0, 0 10, 10 10, 10 0, 0 0)>,
 <LINEARRING (2 2, 2 4, 4 4, 4 2, 2 2)>]
>>> index.tolist()
[0, 1, 1]

r\   r&   r^   T)extract_ringsr   r_   re   s     r>   r   r     sm    z zz("**5H}}X&H}}:;; **84HH&&xtDQGGr=   c                 0    [         R                  " U 40 UD6$ )a  Return 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 <ufuncs.kwargs>` 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

)r   r   r@   s     r>   r   r     s    > !!(5f55r=   c                 0    [         R                  " U 40 UD6$ )a  Get 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 <ufuncs.kwargs>` 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

)r   r   r@   s     r>   r   r   7  s    F X000r=   c                        \ rS rSrSrSrSrSrg)SetPrecisionModei]  r   r&   r'   r-   N)r.   r/   r0   r1   valid_output	pointwisekeep_collapsedr<   r-   r=   r>   rk   rk   ]  s    LINr=   rk   c                 l   [        U[        5      (       a  [        R                  U5      nO&[        R
                  " U5      (       d  [        S5      eU[        R                  :X  a$  [        S:  a  [        R                  " S[        SS9  [        R                  " X[        R                  " U5      40 UD6$ )aP  Return 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 <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 <ufuncs.kwargs>` 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)
<POINT (1 1)>
>>> shapely.set_precision(Point(0.9, 0.9, 0.9), 1.0)
<POINT Z (1 1 0.9)>
>>> shapely.set_precision(LineString([(0, 0), (0, 0.1), (0, 1), (1, 1)]), 1.0)
<LINESTRING (0 0, 0 1, 1 1)>
>>> shapely.set_precision(LineString([(0, 0), (0, 0.1), (0.1, 0.1)]), 1.0, mode="valid_output")
<LINESTRING EMPTY>
>>> shapely.set_precision(LineString([(0, 0), (0, 0.1), (0.1, 0.1)]), 1.0, mode="pointwise")
<LINESTRING (0 0, 0 0, 0 0)>
>>> shapely.set_precision(LineString([(0, 0), (0, 0.1), (0.1, 0.1)]), 1.0, mode="keep_collapsed")
<LINESTRING (0 0, 0 0)>
>>> shapely.set_precision(None, 1.0) is None
True

zmode only accepts scalar values)r(   
   r   z+'pointwise' is only supported for GEOS 3.10r'   )
stacklevel)
isinstancestrrk   	get_valuerH   isscalar	TypeErrorrm   r   warningswarnUserWarningr   r!   rI   )rA   	grid_sizemoderB   s       r>   r!   r!   c  s    d $))$/[[9::)))lZ.G9	

 X"''$-J6JJr=   c                 0    [         R                  " U 40 UD6$ )az  Force 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 <ufuncs.kwargs>` for other keyword arguments.

Examples
--------
>>> import shapely
>>> from shapely import LineString, Point, Polygon, from_wkt
>>> shapely.force_2d(Point(0, 0, 1))
<POINT (0 0)>
>>> shapely.force_2d(Point(0, 0))
<POINT (0 0)>
>>> shapely.force_2d(LineString([(0, 0, 0), (0, 1, 1), (1, 1, 2)]))
<LINESTRING (0 0, 0 1, 1 1)>
>>> shapely.force_2d(from_wkt("POLYGON Z EMPTY"))
<POLYGON EMPTY>
>>> shapely.force_2d(None) is None
True

)r   r   r@   s     r>   r   r     s    6 <<+F++r=   c                     [         R                  " U5      R                  5       (       a  [        S5      e[        R
                  " X40 UD6$ )a  Force 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 <ufuncs.kwargs>` for other keyword arguments.

Examples
--------
>>> import shapely
>>> from shapely import LineString, Point
>>> shapely.force_3d(Point(0, 0), z=3)
<POINT Z (0 0 3)>
>>> shapely.force_3d(Point(0, 0, 0), z=3)
<POINT Z (0 0 0)>
>>> shapely.force_3d(LineString([(0, 0), (0, 1), (1, 1)]))
<LINESTRING Z (0 0 0, 0 1 0, 1 1 0)>
>>> shapely.force_3d(None) is None
True

z1It is not allowed to set the Z coordinate to NaN.)rH   isnananyrd   r   r   )rA   zrB   s      r>   r   r     s:    B 
xx{LMM<<.v..r=   )F)rl   )g        ))rw   enumr   numpyrH   shapelyr   r   r   shapely._enumr   shapely.decoratorsr   r	   r
   __all__r   r   r   r   r   r   r"   r   r   r    r   r   r   r   r   r   r   r   r   r   r   rk   r!   r   r   r-   r=   r>   <module>r      sE      8 8 # 87 " $/ $/N  2  2F !< !<H 7 7B , ,> ; ;D & &6 & &6 & &< x&  &J )= )=X 2 2H 5 56 &E &ER $: $:T &@ &@d ~&':4 (:4L ~&'EH (EHP 6 6B "1 "1Jy  [K [K| , ,: "/ "/r=   