a 1$e@sddlZddlZddlZddlZddlmZddlZddlZ ddl m Z m Z m Z ddlZddlZddlZddlmZddlZddlZddlmZmZddlmZddlmZz ddlZWneydZYn0dd lm Z!dd lm"Z#dd l$m%Z%eej&Z'Gd d d e Z(e e(d,ddZ)d-ddZ*ddZ+ddZ,ddZ-d.ddZ.ddZ/d/ddZ0d0d d!Z1d1d"d#Z2d$d%Z3d2d&d'Z4Gd(d)d)e Z5d*d+Z6dS)3N) lru_cache)ExtensionArrayExtensionDtyperegister_extension_dtype) BaseGeometry)CRS Transformer)AreaOfInterest)query_utm_crs_info)_compat) _vectorized)_get_sindex_classc@s2eZdZeZdZejZe ddZ e ddZ dS) GeometryDtypegeometrycCsDt|tstdt|n"||jkr.|Std|j|dS)Nz0'construct_from_string' expects a string, got {}z!Cannot construct a '{}' from '{}') isinstancestr TypeErrorformattypename__name__)clsstringrK/var/www/html/django/DPS/env/lib/python3.9/site-packages/geopandas/array.pyconstruct_from_string*s    z#GeometryDtype.construct_from_stringcCstSN) GeometryArray)rrrrconstruct_array_type7sz"GeometryDtype.construct_array_typeN) r __module__ __qualname__rrrnpnanna_value classmethodrrrrrrr%s rFcCs(|r|jr|jsdS|j|jks$dSdS)z{ 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_nonerrr _check_crs?s   r+cCs|jr8|j}t|dkr |nd|dddg}nd}|jrt|j}t|dkr\|nd|dddg}nd}tjd||t|ddS)zP Raise 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: {0} Right CRS: {1}  stacklevel)r' to_stringlenjoinwarningswarnr UserWarning)r(r)r0Zleft_srsZ right_srsrrr_crs_mismatch_warnMs  ( $r7cCs"tjr |Stjs|St|SdS)zY Convert internal representation (PyGEOS or Shapely) to external Shapely object. N)compatUSE_SHAPELY_20 USE_PYGEOS vectorizedZ_pygeos_to_shapelygeomrrr_geom_to_shapelyps r>cCs"tjr |Stjs|St|SdS)zY Convert external Shapely object to internal representation (PyGEOS or Shapely). N)r8r9r:r;Z_shapely_to_pygeosr<rrr_shapely_to_geom|s r?cCs$tjrt|tjtfSt|tSdSr)r8r:rpygeosZGeometryrr<rrr_is_scalar_geometrysrAcCstt||dS)a Convert a list or array of shapely objects to a GeometryArray. Validates the elements. Parameters ---------- data : array-like list or array of shapely 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. r&)rr; from_shapelydatar'rrrrBsrBcCst|tstdt|jS)zI Convert GeometryArray to numpy object array of shapely objects. 'geoms' must be a GeometryArray)rr ValueErrorr; to_shapely_data)geomsrrrrGs rGcCstt||dS)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. r&)rr;from_wkbrCrrrrJsrJcKs*t|tstdtj|jfd|i|S)zG Convert GeometryArray to a numpy object array of WKB objects. rEhex)rrrFr;to_wkbrH)rIrKkwargsrrrrLs rLcCstt||dS)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. r&)rr;from_wktrCrrrrNsrNcKs&t|tstdtj|jfi|S)zG Convert GeometryArray to a numpy object array of WKT objects. rE)rrrFr;to_wktrH)rIrMrrrrOs rOcCstt||||dS)aN 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 r&)rr;points_from_xy)xyzr'rrrrPs+rPcsNeZdZdZeZdddZeddZeddZ ed d Z ed d Z e j d d Z ddZ eddZddZddZddZddZddZeddZeddZed d!Zed"d#Zed$d%Zed&d'Zed(d)Zed*d+Zed,d-Zed.d/Zed0d1Zd2d3Zed4d5Z d6d7Z!ed8d9Z"d:d;Z#edd?Z%ddCdDZ&edEdFZ'ddHdIZ(dJdKZ)dLdMZ*dNdOZ+dPdQZ,dRdSZ-dTdUZ.dVdWZ/e0dXdYZ1dZd[Z2d\d]Z3d^d_Z4d`daZ5dbdcZ6dddeZ7dfdgZ8dhdiZ9djdkZ:dldmZ;dndoZdtduZ?dvdwZ@dxdyZAdzd{ZBd|d}ZCd~dZDddZEddZFdddZGdddZHdddZIdddZJddZKddZLddZMdddZNdddZOdddZPdddZQdddZRdddZSeddZTeddZUeddZVeddZWeddZXeddZYeddZZeddZ[ddZ\dddZ]ddZ^dddZ_dddZ`ddZadebddd„ZcddĄZdeddƄZedfddɄ Zfegddd˄Zhdd̈́ZiegddτZjddфZkdebeldҜddԄZmdebeldҜddքZnddd؄ZoegddڄZpddd܄ZqdddބZrddZsddZtddZuddZvZwS)rzj Class wrapping a numpy array of Shapely objects and holding the array-based implementations. NcCsbt||jr|s|j}|j}n(t|tjs4tdn|jdksFtd||_d|_ ||_d|_ dS)Nzx'data' should be array of geometry objects. Use from_shapely, from_wkb, from_wkt functions to construct a GeometryArray.r z;'data' should be a 1-dimensional array of geometry objects.) r __class__r'rHr"ndarrayrndimrF_crs_sindex)selfrDr'rrr__init__s    zGeometryArray.__init__cCstjdtdd|jS)NaAccessing the underlying geometries through the `.data` attribute is deprecated and will be removed in GeoPandas 1.0. You can use `np.asarray(..)` or the `to_numpy()` method instead. Note that if you are using PyGEOS and using this attribute to get an array of PyGEOS geometries, those other methods will always return an array of Shapely geometries. Accessing the underlying PyGEOS geometries directly is deprecated, and you should migrate to use Shapely >= 2.0 instead.r/)r4r5DeprecationWarningrHrYrrrrD)s  zGeometryArray.datacCs|jdurt|j|_|jSr)rXrrHr]rrrsindex9s zGeometryArray.sindexcCs |jduS)a_Check 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)rXr]rrr has_sindex?szGeometryArray.has_sindexcCs|jS)a The 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. )rWr]rrrr'Ws zGeometryArray.crscCs|sdnt||_dS)zSets the value of the crsN)rfrom_user_inputrW)rYvaluerrrr'escCs2|jr.|jjr.tjdtdjt|ddS)zFCheck CRS and warn if the planar operation is done in a geographic CRSzGeometry is in a geographic CRS. Results from '{}' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation. r r/N) r' is_geographicr4r5rinspectstackfunctionr6)rYr0rrrcheck_geographic_crsjs z"GeometryArray.check_geographic_crscCs|jSr)_dtyper]rrrdtypewszGeometryArray.dtypecCs |jdS)Nr)shaper]rrr__len__{szGeometryArray.__len__cCs>t|tjrt|j|Stjj||}t |j||j dSNr&) rnumbersIntegralr>rHpdapiindexerscheck_array_indexerrr')rYidxrrr __getitem__~s zGeometryArray.__getitem__cCsRtjj||}t|tjr"|j}t|tjr8|j}t|t t j frPt |}t|t r|t|tjrntd|j|j|<nt|tst|r8t|rd}n$t|trt |gjd}ntdt|tt t j fr,t jdtd}t|g|dd<Wdn1s0Y||j|<n ||j|<ntdt|d|_dS)Nz)cannot set a single element with an arrayrzshould be valid geometryr rhz5Value should be either a BaseGeometry or None, got %s)rnrorprqrSeriesvalues DataFrameflattenlistr"rUrBrrlrmrFrHrr;isnarsliceemptyobjectr8ignore_shapely2_warningsrrX)rYkeyraZ value_arrayrrr __setitem__s6        .   zGeometryArray.__setitem__cCs:tjrt|j|jfStjr0t|j|jfS|jSdSr) r8r9shapelyrLrHrWr:r@__dict__r]rrr __getstate__s zGeometryArray.__getstate__cCst|ts6t|d}|d|_d|_||_d|_nJd|vrL|d|d<t j rdt |d|d<d|vrtd|d<|j |dS)Nrr rDrHrW)rdictr;rJrWrXrHbasepopr8r:rBrupdate)rYstaterIrrr __setstate__s  zGeometryArray.__setstate__cCs t|jSr)r;is_validrHr]rrrrszGeometryArray.is_validcCs t|jSr)r;is_emptyrHr]rrrrszGeometryArray.is_emptycCs t|jSr)r; is_simplerHr]rrrrszGeometryArray.is_simplecCs t|jSr)r;is_ringrHr]rrrrszGeometryArray.is_ringcCs t|jSr)r; is_closedrHr]rrrrszGeometryArray.is_closedcCs t|jSr)r;has_zrHr]rrrrszGeometryArray.has_zcCs t|jSr)r; geom_typerHr]rrrrszGeometryArray.geom_typecCs|jddt|jSNr/)rfr;arearHr]rrrrs zGeometryArray.areacCs|jddt|jSr)rfr;lengthrHr]rrrrs zGeometryArray.lengthcCstt|j|jdSrk)rr;boundaryrHr'r]rrrrszGeometryArray.boundarycCs"|jddtt|j|jdS)Nrr/r&)rfrr;centroidrHr'r]rrrrs zGeometryArray.centroidcCstj|j||dS)N)ratio allow_holes)r; concave_hullrH)rYrrrrrrszGeometryArray.concave_hullcCstt|j|jdSrk)rr; convex_hullrHr'r]rrrrszGeometryArray.convex_hullcCstt|j|||jdSrk)rr;delaunay_trianglesrHr')rY toleranceZ only_edgesrrrr sz GeometryArray.delaunay_trianglescCstt|j|jdSrk)rr;enveloperHr'r]rrrrszGeometryArray.envelopecCstt|j|jdSrk)rr;minimum_rotated_rectanglerHr'r]rrrrsz'GeometryArray.minimum_rotated_rectanglecCstt|j|jdSrk)rr;exteriorrHr'r]rrrrszGeometryArray.exteriorcCstt|j|jdSrk)rr;extract_unique_pointsrHr'r]rrrrsz#GeometryArray.extract_unique_pointsround@cCs ttj|j||||d|jdS)N) quad_segs join_style mitre_limitr&)rr; offset_curverHr')rYdistancerrrrrrr"szGeometryArray.offset_curvecCs t|jSr)r; interiorsrHr]rrrr.szGeometryArray.interiorscCsttj|j|d|jdS)Nrr&)rr;remove_repeated_pointsrHr')rYrrrrr3sz$GeometryArray.remove_repeated_pointscCstt|j|jdSrk)rr;representative_pointrHr'r]rrrr9sz"GeometryArray.representative_pointcCstt|j|jdSrk)rr;minimum_bounding_circlerHr'r]rrrr<sz%GeometryArray.minimum_bounding_circlecCs t|jSr)r;minimum_bounding_radiusrHr]rrrrAsz%GeometryArray.minimum_bounding_radiuscCstt|j|jdSrk)rr; normalizerHr'r]rrrrDszGeometryArray.normalizecCstt|j|jdSrk)rr; make_validrHr'r]rrrrGszGeometryArray.make_validcCstt|j|jdSrk)rr;reverserHr'r]rrrrJszGeometryArray.reversecCstt|j||jdSrk)rr; segmentizerHr')rYZmax_segment_lengthrrrrMs zGeometryArray.segmentizecKsnt|trTt|t|kr6dt|t|}t|t||sNt||dd|j}tt ||j|fi|S)Nz5Lengths of inputs do not match. Left: {0}, Right: {1}r/) rrr2rrFr+r7rHgetattrr;)opr(r)rMmsgrrr_binary_methodWs   zGeometryArray._binary_methodcCs|d||S)NcoversrrYotherrrrreszGeometryArray.coverscCs|d||S)N covered_byrrrrrrhszGeometryArray.covered_bycCs|d||S)NcontainsrrrrrrkszGeometryArray.containscCs|d||S)NcrossesrrrrrrnszGeometryArray.crossescCs|d||S)NdisjointrrrrrrqszGeometryArray.disjointcCs|d||S)Nequalsrrrrr geom_equalstszGeometryArray.geom_equalscCs|d||S)N intersectsrrrrrrwszGeometryArray.intersectscCs|d||S)NoverlapsrrrrrrzszGeometryArray.overlapscCs|d||S)Ntouchesrrrrrr}szGeometryArray.touchescCs|d||S)NwithinrrrrrrszGeometryArray.withincCs|jd|||dS)NZ equals_exactrr)rYrrrrrgeom_equals_exactszGeometryArray.geom_equals_exactcCs&tjdtdd||dd| S)NzThe 'geom_almost_equals()' method is deprecated because the name is confusing. The 'geom_equals_exact()' method should be used instead.r[r/g? )r4r5 FutureWarningr)rYrdecimalrrrgeom_almost_equalss z GeometryArray.geom_almost_equalscCstt|j|||||jdSrk)rr; clip_by_rectrHr')rYxminyminxmaxymaxrrrrszGeometryArray.clip_by_rectcCst|d|||jdS)N differencer&rrr'rrrrrszGeometryArray.differencecCst|d|||jdS)N intersectionr&rrrrrrszGeometryArray.intersectioncCst|d|||jdS)Nsymmetric_differencer&rrrrrrsz"GeometryArray.symmetric_differencecCst|d|||jdS)Nunionr&rrrrrrszGeometryArray.unioncCst|d|||jdS)N shortest_liner&rrrrrrszGeometryArray.shortest_linecCs|jdd|d||S)Nr/rrfrrrrrrs zGeometryArray.distancecKs"|jdd|jd||fi|S)Nrr/hausdorff_distancerrYrrMrrrrs z GeometryArray.hausdorff_distancecKs"|jdd|jd||fi|S)Nrr/frechet_distancerrrrrrs zGeometryArray.frechet_distancecKsFt|ttfr|dks"|jddttj|j|fd|i||jdS)Nrrr/ resolutionr&) rintfloatrfrr;bufferrHr')rYrrrMrrrrs  zGeometryArray.bufferFcCs(|jddttj|j||d|jdS)Nrr/ normalizedr&)rfrr; interpolaterHr')rYrrrrrrs  zGeometryArray.interpolateTcCsttj|j||d|jdS)N)preserve_topologyr&)rr;simplifyrHr')rYrrrrrrs zGeometryArray.simplifycCs6t|trt|}nt|tr$|j}tj|j||dS)Nr)rrr?rrHr;project)rYrrrrrrs    zGeometryArray.projectcCst|tr|j}t|j|Sr)rrrHr;relaterrrrrs zGeometryArray.relatecCs t|jSr)r; unary_unionrHr]rrrrszGeometryArray.unary_unioncCsttd|j||jdS)Naffine_transformr&rr;Z_affinity_methodrHr')rYmatrixrrrrszGeometryArray.affine_transformcCsttd|j||||jdS)N translater&r)rYxoffyoffZzoffrrrrszGeometryArray.translatecentercCs ttjd|j|||d|jdS)Nrotateorigin use_radiansr&r)rYZanglerrrrrrs  zGeometryArray.rotate?c Cs"ttjd|j||||d|jdS)Nscale)rr&r)rYZxfactZyfactZzfactrrrrrs zGeometryArray.scalec Cs"ttjd|j||||d|jdS)Nskewrr&r)rYZxsZysrrrrrr s zGeometryArray.skewcCs~|jdurtd|dur&t|}n|dur:t|}ntd|j|rR|St|j|dd}t|j |j}t ||dS)u Returns a ``GeometryArray`` with all geometries transformed to a new 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 NzICannot transform naive geometries. Please set a crs on the object first.zMust pass either crs or epsg.TZ always_xyr&) r'rFrr` from_epsgZ is_exact_sameTransformerFromCRSr;Z transformrHr)rYr'Zepsg transformerZnew_datarrrto_crssE    zGeometryArray.to_crsWGS 84c Cs|jstd|j\}}}}|jjrBt||g}t||g}ntt|jddd}|||||\}}}}t||g}||kr|d7}t||g}|ddd}nt||g}t|t ||||dd} zt | d j WSt ytd Yn0d S) uvReturns 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 z$crs must be set to estimate UTM CRS.z EPSG:4326Trih)Zwest_lon_degreeZsouth_lat_degreeZeast_lon_degreeZnorth_lat_degree) datum_nameZarea_of_interestrzUnable to determine UTM CRSN)r' RuntimeError total_boundsrbr"meanrZtransform_boundsr r rrcode IndexError) rYrZminxZminyZmaxxZmaxyZx_centerZy_centerrZ utm_crs_listrrrestimate_utm_crsks:'   zGeometryArray.estimate_utm_crscCst|j|dkrd|j}|rV|}tj|ttjd}t |j |||<|St |j Sn d}t |dS)z8Return the x location of point geometries in a GeoSeriesPointrh fill_valuez5x attribute access only provided for Point geometriesN) rrzallranyr" full_likerr#r;Zget_xrHrFrYr|ZnonemptyZcoordsmessagerrrrQszGeometryArray.xcCst|j|dkrd|j}|rV|}tj|ttjd}t |j |||<|St |j Sn d}t |dS)z8Return the y location of point geometries in a GeoSeriesrrz5y attribute access only provided for Point geometriesN) rrzrrr r"r rr#r;Zget_yrHrFr rrrrRszGeometryArray.ycCst|j|dkrd|j}|rV|}tj|ttjd}t |j |||<|St |j Sn d}t |dS)z8Return the z location of point geometries in a GeoSeriesrrz5z attribute access only provided for Point geometriesN) rrzrrr r"r rr#r;Zget_zrHrFr rrrrSszGeometryArray.zcCs t|jSr)r;boundsrHr]rrrr szGeometryArray.boundsc Cst|dkr&ttjtjtjtjgS|j}tt|dddft|dddft|dddft|dddffS)Nrr r[r,)r2r"arrayr#r ZnanminZnanmax)rYbrrrrs zGeometryArray.total_boundscCs|jjSr)rHsizer]rrrrszGeometryArray.sizecCs|jfSr)rr]rrrri szGeometryArray.shapecCs t|jSr)r2rir]rrrrVszGeometryArray.ndimcOst|j|jdSrk)rrHcopyrW)rYargsrMrrrrszGeometryArray.copycCsddlm}|rL|dus"t|r(d}n$t|tr|||S|S)au 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 If method is specified, this is the maximum number of consecutive NaN values to forward/backward fill. In other words, if there is a gap with more than this number of consecutive NaNs, it will only be partially filled. If method is not specified, this is the maximum number of entries along the entire axis where NaNs 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 supported)NotImplementedErrorrzrr r)rYramethodlimitrmask new_valuesrrrfillna?s  zGeometryArray.fillnacCst|tr|r|S|Snjtjj|rvtjj|svt|}tjj |}t|tj rhtj ||dS|j |ddSt j |||dSdS)a 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. rtF)r)rhrN)rrrrnrotypesZis_string_dtypeZis_object_dtyperOZ pandas_dtypeZ StringDtyperastyper")rYrhrZ string_valuesZpd_dtyperrrr is  zGeometryArray.astypecCsBtjrt|jStjr$t|jStjdd|jDddSdS)zI Boolean NumPy array indicating if each value is missing cSsg|] }|duqSrr).0grrr z&GeometryArray.isna..boolrtN) r8r9rZ is_missingrHr:r@r"rr]rrrrzs   zGeometryArray.isnadropnacCsLt|}ddlm}m}||j|d}|tt|j d|j|_|S)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)Indexrur&N) rLpandasr(ru value_countsrJr"whereindexrz)rYr'rvr(rurrrrr*szGeometryArray.value_countscCsddlm}||\}}|S)zwCompute the ExtensionArray of unique values. Returns ------- uniques : ExtensionArray r) factorize)r)r-)rYr-_Zuniquesrrruniques  zGeometryArray.uniquecCs|jjSr)rHnbytesr]rrrr0szGeometryArray.nbytesr cst||}|j|_|S)aL 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``. )supershiftr')rYZperiodsrZshiftedrTrrr2szGeometryArray.shiftcCst|tr|g}t|S)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 )rrrB)rZscalarsrhrrrr_from_sequences zGeometryArray._from_sequencecCst|}|dfS)asReturn 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)rL)rYvalsrrr_values_for_factorizesz#GeometryArray._values_for_factorizecCst||jdS)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 r&)rJr')rrvoriginalrrr_from_factorizedszGeometryArray._from_factorizedc Csddlm}|jdkr$tjgddS|j}|}||B}|r|}t |d}|jrr||d}nt j dd}|||<n|}|r|j }|d|dd} |d|dd} |t| | | | g7}nd }|||d } |rd| |<| S) aReturn values for sorting. Returns ------- ndarray The transformed values should maintain the ordering between values within the array. See Also -------- ExtensionArray.argsort r)_hilbert_distanceZuint32rtr[g{Gz?r,r N)r)Zgeopandas.tools.hilbert_curver9rr"rrr rzrZnonzerorrrr) rYr9Z mask_emptyZ has_emptyrrIrr=rrrZ distancesrrr_values_for_argsort-s0     z!GeometryArray._values_for_argsort)skipnareturncCs tddSNz%geometries have no minimum or maximumrrYr;rrrargminhszGeometryArray.argmincCs tddSr=r>r?rrrargmaxkszGeometryArray.argmaxcs|rddl}|jjdur|jr8|jjr2dqdnt>tjdtd|| ddj \}}}}Wdn1s0Yd|krd krnnNd|krd krnn6d |krd krnnd |krd krnndndfd d St S)aOFormatting 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``. rNr,rignore)categoryriLriZcstjj|dS)N)Zrounding_precision)rZwktdumpsr< precisionrrr$z*GeometryArray._formatter..) geopandasoptionsZdisplay_precisionr'Z is_projectedr4catch_warnings simplefilterRuntimeWarningrzrrepr)rYZboxedrIrrrrrrFr _formatterns> >     zGeometryArray._formattercCs$tdd|D}t|t|dS)z Concatenate multiple array Parameters ---------- to_concat : sequence of this type Returns ------- ExtensionArray cSsg|] }|jqSr)rHr!garrrr#r$z3GeometryArray._concat_same_type..r&)r"Z concatenater_get_common_crs)rZ to_concatrDrrr_concat_same_types zGeometryArray._concat_same_typecKsH|dks|dkr tt||Stdt|jd|jd|ddS)Nr r'z ' with dtype z does not support reduction ')rrGrrrrh)rYrr;rMrrr_reduceszGeometryArray._reducecCst|S)zj The numpy array interface. Returns ------- values : numpy array )rG)rYrhrrr __array__szGeometryArray.__array__cstfdd}t|tjtjtjfr&tS}||}t|t|krJtdfddt||D}t j |t d}|S)Ncs8t|s&t|ts tjj|r&|}n|gt}|Sr)rArrrnrorZ is_list_liker2)paramZovaluesr]rrconvert_valuess z,GeometryArray._binop..convert_valueszLengths must match to comparecsg|]\}}||qSrr)r!ar)rrrr#r$z(GeometryArray._binop..rt) rrnrur(rwNotImplementedr2rFzipr"Zasarrayr%)rYrrrXZlvaluesZrvaluesresr)rrYr_binops zGeometryArray._binopcCs||tjSr)r]operatoreqrrrr__eq__szGeometryArray.__eq__cCs||tjSr)r]r^nerrrr__ne__szGeometryArray.__ne__cCsHt|r<||jjus,t||jjs,|dur8|SdS||kS)z, Return for `item in self`. NF)r;rzrhr$rrr )rYitemrrr __contains__s    zGeometryArray.__contains__)N)rrr)r)r)F)T)F)rrr)rF)rrrr)rrrF)NN)r)FN)NNNT)T)T)r N)NF)T)T)F)T)N)xrr r!__doc__rrgrZpropertyrDr^r_r'setterrfrhrjrsrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrr staticmethodrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrQrRrSr rrrirVrrrrr rzr%r*r/r0r2r%r4r6r8r:rr@rArOrSrUrVr]r`rbrd __classcell__rrr3rr s&         .                               Z S          * #    '  ; 7  rcCsdd|D}dd|D}dd|D}dd|D}t|dkrHdSt|dkrt|dkrztjd |dd d d |dStd |ddS)NcSs&g|]}|r|jdus|qSr)rzrr'rPrrrr#r$z#_get_common_crs..cSsh|] }|jqSrr&)r!Zarrrrr r$z"_get_common_crs..cSsg|]}|dur|qSrrr!r'rrrr#r$cSsg|] }|jqSr)rrkrrrr#r$rr zJCRS not set for some of the concatenation inputs. Setting output's CRS as z$ (the single non-null crs provided).r[r/z:Cannot determine common CRS for concatenation inputs, got zH. Use `to_crs()` to transform geometries to the same CRS before merging.)r2r4r5rF)Zarr_seqZcrs_setZ crs_not_nonenamesrrrrRs$    rR)F)r,)N)N)F)N)NN)7rlr^r4rc functoolsrnumpyr"r)rnrrrrrZshapely.affinityZshapely.geometryZshapely.geometry.baserZ shapely.opsZ shapely.wktZpyprojrrZ pyproj.aoir Zpyproj.databaser r@ ImportErrorZgeosr r8r r;r^rZfrom_crsrrr+r7r>r?rArBrGrJrLrNrOrPrrRrrrrsb             #      .w