Атмосферная коррекция данных Landsat/ETM+ (COST метод)
- Максим Дубинин
- MindingMyOwnBusiness
- Сообщения: 9129
- Зарегистрирован: 06 окт 2003, 20:20
- Репутация: 748
- Ваше звание: NextGIS
- Откуда: Москва
- Контактная информация:
Атмосферная коррекция данных Landsat/ETM+ (COST метод)
Обсуждение темы:
Атмосферная коррекция данных Landsat/ETM+ (COST метод)
http://gis-lab.info/qa/atcor.html
Атмосферная коррекция данных Landsat/ETM+ (COST метод)
http://gis-lab.info/qa/atcor.html
пристегивайтесь, турбулентность прямо по курсу
-
- Новоприбывший
- Сообщения: 1
- Зарегистрирован: 21 сен 2006, 10:31
- Репутация: 0
Атмосферная коррекция

И получил значения пикселов от 0 до 0.6.
Это что за значения?
Это значения уже спектральной отражательной способности или нет?
И как пересчитать эти значения в привычные коэффициенты спектральной яркости в ERDAS Imagine 8.6 для составления спектральных образов исследуемых объектов и для проведения классификации?
Может кто-нибуть подсказать.
-
- Завсегдатай
- Сообщения: 386
- Зарегистрирован: 07 фев 2004, 14:31
- Репутация: 7
- Откуда: Лозанна
- Контактная информация:
Обновлена модель artor, позволяющая провести атмосферную коррекцию. В COST-ETM7.xls добавлен ряд комментариев, поясняющих алгоритм.
toGidra
Значения могут быть и такими (чаще всего от 0 до 3).
Это конечно не DN - смотрите пояснения в самой статье. Пересчитать эти значения в "привычные коэффициенты спектральной яркости" можно применив алгорим в обратном порядке
toGidra
Значения могут быть и такими (чаще всего от 0 до 3).
Это конечно не DN - смотрите пояснения в самой статье. Пересчитать эти значения в "привычные коэффициенты спектральной яркости" можно применив алгорим в обратном порядке

-
- Новоприбывший
- Сообщения: 5
- Зарегистрирован: 06 дек 2006, 00:24
- Репутация: 0
- Откуда: геофак
- Контактная информация:
-
- Новоприбывший
- Сообщения: 5
- Зарегистрирован: 06 дек 2006, 00:24
- Репутация: 0
- Откуда: геофак
- Контактная информация:

