Как паншарпить LANDSAT 8 в GRASS?

9

Пользуясь QGIS в течение некоторого времени для картографирования, я был «вынужден» копаться в GRASS, так как панхарпинг напрямую не выполняется в QGIS. После более или менее понимания структуры я провел несколько испытаний с i.fusion.brovey. Мне удалось сгенерировать растры с разрешением 15 м, но при загрузке в QGIS я не могу добиться / восстановить хороший цветовой контраст, который присутствует в исходных 30-метровых растрах.

В предыдущем посте кто-то подтвердил это и предложил использовать другую функцию для landsat: «i.his.rgb» (после использования i.landsat.toar, который я не использую). Он начинает конвертировать мин / макс в 1-255 перед использованием i.his.rgb ... У моих оригинальных растров диапазон 0-65535. Это не сработало для меня ...

Я ушел после нескольких дней изо всех сил, чтобы просто выполнить pansharp, я чувствую себя застрявшим ... Почему просто применение i.pansharp не работает для Landsat, как предлагали другие? Может кто-нибудь предложить более простой рабочий процесс для выполнения панорамирования изображений Landsat 8? Ждем несколько советов, заранее спасибо ...

Арно Куппенс
источник

Ответы:

7

Взгляните на i.landsat.rgb - выполняет автоматическую балансировку цветов для изображений LANDSAT, возможно, перед запуском панширования. Вы можете также рассмотреть возможность преобразования цифровых номеров отдельных каналов в яркость или отражательную способность верхней части атмосферы с помощью i.landsat.toar . Смотрите также http://grasswiki.osgeo.org/wiki/LANDSAT

Кстати, наличие диапазона 0-65535 для исходных данных LANDSAT является нечетным, поскольку они доставляются как 8-битные каналы (0-255). Поэтому сначала проверьте ваши входные данные.

Обновление: 10/2013

Для упрощенной процедуры с текущей GRASS GIS 7, см. Страницу введите описание изображения здесьОбработка данных Landsat 8 в GRASS GIS 7: RGB-композиты и панхроматика

markusN
источник
Первоначально я наткнулся на ту статью, на которую вы ссылаетесь с neteler.org, но инструкции не работали вообще (из-за формата с плавающей точкой, созданного с помощью растров i.landsat.toar, поэтому я пришел
Хуан
3
Обратите внимание, что этот совет был дан много лет назад и не отражает новейший код GRASS GIS. Кстати, действительно отличный модуль панширования - это новое дополнение: i.fusion.hpf - объединение панхроматических и низкоразрешающих мультиспектральных данных с высоким разрешением на основе метода сложного
markusN
4

Спасибо Михал и Маркус. Наконец-то мне удалось сделать заточку с помощью указанной функции brov. После выполнения в GRASS я экспортировал растры RGB в gdal GTiff с типом данных Uint16. Маркус, я загружаю и использую несколько изображений Landsat 8 (tiff каждой из 11 полос), получая хорошие визуальные результаты, как в естественном / ложном цвете. Насколько мне известно, оригинальные tiffs landat имеют вид Uint16 со значениями от 0 до 65535 (или я что-то упустил?) ... Спасибо за ссылку GRASS LANDSAT ...

Арно Куппенс
источник
3
Да, изображения Landsat 8 хранятся в формате Uint16.
webrian
3

Обработка резкой привязки FIHS или Brovery Transform довольно проста, если вы следуете формулам из такой статьи, как эта: http://www.sensorsportal.com/HTML/DIGEST/august_2014/Vol_177/P_RP_0183.pdf

Предполагая, что вы используете GRASS GIS и pansharpening Landsat 8.

Установите регион с помощью панорамирования 8 - важно иметь разрешение 15 м, которое будет установлено, если вы сделаете это.

Для Бровери без весов введите это в Map Calculator для синего, зеленого и красного соответственно:

band8 / (band2 + band3 + band4) / 3 * band2

band8 / (band2 + band3 + band4) / 3 * band3

band8 / (band2 + band3 + band4) / 3 * band4

Запустите r.colors для каждого из вышеприведенных выходов и установите выравнивание гистограммы и таблицу цветов на серый.

Запустите r.composite, выберите соответствующие полосы (красный = 4, зеленый = 3, синий = 2)

Проверьте изображение.

Очевидно, что этот процесс мог бы быть сделан с большим вниманием к атмосферной коррекции, растяжению цвета и весам полос, как обсуждалось в статье выше.

Бен
источник
2

Вы также можете попробовать Монтеверди выполнить заточку. Этот инструмент может использовать обычные GeoTIFF в качестве входных данных (однако отдельные многоспектральные полосы должны быть объединены в один многоканальный растр) и работает достаточно хорошо.

Михал Мацкевич
источник
2

После борьбы с этой проблемой я нашел простое решение, которое позволяет вам панорамировать в QGIS (2.18.21 - моя версия) с использованием SCP. Обновление SCP, включающее этот метод, описано здесь https://fromgistors.blogspot.com/2015/07/major-update-semi-automatic-44.html .

Это легко сделать, и это бежало как очарование! Я проверил окончательный результат, и разрешение пикселей составляет 15 м, что соответствует панхроматической 8 полосе. Если кому-то интересно - он использует метод преобразования Брови, как описано в руководстве (см. Https://media.readthedocs.org/pdf/semiautomaticclassificationmanual/latest/semiautomaticclassificationmanual.pdf ).

Шаг 1 для инструмента SCP Pansharpen Шаг 3-4 для SCP pansharpen

Марк Томпсон
источник
1

Добро пожаловать на сайт Janos. В идеале ответы, подобные вашим, должны иметь описание того, почему они являются хорошими методами. ..

Для панорамы снимков Landsat 8 более простым способом является http://www.geosage.com/highview/download.html . Мой совет, попробуйте. Только один клик ....Spectraltransformer

Янош Брезной
источник
-2

Для паншарпенига я использовал спектральный преобразователь с графическим интерфейсом для изображений Landsat 8 введите описание изображения здесь http://www.geosage.com/highview/download.html Очень хорошо. БР Янос

Янош Брезной
источник
2
Добро пожаловать на сайт Janos. В идеале ответы, подобные вашим, должны иметь описание того, почему они являются хорошими методами. Ответы только на ссылки, как правило, перестают существовать, если ссылка исчезает.
Аарон