PrevPrev Go to previous topic
NextNext Go to next topic
Last Post 15 Feb 2017 05:58 AM by  Liam Clark
Problems using FLAASH for Worldview-3 data
 7 Replies
Sort:
You are not authorized to post a reply.
Author Messages

Liam Clark



New Member


Posts:
New Member


--
15 Feb 2017 05:58 AM
    I am attempting to use the FLAASH module on WV3 imagery but get an error saying 'Input image has more bands (8) than expected for the selected sensor (16)...' I have already used the Radiometric Calibration module to convert the raw WV3 imagery, and selected the 'FLAASH options' button when doing so. I'm really stuck on interpreting the error message as FLAASH correctly identifies the imagery as WV3 and automatically fills in the fields, so any help would be much appreciated! There's also the separate issue that I would ideally like to pansharpen the multispectral imagery using the panchromatic band, but FLAASH cannot correct the pan band. Is it acceptable to undertake pansharpening first and to then run FLAASH? I'm wary of doing this as sharpening will inevitably change the imagery.

    Mari Minari



    Veteran Member


    Posts:1968
    Veteran Member


    --
    15 Feb 2017 07:29 AM
    This is a known issue in FLAASH when detecting the sensor type. Try setting the sensor to Unknown, then manually set the sensor altitude. Next, in the Multispectral Parameters, set the filter function to the Worldview3.sli from the \ENVI5.x\classic\filt_func directory. Set the index to start at 9 (where the SWIR bands start). This should allow it to continue to run.

    Mari Minari



    Veteran Member


    Posts:1968
    Veteran Member


    --
    15 Feb 2017 07:38 AM
    With regards to using the single Pan band in FLAASH, single band files are not supported in FLAASH. You could possibly stack the pan band with the VNIR or SWIR bands but it would require you resample one of the data sets up or down. It may also require that you create a custom filter function depending on the bands that are stacked. I would not use ACM on the pan sharpened results if you want to retain the spectral integrity of the MS bands.

    Liam Clark



    New Member


    Posts:
    New Member


    --
    15 Feb 2017 09:53 AM
    Thanks Mari for the quick reply, extremely helpful!

    Liam Clark



    New Member


    Posts:
    New Member


    --
    21 Feb 2017 07:20 AM
    Sorry to bring this back up, but I've tried to follow your instructions this morning and now have different errors. There are three error messages popping up when I click 'Apply': - 'Illegal subscript range: CHANNEL_DATA' Same as above 'Variable is undefined: VIZ'. Any help would be greatly aprpeciated.

    Mari Minari



    Veteran Member


    Posts:1968
    Veteran Member


    --
    21 Feb 2017 07:35 AM
    Can you post the contents of template.txt that is created in the output directory specified in the flaash dialog?

    Liam Clark



    New Member


    Posts:
    New Member


    --
    21 Feb 2017 07:53 AM
    ; ;ENVI FLAASH PARAMETERS TEMPLATE (5.4) ;Written Tue Feb 21 14:20:56 2017 ; ; Project Parameters enviacc.prj.radiance_file = R:\NERC-STEED\WorldView\2016_radiometric_calibration\MS\16AUG14024806-M2AS_R1C1.dat enviacc.prj.reflect_file = R:\NERC-STEED\WorldView\2016_corrected\16AUG14024806-M2AS_R1C1_FLAASH enviacc.prj.filter_func_file = C:\Program Files\Harris\ENVI54\classic\filt_func\worldview3.sli enviacc.prj.filter_func_file_index = 9 enviacc.prj.water_band_choice = 1.13 enviacc.prj.red_channel = 5 enviacc.prj.green_channel = 3 enviacc.prj.blue_channel = 1 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.user_stem_name = 16AUG14024806-M2AS_R1C1_FLAASH_ enviacc.prj.modtran_directory = R:\NERC-STEED\WorldView\2016_corrected\FLAASH_process_files\ ; ; MODTRAN Parameters enviacc.modtran.visvalue = 40.0000 enviacc.modtran.f_resolution = 15.0000 enviacc.modtran.day = 14 enviacc.modtran.month = 8 enviacc.modtran.year = 2016 enviacc.modtran.gmt = 2.8017 enviacc.modtran.latitude = 4.9752 enviacc.modtran.longitude = 117.7841 enviacc.modtran.sensor_altitude = 617.0000 enviacc.modtran.ground_elevation = 0.0000 enviacc.modtran.view_zenith_angle = 180.0000 enviacc.modtran.view_azimuth = 0.0000 enviacc.modtran.atmosphere_model = 1 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 = 4096 enviacc.img.nlines = 4096 enviacc.img.data_type = 4 enviacc.img.margin1 = 0 enviacc.img.margin2 = 0 enviacc.img.nskip = 0 enviacc.img.pixel_size = 1.2000 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.output_diagnostic_files = 0 enviacc.ana.auto_save_template = 1 enviacc.ana.polishing_res = 0 enviacc.ana.aerosol_retrieval = 1 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 = [ 425.0000, 480.0000, 545.0000, 605.0000, 660.0000, 725.0000, 832.5000, 950.0000] enviacc.spc.fwhm = [ 50.000000, 60.000000, 70.000000, 40.000000, 60.000000, 40.000000, 125.000000, 180.000000] enviacc.img.p_input_scale = [ 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000]

    Mari Minari



    Veteran Member


    Posts:1968
    Veteran Member


    --
    21 Feb 2017 08:19 AM
    I am sorry. I understood this error was only occurring with the SWIR bands of WV-3 but it appears it is also happening with the VNIR bands. In this case, do the same steps as before - set the sensor to Unknown and select the WV3.sli but set the index to be 1. I will update the current bug report to reflect that this occurs with the VNIR bands as well as the SWIR.
    You are not authorized to post a reply.


    MOST RECENT

    Democratization of Geospatial Technology Encouraging Active Collaboration – Robert Laudati, Harris Geospatial

    3/21/2018

    Geospatial World – February 2018 – Rob Laudati of Harris explains the ongoing evolution of... more »

    Use Multispectral Imaging to Boost Crop Health: From Data Collection to Analysis

    3/16/2018

    The secret is out! Through the use of multispectral imagery, it is now possible to tailor crop management... more »

    SIGN UP & STAY INFORMED

    Sign up to receive the latest news, events, technologies and special offers.

    SIGN ME UP