(&hjdZddlZddlZddlmZmZmZddlm Z m Z m Z gdZ e dge e ddZe dge e dd Ze d ge e dd Ze dge e dd Ze d ge e dd Ze dge e ddZe dd ge e ddZeZe dZe d ge e ddZe de dZe de dddZdS)z-Set-theoretic operations on geometry objects.N)Geometry GeometryTypelib)deprecate_positionalmultithreading_enabled requires_geos) coverage_unioncoverage_union_all differencedisjoint_subset_uniondisjoint_subset_union_all intersectionintersection_allsymmetric_differencesymmetric_difference_all unary_unionunion union_all grid_size)categoryc |7tj|stdtj|||fi|Stj||fi|S)agReturn 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 ` 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)])) >>> shapely.difference(line, LineString()) >>> 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() >>> box1 = shapely.box(0.1, 0.2, 2.1, 2.1) >>> shapely.difference(box1, box2, grid_size=1) N.grid_size parameter only accepts scalar values)npisscalar ValueErrorrdifference_precr abrkwargss Q/var/www/html/reinick/venv/lib/python3.11/site-packages/shapely/set_operations.pyr r 'sdn{9%% OMNN N"1a==f=== >!Q ) )& ) ))c |7tj|stdtj|||fi|Stj||fi|S)aReturn 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 ` 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)])) >>> box1 = shapely.box(0, 0, 2, 2) >>> box2 = shapely.box(1, 1, 3, 3) >>> shapely.intersection(box1, box2).normalize() >>> box1 = shapely.box(0.1, 0.2, 2.1, 2.1) >>> shapely.intersection(box1, box2, grid_size=1) Nr)rrrrintersection_precrrs r!rrpsed{9%% OMNN N$Q9?????  Aq + +F + ++r"axisc tj|}||}ntj|||j}t j|fi|S)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 ` 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]) >>> shapely.intersection_all([[line1, line2, None]], axis=1).tolist() [] >>> shapely.intersection_all([line1, None]) Nr%start)rasarrayravelrollaxisndimrr geometriesr%r s r!rrs]^J''J |%%'' [$joNNN   5 5f 5 55r"c |7tj|stdtj|||fi|Stj||fi|S)aBReturn 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 ` 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)])) >>> box1 = shapely.box(0, 0, 2, 2) >>> box2 = shapely.box(1, 1, 3, 3) >>> shapely.symmetric_difference(box1, box2).normalize() >>> box1 = shapely.box(0.1, 0.2, 2.1, 2.1) >>> shapely.symmetric_difference(box1, box2, grid_size=1) Nr)rrrrsymmetric_difference_precrrs r!rrsed{9%% OMNN N,Q9GGGGG  #Aq 3 3F 3 33r"c tjdtdtj|}||}ntj|||j}tj |fi|S)aReturn 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 ` 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]) >>> shapely.symmetric_difference_all([[line1, line2, None]], axis=1).tolist() [] >>> shapely.symmetric_difference_all([line1, None]) >>> shapely.symmetric_difference_all([None, None]) 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.) stacklevelNr') warningswarnDeprecationWarningrr)r*r+r,rrr-s r!rr9sn M O  J''J |%%'' [$joNNN  ' = =f = ==r"c |7tj|stdtj|||fi|Stj||fi|S)aMerge 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 ` 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)])) >>> 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() >>> box1 = shapely.box(0.1, 0.2, 2.1, 2.1) >>> shapely.union(box1, box2, grid_size=1) Nr)rrrr union_precrrs r!rrsch{9%% OMNN N~aI88888 9Q $ $V $ $$r"c tj|}||}ntj|||j}t j|tjtj }|6tj |stdt j ||fi|St j |fi|S)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 ` 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]) >>> shapely.union_all([[line1, line2, None]], axis=1).tolist() [] >>> box1 = shapely.box(0, 0, 2, 2) >>> box2 = shapely.box(1, 1, 3, 3) >>> shapely.union_all([box1, box2]).normalize() >>> box1 = shapely.box(0.1, 0.2, 2.1, 2.1) >>> shapely.union_all([box1, box2], grid_size=1) >>> shapely.union_all([None, Point(0, 1)]) >>> shapely.union_all([None, None]) >>> shapely.union_all([]) Nr'r)rr)r*r+r,rcreate_collectionintcrGEOMETRYCOLLECTIONrrunary_union_precr)r.rr%r collectionss r!rrsRJ''J |%%'' [$joNNN 'BGL;<<K{9%% OMNN N#KEEfEEE ?; 1 1& 1 11r"c t||gfi|S)arMerge multiple polygons into one. This is an optimized version of union which assumes the polygons to be non-overlapping. Parameters ---------- a, b : Geometry or array_like Geometry or geometries to merge (union). **kwargs See :ref:`NumPy ufunc docs ` 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() Union with None returns same polygon >>> shapely.coverage_union(polygon_1, None).normalize() )r rrr s r!r r /sB q!f / / / //r"c :tj|}||}n.tjtj|||j}t j|tjtj }t j |fi|S)aReturn 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 ` 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() >>> shapely.coverage_union_all([polygon_1, None]).normalize() >>> shapely.coverage_union_all([None, None]).normalize() Nr') rr)r*r+r,rr:r;rr<r r.r%r r>s r!r r \shJ''J |%%'' [ Jz " "Z_   'BGL;<<K  k 4 4V 4 44r"z3.12.0c t|ts|t|ts|nt|ts|tj||}n\t|ts|tj||}n/t |t |krt dt ||gfddi|S)aMerge 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 ` 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() Union with None returns same polygon: >>> shapely.disjoint_subset_union(polygon_1, None).normalize() Nz(Arrays a and b must have the same lengthr%r) isinstancerr full_likelenrr r@s r!r r sR 1h E191h$-#$9 Ax EAI LA   Ax EAI LA   Q3q66  CDDD $aV > >! >v > >>r")r%c :tj|}||}n.tjtj|||j}t j|tjtj }t j |fi|S)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 ` 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() >>> shapely.disjoint_subset_union_all([polygon_1, None]).normalize() >>> shapely.disjoint_subset_union_all([None, None]).normalize() Nr') rr)r*r+r,rr:r;rr<r rBs r!r r s`J''J |%%'' [ Jz " "Z_   'BGL;<<K  $[ ; ;F ; ;;r")N)NN)__doc__r4numpyrshapelyrrrshapely.decoratorsrrr__all__r6r rrrrrrrr r r r r"r!rNs/33//////////   2{m.@AAA;*;*;*BA;*N{m.@AAA6,6,6,BA6,Dvh);<<<363636=<36~{m.@AAA646464BA64Dvh);<<<B>B>B>=<B>\{m.@AAA8%8%8%BA8%H{F+6HIIIX2X2X2JIX2v  0 0 0Xvh);<<<=5=5=5=<=5@x1?1?1?hx26:<:<:<:<:<:<: