User:Miroff/Orbview-3
Использование снимков Orbview-3 в JOSM
Про public domain вы наверняка уже в все в курсе.
Пререквизиты
1. Интернеты 2. GDAL 3. Nginx или любой другой вебсервер. Пользователи windows могут воспользоваться apache. 4. Горячая голова, чистое сердце и холодные руки.
Добываем
Исходные гранулы добываются через Earth Explorer, так что идем туда.
1. Регистрируемся, бо без регистрации нельзя. 2. Выбираем регион кликая мышкой. 3. Идем во вкладку "Data Sets" и выбираем Orbview-3 4. Идем во вкладку Additional Criteria и выбираем CLoud Level <30%, иначе вместо снимка придется смотреть на облака. 5. Нажимаем Results и выбираем что понравится. 6. Нажмаем Download и ждем.
Привязываем
Снимки Orbview-3 в исходном виде не привязаны и не орторектифицированы. Что это значит и как быть описано тут.
1. Добываем SRTM с гислаба. 2. Склеиваем SRTM в мозаику. 3. Привязываем снимок.
У меня получилась такая команда:
gdalwarp -co TILED=YES -co JPEG_QUALITY=85 -co COMPRESS=JPEG -multi -dstnodata 0 -srcnodata 0 -overwrite -ot byte -t_srs epsg:3857 -rpc -to 'RPC_DEM=srtm_53_02.tif' 3v041007p0000503101a520008301112m_001645635.tif 3v041007p0000503101a520008301112m_001645635.rec.tif
Уменьшаем динамический диапазон
Глубина цвета исходного снимка 16 бит. Нам этого много, потому нужно уменьшить.
Я использовал такую команду:
gdal_translate -scale 30 400 0 255 -ot Byte 3v041007p0000503101a520008301112m_001645635.rec.tif 3v041007p0000503101a520008301112m_001645635.trans.tif
С цифрами 30 и 400 можно поиграться.
Нарезаем
Поскольку JOSM не умеет GeoTIFF такого размера в качестве подложки, займемся нарезкой снимка на тайлы.
Я использовал такую команду:
gdal2tiles.py -w none -n -z 15-18 3v041007p0000503101a520008301112m_001645635.trans.tif tiles
последний аргумент это путь к папке с результатами
Публикуем
Я использовал nginx с примитивнейшим конфигом:
server { root /home/miroff/orbview-3/tiles; server_name localhost; location / { } }
Подключаем к OSM
Запускаем JOSM, идем в настройки и добавляем TMS источник:
tms[18]:http://localhost/{zoom}/{x}/{-y}.png
Затем загружаем нужную территорию, включаем подложку и проверяем привязку по трекам.