像素/数组位置转换为纬度经度gdal Python

Question

我正在尝试将表示.tif的光栅中的位置转换为相应的全局坐标。将整个阵列转换为tif并将其加载到QGIS中，一切都可以参考，但是使用下面的单点计算方法，会有一个轻微的偏移（在结果坐标中为东北偏东。。。。在

在光栅.tif使用ETRS 89 UTM区域32N

有人有主意吗？在

from osgeo import ogr, gdal, osr
import numpy as np


raster = gdal.Open("rasters/raster.tif")
raster_array = np.array(raster.ReadAsArray())


def pixel2coord(x, y):
     xoff, a, b, yoff, d, e = raster.GetGeoTransform()
     xp = a * x + b * y + xoff
     yp = d * x + e * y + yoff
     return(xp, yp)


print(pixel2cood(500,598))

Answer 1

我认为问题可能是xoff和yoff包含最左上像素的左上角的坐标，您需要计算像素中心的坐标。在

def pixel2coord(x, y):
    xoff, a, b, yoff, d, e = raster.GetGeoTransform()

    xp = a * x + b * y + a * 0.5 + b * 0.5 + xoff
    yp = d * x + e * y + d * 0.5 + e * 0.5 + yoff
    return(xp, yp)