ph SSKrSSKrSSKrSSKrSSKJr SSKrSSKr SSK J r J r J r SSKrSSKrSSKrSSKrSSKrSSKJr SSKJrJrJrJrJrJr SSKJr \(aSSKJr \"\R@5r!SSS S S S S SSSS. r"S S1r#1Skr$S S 1r%\RLVs0sHoRN\"URP_M snr)\*"\)RW55r,\RZ"\*"\)R]55\/S9r0"SS\ 5r1\ "\15 S%Sjr2S&Sjr3Sr4Sr5S'Sjr6Sr7S(Sjr8S%Sjr9S(Sjr:Sr;S)S jr<"S!S"\ 5r=S#r>S$r?gs snf)*N) lru_cache)ExtensionArrayExtensionDtyperegister_extension_dtype) BaseGeometry) GEOS_GE_312 HAS_PYPROJ PANDAS_GE_21 PANDAS_GE_22 SHAPELY_GE_21requires_pyproj) SpatialIndex) TransformerPoint LineString LinearRingPolygon MultiPointMultiLineString MultiPolygonGeometryCollection) MISSINGNAGPOINT LINESTRING LINEARRINGPOLYGON MULTIPOINTMULTILINESTRING MULTIPOLYGONGEOMETRYCOLLECTION>rrrdtypec@\rSrSr\rSrSr\S5r \S5r Sr g) GeometryDtype;geometryNc[U[5(d[S[U535eXR:XaU"5$[SUR SUS35e)Nz.'construct_from_string' expects a string, got zCannot construct a 'z' from '') isinstancestr TypeErrortypename__name__)clsstrings gC:\Users\julio\OneDrive\Documentos\Trabajo\IdeasFrscas\Cabanna\env\Lib\site-packages\geopandas/array.pyconstruct_from_string#GeometryDtype.construct_from_string@s]&#&&@fO xx 5L23<<.PQRS Sc[$N) GeometryArray)r1s r3construct_array_type"GeometryDtype.construct_array_typeKsr6) r0 __module__ __qualname____firstlineno__rr.r/na_value classmethodr4r:__static_attributes__r<r6r3r&r&;s9 D DHTTr6r&cU(a#UR(aUR(dgURUR:Xdgg)zo Check if the projection of both arrays is the same. If allow_none is True, empty CRS is treated as the same. TFcrs)leftright allow_nones r3 _check_crsrISs- xxuyy 88uyy  r6cUR(aBURR5n[U5S::aUOSRUSSS/5nOSnUR(aBURR5n[U5S::aUOSRUSSS/5nOSn[R "SUSUS3[ US9 g) zFRaise a CRS mismatch warning with the information on the assigned CRS.2 Nz...zCRS mismatch between the CRS of left geometries and the CRS of right geometries. Use `to_crs()` to reproject one of the input geometries to match the CRS of the other. Left CRS: z Right CRS:   stacklevel)rE to_stringlenjoinwarningswarn UserWarning)rFrGrOleft_srs right_srss r3_crs_mismatch_warnrXas xx88%%'"8}28(3B-QVAW8X yyII'') Y2-I388Ys^U`, such as an authority string (eg "EPSG:4326") or a WKT string. r#Nrfloat64__geo_interface__z&Input must be valid geometry objects: rD)r+r[ndarrayemptyrQobjectr$astypeshapelyis_valid_inputallrappendhasattrr(shaper`r-arrayr9)datarEarroutrbs r3 from_shapelyrus dBJJ ' 'hhs4y/A TaDJJ)3kk&!  ! !# & * * , ,D$ -- 4 233''--d3 4 d 4 "H OPPhhs&)  &&r6cZ[U[5(d [S5eUR$)z?Convert GeometryArray to numpy object array of shapely objects.'geoms' must be a GeometryArray)r+r9 ValueError_data)geomss r3 to_shapelyr{s% e] + +:;; ;;r6c:[[R"XS9US9$)a Convert a list or array of WKB objects to a GeometryArray. Parameters ---------- data : array-like list or array of WKB objects crs : value, optional Coordinate Reference System of the geometry objects. Can be anything accepted by :meth:`pyproj.CRS.from_user_input() `, such as an authority string (eg "EPSG:4326") or a WKT string. on_invalid: {"raise", "warn", "ignore"}, default "raise" - raise: an exception will be raised if a WKB input geometry is invalid. - warn: a warning will be raised and invalid WKB geometries will be returned as None. - ignore: invalid WKB geometries will be returned as None without a warning. - fix: an effort is made to fix invalid input geometries (e.g. close unclosed rings). If this is not possible, they are returned as ``None`` without a warning. Requires GEOS >= 3.11 and shapely >= 2.1.  on_invalidrD)r9rkfrom_wkbrrrEr~s r3rr, ))$FC PPr6c t[U[5(d [S5e[R"U4SU0UD6$)z=Convert GeometryArray to a numpy object array of WKB objects.rwhex)r+r9rxrkto_wkb)rzrkwargss r3rrs4 e] + +:;; >>% 3S 3F 33r6c:[[R"XS9US9$)a Convert a list or array of WKT objects to a GeometryArray. Parameters ---------- data : array-like list or array of WKT objects crs : value, optional Coordinate Reference System of the geometry objects. Can be anything accepted by :meth:`pyproj.CRS.from_user_input() `, such as an authority string (eg "EPSG:4326") or a WKT string. on_invalid : {"raise", "warn", "ignore"}, default "raise" - raise: an exception will be raised if a WKT input geometry is invalid. - warn: a warning will be raised and invalid WKT geometries will be returned as ``None``. - ignore: invalid WKT geometries will be returned as ``None`` without a warning. - fix: an effort is made to fix invalid input geometries (e.g. close unclosed rings). If this is not possible, they are returned as ``None`` without a warning. Requires GEOS >= 3.11 and shapely >= 2.1. r}rD)r9rkfrom_wktrs r3rrrr6c p[U[5(d [S5e[R"U40UD6$)z=Convert GeometryArray to a numpy object array of WKT objects.rw)r+r9rxrkto_wkt)rzrs r3rrs/ e] + +:;; >>% *6 **r6c[R"USS9n[R"USS9nUb[R"USS9n[[R"XU5US9$)a Generate GeometryArray of shapely Point geometries from x, y(, z) coordinates. In case of geographic coordinates, it is assumed that longitude is captured by ``x`` coordinates and latitude by ``y``. Parameters ---------- x, y, z : iterable crs : value, optional Coordinate Reference System of the geometry objects. Can be anything accepted by :meth:`pyproj.CRS.from_user_input() `, such as an authority string (eg "EPSG:4326") or a WKT string. Examples -------- >>> import pandas as pd >>> df = pd.DataFrame({'x': [0, 1, 2], 'y': [0, 1, 2], 'z': [0, 1, 2]}) >>> df x y z 0 0 0 0 1 1 1 1 2 2 2 2 >>> geometry = geopandas.points_from_xy(x=[1, 0], y=[0, 1]) >>> geometry = geopandas.points_from_xy(df['x'], df['y'], df['z']) >>> gdf = geopandas.GeoDataFrame( ... df, geometry=geopandas.points_from_xy(df['x'], df['y'])) Having geographic coordinates: >>> df = pd.DataFrame({'longitude': [-140, 0, 123], 'latitude': [-65, 1, 48]}) >>> df longitude latitude 0 -140 -65 1 0 1 2 123 48 >>> geometry = geopandas.points_from_xy(df.longitude, df.latitude, crs="EPSG:4326") Returns ------- output : GeometryArray rer#rD)r[asarrayr9rkpoints)xyzrEs r3points_from_xyrsRV 1I&A 1I&A} JJq * a0c ::r6c^\rSrSrSr\"5rSSjr\S5r \S5r \S5r \ RS5r S r \S 5rS rS rS rSrSr\S5rSrSSjrSSjr\S5r\S5r\S5r\S5r\S5r\S5r\S5r\S5r\S5r \S5r!Sr"Sr#S r$S!r%S"r&\S#5r'\S$5r(S%r)S&r*\S'5r+\S(5r,S)r-\S*5r.S+r/SS,jr0\S-5r1SS.jr2S/r3S0r4S1r5S2r6S3r7S4r8S5r9SS7jr:SS8jr;S9rSS<jr?SS=jr@SS>jrASS?jrB\CS@5rDSArESBrFSCrGSDrHSErISFrJSGrKSHrLSIrMSJrNSKrOSLrPSMrQSNrRSOrSSPrTSQrUSRrVSSrWSTrXSUrYSVrZSWr[SXr\SYr]SSZjr^SS[jr_SS\jr`SS]jraSS^jrbS_rcS`rdSareSSbjrfScrg\CSd5rhSeriSSfjrjSSgjrkSShjrlSSijrm\nSSjj5ro\nSSkj5rp\Sl5rq\Sm5rr\Sn5rs\So5rt\Sp5ru\Sq5rv\Sr5rw\Ss5rx\St5rySurzSSvjr{SU4Swjjr|SSxjr}SSyjr~SzrSS{\4S|jjrS}r\S~5rSU4Sjjr\SSj5r\SS6S.Sj5rSr\S5rSrSS\S\4SjjrSS\S\4SjjrSSjr\S5rSSjrSSjrSrSrSrSrSrU=r$)r9i;zaClass wrapping a numpy array of Shapely objects. It also holds the array-based implementations. Nc6[XR5(a U(d URnURnOE[U[R 5(d [ S5eURS:Xd [S5eXlSUl X lSUl g)Nzx'data' should be array of geometry objects. Use from_shapely, from_wkb, from_wkt functions to construct a GeometryArray.rz;'data' should be a 1-dimensional array of geometry objects.) r+ __class__rEryr[rgr-ndimrx_crs_sindex)selfrrrEs r3__init__GeometryArray.__init__Cs dNN + +hh::DD"**--M aM    r6chURc[UR5UlUR$)z/Spatial index for the geometries in this array.)rrryrs r3sindexGeometryArray.sindexWs( << ' 3DL||r6cURSL$)aCheck the existence of the spatial index without generating it. Use the `.sindex` attribute on a GeoDataFrame or GeoSeries to generate a spatial index if it does not yet exist, which may take considerable time based on the underlying index implementation. Note that the underlying spatial index may not be fully initialized until the first use. See Also -------- GeoDataFrame.has_sindex Returns ------- bool `True` if the spatial index has been generated or `False` if not. N)rrs r3 has_sindexGeometryArray.has_sindex^s,||4''r6cUR$)aThe Coordinate Reference System (CRS) represented as a ``pyproj.CRS`` object. Returns None if the CRS is not set, and to set the value it :getter: Returns a ``pyproj.CRS`` or None. When setting, the value Coordinate Reference System of the geometry objects. Can be anything accepted by :meth:`pyproj.CRS.from_user_input() `, such as an authority string (eg "EPSG:4326") or a WKT string. )rrs r3rEGeometryArray.crsvsyyr6c[(a&SSKJn U(dSOURU5UlgUb[ R "S[SS9 SUlg)zSet the value of the crs.rCRSNzCannot set the CRS, falling back to None. The CRS support requires the 'pyproj' package, but it is not installed or does not import correctly. The functions depending on CRS will raise an error or may produce unexpected results.rN)r pyprojrfrom_user_inputrrSrTrU)rr_rs r3rErsJ : "$)s/B/B5/IDI  7  DIr6cUR(a[URR(a?[R"S[R "5SR S3[US9 ggg)zGCheck CRS and warn if the planar operation is done in a geographic CRS.z/Geometry is in a geographic CRS. Results from 'rzt' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation. rNN)rE is_geographicrSrTinspectstackfunctionrU)rrOs r3check_geographic_crs"GeometryArray.check_geographic_crssV 88.. MMMMOA&//0199%  /8r6cUR$r8)_dtypers r3r$GeometryArray.dtypes {{r6c URS$)Nr)rprs r3__len__GeometryArray.__len__szz!}r6c[U[R5(aURU$[R R RX5n[URUURS9$NrD) r+numbersIntegralryr]apiindexerscheck_array_indexerr9rE)ridxs r3 __getitem__GeometryArray.__getitem__sV c7++ , ,::c? "ffoo11$<TZZ_$((;;r6c[RRRX5n[ U[R 5(a UR n[ U[R5(aUR R5n[ U[[R-5(a [U5n[ U[5(aD[ U[R5(a [!S5eUR"UR"U'O[ U[$5(d['U5(a['U5(aSnO:[ U[$5(a[U/5R"SnO [)S5e[ U[*[-[R-5(a.[R,"S[.S9nU/USS&X0R"U'OX R"U'O[)SU<35eSUlg)Nz)cannot set a single element with an arrayrzshould be valid geometryrr#z3Value should be either a BaseGeometry or None, got )r]rrrr+Seriesvalues DataFrameflattenlistr[rgrur9rrrxryrr`r-slicerhrir)rkeyr_ value_arrays r3 __setitem__GeometryArray.__setitem__svffoo11$< eRYY ' 'LLE eR\\ * *LL((*E eTBJJ. / / 'E e] + +#w//00 !LMM#kkDJJsO | , ,U E{{E<00$eW-33A6 :;;#ut|bjj899 hhq7 "' A"- 3"' 3EeYO   r6cZ[R"UR5UR4$r8)rkrryrrs r3 __getstate__GeometryArray.__getstate__stzz*DII66r6c[U[5(d8[R"US5nUSUlSUlX lSUlgSU;aURS5US'SU;aSUS'URRU5 g)Nrrrrryr) r+dictrkrrrrybasepop__dict__update)rstaterzs r3 __setstate__GeometryArray.__setstate__s}%&&$$U1X.EaDIDLJDI!&6!2gU" $f MM  'r6cB[R"UR5$r8)rkis_validryrs r3rGeometryArray.is_valid ++r6cB[R"UR5$r8)rkis_valid_reasonryrs r3rGeometryArray.is_valid_reasons&&tzz22r6c[(a [(d [S5e[[R "UR US95$)Nz@Method 'is_valid_coverage' requires shapely>=2.1 and GEOS>=3.12. gap_width)r r ImportErrorboolrkcoverage_is_validryrrs r3is_valid_coverageGeometryArray.is_valid_coverages9 ++R G--djjINOOr6c[(a [(d [S5e[R"UR US9$)NzEMethod 'invalid_coverage_edges' requires shapely>=2.1 and GEOS>=3.12.r)r r rrkcoverage_invalid_edgesryrs r3invalid_coverage_edges$GeometryArray.invalid_coverage_edgess4 ++W --djjINNr6cB[R"UR5$r8)rkis_emptyryrs r3rGeometryArray.is_empty rr6cB[R"UR5$r8)rk is_simpleryrs r3rGeometryArray.is_simple   ,,r6cB[R"UR5$r8)rkis_ringryrs r3rGeometryArray.is_ringstzz**r6cB[R"UR5$r8)rk is_closedryrs r3rGeometryArray.is_closedrr6cB[R"UR5$r8)rkis_ccwryrs r3rGeometryArray.is_ccw~~djj))r6cB[R"UR5$r8)rkhas_zryrs r3rGeometryArray.has_zs}}TZZ((r6cn[(d [S5e[R"UR5$)Nz'has_m' requires shapely>=2.1.)r rrkhas_mryrs r3rGeometryArray.has_m!s&}>? ?}}TZZ((r6c[R"UR5n[[R "[ U5$r8)rk get_type_idrygeometry_type_valuesr[ searchsortedgeometry_type_ids)rress r3 geom_typeGeometryArray.geom_type(s-!!$**-#BOO4Es$KLLr6c`URSS9 [R"UR5$)aReturn the area of the geometries in this array. Raises a UserWarning if the CRS is geographic, as the area calculation is not accurate in that case. Note that the area is calculated in the units of the CRS. Returns ------- np.ndarray of float Area of the geometries. rN)rrkarearyrs r3r GeometryArray.area-s( !!Q!/||DJJ''r6c`URSS9 [R"UR5$)Nr rN)rrklengthryrs r3rGeometryArray.length>s& !!Q!/~~djj))r6cB[R"UR5$r8)rkget_num_coordinatesryrs r3count_coordinatesGeometryArray.count_coordinatesCs**4::66r6cB[R"UR5$r8)rkget_num_geometriesryrs r3count_geometriesGeometryArray.count_geometriesFs))$**55r6cB[R"UR5$r8)rkget_num_interior_ringsryrs r3count_interior_rings"GeometryArray.count_interior_ringsIs--djj99r6cB[R"UR5$r8)rk get_precisionryrs r3rGeometryArray.get_precisionLs$$TZZ00r6c@[R"URUS9$)N)index)rk get_geometryry)rr!s r3r"GeometryArray.get_geometryOs##DJJe<=2.1.rD)r rr9rkconstrained_delaunay_trianglesryrErs r3r0,GeometryArray.constrained_delaunay_trianglesbs8}WX X  2 24:: >DHH  r6cf[[R"UR5URS9$)z7Return the convex hull of the geometries in this array.rD)r9rk convex_hullryrErs r3r3GeometryArray.convex_hulljs$W00<$((KKr6cf[[R"UR5URS9$)z4Return the envelope of the geometries in this array.rD)r9rkenveloperyrErs r3r6GeometryArray.envelopeos$W--djj9txxHHr6cf[[R"UR5URS9$)zEReturn the minimum rotated rectangle of the geometries in this array.rD)r9rkoriented_enveloperyrErs r3minimum_rotated_rectangle'GeometryArray.minimum_rotated_rectanglet"W66tzzBQQr6cf[[R"UR5URS9$r)r9rkget_exterior_ringryrErs r3exteriorGeometryArray.exteriorxr<r6cf[[R"UR5URS9$r)r9rkextract_unique_pointsryrErs r3rB#GeometryArray.extract_unique_points|s"W::4::FDHHUUr6c j[[R"URUUUUS9URS9$)N) quad_segs join_style mitre_limitrD)r9rk offset_curveryrE)rdistancerErFrGs r3rHGeometryArray.offset_curves9   #%'    r6c6Sn/nURHBn[USS5nUbUR[U55 M/SnURS5 MD U(a[R "SSS9 [ R"[U5[S9nX%SS&U$)NF interiorsTzUOnly Polygon objects have interior rings. For other geometry types, None is returned.rrNr#) rygetattrrnrrSrTr[rhrQri)r has_non_poly inner_ringsrbinterior_ring_seqrrs r3rLGeometryArray.interiorss  JJD 'k4 @  ,""4(9#:; $ ""4(  MM4  xxK(7Q r6cd[[R"URUS9URS9$)N tolerancerD)r9rkremove_repeated_pointsryrErrTs r3rU$GeometryArray.remove_repeated_pointss*  * *4:: K  r6cf[[R"UR5URS9$r)r9rkpoint_on_surfaceryrErs r3representative_point"GeometryArray.representative_points"W55djjAtxxPPr6cf[[R"UR5URS9$r)r9rkminimum_bounding_circleryrErs r3r]%GeometryArray.minimum_bounding_circles"W<=2.1.rSrD)r rr9rkmaximum_inscribed_circleryrErVs r3r`&GeometryArray.maximum_inscribed_circles;}QR R  , ,TZZ9 M  r6cB[R"UR5$r8)rkminimum_bounding_radiusryrs r3rc%GeometryArray.minimum_bounding_radiuss..tzz::r6cB[R"UR5$r8)rkminimum_clearanceryrs r3rfGeometryArray.minimum_clearances((44r6c[(d [S5e[[R"UR 5UR S9$)Nz/'minimum_clearance_line' requires shapely>=2.1.rD)r rr9rkminimum_clearance_lineryrErs r3ri$GeometryArray.minimum_clearance_lines3}OP PW;;DJJGTXXVVr6cf[[R"UR5URS9$r)r9rk normalizeryrErs r3rlGeometryArray.normalizes"W..tzz:IIr6Fc[(d [S5e[[R"UR US9UR S9$)Nz('orient_polygons' requires shapely>=2.1.) exterior_cwrD)r rr9rkorient_polygonsryrE)rros r3rpGeometryArray.orient_polygonss:}HI I  # #DJJK Hdhh  r6c0n[(a XS'X#S'OUS:wa [S5e[[R"UR 40UD6UR S9$)Nmethodkeep_collapsedlineworkz:Only the 'linework' method is supported for shapely < 2.1.rD)r rxr9rk make_validryrE)rrsrtrs r3rvGeometryArray.make_valids\ =%8 '5# $# PW// EfE488TTr6cf[[R"UR5URS9$r)r9rkreverseryrErs r3ryGeometryArray.reverses W__TZZ8dhhGGr6ch[[R"URU5URS9$r)r9rk segmentizeryrE)rmax_segment_lengths r3r|GeometryArray.segmentizes+   tzz+= >  r6cf[[R"UR5URS9$r)r9rkforce_2dryrErs r3rGeometryArray.force_2ds"W--djj9txxHHr6cd[[R"URUS9URS9$)N)rrD)r9rkforce_3dryrE)rrs r3rGeometryArray.force_3ds$W--djjA>DHHMMr6cd[[R"URXS9URS9$)N include_zrD)r9rk transformryrE)rtransformationrs r3rGeometryArray.transforms*   djj. N  r6cd[[R"URUS9URS9$)N)directedrD)r9rk line_mergeryrE)rrs r3rGeometryArray.line_merges)   tzzH =488  r6cd[[R"URXS9URS9$)N) grid_sizemoderD)r9rk set_precisionryrE)rrrs r3rGeometryArray.set_precisions*  ! !$** M  r6c 6[U[5(ac[U5[U5:wa%S[U5S[U53n[U5e[ X5(d [ XSS9 UR n[[U5"UR U40UD6$)Nz&Lengths of inputs do not match. Left: z , Right: rN) r+r9rQrxrIrXryrMrk)oprFrGrmsgs r3_binary_methodGeometryArray._binary_methods e] + +4yCJ& YKyU >!o%d**"41=KKEw#DJJ@@@r6c&URSX5$)Ncoversrrothers r3rGeometryArray.covers""8T99r6c&URSX5$)N covered_byrrs r3rGeometryArray.covered_by""<==r6c&URSX5$)Ncontainsrrs r3rGeometryArray.contains"":t;;r6c&URSX5$)Ncontains_properlyrrs r3rGeometryArray.contains_properlys""#6DDr6c&URSX5$)Ncrossesrrs r3rGeometryArray.crosses""9d::r6c&URSX5$)Ndisjointrrs r3rGeometryArray.disjointrr6c&URSX5$)Nequalsrrs r3 geom_equalsGeometryArray.geom_equals rr6c&URSX5$)N intersectsrrs r3rGeometryArray.intersects#rr6c&URSX5$)Noverlapsrrs r3rGeometryArray.overlaps&rr6c&URSX5$)Ntouchesrrs r3rGeometryArray.touches)rr6c&URSX5$)Nwithinrrs r3rGeometryArray.within,rr6cBURSS9 URSXUS9$)NrNdwithin)rIrr)rrrIs r3rGeometryArray.dwithin/s+ !!Q!/""9dH"MMr6c$URSXUS9$)N equals_exactrSrrrrTs r3geom_equals_exactGeometryArray.geom_equals_exact3s"">4)"TTr6cR[(d [S5eURSX5$)Nz.'geom_equals_identical' requires shapely>=2.1.equals_identical)r rrrs r3geom_equals_identical#GeometryArray.geom_equals_identical6s'}NO O""#5tCCr6c j[[R"URXX45URS9$r)r9rk clip_by_rectryrE)rxminyminxmaxymaxs r3rGeometryArray.clip_by_rect?s+  T D$((  r6cJ[URSX5URS9$)N differencerDr9rrErs r3rGeometryArray.differenceDs%    d :  r6cJ[URSX5URS9$)N intersectionrDrrs r3rGeometryArray.intersectionI%    <$((  r6cJ[URSX5URS9$)Nsymmetric_differencerDrrs r3r"GeometryArray.symmetric_differenceNs&    6 D$((  r6cJ[URSX5URS9$)NunionrDrrs r3rGeometryArray.unionSs T00$FDHHUUr6cJ[URSX5URS9$)N shortest_linerDrrs r3rGeometryArray.shortest_lineVs%    =488  r6c H[URSXUS9URS9$)NsnaprSrDrrs r3rGeometryArray.snap[s*   y  Itxx  r6cJ[URSX5URS9$)N shared_pathsrDrrs r3rGeometryArray.shared_paths`rr6cDURSS9 URSX5$)NrrNrIrrs r3rIGeometryArray.distanceis& !!Q!/"":t;;r6c HURSS9 UR"SX40UD6$)NrrNhausdorff_distancerrrrs r3r GeometryArray.hausdorff_distancems, !!Q!/""#7OOOr6c HURSS9 UR"SX40UD6$)NrrNfrechet_distancerrs r3rGeometryArray.frechet_distanceqs, !!Q!/""#5tMfMMr6c [U[[-5(aUS:XdURSS9 [ [ R "URU4SU0UD6URS9$)Nrr rNrErD) r+intrZrr9rkbufferryrE)rrI resolutionrs r3rGeometryArray.bufferusZ8S5[11h!m  % % % 3 NN4::x P: P P  r6cURSS9 [[R"URXS9UR S9$)Nr rN normalizedrD)rr9rkline_interpolate_pointryrE)rrIrs r3 interpolateGeometryArray.interpolate}s; !!Q!/  * *4::x W  r6cd[[R"URXS9URS9$)N)preserve_topologyrD)r9rksimplifyryrE)rrTrs r3rGeometryArray.simplifys/    I    r6c[(a [(d [S5e[[R "UR XS9URS9$)Nz9'simplify_coverage' requires shapely>=2.1 and GEOS>=3.12.)simplify_boundaryrD)r r rr9rkcoverage_simplifyryrE)rrTrs r3simplify_coverageGeometryArray.simplify_coveragesJ ++K   % % I    r6c[U[5(a URn[R"URXS9$)Nr)r+r9ryrkline_locate_point)rrrs r3projectGeometryArray.projects/ e] + +KKE((URRr6c[U[5(a URn[R"URU5$r8)r+r9ryrkrelaters r3rGeometryArray.relates- e] + +KKE~~djj%00r6c[U[5(a URn[R"URX5$r8)r+r9ryrkrelate_pattern)rrpatterns r3rGeometryArray.relate_patterns/ e] + +KKE%%djj%AAr6cV[R"S[SS9 UR5$)NzNThe 'unary_union' attribute is deprecated, use the 'union_all' method instead.rrN)rSrTDeprecationWarning union_allrs r3 unary_unionGeometryArray.unary_unions(  2   ~~r6ctUS:waUb[SUS35eUS:Xa [R"UR5$US:Xa[R"URUS9$US:XaA[ (a [ (d [S5e[R"UR5$[SUS 35e) Nunaryz'grid_size is not supported for method 'z'.coverage)rdisjoint_subsetz=Method 'disjoin_subset' requires shapely>=2.1 and GEOS>=3.12.zMethod 'z?' not recognized. Use 'coverage', 'unary' or 'disjoint_subset'.) rxrkcoverage_union_allryrr r rdisjoint_subset_union_all)rrsrs r3rGeometryArray.union_alls W !6FvhbQR R Z --djj9 9 w $$TZZ9E E ( (!Mkk!S44TZZ@ @6(#%% r6cB[R"UR5$r8)rkintersection_allryrs r3r$GeometryArray.intersection_alls'' 33r6c/nUHOnUbUR(aUnO$[[RU5"U/UQ70UD6nUR U5 MQ [ R "[U5[S9nXGSS&U$)Nr#) rrMrkaffinityrnr[rhrQri)rrFargsrrtrbrrrs r3_affinity_methodGeometryArray._affinity_methodsq D|t}}g..3DJ4J6J JJsO  xxD 0Q r6c`[URSURU5URS9$)Naffine_transformrDr9r)ryrE)rmatrixs r3r,GeometryArray.affine_transforms-  ! !"4djj& I  r6c b[URSURXU5URS9$)N translaterDr-)rxoffyoffzoffs r3r1GeometryArray.translates.  ! !+tzz4t L  r6c ^[URSURXUS9URS9$)Nrotateorigin use_radiansrDr-)rangler9r:s r3r7GeometryArray.rotates8  ! !$**e  "    r6c ^[URSURXX4S9URS9$)Nscale)r9rDr-)rxfactyfactzfactr9s r3r>GeometryArray.scales8  ! !U5 "    r6c ^[URSURXX4S9URS9$)Nskewr8rDr-)rxsysr9r:s r3rDGeometryArray.skews8  ! ! B6 "    r6chSSKJn URc [S5eUbUR U5nO UbUR U5nO [S5eURR U5(aU$[URUSS9n[URUR5n[XQS9$)uTransform all geometries to a different coordinate reference system. Transform all geometries in a GeometryArray to a different coordinate reference system. The ``crs`` attribute on the current GeometryArray must be set. Either ``crs`` or ``epsg`` may be specified for output. This method will transform all points in all objects. It has no notion of projecting entire geometries. All segments joining points are assumed to be lines in the current projection, not geodesics. Objects crossing the dateline (or other projection boundary) will have undesirable behavior. Parameters ---------- crs : pyproj.CRS, optional if `epsg` is specified The value can be anything accepted by :meth:`pyproj.CRS.from_user_input() `, such as an authority string (eg "EPSG:4326") or a WKT string. epsg : int, optional if `crs` is specified EPSG code specifying output projection. Returns ------- GeometryArray Examples -------- >>> from shapely.geometry import Point >>> from geopandas.array import from_shapely, to_wkt >>> a = from_shapely([Point(1, 1), Point(2, 2), Point(3, 3)], crs=4326) >>> to_wkt(a) array(['POINT (1 1)', 'POINT (2 2)', 'POINT (3 3)'], dtype=object) >>> a.crs # doctest: +SKIP Name: WGS 84 Axis Info [ellipsoidal]: - Lat[north]: Geodetic latitude (degree) - Lon[east]: Geodetic longitude (degree) Area of Use: - name: World - bounds: (-180.0, -90.0, 180.0, 90.0) Datum: World Geodetic System 1984 - Ellipsoid: WGS 84 - Prime Meridian: Greenwich >>> a = a.to_crs(3857) >>> to_wkt(a) array(['POINT (111319.490793 111325.142866)', 'POINT (222638.981587 222684.208506)', 'POINT (333958.47238 334111.171402)'], dtype=object) >>> a.crs # doctest: +SKIP Name: WGS 84 / Pseudo-Mercator Axis Info [cartesian]: - X[east]: Easting (metre) - Y[north]: Northing (metre) Area of Use: - name: World - 85°S to 85°N - bounds: (-180.0, -85.06, 180.0, 85.06) Coordinate Operation: - name: Popular Visualisation Pseudo-Mercator - method: Popular Visualisation Pseudo Mercator Datum: World Geodetic System 1984 - Ellipsoid: WGS 84 - Prime Meridian: Greenwich rrzICannot transform naive geometries. Please set a crs on the object first.zMust pass either crs or epsg.T always_xyrD) rrrErxr from_epsg is_exact_sameTransformerFromCRSrryr9)rrEepsgr transformernew_datas r3to_crsGeometryArray.to_crssJ  88 8  ?%%c*C  --%C<= = 88 ! !# & &K(3$G TZZ)>)>?X//r6c SSKJn SSKJn SSKJn UR (d [S5eURupVpxUR R(a/[R"XW/5n [R"Xh/5n O[UR SSS9n U RXVXx5upVpx[R"Xh/5n XW:a(US - n[R"XW/5n U S -S -S - n O[R"XW/5n U"UU"U U U U S 9S 9n URU SR5$![ a [S 5ef=f)u}Return the estimated UTM CRS based on the bounds of the dataset. .. versionadded:: 0.9 .. note:: Requires pyproj 3+ Parameters ---------- datum_name : str, optional The name of the datum to use in the query. Default is WGS 84. Returns ------- pyproj.CRS Examples -------- >>> import geodatasets >>> df = geopandas.read_file( ... geodatasets.get_path("geoda.chicago_commpop") ... ) >>> df.geometry.values.estimate_utm_crs() # doctest: +SKIP Name: WGS 84 / UTM zone 16N Axis Info [cartesian]: - E[east]: Easting (metre) - N[north]: Northing (metre) Area of Use: - name: Between 90°W and 84°W, northern hemisphere between equator and 84°N,... - bounds: (-90.0, 0.0, -84.0, 84.0) Coordinate Operation: - name: UTM zone 16N - method: Transverse Mercator Datum: World Geodetic System 1984 ensemble - Ellipsoid: WGS 84 - Prime Meridian: Greenwich rr)AreaOfInterest)query_utm_crs_infoz$crs must be set to estimate UTM CRS.z EPSG:4326TrIih)west_lon_degreesouth_lat_degreeeast_lon_degreenorth_lat_degree) datum_namearea_of_interestzUnable to determine UTM CRS)rr pyproj.aoirTpyproj.databaserUrE RuntimeError total_boundsrr[meanrMtransform_boundsrKcode IndexError) rr[rrTrUminxminymaxxmaxyx_centery_centerrO utm_crs_lists r3estimate_utm_crsGeometryArray.estimate_utm_crs_s;N -6xxEF F!%!2!2D 88 ! !ww|,Hww|,H-TXX{dSK%0%A%AD& "Dww|,H{ 77D<0%^s2c977D<0)!+ (!) (!)   >==a!5!56 6 ><= = >s D::EcURUR5)S:HR5(aURnUR 5(aRU)n[ R "U[[ RS9n[R"URU5X2'U$[R"UR5$Sn[U5e)z9Return the x location of point geometries in a GeoSeries.rr$ fill_valuez5x attribute access only provided for Point geometries) rr`rmranyr[ full_likerZnanrkget_xryrxrrhnonemptycoordsmessages r3rGeometryArray.x NNDIIK< (G 3 8 8 : :MMEyy{{!6heO#*==H1E#F  }}TZZ00MGW% %r6cURUR5)S:HR5(aURnUR 5(aRU)n[ R "U[[ RS9n[R"URU5X2'U$[R"UR5$Sn[U5e)z9Return the y location of point geometries in a GeoSeries.rroz5y attribute access only provided for Point geometries) rr`rmrrqr[rrrZrsrkget_yryrxrus r3rGeometryArray.yrzr6cURUR5)S:HR5(aURnUR 5(aRU)n[ R "U[[ RS9n[R"URU5X2'U$[R"UR5$Sn[U5e)z9Return the z location of point geometries in a GeoSeries.rroz5z attribute access only provided for Point geometries) rr`rmrrqr[rrrZrsrkget_zryrxrus r3rGeometryArray.zrzr6c[(d [S5eURUR5)S:HR 5(aUR nUR 5(aRU)n[R"U[[RS9n[R"URU5X2'U$[R"UR5$Sn[U5e)z;Return the m coordinate of point geometries in a GeoSeries.z'm' requires shapely>=2.1.rroz5m attribute access only provided for Point geometries)r rrr`rmrrqr[rrrZrsrkget_mryrxrus r3mGeometryArray.ms}:; ; NNDIIK< (G 3 8 8 : :MMEyy{{!6heO#*==H1E#F  }}TZZ00MGW% %r6cB[R"UR5$r8)rkboundsryrs r3rGeometryArray.boundsrr6c [U5S:XaR[R"[R[R[R[R/5$URn[ R "5 [ R"SS[5 [R"[R"USS2S45[R"USS2S45[R"USS2S45[R"USS2S4545sSSS5 $!,(df  g=f)NrignorezAll-NaN slice encounteredrr) rQr[rqrsrrScatch_warningsfilterwarningsRuntimeWarningnanminnanmax)rbs r3r`GeometryArray.total_boundss t9>88RVVRVVRVVRVV<= = KK  $ $ &  # #6 88IIa1g&IIa1g&IIa1g&IIa1g&  ' & &s B"D// D=c.URR$r8)rysizers r3rGeometryArray.sizeszzr6cUR4$r8)rrs r3rpGeometryArray.shapes |r6c,[UR5$r8)rQrprs r3rGeometryArray.ndims4::r6cZ[URR5URS9$r)r9rycopyr)rr(rs r3rGeometryArray.copy"sTZZ__.DII>>r6c&SSKJn U(a<Ub[R"U5(aSnO[ U5(d [ S5eU"UR XUS9nU(aUcSU[R"U5)'[XPRS9$)Nr)takez&provide geometry or None as fill value) allow_fillrprD) pandas.api.extensionsrr]r`rcr-ryrkrlr9rE)rindicesrrprresults r3rGeometryArray.take&sy. !RWWZ%8%8! (44 HIIdjj'ZX *,6:FG**622 3V22r6c n>[(a[(dUbX5S'OX5S'[TU] "SXUS.UD6$)N limit_area)rslimitrr<)r r super_pad_or_backfill)rrsrrrrrs r3rGeometryArray._pad_or_backfill5sB < %(2|$#-< w'XvXQWXXr6c6Ub [S5eUR5nU(aUR5nOUnUR5(dU$Ub;U[ U5:a,UR 5U:nUR5(aSXW'[U5(aS/nOP[ U5(aU/nO<[U[5(aXnO"[S[U5R<35e[R"U[S9nXRU'U$)a Fill NA values with geometry (or geometries) or using the specified method. Parameters ---------- value : shapely geometry object or GeometryArray If a geometry value is passed it is used to fill all missing values. Alternatively, an GeometryArray 'value' can be given. It's expected that the GeometryArray has the same length as 'self'. method : {'backfill', 'bfill', 'pad', 'ffill', None}, default None Method to use for filling holes in reindexed Series pad / ffill: propagate last valid observation forward to next valid backfill / bfill: use NEXT valid observation to fill gap limit : int, default None The maximum number of entries where NA values will be filled. copy : bool, default True Whether to make a copy of the data before filling. If False, then the original should be modified and no new memory should be allocated. Returns ------- GeometryArray Nz)fillna with a method is not yet supportedFzT'value' parameter must be None, a scalar geometry, or a GeoSeries, but you passed a r#)NotImplementedErrorr`rrqrQcumsumrcr+r9r-r.r0r[rriry) rr_rsrrmask new_valuesmodify value_arrs r3fillnaGeometryArray.fillnaAs6  %&QR Ryy{ JJxxzz   T!2[[]U*Fzz||$ ;;FE  ' 'GE } - -KE$$(K$8$8#;= JJuF3 !*r6cd[U[5(aU(aUR5$U$[RR R U5(a[RR RU5(dw[U5n[RR RU5n[U[R5(a[R"X4S9$URUSS9$U(d[R"XS9$[R"XUS9$)az Cast to a NumPy array with 'dtype'. Parameters ---------- dtype : str or dtype Typecode or data-type to which the array is cast. copy : bool, default True Whether to copy the data, even if not necessary. If False, a copy is made only if the old dtype does not match the new dtype. Returns ------- array : ndarray NumPy ndarray with 'dtype' for its dtype. r#F)rr$r)r+r&rr]rtypesis_string_dtypeis_object_dtyper pandas_dtype StringDtyperqrjr[r)rr$r string_valuespd_dtypes r3rjGeometryArray.astype}s$ e] + +yy{" VV\\ ) )% 0 09U9U : : #4LMvv||007H(BNN33xx >> ''E': : zz$44xx==r6cB[R"UR5$)z8Boolean NumPy array indicating if each value is missing.)rk is_missingryrs r3r`GeometryArray.isnas!!$**--r6dropnac [U5nSSKJnJn U"U5R US9nU"[ [ R"URR5SUR555UlU$)z Compute a histogram of the counts of non-null values. Parameters ---------- dropna : bool, default True Don't include counts of NaN Returns ------- pd.Series r)Indexr)rN) rpandasrr value_countsrr[wherer!r`)rrrrrrs r3rGeometryArray.value_countssb$(,,F,;  RXXfll//14F G   r6c&SSKJn U"U5up#U$)zWCompute the ExtensionArray of unique values. Returns ------- uniques : ExtensionArray r) factorize)rr)rr_uniquess r3uniqueGeometryArray.uniques %t_ r6c.URR$r8)rynbytesrs r3rGeometryArray.nbytesszz   r6cH>[TU]X5nURUlU$)a Shift values by desired number. Newly introduced missing values are filled with ``self.dtype.na_value``. Parameters ---------- periods : int, default 1 The number of periods to shift. Negative values are allowed for shifting backwards. fill_value : object, optional (default None) The scalar value to use for newly introduced missing values. The default is ``self.dtype.na_value``. Returns ------- GeometryArray Shifted. Notes ----- If ``self`` is empty or ``periods`` is 0, a copy of ``self`` is returned. If ``periods > len(self)``, then an array of size len(self) is returned, with all values filled with ``self.dtype.na_value``. )rshiftrE)rperiodsrpshiftedrs r3rGeometryArray.shifts#>'-4hh r6cH[U[5(aU/n[U5$)a Construct a new ExtensionArray from a sequence of scalars. Parameters ---------- scalars : Sequence Each element will be an instance of the scalar type for this array, ``cls.dtype.type``. dtype : dtype, optional Construct for this particular dtype. This should be a Dtype compatible with the ExtensionArray. copy : boolean, default False If True, copy the underlying data. Returns ------- ExtensionArray )r+rru)r1scalarsr$rs r3_from_sequenceGeometryArray._from_sequences#* g| , ,iGG$$r6rc[U5$)a Construct a new ExtensionArray from a sequence of strings. Parameters ---------- strings : Sequence Each element will be an instance of the scalar type for this array, ``cls.dtype.type``. dtype : dtype, optional Construct for this particular dtype. This should be a Dtype compatible with the ExtensionArray. copy : bool, default False If True, copy the underlying data. Returns ------- ExtensionArray )r)r1stringsr$rs r3_from_sequence_of_strings'GeometryArray._from_sequence_of_stringss*  r6c [U5nUS4$)aReturn an array and missing value suitable for factorization. Returns ------- values : ndarray An array suitable for factorization. This should maintain order and be a supported dtype (Float64, Int64, UInt64, String, Object). By default, the extension array is cast to object dtype. na_value : object The value in `values` to consider missing. This will be treated as NA in the factorization routines, so it will be coded as `na_sentinal` and not included in `uniques`. By default, ``np.nan`` is used. N)r)rvalss r3_values_for_factorize#GeometryArray._values_for_factorize0s d|Tzr6c([XRS9$)a$ Reconstruct an ExtensionArray after factorization. Parameters ---------- values : ndarray An integer ndarray with the factorized values. original : ExtensionArray The original ExtensionArray that factorize was called on. See Also -------- pandas.factorize ExtensionArray.factorize rD)rrE)r1roriginals r3_from_factorizedGeometryArray._from_factorizedCs"LL11r6czSSKJn URS:Xa[R"/SS9$UR nUR 5nUR5U-nUR 5(ahUR5n[R"U)5SnUR(aXSnO [RRSS5nXuU'OUnU(aGURnUSUS- S-n USUS- S-n U[R"U *U *X/5- nOS nU"XXS 9n U(aSX'U $) zReturn values for sorting. Returns ------- ndarray The transformed values should maintain the ordering between values within the array. See Also -------- ExtensionArray.argsort r)_hilbert_distanceuint32r#rg{Gz?rrN)r`)geopandas.tools.hilbert_curverrr[rqrrqr`rnonzerorkr(rr`) rr mask_empty has_emptyrrzrrbr`r2r3 distancess r3_values_for_argsort!GeometryArray._values_for_argsortVs- D 99>88Bh/ /]] NN$ yy{Z' 88::IIKEjj$'*G||AJ'''--a3$KE !--L Ol1o5=D Ol1o5=D BHHtedUD%?@ @LL%eG %&I !r6skipnareturnc[S5eNz%geometries have no minimum or maximumr-rrs r3argminGeometryArray.argmin?@@r6c[S5errrs r3argmaxGeometryArray.argmaxrr6c(^U(aSSKnURRmTcUR(a!URR(aSmOSmO[ R "5 [ R"S[S9 XR5)SSRup4pVSSS5 SWs=::aS ::a6O O3SWs=::aS ::a&O O#S Ws=::aS ::aO OS Ws=::aS ::aO OSmOSmU4S j$[$!,(df  N_=f) aReturn a formatting function for scalar values. This is used in the default '__repr__'. The returned formatting function receives instances of your scalar type. Parameters ---------- boxed: bool, default False An indicated for whether or not your array is being printed within a Series, DataFrame, or Index (True), or just by itself (False). This may be useful if you want scalar values to appear differently within a Series versus on its own (e.g. quoted or not). Returns ------- Callable[[Any], str] A callable that gets instances of the scalar type and returns a string. By default, :func:`repr` is used when ``boxed=False`` and :func:`str` is used when ``boxed=True``. rNrr r)category iLrViZc.>[R"UTS9$)N)rounding_precision)rkr)rb precisions r3*GeometryArray._formatter..st Rr6) geopandasoptionsdisplay_precisionrE is_projectedrSr simplefilterrr`r`repr)rboxedrrrrrrs @r3 _formatterGeometryArray._formatters.  !));;I 88xx,,$% $% "002 --hP15yy{l1CCR1H1U1U.D3,,!T0S0 D.B. D.B.%& %& R R !32s -=D Dc[R"UVs/sHo"RPM sn5n[U[ U5S9$s snf)zuConcatenate multiple array. Parameters ---------- to_concat : sequence of this type Returns ------- ExtensionArray rD)r[ concatenateryr9_get_common_crs)r1 to_concatgarrs r3_concat_same_typeGeometryArray._concat_same_types:~~)<)Bxx)<=Ty'ABB=sAc US;a[URU5"US9$[S[U5RSUR SUS35e)N)rqrm)keepdimsr*z ' with dtype z does not support reduction ')rMryr-r.r0r$)rr/rrrs r3_reduceGeometryArray._reduces\ > !4::t,h? ?T ##$M$**>++/& 3  r6cU(a7UbU[R"S5:XaURR5$UR$)zlReturn the data as a numpy array. This is the numpy array interface. Returns ------- values : numpy array ri)r[r$ryr)rr$rs r3 __array__GeometryArray.__array__s6 U]erxx/A&A::??$ $zzr6c^U4Sjn[U[R[R-[R-5(a[ $TnU"U5n[ U5[ U5:wa [S5e[XE5VVs/sH upgU"Xg5PM nnn[R"U[S9nU$s snnf)Nc>[U5(dG[U[5(d.[RR R U5(aUnU$U/[T5-nU$r8)rcr+rr]rr is_list_likerQ)paramovaluesrs r3convert_values,GeometryArray._binop..convert_valuessV&u--5.11RVV\\5N5Nu5U5UN!'CI-Nr6zLengths must match to comparer#) r+r]rrrNotImplementedrQrxzipr[rr) rrrrlvaluesrvaluesarrs ` r3_binopGeometryArray._binops  eRYY1BLL@ A A! ! ' w<3w< '<= ='*'&;<&;FQr!x&;<jjD) =s B<cBURU[R5$r8)r&operatoreqrs r3__eq__GeometryArray.__eq__{{5(++..r6cBURU[R5$r8)r&r)ners r3__ne__GeometryArray.__ne__r-r6c[U5(a^XRRLd'[XRR5(dUcUR5R 5$gX:HR 5$)zReturn for `item in self`.F)r`r$r@r+r.rq)ritems r3 __contains__GeometryArray.__contains__sZ :: +++dJJOO44<yy{(( !!##r6)rryrrrEr8))roundg@F)ruT)r) valid_output))T)rN)r6r6r6)centerF)?r=r=r<)r6r6r<FNN)zWGS 84)FN)NNT)NNNT)rN)NF)TF)r0r=r>r?__doc__r&rrpropertyrrrEsetterrr$rrrrrrrrrrrrrrrrrr rrrrrr"r%r(r-r0r3r6r:r?rBrHrLrUrZr]r`rcrfrirlrprvryr|rrrrr staticmethodrrrrrrrrrrrrrrrrrrrrrrrrIrrrrrr rrrrrr$r)r,r1r7r>rDrrQrlrrrrrr`rrprrrrrrjr`rrrrrrArrrrrrrrr rrrr&r+r0r4rB __classcell__)rs@r3r9r9;s _F( ((.   ZZ$ <#Z7 (&,,3PO,,--++--**)))) MM(( **76:1=IIIIV LLIIRRRV  0 QX ;5W J  UH IN    A A:><E;<:><;:NUD    V   <PN     S 1 B &4      Y0Y0vP>P>l & & & & & &&&$**0? 3 9= Y:x'>R. > !!!N%%09=E!!,&22$9vATASAATASA5n C C  6// $ $r6r9c4UVs/sH6oR5R5(aURcM4UPM8 nn/nUH0nURU;dMURUR5 M2 UVs/sH oDcMUPM nnUVs/sHoDRPM nn[ U5S:Xag[ U5S:Xa0[ U5S:wa[ R"SUSS3SS9 US$[SUS35es snfs snfs snf) NrrzJCRS not set for some of the concatenation inputs. Setting output's CRS as z$ (the single non-null crs provided).rrNz:Cannot determine common CRS for concatenation inputs, got zH. Use `to_crs()` to transform geometries to the same CRS before merging.) r`rmrErnr/rQrSrTrx)arr_seqr unique_crsrsrE crs_not_nonenamess r3rr#s$PGbGGIMMOOrGGPJ 77* $   cgg &$.A:CC:LA!- .#XXE . <A <A z?a  MM++08*566   A  DUGLQ Q /QB .s3D D D DDcH[R"U5n[R"U5n[R"X)SS9nU"USS2S4USS2S45n[R "X)R 5[R"U5R5X2)'[R"XSS9nU"USS2S4USS2S4USS2S45n[R "XR 5[R"U5R5X2'U$)NFrrrTr) rkrr[ empty_likeget_coordinatesset_coordinatesrrqT)rrfuncrrrw new_coords_zcoords_zs r3rrCs MM$ E ]]4 F  $ $T&\U CFq!t fQTl3L,, V RXXl355F6N&&t{dCHAAAGL++ BHH\244FM Mr6r9)rr8)Nraiser>)@rrr)rS functoolsrnumpyr[rr]rrrrrkshapely.affinityshapely.geometry shapely.ops shapely.wktshapely.geometry.baser_compatr r r r r rrrrrfrom_crsrM_namesPOLYGON_GEOM_TYPESLINE_GEOM_TYPESPOINT_GEOM_TYPES GeometryTyper_r/ type_mappingrkeysrrqrrirr&rIrXr`rcrur{rrrrrr9rr)ps0r3rcsr .! "";#7#78  (".  0A\*181E1EF1EA'1EF **,-xx\%8%8%: ;6JN*' 8,*''TQ24Q2+0;fe$Ne$P/@[8Gs#E