Паншарпенинг в QGIS с использованием Orfeo ToolBox: различия между версиями

Материал из GIS-Lab
Перейти к навигации Перейти к поиску
Нет описания правки
Нет описания правки
Строка 3: Строка 3:


==Введение==
==Введение==
При работе с данными дистанционного зондирования Земли зачастую приходится решать задачу улучшения снимков с низким разрешением за счёт снимков с высоким. В QGIS нет встроенного  
При работе с данными дистанционного зондирования Земли зачастую приходится решать задачу улучшения снимков с низким разрешением за счёт снимков с высоким. В QGIS нет встроенного инструмента паншарпенинга, но есть возможность подключения инструментов [https://www.orfeo-toolbox.org/ Orfeo ToolBox] (OTB), где они есть.


==Немного теории==
==Немного теории==
Строка 15: Строка 15:


где i и j - индексы пикселей.
где i и j - индексы пикселей.
==Паншарпенинг в OTB и QGIS==
С программной точки зрения процесс паншарпенинга делится на две стадии: 1) подготовительную, во время которой разрешение и экстент мультиспектрального растра приводится в соответствие с экстентом и разрешением панхроматического растра (это необходимо для проведения операций растровой алгебры, вовлекающих оба растра); 2) непосредственно паншарпенинг.
===Использование утилиты OTB===
OTB по своей сути является набором консольных утилит. За паншарпенинг отвечает "приложение" [https://www.orfeo-toolbox.org//CookBook/CookBooksu36.html#x57-770004.2.2 otbcli_BundleToPerfectSensor]. Оно совмещает в себе обе стадии, описанные выше, и запускается следующей консольной командой:
  otbcli_BundleToPerfectSensor -inp pan_image -inxs xs_image -out output_image
где pan_image, xs_image и output_image соответсвенно - пути к панхроматическому, мультиспектральному и результирующему растрам. О дополнительных опциях читайте в документации.
===Запуск паншарпенинга из QGIS===
В QGIS инструменты OTB подключаются через модуль Processing.
[[Файл:Enable OTB.png|center]]
==Ссылки по теме==
==Ссылки по теме==
[http://gis-lab.info/qa/pan-sharpening-erdas.html Паншарпенинг в ERDAS]
[http://gis-lab.info/qa/pan-sharpening-erdas.html Паншарпенинг в ERDAS]


[http://gis-lab.info/qa/grass7-landsat8-processing.html Обработка и интерпретация данных Landsat 8 (OLI) средствами GRASS GIS 7]
[http://gis-lab.info/qa/grass7-landsat8-processing.html Обработка и интерпретация данных Landsat 8 (OLI) средствами GRASS GIS 7]

Версия от 10:53, 4 апреля 2015

Эта страница является черновиком статьи.


В данной статье рассказывается о том как произвести паншарпенинг в QGIS посредством инструментария Orfeo ToolBox.

Введение

При работе с данными дистанционного зондирования Земли зачастую приходится решать задачу улучшения снимков с низким разрешением за счёт снимков с высоким. В QGIS нет встроенного инструмента паншарпенинга, но есть возможность подключения инструментов Orfeo ToolBox (OTB), где они есть.

Немного теории

На удивление сложно найти определение для слова "паншарпенинг". Попробуем сформулировать его своими словами взяв за основу описание из книги Шовергердт Р.А. Дистанционное зондирование. Модели и методы обработки изображений - М.: Техносфера, 2010. (с. 416):

Паншарпенинг (от анг. Panchromatic sharpening) - это процесс объединения изображений в пространственной области основная задача которого заключается в передаче высокочастотного содержания изображения с высоким разрешением (обычно панхроматического) изображению с низким разрешением (обычно мультиспектральному).

Одним из наиболее простых алгоритмов паншарпенинга является модуляция высоких частот. В этом алгоритме для получения улучшенного изображения в канале PXS панхроматический снимок PAN попиксельно умножается на изображение с низким разрешением XS, после чего полученный результат нормируется на низкочастотную компоненту панхроматического снимка PANsmooth (панхроматический снимок, обработанный сглаживающим фильтром с окном, соответствующим размеру пикселя изображения с низким разрешением):

Formula OTB pansharpening.png

где i и j - индексы пикселей.

Паншарпенинг в OTB и QGIS

С программной точки зрения процесс паншарпенинга делится на две стадии: 1) подготовительную, во время которой разрешение и экстент мультиспектрального растра приводится в соответствие с экстентом и разрешением панхроматического растра (это необходимо для проведения операций растровой алгебры, вовлекающих оба растра); 2) непосредственно паншарпенинг.

Использование утилиты OTB

OTB по своей сути является набором консольных утилит. За паншарпенинг отвечает "приложение" otbcli_BundleToPerfectSensor. Оно совмещает в себе обе стадии, описанные выше, и запускается следующей консольной командой:

 otbcli_BundleToPerfectSensor -inp pan_image -inxs xs_image -out output_image

где pan_image, xs_image и output_image соответсвенно - пути к панхроматическому, мультиспектральному и результирующему растрам. О дополнительных опциях читайте в документации.

Запуск паншарпенинга из QGIS

В QGIS инструменты OTB подключаются через модуль Processing.

Enable OTB.png

Ссылки по теме

Паншарпенинг в ERDAS

Обработка и интерпретация данных Landsat 8 (OLI) средствами GRASS GIS 7