Comprensión de las propiedades de creación de Map Tile

Question

Estoy intentando crear un mapa base de azulejos. Tengo una imagen de la que quiero crear el azulejo y estoy usando GDAL.

He conseguido crear las baldosas pero cuando veo el resultado parece que el mapa de baldosas está:

1. Se voltea
2. Las proporciones no son correctas

Esta es la imagen superpuesta del mapa:

Este es el resultado del alicatado:

He comprobado muchas veces que las coordenadas son buenas comparando las esquinas de la imagen en una ortofoto.

Otra cosa es que en google earth con las mismas coordenadas la imagen del mapa se ajusta con precisión al mapa de google.

EDITADO

gdalinfo map_overlay.png

con esta salida:

Conductor: PNG/Portable Network Graphics Archivos: map_overlay.png El tamaño es 3527, 2494 El sistema de coordenadas es `' Metadatos de la estructura de la imagen:
INTERLEAVE=PIXEL Coordenadas de esquina: Superior izquierda ( 0.0, 0.0) Superior izquierda ( 0,0, 2494,0) Superior derecha ( 3527,0, 0,0) Inferior derecha ( 3527,0, 2494,0) Derecha ( 3527,0, 2494,0) Centro ( 1763,5, 1247,0) Banda 1 Block=3527x1 Type=Byte, ColorInterp=Red Mask Flags: PER_DATASET ALPHA Band 2 Block=3527x1 Type=Byte, ColorInterp=Green Mask Flags: PER_DATASET ALPHA Band 3 Block=3527x1 Type=Byte, ColorInterp=Azul
Banderas de la máscara: PER_DATASET ALPHA Band 4 Block=3527x1 Type=Byte, ColorInterp=Alpha

Entonces traduce:

gdal_translate -de VRT -a_srs EPSG:4326 -gcp 0 0 31.7431761644 35.1680410195 -gcp 3527 0 31.7493769674 35.1784535489 -gcp 3527 2492 31.7431011291 35.1784951643 -gcp 0 2494 31.7431761644 35.1680410195 map_overlay.png map_overlay.vrt

con esta salida:

El tamaño del archivo de entrada es 3527, 2494

y por último:

gdal2tiles.py -z 14-21 map_overlay.vrt

La superposición es una simple imagen .png sin datos especiales.

Agradeceré cualquier ayuda.

Shani

Answer 1

Intenta utilizar la imagen coordinada con la métrica antes de alicatarla como:

   gdalwarp -of GTiff -t_srs EPSG:3857 input.tif output.tif 

---

ACTUALIZACIÓN

Nota: aunque especifiquemos los gcp's, gdal_translate no especificaría las coordenadas de las esquinas del tiff.

Además, no es necesario crear un conjunto de datos virtuales (.vrt) porque sólo hay un archivo ...

utilice el siguiente código:

  gdal_translate -of  GTiff  -gcp 0 0 31.7431761644 35.1680410195 -gcp 3527 0 
       31.7493769674 35.1784535489 -gcp 3527 2492 31.7431011291 35.1784951643
      -gcp 0 2494 31.7431761644 35.1680410195 map_overlay.png map_overlay.tif

entonces

  gdalwarp -s_srs epsg:4326 -t_srs epsg:4326  map_overlay.tif warped_map_overlay.tif

y ahora intente usar gdalinfo... si no ve ningún problema, ejecute este comando:

  gdal2tiles.py -z 14-21 warped_map_overlay.tif

si sus baldosas no encajan, ejecute el código anterior antes de actualizar, he dado.

Puede obtener más información aquí

espero que te ayude....