Soy nuevo en gdal y tengo una simple pregunta. Quiero transformar una imagen de la proyección OpenStreetMap a la proyección Nasa Blue Marble. Lo que entiendo es que la proyección OpenStreetMap es EPSG:3857 y la proyección Blue Marble es EPSG:4326 (http://www.paulillsley.com/gia/index.html).
Así que simplemente descargué un azulejo de muestra de 256*256 como este: http://a.tah.openstreetmap.org/Tiles/tile/0/0/0.png y aplicar el siguiente comando:
gdalwarp -s_srs EPSG:3857 -t_srs EPSG:4326 "0.png" "o_bluemarble.png"
He recibido el siguiente error: "Copiando la tabla de colores de 0.png al nuevo archivo. ERROR 1: No se ha podido calcular una transformación entre coordenadas de píxeles/líneas y coordenadas georeferenciadas para 0.png. No hay transformación afín ni GCPs".
¿Por dónde debo empezar?
Actualización 1
Lo tengo. Necesito definir los límites. En primer lugar, tengo los límites:
gdaltransform -s_srs EPSG:4326 -t_srs EPSG:3857
-180 85.05
-20037508.3427892 20036051.9193368 0Utilicé 85,05 ya que es la aproximación al infinito en Spherical Mercator que utiliza Bing. Luego, utilicé los siguientes comandos:
./gdal_translate -of Gtiff -co "tfw=yes" -a_ullr -20037508.3427892 20036051.9193368 20037508.3427892 -20036051.9193368 -a_srs "EPSG:3857" "input.png" "input_tfw.tiff"
./gdalwarp -s_srs EPSG:3857 -t_srs EPSG:4326 -ts 256 128 "input_tfw.tiff" "output.tif"El problema ahora es que al comparar una imagen descargada de openstreetmap con Blue Marble, los polos no están bien alineados. Me refiero aquí a Groenlandia, no al polo lejano. ¿Cuál crees que es el problema?
