如何解决在Python中用geopandas重叠两个shapefile
我有这两个shapefile。 [第一个链接上带有省的shapefile] [带有链接2上的区的第二个shapefile]
我需要加入/合并这两个shapefile以返回如下图: [![莫桑比克各区如下图所示]
moz_admin = 1:https://drive.google.com/file/d/1MYNkuKFXP0dt76G9OgeBotjF5UmiRqEl/view?usp=sharing
moz_admin_district = [2]:https://drive.google.com/file/d/1idf5VKgN8PZgdoAcBVEcY9pDa1WYpg7-/view?usp=sharing
到目前为止我所做的:
import os
import geopandas as gpd
file = os.listdir(r'pathtofilefolder')
path = [os.path.join(r'folder',i) for i in file if ".shp" in i]
gdf = gpd.GeoDataFrame(pd.concat([gpd.read_file(i) for i in path],ignore_index=True),crs=gpd.read_file(path[0]).crs)
和
moz_admin.to_crs(moz_admin_district.crs,inplace=True) #Change projections
gpd.sjoin(moz_admin_district,moz_admin,how='union',op='within').shape #Join 2 shapefiles
gpd.sjoin(moz_admin_district,op='within').plot()
和
diffs = []
gdfs = [moz_admin,moz_admin_district]
for idx,gdf in enumerate(gdfs):
if idx < 2:
diffs.append(gdf.symmetric_difference(gdfs[idx+1]).iloc[0])
diffs.append(moz_admin_district.iloc[0].geometry)
和
gdf = gpd.GeoDataFrame(pd.concat([moz_admin,moz_admin_district]))
和
from shapely.geometry import Polygon
res_union = gpd.overlay(moz_admin,moz_admin_district,how='union')
和
so = moz_admin.merge(moz_admin_district,on='ADM1_PT',how='inner')
df = tab_df.merge(spatial_df,on='mukey',how='right')
sof = gpd.GeoDataFrame(so)
和
merged_master = gpd.GeoDataFrame(pd.merge(moz_admin,how='left',left_on="ADM1_PT",right_on="ADM1_PT"))
merged = merged_master[['ADM1_PT','ADM2_PT','geometry']]
和
so = gpd.GeoDataFrame(pd.concat([moz_admin,moz_admin_district],ignore_index=True))
和
sjoined_listings = gpd.sjoin(moz_admin_district,op='intersects')
和
merged=sjoin(moz_admin_district,op='intersects')
没有成功!
解决方法
#将两个形状文件重叠为一个 ax = moz_admin.plot(颜色='none',edgecolor ='black',线宽= 1.5) moz_admin_district.plot(ax = ax,color ='none',edgecolor ='grey',线宽= 0.5)
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。