ToB企服应用市场:ToB评测及商务社交产业平台

标题: geopandas空间连接 [打印本页]

作者: 刘俊凯    时间: 2022-12-2 02:59
标题: geopandas空间连接

 
 
 
import geopandas as gpb
Berlin = gpb.read_file(r'E:\项目\市.shp')
Berlin.head()
a = gpb.read_file(r'G:\空间单元WGS84.shp')
a.head()
 
b =gpb.sjoin(left_df=a[['grid_id','geometry']], right_df=Berlin[['fnid','geometry']], op='intersects',how = 'right')
b.dropna(axis=0,how='any',inplace = True)
len(b)
 
e = b[['grid_id','fnid']]
e.to_csv(r'E:\项目\20211013daas平台\01数据\20220621联通数据改版\连接后shp\upload_grid_geo2022.csv')
#data.plot()
#针对一对多的情况
 

免责声明:如果侵犯了您的权益,请联系站长,我们会及时删除侵权内容,谢谢合作!




欢迎光临 ToB企服应用市场:ToB评测及商务社交产业平台 (https://dis.qidao123.com/) Powered by Discuz! X3.4