8.17. ASTER

Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) is a Japanese instrument. ASP can process ASTER Level 1A VNIR images. These are acquired with a stereo rig consisting of two cameras, pointing nadir and back. The orbit is sun-synchronous, at an elevation of 705 km. The ground sample distance is 15 meters/pixel.

See a ready-made ASTER example. It has the input images and cameras, ASP outputs, and instructions for how to run it. Also see a workbook with illustrations.

ASP can correct for the jitter in these cameras (Section 16.37.8).

8.17.1. Fetching the data

ASTER satellite images are freely available from:

When visiting that page, select a region on the map, search for AST_L1A, and choose ASTER L1A Reconstructed Unprocessed Instrument Data V003. (The same interface can be used to obtain pre-existing ASTER DEMs.) If too many results are shown, narrow down the choices by using a range in time or deselecting unwanted items manually. Examining the data thumbnails is helpful, to exclude those with clouds, etc. Then click to download.

It is very important that, at the very last step, when finalizing the order options, choose GeoTIFF as the data format, rather than HDF-EOS. This way the images and metadata will come already extracted from the HDF file.

ASTER L1B images are also available. These are produced by projecting L1A images onto the WGS84 ellipsoid at zero elevation. ASTER L1B images can be processed with ASP by using the mapprojection workflow (Section 6.1.7). The user should invoke parallel_stereo with the L1B images (already mapprojected), L1A cameras, output prefix, and a DEM with zero height above the WGS84 datum. The results are nearly the same as obtained with L1A images.

8.17.2. Data preparation

In this example will use the dataset AST_L1A_00307182000191236_20160404141337_21031 near San Luis Reservoir in Northern California. This dataset will contain TIFF images and meta-information as text files. We use the tool aster2asp to parse it:

aster2asp dataDir -o out

Here, dataDir is the directory containing the *VNIR*tif files produced by unzipping the ASTER dataset (sometimes this creates a new directory, and sometimes the files are extracted in the current one).

This command will create 4 files, named:

out-Band3N.tif out-Band3B.tif out-Band3N.xml out-Band3B.xml

We refer again to the tool’s documentation page regarding details of how these files were created.

Open the images in stereo_gui (Section 16.67) as:

stereo_gui out-Band3N.tif out-Band3B.tif

and ensure that they are of good quality, or else get another dataset.

8.17.3. Stereo with raw images

Run parallel_stereo (Section 16.50):

parallel_stereo -t aster         \
  --stereo-algorithm asp_mgm     \
  --subpixel-mode 9              \
   out-Band3N.tif out-Band3B.tif \
   out-Band3N.xml out-Band3B.xml \
   out_stereo/run

This used the asp_mgm algorithm, which is the most accurate algorithm ASP has. One can also try the option --subpixel-mode 2 which will be much slower but produce better results.

See Section 6 for a discussion about various stereo algorithms and speed-vs-quality choices.

This is followed by DEM creation with point2dem:

point2dem -r earth --stereographic --auto-proj-center \
  out_stereo/run-PC.tif

This will create a DEM named out_stereo/run-DEM.tif using an auto-guessed local stereographic projection with auto-guessed resolution (about 15 m / pixel, the image ground sample distance). See the point2dem documentation in Section 16.55 for more details about how to run this program.

Visualize the DEM with stereo_gui (Section 16.67):

stereo_gui --hillshade out_stereo/run-DEM.tif

8.17.4. Stereo with mapprojected images

To improve the results for steep terrain, one may consider doing stereo as before, followed by mapprojection onto a coarser and smoother version of the obtained DEM, and then redoing stereo with mapprojected images (per the suggestions in Section 6.1.7).

# Initial stereo
parallel_stereo -t aster         \
  --stereo-algorithm asp_mgm     \
  --subpixel-mode 9              \
   out-Band3N.tif out-Band3B.tif \
   out-Band3N.xml out-Band3B.xml \
   out_stereo/run

# Create a low-resolution smooth DEM at 200 meters/pixel
point2dem -r earth --stereographic --auto-proj-center \
  --tr 200 out_stereo/run-PC.tif -o out_stereo/run-200m

# Mapproject onto this DEM at 15 meters/pixel
mapproject --tr 15 out_stereo/run-200m-DEM.tif        \
  out-Band3N.tif out-Band3N.xml out-Band3N_proj.tif
mapproject --tr 15 out_stereo/run-200m-DEM.tif        \
  out-Band3B.tif out-Band3B.xml out-Band3B_proj.tif

# Run parallel_stereo with the mapprojected images
parallel_stereo -t aster                  \
  --stereo-algorithm asp_mgm              \
  --subpixel-mode 9                       \
  out-Band3N_proj.tif out-Band3B_proj.tif \
  out-Band3N.xml out-Band3B.xml           \
  out_stereo_proj/run                     \
  out_stereo/run-200m-DEM.tif

# Create the final DEM
point2dem -r earth --stereographic --auto-proj-center \
 out_stereo_proj/run-PC.tif

It is very important to use the same resolution (option --tr) for both images when mapprojecting. That helps making the resulting images more similar and reduces the processing time (Section 6.1.7.4).

One could consider mapprojecting at a higher resolution, for example, at 10 meters/pixel.

It is suggested to also create and inspect the intersection error image (Section 16.55). If it is large (comparable to ground sample distance), the cameras should be bundle-adjusted first (Section 16.5).

See Fig. 16.14 for an illustration.

8.17.5. Using the CSM model

An ASTER camera model consists of a sequence of satellite position samples and a set of camera directions (sight vectors, in world coordinates), sampled at about a dozen image rows and columns. Interpolation is used in-between.

ASP can, in addition, fit a CSM linescan model (Section 8.12) on-the-fly to the ASTER model. This has the advantage that instead of a set of directions on a grid, there is one camera orientation at each satellite position sample. These will be used to solve for jitter in ASTER cameras (Section 16.37.8).

This functionality can be turned on with the option --aster-use-csm in stereo, bundle adjustment, mapprojection, and cam_test (Section 16.9). This option is implicitly assumed when solving for jitter, as that tool only works with CSM cameras.

The CSM model is produced by optimizing the optical center, focal length, and camera orientations, to fit best the provided ASTER sight vectors. No ground information is used, or stereo pair knowledge. The satellite positions do not change. This model results in a triangulated surface that is different by an average of 1 m or so vertically from the one obtained with the original cameras, but this is very small given the ground sample distance of 15 meters, and is not noticeable when taking the difference with a prior terrain model.

The cam_test documentation also describes how to compare the existing ASTER and new CSM-based implementations.

The bundle adjustment program (Section 16.5) will optimize and save the produced CSM models (Section 8.12.6), if invoked with this switch. To save the best-fit CSM models with no further refinement, invoke this tool with zero iterations.

The CSM model may be further refined by tying together multiple datasets and using ground constraints (Section 12.2.2).