Использование Landsat 7 для перехода от NDVI к коэффициенту излучения

Question

Я использую Landsat 7 для расчета температуры поверхности земли.
Я понимаю основные понятия, лежащие в основе преобразования, однако я смущен тем, как влиять на коэффициент излучения в мою модель.
Я использую построитель моделей для своих расчетов и создал несколько модулей, в которых используются корректирующие переменные Gain, Bias Offset, Landsat K1 и Landsat K2.
Я преобразовал DN в значения сияния.

Теперь мне нужно учитывать последнюю и, вероятно, самую запутанную (для меня) часть: Коэффициент излучения.
Я бы хотел рассчитать коэффициент излучения с помощью NDVI.
У меня есть модельная процедура, построенная для расчета уровня NDVI (band4-band3)/(band4 + band3).
Я также вычислил Pv, который представляет собой долю растительности, рассчитанную по: [NDVI - NDVI_min]/[NDVI_max-NDVI_min]^2.
Теперь, используя метод растительного покрытия, мне нужно только Ev и Eg.
Я не понимаю, как найти эти значения для вычисления значения общей излучательной способности на ячейку.

Кто-нибудь имеет представление о том, как я могу включить коэффициент излучения в мою формулировку?
Я немного запутался в том, как получить это значение ...

Answer 1

Я считаю, что коэффициент излучения часто включается как часть набора данных. В качестве альтернативы базы данных коэффициентов выбросов существуют (например, база данных ASTER здесь: https://lpdaac.usgs.gov/about/news_archive/aster_global_emissivity_database_ged_product_release и другие обычно поддерживаются академическими отделами.)

Значения Ev = 0.99 и Eg = 0.97 используются, а метод выбора обсуждается на p , 436 здесь: ftp://atmosfera.cl/pub/elias/Paula/2004_Sobrino_RSE.pdf (J.A. Sobrino et al. Поиск температуры поверхности земли от LANDSAT TM 5, дистанционное зондирование окружающей среды 90, 2004, стр. 434-440).

Другой подход взят здесь: http://fromgistors.blogspot.com/2014/01/estimation-of-land-surface-temperature.html

Оценка поверхности земли температура Есть несколько исследований о расчете температуры поверхности земли. Например, используя NDVI для оценки эмиссионной способности поверхности земли (Sobrino, et al., 2004) или используя классификацию наземного покрова для определения коэффициента излучения поверхности земли для каждого класса (Weng, et al., 2004).

Например, значения излучательной способности (e) различных типов растительного покрова приведены в следующей таблице (от Mallick, et al., 2012).

 
Soil: 0.928 
Grass: 0.982 
Asphalt: 0.942 
Concrete: 0.937 

Таким образом, температура поверхности земли может быть рассчитана как (Weng, и др 2004).:

Т = ТБ/[1 + (? * ТБ /) Lne]
где:
? = длина волны излучаемого излучения
? = Н * с/с (1,438 * 10^-2 м К)
ч = постоянная Планка (6,626 * 10^-34 Js)
з = постоянная Больцмана (1,38 * 10^-23 Дж/К)
c = скорость света (2.998 * 10^8 м/с)

Значения?для тепловых полос Landsat setellites перечислены в следующей таблице:

 
|  Satellite   | Band | Center wavelength (µm) | 
| Landsat 4, 5, and 7  | 6 | 11.45     | 
| Landsat 8    | 10 | 10.8     | 
| Landsat 8    | 11 | 12      | 

Дальнейшее чтение по выбору теплоотдачей, смотрите раздел 2.3, Лучеиспускаемость Retrieval, здесь: https://books.google.com/books?id=XN4uAYlexnsC&lpg=PA51&ots=YQrmDa2S1G&dq=vegetation%20and%20bare%20soil%20emissivity&pg=PA50#v=onepage&q&f=false