Установка программы GDAL

  1. Скачайте и запустите osgeo4w-setup-x86.exe или osgeo4w-setup-x86_64.exe для x64 по ссылке http://download.osgeo.org/osgeo4w/
  2. Программу нужно установить в стандартном режиме ("Стандартная Desktop установка").

  3. На втором шаге выберите сайт откуда загружать установку

  4. Оставить выбранным только пакет GDAL и установить

  5. Опять запустить установку в режиме "Расширенная установка" 

  6. На втором шаге - "Установить из интернет"

  7. На следующих шагах нажимаем "Далее"

  8. На этапе выбора пакета, раскройте ветку "Libs"

  9. Найдите пункт "2.149k python-numpy: Python NumPy (Numerical) Extension" и переключите его в "2.7-1.5.1-1" (в новых версиях может отличаться наименование и версия)

  10. Если программа выдаёт ошибку, что не хватает библиотеки MSVCR100.DLL, то возьмите два файла MSVCR100.DLL, MSVCP100.DLL и скопируйте их в папку C:\WINDOWS\SYSTEM32

Решение проблемы с нумерацией тайлов

Утилита gdal2tiles по умолчанию создаёт индексацию по стандарту OSGeo TMS (Tile Map Specification), чтобы формировать тайлы по стандарту Google OSM:

  1. Скачайте подправленный скрипт для правильной нумерации тайлов gdal2tiles.py (Google/OSM)

  2. Замените gdal2tiles в папке С:\OSGeo4W64\bin\

Создание тайлов (индексация)

Команда склеивания всех растров в папке test

gdalbuildvrt test.vrt -vrtnodata "255 255 255" -hidenodata -resolution highest test/*.jpg

Команда преобразования TIF файлов в RGB (если в файлах другая цветовая модель)

gdal_translate -of vrt -expand rgba test.tif test.vrt

Обратите внимание! Наименования файлов должны содержать только латинские символы. Файлы привязки растров должны иметь расширение tfw, jpw, sdw.

Команда преобразования растра в систему координат EPSG:4326

СК 63 - gdalwarp.exe -s_srs "+proj=tmerc +ellps=krass +lon_0=35d30E +lat_0=0d0N +k=1 +x_0=5300000 +y_0=-9214.69 +towgs84=+25,-141,-78.5,0.0,0.35,0.736,0.0" -t_srs EPSG:4326 -r bilinear test.vrt test.tif

EPSG:4326 - gdalwarp.exe -t_srs EPSG:4326 -r bilinear test.vrt test.tif (или можно пропустить этот шаг)

Команда построения тайлов. Результат будет сохранён в папке test

gdal2tiles -p mercator --s_srs EPSG:4326 test.tif

Отображение тайлов в Leaflet

Для отображения тайлов с использованием Leaflet:

TMS - L.tileLayer('http://{s}.somedomain.com/{foo}/{z}/{x}/{y}.png', {tms: true});

OSM - L.tileLayer('http://{s}.somedomain.com/{foo}/{z}/{x}/{y}.png', {});