Python shp文件的地理坐标转换

Python 将shp文件的地理坐标转换为经纬度坐标

引言

我们在绘制地图的过程中，经常需要使用shp文件绘制某个区域的边界、轮廓等。然而，我们下载的某些shp文件的地理坐标不是经纬度坐标，导致直接读取shp文件绘图时出现问题，所以，我们需要将该文件的地理坐标转换为经纬度坐标。

安装库

pip install geopandas

完整示例

import geopandas as gpd

# 读取Shape文件
shapefile = gpd.read_file('../data_recon/borders/fr_data/fr.shp')

# 定义目标坐标转换器
target_crs = {'init': 'epsg:4326'}

# 坐标转换
shapefile.to_crs(target_crs, inplace=True)

# 保存转换后的Shape文件
shapefile.to_file('../data_recon/borders/fr_data/loess.shp')

这样就可以使用新的shp文件去绘图啦！或者可以不保存文件，在绘图代码中直接使用shapefile对象绘图。

参考资料

https://blog.51cto.com/u_16175489/8591450