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.
--
https://mail.python.org/mailman/listinfo/python-list
