Конвертация данных TM,ETM+ в показатели излучения на сенсоре
- Максим Дубинин
- MindingMyOwnBusiness
- Сообщения: 9129
- Зарегистрирован: 06 окт 2003, 20:20
- Репутация: 747
- Ваше звание: NextGIS
- Откуда: Москва
- Контактная информация:
Конвертация данных TM,ETM+ в показатели излучения на сенсоре
Обсуждение темы:
Конвертация данных TM, ETM+ в показатели излучения на сенсоре
http://gis-lab.info/qa/dn2radiance.html
Конвертация данных TM, ETM+ в показатели излучения на сенсоре
http://gis-lab.info/qa/dn2radiance.html
Последний раз редактировалось Максим Дубинин 28 мар 2007, 20:02, всего редактировалось 1 раз.
пристегивайтесь, турбулентность прямо по курсу
Кого интересует расчет температуры по данным теплового диапазона 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)
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)
-
- Новоприбывший
- Сообщения: 1
- Зарегистрирован: 18 май 2006, 13:44
- Репутация: 0
-
- Новоприбывший
- Сообщения: 5
- Зарегистрирован: 30 май 2006, 16:41
- Репутация: 0
- Откуда: Moscow
Очень нужно пересчитать полученные значения для Landsat7/ETM+ в абсолютные безразмерные значения отражения (reflectance).
в целом понятно что пересчет проводить по формуле
![Изображение](http://ltpwww.gsfc.nasa.gov/IAS/handbook/handbook_htmls/chapter11/images/reflect_formula.gif)
из которой не совсем понятно как выбирать
? данные есть в таблице, а что делать с промежуточными датами неясно.
с остальным вроде как все понятно.
Если кто-нибудьподелиться своим опытом по данному вопросу будуочень признателен.
Заранее спасибо!
в целом понятно что пересчет проводить по формуле
![Изображение](http://ltpwww.gsfc.nasa.gov/IAS/handbook/handbook_htmls/chapter11/images/reflect_formula.gif)
из которой не совсем понятно как выбирать
![Изображение](http://ltpwww.gsfc.nasa.gov/IAS/handbook/handbook_htmls/chapter11/images/earth_sun.gif)
с остальным вроде как все понятно.
Если кто-нибудьподелиться своим опытом по данному вопросу будуочень признателен.
Заранее спасибо!
-
- Новоприбывший
- Сообщения: 5
- Зарегистрирован: 02 дек 2006, 20:16
- Репутация: 0
- Откуда: Москва
В этой формуле пересчета 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 кажется не зависит от места съемки, но снимки все на ЕТР (Тверская, Калужская области)
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 кажется не зависит от места съемки, но снимки все на ЕТР (Тверская, Калужская области)
-
- Гуру
- Сообщения: 534
- Зарегистрирован: 30 ноя 2006, 13:31
- Репутация: 116
- Откуда: Moscow
со значением d уже давно разобрался
Эта величина не постоянная на сколько я понял и зависит от времени съемки и координат места, расчитать ее можно на сайте http://ssd.jpl.nasa.gov/horizons.cgi в статьях есть подробное описание как это сделать http://gis-lab.info/qa/earthsundist.html
![Very Happy :D](./images/smilies/icon_biggrin.gif)
Эта величина не постоянная на сколько я понял и зависит от времени съемки и координат места, расчитать ее можно на сайте http://ssd.jpl.nasa.gov/horizons.cgi в статьях есть подробное описание как это сделать http://gis-lab.info/qa/earthsundist.html
-
- Новоприбывший
- Сообщения: 5
- Зарегистрирован: 02 дек 2006, 20:16
- Репутация: 0
- Откуда: Москва
Меня сейчас волнует такая проблема: там в формуле расчета отраженного солнечного коротковолнового излучения и собственного длинноволнового излучения земной поверхности (ну где L=((Lmax-Lmin)/(Qcalmax - Qcalmin))*(Qcal-Qcalmin)+Lmin значение Qcalmax может быть как=0 или как=1. На форуме написано очень расплывчато о том как определить чему равно это значение. меня интересуют как раз значения Qcalmax для 27.04.2000, 27.09.2000, 03.05.1990 (последняя дата - Landsat5) [/img]
-
- Гуру
- Сообщения: 534
- Зарегистрирован: 30 ноя 2006, 13:31
- Репутация: 116
- Откуда: Moscow
А мета файле разве нету описания чему они равны? покрайне мере для первых 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
вот кусочек из мета к 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
-
- Новоприбывший
- Сообщения: 5
- Зарегистрирован: 02 дек 2006, 20:16
- Репутация: 0
- Откуда: Москва
Точно, конечно есть,я забыл совсем про метафайл. Только вот на 22.03.2001 у меня почему-то нет метафайла,пытался найти параметры калибровки на http://landsat.usgs.gov/, но запутался в обозначениях.....
-
- Новоприбывший
- Сообщения: 5
- Зарегистрирован: 06 дек 2006, 00:24
- Репутация: 0
- Откуда: геофак
- Контактная информация:
-
- Новоприбывший
- Сообщения: 5
- Зарегистрирован: 02 дек 2006, 20:16
- Репутация: 0
- Откуда: Москва
- Максим Дубинин
- MindingMyOwnBusiness
- Сообщения: 9129
- Зарегистрирован: 06 окт 2003, 20:20
- Репутация: 747
- Ваше звание: NextGIS
- Откуда: Москва
- Контактная информация:
-
- Новоприбывший
- Сообщения: 4
- Зарегистрирован: 30 май 2009, 21:16
- Репутация: 5
Re: Конвертация данных TM,ETM+ в показатели излучения на сенсоре
Следует отметить, что есть новые уточненные параметры для пересчета значений 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
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
-
- Гуру
- Сообщения: 534
- Зарегистрирован: 30 ноя 2006, 13:31
- Репутация: 116
- Откуда: Moscow
Re: Конвертация данных TM,ETM+ в показатели излучения на сенсоре
изменилась ссылка на handbook
http://landsathandbook.gsfc.nasa.gov/ha ... k_toc.html
http://landsathandbook.gsfc.nasa.gov/ha ... k_toc.html
-
- Новоприбывший
- Сообщения: 2
- Зарегистрирован: 04 апр 2010, 20:20
- Репутация: 0
Re: Конвертация данных TM,ETM+ в показатели излучения на сенсоре
Добрый день,
А не подскажете, где можно взять данные для конвертации для снимков QuickBird и Ikonos?
А не подскажете, где можно взять данные для конвертации для снимков QuickBird и Ikonos?
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 0 гостей