Ch%zSrSSKrSSKrSSKrSSKJrJrJ r SSK J r J r SSK JrJrJrJrJr SrS SjrS Sjrg) z+Module to clip vector data using GeoPandas.N) MultiPolygonPolygonbox) GeoDataFrame GeoSeries)LINE_GEOM_TYPESPOINT_GEOM_TYPESPOLYGON_GEOM_TYPES _check_crs_crs_mismatch_warnc[RRRU5=(a* [ U[ [ -[-[-5(+$)at 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 gC:\Users\julio\OneDrive\Documentos\Trabajo\Ideas Frescas\venv\Lib\site-packages\geopandas/tools/clip.py_mask_is_list_like_rectanglersF ::   ( ( . z lY&0<?84c [U5nU(a [U6nOUnURURR USUS9nUR S:gnUR 5(dU$[U[5(aUR5nU(a?URRUR"U6URXgR4'OURRURU5URXgR4'OVUR5nU(aURUR"U6Xv'O URURU5Xv'U(aXwR )nU$)a+ 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#hh -DQG &&'1N     '<((,,.   ''7DDdK KK(E(EE F   ''7DDTJ KK(E(EE F ,,. &-nn^&D&Q&QSW&XG #&-nn^&D&Q&QRV&WG #+++, NrcX[U[[-5(d[S[ U535e[ U5n[U[[-[ -[-5(dU(d[S[ U535eU(a[U5S:wa [S5e[U[[-5(a[X5(d [XSS9 [U[[-5(a URnO;U(aUnO1UR(d URO[R4S-nURnUSUS::a$USUS::aUS US::a US US::dUR S S$[U[[-5(aUR"R%5nOUn['XUS 9nU(GaRUR(S :HR+5n UR(S :HR+5n U =(a U (+n U (a[,R."S SS9 U$[1UR(R3[45R+5UR(R3[65R+5UR(R3[85R+5/5n [1UR(R3[45R+5UR(R3[65R+5UR(R3[85R+5/5n X:nU S :a[,R."SSS9 U$U (dU(aUR(R SnU (aUR;SS9nU[4;a.UR<UR(R3[45nU$U[6;a,UR<UR(R3[65nU$)a 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 rcliprPlsyB c<)3 4 4I$s)UVV8> t\I5?,N O O% 004T | =  Ta W  $ y011#$$ sQ 7$ y011$$  '+mm4;;"&&QG 1+ #'!* *CqkWQZ'gajHQK.Gxx|$ y011 //1  !#4@G",,0DDIIK*>>CCE (>->  MM8 R NE #MM&&'9:>>@MM&&7;;=MM&&'78<<>   #%%**+=>BBD%%**?;??A%%**+;<@@B  *rXs6177-*@F_r