Как построить тайлы для отображения на Leaflet.js картах. Решение проблемы с нумерацией тайлов
Установка программы GDAL
- Скачайте и запустите osgeo4w-setup-x86.exe или osgeo4w-setup-x86_64.exe для x64 по ссылке http://download.osgeo.org/osgeo4w/
-
Программу нужно установить в стандартном режиме ("Стандартная Desktop установка").
-
На втором шаге выберите сайт откуда загружать установку
-
Оставить выбранным только пакет GDAL и установить
-
Опять запустить установку в режиме "Расширенная установка"
-
На втором шаге - "Установить из интернет"
-
На следующих шагах нажимаем "Далее"
-
На этапе выбора пакета, раскройте ветку "Libs"
-
Найдите пункт "2.149k python-numpy: Python NumPy (Numerical) Extension" и переключите его в "2.7-1.5.1-1" (в новых версиях может отличаться наименование и версия)
-
Если программа выдаёт ошибку, что не хватает библиотеки MSVCR100.DLL, то возьмите два файла MSVCR100.DLL, MSVCP100.DLL и скопируйте их в папку C:\WINDOWS\SYSTEM32
Решение проблемы с нумерацией тайлов
Утилита gdal2tiles по умолчанию создаёт индексацию по стандарту OSGeo TMS (Tile Map Specification), чтобы формировать тайлы по стандарту Google OSM:
- Скачайте подправленный скрипт для правильной нумерации тайлов gdal2tiles.py (Google/OSM)
- Замените 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', {});