Quiero convertir archivos de datos de cuadrícula ascii de Deutscher Wetterdienst a GeoTIFF utilizando GDAL en Python, pero tengo problemas para añadir el sistema de coordenadas correcto (ya que DWD utiliza su propia cuadrícula RADOLAN es un poco problemático).
Estoy utilizando CreateCopy() que por lo que he leído utiliza el comando gdal_translate. En gdal_translate simplemente añadiría los valores correctos con appending gdal_translate -a_srs="+proj=stere +lat_0=90.0 +lon_0=10.0 +lat_ts=60.0 +a=6370040 +b=6370040 +units=m"
pero no encuentro la forma de hacerlo en la API de Python.
¿Hay alguna forma de añadir estos valores mientras se traduce, o de añadirlos después de crear el archivo utilizando el conjunto de datos devuelto por .CreateCopy()?