Паншерпенинг Ландсата 8
Добавлено: 26 ноя 2013, 15:11
У кого-нибудь получилось в GRASS сделать более или менее сносный паршерпенинг снимков Landsat 8?
Форумы GIS-Lab.info
Геоинформационные системы (ГИС) и Дистанционное зондирование Земли
https://gis-lab.info/forum/
Паншерпенинг Ландсата 8
Страница 1 из 4
Re: Паншерпенинг Ландсата 8
Добавлено: 26 ноя 2013, 15:35
Да, с пом. какого-то стандартного алгоритма, сейчас, сходу не могу найти куда записал команду
Re: Паншерпенинг Ландсата 8
Добавлено: 26 ноя 2013, 15:37
Снимки на вход были после атмосферной коррекции? Там три алгоритма: brovey, PCA, IHS.
Re: Паншерпенинг Ландсата 8
Добавлено: 26 ноя 2013, 17:07
Использовал http://grass.osgeo.org/grass70/manuals/ ... arpen.html Метод IHS.rhot писал(а):Снимки на вход были после атмосферной коррекции? Там три алгоритма: brovey, PCA, IHS.
Без коррекции, т.к. это все равно только для визуадизации можно использовать.
Re: Паншерпенинг Ландсата 8
Добавлено: 26 ноя 2013, 17:12
Там только для начала нужно прогнать Landsat 8 через GDAL, чтобы получить 8bit.
Re: Паншерпенинг Ландсата 8
Добавлено: 26 ноя 2013, 17:15
Странно...
Без коррекции вообще пустые растры выдаёт.
Версия 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 Re: Паншерпенинг Ландсата 8
Добавлено: 26 ноя 2013, 17:22
Не подскажите команду?bolotoved писал(а):Там только для начала нужно прогнать Landsat 8 через GDAL, чтобы получить 8bit.
Преобразование 16-bit d 8 bit b паншарп Landsat-8
Добавлено: 26 ноя 2013, 18:01
Сорри, похоже, что я все делал с пом. 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.tif2) Склеиваем RGB-композит (4, 5 и 6 каналы в моем случае)
3)Теперь сам паншарп по 8-му:
Код: Выделить всё
gdal_landsat_pansharp -rgb LC81670152013184LGN00_8bit.TIF -pan LC81670152013184LGN00_B8_8bit.tif -o out3.tifRe: Паншерпенинг Ландсата 8
Добавлено: 26 ноя 2013, 18:50
Мжно в кугисе калькулятором растров, экспериментировал для того, чтобы далеким от грасса людям обьяснить тупейший способ.
[ Сообщение с мобильного устройства ]
[ Сообщение с мобильного устройства ]
Re: Паншерпенинг Ландсата 8
Добавлено: 26 ноя 2013, 19:09
Поясните, пож-та, недалёким (от грасса) людям, как.gimran писал(а):Мжно в кугисе калькулятором растров, экспериментировал для того, чтобы далеким от грасса людям обьяснить тупейший способ.
Re: Паншерпенинг Ландсата 8
Добавлено: 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.
Если приемлемо преобразование Бровея, то после допиливания gdal_calc.py по рецепту Voltrona (http://blogs.gis-lab.info/apple-touch-i ... /?paged=45), этот метод можно также применять в GDAL вместе с dans-gdal-scripts.
Re: Паншерпенинг Ландсата 8
Добавлено: 26 ноя 2013, 19:20
В GRASS никак нельзя преобразовать 16-бит в 8-бит?
Re: Паншерпенинг Ландсата 8
Добавлено: 26 ноя 2013, 19:32
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 в числитель, соответственно.
А потом полученные растры в композит.
[ Сообщение с мобильного устройства ]
Re: Паншерпенинг Ландсата 8
Добавлено: 26 ноя 2013, 19:34
То есть первый = (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)
[ Сообщение с мобильного устройства ]
Re: Паншерпенинг Ландсата 8
Добавлено: 27 ноя 2013, 01:37
Можно: http://grass.osgeo.org/grass70/manuals/r.rescale.htmlrhot писал(а):В GRASS никак нельзя преобразовать 16-бит в 8-бит?
Но нужно разбираться как правильно настроить гистограмму.