Форумы GIS-Lab.info

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

http://gis-lab.info/qa/dn2radiance.html

Кого интересует расчет температуры по данным теплового диапазона 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)

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

Очень нужно пересчитать полученные значения для Landsat7/ETM+ в абсолютные безразмерные значения отражения (reflectance).
в целом понятно что пересчет проводить по формуле

из которой не совсем понятно как выбирать ? данные есть в таблице, а что делать с промежуточными датами неясно.
с остальным вроде как все понятно.
Если кто-нибудьподелиться своим опытом по данному вопросу будуочень признателен.
Заранее спасибо!

В этой формуле пересчета 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 кажется не зависит от места съемки, но снимки все на ЕТР (Тверская, Калужская области)

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

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

А мета файле разве нету описания чему они равны? покрайне мере для первых 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

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

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

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

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

Следует отметить, что есть новые уточненные параметры для пересчета значений 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

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

Добрый день,

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