Все методы – абсолютно различны. Вас, как ученого, это не смущает? Где именно на форуме это было объяснено? Не было ничего и нигде. И никогда. Я использовала в прошлом проекте AST_08, поскольку тогда не смогла разобраться с Landsat и понять, как по нему рассчитать температуру.
Количество методов обусловлено тем, что значения, которые показывает тепловой канал (после пересчета - brightness temperature) и температура земной поверхности (surface temperature) связаны весьма косвенно. Представьте - у Вас нет никакой информации о состояние поверхности и атмосферы, есть только излучение зафиксированное съемочной аппаратурой и Вы должны определить как это излучение исказилось около земной поверхности и пройдя через атмосферу, прежде чем было зафиксировано аппратурой. Отсюда и множество подходов и методов и подключение дополнительных источников информации о состоянии поверхности и атмосферы. Одни более простые и требуют меньше сторонних данных, другие более сложные.
Качеством получаемых результатов (surface temperature) в абсолютных значениях тоже далеки от заявленных 0.1-1 градуса чувствительности аппаратуры. Лет 5 назад Вы поднимали эту тему.
Какая-то классификация поверхности там используется. Бред какой-то.
Растр landcover используется для задания излучательной способности (emissivity)
what-is-emissivity-and-why-is-it-important
7-8 листы pdf
В инструменте GEOBIA in ArcGIS, который Вы приводите в соседней теме, для оценки emissivity используется NDVI
12 лист pdf .
Статьи, которые без софта вообще бесполезны.
Как раз наоборот.
Скажите честно – ВЫ умеете строить температуру по ландсату?
Задачи определения surface temperature у меня не возникало. Для тех задач что были (лет 10 назад) хватало brightness temperature. Проблем с brightness temperature нет никаких, все решается 2 формулами. Да, можно сказать,что Я не умею строить температуру поверхности по Ландсат (хотя аппарат тут значения не имеет). С surface temperature чувствуется элемент шаманизма, поэтому не интересно пробовать
.
Упрощенный порядок действий на любом растровом калькуляторе:
1)Берете любой из тепловых каналов, например 10.
2)Находите в метаданных параметр RADIANCE_MULT_BAND_10 и RADIANCE_ADD_BAND_10, обычно эти параметры равны 3.3420E-04 и 0.10000 и считаете по формуле:
Код: Выделить всё
Lλ=RADIANCE_MULT_BAND∗BAND10+RADIANCE_ADD_BAND
where:
Lλ-TOA spectral radiance (Watts/(m2*srad*µm))
3)Находите в метаданных параметр K1_CONSTANT_BAND_10 и K2_CONSTANT_BAND_10, обычно эти параметры равны 774.8853 и 1321.0789.Теперь Вы можете посчитать TOA brightness temperature по формуле:
where:
TB-TOA brightness temperature (K)
Lλ-TOA spectral radiance (Watts/(m2*srad*µm))
K1-Band-specific thermal conversion constant from the metadata (K1_CONSTANT_BAND_x, where x is the thermal band number)
K2-Band-specific thermal conversion constant from the metadata (K2_CONSTANT_BAND_x, where x is the thermal band number)
Переход от DN к TOA brightness temperature делается только
так, никаких других вариантов быть не может. Для уменьшения числа операций можно объединить формулы
Код: Выделить всё
TB=K2/ln[(K1/(RADIANCE_MULT_BAND∗BAND10+RADIANCE_ADD_BAND))+1]
4)Готовимся к переходу к surface temperature и считаем по формуле:
where:
λ-wavelength of emitted radiance, для 10 канала - значение 10.8
c2=h∗c/s=1.4388∗10−2 m K = 14388 µm K
h = Planck’s constant = 6.626∗10−34 J s
s = Boltzmann constant = 1.38∗10−23 J/K
c = velocity of light = 2.998∗108 m/s
Подготовка состоит в том, чтобы определить, что Вы будете использовать в качестве значения emissivity - переклассифированный landcover, значение на основе NDVI, константу (если например интересующий вас участок однороден).
Если emissivity вводите в виде константы, например равно 0.98, и добавляем перевод в градусы Цельсия, то формула получается:
where:
TB-TOA brightness temperature (K)
TS-surface temperature (C)