и все(((((((((((.
Слышал что каналы Ландсат как-то корректируют на среднюю прозрачность атмосферы (коэф. 0.6) как раз на этом этапе (L1G)
но ссылку найти не могу конкретную.

...чувствовал себя дикарём среди миклухо-маклаев...
- Максим Дубинин
- MindingMyOwnBusiness
- Сообщения: 9129
- Зарегистрирован: 06 окт 2003, 20:20
- Репутация: 748
- Ваше звание: NextGIS
- Откуда: Москва
- Контактная информация:
Да, вы правы, атамосферную коррекцию L1G не проходят. Про среднюю прозрачность я даже не слышал и в Landsat Handbook про это ничего не сказано, поэтому я бы сказал, что этого для L1G скорее не делается.
А фраза эта относилась скорее к продуктам второго уровня обработки, но вы правы, лучше ее переписать, чтобы не возникало путаницы.
А фраза эта относилась скорее к продуктам второго уровня обработки, но вы правы, лучше ее переписать, чтобы не возникало путаницы.
пристегивайтесь, турбулентность прямо по курсу
-
- Новоприбывший
- Сообщения: 5
- Зарегистрирован: 06 дек 2006, 00:24
- Репутация: 0
- Откуда: геофак
- Контактная информация:
Спасибо, похоже я принял желаемое за действительное на уровне L1G. ничего атмосферного нет(((((, а про коффициент видимо это была какая-то частная работа, чтото в духе этого -
http://d902.iki.rssi.ru/theses-cgi/thesis.pl?id=559
http://d902.iki.rssi.ru/theses-cgi/thesis.pl?id=559
...чувствовал себя дикарём среди миклухо-маклаев...
-
- Интересующийся
- Сообщения: 33
- Зарегистрирован: 17 янв 2008, 16:14
- Репутация: 0
-
- Новоприбывший
- Сообщения: 10
- Зарегистрирован: 30 июл 2008, 18:09
- Репутация: 0
Пытаемся применить этот метод для снимков TERRA.
Не нашли значений для Lmin Lmax - это вроде количество приходящего излучения которое максимально м минимально может зафиксировать сенсор для каждого канала. не совсем понятно как это связано с длиной волны(если связано) и почему на сайте TERRA дается только одно значение для "Spectral Radiance" для каждого канала http://modis.gsfc.nasa.gov/about/specifications.php хотя вроде должно быть два. минимум)
может кто хоть ссылку даст на объяснение как это(атмосферную коррекцию) надо делать для TERRA
Не нашли значений для Lmin Lmax - это вроде количество приходящего излучения которое максимально м минимально может зафиксировать сенсор для каждого канала. не совсем понятно как это связано с длиной волны(если связано) и почему на сайте TERRA дается только одно значение для "Spectral Radiance" для каждого канала http://modis.gsfc.nasa.gov/about/specifications.php хотя вроде должно быть два. минимум)
может кто хоть ссылку даст на объяснение как это(атмосферную коррекцию) надо делать для TERRA
-
- Новоприбывший
- Сообщения: 3
- Зарегистрирован: 05 май 2009, 22:47
- Репутация: 0
- Откуда: Краснодар
Re: Атмосферная коррекция данных Landsat/ETM+ (COST метод)
Спасибо огромное за столь полезную статью! Снимки после данной атмосферной коррекции получаются очень четкими и яркими.
Некоторые дополнения к статье и к модели:
1. Здесь в расчетах присутствует не "зенит", а "зенитное расстояние".
Зенит - точка пересечения небесной сферы с отвесной линией над головой наблюдателя.
Зенитное расстояние - это дуга вертикала между зенитом и светилом.
<<В горизонтальной системе основным кругом служит математический, или истинный, горизонт NESW, полюсом — зенит Z места наблюдения. Для определения положения светила s проводят через него и Z большой круг, который называется кругом высоты, или вертикалом, данного светила. Дуга Zs вертикала от зенита до светила называется его зенитным расстоянием z и является первой координатой; z может иметь любое значение от 0° (для зенита Z) до 180° (для надира Z'). Вместо z пользуются также высотой светила h, равной дуге круга высоты от горизонта до светила. Высота отсчитывается в обе стороны от горизонта от 0° до 90° и считается положительной, если светило находится над горизонтом, и отрицательной — если светило под горизонтом. При таком условии всегда справедливо соотношение z + h = 90°>>.
http://www.cultinfo.ru/fulltext/1/001/008/080/663.htm
В статье зенитное расстояние z - это Q (SolarZenithAngle), высота светила h - это SUN_ELEVATION.
2. Входной снимок в модели должен из себя представлять единый, многоканальный файл (каналы 1, 2, 3, 4, 5, 7, объединенные с помощью команды Layer Stack).
3. Значения минимального значения DN (номер строки) были взяты без сортировки столбца Histogram.
4. Формулы из EXCEL нужно копировать в ERDAS без "MODEL GMD = ". Необходимо проверить, чтобы вставленные в формулы новые значения имели разделитель целой и дробной части числа в виде точки, а не запятой.
5. В графической модели после вставки каждой формулы необходимо заменить переменные $n1_phnx(1) и др. на переменные соответствующего канала для снимка. То есть, если корректируемый снимок называется Landsat-173-28, то в формуле для 1 канала переменная $n1_phnx(1) должна быть заменена на переменную $n1_Landsat-173-28(1), в формуле для 2 канала переменная $n1_phnx(2) должна быть заменена на переменную $n1_Landsat-173-28(2), и т.д.

Некоторые дополнения к статье и к модели:
1. Здесь в расчетах присутствует не "зенит", а "зенитное расстояние".
Зенит - точка пересечения небесной сферы с отвесной линией над головой наблюдателя.
Зенитное расстояние - это дуга вертикала между зенитом и светилом.
<<В горизонтальной системе основным кругом служит математический, или истинный, горизонт NESW, полюсом — зенит Z места наблюдения. Для определения положения светила s проводят через него и Z большой круг, который называется кругом высоты, или вертикалом, данного светила. Дуга Zs вертикала от зенита до светила называется его зенитным расстоянием z и является первой координатой; z может иметь любое значение от 0° (для зенита Z) до 180° (для надира Z'). Вместо z пользуются также высотой светила h, равной дуге круга высоты от горизонта до светила. Высота отсчитывается в обе стороны от горизонта от 0° до 90° и считается положительной, если светило находится над горизонтом, и отрицательной — если светило под горизонтом. При таком условии всегда справедливо соотношение z + h = 90°>>.
http://www.cultinfo.ru/fulltext/1/001/008/080/663.htm
В статье зенитное расстояние z - это Q (SolarZenithAngle), высота светила h - это SUN_ELEVATION.
2. Входной снимок в модели должен из себя представлять единый, многоканальный файл (каналы 1, 2, 3, 4, 5, 7, объединенные с помощью команды Layer Stack).
3. Значения минимального значения DN (номер строки) были взяты без сортировки столбца Histogram.
4. Формулы из EXCEL нужно копировать в ERDAS без "MODEL GMD = ". Необходимо проверить, чтобы вставленные в формулы новые значения имели разделитель целой и дробной части числа в виде точки, а не запятой.
5. В графической модели после вставки каждой формулы необходимо заменить переменные $n1_phnx(1) и др. на переменные соответствующего канала для снимка. То есть, если корректируемый снимок называется Landsat-173-28, то в формуле для 1 канала переменная $n1_phnx(1) должна быть заменена на переменную $n1_Landsat-173-28(1), в формуле для 2 канала переменная $n1_phnx(2) должна быть заменена на переменную $n1_Landsat-173-28(2), и т.д.
-
- Завсегдатай
- Сообщения: 410
- Зарегистрирован: 24 авг 2006, 17:35
- Репутация: 10
Re: Атмосферная коррекция данных Landsat/ETM+ (COST метод)
Народ, а мог бы кто-нить выложить примеры снимков после коррекции? Заранее спасибо!
В мета-файлах, которые идут в комплекте со снимками, выкачиваемыми через USGS GLOVIS, вся необходимая инфа для рассчета есть?
В мета-файлах, которые идут в комплекте со снимками, выкачиваемыми через USGS GLOVIS, вся необходимая инфа для рассчета есть?
-
- Гуру
- Сообщения: 534
- Зарегистрирован: 30 ноя 2006, 13:31
- Репутация: 116
- Откуда: Moscow
Re: Атмосферная коррекция данных Landsat/ETM+ (COST метод)
Да, вся информация есть в файле MTL.
-
- Участник
- Сообщения: 60
- Зарегистрирован: 12 май 2009, 13:03
- Репутация: 0
- Откуда: Минск
Re: Атмосферная коррекция данных Landsat/ETM+ (COST метод)
Есть ошибка в формуле в разделе 2 пункт е) - ESUN должен быть в числителе, d2 - в знаменателе.
-
- Новоприбывший
- Сообщения: 1
- Зарегистрирован: 14 янв 2010, 17:04
- Репутация: 0
- Контактная информация:
Re: Атмосферная коррекция данных Landsat/ETM+ (COST метод)
Помимо ошибки, найденной Aleksandr'om, есть еще одна, более фундаментальная ошибка в описании.
Если взглянуть на формулу по вычислению Value, то возникает вопрос о размерности в следующем месте формулы: (DN-Lhaze). Т.е. мы от DigitalNumbers, которые по сути просто градации яркости, вычитаем Lhaze, значение дымки размерностью [W m^(-2) sr^(-1) mkm^(-1)], что в корне неверно, это все равно, что мы из человек вычитаем лопаты. Поэтому остается открытым вопрос о размерности не только Value, но и о том, что это вообще за значения.
Если взглянуть на статью 1989г, то на стр. 4 в формуле (5) вместо DN находится RAD, т.е. показатели излучения на сенсоре. Так же это подтвержает и статья 1996г: стр.3, формула (2). По сути (2) - это тоже самое, что и (5) из 1989г.
Поэтому, я несколько не доверяю предложенной модели.
Возможно, я ошибаюсь, но то, что мы не можем из DN вычитать Lhaze, на мой взгляд - это факт.
Очень интересно было бы ознакомиться со статьей Chavez'а от 1975г, на которую ссылаются при описании dark-object subtraction, но я ее, увы, не нашла. Если у кого-то она есть, поделитесь, пожалуйста.
Если взглянуть на формулу по вычислению Value, то возникает вопрос о размерности в следующем месте формулы: (DN-Lhaze). Т.е. мы от DigitalNumbers, которые по сути просто градации яркости, вычитаем Lhaze, значение дымки размерностью [W m^(-2) sr^(-1) mkm^(-1)], что в корне неверно, это все равно, что мы из человек вычитаем лопаты. Поэтому остается открытым вопрос о размерности не только Value, но и о том, что это вообще за значения.
Если взглянуть на статью 1989г, то на стр. 4 в формуле (5) вместо DN находится RAD, т.е. показатели излучения на сенсоре. Так же это подтвержает и статья 1996г: стр.3, формула (2). По сути (2) - это тоже самое, что и (5) из 1989г.
Поэтому, я несколько не доверяю предложенной модели.
Возможно, я ошибаюсь, но то, что мы не можем из DN вычитать Lhaze, на мой взгляд - это факт.
Очень интересно было бы ознакомиться со статьей Chavez'а от 1975г, на которую ссылаются при описании dark-object subtraction, но я ее, увы, не нашла. Если у кого-то она есть, поделитесь, пожалуйста.
-
- Интересующийся
- Сообщения: 40
- Зарегистрирован: 19 янв 2012, 10:03
- Репутация: 3
Re: Атмосферная коррекция данных Landsat/ETM+ (COST метод)
Посмотрел статью, появились вопросы:
1.
- Полностью согласен. Те DN надо еще пересчитать в L по формуле в пункте д.
2. Подскажите, пожалуйста, а как полученные значения соотносятся с радиометрическим разрешением канала?
Интересует вопрос калибровки данных разных сенсоров. Или вообще не нужно смотреть на это разрешение, если, говорится, что мы получили абсолютные безразмерные значения отражения?
1.
barash писал(а):Если взглянуть на статью 1989г, то на стр. 4 в формуле (5) вместо DN находится RAD, т.е. показатели излучения на сенсоре.
- Полностью согласен. Те DN надо еще пересчитать в L по формуле в пункте д.
2. Подскажите, пожалуйста, а как полученные значения соотносятся с радиометрическим разрешением канала?
Интересует вопрос калибровки данных разных сенсоров. Или вообще не нужно смотреть на это разрешение, если, говорится, что мы получили абсолютные безразмерные значения отражения?
Последний раз редактировалось Cocodrilo 19 сен 2012, 14:53, всего редактировалось 1 раз.
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 1 гость