geopandas bug ?
Peter Otten
__peter__ at web.de
Sat Jan 14 08:42:23 EST 2017
Xristos Xristoou wrote:
> i want to create a simple spatial joing using geopandas but i thing so
> geopandas has bug ?
Have you tried the examples on <http://geopandas.org/mergingdata.html>?
Do they work? If yes, inspect your data, does it have the same format?
> geopandas code :
>
> from geopandas import gpd
> import geopandas
> points = geopandas.GeoDataFrame.from_file('points.shp') # or geojson etc
> polys = geopandas.GeoDataFrame.from_file('polygons.shp')
> pointInPoly = gpd.sjoin(points, polys, how='left',op='within')
>
> error :
>
> Traceback (most recent call last):
> File "/home/sarantis/testshapely/sumpointsinsidepolygon/testgeo.py",
> line 7, in <module>
> pointInPoly = gpd.sjoin(points, polys, how='left',op='within')
> File "/usr/local/lib/python2.7/dist-packages/geopandas/tools/sjoin.py",
> line 57, in sjoin
> r_idx = np.concatenate(idxmatch.values)
> ValueError: need at least one array to concatenate
>
> and if i change the imports with the some code :
>
> import geopandas
> import pandas as pd
> import geopandas as gpd
> from geopandas import GeoDataFrame, read_file
> from geopandas.tools import sjoin
> from shapely.geometry import Point, mapping,shape
> import pandas as gpd
>
> i take that error :
>
> pointInPoly = gpd.sjoin(points, polys, how='left',op='within')
> AttributeError: 'module' object has no attribute 'sjoin'
>
>
> any idea why ?
My crystal ball says that either points or polys is empty ;)
More information about the Python-list
mailing list