i.pansharpengrass - Интернет в облаке

ПРОГРАММА:

ИМЯ

я.паншарпен - Алгоритмы объединения изображений для повышения резкости мультиспектральных изображений с высоким разрешением
панхроматические каналы

КЛЮЧЕВЫЕ СЛОВА

изображения, слияние, резкость, Брови, IHS, HIS, PCA

СИНТАКСИС

я.паншарпен
я.паншарпен --Помогите
я.паншарпен [-sl] red=имя зеленый=имя синий=имя кастрюля=имя выходной=базовое имя метод=string
[-затирать] [-помощь] [-подробный] [-тихий] [-ui]

Флаги:
-s
Последовательная обработка, а не параллельная обработка

-l
Перебалансируйте синий канал для LANDSAT

- перезаписать
Разрешить выходным файлам перезаписывать существующие файлы

--Помогите
Распечатать сводку использования

--подробный
Подробный вывод модуля

--тихий
Тихий выход модуля

--уи
Принудительный запуск диалогового окна GUI

Параметры:
red=имя [требуется]
Имя растровой карты, которая будет использоваться для

зеленый=имя [требуется]
Имя растровой карты, которая будет использоваться для

синий=имя [требуется]
Имя растровой карты, которая будет использоваться для

кастрюля=имя [требуется]
Имя растровой карты, которая будет использоваться для панхроматического канала высокого разрешения

выходной=базовое имя [требуется]
Имя для выходной растровой карты (ов) с базовым именем

метод=string [требуется]
Метод заточки сковороды
Опции: Брови ихс, PCA
По умолчанию: IHS

ОПИСАНИЕ

я.паншарпен использует панхроматический диапазон высокого разрешения от мультиспектрального изображения до
резкость 3 полосы с более низким разрешением. Затем 3 полосы с более низким разрешением можно объединить в
цветное изображение RGB с более высоким (более детальным) разрешением, чем это возможно при использовании
оригинальные 3 полосы. Например, Landsat ETM имеет спектральные полосы низкого разрешения 1 (синий), 2
(зеленый), 3 (красный), 4 (ближний ИК), 5 (средний ИК) и 7 (средний ИК) с разрешением 30 м, а также высокий
разрешение панхроматический диапазон 8 при разрешении 15 м. Панельная заточка позволяет использовать ленты 3-2-1 (или
другие комбинации диапазонов разрешения 30 м, например 4-3-2 или 5-4-2) для объединения в 15-метровый
разрешение цветного изображения.
i.pansharpen предлагает на выбор три различных алгоритма панорамирования: IHS, Brovey,
и СПС.
Что касается IHS кастрюля заострение, исходные 3 полосы с более низким разрешением, выбранные как красный, зеленый и
синие каналы для создания составного изображения RGB преобразуются в IHS (интенсивность,
оттенок и насыщенность) цветовое пространство. Затем панхроматическая полоса заменяется на
канал интенсивности (I) в сочетании с исходными каналами оттенка (H) и насыщенности (S), и
преобразован обратно в цветовое пространство RGB с более высоким разрешением панхроматической полосы. В
алгоритм для этого может быть представлен как: RGB -> IHS -> [pan] HS -> RGB.
С Брови кастрюля заострение, каждая из 3 полос с более низким разрешением и панхроматическая полоса
объединяются с использованием следующего алгоритма для расчета 3 новых полос на более высоком
разрешение (пример для диапазона 1):
band1
новый band1 = ----------------------- * panband
группа1 + группа2 + группа3
In PCA кастрюля заострение, анализ главных компонент выполняется на исходных 3 нижних
полосы разрешения для создания трех основных компонентных изображений (ПК3, ПК1 и ПК2) и их
ассоциированные собственные векторы (EV), такие что:
группа1 группа2 группа3
PC1: EV1-1 EV1-2 EV1-3
PC2: EV2-1 EV2-2 EV2-3
PC3: EV3-1 EV3-2 EV3-3
и
PC1 = EV1-1 * band1 + EV1-2 * band2 + EV1-3 * band3 - среднее (полосы 1,2,3)
Затем выполняется обратный PCA, при котором панхроматическая полоса заменяется на PC1. Сделать это,
матрица собственных векторов инвертирована (в данном случае транспонирована), изображения ПК
умноженный на собственные векторы с панхроматической полосой, замененной PC1, и средним значением
каждая полоса добавляется к каждой полосе преобразованного изображения с использованием следующего алгоритма (пример
для диапазона 1):
диапазон1 '= панорама * EV1-1 + PC2 * EV2-1 + PC3 * EV3-1 + среднее значение (диапазон1)
Назначение каналов зависит от спутника. Примеры спутниковых снимков
с панхроматическими полосами высокого разрешения, а спектральные полосы с более низким разрешением включают
Landsat 7 ETM, QuickBird и SPOT.

