坐标系转换（wgs84tomercator3857）

原创于 2025-07-23 10:24:04 发布 · 361 阅读

5 ·

CC 4.0 BY-SA版权

文章标签：

#python

# WGS84转墨卡托投影
import arcpy
import os

# 设置输入和输出文件夹路径
input_folder = r"F:\bianjie"  # 替换为实际路径
output_folder = r"F:\bianjie"  # 替换为实际路径

# 确保输出文件夹存在
if not os.path.exists(output_folder):
    os.makedirs(output_folder)

# 定义输出坐标系（Web Mercator，EPSG:3857）
output_coordinate_system = arcpy.SpatialReference(3857)  # 或使用名称："WGS 1984 Web Mercator (auxiliary sphere)"

# 遍历输入文件夹中的所有SHP文件
arcpy.env.workspace = input_folder
shp_files = arcpy.ListFeatureClasses("*.shp")  # 获取所有SHP文件

# 批量处理
for shp in shp_files:
    try:
        # 构造输入和输出路径
        input_shp = os.path.join(input_folder, shp)
        output_shp = os.path.join(output_folder, f"{os.path.splitext(shp)[0]}_mercator3857.shp")

        # 执行投影变换
        arcpy.Project_management(
            in_dataset=input_shp,
            out_dataset=output_shp,
            out_coor_system=output_coordinate_system,
            transform_method="",  # WGS84转Web Mercator无需地理变换
            preserve_shape="PRESERVE_SHAPE"
        )
        print(f"成功转换: {shp} → {os.path.basename(output_shp)}")

    except Exception as e:
        print(f"转换失败: {shp} | 错误: {str(e)}")