[Gfoss] landsat metadata

Local Chapters GFOSS.it
1

Ciao,

partendo dai dati scaricabili su :
http://glcf.umiacs.umd.edu/

avrei una domanda riguardo i metadata delle landsat :

mentre per le tm5 i valori di gain e bias
sono esplicitamente indicati nel file .hdr >

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

nel file di metadata delle etm+ invece ho :

GROUP = L1G_PRODUCT_METADATA
    BAND_COMBINATION = "123456678"
    CPF_FILE_NAME = "L7CPF20000719_20000930_09"
    GROUP = MIN_MAX_RADIANCE
      LMAX_BAND1 = 191.600
      LMIN_BAND1 = -6.200
      LMAX_BAND2 = 196.500
      LMIN_BAND2 = -6.400
      LMAX_BAND3 = 152.900
      LMIN_BAND3 = -5.000
      LMAX_BAND4 = 241.100
      LMIN_BAND4 = -5.100
      LMAX_BAND5 = 31.060
      LMIN_BAND5 = -1.000
      LMAX_BAND61 = 17.040
      LMIN_BAND61 = 0.000
      LMAX_BAND62 = 12.650
      LMIN_BAND62 = 3.200
      LMAX_BAND7 = 10.800
      LMIN_BAND7 = -0.350
      LMAX_BAND8 = 243.100
      LMIN_BAND8 = -4.700
    END_GROUP = MIN_MAX_RADIANCE
    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

leggendo il cap.11 del prodotto landsat :

http://landsathandbook.gsfc.nasa.gov/handbook/handbook_htmls/chapter11/chapter11.html

alla sezione : 11.3.1

trovo :

Lλ = "gain" * QCAL + "offset"

offset sarebbe il mio "bias" giusto ?

Lλ = ((LMAXλ - LMINλ)/(QCALMAX-QCALMIN)) * (QCAL-QCALMIN) + LMINλ

cio significa che per le ETM+ posso considerare :

gain : ((LMAXλ - LMINλ)/(QCALMAX-QCALMIN))

bias : LMINλ

sto sbaglando ?

ancora una domanda ... una delle scene che mi serve ... path-row : p189r32 ... tm5
ha il satellite molto inclinato ... cio mi dice che la regola del coseno non è adatta per una correzione topografica basata sul dtm ...

avete qualche link - pdf o doc vari che spieghino la regola di minnaert ?
sono interessato soprattutto a capire come e in base a cosa ... regolare i parametri presenti nella formula...

grazie!!!

buon lavoro.......
... e per chi può ...
.. happy dives!

2

http://www.isprs.org/istanbul2004/comm3/papers/347.pdf
http://www.asf.alaska.edu/educational/seminarpdf/2004/corrections.pdf
http://www.gmes-geoland.info/upload/files/CSL-0350-TN-02-BestPractice-TopographicNormalisation-I1.0.pdf

Dall’ultimo documento:

“The values of k vary between 0 and 1 for different land cover categories. Smaller values weaken the influence of the quotient in the Minnaert equation. Thus the overcorrecting in the cosine law for areas of low cos(i) is counteracted. The value of the Minnaert constant is determined empirically by linearising the Minnaert equation logarithmically and estimating the slope by linear regression. The value of the slope is the value of the Minnaert constant. As the Minnaert constant depends on land cover category, it is recommended to use forested areas with the same characteristics (forest type / crown cover / age class) for regression input (implemented in JR software IMPACT). Thus, an average regression function can be calculated, even though non-forested areas will be under- or over-corrected. Further research work still remains in order to improve the Minnaert correction and make it independent from pre-knowledge of land cover.”

La tua interpretazione di gain e bias per etm+ mi sembra corretta per QCALMIN = 0.
Per QCALMIN diversa da 0 hai invece:

bias = LMINλ - QCALMIN((LMAXλ - LMINλ)/(QCALMAX-QCALMIN))

…ma qui non capisco io come mai riscalino da 1 a 255…

Attenzione: ogni mia vaccata è attribuibile al lunedì mattina, nonché al fatto che è troppo tempo che non faccio preprocessing.

M

Il 06/07/07, massimo di stefano <massimodisasha@yahoo.it> ha scritto:

Ciao,

partendo dai dati scaricabili su :
http://glcf.umiacs.umd.edu/

avrei una domanda riguardo i metadata delle landsat :

mentre per le tm5 i valori di gain e bias
sono esplicitamente indicati nel file .hdr >

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

nel file di metadata delle etm+ invece ho :

GROUP = L1G_PRODUCT_METADATA
BAND_COMBINATION = “123456678”
CPF_FILE_NAME = “L7CPF20000719_20000930_09”
GROUP = MIN_MAX_RADIANCE
LMAX_BAND1 = 191.600
LMIN_BAND1 = -6.200
LMAX_BAND2 = 196.500
LMIN_BAND2 = -6.400
LMAX_BAND3 = 152.900
LMIN_BAND3 = -5.000
LMAX_BAND4 = 241.100
LMIN_BAND4 = - 5.100
LMAX_BAND5 = 31.060
LMIN_BAND5 = -1.000
LMAX_BAND61 = 17.040
LMIN_BAND61 = 0.000
LMAX_BAND62 = 12.650
LMIN_BAND62 = 3.200
LMAX_BAND7 = 10.800
LMIN_BAND7 = -0.350
LMAX_BAND8 = 243.100
LMIN_BAND8 = - 4.700
END_GROUP = MIN_MAX_RADIANCE
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

leggendo il cap.11 del prodotto landsat :

http://landsathandbook.gsfc.nasa.gov/handbook/handbook_htmls/
chapter11/chapter11.html

alla sezione : 11.3.1

trovo :

Lλ = “gain” * QCAL + “offset”

offset sarebbe il mio “bias” giusto ?

Lλ = ((LMAXλ - LMINλ)/(QCALMAX-QCALMIN)) * (QCAL-QCALMIN) + LMINλ

cio significa che per le ETM+ posso considerare :

gain : ((LMAXλ - LMINλ)/(QCALMAX-QCALMIN))

bias : LMINλ

sto sbaglando ?

ancora una domanda … una delle scene che mi serve … path-row :
p189r32 … tm5
ha il satellite molto inclinato … cio mi dice che la regola del
coseno non è adatta per una correzione topografica basata sul dtm …

avete qualche link - pdf o doc vari che spieghino la regola di
minnaert ?
sono interessato soprattutto a capire come e in base a cosa …
regolare i parametri presenti nella formula…

grazie!!!

buon lavoro…
… e per chi può …
… happy dives!

Gfoss mailing list: 220 iscritti (28-05-2007)
Gfoss@faunalia.com
http://www.faunalia.com/cgi-bin/mailman/listinfo/gfoss