[R-sig-Geo] Three-dimensional interpolation of temperature profiles using gstat?

2014-06-03 Thread Trevor Joyce
Hi, I'm working on a problem in which I'd like to interpolate a set of oceanographic temperature profiles over a three-dimensional matrix of unsampled x, y, & z locations based on a SpatialPointsDataFrame (WOD) in which I've stored temperature data (WOD$TEMP) in association with coordinate informa

Re: [R-sig-Geo] Export object class SpatialPixelsDataFrame to GIS

2014-06-03 Thread Sarah Goslee
Hi, You can use writeGDAL() from the rgdal package. Sarah On Tue, Jun 3, 2014 at 4:24 PM, Hazlehurst, Jenny A wrote: > ?Thanks in advance. Using the package adehabitatHR for kernel analysis of > bird home ranges. > > > Looking for how to get a specific kernel, which is of class > SpatialPixel

[R-sig-Geo] Export object class SpatialPixelsDataFrame to GIS

2014-06-03 Thread Hazlehurst, Jenny A
?Thanks in advance. Using the package adehabitatHR for kernel analysis of bird home ranges. Looking for how to get a specific kernel, which is of class SpatialPixelsDataFrame, into a format use-able by GIS platform, mostly QGIS, so that I can work within specific kernels in the field. Thanks

Re: [R-sig-Geo] Adding smoothing factor to terrain function in raster

2014-06-03 Thread Alex Zvoleff
Whoops - somehow a double comma worked its way into the ."Call" lines in the first patch I sent. The attached patch is updated to address that issue. Alex On Fri, May 30, 2014 at 10:48 AM, Alex Zvoleff wrote: > Dear list (and Robert Hijmans in particular), > > I am using the terrain function in

Re: [R-sig-Geo] Failure during raster IO

2014-06-03 Thread Ariel Ortiz-Bobea
Ran into a similar problem and was getting "Failure during raster IO" when using extract on a very large raster with a small polygon. My solution was to crop the large raster before doing the extract. Maybe this can help other folks in the future. Cheers. Ariel - Ariel Ortiz-Bobea Fellow at

Re: [R-sig-Geo] Cut SpatialLines into two Lines when intersecting with SpatialPolyon

2014-06-03 Thread Roger Bivand
On Tue, 3 Jun 2014, Wouter Steenbeek wrote: I'm sorry I didn't provide a toy example. I figured it was a simple question and I just couldn't find some well-known function. I am not very familiar with readWKT(), but I used the sp_intro text to generate the following example: library(sp) Sr1

Re: [R-sig-Geo] Cut SpatialLines into two Lines when intersecting with SpatialPolyon

2014-06-03 Thread Wouter Steenbeek
I'm sorry I didn't provide a toy example. I figured it was a simple question and I just couldn't find some well-known function. I am not very familiar with readWKT(), but I used the sp_intro text to generate the following example: library(sp) Sr1 = Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2))) Sr2

Re: [R-sig-Geo] Error in extract function in Raster package on Mac OS X

2014-06-03 Thread Justin Michell
All fixed, thanks. I needed to install the "Spring Dance” version of R. Not quite sure what was wrong previously, but it works now! Cheers Justin On Jun 3, 2014, at 3:20 PM, Alaba Boluwade wrote: >> r <- raster(ncol=36, nrow=18) >> r[] <- 1:ncell(r) >> >> ### >

[R-sig-Geo] Error in extract function in Raster package on Mac OS X

2014-06-03 Thread Justin Michell
Dear List I have encountered a problem with the ‘extract’ function in Raster package, previously it worked for me but not anymore. Here is an example from the Raster package documentation. It does not work on my Mac OS X 10.9.3 running R version 3.0.2 (2013-09-25) -- "Frisbee Sailing” (64 bit).

Re: [R-sig-Geo] Cut SpatialLines into two Lines when intersecting with SpatialPolyon

2014-06-03 Thread Roger Bivand
On Tue, 3 Jun 2014, Wouter Steenbeek wrote: [My apologies if this is a double post. I wrote the first post in HTML, and I think it may not have gotten through.] Thanks for respecting our guidelines - it is much harder to infect plain text, which is why HTML is not welcome. I have a Spatia

[R-sig-Geo] Cut SpatialLines into two Lines when intersecting with SpatialPolyon

2014-06-03 Thread Wouter Steenbeek
[My apologies if this is a double post. I wrote the first post in HTML, and I think it may not have gotten through.] I have a SpatialLines(DataFrame) (e.g. streets) and SpatialPolygons(DataFrame) (e.g. neighborhoods). When a Line intersects with a polygon boundary, I want to split this Line int

Re: [R-sig-Geo] (Linear) interpolation in three (or more) dimensions

2014-06-03 Thread Edzer Pebesma
package gstat offers simple 3D interpolation (idw, kriging), e.g. library(gstat) example(gstat3D) and also spatio-temporal kriging; some memory limitations there may be dealt with by using local kriging (see ?krigeST). On 06/03/2014 08:27 AM, Mark Payne wrote: > Hi, > > I have a dataset that co