Паншерпенинг Ландсата 8
- rhot
- Гуру
- Сообщения: 1727
- Зарегистрирован: 25 янв 2011, 17:50
- Репутация: 194
- Ваше звание: доктор
- Откуда: Архангельск
Паншерпенинг Ландсата 8
У кого-нибудь получилось в GRASS сделать более или менее сносный паршерпенинг снимков Landsat 8?
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________
-
- Гуру
- Сообщения: 920
- Зарегистрирован: 30 дек 2008, 14:11
- Репутация: 236
- Откуда: Ханты-Мансийск
- Контактная информация:
Re: Паншерпенинг Ландсата 8
Да, с пом. какого-то стандартного алгоритма, сейчас, сходу не могу найти куда записал команду
- Вложения
-
- Selection_025.png (72.19 КБ) 11014 просмотров
- rhot
- Гуру
- Сообщения: 1727
- Зарегистрирован: 25 янв 2011, 17:50
- Репутация: 194
- Ваше звание: доктор
- Откуда: Архангельск
Re: Паншерпенинг Ландсата 8
Снимки на вход были после атмосферной коррекции? Там три алгоритма: brovey, PCA, IHS.
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________
-
- Гуру
- Сообщения: 920
- Зарегистрирован: 30 дек 2008, 14:11
- Репутация: 236
- Откуда: Ханты-Мансийск
- Контактная информация:
Re: Паншерпенинг Ландсата 8
Использовал http://grass.osgeo.org/grass70/manuals/ ... arpen.html Метод IHS.rhot писал(а):Снимки на вход были после атмосферной коррекции? Там три алгоритма: brovey, PCA, IHS.
Без коррекции, т.к. это все равно только для визуадизации можно использовать.
-
- Гуру
- Сообщения: 920
- Зарегистрирован: 30 дек 2008, 14:11
- Репутация: 236
- Откуда: Ханты-Мансийск
- Контактная информация:
Re: Паншерпенинг Ландсата 8
Там только для начала нужно прогнать Landsat 8 через GDAL, чтобы получить 8bit.
- rhot
- Гуру
- Сообщения: 1727
- Зарегистрирован: 25 янв 2011, 17:50
- Репутация: 194
- Ваше звание: доктор
- Откуда: Архангельск
Re: Паншерпенинг Ландсата 8
Странно...
Без коррекции вообще пустые растры выдаёт.
Версия GRASS какая? У меня:
Без коррекции вообще пустые растры выдаёт.
Версия 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
Не подскажите команду?bolotoved писал(а):Там только для начала нужно прогнать Landsat 8 через GDAL, чтобы получить 8bit.
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________
-
- Гуру
- Сообщения: 920
- Зарегистрирован: 30 дек 2008, 14:11
- Репутация: 236
- Откуда: Ханты-Мансийск
- Контактная информация:
Преобразование 16-bit d 8 bit b паншарп Landsat-8
Сорри, похоже, что я все делал с пом. GDAL, т.ч. не удивительно, что я не смог найти грассовского скрипта.
Для начала нужно установить спец. утилиту:
https://github.com/gina-alaska/dans-gdal-scripts
1) Преобразуем 16-bit в 8-bit, c выравниванием гистограммы:
2) Склеиваем RGB-композит (4, 5 и 6 каналы в моем случае)
3)Теперь сам паншарп по 8-му:
Для начала нужно установить спец. утилиту:
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
Мжно в кугисе калькулятором растров, экспериментировал для того, чтобы далеким от грасса людям обьяснить тупейший способ.
[ Сообщение с мобильного устройства ]
[ Сообщение с мобильного устройства ]
-
- Гуру
- Сообщения: 5173
- Зарегистрирован: 26 сен 2009, 16:26
- Репутация: 793
- Ваше звание: званий не имею
- Откуда: Москва
Re: Паншерпенинг Ландсата 8
Поясните, пож-та, недалёким (от грасса) людям, как.gimran писал(а):Мжно в кугисе калькулятором растров, экспериментировал для того, чтобы далеким от грасса людям обьяснить тупейший способ.
Редактор материалов, модератор форума
-
- Интересующийся
- Сообщения: 31
- Зарегистрирован: 23 сен 2013, 18:45
- Репутация: 7
Re: Паншерпенинг Ландсата 8
В 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.
Если приемлемо преобразование Бровея, то после допиливания 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
В GRASS никак нельзя преобразовать 16-бит в 8-бит?
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________
- gimran
- Гуру
- Сообщения: 1902
- Зарегистрирован: 07 июл 2010, 15:43
- Репутация: 242
- Откуда: Уфа
Re: Паншерпенинг Ландсата 8
DN_b1
DN_fused = ------------------------------- * DN_pan
DN_b1 + DN_b2 + DN_b3
Где b1, b2, b3 каналы rgb, pan - панхром
Это для первого канала.
Для второго и третьего тоже самое, только вместо b1 ставим b2 и b3 в числитель, соответственно.
А потом полученные растры в композит.
[ Сообщение с мобильного устройства ]
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
То есть первый = (b1 × b_pan)/(b1+b2+b3)
Второй = (b2 × b_pan)/(b1+b2+b3)
Третий = (b3 × b_pan)/(b1+b2+b3)
[ Сообщение с мобильного устройства ]
Второй = (b2 × b_pan)/(b1+b2+b3)
Третий = (b3 × b_pan)/(b1+b2+b3)
[ Сообщение с мобильного устройства ]
-
- Гуру
- Сообщения: 920
- Зарегистрирован: 30 дек 2008, 14:11
- Репутация: 236
- Откуда: Ханты-Мансийск
- Контактная информация:
Re: Паншерпенинг Ландсата 8
Можно: http://grass.osgeo.org/grass70/manuals/r.rescale.htmlrhot писал(а):В GRASS никак нельзя преобразовать 16-бит в 8-бит?
Но нужно разбираться как правильно настроить гистограмму.
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 7 гостей