Comprender las propiedades de creación de mosaico de mapas

Estoy tratando de crear un mapa base de mosaico. Tengo una imagen desde la que quiero crear el mosaico y estoy usando GDAL.

Logré crear los mosaicos, pero cuando veo el resultado parece que el mapa en mosaico es:

1. Se voltea
2. Las proporciones no son correctas

Esta es la imagen de superposición del mapa:

Este es el resultado del mosaico:

Verifiqué 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 precisamente al mapa de Google.

EDITADO

gdalinfo map_overlay.png

con esta salida:

Controlador: PNG / Archivos de gráficos de red portátiles: map_overlay.png El tamaño es 3527, 2494 El sistema de coordenadas es `` Metadatos de estructura de imagen:
INTERLEAVE = PIXEL Coordenadas de esquina: Superior izquierda (0.0, 0.0) Inferior izquierda (0.0, 2494.0) Superior derecha (3527.0 , 0.0) Inferior derecha (3527.0, 2494.0) Centro (1763.5, 1247.0) Bloque de banda 1 = 3527x1 Tipo = Byte, ColorInterp = Indicadores de máscara rojos: PER_DATASET ALPHA Bloque de banda 2 = 3527x1 Tipo = Byte, ColorInterp = Indicadores de máscara verdes: PER_DATASET ALPHA Bloque de banda 3 = 3527x1 Tipo = Byte, ColorInterp =
Banderas de máscara azul : PER_DATASET ALPHA Bloque de banda 4 = 3527x1 Tipo = Byte, ColorInterp = Alfa

Luego traduce:

gdal_translate -of 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.7436164.png_44191910.png

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 imagen .png simple sin ningún dato especial.

Agradeceré cualquier ayuda.

Shani

Intente utilizar la imagen coordinada métrica antes de colocarla en mosaico como:

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

ACTUALIZAR

Nota: incluso si especificamos los gcp, gdal_translateno especificaría las coordenadas de esquina del tiff.

Además de esto, no necesita crear un conjunto de datos virtual (.vrt) debido a que solo tiene un archivo ...

usa 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 intenta usar gdalinfo ... si no ves ningún problema, ejecuta este comando:

  gdal2tiles.py -z 14-21 warped_map_overlay.tif

si sus mosaicos no encajan, ejecute el código anterior antes de la actualización, lo he dado.

Puedes obtener más información aquí

Espero que te ayude....