Convertir un geotiff web-mercator en una proyección azimutal equidistante

Question

Tengo un archivo tiff georeferenciado, en la proyección web-mercator; y necesito convertirlo a una proyección azimutal equidistante. Soy nuevo en el GIS y no puedo averiguar cómo hacer esto.

Estoy abierto a usar cualquier programa para cualquier sistema operativo, siempre que haga el trabajo.

Y en lo que respecta a los sistemas operativos específicos, tengo computadoras con..: Linux-Debian, Mac OX 10.6.8, Windows 7 y Windows XP.

Answer 1

Su fuente tif está en Pseudo mercator, pero la extensión almacenada dentro del archivo está en grados. Esto no puede ser interpretado correctamente por gdalwarp, por lo que entrega el archivo fuente intacto. Puedes obtener la extensión correcta en las coordenadas de Pseudo Mercator de https://github.com/mapnik/mapnik/wiki/XMLConfigReference : -20037508.34, -20037508.34, 20037508.34, 20037508.34 . Con eso, puedes establecer la extensión con gdal_traducir manualmente antes de reproyectar.

La proyección del objetivo no está cubierta por el EPSG, pero ESRI ha inventado un código 102016, que está incluido en GDAL y QGIS como EPSG:102016 :

+proj=aeqd +lat_0=90 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs

Recibirás desagradables artefactos a lo largo del meridiano de 180°, así que una opción extra SOURCE_EXTRA es necesario. Un punto ciego en el polo norte permanecerá, porque el pseudo-mercader no lo cubre:

gdal_translate -a_srs EPSG:3857 -of GTiff -a_ullr -20037508.34 20037508.34 20037508.34 -20037508.34 F:/Download/63575605298.3c.tif F:/Karten/merc_corrected.tif
gdalwarp -overwrite -s_srs EPSG:3857 -t_srs EPSG:102016 -ts 8192 8192 -wo SOURCE_EXTRA=1000 -of GTiff F:/Karten/merc_corrected.tif F:/Karten/world_aeqd.tif

Answer 2

1. Cargue su capa de trama en el QGIS. Haga clic con el botón derecho del ratón en el nombre de la capa rasterizada en la lista de capas y seleccione Salve a As... del menú emergente. En el cuadro de diálogo Guardar como, seleccione un formato de destino, un nombre y otro CRS.

2. También puedes usar gdal_warp (Utilidades GDAL) desde la línea de comandos:

gdalwarp -s_srs +init=epsg: yyyy -t_srs +init=epsg: xxxx input_raster output_raster

s_srs es la proyección fuente, t_srs es la proyección objetivo.