Тайлы из Landsat 8: различия между версиями
Zverik (обсуждение | вклад) (→Паншарпенинг: вам пакет) |
Zverik (обсуждение | вклад) (→Паншарпенинг: коэфф. lum) |
||
| Строка 17: | Строка 17: | ||
У нас есть 3 канала с разрешением 15 м/пк и один — 30 м/пк. Паншерпенинг — это объединение их всех в один цветной снимок максимального разрешения. Раньше это делали длинным цепочками команд GDAL или ImageMagick, но нам повезло. Умные люди из GINA (университет Аляски) написали [https://github.com/gina-alaska/dans-gdal-scripts/ скрипт], и теперь достаточно одной команды. Набор аляскинских скриптов есть в большинстве дистрибутивов, поищите пакет ''dans-gdal-scripts''. Команда такая: | У нас есть 3 канала с разрешением 15 м/пк и один — 30 м/пк. Паншерпенинг — это объединение их всех в один цветной снимок максимального разрешения. Раньше это делали длинным цепочками команд GDAL или ImageMagick, но нам повезло. Умные люди из GINA (университет Аляски) написали [https://github.com/gina-alaska/dans-gdal-scripts/ скрипт], и теперь достаточно одной команды. Набор аляскинских скриптов есть в большинстве дистрибутивов, поищите пакет ''dans-gdal-scripts''. Команда такая: | ||
gdal_landsat_pansharp -pan B8.TIF -ndv 0 -o result.tif -rgb B4.TIF -rgb B3.TIF -rgb B2.TIF | gdal_landsat_pansharp -pan B8.TIF -ndv 0 -o result.tif \ | ||
-rgb B4.TIF -rgb B3.TIF -rgb B2.TIF -lum B3.TIF 0.5 -lum B4.TIF 0.5 | |||
Я опустил длинные идентификаторы снимков, вида «LC81850172014281LGN00». С ними будет длиннее и менее понятно: | Я опустил длинные идентификаторы снимков, вида «LC81850172014281LGN00». С ними будет длиннее и менее понятно: | ||
| Строка 24: | Строка 24: | ||
gdal_landsat_pansharp -pan LC81850172014281LGN00_B8.TIF -ndv 0 -o result.tif \ | gdal_landsat_pansharp -pan LC81850172014281LGN00_B8.TIF -ndv 0 -o result.tif \ | ||
-rgb LC81850172014281LGN00_B4.TIF -rgb LC81850172014281LGN00_B3.TIF -rgb LC81850172014281LGN00_B2.TIF \ | -rgb LC81850172014281LGN00_B4.TIF -rgb LC81850172014281LGN00_B3.TIF -rgb LC81850172014281LGN00_B2.TIF \ | ||
-lum LC81850172014281LGN00_B3.TIF 0.5 -lum LC81850172014281LGN00_B4.TIF 0.5 | |||
Что тут происходит? Первые три параметра передают панхроматический канал (восьмой), значение «нет данных» и имя выходного файла. Следующая тройка с <code>-ref</code> — каналы для красного, зелёного и синего цветов. В Landsat 8 им соответствуют 4, 3 и 2, но можно выбрать другие, из [http://gis-lab.info/qa/landsat-bandcomb.html невидимых областей спектра]. | Что тут происходит? Первые три параметра передают панхроматический канал (восьмой), значение «нет данных» и имя выходного файла. Следующая тройка с <code>-ref</code> — каналы для красного, зелёного и синего цветов. В Landsat 8 им соответствуют 4, 3 и 2, но можно выбрать другие, из [http://gis-lab.info/qa/landsat-bandcomb.html невидимых областей спектра]. | ||
Наконец, параметры <code>-lum</code> корректируют наложение цвета на панхроматический канал. [https://github.com/gina-alaska/dans-gdal-scripts/wiki/Gdal_landsat_pansharp Описание алгоритма] слишком сложное, мы | Наконец, параметры <code>-lum</code> корректируют наложение цвета на панхроматический канал. [https://github.com/gina-alaska/dans-gdal-scripts/wiki/Gdal_landsat_pansharp Описание алгоритма] слишком сложное, мы взяли коэффициенты из [https://gist.github.com/dwtkns/d18828165160384d0fe1 этих функций] Дерека Уоткинса. | ||
== Загрузка в QGIS == | == Загрузка в QGIS == | ||
Версия от 15:59, 1 февраля 2015
Главное свойство ДЗЗ Landsat 8 — свободная лицензия, которая позволяет обрисовывать снимки для, в частности, OpenStreetMap. Эта статья не вдаётся в теорию, но содержит шаги для поиска и обработки снимков Landsat 8 и их превращения в тайлы для публикации или обрисовки в JOSM.
Получение снимка
Раньше нужно было шоркаться по EarthExplorer, вводить параметры поиска, что-то выбирать. Спасибо DevelopmentSeed, много лет работающей с Landsat, теперь выбрать нужный снимок проще. Заходите, приближайте свою область, кликайте в точку, фильтруйте и скачивайте:
http://libra.developmentseed.org
Напоминаем, что каждый снимок весит 500-800 мегабайт в архиве (2 гигабайта вне его), и для обработки потребуются ещё гигабайта четыре.
В архиве вы получите 11 каналов и слой QA, на котором отмечены точки без информации о земной поверхности (например, облака). Мы будем работать с каналами 2, 3, 4 (синий, зелёный, красный) и 8 (ч/б с 15-метровым разрешением).
Паншарпенинг
У нас есть 3 канала с разрешением 15 м/пк и один — 30 м/пк. Паншерпенинг — это объединение их всех в один цветной снимок максимального разрешения. Раньше это делали длинным цепочками команд GDAL или ImageMagick, но нам повезло. Умные люди из GINA (университет Аляски) написали скрипт, и теперь достаточно одной команды. Набор аляскинских скриптов есть в большинстве дистрибутивов, поищите пакет dans-gdal-scripts. Команда такая:
gdal_landsat_pansharp -pan B8.TIF -ndv 0 -o result.tif \ -rgb B4.TIF -rgb B3.TIF -rgb B2.TIF -lum B3.TIF 0.5 -lum B4.TIF 0.5
Я опустил длинные идентификаторы снимков, вида «LC81850172014281LGN00». С ними будет длиннее и менее понятно:
gdal_landsat_pansharp -pan LC81850172014281LGN00_B8.TIF -ndv 0 -o result.tif \ -rgb LC81850172014281LGN00_B4.TIF -rgb LC81850172014281LGN00_B3.TIF -rgb LC81850172014281LGN00_B2.TIF \ -lum LC81850172014281LGN00_B3.TIF 0.5 -lum LC81850172014281LGN00_B4.TIF 0.5
Что тут происходит? Первые три параметра передают панхроматический канал (восьмой), значение «нет данных» и имя выходного файла. Следующая тройка с -ref — каналы для красного, зелёного и синего цветов. В Landsat 8 им соответствуют 4, 3 и 2, но можно выбрать другие, из невидимых областей спектра.
Наконец, параметры -lum корректируют наложение цвета на панхроматический канал. Описание алгоритма слишком сложное, мы взяли коэффициенты из этих функций Дерека Уоткинса.
Загрузка в QGIS
http://gis-lab.info/qa/qtiles.html
Коррекция цвета
Тайлы
- gdal2tiles
- подключение в JOSM