Паншерпенинг Ландсата 8

Аватара пользователя
rhot
Гуру
Сообщения: 1727
Зарегистрирован: 25 янв 2011, 17:50
Репутация: 194
Ваше звание: доктор
Откуда: Архангельск

Паншерпенинг Ландсата 8

Сообщение rhot » 26 ноя 2013, 15:11

У кого-нибудь получилось в GRASS сделать более или менее сносный паршерпенинг снимков Landsat 8?
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________

bolotoved
Гуру
Сообщения: 920
Зарегистрирован: 30 дек 2008, 14:11
Репутация: 236
Откуда: Ханты-Мансийск
Контактная информация:

Re: Паншерпенинг Ландсата 8

Сообщение bolotoved » 26 ноя 2013, 15:35

Да, с пом. какого-то стандартного алгоритма, сейчас, сходу не могу найти куда записал команду :cry:
Вложения
Selection_025.png
Selection_025.png (72.19 КБ) 11014 просмотров

Аватара пользователя
rhot
Гуру
Сообщения: 1727
Зарегистрирован: 25 янв 2011, 17:50
Репутация: 194
Ваше звание: доктор
Откуда: Архангельск

Re: Паншерпенинг Ландсата 8

Сообщение rhot » 26 ноя 2013, 15:37

Снимки на вход были после атмосферной коррекции? Там три алгоритма: brovey, PCA, IHS.
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________

bolotoved
Гуру
Сообщения: 920
Зарегистрирован: 30 дек 2008, 14:11
Репутация: 236
Откуда: Ханты-Мансийск
Контактная информация:

Re: Паншерпенинг Ландсата 8

Сообщение bolotoved » 26 ноя 2013, 17:07

rhot писал(а):Снимки на вход были после атмосферной коррекции? Там три алгоритма: brovey, PCA, IHS.
Использовал http://grass.osgeo.org/grass70/manuals/ ... arpen.html Метод IHS.
Без коррекции, т.к. это все равно только для визуадизации можно использовать.

bolotoved
Гуру
Сообщения: 920
Зарегистрирован: 30 дек 2008, 14:11
Репутация: 236
Откуда: Ханты-Мансийск
Контактная информация:

Re: Паншерпенинг Ландсата 8

Сообщение bolotoved » 26 ноя 2013, 17:12

Там только для начала нужно прогнать Landsat 8 через GDAL, чтобы получить 8bit.

Аватара пользователя
rhot
Гуру
Сообщения: 1727
Зарегистрирован: 25 янв 2011, 17:50
Репутация: 194
Ваше звание: доктор
Откуда: Архангельск

Re: Паншерпенинг Ландсата 8

Сообщение rhot » 26 ноя 2013, 17:15

Странно...
Без коррекции вообще пустые растры выдаёт.
Версия GRASS какая? У меня:

Код: Выделить всё

GRASS version: 7.0.svn                                                          
GRASS SVN Revision: 58304                                                       
Build Date: 2013-11-25                                                          
GIS Library Revision: 58048 (2013-10-18)                                        
GDAL/OGR: 1.9.0                                                                 
PROJ.4: 4.7.0                                                                   
GEOS:                                                                           
SQLite: 3.7.13                                                                  
Python: 2.7.3                                                                   
wxPython: 2.8.12.1                                                              
Platform: Linux-3.2.0-4-amd64-x86_64-with-debian-7.2  
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________

Аватара пользователя
rhot
Гуру
Сообщения: 1727
Зарегистрирован: 25 янв 2011, 17:50
Репутация: 194
Ваше звание: доктор
Откуда: Архангельск

Re: Паншерпенинг Ландсата 8

Сообщение rhot » 26 ноя 2013, 17:22

bolotoved писал(а):Там только для начала нужно прогнать Landsat 8 через GDAL, чтобы получить 8bit.
Не подскажите команду? :)
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________

bolotoved
Гуру
Сообщения: 920
Зарегистрирован: 30 дек 2008, 14:11
Репутация: 236
Откуда: Ханты-Мансийск
Контактная информация:

Преобразование 16-bit d 8 bit b паншарп Landsat-8

Сообщение bolotoved » 26 ноя 2013, 18:01

Сорри, похоже, что я все делал с пом. GDAL, т.ч. не удивительно, что я не смог найти грассовского скрипта.

Для начала нужно установить спец. утилиту:
https://github.com/gina-alaska/dans-gdal-scripts

1) Преобразуем 16-bit в 8-bit, c выравниванием гистограммы:

Код: Выделить всё

