Landsat8, температура поверхности через ENVI 5.1 \ ArcGIS

[Natalia Novoselova]

Необходимо разобраться в вопросе - какой алгоритм создания карты температур поверхности по Landsat8 в ENVI
(в моем случае ENVI 5.1, вроде бы, уже поддерживающий новые мета-файлы ландсатов).

Поджимает время (нужно для проверки идей по моему будущему проекту на полевой практике, к которой сейчас готовлюсь), потому, буду очень признательна за рекомендации где искать такой алгоритм, может быть, уже технические статьи об этом. Или ваш опыт.

Когда сама научусь, весь алгоритм по пунктам - опишу в данной теме.

[Natalia Novoselova]

Прибавила в теме ENVI, ArcGIS

Просто я пока не владею иным софтом. Но тема, как мне кажется, общей необходимости. Если кому-то известны другие методы, можно тоже дать ссылки.

По итогам, как минимум, в теме - опишу алгоритм.

[Эдуард Казаков]

http://blogs.esri.com/esri/arcgis/2014/ ... ands-tirs/

http://geohackers.in/2013/08/using-data ... mperature/

http://www.youtube.com/watch?v=ST-tclrTMzc

[Эдуард Казаков]

Статейка
Land Surface Temperature Retrieval Methods From Landsat-8 Thermal Infrared Sensor Data

[Natalia Novoselova]

Спасибо, буду изучать. Ценные видео - показывают ряд аналогичных.
Да, вот так и надо учиться.

Склонясь теперь к использованию ArcGIS, или QGIS.. в общем, мне бы хоть как-то.

[Эдуард Казаков]

Вот ещё одна статья. Но это всё больше про концептуальные вещи, а не реализацию в конкретных продуктах.
Derivation of Land Surface Temperature for Landsat-8 TIRS Using a Split Window Algorithm

Когда-то решал задачу поиска тепловых аномалий с некоторыми условиями по восьмому ландсату, не имея в распоряжении никаких пакетов типа ENVI, писал сам на питоне с gdal, по подобным материалам. Но это была глубоко побочная задача, поэтому так толком и не получилось ничего, да и я быстро пришел к выводу, по получавшимся результатам, что без полевой верификации и калибровки это практически пустая трата времени, особенно если хочется получать "настоящий абсолютный цельсий".
В общем, сомневаюсь, что задача решается кнопками без отдельной изыскательской работы.

[Natalia Novoselova]

del

[Natalia Novoselova]

, особенно если хочется получать "настоящий абсолютный цельсий".
В общем, сомневаюсь, что задача решается кнопками без отдельной изыскательской работы.

Какого рода исзыскательской?

Мне так казалось, что достаточно точна должны быть методология, да и техника уже лучше, чем раньше.
В общем, моя цель использовать параметр температуры для выявления качества лесного покрова, в совокупности с NDVI. Может быть, еще с чем-то.
То есть, так думаю, если все территории обработать по одному методу, то сравнить их можно будет, сам чистый цельсий не так важен?

Ну так я пока думаю..

[Natalia Novoselova]

Использовала метод конвертирования Landsat8 в температуру, описанный в блоге ArcGIS:

http://blogs.esri.com/esri/arcgis/2014/ ... ands-tirs/

Я выполнила все шаги, но не могу понять, правильно ли я получила температуру на выходе?


Получается растр, значения которого явно не градусы Цельсий : 28122, 30472
Хотя похоже на цельсий умноженный на 1000
Однако про это в алгоритме не сказано (было бы, как мне кажется, там всё очень детально).

[Natalia Novoselova]

Вот, вложением - построенная по этому методу цепочка. Ее можно применить к любому Landsat-8 и попробовать перевести его в температуру. Или посмотреть - что там не так? Почему значения больше в 1000 раз от того, что похоже на температуру..


Также в этом методе не говорится о необходимости проводить атмосферную коррекцию. Я ее не проводила.

[nadiopt]

что за невероятное расширение файла?

[Natalia Novoselova]

что за невероятное расширение файла?

Это положено в zip.

Расширение делала, как там описано:

Step 5 – Applying this function chain to any Landsat 8 image

When you’re finished, save this as a raster function template.xml file. Then, open the Edit Raster Function geoprocessing tool. Select the Landsat 8 thermal raster product that you want to convert to degrees Fahrenheit as the input raster. Choose the Raster Function Template that you just saved, and leave the Function Name blank. Repeat this for as many Landsat 8 scenes as you want.


И на приведенной картинке (того блога о котором идет речь) под Шагом 5 - такой же файл (я сделала по этому образцу):
L8 thermal conversion to F.rft.xml

[ericsson]

А откуда у вас там в цепочке расширений взялось .rtf. ? Вы что, из Word или Wordpad это сохраняете?

[Natalia Novoselova]

А откуда у вас там в цепочке расширений взялось .rtf. ? Вы что, из Word или Wordpad это сохраняете?

Да, вы это правильно заметили. Как бы не в этом была причина.

Посмотрела: в примере в блоге файл называется ***.rft.xml

А у меня, при сохранении - автоматом дается расширение .rft
Хотя тип файла программа обозначает Raster Function Template XML Files

То есть она у меня сохраняет не правильно, хотя при этом - открывает сама этот файл, и использует (не выдает ошибку) Edit Raster Function Tool

А как это исправить?

[Natalia Novoselova]

del