Паншарпенинг данных Landsat 8 средствами свободного ПО

Материал из GIS-Lab
Перейти к навигации Перейти к поиску
Эта страница является черновиком статьи.


Здесь будет описано (если будет), как можно сделать паншарпенинг данных Landsat 8 с помощью инструментов различного свободного ПО.

По мотивам этой темы на форуме.

Что такое "паншарпенинг"?

Паншарпенинг (от англ. "panchromatic sharpening" или "pan sharpening") — процесс, позволяющий получить из панхроматического и мультиспектрального каналов двух продуктов одно изображение. Панхроматический канал, как правило, имеет более высокое пространственное разрешение, мультиспектральный — более низкое. В результате слияния каналов получается цветное изображение с высоким разрешением [1].

Теория о методе, когда впервые появился, кто создатель - мне кажется, это лишнее (А.М.)

Методы паншарпенинга - лучше кратко (А.М.)

  • brovey
  • IHS
  • PCA


Кратко о Landsat 8

Тестовые данные

В качестве тестовых данных используется сцена Landsat 8 от 29.06.2013 (Entity ID: LC81710172013180LGN00). Получить снимки можно с помощью [EarthExplorer http://earthexplorer.usgs.gov/]. На вкладке "Search Criteria" в графе "Path/Row" вводим соответственно "170" для "Path" и "17" для "Row", в графе "Data Range" для "Search from:" и "to:" вписываем дату снимка "06/29/2013" ; затем на вкладке "Data Sets" выбираем "Landsat Archive" -> "L8 OLI/TIRS"; выбор закончен, жмём "Results>>". На странице результатов жмём на иконку с зелёной стрелкой "Download Options", в появившемся окне выбираем нижний вариант (Level 1 GeoTIFF Data Product (821.9 MB)) -> "Select Download Option" -> "Download". Внимание: загружаемый архив весит около 822 Мб, а в разархирированном виде сцена занимает на диске около 2 Гб.

Условия выполнения

Данные Landsat 8 поставляются в 16-битном растровом формате (Tiff). Алгоритмы GDAL и GRASS работают пока только с 8-битными форматами данных. Поэтому для получения растров повышенного разрешения необходимо сначала преобразовать входные данные в 8-битовый формат. Это точно? Или речь идет о каких-то конкретных алгоритмах? (Просто я постоянно работаю с 16-битными данными и никаких проблем не встречал -- КД)

Паншарпенинг средствами свободном ПО

GRASS

См. http://courses.neteler.org/processing-landsat-8-data-in-grass-gis-7-rgb-composites-and-pan-sharpening/

  1. r.rescale.eq
  2. i.pansharpen

GDAL

  • допиленный gdal_calc.py + gdal_contrast_stretch (dans-gdal-scripts)

Orfeo Toolbox / Monteverdi

http://live.osgeo.org/ru/overview/otb_overview.html

http://www.orfeo-toolbox.org/Applications/Pansharpening.html otbcli_Pansharpening

Источник — https://wiki.gis-lab.info/index.php?title=Паншарпенинг_данных_Landsat_8_средствами_свободного_ПО&oldid=16489