Changes between Version 30 and Version 31 of Procedures/EagleHawkProcessing

Jul 4, 2011, 10:58:58 AM (13 years ago)



  • Procedures/EagleHawkProcessing

    v30 v31  
    33This guide describes how to process hyperspectral data using the apl suite. This should be used for 2011 data onwards. If you are processing earlier data, see [wiki:Procedures/AZEagleHawkProcessing here] for instructions for processing with the az suite.
    5 The apl suite consists of the following programs: apl cal, apl nav, apl tran, apl map.
     5The apl suite consists of the following programs: aplcal, aplnav, apltran, aplmap.
    77== Creating config file ==
    4242=== Individual processing stages: ===
    44 ==== Stage 1: Radiometric Calibration ====
    45 The software that performs this is called aplcal. It uses the cal files to calibrate the raw data to level 1b. It strips off the FODIS region (also calibrated to level 1b), "flips" the data such that:
    46  *  eagle runs low->high wavelengths
    47  *  hawk mirrored so that it is comparable to eagle (i.e. things on the left of the aircraft are on the left of the data file)
    48 It also applies smear correction to the eagle data. Rather than mask out bad pixels it creates a separate mask file that contains flags for the following:
    49  *  overflows
    50  *  underflows
    51  *  smear affected
    52  *  bad ccd pixel
    53  *  dropped scans
    54 Example command to calibrate a file:
    55    `aplcal -input VNIR11092-1.raw -output e09201b.bil -calfile calibration/2011/eagle/SN001 `
    56 Example command to calibrate first 2501 lines of file:
    57    `aplcal -input VNIR11092-1.raw -output e09201b.bil -calfile calibration/2011/eagle/SN001 -lines 0 2500`
     44You shouldn't have to worry about this unless something goes wrong. However something often does! Detailed explanation of each step is explained [wiki:Procedures/AplSuiteDetails here]
    59 ==== Stage 2: Navigation syncing ====
    60 This is performed using the aplnav software. It can use navigation data either from an SBET or specim nav file, depending on whether you want to use post-processed or real-time data.
    61 It is possible to add a timing offset to shift the navgation data w.r.t the scan lines, or to shift the position w.r.t to attitude. It would be rare to require the latter as this suggests a problem with the SBET file.
    63 Lever arms and boresight values are entered here to offset the position and attitude to the eagle/hawk sensor view direction. It is possible to smooth the data (using a triangle filter) and interpolate the scan time data using either linear or spline methods.
    64 BIL files containing quality flags for the navigation data are output also.
    66 '''IMPORTANT NOTE:'''
    68  * Lever arms are entered as X Y Z in the aircraft coordinates (X +ve forward, Y +ve to starboard, Z +ve down). This is different to using aznav, which used X +ve forward, Y +ve port, Z +ve up.
    69  * Specim have previously said that there is a -0.055 second offset to apply to the sync message. Rather than hard code this into aplnav, you should add it onto the scantimeoffset parameter. This is automatically done using the processing scripts by the sct_global_offset parameter.
    71 Example command to sync post-processed navigation data:
    72    `aplnav -nav VNIR092-11-1.nav -sbet sbet_2011092.out -lev1 e092011b.bil -interp Spline -leverarm 0.559 0.015 1.543 -boresight -0.28 0.16 0.38 -scantimeoffset 0.045 -output e092011b_p_sct0.1_nav_post_processed.bil -qualityfile e092011b_nav_post_processed_quality.bil`
    74 ==== Stage 3: Geocorrection ====
    75 The first step of geocorrection is to generate the IGM file using aplcorr. This takes an optional 1-band BIL / BSQ DEM (in WGS-84 Lat/Lon) and the instrument FOV vector file.
    77 Example command:
    78    `aplcorr -lev1file e092011b.bil -vvfile eagle_fov_fullccd_vectors.bil -navfile e092011b_p_sct0.1_nav_post_processed.bil -dem ASTER_resampled_0.1667_arcsec.dem  -igmfile e092013b_p_sct0.1.igm`
    80 The second step of geocorrection is to re-project the IGM file into the chosen map projection, using apltran.
    82 Example commands:
    83    `apltran -inproj latlong WGS84 -igm e092013b_p_sct0.1.igm -outproj utm_wgs84N 31 -output e092013b_p_sct0.1_utm31n.igm`
    85    `apltran -inproj latlong WGS84 -igm e092013b_p_sct0.1.igm -outproj osng /users/rsg/arsf/dems/ostn02/OSTN02_NTv2.gsb -output e092013b_p_sct0.1_utm31n.igm`
    87 The third step is then to map it.
    89 Example commands:
    90    `aplmap -igm e092013b_p_sct0.1_utm31n.igm -lev1 e092011b.bil -mapname e092013b_utm31.mapped -bandlist 30 15 7`
    92    `aplmap -igm e092013b_p_sct0.1_utm31n.igm -lev1 e092011b.bil -mapname e092013b_utm31.mapped -bandlist 30 15 7` -pixelsize 2 2
    94    `aplmap -igm e092013b_p_sct0.1_utm31n.igm -lev1 e092011b.bil -mapname e092013b_utm31.mapped -bandlist 30 15 7` -oversample 3 3
    97 === Making a Delivery ===
     46=== Making a delivery ===
    9948Use the script to make the delivery directory. Run it from within the main project directory. By default it runs in dry run mode.