Changes between Initial Version and Version 1 of Procedures/AplSuiteDetails

Jul 4, 2011, 10:56:57 AM (11 years ago)



  • Procedures/AplSuiteDetails

    v1 v1  
     1=== Details on individual processing stages: ===
     3==== Stage 1: Radiometric Calibration ====
     4The 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:
     5 *  eagle runs low->high wavelengths
     6 *  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)
     7It 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:
     8 *  overflows
     9 *  underflows
     10 *  smear affected
     11 *  bad ccd pixel
     12 *  dropped scans
     13Example command to calibrate a file:
     14   `aplcal -input VNIR11092-1.raw -output e09201b.bil -calfile calibration/2011/eagle/SN001 `
     15Example command to calibrate first 2501 lines of file:
     16   `aplcal -input VNIR11092-1.raw -output e09201b.bil -calfile calibration/2011/eagle/SN001 -lines 0 2500`
     18==== Stage 2: Navigation syncing ====
     19This 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.
     20It 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.
     22Lever 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.
     23BIL files containing quality flags for the navigation data are output also.
     25'''IMPORTANT NOTE:'''
     27 * 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.
     28 * 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.
     30Example command to sync post-processed navigation data:
     31   `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`
     33==== Stage 3: Geocorrection ====
     34The 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.
     36Example command:
     37   `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`
     39The second step of geocorrection is to re-project the IGM file into the chosen map projection, using apltran.
     41Example commands:
     42   `apltran -inproj latlong WGS84 -igm e092013b_p_sct0.1.igm -outproj utm_wgs84N 31 -output e092013b_p_sct0.1_utm31n.igm`
     44   `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`
     46The third step is then to map it.
     48Example commands:
     49   `aplmap -igm e092013b_p_sct0.1_utm31n.igm -lev1 e092011b.bil -mapname e092013b_utm31.mapped -bandlist 30 15 7`
     51   `aplmap -igm e092013b_p_sct0.1_utm31n.igm -lev1 e092011b.bil -mapname e092013b_utm31.mapped -bandlist 30 15 7` -pixelsize 2 2
     53   `aplmap -igm e092013b_p_sct0.1_utm31n.igm -lev1 e092011b.bil -mapname e092013b_utm31.mapped -bandlist 30 15 7` -oversample 3 3