ASTER Processing

Description

This procedure is designed to help the user process an ASTER VIS/SWIR/TIR data and run it through "standard processing", correcting for various noises and producing an emissivity array.

For a brief description of the standard processing steps, including when and how they should be used please see the paper:

Cooper, F. J., B. A. Adams, C. S. Edwards, and K. V. Hodges (2012), Large normal-sense displacement on the South Tibetan fault system in the eastern Himalaya, Geology, doi: 10.1130/g33318.1.

Procedure

Due to the format HDF4 that standard ASTER products from archives (e.g. glovis.usgs.gov) the data must be read into davinci specially using the load_aster() function which relies on GDAL to do the data translation. By default we are going to process ASTER L1B data that contain all VIS/SWIR/TIR bands. Other data types can be processed (e.g. 09, 05 etc.) but the details for these different products are not included in this discussion.

dv>a = load_aster("AST_L1B_00303202005145311_20121120130103_24164.hdf")
struct, 7 elements
  anc: struct, 60 elements
  ul: 2x1x1 array of double, bsq format [16 bytes]
  ur: 2x1x1 array of double, bsq format [16 bytes]
  ll: 2x1x1 array of double, bsq format [16 bytes]
  lr: 2x1x1 array of double, bsq format [16 bytes]
  bandlist: Text Buffer with 15 lines of text
  data: struct, 4 elements
       vnir: 4980x4200x3 array of float, bsq format [250,992,000 bytes]
       vnir3B: 4980x4600x1 array of float, bsq format [91,632,000 bytes]
       swir: 2490x2100x6 array of float, bsq format [125,496,000 bytes]
       ir: 830x700x5 array of float, bsq format [11,620,000 bytes]

Continue adding more sections that correspond to clean up (aster_deplaid, rrmnoise_pca, thresholding, NDVI, etc.)

Editing tricks Headers look like this

boxes for code examples look like this
you use a single indent

Images can be added and their size specified like this..then click the link to upload the file. 300px

internal links look like this I LOVE ASTER

external links look like this mars.asu.edu

And equal signs need to be coded in HTML (really annoying).

The code for that is =