Как паншарпить 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-композиты и панхроматика


Первоначально я наткнулся на ту статью, на которую вы ссылаетесь с 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. В идеале ответы, подобные вашим, должны иметь описание того, почему они являются хорошими методами. Ответы только на ссылки, как правило, перестают существовать, если ссылка исчезает.
Аарон
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.