Конвертация данных Landsat в значения температуры - Практика

[Максим Дубинин]

Обсуждение статьи "Конвертация данных Landsat TM/ETM+ в значения температуры - Практика"

http://gis-lab.info/qa/dn2temperature-sa.html

[Максим Дубинин]

Уточнены формулы, добавлено про разницу в Qcal для LPGS и NLAPS сцен.

[jeny]

Подскажите, пожалуйста, а в этих двух формулах не может быть ошибки:
Для High Gain: (1282.71/((666.09/ ((((12.65-3.2)/Qcal) *[Value]) +3.2) )+1 ).ln) - 273
Для Low Gain: (1282.71/((666.09/ ((((17.04)/Qcal) *[Value])))+1 ).ln) - 273?

Ведь если Qcal - калиброванное значение (DN), то из формулы
L=Lmax-Lmin/(Qcalmax-Qcalmin) * (Qcal-Qcalmin) + Lmin
в случае предобработки LPGS (Qcalmin=1), вроде бы, должна получиться формула
(1282.71/((666.09/ ((((12.65-3.2)/254) *([Value]-1)) +3.2) )+1 ).ln) – 273
(1282.71/((666.09/ ((((17.04)/254) *([Value]-1))))+1 ).ln) – 273
или это я сама где-то ошиблась?

[alexandr cherepanov]

да все правильно, для (Qcalmin=1) - [Value]-1
, для (Qcalmin=0) - [Value]

[kog9]

Вопрос. Пытаюсь то же самое сделать в ArcMap, выкидывает ошибки.
Сначала деление на ноль (исправляется удалением строки с нулевым значением), после выкидывает ошибку "требуется объект с непонятным номером". В чем может быть проблема?

[novia]

скриншот в студию, инструмент(параметры) и данные, с которыми запускаете.

[AlexKostin]

Вот какая интересная штука:
В этих двух формулах есть .ln.
Для High Gain: (1282.71/((666.09/ ((((12.65-3.2)/Qcal) *[Value]) +3.2) )+1 ).ln) - 273
Для Low Gain: (1282.71/((666.09/ ((((17.04)/Qcal) *[Value])))+1 ).ln) - 273
Из-за этой компоненты ArcGIS 10 не производит вычислений, а выдает сообщение об ошибке. Убираю из формулы ".ln" - вычисление происходит (понятно, что вычисляется что-то не то).
Мне кажется, что тут какая-то ошибка в написании формулы или может я чего-то недопонимаю. Какой-то кусочек этой строки логарифмируется. Разъясните если можно, больно уж интересная тема "температурные поля".

Перепробовал по разному, не вычисляется ничего. Наличие вот такого в формуле ".ln" не работает ни в ArcView 3.2 ни в ArcGIS 10. Кто-то может сам на практике пробовал? Или есть ссылка на первоисточник?
В ссылке на файл http://earthsys.ag.ohio-state.edu/nr797/Exercise3.doc нарисована таже формула почемуто. Неужели у ее автора она сработала - странно.

[aspirant]

Хоть у кого-то эта штука сработала?

[gamm]

вот тут все написано: http://serc.carleton.edu/NAGTWorkshops/ ... 48433.html

детали тут:
http://serc.carleton.edu/files/NAGTWork ... ing_te.pdf

и вообще http://serc.carleton.edu - знатный ресурс в части ГИС-преподавания и обучения ...