¿Cómo convierto coordenadas afines a lat / lng?

Soy extremadamente nuevo en SIG.

Estoy usando gdalpara leer en un mapa de uso / cobertura del suelo y necesito seleccionar el lat / lng de ciertos tipos de cobertura del suelo para indexar en un conjunto de datos diferente que se expresa solo en lat / lng. Desafortunadamente, no entiendo la forma de las coordenadas X e Y que me ha dado de la transformación geográfica, en concreto el originXy originYa continuación:

geotransform = dataset.GetGeoTransform()
originX = geotransform[0]
originY = geotransform[3]

Imprimir estos valores me da coordenadas como (447466.693808, 4952570.40529). ¿Cómo se relacionan con la latitud y longitud originales?

Editar:

Aquí hay un ejemplo simple de Python que me dio lo que estaba buscando:

srs = osr.SpatialReference()
srs.ImportFromWkt(dataset.GetProjection())

srsLatLong = srs.CloneGeogCS()
ct = osr.CoordinateTransformation(srs,srsLatLong)
print ct.TransformPoint(originX,originY)

Robado de: tolatlong.py

gdal_translate reproyectará sus datos de cualquier proyección en cualquier otra (en este caso desea EPSG: 4326) usando:

gdal_translate -a_srs epsg:4326 srcfile new_file 

o podría usar gdaltrasform para convertir los puntos (y estoy seguro de que también puede acceder a eso desde Python (?))

La geotransformación se documenta en https://gdal.org/user/raster_data_model.html . La idea es que tome coordenadas (x, y) directamente del conjunto de datos, aplique una transformación lineal para obtener (u, v) con

u = a*x + b*y
v = c*x + d*y

(puede tomar esto como la definición de una transformación lineal), luego cambie el origen agregando geotransform [0] a u y geotransform [3] a v. Eso da la "transformación afín" de (x, y). Realmente tiene la intención de rotar, cambiar la escala, quizás corregir un poco para algunos errores de sesgo y reubicar las coordenadas específicas de los datos (x, y) para que coincidan con un sistema de coordenadas conocido. Se supone que el resultado produce coordenadas proyectadas. Esto simplemente significa que hay un procedimiento matemático que toma (longitud, latitud) y los convierte en las coordenadas conocidas: esto se llama "proyección". "Desproyectar" es hacer lo contrario; entonces, si sabes qué proyección es necesaria, la aplicas a las coordenadas transformadas afines (x, y) para obtener la latitud y longitud.

Por cierto, los valores de las constantes a, b, c, d están dados por las entradas 1, 2, 4 y 5 en la matriz de geotransformación.

Puedes usar lo siguiente:

import gdal, numpy as n
def coord(file):
    padfTransform = file.GetGeoTransform()
    indices = n.indices(file.ReadAsArray().shape)
    xp = padfTransform[0] + indices[1]*padfTransform[1] + indices[1]*padfTransform[2]   
    yp = padfTransform[3] + indices[0]*padfTransform[4] + indices[0]*padfTransform[5]  
    return xp,yp
file = gdal.Open('Yourfile.tif') # A GeoTiff file
x,y = coord(file)

coord devolverá la longitud (x) y la latitud (y) de todos los píxeles. Tenga en cuenta que las coordenadas son de la esquina izquierda de un píxel