
    \i-k                        S r SSKrSSKrSSKJrJrJr  SSKJ	r	J
r
Jr  / SQr\	" S/\S9\
SS j5       5       r\	" S/\S9\
SS	 j5       5       r\	" S
/\S9\
SS j5       5       r\	" S/\S9\
SS j5       5       r\	" S
/\S9\
SS j5       5       r\	" S/\S9\
SS j5       5       r\	" SS
/\S9\
SS j5       5       r\r\
S 5       r\	" S
/\S9\
SS j5       5       r\" S5      \
S 5       5       r\" S5      \
SS.S j5       5       rg)z-Set-theoretic operations on geometry objects.    N)GeometryGeometryTypelib)deprecate_positionalmultithreading_enabledrequires_geos)coverage_unioncoverage_union_all
differencedisjoint_subset_uniondisjoint_subset_union_allintersectionintersection_allsymmetric_differencesymmetric_difference_allunary_unionunion	union_all	grid_size)categoryc                     Ub>  [         R                  " U5      (       d  [        S5      e[        R                  " XU40 UD6$ [        R
                  " X40 UD6$ )a  Return the part of geometry A that does not intersect with geometry B.

If grid_size is nonzero, input coordinates will be snapped to a precision
grid of that size and resulting coordinates will be snapped to that same
grid.  If 0, this operation will use double precision coordinates.  If None,
the highest precision of the inputs will be used, which may be previously
set using set_precision.  Note: returned geometry does not have precision
set unless specified previously by set_precision.

Parameters
----------
a : Geometry or array_like
    Geometry or geometries to subtract b from.
b : Geometry or array_like
    Geometry or geometries to subtract from a.
grid_size : float, optional
    Precision grid size; will use the highest precision of the inputs by default.
**kwargs
    See :ref:`NumPy ufunc docs <ufuncs.kwargs>` for other keyword arguments.

Notes
-----

.. deprecated:: 2.1.0
    A deprecation warning is shown if ``grid_size`` is specified as a
    positional argument. This will need to be specified as a keyword
    argument in a future release.

See Also
--------
set_precision

Examples
--------
>>> import shapely
>>> from shapely import LineString
>>> line = LineString([(0, 0), (2, 2)])
>>> shapely.difference(line, LineString([(1, 1), (3, 3)]))
<LINESTRING (0 0, 1 1)>
>>> shapely.difference(line, LineString())
<LINESTRING (0 0, 2 2)>
>>> shapely.difference(line, None) is None
True
>>> box1 = shapely.box(0, 0, 2, 2)
>>> box2 = shapely.box(1, 1, 3, 3)
>>> shapely.difference(box1, box2).normalize()
<POLYGON ((0 0, 0 2, 1 2, 1 1, 2 1, 2 0, 0 0))>
>>> box1 = shapely.box(0.1, 0.2, 2.1, 2.1)
>>> shapely.difference(box1, box2, grid_size=1)
<POLYGON ((2 0, 0 0, 0 2, 1 2, 1 1, 2 1, 2 0))>

.grid_size parameter only accepts scalar values)npisscalar
ValueErrorr   difference_precr   abr   kwargss       Y/var/www/html/kml_chatgpt/mouzaenv/lib/python3.13/site-packages/shapely/set_operations.pyr   r   '   sR    n {{9%%MNN""1=f==>>!)&))    c                     Ub>  [         R                  " U5      (       d  [        S5      e[        R                  " XU40 UD6$ [        R
                  " X40 UD6$ )a  Return the geometry that is shared between input geometries.

If grid_size is nonzero, input coordinates will be snapped to a precision
grid of that size and resulting coordinates will be snapped to that same
grid.  If 0, this operation will use double precision coordinates.  If None,
the highest precision of the inputs will be used, which may be previously
set using set_precision.  Note: returned geometry does not have precision
set unless specified previously by set_precision.

Parameters
----------
a, b : Geometry or array_like
    Geometry or geometries to intersect with.
grid_size : float, optional
    Precision grid size; will use the highest precision of the inputs by default.
**kwargs
    See :ref:`NumPy ufunc docs <ufuncs.kwargs>` for other keyword arguments.

Notes
-----

.. deprecated:: 2.1.0
    A deprecation warning is shown if ``grid_size`` is specified as a
    positional argument. This will need to be specified as a keyword
    argument in a future release.

See Also
--------
intersection_all
set_precision

Examples
--------
>>> import shapely
>>> from shapely import LineString
>>> line = LineString([(0, 0), (2, 2)])
>>> shapely.intersection(line, LineString([(1, 1), (3, 3)]))
<LINESTRING (1 1, 2 2)>
>>> box1 = shapely.box(0, 0, 2, 2)
>>> box2 = shapely.box(1, 1, 3, 3)
>>> shapely.intersection(box1, box2).normalize()
<POLYGON ((1 1, 1 2, 2 2, 2 1, 1 1))>
>>> box1 = shapely.box(0.1, 0.2, 2.1, 2.1)
>>> shapely.intersection(box1, box2, grid_size=1)
<POLYGON ((2 2, 2 1, 1 1, 1 2, 2 2))>

r   )r   r   r   r   intersection_precr   r   s       r!   r   r   p   sT    d {{9%%MNN$$Q9???A+F++r"   axisc                     [         R                  " U 5      n Uc  U R                  5       n O[         R                  " XU R                  S9n [
        R                  " U 40 UD6$ )a  Return the intersection of multiple geometries.

This function ignores None values when other Geometry elements are present.
If all elements of the given axis are None, an empty GeometryCollection is
returned.

Parameters
----------
geometries : array_like
    Geometries to calculate the intersection of.
axis : int, optional
    Axis along which the operation is performed. The default (None)
    performs the operation over all axes, returning a scalar value.
    Axis may be negative, in which case it counts from the last to the
    first axis.
**kwargs
    See :ref:`NumPy ufunc docs <ufuncs.kwargs>` for other keyword arguments.

Notes
-----

.. deprecated:: 2.1.0
    A deprecation warning is shown if ``axis`` is specified as a
    positional argument. This will need to be specified as a keyword
    argument in a future release.

See Also
--------
intersection

Examples
--------
>>> import shapely
>>> from shapely import LineString
>>> line1 = LineString([(0, 0), (2, 2)])
>>> line2 = LineString([(1, 1), (3, 3)])
>>> shapely.intersection_all([line1, line2])
<LINESTRING (1 1, 2 2)>
>>> shapely.intersection_all([[line1, line2, None]], axis=1).tolist()
[<LINESTRING (1 1, 2 2)>]
>>> shapely.intersection_all([line1, None])
<LINESTRING (0 0, 2 2)>

r%   start)r   asarrayravelrollaxisndimr   r   
geometriesr%   r    s      r!   r   r      sP    ^ J'J|%%'
[[jooN

5f55r"   c                     Ub>  [         R                  " U5      (       d  [        S5      e[        R                  " XU40 UD6$ [        R
                  " X40 UD6$ )a  Return the geometry with the portions of input geometries that do not intersect.

If grid_size is nonzero, input coordinates will be snapped to a precision
grid of that size and resulting coordinates will be snapped to that same
grid.  If 0, this operation will use double precision coordinates.  If None,
the highest precision of the inputs will be used, which may be previously
set using set_precision.  Note: returned geometry does not have precision
set unless specified previously by set_precision.

Parameters
----------
a, b : Geometry or array_like
    Geometry or geometries to evaluate symmetric difference with.
grid_size : float, optional
    Precision grid size; will use the highest precision of the inputs by default.
**kwargs
    See :ref:`NumPy ufunc docs <ufuncs.kwargs>` for other keyword arguments.

Notes
-----

.. deprecated:: 2.1.0
    A deprecation warning is shown if ``grid_size`` is specified as a
    positional argument. This will need to be specified as a keyword
    argument in a future release.

See Also
--------
symmetric_difference_all
set_precision

Examples
--------
>>> import shapely
>>> from shapely import LineString
>>> line = LineString([(0, 0), (2, 2)])
>>> shapely.symmetric_difference(line, LineString([(1, 1), (3, 3)]))
<MULTILINESTRING ((0 0, 1 1), (2 2, 3 3))>
>>> box1 = shapely.box(0, 0, 2, 2)
>>> box2 = shapely.box(1, 1, 3, 3)
>>> shapely.symmetric_difference(box1, box2).normalize()
<MULTIPOLYGON (((1 2, 1 3, 3 3, 3 1, 2 1, 2 2, 1 2)), ((0 0, 0 2, 1 2, 1 1, ...>
>>> box1 = shapely.box(0.1, 0.2, 2.1, 2.1)
>>> shapely.symmetric_difference(box1, box2, grid_size=1)
<MULTIPOLYGON (((2 0, 0 0, 0 2, 1 2, 1 1, 2 1, 2 0)), ((2 2, 1 2, 1 3, 3 3, ...>

r   )r   r   r   r   symmetric_difference_precr   r   s       r!   r   r      sT    d {{9%%MNN,,Q9GGG##A3F33r"   c                     [         R                  " S[        SS9  [        R                  " U 5      n Uc  U R                  5       n O[        R                  " XU R                  S9n [        R                  " U 40 UD6$ )aM  Return the symmetric difference of multiple geometries.

This function ignores None values when other Geometry elements are present.
If all elements of the given axis are None an empty GeometryCollection is
returned.

.. deprecated:: 2.1.0

    This function behaves incorrectly and will be removed in a future
    version. See https://github.com/shapely/shapely/issues/2027 for more
    details.

Parameters
----------
geometries : array_like
    Geometries to calculate the combined symmetric difference of.
axis : int, optional
    Axis along which the operation is performed. The default (None)
    performs the operation over all axes, returning a scalar value.
    Axis may be negative, in which case it counts from the last to the
    first axis.
**kwargs
    See :ref:`NumPy ufunc docs <ufuncs.kwargs>` for other keyword arguments.

Notes
-----

.. deprecated:: 2.1.0
    A deprecation warning is shown if ``axis`` is specified as a
    positional argument. This will need to be specified as a keyword
    argument in a future release.

See Also
--------
symmetric_difference

Examples
--------
>>> import shapely
>>> from shapely import LineString
>>> line1 = LineString([(0, 0), (2, 2)])
>>> line2 = LineString([(1, 1), (3, 3)])
>>> shapely.symmetric_difference_all([line1, line2])
<MULTILINESTRING ((0 0, 1 1), (2 2, 3 3))>
>>> shapely.symmetric_difference_all([[line1, line2, None]], axis=1).tolist()
[<MULTILINESTRING ((0 0, 1 1), (2 2, 3 3))>]
>>> shapely.symmetric_difference_all([line1, None])
<LINESTRING (0 0, 2 2)>
>>> shapely.symmetric_difference_all([None, None])
<GEOMETRYCOLLECTION EMPTY>

zThe symmetric_difference_all function behaves incorrectly and will be removed in a future version. See https://github.com/shapely/shapely/issues/2027 for more details.   )
stacklevelr'   )
warningswarnDeprecationWarningr   r)   r*   r+   r,   r   r   r-   s      r!   r   r   9  sl    n MM	O 	 J'J|%%'
[[jooN
''
=f==r"   c                     Ub>  [         R                  " U5      (       d  [        S5      e[        R                  " XU40 UD6$ [        R
                  " X40 UD6$ )a  Merge geometries into one.

If grid_size is nonzero, input coordinates will be snapped to a precision
grid of that size and resulting coordinates will be snapped to that same
grid.  If 0, this operation will use double precision coordinates.  If None,
the highest precision of the inputs will be used, which may be previously
set using set_precision.  Note: returned geometry does not have precision
set unless specified previously by set_precision.

Parameters
----------
a, b : Geometry or array_like
    Geometry or geometries to merge (union).
grid_size : float, optional
    Precision grid size; will use the highest precision of the inputs by default.
**kwargs
    See :ref:`NumPy ufunc docs <ufuncs.kwargs>` for other keyword arguments.

Notes
-----

.. deprecated:: 2.1.0
    A deprecation warning is shown if ``grid_size`` is specified as a
    positional argument. This will need to be specified as a keyword
    argument in a future release.

See Also
--------
union_all
set_precision

Examples
--------
>>> import shapely
>>> from shapely import LineString
>>> line = LineString([(0, 0), (2, 2)])
>>> shapely.union(line, LineString([(2, 2), (3, 3)]))
<MULTILINESTRING ((0 0, 2 2), (2 2, 3 3))>
>>> shapely.union(line, None) is None
True
>>> box1 = shapely.box(0, 0, 2, 2)
>>> box2 = shapely.box(1, 1, 3, 3)
>>> shapely.union(box1, box2).normalize()
<POLYGON ((0 0, 0 2, 1 2, 1 3, 3 3, 3 1, 2 1, 2 0, 0 0))>
>>> box1 = shapely.box(0.1, 0.2, 2.1, 2.1)
>>> shapely.union(box1, box2, grid_size=1)
<POLYGON ((2 0, 0 0, 0 2, 1 2, 1 3, 3 3, 3 1, 2 1, 2 0))>

r   )r   r   r   r   
union_precr   r   s       r!   r   r     sP    h {{9%%MNN~~aI88899Q$V$$r"   c                    [         R                  " U 5      n Uc  U R                  5       n O[         R                  " XU R                  S9n [
        R                  " U [         R                  " [        R                  5      5      nUb=  [         R                  " U5      (       d  [        S5      e[
        R                  " XA40 UD6$ [
        R                  " U40 UD6$ )a  Return the union of multiple geometries.

This function ignores None values when other Geometry elements are present.
If all elements of the given axis are None an empty GeometryCollection is
returned.

If grid_size is nonzero, input coordinates will be snapped to a precision
grid of that size and resulting coordinates will be snapped to that same
grid.  If 0, this operation will use double precision coordinates.  If None,
the highest precision of the inputs will be used, which may be previously
set using set_precision.  Note: returned geometry does not have precision
set unless specified previously by set_precision.

`unary_union` is an alias of `union_all`.

Parameters
----------
geometries : array_like
    Geometries to merge/union.
grid_size : float, optional
    Precision grid size; will use the highest precision of the inputs by default.
axis : int, optional
    Axis along which the operation is performed. The default (None)
    performs the operation over all axes, returning a scalar value.
    Axis may be negative, in which case it counts from the last to the
    first axis.
**kwargs
    See :ref:`NumPy ufunc docs <ufuncs.kwargs>` for other keyword arguments.

Notes
-----

.. deprecated:: 2.1.0
    A deprecation warning is shown if ``grid_size`` or ``axis`` are
    specified as positional arguments. In a future release, these will
    need to be specified as keyword arguments.

See Also
--------
union
set_precision

Examples
--------
>>> import shapely
>>> from shapely import LineString, Point
>>> line1 = LineString([(0, 0), (2, 2)])
>>> line2 = LineString([(2, 2), (3, 3)])
>>> shapely.union_all([line1, line2])
<MULTILINESTRING ((0 0, 2 2), (2 2, 3 3))>
>>> shapely.union_all([[line1, line2, None]], axis=1).tolist()
[<MULTILINESTRING ((0 0, 2 2), (2 2, 3 3))>]
>>> box1 = shapely.box(0, 0, 2, 2)
>>> box2 = shapely.box(1, 1, 3, 3)
>>> shapely.union_all([box1, box2]).normalize()
<POLYGON ((0 0, 0 2, 1 2, 1 3, 3 3, 3 1, 2 1, 2 0, 0 0))>
>>> box1 = shapely.box(0.1, 0.2, 2.1, 2.1)
>>> shapely.union_all([box1, box2], grid_size=1)
<POLYGON ((2 0, 0 0, 0 2, 1 2, 1 3, 3 3, 3 1, 2 1, 2 0))>
>>> shapely.union_all([None, Point(0, 1)])
<POINT (0 1)>
>>> shapely.union_all([None, None])
<GEOMETRYCOLLECTION EMPTY>
>>> shapely.union_all([])
<GEOMETRYCOLLECTION EMPTY>

r'   r   )r   r)   r*   r+   r,   r   create_collectionintcr   GEOMETRYCOLLECTIONr   r   unary_union_precr   )r.   r   r%   r    collectionss        r!   r   r     s    R J'J|%%'
[[jooN
 ''BGGL;;<K {{9%%MNN##KEfEE??;1&11r"   c                     [        X/40 UD6$ )a  Merge multiple polygons into one.

This is an optimized version of union which assumes the polygons to be
non-overlapping.
If this assumption is not met, the exact result is not guaranteed
(depending on the GEOS version, it may return the input unchanged or raise
an error).

Parameters
----------
a, b : Geometry or array_like
    Geometry or geometries to merge (union).
**kwargs
    See :ref:`NumPy ufunc docs <ufuncs.kwargs>` for other keyword arguments.

See Also
--------
coverage_union_all

Examples
--------
>>> import shapely
>>> from shapely import Polygon
>>> polygon_1 = Polygon([(0, 0), (0, 1), (1, 1), (1, 0), (0, 0)])
>>> polygon_2 = Polygon([(1, 0), (1, 1), (2, 1), (2, 0), (1, 0)])
>>> shapely.coverage_union(polygon_1, polygon_2).normalize()
<POLYGON ((0 0, 0 1, 1 1, 2 1, 2 0, 1 0, 0 0))>

Union with None returns same polygon

>>> shapely.coverage_union(polygon_1, None).normalize()
<POLYGON ((0 0, 0 1, 1 1, 1 0, 0 0))>

)r
   r   r   r    s      r!   r	   r	   /  s    H qf///r"   c                 \   [         R                  " U 5      n Uc  U R                  5       n O3[         R                  " [         R                  " U 5      XR                  S9n [
        R                  " U [         R                  " [        R                  5      5      n[
        R                  " U40 UD6$ )a  Return the union of multiple polygons of a geometry collection.

This is an optimized version of union which assumes the polygons
to be non-overlapping.

This function ignores None values when other Geometry elements are present.
If all elements of the given axis are None, an empty GeometryCollection is
returned (before GEOS 3.12 this was an empty MultiPolygon).

Parameters
----------
geometries : array_like
    Geometries to merge/union.
axis : int, optional
    Axis along which the operation is performed. The default (None)
    performs the operation over all axes, returning a scalar value.
    Axis may be negative, in which case it counts from the last to the
    first axis.
**kwargs
    See :ref:`NumPy ufunc docs <ufuncs.kwargs>` for other keyword arguments.

Notes
-----

.. deprecated:: 2.1.0
    A deprecation warning is shown if ``axis`` is specified as a
    positional argument. This will need to be specified as a keyword
    argument in a future release.

See Also
--------
coverage_union

Examples
--------
>>> import shapely
>>> from shapely import Polygon
>>> polygon_1 = Polygon([(0, 0), (0, 1), (1, 1), (1, 0), (0, 0)])
>>> polygon_2 = Polygon([(1, 0), (1, 1), (2, 1), (2, 0), (1, 0)])
>>> shapely.coverage_union_all([polygon_1, polygon_2]).normalize()
<POLYGON ((0 0, 0 1, 1 1, 2 1, 2 0, 1 0, 0 0))>
>>> shapely.coverage_union_all([polygon_1, None]).normalize()
<POLYGON ((0 0, 0 1, 1 1, 1 0, 0 0))>
>>> shapely.coverage_union_all([None, None]).normalize()
<GEOMETRYCOLLECTION EMPTY>

r'   )r   r)   r*   r+   r,   r   r:   r;   r   r<   r	   r.   r%   r    r>   s       r!   r
   r
   _  s    h J'J|%%'
[[JJz"__

 ''BGGL;;<K k4V44r"   z3.12.0c                    [        U [        5      (       d  U c  [        U[        5      (       d  Uc  O[        U [        5      (       d  U c  [        R                  " X5      n OR[        U[        5      (       d  Uc  [        R                  " X5      nO#[	        U 5      [	        U5      :w  a  [        S5      e[        X/4SS0UD6$ )ak  Merge multiple polygons into one using algorithm optimised for subsets.

This is an optimized version of union which assumes inputs can be
divided into subsets that do not intersect.

If there is only one such subset, performance can be expected to be worse than
:func:`union`. As such, it is recommeded to use ``disjoint_subset_union`` with
GeometryCollections rather than individual geometries.

.. versionadded:: 2.1.0

Parameters
----------
a, b : Geometry or array_like
    Geometry or geometries to merge (union).
**kwargs
    See :ref:`NumPy ufunc docs <ufuncs.kwargs>` for other keyword arguments.

See Also
--------
union
coverage_union
disjoint_subset_union_all

Examples
--------
>>> import shapely
>>> from shapely import Polygon
>>> polygon_1 = Polygon([(0, 0), (0, 1), (1, 1), (1, 0), (0, 0)])
>>> polygon_2 = Polygon([(1, 0), (1, 1), (2, 1), (2, 0), (1, 0)])
>>> shapely.disjoint_subset_union(polygon_1, polygon_2).normalize()
<POLYGON ((0 0, 0 1, 1 1, 2 1, 2 0, 1 0, 0 0))>

Union with None returns same polygon:

>>> shapely.disjoint_subset_union(polygon_1, None).normalize()
<POLYGON ((0 0, 0 1, 1 1, 1 0, 0 0))>
z(Arrays a and b must have the same lengthr%   r   )
isinstancer   r   	full_likelenr   r   r@   s      r!   r   r     s    R 	1h191h19	Ax	 	 AILL	Ax	 	 AILL	Q3q6	CDD$aV>!>v>>r"   )r%   c                \   [         R                  " U 5      n Uc  U R                  5       n O3[         R                  " [         R                  " U 5      XR                  S9n [
        R                  " U [         R                  " [        R                  5      5      n[
        R                  " U40 UD6$ )a  Return the union of multiple polygons.

This is an optimized version of union which assumes inputs can be divided into
subsets that do not intersect.

If there is only one such subset, performance can be expected to be worse than
:func:`union_all`.

This function ignores None values when other Geometry elements are present.
If all elements of the given axis are None, an empty GeometryCollection is
returned.

.. versionadded:: 2.1.0

Parameters
----------
geometries : array_like
    Geometries to union.
axis : int, optional
    Axis along which the operation is performed. The default (None)
    performs the operation over all axes, returning a scalar value.
    Axis may be negative, in which case it counts from the last to the
    first axis.
**kwargs
    See :ref:`NumPy ufunc docs <ufuncs.kwargs>` for other keyword arguments.

See Also
--------
coverage_union_all
union_all
disjoint_subset_union

Examples
--------
>>> import shapely
>>> from shapely import Polygon
>>> polygon_1 = Polygon([(0, 0), (0, 1), (1, 1), (1, 0), (0, 0)])
>>> polygon_2 = Polygon([(1, 0), (1, 1), (2, 1), (2, 0), (1, 0)])
>>> shapely.disjoint_subset_union_all([polygon_1, polygon_2]).normalize()
<POLYGON ((0 0, 0 1, 1 1, 2 1, 2 0, 1 0, 0 0))>
>>> shapely.disjoint_subset_union_all([polygon_1, None]).normalize()
<POLYGON ((0 0, 0 1, 1 1, 1 0, 0 0))>
>>> shapely.disjoint_subset_union_all([None, None]).normalize()
<GEOMETRYCOLLECTION EMPTY>
r'   )r   r)   r*   r+   r,   r   r:   r;   r   r<   r   rB   s       r!   r   r     s    ` J'J|%%'
[[JJz"__

 ''BGGL;;<K $$[;F;;r"   )N)NN)__doc__r4   numpyr   shapelyr   r   r   shapely.decoratorsr   r   r   __all__r6   r   r   r   r   r   r   r   r   r	   r
   r   r    r"   r!   <module>rN      s   3   / / 2 {m.@A;*  B;*N {m.@A6,  B6,D vh);<36  =36~ {m.@A64  B64D vh);<B>  =B>\ {m.@A8%  B8%H {F+6HIX2  JX2v  #0 #0^ vh);<=5  ==5@ x1?  1?h x26 :<  :<r"   