多个geojson经过坐标系转换后合并为一个shp数据

ylfmsn

已于 2023-06-15 13:43:59 修改

阅读量701

点赞数

分类专栏： python 文章标签： python pandas

于 2023-06-15 10:18:51 首次发布

本文链接：https://blog.csdn.net/ylfmsn/article/details/131222196

版权

python 专栏收录该内容

12 篇文章 3 订阅

订阅专栏

该代码示例展示了如何使用geopandas和pyproj库将多个WGS84坐标系统的GeoJSON文件转换到EPSG2363西安80高斯投影，并将它们合并成一个单一的Shapefile（shp）数据。首先，创建一个空的GeoDataFrame，然后逐个读取GeoJSON文件，进行坐标转换，并添加文件名作为特征。转换后的GeoDataFrames被合并，过滤掉面积小于100的要素，最后将结果保存为Shapefile。

摘要由CSDN通过智能技术生成

利用geopandas和pyproj将多个geojson数据经过坐标转换后合并成一个shp数据：

import geopandas as gpd
from pyproj import CRS
import pandas as pd


BASE_DIR = 'D:/demo'

# 输入的WGS84 GeoJSON文件列表
geojson_files = [f"{BASE_DIR}/file1.geojson", f"{BASE_DIR}/file2.geojson", f"{BASE_DIR}/file3.geojson"]

# 目标坐标系的EPSG代码，使用一个西安80高斯投影举例
target_epsg = 2363

# 创建一个空的GeoDataFrame
merged_gdf = gpd.GeoDataFrame()

# 定义WGS84和目标坐标系之间的坐标转换
wgs84_crs = CRS.from_epsg(4326)
target_crs = CRS.from_epsg(target_epsg)

# 逐个读取、转换和合并GeoJSON文件
for geojson_file in geojson_files:
    gdf = gpd.read_file(geojson_file, crs=wgs84_crs)
    # 为每个geodataframe新添加字符类型的字段name，并将文件名赋值给字段name
    gdf['name'] = geojson_file[-13:-8]
    gdf = gdf.to_crs(target_crs)
    merged_gdf = pd.concat([merged_gdf, gdf], ignore_index=True)

# 过滤掉面积小于100的要素
gdf_filtered = merged_gdf[merged_gdf["area"] >= 100]
# 保存为shp文件
gdf_filtered.to_file(output_shp, driver='ESRI Shapefile', encoding='utf-8')

geopandas版本：0.12.0，pandas版本：1.5.3，pyproj版本：3.4.0