X
PrevPrev Go to previous topic
NextNext Go to next topic
Last Post 27 Feb 2019 12:23 PM by  MariM
batch reproject GLT with bowtie correction
 3 Replies
Sort:
You are not authorized to post a reply.
Author Messages

Ben



New Member


Posts:2
New Member


--
27 Feb 2019 09:04 AM
    Hi,

    I was wondering if anyone has experience in using DIL script to batch "Reproject GLT with Bowtie Correction" function? specifically dealing with GOES-16 L1B radiance data?
    I downloaded a bunch of GOES-16 ABI L1B radiance data (CONUS) from NOAA (OR_ABI_L1b_RadC_*.nc). In ENVI , I was able to open the files one by one and reproject the raster with GLT information. But have not had much luck to get a script working in IDL to batch process many files. Tried the following page but did not work:
    https://www.harrisgeospatial.com/docs/envireprojectglttask.html

    any advice will be appreciated!


    MariM



    Veteran Member


    Posts:2396
    Veteran Member


    --
    27 Feb 2019 10:47 AM
    What did not work? What error did you receive?

    Ben



    New Member


    Posts:2
    New Member


    --
    27 Feb 2019 11:48 AM
    I tried to follow the example from
    https://www.harrisgeospatial.com/docs/envireprojectglttask.html

    e = ENVI()

    infile = 'OR_ABI-L1b-RadC-M3C01_G16_s20181062057233_e20181062100006_c20181062100050.nc'
    Raster = e.OpenRaster(infile)

    Task = ENVITask('ReprojectGLT')

    Task.Input_Raster = Raster[0]
    Task.Latitude_Raster = Raster[1]
    Task.Longitude_Raster = Raster[2]

    then the error message appears
    % Attempt to subscript RASTER with <INT ( 2)> is out of range.

    if you can provide an example code/script showing how to reproject GLT with GOES-16 image (instead of VIIRS data as on the webpage above), that would be great

    MariM



    Veteran Member


    Posts:2396
    Veteran Member


    --
    27 Feb 2019 12:23 PM
    Specifying the Raster[0], [1] and [2] depends entirely on how many rasters are opened when you use OpenRaster. You would need to do a help on the raster object and determine if indeed, raster[0] is your radiance data, raster[1] is latitude, raster[2] is longitude, etc.

    Were you able to open the .nc file with Open As->Scientific Formats which opens the Scientific Browser? If so, then you can save this as a template and when you use OpenRaster, use your template to open the file so you know which rasters are in your file and you are passing the correct information to the GLT.
    You are not authorized to post a reply.