Установка программы 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

Команда преобразования растра в систему координат 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', {});