Конвертация данных Landsat TM/ETM+ в значения температуры - Практика
по адресу http://gis-lab.info/qa/dn2temperature-sa.html
Пример того как можно провести пересчет для данных ETM+ с использованием Arcview Spatial Analyst.
Откройте новый проект ArcView. Подключите два дополнительных модуля – Image Analysis и Spatial Analyst («File → Extensions») и установите рабочий каталог в меню Вида («File → Set Working Directory»). Это важно сделать вначале, иначе результаты обработки будут по умолчанию сохраняться в папку c:\windows\temp.
2. Загрузите оба тепловых канала (полученных в режиме работы high gain и low gain) в Вид ArcView (при загрузке Data source type необходимо указать Image Analyst Data Source)
3. Поочередно сконвертируйте изображения в гриды. Для этого в главном меню Вида выберите «Theme → Save Image As» и укажите List Files of Type – GRID. Когда AV предложит добавить созданный только что грид в Вид, необходимо нажать No (так как, по умолчанию он будет загружать новую тему не как грид, а как image analyst изображение).
4. Теперь необходимо загрузить созданные гриды в Вид. Для это в главном меню выбираете «View → Add Theme». Не забудьте указать тип загружаемых данных как Grid Data source.
5. Выберите одно из конвертированных изображений, откройте его атрибутивную таблицу, включите режим редактирования таблицы («Table → Start Editing»), добавьте новое поле («Edit → Add Field»). Укажите поле Name – Temperature, Type – Number, Width – 16, Decimal Place- 2. Нажмите OK.
6. Далее выберите кнопку Field Calculator и наберите в нем такое выражение:
Для 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 = Qcalmax - Qcalmin = 254, если Qcalmin = 1 (предобработка LPGS)
Qcal = Qcalmax - Qcalmin = 255, если Qcalmin = 0 (предобработка NLAPS).
Значения параметров необходимо уточнить в метаданных сцены. Подробнее о разнице между LPGS и NLAPS.
*Если съемка произведена в период с 24.11.1999 по 8.12.1999, то результаты вычисления для High Gain и Low Gain необходимо поменять местами, т.е. High Gain в этот период = Low Gain, и наоборот.
7. Нажмите ОК. Таким образом, в новом столбце Temperature значения DN будут пересчитаны в значения температуры в градусах Цельсия.
8. Сохраните результат «Table → Stop Editing»
9. Теперь в Виде в редакторе легенды (Legend Editor) укажите новое поле классификации (Classification Field) – Temperature. Нажмите применить - Apply.
Обратите внимание, что
- Эта формула применима для данных полученных со спутника Landsat-7, после 20 декабря 2000 года
- Результатом будут являться значения температуры, выраженные в градусах Цельсия. Чтобы получить значения в градусах Кельвина - не пишите в конце выражения (–273).
Подробную информацию о том, какие параметры необходимы для конвертации данных Landsat 4 и 5 в значения температур, можно найти здесь.