Set of shell scripts used to create prototype terrain mosaic layers from CGIR SRTM and ASTER GDEM data tiles.
Two prototypes generated: 3 ArcSecond and 30 ArcSecond resolution. There are some minor inconsistencies in these scripts; however, the mosaic layer components that they generate appear to be valid, except for some...

First version of this script; dKc constant is not optimally set yet. Need more experimentation.

Version of pixel pair plotting script that adds 'Lowess Lines' to center of each scatterplot.

Note: Loess lines feature has only been tested with the input file "pixelPairs36000_5_8EvenSortCol1.csv", used in the script.
Some table column names are hard-coded into the script.

This version has improved, vertically-stacked plots.

First version of 'image blending' solution that employs a linear decay function to blend ASTER and SRTM pixel values south of 60 degree N Latitude line.

First version of ASTER/SRTM boundary adjustment R script that tests a prototype 'exponential decay' function on a sub image of the DEM imaage mosaic.

First version of R script that adjusts the boundary between ASTER and SRTM components of DEM image mosaic created for Environment and Organisms project. This prototype version uses a set of 'narrow' test images, applies a linear ramp function. Later versions will use more sophisticated functions.

this table created using the latest and best-aligned versions of the 'boundary validation' (ASTER/SRTM/CDEM) subimages. Input to the script: CreateImageDiffPlotsOverlayNewPlots.r

Added option: write individual plots (not trplets - yet) to separate JPG files. New command line argument: JpgPlotFileFlag

Table used by CreateImageDiffPlotsOverlayNewPlots.r
to create latest set of (normalized) difference scatterplots

This version used to create the large (36000 row) table used to generate latest bounday analysis scatterplots.

This version creates the (normalized) difference image plots
that Rick, Jim, and Mark are reviewing as part of the final 'ASTER/SRTM boundary analysis' work.
The input file pixelPairs36000_5_8EvenSortCol1.csv also checked in.

The first version of 'CreateImageDiffPlotsOverlay.r' that produces the 'three-panel' plot containing the boundary differences for ASTER / Border / SRTM.

Latest version of Elevation Difference Image Analysis plot generator: Plots remain on screen until end of script (user presses 'final' <cr>

TODO: Add legend to plots.

correct CDEM image values

With correct CEDM input image filename.

Uses correct, even boundary input images

made with latest version of makeImagePairTable w/correct images.

created with latest version of code 5/15

Tested Sunday 5/15 - produces table tableForRic2000_5_8.csv

In header, better 'to run' instructions.

input table for the R script: MosaicCdemDifferencePlots.r
Created by the R script: makeImagePairFiles.

First working version - very simple - just 'source' the file and press <CR> to walk through the plots.

Next version: 'panel plots' and improved plot axes.

Very first version of plot making script. Plots need consolidation into 'panel plots', and plot labels need to be improved.

Uses input file:
tableForMark4000_5_8_SortColID.csv, also checked in.

These scripts are complete but may not run correctly due to conflicts in the extents of the input raster files. Recheck inputs.

Version runs correctly, though the CDEM edge input image might not be the latest and best. Recheck inputs

ReadMe file describes how to use shell scripts (also checked in) to create example Srtm/Aster DEM image mosaic of Western Canada region spanning 60 deg N Latitude boundary

Sample output file from R script makeMarkTable.r - sample of 4000 random columns, row 5 through row 8 straddling 60 degree boundary.

Producing valid table, but from column samples taken entirely from mosaic. Next step: take them from Aster and Cgiar images.

First version, has issues with column indexing, I am fixing tonight.

Added collection of Linux shell scripts and R scripts that I used in my initial analysis of SRTM/CGIAR, CDEM, and ASTER GDEM digital terrain data sets.

added empty dir structure for terrestrial work