Ch3BSSKrSSKJr SSKrSSKrSSKJrJ r SSK J r SSK J r JrJrJrJr SrSrSrS rS rS rSS jrS rg)N)reduce) GeoDataFrame GeoSeries) PANDAS_GE_30)LINE_GEOM_TYPESPOINT_GEOM_TYPESPOLYGON_GEOM_TYPES _check_crs_crs_mismatch_warncURS:Xdo[(a3SUR;aURSSS9nUR S5nU$SUR;aURSSSS9 UR SSS9 U$)ztEnsure that the geometry column is called 'geometry'. If another column with that name exists, it will be dropped. geometryaxisTrinplace)r)_geometry_column_namercolumnsdroprename_geometry)dfs jC:\Users\julio\OneDrive\Documentos\Trabajo\Ideas Frescas\venv\Lib\site-packages\geopandas/tools/overlay.py_ensure_geometry_columnrs  # #z 1 <RZZ'WWZaW0##J/B IRZZ' D9   z4  8 IcURRURSSS9up#URS:GaUURS:GaDURR U5nUR SSS9 URR U5nUR SSS9 UR U5nURR[5nXgR5URU'[R"X#S.5nUn UR SS9nUR SS9nURURUR SS 9S SS 9n U RURUR SS 9S SS S9n [#XUR$S9$UR&SSRUR&SSRURR(SS 9SSS S9n [*R,U S '[*R,U S 'U U R.RURR(5R15URR(/-$)z8Overlay Intersection operation used in overlay function. intersectsT predicatesortrrr)__idx1__idx2rrrr!)left_on right_indexr"_1_2)r$r%suffixesr crsN left_indexr%r))sindexqueryr sizetake reset_index intersection geom_typeisinr make_validlocpd DataFramemergerrrr+ilocnamenpnanrtolist) df1df2idx1idx2leftright intersectionspoly_ixpairs_intersectgeom_intersectdfinterresults r_overlay_intersectionrL#s.!!#,,,T!RJD yy1}Q||  & dD1 !!$' tT2))%0 ))../AB%2%;%F%F%H '",,$'GH&oo4o(oo4o(!'' HHS..QH 7(  -- HHS..QH 7!   G#''JJ"1## HHRaL  cll//a  8! $  66x66x NN   1 1 2 9 9 ;s||?P?P>Q Q  rc XURRURSSS9up#[R"USS9upE[R "X5SS5n[ URR5Vs/sHnXt;aURS5O/PM nn/n [URU5HGup[SU /[URRU 5-5n U RU 5 MI [XRUR S 9n U R"R%[&5nXR)5U R*U'XR,)R/5nX R,)R/5nUUUR0'U$s snf) z6Overlay Difference operation used in overlay function.rTr) return_indexrNrc$URU5$)N) difference)xys r%_overlay_difference..as ar)indexr+)r.r/r r=uniquesplitranger0popziprlistr;appendrrUr+r4r5r r6r7is_emptycopyr)r@rArBrC idx1_uniqueidx1_unique_indices idx2_splitidxsidxnew_ggeom neighboursnew differencesrG geom_diffdfdiffs r_overlay_differencerkSs}!!#,,,T!RJD')yyD'I$K$AB 78J**+ +C!/ qR7+   E d3 (4&4 8I8I*8U3V*V   S 4 E@K##(();>@KOOG1112779I &&& ' , , .F+4F6 ' '( M# s9"F'ct[X5n[X5n[U5nURVs/sHnXBR;aUOUS3PM nnXSl[R "X#/SSS9n[ UR5nURS5 URS5 URUS9$s snf)z4Overlay Identity operation used in overlay function.r'TF ignore_indexrr r) rLrkrrr8concatr[remover\reindex)r@rAdfintersection dfdifferencecol new_columnsrKrs r_overlay_identityrwms*34N&s0L*<8L  '''C,,,SE*<''YY5Du UF>))*G NN: NN: >>'> **s B5cz[X5n[X5n[[U55US'[[U55US'[RUS'[RUS'[ U5n[ U5nUR USS/SSS9nURR5nSUl URURR5S4URURR5'URSS/S S S 9 URS S S 9 [XEURS 9nU$)z@Overlay Symmetric Difference operation used in overlay function.r!r"outerr&)onhowr)r geometry_2 geometry_1rTrr r*)rkrXlenr=r>rr:r}r^r<r7isnullrr2rr+)r@rAdfdiff1dfdiff2dfsymr s r_overlay_symmetric_diffrs(!#+G!#+Gc'l+GHc'l+GHGHGH%g.G%g.G MMXx(g   E$$&HHM.3ii !<//HLL!!((*+ JJ l+!TJB 4. sww ?E Lrc[X5n[X5n[R"X#/SSS9n[ UR 5nUR S5 URS5 URUS9$)z1Overlay Union operation used in overlay function.TFrmr ro) rLrr8rpr[rrqr\rr)r@rArJrdfunionrs r_overlay_unionrsb#C-G #C -Eii(t%HG7??#G NN: NN: ??7? ++rct^/SQnX%;a[SUSU35e[U[5(d[U[5(a [S5e[ X5(d [ XSS9 UcSnSnOS n[ X/5HupxURR[5R5n URR[5R5n URR[5R5n [XU /5S :dM[S US -S 35e US :XaURn URn U SU S::a$U SU S::aU S U S::a U S U S::dURSSR!URSSR#UR$R&S S9SSSS9nUUR(R#UR$R&5R+5UR$R&/-$U4SjnU(aURRSnU"U5nU"U5n[,R."5 [,R0"SSS9 US:Xa [3X5nOGUS :Xa [5X5nO5US:Xa [7X5nO#US:Xa [9X5nOUS:Xa [;X5nUS;aWR#SS/S SS9 SSS5 U(a [=WWU5nWR?SSS9 U$!,(df  N4=f) aFPerform spatial overlay between two GeoDataFrames. Currently only supports data GeoDataFrames with uniform geometry types, i.e. containing only (Multi)Polygons, or only (Multi)Points, or a combination of (Multi)LineString and LinearRing shapes. Implements several methods that are all effectively subsets of the union. See the User Guide page :doc:`../../user_guide/set_operations` for details. Parameters ---------- df1 : GeoDataFrame df2 : GeoDataFrame how : string Method of spatial overlay: 'intersection', 'union', 'identity', 'symmetric_difference' or 'difference'. keep_geom_type : bool If True, return only geometries of the same geometry type as df1 has, if False, return all resulting geometries. Default is None, which will set keep_geom_type to True but warn upon dropping geometries. make_valid : bool, default True If True, any invalid input geometries are corrected with a call to make_valid(), if False, a `ValueError` is raised if any input geometries are invalid. Returns ------- df : GeoDataFrame GeoDataFrame with new set of polygons and attributes resulting from the overlay Examples -------- >>> from shapely.geometry import Polygon >>> polys1 = geopandas.GeoSeries([Polygon([(0,0), (2,0), (2,2), (0,2)]), ... Polygon([(2,2), (4,2), (4,4), (2,4)])]) >>> polys2 = geopandas.GeoSeries([Polygon([(1,1), (3,1), (3,3), (1,3)]), ... Polygon([(3,3), (5,3), (5,5), (3,5)])]) >>> df1 = geopandas.GeoDataFrame({'geometry': polys1, 'df1_data':[1,2]}) >>> df2 = geopandas.GeoDataFrame({'geometry': polys2, 'df2_data':[1,2]}) >>> geopandas.overlay(df1, df2, how='union') df1_data df2_data geometry 0 1.0 1.0 POLYGON ((2 2, 2 1, 1 1, 1 2, 2 2)) 1 2.0 1.0 POLYGON ((2 2, 2 3, 3 3, 3 2, 2 2)) 2 2.0 2.0 POLYGON ((4 4, 4 3, 3 3, 3 4, 4 4)) 3 1.0 NaN POLYGON ((2 0, 0 0, 0 2, 1 2, 1 1, 2 1, 2 0)) 4 2.0 NaN MULTIPOLYGON (((3 4, 3 3, 2 3, 2 4, 3 4)), ((4... 5 NaN 1.0 MULTIPOLYGON (((2 3, 2 2, 1 2, 1 3, 2 3)), ((3... 6 NaN 2.0 POLYGON ((3 5, 5 5, 5 3, 4 3, 4 4, 3 4, 3 5)) >>> geopandas.overlay(df1, df2, how='intersection') df1_data df2_data geometry 0 1 1 POLYGON ((2 2, 2 1, 1 1, 1 2, 2 2)) 1 2 1 POLYGON ((2 2, 2 3, 3 3, 3 2, 2 2)) 2 2 2 POLYGON ((4 4, 4 3, 3 3, 3 4, 4 4)) >>> geopandas.overlay(df1, df2, how='symmetric_difference') df1_data df2_data geometry 0 1.0 NaN POLYGON ((2 0, 0 0, 0 2, 1 2, 1 1, 2 1, 2 0)) 1 2.0 NaN MULTIPOLYGON (((3 4, 3 3, 2 3, 2 4, 3 4)), ((4... 2 NaN 1.0 MULTIPOLYGON (((2 3, 2 2, 1 2, 1 3, 2 3)), ((3... 3 NaN 2.0 POLYGON ((3 5, 5 5, 5 3, 4 3, 4 4, 3 4, 3 5)) >>> geopandas.overlay(df1, df2, how='difference') geometry df1_data 0 POLYGON ((2 0, 0 0, 0 2, 1 2, 1 1, 2 1, 2 0)) 1 1 MULTIPOLYGON (((3 4, 3 3, 2 3, 2 4, 3 4)), ((4... 2 >>> geopandas.overlay(df1, df2, how='identity') df1_data df2_data geometry 0 1 1.0 POLYGON ((2 2, 2 1, 1 1, 1 2, 2 2)) 1 2 1.0 POLYGON ((2 2, 2 3, 3 3, 3 2, 2 2)) 2 2 2.0 POLYGON ((4 4, 4 3, 3 3, 3 4, 4 4)) 3 1 NaN POLYGON ((2 0, 0 0, 0 2, 1 2, 1 1, 2 1, 2 0)) 4 2 NaN MULTIPOLYGON (((3 4, 3 3, 2 3, 2 4, 3 4)), ((4... See Also -------- sjoin : spatial join GeoDataFrame.overlay : equivalent method Notes ----- Every operation in GeoPandas is planar, i.e. the potential third dimension is not taken into account. )r3unionidentitysymmetric_differencerPz `how` was 'z' but is expected to be in z4overlay currently only implemented for GeoDataFrames stacklevelNTFrrz contains mixed geometry types.r3rrr&r,c>UR5nURR[5R 5(aUR R )nURnT(aNURX4R5URX4'UR5(a [USSS9nU$UR5(a[SUR5S35eU$)NPolygonF)r4keep_geom_type_warningz,You have passed make_valid=False along with zo invalid input geometries. Use make_valid=True or make sure that all geometries are valid before using overlay.)r^r4r5r allr is_validrr7r6any_collection_extract ValueErrorsum)rmaskrur6s r _make_validoverlay.._make_valid8s WWY <<  / 0 4 4 6 6KK(((D**C$&FF49$5$@$@$Bty!88::,iB  Bxxzl#66  rignorezCRS mismatch between the CRS)messagerPrrr)r3rrrr!r"rr ) r isinstancerNotImplementedErrorr r enumerater4r5r rrrr total_boundsr;r:rr r<rr?warningscatch_warningsfilterwarningsrkrLrrrwrr2)r@rAr{keep_geom_typer6 allowed_howsrir poly_check lines_check points_checkbox_gdf1box_gdf2rKrr4s ` roverlayrsrL ;se+F|nUVV#y!!ZY%?%?! B   c  32!%!&C:&\\&&'9:>>@ ll''8<<> ||(()9:>>@  6 7! ;%1q5'1P&QR R ' n####qkXa[(x{hqk/I1+!,8A;(1+3MXXbq\''! !!#,,"3"3!!< % (F ##CLL$5$56==?3<&B&B&D D %%%3 .7.A.H.H =*+!" \\!_N ZZ((--j9 :F <<?2Ka1A5;Q ;345     Mr)r3NT)r functoolsrnumpyr=pandasr8 geopandasrrgeopandas._compatrgeopandas.arrayrrr r r rrLrkrwrrrrrrrsL-*"- `4+06 ,}@;r