То есть классификация landcover выступает аналогом NDVI.
Нет, цель максимально корректно задать параметр emissivity, потому что ошибка в нем существенно сказывается на погрешности определения температура, приводятся числа порядка 0.4-0.9 градуса за 0.01 emissivity. Но наверное это и не так страшно, если ошибка будет, то близкая для всех и Ваша неоднородность не пострадает. Если Вам интересна только температура где есть растительность, а остальные территории не интересны, тогда да, использовать landcover бессмысленно и проще всего обойтись оценкой emissivity по NDVI или константой.
Тогда смотрите работы, которые используют NDVI для этой цели.
Например очень много ссылок на
Sobrino, J.A.; Jiménez-Muñoz, J.C.; Paolini, L. Land surface temperature retrieval from LandsatTM 5. Remote Sens. Environ. 2004, 90, 434–440.
Если взять код как это делается в GEOBIA in ArcGIS, то там все очень весело
Есть модель описывающая связь между NDVI и emissivity (Van De Griend and Ove, 1993)
Понятно, что все значения NDVI эта модель описывать не может и будет подходить только для растительности.
Далее на основе NDVI растр разбивается на 4 категории
Код: Выделить всё
-1.0,-0.18,1 класс
-0.18,0.145,2 класс
0.145,0.728,3 класс
0.728,1.0,4 класс
По сути тот же landcover, только урезанный
А далее еще веселее
Код: Выделить всё
1 класс, использует константу emissivity 0.985
2 класс, использует константу emissivity 0.919
3 класс, должен был использовать 1.0094+0.047*(Ln(NDVI)), но по факту
использует константу NDVI=0.5 и как результат emissivity 0.976
4 класс, использует константу emissivity 0.994
И как тогда измениться последняя формула?
Ваша задача получить растр emissivity и использовать его в калькуляторе
Код: Выделить всё
TS=(TB/[1+(10.8∗TB/14388)∗ln(emissivity ))-273.15
Я понимаю методику так.
Да, для расчета NDVI, каналы красный и ближний ИК. Перевод в radiance и атмосферная коррекция (DOS - проще, 6S-лучше). Можно сразу посмотреть Landsat Surface Reflectance Level-2 Science Products, ссылка на него уже была.
Для тепловых
сложнее.
Barsi_AtmCorr_SPIE05.pdf
И в цепочку перед расчетом TOA brightness temperature добавляется
where LT is the radiance of a blackbody target of kinetic temperature T;
τ is the atmospheric transmission;
ε is the emissivity of the surface, specific to the target type;
Lu is the upwelling or atmospheric path radiance;
Ld is the downwelling or sky radiance;
Lλ-TOA spectral radiance measured by the instrument (Watts/(m2*srad*µm))
Эти параметры скорее всего можно получить
тут.