Понимание свойств создания Map Tile

Я пытаюсь создать карту основы плитки. У меня есть изображение, из которого я хочу создать плитку, и я использую GDAL.

Мне удалось создать плитки, но когда я вижу результат, кажется, что мозаичная карта:

1. Это перевернуто
2. Пропорции не верны

Это изображение наложения карты:

Это результат тайлинга:

Я много раз проверял, что координаты хороши, сравнивая углы изображения на ортофото.

Другое дело, что в Google Earth с такими же координатами изображение карты точно соответствует карте Google.

отредактированный

gdalinfo map_overlay.png

с этим выводом:

Драйвер: PNG / Portable Network Graphics Файлы: map_overlay.png Размер 3527, 2494 Система координат is '' Структура изображения Метаданные:
INTERLEAVE = PIXEL Угловые координаты: Верхний левый (0.0, 0.0) Нижний левый (0.0, 2494.0) Верхний правый (3527.0 , 0,0) Нижний правый (3527,0, 2494,0) Центр (1763,5, 1247,0) Блок 1 диапазона = 3527x1 Тип = байт, ColorInterp = Красные флаги маски: PER_DATASET ALPHA Блок 2 блока = 3527x1 Тип = байт, ColorInterp = Зеленые флаги маски: PER_DATASET ALPHA Блок 3 блока = 3527x1 тип = байт, ColorInterp =
флаги с синей маской: PER_DATASET ALPHA Блок 4 блока = 3527x1 тип = байт, ColorInterp = альфа

Затем переведите:

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 -ggp 0 2494 31.780804.

с этим выводом:

Размер входного файла 3527, 2494

И наконец:

gdal2tiles.py -z 14-21 map_overlay.vrt

Наложение - это простое изображение .png без каких-либо специальных данных.

Я буду признателен за любую помощь.

Шани

Попробуйте использовать метрическую скоординированную картинку, прежде чем использовать ее как:

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

ОБНОВИТЬ

Примечание: даже если мы укажем gcp, gdal_translateне будет указывать угловые координаты tiff.

Кроме того, вам не нужно создавать виртуальный набор данных (.vrt) из-за наличия только одного файла ...

используйте следующий код:

  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

тогда

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

и теперь попробуйте использовать gdalinfo ... если вы не видите никаких проблем, выполните эту команду:

  gdal2tiles.py -z 14-21 warped_map_overlay.tif

если ваши плитки не подходят, запустите код выше перед обновлением, я дал.

Вы можете получить больше информации здесь

я надеюсь, что это поможет вам ....