GIS-LAB

Географические информационные системы и дистанционное зондирование


Форумы GIS-Lab.info

Геоинформационные системы (ГИС) и Дистанционное зондирование Земли


Атмосферная коррекция данных Landsat/ETM+ (COST метод)

Обсуждение материалов сайта: вопросы, замечания, предложения

Атмосферная коррекция данных Landsat/ETM+ (COST метод)

Сообщение Максим Дубинин » 07 мар 2006, 02:56

Обсуждение темы:
Атмосферная коррекция данных Landsat/ETM+ (COST метод)

http://gis-lab.info/qa/atcor.html
  • 2

пристегивайтесь, турбулентность прямо по курсу
Аватара пользователя
Максим Дубинин
NextGIS
 
Зарегистрирован: 06 окт 2003
Откуда: Москва (на карте)
 
Сообщения: 8855
Репутация: 613
Статьи: 231
Проекты: 12/6

Атмосферная коррекция

Сообщение Gidra » 21 сен 2006, 11:51

:D Сделал я атмосферную коррекцию для снимка ETM+/Landsat 7 по методу COST.
И получил значения пикселов от 0 до 0.6.
Это что за значения?
Это значения уже спектральной отражательной способности или нет?

И как пересчитать эти значения в привычные коэффициенты спектральной яркости в ERDAS Imagine 8.6 для составления спектральных образов исследуемых объектов и для проведения классификации?
Может кто-нибуть подсказать.
  • 0

Gidra
Новоприбывший
 
Зарегистрирован: 21 сен 2006
 
Сообщения: 1
Репутация: 0

Сообщение Анна » 26 сен 2006, 22:08

Обновлена модель artor, позволяющая провести атмосферную коррекцию. В COST-ETM7.xls добавлен ряд комментариев, поясняющих алгоритм.

toGidra
Значения могут быть и такими (чаще всего от 0 до 3).
Это конечно не DN - смотрите пояснения в самой статье. Пересчитать эти значения в "привычные коэффициенты спектральной яркости" можно применив алгорим в обратном порядке ;)
  • 0

Анна
Завсегдатай
 
Зарегистрирован: 07 фев 2004
Откуда: Лозанна (на карте)
 
Сообщения: 386
Репутация: 7
Статьи: 15
Проекты: 5

Сообщение Robert » 04 май 2007, 00:58

"В настоящее время, большая часть данных попадающих в руки конечного пользователя уже атмосферно скорректирована"
А где для ландсата можно почитать про енто? что они там с ним делают? шум удаляют или чего?
  • 0

...чувствовал себя дикарём среди миклухо-маклаев...
Robert
Новоприбывший
 
Зарегистрирован: 06 дек 2006
Откуда: геофак (на карте)
 
Сообщения: 5
Репутация: 0

Сообщение Robert » 08 май 2007, 02:28

