ArcMap中使用ArcPy InsertCursor插入的GPS点位置偏移，坐标系核对无误仍存在异常

我太懂你这种反复核对坐标系却还是出问题的抓狂感了！你本来想每月用Python脚本把野外采集的GPS井位更新到ArcMap的shapefile里，属性和新点的相对位置都对，但一更新现有文件，新点就跑到完全离谱的地方——哪怕用和已有点一模一样的坐标也不行。你已经翻来覆去检查了shapefile、边界图层和数据框的坐标系（GCS_WGS_1984和WGS_1984_UTM_Zone_16N），但问题还是没解决，核心疑问应该就是：Python里的空间参考是不是没定义对？

其实你代码里有两个关键问题，直接导致了点的偏移：

问题1：空间参考只创建了对象，却没真正用到 #

你写了arcpy.SpatialReference(32616)，但没有把这个对象赋值给变量，也没有关联到要插入的几何上，这行代码等于白写，完全没起到作用。

问题2：坐标顺序和坐标系不匹配，直接插错了单位 #

你的GPS是十进制经纬度（纬度43.20779，经度-89.78513），对应GCS_WGS_1984（WKID4326，度单位），但目标shapefile用的是UTM16N（WKID32616，米单位）。你之前直接把(lat, lon)传给SHAPE@XY，相当于把纬度当X（经度）、经度当Y（纬度），还把度数值当成米来插入，这肯定会让点跑到完全错误的位置！

修正后的代码及说明 #

下面是调整后的脚本，我会标注关键细节：

import arcpy

# 测试用的GPS坐标（格式：(Well_ID, (纬度, 经度))）
coords = [(1, (43.20779, -89.78513))]

# 目标shapefile的路径
shapefile = r'M:\ARC\gis_users\Missing_Wells\Missing_Wells.shp'

# 明确定义两个空间参考：输入的地理坐标系（GPS用的WGS84）和目标的投影坐标系（UTM16N）
input_sr = arcpy.SpatialReference(4326)
target_sr = arcpy.SpatialReference(32616)

# 使用SHAPE@而不是SHAPE@XY，这样可以传入带空间参考的几何对象
with arcpy.da.InsertCursor(shapefile, ['Well_ID', 'Y_Coord', 'X_Coord', 'SHAPE@']) as cursor:
    for Well_ID, (lat, lon) in coords:
        # 重点：ArcPy的Point是(x, y)，对应地理坐标系的「经度，纬度」，必须调换顺序！
        point = arcpy.Point(lon, lat)
        # 先创建带输入坐标系的PointGeometry，再投影到目标shapefile的坐标系
        projected_point = arcpy.PointGeometry(point, input_sr).projectAs(target_sr)
        # 插入行：保留原始经纬度到属性字段，同时插入投影后的几何
        cursor.insertRow((Well_ID, lat, lon, projected_point))
        print(f'已成功添加Well_{Well_ID}')

额外验证小技巧 #

你可以打开目标shapefile的属性表，看看已有点的X/Y坐标值——应该是UTM的米数（比如X大概在50万左右，Y大概在478万左右）。而你之前插入的是经纬度数值（X=-89，Y=43），数值差了好几个数量级，肯定会偏移到地图视野外，你可以试试缩小地图到最小，说不定能找到这些“失踪”的点。

备注：内容来源于stack exchange，提问作者Publicchee