मैं लाखों जीपीएस बिंदुओं के लिए एक देश (और कभी-कभी राज्य) खोजने के लिए नीचे दिए गए कोड का उपयोग कर रहा हूं। वर्तमान में कोड प्रति सेकंड एक सेकंड लगता है, जो अविश्वसनीय रूप से धीमा है। शेपफाइल 6 एमबी है।
मैंने पढ़ा कि जियोपैन्ड्स स्थानिक जुड़ावों के लिए rtrees का उपयोग करता है, जिससे वे अविश्वसनीय रूप से कुशल हो जाते हैं, लेकिन यह यहां काम नहीं करता है। मैं क्या गलत कर रहा हूं? मैं प्रति सेकंड एक हजार अंक की उम्मीद कर रहा था।
शेपफाइल और सीएसवी को यहां डाउनलोड किया जा सकता है (5MB): https://www.dropbox.com/s/gdkxtpqupj0sidm/SpatialJoin.zip?dl=0
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 time
#parameters
shapefile="K:/.../Shapefiles/Used/World.shp"
df=pd.read_csv("K:/.../output2.csv",index_col=None,nrows=20)# Limit to 20 rows for testing
if __name__=="__main__":
start=time.time()
df['geometry'] = df.apply(lambda z: Point(z.Longitude, z.Latitude), axis=1)
PointsGeodataframe = gpd.GeoDataFrame(df)
PolygonsGeodataframe = gpd.GeoDataFrame.from_file(shapefile)
PointsGeodataframe.crs = PolygonsGeodataframe.crs
print time.time()-start
merged=sjoin(PointsGeodataframe, PolygonsGeodataframe, how='left')
print time.time()-start
merged.to_csv("K:/01. Personal/04. Models/10. Location/output.csv",index=None)
print time.time()-start
आपका डाटा लिंक 404 है
—
हारून