Конвертация данных TM,ETM+ в показатели излучения на сенсоре

Обсуждение материалов сайта: вопросы, замечания, предложения
Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9128
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 747
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:

Конвертация данных TM,ETM+ в показатели излучения на сенсоре

Сообщение Максим Дубинин » 30 окт 2004, 00:44

Обсуждение темы:
Конвертация данных TM, ETM+ в показатели излучения на сенсоре

http://gis-lab.info/qa/dn2radiance.html
Последний раз редактировалось Максим Дубинин 28 мар 2007, 20:02, всего редактировалось 1 раз.
пристегивайтесь, турбулентность прямо по курсу

Andrey

Сообщение Andrey » 30 окт 2004, 07:01

Кого интересует расчет температуры по данным теплового диапазона Landsat TM/ETM (band 6) рекомендую ознакомиться с докой
http://earthsys.ag.ohio-state.edu/nr797/Exercise3.doc
Также алгоритмы расчета можно скачать в виде моделей к ERDAS со страницы доступных моделей http://www.gis.leica-geosystems.com/Sup ... fault.aspx
или по прямым ссылкам
Convert Landsat 7 spectral radiance to temperature.
http://www.gis.leica-geosystems.com/Sup ... nload_ID=3
Converting spectral radiance from band 6 of Landsat 4 and Landsat 5 data to temperature values in degrees Celsius.
http://www.gis.leica-geosystems.com/Sup ... nload_ID=7

Могу предложить также алгоритм расчета температуры по данным band 6 Landsat 5, заимствованный из вышеупомянутой модели ERDAS, но в виде ENVI EXPRESSION (содержание файла *.exp):
ENVI EXPRESSIONS
1.238+(15.600-1.238)*b6/255
1260.56/(607.76/b1+1)

Процедура состоит из двух этапов:
1. Convert Calibrated Digital numbers to Spectral Radiance
1.238+(15.600-1.238)*b6/255
(b6 = 6band landsat5)
2. Convert Spectral Radiance to Temperature
1260.56/(607.76/b1+1)
(b1 = Spectral Radiance)

Tanya
Новоприбывший
Сообщения: 1
Зарегистрирован: 18 май 2006, 13:44
Репутация: 0

Сообщение Tanya » 18 май 2006, 13:49

Нигде не могу найти калибровочные файлы для Landsat 4. Подскажите, пожалуйста!

AlexChe
Новоприбывший
Сообщения: 5
Зарегистрирован: 30 май 2006, 16:41
Репутация: 0
Откуда: Moscow

Сообщение AlexChe » 31 май 2006, 14:13

Очень нужно пересчитать полученные значения для Landsat7/ETM+ в абсолютные безразмерные значения отражения (reflectance).
в целом понятно что пересчет проводить по формуле
Изображение
из которой не совсем понятно как выбирать Изображение? данные есть в таблице, а что делать с промежуточными датами неясно.
с остальным вроде как все понятно.
Если кто-нибудьподелиться своим опытом по данному вопросу будуочень признателен.
Заранее спасибо!

GeoFaker
Новоприбывший
Сообщения: 5
Зарегистрирован: 02 дек 2006, 20:16
Репутация: 0
Откуда: Москва

Сообщение GeoFaker » 02 дек 2006, 20:27

