import geopandas as gpd
import glob
file_name = "****.shp"
data = gpd.read_file(file_name, encoding="utf-8")
data.crs # {'init': 'epsg:4326'}
reprojeted_data = data.to_crs(crs="+proj =utm +zone=48 +ellps=WGS84 +datum=WGS84 +units=m +no_defs")
reprojected_data.to_file("result.shp", encoding="utf-8")
geopandas坐标转换
最新推荐文章于 2025-04-01 15:27:42 发布