Land Surface Temperature plugin (Расчет NDVI, Surface Temperature из Landsat, ASTER)

[Natalia Novoselova]

Одним (наконец то давшим результат!) итогом моих поисков надежных автоматизированных методов расчета температуры поверхности по Landasat8 (по скачиваемому продукту Landsat8 OLI/TIRS C1 Level-1 )
Нашла Land Surface Temperature (LST) QGIS plugin
http://qgislstplugin.blogspot.com

сейчас в списке плагинов QGIS он обозначен как «экспериментальный»

Может быть, сообществу будет интересно на него посмотреть, дать свою оценку.
Как уже говорила в другой теме, одно из того, что намереваюсь сделать, провести верификацию данных со спутника, методов расчета и измерений температуры реальных. Интересно же сравнить разные методы оценки температуры по космоснимку и что в реальности показывает термометр.
Измерения планирую проводить в указанных точках обзора.


В отличие от других модулей – там все понятно и наглядно. Поэтапно вычисляешь все параметры. Расчет emissivity идет на основе NDVI.
Получилась вполне правдоподобная карта температуры. В центре – мой фрагмент тропического леса. Как и ожидалось, его температура намного ниже соседних полей, и внутри, о счастье, также есть неоднородность (которая мне и нужна).

Единственное, что мне пока не ясно. Сам этот плагин проводит атмосферную коррекцию? Или все-таки ее надо проводить отдельно? Для ландсатов (Последний шаг расчета "Land Surface Temperature Algorithm > во вкладке "Single Channel Algorithm (Landsat)") можно указать параметры атмосферы.

Это

Atmospheric Water Vapour
At Sensor Radiance


At Sensor Radiance – вроде бы это то, что считается в первом шаге (т.е. полученный растр “Radiance”).

Atmospheric Water Vapour – наверное, можно узнать для своей территории из данных метеонаблюдений и внести.

И мой вопрос.

Весь алгоритм расчета плагина полный? Если внести «Atmospheric Water Vapour» то это будет атмосферная коррекция? Или ее нужно делать еще отдельно для повышения точности оценки температуры земной поверхности?

[alexandr cherepanov]

Сам этот плагин проводит атмосферную коррекцию? Или все-таки ее надо проводить отдельно?

Все необходимое для обработки кроме самих значений параметров есть.

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

The Mono Window Algorithm (MWA)
The Radiative Transfer Equation (RTE)
The Single Channel Algorithm (SCA)

Все эти алгоритмы с учетом влияния атмосферы

[Natalia Novoselova]

Большое спасибо! Не знала, что атмосферные параметры можно так взять!

Было случайно удалено. Вот это посоветовали для расчета атмосферных параметров

Atmospheric Correction Parameter Calculator
https://atmcorr.gsfc.nasa.gov

Получается, что в данной теме сейчас собрано все необходимое, чтобы полу-автоматически рассчитать температуру поверхности из Landsat-8 с учетом атмосферной коррекции. QGIS и его сообщество как всегда, оказались на высоте.. :

Перенесла ваш дополняющий ответ из другой темы по этому плагину.

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

LT=(Lλ-Lu-τ(1−ε)Ld)/τε

Это The Radiative Transfer Equation в плагине.
Первый вариант без атмосферной коррекции тепла - The Planck Function.
NDVI считается по тому что ему дадите на вход, сам данные не готовит. Будут сырые - посчитает по сырым как есть.
С автоматизацией могут быть проблемы, по скольку параметры нужно находить и задавать.

Важные пояснения с этого поста: viewtopic.php?t=24286&start=15#p164249

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

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

[ericsson]

необходимо понимать теорию

Всегда. О чем вам и твердят уже не первый год.

[Natalia Novoselova]

Здесь про атмосферную коррекцию, которую можно использовать:

Как выбрать параметры для атмосферной коррекции Landsat-8 в калькуляторе
viewtopic.php?f=20&t=24374