ПРИМЕЧАНИЯ

В настоящее время модуль работает только с 8-битными изображениями.
Команда временно изменяет расчетную область на высокое разрешение
панхроматическая полоса во время расчета резкости, затем восстанавливает предыдущую область
настройки. Текущие координаты региона (и нулевые значения) соблюдаются. Высота
Панхроматическое изображение с разрешением - это гистограмма, соответствующая полосе, которую она заменяет до
подстановка (т.е. канал интенсивности для повышения резкости IHS, выбранная полоса низкого разрешения
для каждого цветового канала с повышением резкости по Брови и изображение PC1 для повышения резкости PCA).
По умолчанию команда будет пытаться использовать параллельную обработку, используя до 3 ядер.
одновременно. Флаг -s отключит параллельную обработку, но использует оптимизированный
Выражение r.mapcalc для уменьшения дискового ввода-вывода.
Три выходных канала с панорамированием резкости можно комбинировать с d.rgb or р.композит. Цвета
может быть дополнительно оптимизирован с помощью i.colors.enhance. В то время как полученное цветное изображение будет
быть с более высоким разрешением во всех случаях, 3 алгоритма панорамирования различаются по терминам
спектрального отклика.

ПРИМЕРЫ

Панорамирование изображения Landsat из Боулдера, Колорадо, США:
# R, G, B композит на 30 м
g.region raster=p034r032_7dt20010924_z13_10 -p
d.rgb b=p034r032_7dt20010924_z13_10 g=lp034r032_7dt20010924_z13_20
r=p034r032_7dt20010924_z13_30
# i.pansharpen с алгоритмом IHS
i.pansharpen red=p034r032_7dt20010924_z13_30 green=p034r032_7dt20010924_z13_20
blue=p034r032_7dt20010924_z13_10 pan=p034r032_7dp20010924_z13_80
output = ihs321 method = ihs
# дисплей на 15м
g.region raster = ihs321_blue -p
d.rgb b = ihs321_blue g = ihs321_green r = ihs321_red
Результаты:

R, G, B композитный of Landsat at 30m R, G, B композитный of Брови заточены изображение at 15m

R, G, B композитный of IHS заточены изображение at 15m R, G, B композитный of PCA заточены изображение at 15 м "

Пример: LANDSAT ETM + (Landsat 7), примерный набор данных Северная Каролина:
# оригинал на 28м
g.region raster = lsat7_2002_10 -p
д.мон wx0
d.rgb b=lsat7_2002_10 g=lsat7_2002_20 r=lsat7_2002_30
# i.pansharpen с алгоритмом IHS
i.pansharpen red = lsat7_2002_30 @ ПОСТОЯННЫЙ \
green=lsat7_2002_20 blue=lsat7_2002_10 \
pan = lsat7_2002_80 method = ihs \
output = lsat7_2002_ihs
# дисплей на 14.25м
g.region raster = lsat7_2002_ihs_red -p
г. стереть
d.rgb r = lsat7_2002_ihs_red g = lsat7_2002_ihs_green b = lsat7_2002_ihs_blue
# сравнить до / после (поддержка RGB в "Advanced"):
g.gui.mapswipe
# по желанию балансировка цвета:
i.colors.enhance r = lsat7_2002_ihs_red g = lsat7_2002_ihs_green b = lsat7_2002_ihs_blue

Используйте i.pansharpengrass в Интернете с помощью сервисов onworks.net