[R] smoothing bathymetry
currie.jens at gmail.com
Wed Apr 6 20:01:07 CEST 2011
Dear R Users,
Using the following R-script I created the first image
> dep <- interp(long, lat, depth, xo=seq(1,990,10), yo=seq(1,990,10),
+ extrap=FALSE, ncp=0,duplicate = "mean", dupfun = NULL)
Where "long" are x-coordinates between 1 and 1000, "lat" are y-cordinates
between 1 and 1000, and depth are depth values corresponding to the x and y
coordinates. All data are in vector form.
I would like to extrapolate the data points beyond the data set, by setting
extrap=TRUE, filling out the entire 1000x1000 grid.
I attempted this using the following script, which generated the second
>dep <- interp(long, lat, depth, xo=seq(1,990,10), yo=seq(1,990,10),
+ linear=FALSE,extrap=TRUE, ncp=NULL,duplicate = "mean", dupfun = NULL)
However, as seen in the image, this results in a messy looking bathymetry,
with erros. Any suggestions on what went wrong or other packages or
functions I could try to would be greatly appreciated?
I played around with the krig.image function in the "fields" package, but
didnt have much luck.
View this message in context: http://r.789695.n4.nabble.com/smoothing-bathymetry-tp3431391p3431391.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help