X
PrevPrev Go to previous topic
NextNext Go to next topic
Last Post 22 Jan 2018 06:14 AM by  MariM
ASTER - FLAASH
 5 Replies
Sort:
You are not authorized to post a reply.
Author Messages

Ludovic Legros



New Member


Posts:29
New Member


--
19 Jan 2018 10:27 AM
    Hello,

    I tried an atmospherical correction with FLAASH on ASTER L1T data (The ASTER Level 1 Precision Terrain Corrected Registered At-Sensor Radiance (AST_L1T) data contains calibrated at-sensor radiance, which corresponds with the ASTER Level 1B (AST_L1B), that has been geometrically corrected, and rotated to a north-up UTM projection), but my values are really odd.

    I have strong positive and negative values (+ 32 010 to -32 752). I tried to follow step by step the tutorial here:
    http://www.harrisgeospati...nvi/ASTER_FLAASH.pdf

    But it doesn't work well. The tutorial is about ASTER L1A data, do you think the problems comes from the type of data ?

    Here's the template file:

    ;
    ;ENVI FLAASH PARAMETERS TEMPLATE (5.4.1)
    ;Written Fri Jan 19 12:16:18 2018
    ;
    ; Project Parameters
    enviacc.prj.radiance_file = G:\LL_WORKSPACE\OTR-SI411_GoldSpot_Arizona\02_Rasters\01_Images_base\ASTER\HDF\FLAASH\ASTER_GoldSpot_Arizona_VNIR-SWIR_FLAASH
    enviacc.prj.reflect_file = G:\LL_WORKSPACE\OTR-SI411_GoldSpot_Arizona\02_Rasters\01_Images_base\ASTER\HDF\FLAASH\ASTER_GoldSpot_Arizona_VNIR-SWIR_FLAASH
    enviacc.prj.filter_func_file = G:\ENVI_5.4.1_Trial\ENVI54\classic\filt_func\aster.sli
    enviacc.prj.filter_func_file_index = 0
    enviacc.prj.water_band_choice = 1.13
    enviacc.prj.red_channel = 2
    enviacc.prj.green_channel = 1
    enviacc.prj.blue_channel = 0
    enviacc.prj.water_abs_channel = 0
    enviacc.prj.water_ref_channel = 0
    enviacc.prj.kt_upper_channel = 0
    enviacc.prj.kt_lower_channel = 2
    enviacc.prj.kt_cutoff = 0.0800
    enviacc.prj.kt_ratio = 0.5000
    enviacc.prj.cirrus_channel = 0
    enviacc.prj.water_retrieval = 0
    enviacc.prj.modtran_directory = G:\LL_WORKSPACE\OTR-SI411_GoldSpot_Arizona\02_Rasters\01_Images_base\ASTER\HDF\FLAASH\Fichier_FLAASH\
    ;
    ; MODTRAN Parameters
    enviacc.modtran.visvalue = 40.0000
    enviacc.modtran.f_resolution = 15.0000
    enviacc.modtran.day = 12
    enviacc.modtran.month = 9
    enviacc.modtran.year = 2003
    enviacc.modtran.gmt = 18.2408
    enviacc.modtran.latitude = 33.6121
    enviacc.modtran.longitude = -110.7300
    enviacc.modtran.sensor_altitude = 705.0000
    enviacc.modtran.ground_elevation = 1.7000
    enviacc.modtran.view_zenith_angle = 180.0000
    enviacc.modtran.view_azimuth = 0.0000
    enviacc.modtran.atmosphere_model = 2
    enviacc.modtran.aerosol_model = 1
    enviacc.modtran.multiscatter_model = 2
    enviacc.modtran.disort_streams = 8
    enviacc.modtran.co2mix = 390.0000
    enviacc.modtran.water_column_multiplier = 1.0000
    ;
    ; Image Parameters
    enviacc.img.nspatial = 5630
    enviacc.img.nlines = 4994
    enviacc.img.data_type = 2
    enviacc.img.margin1 = 0
    enviacc.img.margin2 = 0
    enviacc.img.nskip = 0
    enviacc.img.pixel_size = 15.0000
    enviacc.img.sensor_name = ASTER
    ;
    ; Analysis Parameters
    enviacc.ana.aerosol_scaleht = 1.5000
    enviacc.ana.use_adjacency = 1
    enviacc.ana.output_scale = 10000.0000
    enviacc.ana.output_diagnostic_files = 0
    enviacc.ana.auto_save_template = 1
    enviacc.ana.polishing_res = 0
    enviacc.ana.aerosol_retrieval = 0
    enviacc.ana.calc_wl_correction = 0
    enviacc.ana.reuse_modtran_calcs = 0
    enviacc.ana.use_square_slit_function = 0
    enviacc.ana.convolution_method = fft
    enviacc.ana.use_tiling = 1
    enviacc.ana.tile_size = 100.0000
    ;
    ; GParameters Parameters
    enviacc.gparameters.gfact = 0
    enviacc.gparameters.gval = 0.0000
    enviacc.gparameters.use_lr_method = 0
    enviacc.gparameters.multikernel = 1
    enviacc.gparameters.autoroi = auto
    enviacc.gparameters.selarg = a
    ;
    ; Spectral Parameters
    enviacc.spc.wavelength_units = micron
    enviacc.spc.lambda = [
    0.5560, 0.6610, 0.8070, 1.6560, 2.1670, 2.2090, 2.2620,
    2.3360, 2.4000]
    enviacc.spc.fwhm = [
    0.080000, 0.060000, 0.100000, 0.100000, 0.040000, 0.040000, 0.050000,
    0.070000, 0.070000]
    enviacc.img.p_input_scale = [
    10.0000, 10.0000, 10.0000, 10.0000, 10.0000, 10.0000, 10.0000,
    10.0000, 10.0000]


    MariM



    Veteran Member


    Posts:2396
    Veteran Member


    --
    19 Jan 2018 11:50 AM
    When the results from FLAASH hit the maximum and minimum like this, it usually means that the scale factor used on import was incorrect. Did you first use Radiometric Calibration with 'FLAASH Settings" prior to FLAASH? If so, then the scale factor you would use would be 1.0. I would do the following:
    1. Process the VNIR and SWIR bands using Radiometric Calibration with FLAASH Settings
    2. Layer Stack the VNIR+SWIR
    3. Convert interleave back to BIL
    4. Process cube in FLAASH

    Ludovic Legros



    New Member


    Posts:29
    New Member


    --
    19 Jan 2018 01:27 PM
    As the ASTER Level 1T seemed to be already corrected, I thought it wasn't necessary to run a radiometric calibration...

    So I did a Radiometric calibration and then I ran FLAASH with a scale factor of 1.0 and it works.

    I have a question, the reflectance values shouldn't be from 0.0 to 1.0 (0 to 100% of reflectance) ?
    Because I have values from 300 to 1700.

    Thanks.

    MariM



    Veteran Member


    Posts:2396
    Veteran Member


    --
    19 Jan 2018 02:12 PM
    FLAASH outputs scaled reflectance stored as integers where the scaling factor is 10000. You can use Band Math to convert the scaled reflectance to percent using:

    float(b1) / 10000
    and use the Map to File option to select the entire reflectance cube. However, those values look fairly low if 1700 is the maximum, unless it is an image of mostly water.

    Ludovic Legros



    New Member


    Posts:29
    New Member


    --
    19 Jan 2018 02:37 PM
    I made a mistake, 1700 was a visual estimation of the maximum for band 1.

    Here are the stats:

    Basic Stats Min Max Mean Stdev
    Band 1 -720 4324 542.532637 502.969436
    Band 2 -21 5235 703.287215 660.983475
    Band 3 -325 4253 1274.252536 1050.116195
    Band 4 -68 4922 1422.767457 1196.905315
    Band 5 -28 4196 1050.507477 897.434710
    Band 6 -27 3989 1033.907243 887.104901
    Band 7 -26 4387 1045.919617 893.798647
    Band 8 -23 3838 895.856264 772.731708
    Band 9 -22 4319 917.721718 778.632051

    It seems ok ?

    MariM



    Veteran Member


    Posts:2396
    Veteran Member


    --
    22 Jan 2018 06:14 AM
    Yes, those look reasonable at 30-50% reflectance. It is not uncommon to have some negative values in pixels that are over water or shadow. They do not model well and are too near %0 reflectance/radiance.
    You are not authorized to post a reply.