gdal_contrast_stretch -histeq 32768 LC81670152013184LGN00_B4.TIF LC81670152013184LGN00_B4_8bit.tif 
gdal_contrast_stretch -histeq 32768 LC81670152013184LGN00_B5.TIF LC81670152013184LGN00_B5_8bit.tif 
gdal_contrast_stretch -histeq 32768 LC81670152013184LGN00_B6.TIF LC81670152013184LGN00_B6_8bit.tif 
gdal_contrast_stretch -histeq 32768 LC81670152013184LGN00_B8.TIF LC81670152013184LGN00_B8_8bit.tif


2) Склеиваем RGB-композит (4, 5 и 6 каналы в моем случае)

3)Теперь сам паншарп по 8-му:

Код: Выделить всё

gdal_landsat_pansharp -rgb LC81670152013184LGN00_8bit.TIF -pan LC81670152013184LGN00_B8_8bit.tif -o out3.tif

Аватара пользователя
gimran
Гуру
Сообщения: 1902
Зарегистрирован: 07 июл 2010, 15:43
Репутация: 242
Откуда: Уфа

Re: Паншерпенинг Ландсата 8

Сообщение gimran » 26 ноя 2013, 18:50

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

[ Сообщение с мобильного устройства ]

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 793
Ваше звание: званий не имею
Откуда: Москва

Re: Паншерпенинг Ландсата 8

Сообщение Александр Мурый » 26 ноя 2013, 19:09

gimran писал(а):Мжно в кугисе калькулятором растров, экспериментировал для того, чтобы далеким от грасса людям обьяснить тупейший способ.
Поясните, пож-та, недалёким (от грасса) людям, как. :)
Редактор материалов, модератор форума

Strix
Интересующийся
Сообщения: 31
Зарегистрирован: 23 сен 2013, 18:45
Репутация: 7

Re: Паншерпенинг Ландсата 8

Сообщение Strix » 26 ноя 2013, 19:11

В gdal_contrast_stretch со значением параметра -histeq 32768 режет гистограмму (по крайней мере у меня); экспериментальным путем установил, что хорошие результаты дает при значении 80-100.
Если приемлемо преобразование Бровея, то после допиливания gdal_calc.py по рецепту Voltrona (http://blogs.gis-lab.info/apple-touch-i ... /?paged=45), этот метод можно также применять в GDAL вместе с dans-gdal-scripts.

Аватара пользователя
rhot
Гуру
Сообщения: 1727
Зарегистрирован: 25 янв 2011, 17:50
Репутация: 194
Ваше звание: доктор
Откуда: Архангельск

Re: Паншерпенинг Ландсата 8

Сообщение rhot » 26 ноя 2013, 19:20

В GRASS никак нельзя преобразовать 16-бит в 8-бит?
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________

Аватара пользователя
gimran
Гуру
Сообщения: 1902
Зарегистрирован: 07 июл 2010, 15:43
Репутация: 242
Откуда: Уфа

Re: Паншерпенинг Ландсата 8

Сообщение gimran » 26 ноя 2013, 19:32

DN_b1
DN_fused = ------------------------------- * DN_pan
DN_b1 + DN_b2 + DN_b3

Где b1, b2, b3 каналы rgb, pan - панхром
Это для первого канала.
Для второго и третьего тоже самое, только вместо b1 ставим b2 и b3 в числитель, соответственно.
А потом полученные растры в композит.

[ Сообщение с мобильного устройства ]

Аватара пользователя
gimran
Гуру
Сообщения: 1902
Зарегистрирован: 07 июл 2010, 15:43
Репутация: 242
Откуда: Уфа

Re: Паншерпенинг Ландсата 8

Сообщение gimran » 26 ноя 2013, 19:34

То есть первый = (b1 × b_pan)/(b1+b2+b3)
Второй = (b2 × b_pan)/(b1+b2+b3)
Третий = (b3 × b_pan)/(b1+b2+b3)

[ Сообщение с мобильного устройства ]

bolotoved
Гуру
Сообщения: 920
Зарегистрирован: 30 дек 2008, 14:11
Репутация: 236
Откуда: Ханты-Мансийск
Контактная информация:

Re: Паншерпенинг Ландсата 8

Сообщение bolotoved » 27 ноя 2013, 01:37

rhot писал(а):В GRASS никак нельзя преобразовать 16-бит в 8-бит?
Можно: http://grass.osgeo.org/grass70/manuals/r.rescale.html
Но нужно разбираться как правильно настроить гистограмму.

Ответить

Вернуться в «GRASS»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 7 гостей