rbi%xdZddlZddlZddlZddlmZmZm Z ddl m Z m Z ddl mZmZmZmZmZdZd dZd d ZdS) z+Module to clip vector data using GeoPandas.N) MultiPolygonPolygonbox) GeoDataFrame GeoSeries)LINE_GEOM_TYPESPOINT_GEOM_TYPESPOLYGON_GEOM_TYPES _check_crs_crs_mismatch_warnctjj|o-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 d/var/www/html/mdtn/previsions/meteo_cartes/venv/lib/python3.11/site-packages/geopandas/tools/clip.py_mask_is_list_like_rectanglersI :  ( ( . . z lY&0<?884Fct|}|r t|}n|}|j|j|d|}|jdk}|s|St|tru| }|r*|j j |j ||j ||jf<n|j j |||j ||jf<nR| }|r|j |j |||<n#|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/)s~29>>$"Dz#h -DQQG &'1N     '<((X,,..  D '7DdK K(EE F F  '7DDTJJ K(EE F F ,,..  X&Qgn^&D&QSW&XGN # #&-n^&D&Q&QRV&W&WGN #-7++, Nrct|ttzstdt |t |}t|ttzt ztzs!|stdt ||r"t|dkrtdt|ttzr"t||st||dt|ttzr|j }n#|r|}n|j s|j ntjfdz}|j }|d|dkr6|d|dkr$|d |dkr|d |dks|jd dSt|ttzr|j}n|}t'||| }|rO|jd k} |jd k} | o| } | rt-jd dnt1|jt4|jt6|jt8g} t1|jt4|jt6|jt8g} | | k}| d krt-jddn| s|r|jjd}| r|d}|t4vr+|j|jt4}n3|t6vr*|j|jt6}|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 rcliprPls9B c<)3 4 4WU$s))UUVVV8>> t\I5?,N O O %  =04T  = =    Ta W   $ y0118#t$$ 8 sDQ 7 7 7 7$ y011 G$ G '+mF4;;"&QG 1+ # #'!* *C*CqkWQZ''gajHQK.G.Gx|$ y011 //11  !#}4@@@G.S",0DDIIKK*>>CCEE (>-> ( S M8      #M&&'9::>>@@M&&77;;==M&&'788<<>>    #%**+=>>BBDD%**?;;??AA%**+;<<@@BB   *,<rXs117777777777--------*@@@@F______r