X
PrevPrev Go to previous topic
NextNext Go to next topic
Last Post 12 Mar 2009 12:22 PM by  anon
Landsat L1T (and perhaps other formats) Pre-processing Problem
 2 Replies
Sort:
You are not authorized to post a reply.
Author Messages

anon



New Member


Posts:
New Member


--
12 Mar 2009 12:22 PM
    I've think I've discovered a problem with the Landsat TM pre-processing functionality built into ENVI.  As accessed through the ENVI menu system, the tmcal_doit function does not correctly adjust for the gain/bias as used in the L1T processing for Landsat 5 or 7.  It appears as though it does not know that the minimum QCal DN value is 1, not 0, and miscalculates both the gain and bias when converting to at-sensor-radiance.  This causes errors in both calculations of at-sensor-radiance and top-of-atmosphere reflectance. For that reason, I've created an IDL script to parse the MTL metadata file and for each of the visual, nir- and mid-IR bands do a radiance conversion with gains/biases calculated from the metadata file and then stacks the resulting calibrated data.  It can process L1T data from both Landsat 5 and 7 (and I assume 4 but haven't tested it), and generates correct at-sensor-radiance ranges for non-zero QCal DNs. Is this problem on ITT's radar and will it be fixed in future versions or would there be interest in me cleaning up my processing script to release to the community?

    MariM



    Veteran Member


    Posts:2396
    Veteran Member


    --
    12 Mar 2009 01:33 PM
    We have a patch available for 4.6 now that will address this issue.  The new routine will open the geotiff files using the metadata (_MTL.txt) and group them into the appropriate files (VIS-SWIR, TIR, and Pan, if available).  It also uses the QCALMIN/QCALMAX and LMIN/LMAX from the metadata to calibrate.  The update allows the user to manually check and edit the LMIN/LMAX or gain/bias values.  For the case of Geotiff data that does not have metadata, you can read in the single band TIFF and based on the date you enter, ENVI will use hard-coded LMIN/LMAX values of the appropriate sensor for you (as published on the USGS website), or you can manually enter/change them.  With this update, the 'Get calibration parameters from the web' is no longer available.  Please contact Technical Support if you are interested in this patch.

    MariM



    Veteran Member


    Posts:2396
    Veteran Member


    --
    25 Feb 2010 10:01 AM
    The ENVI routine for it is the same TMCAL_DOIT. 
    You are not authorized to post a reply.