|
1
|
# 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
|
src.dir <- "/home/regetz/git/environmental-layers/terrain/research/oregon"
|
|
12
|
data.dir <- "/home/layers/commons/topo/v2/experimental"
|
|
13
|
|
|
14
|
source(file.path(src.dir, "multiscalesmooth.R"))
|
|
15
|
|
|
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)
|
|
23
|
|
