Обсудить в форуме Комментариев 34Редактировать в вики
В статье описывается упрощённая технология коррекции многозональных снимков Landsat: компенсация радиометрической неоднородности режима съёмки, а также пересчёт значений теплового канала в температуру
Содержание |
Материалы Landsat в том виде, как они предоставляются пользователю, уже в значительной мере скорректированы, например, с учетом кривизны поверхности Земли, особенностей рельефа снимаемой территории, применено масштабирование в пределах возможных значений регистрируемых величин и т.п. Однако все эти коррекции направлены, главным образом, на обеспечение надёжного визуального дешифрирования наземных объектов.
Если же главным при анализе космоснимков становится не форма, рисунок и взаимное расположение объектов, а их индивидуальные спектральные характеристики, то требуется дополнительная коррекция изображений. Кроме того, тепловой канал, содержащийся в материалах Landsat, начиная с № 3, позволяет достаточно точно определить температуру подстилающей поверхности. Для этого также потребуется произвести определённые вычисления с исходными значениями.
Типичным примером такого анализа является вычисление вегетационных индексов по многозональным материалам, в первую очередь – нормализованного разностного вегетационного индекса NDVI. При интерпретации схем распределения NDVI крайне важными становятся не относительные изменения значений по одной сцене, а соотношения NDVI по одной сцене или разновременным сценам в абсолютном выражении.
Например, возьмём две сцены:
Эти сцены имеют взаимное перекрытие, поэтому можно сравнить, как выглядит на них один и тот же участок территории возле Кольской АЭС. Правомерно предположить, что за одни сутки (с погрешностью в 12 минут, судя по времени регистрации) в состоянии растительных сообществ не должно было произойти никаких значимых изменений. Но схемы NDVI, построенные для этих сцен, демонстрируют нечто совсем иное (рис. 1).
Сравнивая эти схемы, можно решить, что северная тайга в последний день июля 2004 года испытала быстрый расцвет. Судя по значениям NDVI, которые в среднем изменились на 18%, растительность повсеместно из категории угнетённой разреженной перешла в разряд пышной сомкнутой.
Этот курьёзный случай может быть объяснён только одним образом – нельзя сравнивать несравнимое. Прежде чем проводить вычисления NDVI, надо было выполнить коррекцию значений каналов каждой сцены в соответствии с калибровочными коэффициентами каждой из них.
Дело в том, что в данном случае сцены были засняты разными спутниками, аппаратура которых была настроена по-разному. Следовательно, зарегистрированные ими изображения без дополнительной радиометрической коррекции не годятся для тонкого спектрального анализа.
Дистанционное определение температуры поверхности Земли находит применение при геотермическом картировании, которое используется для решения широкого круга фундаментальных и прикладных задач:
Значения, содержащиеся в файлах каналов многозонального изображения, представляют собой безразмерную величину. Она пропорциональна интенсивности излучения, достигающего орбиты, на которой находится спутник Landsat[1]. То, что мы получаем из файла – дискретное калиброванное значение пикселя , – номинировано в условных относительных единицах DN (Digital Numbers – числовые значения).
Требуется получить из значения отражающей способности (альбедо)
наземных объектов, видимых на космоснимке.
Для материалов Landsat до 7 включительно, полученных из архива USGS, это производится следующим образом (1):
где – интенсивность излучения от объекта, достигшего орбиты Landsat,
– расстояние между Землёй и Солнцем,
– коэффициент светимости для каждого канала,
– высота стояния Солнца над горизонтом в момент съёмки.
Интенсивность излучения объекта, достигающего орбиты , вычисляется по формуле (2):
Калибровочные коэффициенты и
можно найти в файле мета-данных с именем «*_mtl.txt», который поставляется вместе со сценой в одном архивном файле. Из списка параметров, приведённых на рис. 2, должны быть выбраны RADIANCE_MULT (
) и RADIANCE_ADD (
), соответствующие нужному каналу.
Значение интенсивности излучения на сенсоре имеет размерность Вт/м²/ср/мкм, т.е. мощность излучения, падающая на единицу площади земной поверхности сквозь телесный угол в 1 стерадиан, взятого относительно единицы длины волны излучения.
Расстояние между Землёй и Солнцем может быть вычислено по приближённой формуле:
где – порядковый номер в году дня получения изображения. Его можно определить по дате регистрации сцены, которая также находится в файле метаданных (рис. 3):
Если для вычислений использовать MS Excel, то в нём для получения значения удобно использовать формулу
=Ref-ДАТА(ГОД(Ref);1;1)
, где Ref – ссылка на ячейку, где записана дата регистрации.
Величина измеряется в астрономических единицах (а.е.).
Значения светимости для каждого канала могут быть взяты из следующей таблицы (табл. 1):
Табл. 1. Коэффициенты светимости для каналов Landsat, Вт/м²/мкм
Каналы | 1 | 2 | 3 | 4 | 5 | 7 | 8 |
---|---|---|---|---|---|---|---|
Landsat 4 | 1957 | 1825 | 1557 | 1033 | 214.9 | 80.72 | – |
Landsat 5 | 1957 | 1826 | 1554 | 1036 | 215.0 | 80.67 | – |
Landsat 7 | 1970 | 1842 | 1547 | 1044 | 225.7 | 82.06 | 1369 |
И последний, необходимый для вычисления альбедо, параметр – высота стояния Солнца над горизонтом в момент съёмки – может быть взят из файла с метаданными сцены (рис. 4). Значение это в файле приведено в градусах. Перед подстановкой в формулу его надо перевести в радианы.
Для материалов Landsat 8, взятых из архива USGS, альбедо вычисляется проще:
Калибровочные коэффициенты и
для сцены Landsat 8 также можно найти в файле метаданных (рис. 5).
После вычисления по представленным формулам значения альбедо по разным каналам могут быть использованы, например, для определения величины NDVI:
Здесь NIR – номер канала ближнего инфракрасного диапазона, RED – номер канала красного диапазона. Определить, какие номера каналов соответствует этим диапазонам, можно с помощью табл. 2.
Табл. 2. Соответствие между диапазонами электромагнитного излучения и номерами каналов спутников разных поколений миссии Landsat[2][3][4][5][6]
Диапазон | Landsat 1 | Landsat 2 | Landsat 3 | Landsat 4 | Landsat 5 | Landsat 7 | Landsat 8 |
---|---|---|---|---|---|---|---|
RED | 5 | 5 | 5 | 3 | 3 | 3 | 4 |
NIR | 6 | 6 | 6 | 4 | 4 | 4 | 5 |
Тепловой | – | – | 8 | 6 | 6 | 61, 62 | 10, 11 |
По значениям тепловых каналов (табл. 2) можно определить радиояркостную температуру подстилающей поверхности. Теоретически точность оценки температуры около 0,5°С, однако, дымка в атмосфере занижает значения на несколько градусов.
Исходными данными для определения температуры служат значения интенсивности излучения , пришедшего на сенсор спутника и зарегистрированного соответствующим тепловым каналом. Величина
определяется в соответствии с методикой, описанной в предыдущем разделе (формула 2).
Для Landsat 4-7 температура вычисляется по следующей формуле:
Здесь и
– калибровочные константы, значения которых можно взять из табл. 3.
Табл. 3. Калибровочные константы Landsat 4-7 для вычисления температуры, Вт/м²/мкм
Диапазон | Landsat 4 | Landsat 5 | Landsat 7 |
---|---|---|---|
![]() |
671.62 | 607.76 | 666.09 |
![]() |
1284.30 | 1260.56 | 1282.71 |
Материалы Landsat 7 содержат два тепловых канала: 61 и 62, которые отличаются настройками усиления сигнала. Для вычисления температуры можно использовать в формуле (4) любой из них или оба – они дают почти одинаковые значения (±0,1°С).
Почти также[7] и с помощью той же формулы (4) производится вычисление температуры по данным Landsat 8, за исключением того, что значения констант и
надо брать отдельно для каждого теплового канала данной сцены из её файла метаданных (рис. 5).
Значения температуры, определённые по каналам 10 и 11 (они различаются охватываемыми интервалами теплового диапазона) Landsat 8, отличаются друг от друга на 1,5-3°С. В общем случае, можно усреднить их.
Пример результатов определения значений температуры подстилающей поверхности по данным теплового канала показан на рис. 6.
На нём видно, что городская территория уже хорошо прогревается тёплым весенним солнцем, в то время как за городом ещё лежит холодный снег. Одновременно выделяются интенсивные положительные аномалии от крупных заводов, ТЭЦ и очистных сооружений.
Константин Силкин 22:34, 24 сентября 2015 (MSD)
Обсудить в форуме Комментариев 34Редактировать в вики
Последнее обновление: 2015-12-15 05:49
Дата создания: 16.09.2015
Автор(ы): Константин Силкин
© GIS-Lab и авторы, 2002-2021. При использовании материалов сайта, ссылка на GIS-Lab и авторов обязательна. Содержание материалов - ответственность авторов. (подробнее).