Τη Σάββατο, 14 Ιανουαρίου 2017 - 4:30:48 μ.μ. UTC+2, ο χρήστης Peter Otten έγραψε: > Xristos Xristoou wrote: > > > Τη Σάββατο, 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? > > Looks like you chose to ignore the hard part. > > >> > 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 > > >> > 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. > > Can I download those files somewhere? > > > what is the error ? > > No idea. Without those files I cannot run your code.
yes i sent you -- https://mail.python.org/mailman/listinfo/python-list