В этой формуле пересчета d - расстояние между Солнцем и Землей в астрономических единицах. Эти данные есть на официальном сайте Landsat (http://landsat.usgs.gov/) и скорее всего в метафайлах (файл,поставляющийся вместе со снимком и имеющий расширение *.met, его можно открыть текстовым редактором). У меня есть значения d для следующих сроков съемки:
Landsat 7 22.03.2001 d=0.9965
L7 27.04.2000 d=1.0067
L7 20.06.2002 d=1.0161
L7 27.09.2000 d=1.0002
L5 03.05.1990 d=1.0081

Это значение d кажется не зависит от места съемки, но снимки все на ЕТР (Тверская, Калужская области)

alexandr cherepanov
Гуру
Сообщения: 534
Зарегистрирован: 30 ноя 2006, 13:31
Репутация: 116
Откуда: Moscow

Сообщение alexandr cherepanov » 02 дек 2006, 20:36

со значением d уже давно разобрался :D
Эта величина не постоянная на сколько я понял и зависит от времени съемки и координат места, расчитать ее можно на сайте http://ssd.jpl.nasa.gov/horizons.cgi в статьях есть подробное описание как это сделать http://gis-lab.info/qa/earthsundist.html

GeoFaker
Новоприбывший
Сообщения: 5
Зарегистрирован: 02 дек 2006, 20:16
Репутация: 0
Откуда: Москва

Сообщение GeoFaker » 03 дек 2006, 15:59

Меня сейчас волнует такая проблема: там в формуле расчета отраженного солнечного коротковолнового излучения и собственного длинноволнового излучения земной поверхности (ну где L=((Lmax-Lmin)/(Qcalmax - Qcalmin))*(Qcal-Qcalmin)+Lmin значение Qcalmax может быть как=0 или как=1. На форуме написано очень расплывчато о том как определить чему равно это значение. меня интересуют как раз значения Qcalmax для 27.04.2000, 27.09.2000, 03.05.1990 (последняя дата - Landsat5) [/img]

alexandr cherepanov
Гуру
Сообщения: 534
Зарегистрирован: 30 ноя 2006, 13:31
Репутация: 116
Откуда: Moscow

Сообщение alexandr cherepanov » 04 дек 2006, 08:37

А мета файле разве нету описания чему они равны? покрайне мере для первых 2 снимков должны быть
вот кусочек из мета к p160r021_7t20020728_z41

GROUP = MIN_MAX_PIXEL_VALUE
QCALMAX_BAND1 = 255.0
QCALMIN_BAND1 = 1.0
QCALMAX_BAND2 = 255.0
QCALMIN_BAND2 = 1.0
QCALMAX_BAND3 = 255.0
QCALMIN_BAND3 = 1.0
QCALMAX_BAND4 = 255.0
QCALMIN_BAND4 = 1.0
QCALMAX_BAND5 = 255.0
QCALMIN_BAND5 = 1.0
QCALMAX_BAND61 = 255.0
QCALMIN_BAND61 = 1.0
QCALMAX_BAND62 = 255.0
QCALMIN_BAND62 = 1.0
QCALMAX_BAND7 = 255.0
QCALMIN_BAND7 = 1.0
QCALMAX_BAND8 = 255.0
QCALMIN_BAND8 = 1.0
END_GROUP = MIN_MAX_PIXEL_VALUE

GeoFaker
Новоприбывший
Сообщения: 5
Зарегистрирован: 02 дек 2006, 20:16
Репутация: 0
Откуда: Москва

Сообщение GeoFaker » 04 дек 2006, 12:14

Точно, конечно есть,я забыл совсем про метафайл. Только вот на 22.03.2001 у меня почему-то нет метафайла,пытался найти параметры калибровки на http://landsat.usgs.gov/, но запутался в обозначениях.....

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

Сообщение Robert » 06 дек 2006, 00:32

Обычно для каждого снимка метафайл есть в комплекте: файл с окончанием MTL: например L71181021_02120010322_MTL. Так что надо спрашивать там, где брал снимок ;).
У меня как раз заволялся... прислать?)))
...чувствовал себя дикарём среди миклухо-маклаев...

GeoFaker
Новоприбывший
Сообщения: 5
Зарегистрирован: 02 дек 2006, 20:16
Репутация: 0
Откуда: Москва

Сообщение GeoFaker » 06 дек 2006, 00:39

Давай Роб, пришли, ага. А то по другому видимо никогда мне не пришлешь......Штеф.

Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9128
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 747
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:

Сообщение Максим Дубинин » 28 мар 2007, 20:01

обновлена таблица с коэфициентами для L5-TM, добавлен 6 канал и источник таблицы
пристегивайтесь, турбулентность прямо по курсу

olgatut
Новоприбывший
Сообщения: 4
Зарегистрирован: 30 май 2009, 21:16
Репутация: 5

Re: Конвертация данных TM,ETM+ в показатели излучения на сенсоре

Сообщение olgatut » 30 май 2009, 21:24

Следует отметить, что есть новые уточненные параметры для пересчета значений Landsat 5 TM, отличия могут достигать 15% для снимков 1984-1991 гг.. См.:

Gyanesh Chander, Member, IEEE, Brian L. Markham, and Julia A. Barsi. 2007. Revised Landsat-5 Thematic Mapper Radiometric Calibration IEEE GEOSCIENCE AND REMOTE SENSING LETTERS, VOL. 4, NO. 3, JULY 2007
http://neptune.gsfc.nasa.gov/publicatio ... Mapper.pdf

alexandr cherepanov
Гуру
Сообщения: 534
Зарегистрирован: 30 ноя 2006, 13:31
Репутация: 116
Откуда: Moscow

Re: Конвертация данных TM,ETM+ в показатели излучения на сенсоре

Сообщение alexandr cherepanov » 11 ноя 2009, 10:39

изменилась ссылка на handbook
http://landsathandbook.gsfc.nasa.gov/ha ... k_toc.html

timorlov
Новоприбывший
Сообщения: 2
Зарегистрирован: 04 апр 2010, 20:20
Репутация: 0

Re: Конвертация данных TM,ETM+ в показатели излучения на сенсоре

Сообщение timorlov » 04 апр 2010, 20:46

Добрый день,

А не подскажете, где можно взять данные для конвертации для снимков QuickBird и Ikonos?

Ответить

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

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

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