Просмотр исходного текста страницы Тайлы из Landsat 8

== Паншарпенинг ==

У нас есть 3 канала с разрешением 30 м/пк и один — 15 м/пк. Паншарпенинг — это объединение их всех в один цветной снимок максимального разрешения. Раньше его делали длинным цепочками команд 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 -lum B3.TIF 0.5 -lum B4.TIF 0.5

Здесь опущены длинные идентификаторы снимков, вида «LC81850172014281LGN00». Для удобства пригодится такой скрипт:

 BASE=LC81850172014281LGN00
 gdal_landsat_pansharp -pan ${BASE}_B8.TIF -ndv 0 -o result.tif \
  -rgb ${BASE}_B4.TIF -rgb ${BASE}_B3.TIF -rgb ${BASE}_B2.TIF \
  -lum ${BASE}_B3.TIF 0.5 -lum ${BASE}_B4.TIF 0.5

Что тут происходит? Первые три параметра передают панхроматический канал (восьмой), значение «нет данных» и имя выходного файла. Следующая тройка с <code>-rgb</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 Описание алгоритма] слишком сложное, мы взяли коэффициенты из [https://gist.github.com/dwtkns/d18828165160384d0fe1 этих функций] Дерека Уоткинса.