X
PrevPrev Go to previous topic
NextNext Go to next topic
Last Post 26 Jun 2016 07:40 PM by  anon
FLAASH error: acc_called_t5set
 7 Replies
Sort:
You are not authorized to post a reply.
Author Messages

anon



New Member


Posts:
New Member


--
26 Jun 2016 07:40 PM
    Hi, I'm processing GOCI (Geostationary Ocean Color Imager launched by Korea) radiance data, and tried to use FLAASH to correct atmospheric effects.I know there are some requirements of the data format by FLAASH, so I converted the radiance data type to float, convert BSQ to BIL, correct the image coordinates. FLAASH failed and the following error poped up: ACC error: acc_called_t5set:read tape5 back in to check format. IDL error: Unable to apply format code F to input: "**********". 2016/6/27 ACC_CALLED_T5TEST ACC_HYPER2 ACC_ACC ENVI_ACC_EVENT ENVI_FLAASH_PROCESS_MANAGER $MAIN$ Here is the template.txt: ; ;ENVI FLAASH PARAMETERS TEMPLATE (5.1) ;Written Sun Jun 26 18:57:43 2016 ; ; Project Parameters enviacc.prj.radiance_file = D:\GOCI\taihu_rad enviacc.prj.reflect_file = D:\GOCI\taihu_flaash enviacc.prj.filter_func_file = D:\GOCI\GOCI_SRF enviacc.prj.filter_func_file_index = 0 enviacc.prj.water_band_choice = 1.13 enviacc.prj.red_channel = 5 enviacc.prj.green_channel = 4 enviacc.prj.blue_channel = 2 enviacc.prj.water_abs_channel = 0 enviacc.prj.water_ref_channel = 0 enviacc.prj.kt_upper_channel = 0 enviacc.prj.kt_lower_channel = 5 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 = D:\GOCI\ ; ; MODTRAN Parameters enviacc.modtran.visvalue = 40.0000 enviacc.modtran.f_resolution = 5.0000 enviacc.modtran.day = 17 enviacc.modtran.month = 8 enviacc.modtran.year = 2015 enviacc.modtran.gmt = 0.0000 enviacc.modtran.latitude = 31.2434 enviacc.modtran.longitude = 120.2392 enviacc.modtran.sensor_altitude = 35786.0000 enviacc.modtran.ground_elevation = 0.0500 enviacc.modtran.view_zenith_angle = 180.0000 enviacc.modtran.view_azimuth = 0.0000 enviacc.modtran.atmosphere_model = 2 enviacc.modtran.aerosol_model = 5 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 = 133 enviacc.img.nlines = 138 enviacc.img.data_type = 4 enviacc.img.margin1 = 0 enviacc.img.margin2 = 0 enviacc.img.nskip = 0 enviacc.img.pixel_size = 500.0000 enviacc.img.sensor_name = UNKNOWN-MSI ; ; Analysis Parameters enviacc.ana.aerosol_scaleht = 1.5000 enviacc.ana.use_adjacency = 1 enviacc.ana.output_scale = 10000.0000 enviacc.ana.polishing_res = 9 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 = nanometer enviacc.spc.lambda = [ 412.0000, 443.0000, 490.0000, 555.0000, 660.0000, 680.0000, 745.0000, 865.0000] enviacc.spc.fwhm = [ 20.000000, 20.000000, 20.000000, 20.000000, 20.000000, 10.000000, 20.000000, 40.000000] enviacc.img.p_input_scale = [ 10.0000, 10.0000, 10.0000, 10.0000, 10.0000, 10.0000, 10.0000, 10.0000] Please help. Thanks!

    MariM



    Veteran Member


    Posts:2396
    Veteran Member


    --
    27 Jun 2016 07:10 AM
    It could be the altitude: enviacc.modtran.sensor_altitude = 35786.0000 In FLAASH, this is defined in km so this would be 35786km. Did you enter this in meters?

    Deleted User



    New Member


    Posts:
    New Member


    --
    27 Jun 2016 08:26 PM
    Thank you for your help! I’m sure the sensor altitude is 35786km because the GOCI sensor is equipped in the geostationary orbit satellite. Here is the wiki about GOCI https://en.wikipedia.org/wiki/Geostationary_Ocean_Color_Imager Is there any limitation about the sensor altitude in FLAASH?

    MariM



    Veteran Member


    Posts:2396
    Veteran Member


    --
    28 Jun 2016 11:24 AM
    No, I do not believe it is a limit in altitude and it does look correct. This error sounds like it is unable to read something from a file. How does your spectral response (filter) function look? Can you also post the error from the journal file to see the last few processes that took place before the error?

    Deleted User



    New Member


    Posts:
    New Member


    --
    29 Jun 2016 03:30 AM
    I am sure the spectral response function is the same as the image in this link :http://222.236.46.36/page...201301/20130104.html And I convert the ASCII file to the .sli file which can be opened by the spectral library viewer in ENVI. Here is the journal.txt: Warning: first requested MODTRAN run calls for existing statistics; changing to no-statistics *****************FLAASH-SLANT**************** Spectral Sciences Atmospheric Correction Code 1998-2008 Version 01302013 ** See Instructions for Use ** ************* MAIN MENU *************** Commands: m = set up MODTRAN input file r = Run MODTRAN h = Hyperspectral data analysis a = About FLAASH (Instructions) c = Convert existing MODTRAN tape7 w = Read Warfighter run file, create initial tape5 q = Quit script command h Path for output files (except output reflectance cube) is D:\GOCI\ Enter u = Use this path or a = select Another path : script command u ************ DATA ANALYSIS MENU ************* Setup Commands: S = I/O Setup, interactive s = I/O Setup via serial queries d = new radiance Data cube a = Adjacency algorithm setup (create averaged image) q = Quit to main menu Processing Commands: i = simulate radiance Image from reflectance M = prepare MODTRAN run R = run MODTRAN WC = Channel Wavelength Correction Utility Commands: ? = menu help p = Plot pixel spectra e = display Existing reflectance file script command d Reading header (if any) from image file D:\GOCI\taihu_rad Enter u = Use this file, a = select Another file (dialog), or d = select another file (Direct filename entry) : script command u Select method for representing instrument spectral filter function. Type 1 for calculated function, 2 for a function defined by a file, 3 for external channel definitions passed in. ...type 99 to accept default from call. script command 99 ************ DATA ANALYSIS MENU ************* Setup Commands: S = I/O Setup, interactive s = I/O Setup via serial queries d = new radiance Data cube l = Load and display radiance data cube a = Adjacency algorithm setup (create averaged image) q = Quit to main menu Processing Commands: i = simulate radiance Image from reflectance v = aerosol Visibility determination k = Kaufman (ratio) visibility retrieval M = prepare MODTRAN run R = run MODTRAN WC = Channel Wavelength Correction Utility Commands: ? = menu help p = Plot pixel spectra e = display Existing reflectance file script command l Channel definitions: type 1 to use passed-in default, 2 to select manually between automatic and external-file channel definitions. script command 1 Warning: KT calculation not defined in channel selection file D:\GOCI\channel_def.txt Warning: vegetation parameters not found in channel selection file D:\GOCI\channel_def.txt Loading spectral data... 9045 blank pixels out of 18354 were found. *Not displaying RGB radiance image, cirrus cloud image and last-line spectra ************ DATA ANALYSIS MENU ************* Setup Commands: S = I/O Setup, interactive s = I/O Setup via serial queries d = new radiance Data cube l = Load and display radiance data cube a = Adjacency algorithm setup (create averaged image) q = Quit to main menu Processing Commands: i = simulate radiance Image from reflectance v = aerosol Visibility determination k = Kaufman (ratio) visibility retrieval c = Column water retrieval MSc = User-specified Column water M = prepare MODTRAN run R = run MODTRAN WC = Channel Wavelength Correction Utility Commands: ? = menu help p = Plot pixel spectra e = display Existing reflectance file script command M Select option: g = standard GUI interface v = visibility setup script command call Estimated MODTRAN run time: 0.035526826 min ACC error... IDL error: -361 Unable to apply format code F to input: "**********". in context: acc_called_t5set:read tape5 file back in to check format. calls: ACC_CALLED_T5SET ACC_HYPER2 ACC_ACC ENVI_ACC_EVENT ENVI_FLAASH_PROCESS_MANAGER IDLRTMAIN $MAIN$

    MariM



    Veteran Member


    Posts:2396
    Veteran Member


    --
    29 Jun 2016 09:55 AM
    I agree that the spectral filter function looks reasonable. I asked one of our engineers and they agreed that the error is typically caused by reading an unexpected ascii string. The only thing that might be unexpected that I can see is the sensor elevation. Unfortunately, I do not have access to any GOCI data to test myself. If you set the sensor to something like 700km, do you get the same error?

    Deleted User



    New Member


    Posts:
    New Member


    --
    29 Jun 2016 07:56 PM
    No,this time it works well when I set the sensor altitude to 700km though the spectrum curve is absolutely wrong.The reflectance is too high in the blue(412nm and 443nm).

    MariM



    Veteran Member


    Posts:2396
    Veteran Member


    --
    30 Jun 2016 10:24 AM
    The developers thought that using the lower altitude should not make much difference since there is little atmosphere beyond 100km. You can also turn off adjacency correction as it is negligible with large pixel sizes. I was able to find a GOCI dataset and process it in FLAASH. I did not think the spectra looked unusual: How did you convert to radiance?
    You are not authorized to post a reply.