enabled overwrite and fixed input check typo
added GRASS/Python implementation of terrain noise calculation
simplified code by using 'a' flag for r.neighbors (improves 3abda63)
reorganized AML terrain scripts to reduce superficial redundancy
renamed script to multiscalesmooth.R
numerous comment, format, and variable name changes
now actually refining before smoothing (GRASS version)
reversed order of refine-and-smooth step
added GRASS/Python implementation of multiscale smoother
added R script to run multiscalesmooth example on Oregon SRTM
switched to more accurate chisq test
enhanced smoother to accept sd as a raster, not just a constant
tweaked var names to match Gallant pub; added comments
fixed vwg step
continued R translation of multiscale smoother
initial R translation of multiscale smoother
added old Oregon topo AML scripts (Tien Ming Lee, John Gallant)
finished organisms->layers content migration
folded Nunokawa backup file changes into vrt script file
cleaned up file permissions
added in various previously uncommitted scripts
reorganized terrain-related files
removed obsolete/preliminary DEM 60N assessment scripts
removed obsolete/preliminary DEM 60N mosaicing scripts
removed obsolete/preliminary DEM regional mosaicing scripts
added old R script to assess our SRTM holdings
tweaked aster-check script to work from anywhere
added R script to compare ASTER tile names to locations
added PDF of boundary assessment doc sent to the group
final minor edits before posting to the group
just minor tweaks to code comments
added histogram plot of ASTER-SRTM differences
removed placeholder cruft from flowdir correlation plots
spell-checked boundary analysis doc (aspell)
more boundary analysis updates, especially discussion of stats
added correlations/RMSE discussion, plus many new edits
numerous changes and revisions to boundary assessment doc
switched plot to bin by SRTM elev rather than ASTER
added ASTER vs SRTM elevation comparision plots
revamped aspect/flowdir boundary assessments to use circular statistics
split the dem-fusion code snippets into separate script files
minor cleanup of code comments
added multiresolution spline fused layer as terraflow input
added prose text to boundary document, and (partly) refactored figures
early version of 60N boundary assessment report
updated boundary assessment code to use multires spline results
fixed critical typo in post-processing enblended layer
added GRASS commands to run terraflow on Canada boundary DEMs
added R code to create various boundary assessment figures
added R code to 'enblend' DEMs (multires splines, Burt & Adelson, 1983)
superficial tweaks to DEM fusion R code
added code to profile aspect near 60N boundary (just like slope code)
separated out code to profile slope near 60N boundary
added GRASS/R code to run/compare flow analysis on Oregon test area
added GRASS terraflow code plus basic R diagnostics for 60N test case
added slope boundary descrepancy plotting code
added GDAL/R code for fusing and correcting DEMs near Canadian border
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.rto 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 plotsthat 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.rCreated 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.