geopandas bug ?
Xristos Xristoou
saxri89 at gmail.com
Sat Jan 14 09:06:15 EST 2017
Τη Σάββατο, 14 Ιανουαρίου 2017 - 3:43:10 μ.μ. UTC+2, ο χρήστης Peter Otten έγραψε:
> 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 ;)
is not empty and yes i have two shapefiles from qgis.what is the error ?
More information about the Python-list
mailing list