h%vdZddlZddlZddlZddlmZmZm Z ddl m Z m Z ddl mZmZmZmZmZdZd dZd dZy) z+Module to clip vector data using GeoPandas.N) MultiPolygonPolygonbox) GeoDataFrame GeoSeries)LINE_GEOM_TYPESPOINT_GEOM_TYPESPOLYGON_GEOM_TYPES _check_crs_crs_mismatch_warnctjjj|xr&t |t t ztztz S)a Check if the input mask is list-like and not an instance of specific geometric types. Parameters ---------- mask : GeoDataFrame, GeoSeries, (Multi)Polygon, list-like Polygon vector layer used to clip ``gdf``. Returns ------- bool True if `mask` is list-like and not an instance of `GeoDataFrame`, `GeoSeries`, `Polygon`, or `MultiPolygon`, otherwise False. ) pandasapitypes is_list_like isinstancerrrr)masks G/var/www/html/immo/lib/python3.12/site-packages/geopandas/tools/clip.py_mask_is_list_like_rectanglersF ::   ( ( . z lY&0<?84ct|}|r t|}n|}|j|jj |d|}|j dk7}|j s|St|tr|j}|r@|jj|j||j||jf<n|jj|j||j||jf<nS|j}|r |j|j|||<n!|j|j|||<|r||j }|S)as Clip geometry to the polygon/rectangle extent. Clip an input GeoDataFrame to the polygon extent of the polygon parameter. Parameters ---------- gdf : GeoDataFrame, GeoSeries Dataframe to clip. mask : (Multi)Polygon, list-like Reference polygon/rectangle for clipping. sort : boolean, default False If True, the results will be sorted in ascending order using the geometries' indexes as the primary key. Returns ------- GeoDataFrame The returned GeoDataFrame is a clipped subset of gdf that intersects with polygon/rectangle. intersects) predicatesortPoint)rrilocsindexquery geom_typeanyrrcopygeometryvalues clip_by_rectloc_geometry_column_name intersectionis_empty)gdfrrclipping_by_rectangleintersection_polygongdf_subnon_point_maskclippeds r_clip_gdf_with_maskr/)sf29>"Dz#hh -DQG &&'1N    '<(,,. D  ''7DDdK KK(E(EE F   ''7DDTJ KK(E(EE F ,,. &Qgnn^&D&Q&QSW&XGN #&-nn^&D&Q&QRV&WGN #7+++, Nrct|ttzstdt |t |}t|ttzt ztzs|stdt ||rt|dk7r tdt|ttzrt||st||dt|ttzr |j}n1|r|}n,|js |jntjfdz}|j}|d|dkr!|d|dkr|d |dkr |d |dks|j d dSt|ttzr|j"j%}n|}t'||| }|r8|j(d k(j+} |j(d k(j+} | xr| } | rt-j.d d|St1|j(j3t4j+|j(j3t6j+|j(j3t8j+g} t1|j(j3t4j+|j(j3t6j+|j(j3t8j+g} | | k}| d kDrt-j.dd|S| s|r|j(j d}| r|j;d}|t4vr.|j<|j(j3t4}|S|t6vr,|j<|j(j3t6}|S)as Clip points, lines, or polygon geometries to the mask extent. Both layers must be in the same Coordinate Reference System (CRS). The ``gdf`` will be clipped to the full extent of the clip object. If there are multiple polygons in mask, data from ``gdf`` will be clipped to the total boundary of all polygons in mask. If the ``mask`` is list-like with four elements ``(minx, miny, maxx, maxy)``, a faster rectangle clipping algorithm will be used. Note that this can lead to slightly different results in edge cases, e.g. if a line would be reduced to a point, this point might not be returned. 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 ---------- gdf : GeoDataFrame or GeoSeries Vector layer (point, line, polygon) to be clipped to mask. mask : GeoDataFrame, GeoSeries, (Multi)Polygon, list-like Polygon vector layer used to clip ``gdf``. The mask's geometry is dissolved into one geometric feature and intersected with ``gdf``. If the mask is list-like with four elements ``(minx, miny, maxx, maxy)``, ``clip`` will use a faster rectangle clipping (:meth:`~GeoSeries.clip_by_rect`), possibly leading to slightly different results. keep_geom_type : boolean, default False If True, return only geometries of original type in case of intersection resulting in multiple geometry types or GeometryCollections. If False, return all resulting geometries (potentially mixed-types). sort : boolean, default False If True, the results will be sorted in ascending order using the geometries' indexes as the primary key. Returns ------- GeoDataFrame or GeoSeries Vector data (points, lines, polygons) from ``gdf`` clipped to polygon boundary from mask. See Also -------- GeoDataFrame.clip : equivalent GeoDataFrame method GeoSeries.clip : equivalent GeoSeries method Examples -------- Clip points (grocery stores) with polygons (the Near West Side community): >>> import geodatasets >>> chicago = geopandas.read_file( ... geodatasets.get_path("geoda.chicago_health") ... ) >>> near_west_side = chicago[chicago["community"] == "NEAR WEST SIDE"] >>> groceries = geopandas.read_file( ... geodatasets.get_path("geoda.groceries") ... ).to_crs(chicago.crs) >>> groceries.shape (148, 8) >>> nws_groceries = geopandas.clip(groceries, near_west_side) >>> nws_groceries.shape (7, 8) z/'gdf' should be GeoDataFrame or GeoSeries, got zJ'mask' should be GeoDataFrame, GeoSeries,(Multi)Polygon or list-like, got zIIf 'mask' is list-like, it must have four values (minx, miny, maxx, maxy)) stacklevelrN)rGeometryCollectionzKkeep_geom_type can not be called on a GeoDataFrame with GeometryCollection.z>keep_geom_type can not be called on a mixed type GeoDataFrame.F) index_parts)rrr TypeErrortyperrrlenr r total_boundsr(boundsnpnanrr" union_allr/rr warningswarnsumisinr rr exploder%)r)rkeep_geom_typerr*box_maskbox_gdf combined_maskr.geomcoll_concat geomcoll_orignew_collectionorig_types_totalclip_types_total more_types orig_types rcliprPlspB c<)3 4I$s)UVV8> t\I5?,N O% 004T | =  Ta W  $ y01#t$ sDQ 7$ y01$$  '+mm4;;"&&QG 1+ #'!* *CqkWQZ'gajHQK.Gxx|$ y01 //1  !#}4@G",,0DDIIK*>>CCE (>->  MM8 R NE #MM&&'9:>>@MM&&7;;=MM&&'78<<>   #%%**+=>BBD%%**?;??A%%**+;<@@B  *,<rXs6177-*@F_r