Форумы GIS-Lab.info

вопрос по формуле конвертации значий ячеек снимка в температуру для Landsat 7:

Для High Gain: (1282.71/((666.09/ ((((12.65-3.2)/255) *[Value]) +3.2) )+1 ).ln) - 273
Для Low Gain: (1282.71/((666.09/ ((((17.04)/255) *[Value])))+1 ).ln) - 273

предполагается, что цифра "255" - это "QCALMAX - QCALMIN" для сенсора ETM+, однако в .met файле на снимок от 19 августа 2000 года пишется:
GROUP = MIN_MAX_PIXEL_VALUE
...
QCALMAX_BAND61 = 255.0
QCALMIN_BAND61 = 1.0
QCALMAX_BAND62 = 255.0
QCALMIN_BAND62 = 1.0
...
т.е. цифра должна быть "254"?


и сразу же еще вопрос, где взять значения QCALMAX, QCALMIN для сенсора ТМ (Landsat 5)?

потому что в файле .hdr для снимка от 16 августа 1990 я не нашел эти данные:
PRODUCT =10110800049 WRS =180/028F ACQUISITION DATE =19900816 SATELLITE =L5 INSTRUMENT =TM10 PRODUCT TYPE =ORBIT ORIENTED PRODUCT SIZE =FULL SCENE TYPE OF GEODETIC PROCESSING =SYSTEMATIC RESAMPLING =CC RAD GAINS/BIASES = +0.6024/-1.5000 +1.1749/-2.8000 +0.8059/-1.2000 +0.8145/-1.5000 +0.1081/-0.3700 +0.0552/+1.2378 +0.0570/-0.1500 TAPE SPANNING FLAG=1/1 START LINE #= 1 LINES PER VOL= 5960 ORIENTATION =-13.58 PROJECTION =UTM 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 EARTH ELLIPSOID =WGS84 SEMI-MAJOR AXIS =6378137.000 SEMI-MINOR AXIS =6356752.314 PIXEL SIZE =30.00 PIXELS PER LINE= 6920 LINES PER IMAGE= 5960 UL 0292713.2962E 470014.0625N 230421.594 5211703.500 UR 0320651.5753E 463700.2957N 432186.469 5162957.000 LR 0313620.7243E 450251.4553N 390203.375 4989187.000 LL 0290102.3268E 452531.6040N 188438.469 5037933.500 BANDS PRESENT =1234567 BLOCKING FACTOR = 1 RECORD LENGTH = 6920 SUN ELEVATION =49 SUN AZIMUTH =132 CENTER 0303256.4747E 460152.0742N 310323.469 5100427.000 3460 2980 OFFSET= 0 REVB

Да, если в метаданных Qcalmin = 1, то надо использовать 254 а не 255. В целом, то что нужно использоватть зависит от процедуры предобработки и когда была сделана сцена.
Есть два технологических процесса, которыми может быть создана сцена: LPGS и NLAPS.

После 5 апреля 2004 для них обоих, Qcalmin = 1, Qcalmax = 255
До, для NLAPS 0-255, для LPGS 1-255.
Подробнее об этих и других отличиях здесь

По второму вопросу: http://gis-lab.info/qa/dn2radiance.html

Статью про температуру поправил с учетом этого.

к сожалению в http://gis-lab.info/qa/dn2radiance.html не указан QCALMAX для Landsat 5 TM, а в таблице с коэфициентами Lmin и Lmax нет 6 канала.

Из статьи: Qcalmax - максимальное калиброванное значение DN (255)

Таблицу обновил, значения теперь есть.