Nikos wrote: > > While i.landsat.toar works just fine with all bands in converting DNs to > > Radiance/Reflectance (with some method=dos?), it seems, at least in once > > case (scene LE71840312002041SGS00), it fails to convert bands 61 and 62 in > > temperature!
> > Same happens (as described below in both G64 and G70). > > Can anyone confirm or reject? > [deleted stuff] > > Using the old metadata file (in this case: LE71840312002041SGS00_MTLold.txt) > works fine. > > GRASS 6.4.3svn (wgs84_utm_z34_test):~ > r.info -r B.DOS1.ToCR.61 > min=158.263805493393 > max=338.623091571347 > > GRASS 6.4.3svn (wgs84_utm_z34_test):~ > r.info -r B.DOS1.ToCR.62 > min=240.06999848884 > max=322.080084447704 > I guess it is a bug (something that got away in latest updates of > i.landsat.toar towards the newer metadata file). Will try to check > further... not today though, deadlines approaching. This post passed unnoticed. Currently, newer MTL files contain for temperature data file names the following RADIANCE_MAXIMUM_BAND_6_VCID_1 = 17.040 RADIANCE_MINIMUM_BAND_6_VCID_1 = 0.000 RADIANCE_MAXIMUM_BAND_6_VCID_2 = 12.650 RADIANCE_MINIMUM_BAND_6_VCID_2 = 3.200 .. QUANTIZE_CAL_MAX_BAND_6_VCID_1 = 255 QUANTIZE_CAL_MIN_BAND_6_VCID_1 = 1 QUANTIZE_CAL_MAX_BAND_6_VCID_2 = 255 QUANTIZE_CAL_MIN_BAND_6_VCID_2 = 1 as opposed to older MTL files which instruct LMAX_BAND61 = 17.040 LMIN_BAND61 = 0.000 LMAX_BAND62 = 12.650 LMIN_BAND62 = 3.200 .. QCALMAX_BAND61 = 255.0 QCALMIN_BAND61 = 1.0 QCALMAX_BAND62 = 255.0 QCALMIN_BAND62 = 1.0 Can you give me a hint (I want to --at least-- try to understand and maybe fix this!) on how to adapt the landsat_met.c file so as to account for RADIANCE_MAXIMUM_BAND_6_VCID_%d RADIANCE_MINIMUM_BAND_6_VCID_%d QUANTIZE_CAL_MAX_BAND_6_VCID_%d QUANTIZE_CAL_MIN_BAND_6_VCID_%d which are in newer MTL files? Or is it to complex for C-ignorants like me? Thanks, N _______________________________________________ grass-dev mailing list grass-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-dev