| 1 |
13367984
|
Jim Regetz
|
# Experimental R procedural code to run multiscalesmooth function on
|
| 2 |
|
|
# srtm_oregon binary grid produced by Ming
|
| 3 |
|
|
#
|
| 4 |
|
|
# Jim Regetz
|
| 5 |
|
|
# NCEAS
|
| 6 |
|
|
# Crated on 27-Mar-2012
|
| 7 |
|
|
|
| 8 |
|
|
library(raster)
|
| 9 |
|
|
|
| 10 |
|
|
# set up some local paths on atlas
|
| 11 |
5af6ef96
|
Jim Regetz
|
src.dir <- "/home/regetz/git/environmental-layers/terrain/research/oregon"
|
| 12 |
13367984
|
Jim Regetz
|
data.dir <- "/home/layers/commons/topo/v2/experimental"
|
| 13 |
|
|
|
| 14 |
5af6ef96
|
Jim Regetz
|
source(file.path(src.dir, "multiscalesmooth.R"))
|
| 15 |
13367984
|
Jim Regetz
|
|
| 16 |
|
|
srtm_oregon <- raster(file.path(data.dir, "srtm_oregon"))
|
| 17 |
|
|
|
| 18 |
|
|
# create a subset of data to work with
|
| 19 |
|
|
srtm.sub <- crop(srtm_oregon, extent(srtm_oregon, r1=4001, r2=6000,
|
| 20 |
|
|
c1=4001, c2=6000))
|
| 21 |
|
|
|
| 22 |
|
|
srtm.sub.sm <- multiscalesmooth(srtm.sub, sd=0.001)
|
