Más

Convierta un mercator web geotiff en una proyección equidistante azimutal

Convierta un mercator web geotiff en una proyección equidistante azimutal


Tengo un archivo tiff georreferenciado, en la proyección web-mercator; y lo necesito convertido en una proyección equidistante azimutal. Soy nuevo en SIG y no sé 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 sistemas operativos específicos, tengo computadoras con: Linux-Debian, Mac OX 10.6.8, Windows 7 y Windows XP.


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. Puede obtener la extensión correcta en coordenadas de pseudo mercator en https://github.com/mapnik/mapnik/wiki/XMLConfigReference:-20037508.34, -20037508.34, 20037508.34, 20037508.34. Con eso, puede establecer la extensión con gdal_translate manualmente antes de reproyectar.

La proyección de destino no está cubierta por EPSG, pero ESRI ha inventado un código 102016, que se incluye en GDAL y QGIS comoEPSG: 102016:

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

Obtendrá artefactos desagradables a lo largo del meridiano de 180 °, por lo que es una opción adicionalSOURCE_EXTRAes necesario. Quedará un punto ciego en el polo norte, porque el pseudo mercator 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 -Grst_20 -s -s 38 8192 -dos SOURCE_EXTRA = 1000 -of GTiff F: /Karten/merc_corrected.tif F: /Karten/world_aeqd.tif


  1. Cargue su capa ráster en QGIS. Haga clic derecho en el nombre de la capa ráster en la lista de capas y seleccione Guardar como… en el menú emergente. En el cuadro de diálogo Guardar como, seleccione un formato de destino, un nombre y otro CRS.

  2. También puede usar gdal_warp (GDAL Utilities) desde la línea de comando:

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

s_srs es la proyección de origen, t_srs es la proyección de destino.


Ver el vídeo: Tutorial de Representación Proyecciones en Mapas con QGIS