Альбедо по снимкам Landsat в ENVI

[А.А.Карасев]

Добрый день. Подскажите пожалуйста, как рассчитать альбедо по снимкам Landsat в ENVI?
Я обращался с этим вопросом к дистрибьюторам ENVI и получил такой ответ: "Для перевода значений пикселей DN съемочной системы Landsat 5 в единицы отражения в ПК ENVI Вам необходимо воспользоваться инструментом LandsatCalibration. Программа позволяет переводить значения DN в единицы яркости (Radiance) и в единицы отражения (Reflectance). Вам необходим перевод значений в единицы отражения Reflectance. Также целесообразно для получения корректных данных после перевода значений в единицы отражения выполнить атмосферную коррекцию."
Значения, которые у меня получились - достигают 3 единиц. Выходит что это не может быть альбедо в классическом понимании этого термина.
Подскажите, как рассчитать альбедо, чтоб показатели были от 0 до 1 (или до 100%)?

[nadiopt]

так а вы все сделали как написано и правильно?

[А.А.Карасев]

так а вы все сделали как написано и правильно?

Да, я сделал как было описано, но показатели достигают значения 3.

[alexandr cherepanov]

Версия 5.0.2 и выше. Classic
File -> Open External File -> Landsat - > GeoTIFF with Metadata указываете свой *_MTL.txt
Basic Tools -> Preprocessing -> Calibration Utilities -> Landsat Calibration выбираете свой *_MTL.txt ( тот что с 6 каналами, для пятого ландсата он должен быть вторым в списке)
Все параметры в окне должны выставиться на автомате. Reflectance стоит по умолчанию.

[gamm]

я сделал как было описано, но показатели достигают значения 3.

непонятно, что конкретно вы делали - посмотрите, какие значения у вас получатся без атмосферной коррекции. Если опять больше 3, то неверно пересчитаны DN, если нет, то что-то неладно с атмосферной коррекцией

[А.А.Карасев]

Версия 5.0.2 и выше. Classic
File -> Open External File -> Landsat - > GeoTIFF with Metadata указываете свой *_MTL.txt
Basic Tools -> Preprocessing -> Calibration Utilities -> Landsat Calibration выбираете свой *_MTL.txt ( тот что с 6 каналами, для пятого ландсата он должен быть вторым в списке)
Все параметры в окне должны выставиться на автомате. Reflectance стоит по умолчанию.

Я проделываю Ваш алгоритм, визуализирую конечный результат, и вижу, что в 3 канале у меня reflectance 2.78.
С сайта Экселис взято, что: "Reflectance is the proportion of the radiation striking a surface to the radiation reflected off of it".
Из этого понятия следует, что показатель reflectance противоположен альбедо (где показано отношение отраженной к приходящей и естественно оно не может быть больше 1).
Будет ли это правильным что после подсчетов reflectance я в математике каналов пропишу: 1/b1-7 и получу альбедо?

Хотел бы еще добавить что показатель reflectance, к примеру для лесных территорий меньше 1 во 2 канале, т.е. выходит что исходящей радиации этого спектра больше чем приходящей? Я еще понимаю когда такое может быть в тепловом канале, но здесь...

[alexandr cherepanov]

Нет, такого быть не может. Название сцены напишите.

[А.А.Карасев]

LT51760271986189XXX04

[alexandr cherepanov]

Исходник:

Код: Выделить всё

Filename: F:\LT51760271986189XXX04\LT51760271986189XXX04_MTL.txt
Dims: Full Scene (59,572,051 points)

Basic Stats    Min    Max       Mean      Stdev
     Band 1      0    255  62.769479  46.773652
     Band 2      0    254  26.604165  20.100512
     Band 3      0    255  24.373000  19.557188
     Band 4      0    255  40.070921  35.521895
     Band 5      0    255  48.645118  43.286436
     Band 6      0    255  22.825697  21.525488

Reflectance

Код: Выделить всё

Filename: [Memory2] (8071x7381x6)
Dims: Full Scene (59,572,051 points)

Basic Stats       Min       Max      Mean     Stdev
     Band 1  0.000000  0.333216  0.080279  0.059883
     Band 2  0.000000  0.722058  0.070690  0.053667
     Band 3  0.000000  0.672009  0.061108  0.049524
     Band 4  0.000000  0.838100  0.127242  0.114460
     Band 5  0.000000  0.536718  0.098387  0.089188
     Band 6  0.000000  0.773165  0.063550  0.062440

Все там нормально.

[А.А.Карасев]

Уважаемый Александр Сергеевич (если я правильно понял), огромное Вам спасибо за помощь. Я действительно неверно рассчитал и отнял у Вас время. Еще раз спасибо за отзывчивость!

[gornak]

Объясните тогда пожалуйста, кто разобрался, такой момент:
сцена LC81810202013154LGN00
калибровочные коэффициенты из
_MTL.txt
.........
GROUP = MIN_MAX_RADIANCE
RADIANCE_MAXIMUM_BAND_1 = 738.68207
RADIANCE_MINIMUM_BAND_1 = -61.00061

.........

GROUP = MIN_MAX_REFLECTANCE
REFLECTANCE_MAXIMUM_BAND_1 = 1.210700
REFLECTANCE_MINIMUM_BAND_1 = -0.099980
........
(примерно то же во всех каналах для всех сцен данной местности)
т.е. для откалиброванных данных априори возможны отрицательные значения (что и получается- использую Arcgis калькулятор растра- ввожу формулу из handbook -а), а также значения reflectance>1,
каков же тогда физический смысл таких значений и параметра reflectance вообще