Python 将utm坐标转换为参照相邻分区的坐标

我有成千上万的NAD83 UTM 13N光栅。我试图使用arcpy.GetCellValue_管理（raster.tif，point）按点提取数据，但数据的最西侧位于UTM 12N区域。有没有办法从12N获得坐标，但参考13N？项目要求所有数据均为UTM 13N格式，即使是全州项目。我知道这很愚蠢

这可以通过GDAL实现。将dataPoints.shp保存在所需的UTM网格（即UTM 13N）中，然后使用GDAL加载点图层、获取字段、获取几何图形、获取边界坐标、geotransform、光栅带、点坐标（UTM 13N）并将光栅读取为阵列。在所有光栅上构建一个循环，它工作得非常快。感谢您提供的详细信息

from osgeo import gdal, ogr
shp_filename = 'C:\\Path\\dataPoints_UTM13.shp'
ds = ogr.Open(shp_filename)
lyr = ds.GetLayer()
for feat in lyr:
    point_id_obj = feat.GetField("Sample")
    name = feat.GetField("Location_D")
    geom = feat.GetGeometryRef()
    mx, my = geom.GetX(), geom.GetY()
    path = 'C:\\RasterPath'
    raster = 'myraster'
    ras_open = gdal.Open('{a}\\{b}.tif'.format(a=path, b=raster))
    gt = aws_open.GetGeoTransform()
    rb = aws_open.GetRasterBand(1)
    px = abs(int((mx - gt[0]) / gt[1]))
    py = int((my - gt[3]) / gt[5])
    ras_obj = rb.ReadAsArray(px, py, 1, 1)
    print point_id_obj
    print name
    print mx, my