:arrow: Я тут почитал документацию - про уровень обработки L1G - там вобще ни слова про атмосферную коррекцию, там только геометрия и радиометрия (радиометрическая - насколько я понимаю - это устраннеие недостатков при юстировке камеры, всякие там биения и прочая "сенсорика", геометрия - это типа они так привязку называют)
и все(((((((((((.
Слышал что каналы Ландсат как-то корректируют на среднюю прозрачность атмосферы (коэф. 0.6) как раз на этом этапе (L1G)
но ссылку найти не могу конкретную.
:( Может кто знает поточнее? а то ссылаться надо куда-то, не хочется людей обманывать.
  • 0

...чувствовал себя дикарём среди миклухо-маклаев...
Robert
Новоприбывший
 
Зарегистрирован: 06 дек 2006
Откуда: геофак (на карте)
 
Сообщения: 5
Репутация: 0

Сообщение Максим Дубинин » 09 май 2007, 17:39

Да, вы правы, атамосферную коррекцию L1G не проходят. Про среднюю прозрачность я даже не слышал и в Landsat Handbook про это ничего не сказано, поэтому я бы сказал, что этого для L1G скорее не делается.

А фраза эта относилась скорее к продуктам второго уровня обработки, но вы правы, лучше ее переписать, чтобы не возникало путаницы.
  • 0

пристегивайтесь, турбулентность прямо по курсу
Аватара пользователя
Максим Дубинин
NextGIS
 
Зарегистрирован: 06 окт 2003
Откуда: Москва (на карте)
 
Сообщения: 8855
Репутация: 613
Статьи: 231
Проекты: 12/6

Сообщение Robert » 18 май 2007, 17:16

Спасибо, похоже я принял желаемое за действительное на уровне L1G. ничего атмосферного нет(((((, а про коффициент видимо это была какая-то частная работа, чтото в духе этого -
http://d902.iki.rssi.ru/theses-cgi/thesis.pl?id=559
  • 0

...чувствовал себя дикарём среди миклухо-маклаев...
Robert
Новоприбывший
 
Зарегистрирован: 06 дек 2006
Откуда: геофак (на карте)
 
Сообщения: 5
Репутация: 0

Сообщение tim_orlov » 17 янв 2008, 17:52

Два вопроса:
1. Надо ли сортировать гистограмму при выборе минимального значение DN (digital number).
2. Что надо вводить дальше номер строки? или соответсвеующее тому номеру значение value?
  • 0

tim_orlov
Интересующийся
 
Зарегистрирован: 17 янв 2008
 
Сообщения: 31
Репутация: 0

Сообщение Programmer » 01 авг 2008, 18:25

Пытаемся применить этот метод для снимков TERRA.
Не нашли значений для Lmin Lmax - это вроде количество приходящего излучения которое максимально м минимально может зафиксировать сенсор для каждого канала. не совсем понятно как это связано с длиной волны(если связано) и почему на сайте TERRA дается только одно значение для "Spectral Radiance" для каждого канала http://modis.gsfc.nasa.gov/about/specifications.php хотя вроде должно быть два. минимум)

может кто хоть ссылку даст на объяснение как это(атмосферную коррекцию) надо делать для TERRA
  • 0

Programmer
Новоприбывший
 
Зарегистрирован: 30 июл 2008
 
Сообщения: 10
Репутация: 0

Re: Атмосферная коррекция данных Landsat/ETM+ (COST метод)

Сообщение Svetla_GIS » 08 май 2009, 00:17

Спасибо огромное за столь полезную статью! Снимки после данной атмосферной коррекции получаются очень четкими и яркими. :D

Некоторые дополнения к статье и к модели:
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), и т.д.
  • 0

Svetla_GIS
Новоприбывший
 
Зарегистрирован: 05 май 2009
Откуда: Краснодар (на карте)
 
Сообщения: 3
Репутация: 0

Re: Атмосферная коррекция данных Landsat/ETM+ (COST метод)

Сообщение manisha » 25 окт 2009, 23:23

Народ, а мог бы кто-нить выложить примеры снимков после коррекции? Заранее спасибо!
В мета-файлах, которые идут в комплекте со снимками, выкачиваемыми через USGS GLOVIS, вся необходимая инфа для рассчета есть?
  • 0

manisha
Завсегдатай
 
Зарегистрирован: 24 авг 2006
 
Сообщения: 410
Репутация: 10
Статьи: 2
Проекты: 1

Re: Атмосферная коррекция данных Landsat/ETM+ (COST метод)

Сообщение alexandr cherepanov » 26 окт 2009, 11:18

Да, вся информация есть в файле MTL.
  • 0

alexandr cherepanov
Завсегдатай
 
Зарегистрирован: 30 ноя 2006
Откуда: Moscow (на карте)
 
Сообщения: 469
Репутация: 59
Статьи: 3
Проекты: 1

Re: Атмосферная коррекция данных Landsat/ETM+ (COST метод)

Сообщение Aleksander » 18 дек 2009, 16:35

Есть ошибка в формуле в разделе 2 пункт е) - ESUN должен быть в числителе, d2 - в знаменателе.
  • 0

Aleksander
Участник
 
Зарегистрирован: 12 май 2009
Откуда: Минск (на карте)
 
Сообщения: 60
Репутация: 0

Re: Атмосферная коррекция данных Landsat/ETM+ (COST метод)

Сообщение barash » 14 янв 2010, 20:39

Помимо ошибки, найденной 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, но я ее, увы, не нашла. Если у кого-то она есть, поделитесь, пожалуйста.
  • 0

barash
Новоприбывший
 
Зарегистрирован: 14 янв 2010
 
Сообщения: 1
Репутация: 0

Re: Атмосферная коррекция данных Landsat/ETM+ (COST метод)

Сообщение Cocodrilo » 11 сен 2012, 16:18

Посмотрел статью, появились вопросы:

1.
barash писал(а):Если взглянуть на статью 1989г, то на стр. 4 в формуле (5) вместо DN находится RAD, т.е. показатели излучения на сенсоре.

- Полностью согласен. Те DN надо еще пересчитать в L по формуле в пункте д.

2. Подскажите, пожалуйста, а как полученные значения соотносятся с радиометрическим разрешением канала?
Интересует вопрос калибровки данных разных сенсоров. Или вообще не нужно смотреть на это разрешение, если, говорится, что мы получили абсолютные безразмерные значения отражения?
  • 0

Последний раз редактировалось Cocodrilo 19 сен 2012, 15:53, всего редактировалось 1 раз.
Cocodrilo
Интересующийся
 
Зарегистрирован: 19 янв 2012
 
Сообщения: 38
Репутация: 3

След.

Вернуться в Материалы сайта

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 2


(Геокруг)

© GIS-Lab и авторы, 2002-2013. При использовании материалов сайта, ссылка на GIS-Lab и авторов обязательна. Содержание материалов - ответственность авторов (подробнее).