Як зафіксувати LANDSAT 8 в GRASS?

9

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

У попередньому дописі хтось підтвердив це і запропонував використовувати іншу функцію для landsat: 'i.his.rgb' (після використання i.landsat.toar, яку я не використовую). Він починає конвертувати min / max до 1-255, перш ніж використовувати i.his.rgb ... Мої оригінальні растри мають діапазон 0-65535. Мені це не вийшло ...

Я кинув після декількох днів, що намагаються просто виконати паншарп, я відчуваю, що застряг ... Чому просто застосування i.pansharp не працює для Landsat, як пропонують інші? Чи може хтось запропонувати більш простий робочий процес для виконання панормотування зображень Landsat 8? З нетерпінням чекаю поради, заздалегідь дякую ...

— Арноуд Купенс
джерело

7

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

BTW: діапазон 0-65535 для вихідних даних LANDSAT є дивним, оскільки вони передаються у вигляді 8-бітових каналів (0-255). Тому спочатку перевірте свої вхідні дані.

Оновлення: 10/2013

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

— markusN
джерело

Я спочатку натрапив на цю статтю, з якою ви посилаєтесь з neteler.org, але інструкції взагалі не працювали (через формат растрової плаваючої точки, що генерується i.landsat.toar, тому я в кінцевому підсумку дотримувався поради звідси: gis.stackexchange.com / a / 39239/9857

— Хуан

3

Зауважте, що ця порада була надана багато років тому і не відображає останнього GIS-коду GRASS. BTW: по-справжньому чудовий модуль панормотування - це новий додаток: i.fusion.hpf - З'єднання панкроматичних даних із високою роздільною здатністю та низькою роздільною здатністю, засновані на техніці доповнення до фільтру високих частот, див. Trave.osgeo.org/grass70/manuals/ addons / i.fusion.hpf.html

— markusN

4

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

— Арноуд Купенс
джерело

3

Так, зображення Landsat 8 зберігаються як типи Uint16.

— Вебріан

3

Панель різкості FIHS або Brovery Transform є досить простою, якщо ви дотримуєтесь формул із такої статті: http://www.sensorsportal.com/HTML/DIGEST/august_2014/Vol_177/P_RP_0183.pdf

Якщо припустити, що ви використовуєте ГРІС ГІС та паніруючий Landsat 8.

Встановіть свій регіон, використовуючи панораму 8 tif - важливо мати роздільну здатність 15 м, яка буде встановлена, якщо ви це зробите.

Для Brovery без ваг введіть це в Калькулятор карт синього, зеленого та червоного кольорів відповідно:

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-смуговій. Якщо хтось цікавий - він використовує метод Brovey Transform, як описано в посібнику (див. Https://media.readthedocs.org/pdf/semiautomaticclassificationmanual/latest/semiautomaticclassificationmanual.pdf ).

— Марк Томпсон
джерело

1

Ласкаво просимо на сайт Janos. В ідеалі відповіді на кшталт вашої повинні мати опис того, чому вони є хорошими методами. ..

Для панормотування зображень Landsat 8 простішим способом є http://www.geosage.com/highview/download.html . Моя порада, спробуйте. Лише одним клацанням…. Спектралтрансформатор

— Янош Брезнай
джерело

-2

Для pansharpenig я використав спектральний трансформатор GUI для фотографій Landsat 8 введіть тут опис зображення http://www.geosage.com/highview/download.html Дуже добре. BR Янош

— Янош Брезнай
джерело

2

Ласкаво просимо на сайт Janos. В ідеалі відповіді на кшталт вашої повинні мати опис того, чому вони є хорошими методами. Відповіді лише на посилання стають непрацездатними, якщо посилання відійде.

— Аарон