Τη Σάββατο, 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 ? -- https://mail.python.org/mailman/listinfo/python-list
