[R-sig-Geo] apply a formula by matrix

2011-07-07 Thread Alfredo Alessandrini
Hi, I've to apply a formula to a raster in relation to an index table created from raster stack. I'm trying to explain the problem: I've three raster with same resolution, that I join in a raster RasterStack, and from it I can obtain a table like this: > data <- stack(dem,aer550,aer550_type)

[R-sig-Geo] how to get a SpatialGrid object from a .gri file ?

2011-07-07 Thread camille garcin
Hello I got a .gri file (data about density of population in Ireland from http://www.diva-gis.org/gdata) that I would like to read with the readGDAL function from the *rgdal* package. So have I to get a SpatialGrid object from this file. Could you tell me if it's possible to do that with the * rgd

[R-sig-Geo] spatial modelling Lectureship at the University of Bristol

2011-07-07 Thread Malcolm Fairbrother
Dear list, Please see below an announcement about a job opening in the School of Geographical Sciences at the University of Bristol. The post is only just being advertised today, but there is a need to have someone in post by October -- so there is a very short time before the deadline. Subscr

[R-sig-Geo] plot3D raster coordinates

2011-07-07 Thread Etienne B. Racine
Hi list, I'd like to plot3D() a raster layer using it's own coordinates instead of the row-col system so I can easily add points. Here's a demo library(raster) library(rgl) r <- raster(matrix(1, nrow = 5, ncol = 5)) extent(r) # class : Extent # xmin: 0 # xmax: 1 # ymin

Re: [R-sig-Geo] plot3D raster coordinates

2011-07-07 Thread Etienne B. Racine
I don't know if this will clarify or add to the confusion as I don't know if it's the same issue or another one with the same effect, but when using a projected raster, I get some weird coordinates in 3d (really large). I can't see how it is linked to my raster coordinates. Example : library(raste

[R-sig-Geo] Error in anyDuplicated.default(coords) message with trinb

2011-07-07 Thread Stratford, Jeffrey
Hi everyone, I'm attempting to do some spatial analyses with a relatively simple data set and I'm running into some issues. Here's my code: sosp09 <- read.csv("f:\\sosp\\2009\\sosp09.means.csv", header=T) sp.point <- cbind(sosp09$x, sosp09$y) colnames(sp.point) <- c("LONG","LAT") # Pro

Re: [R-sig-Geo] Error in anyDuplicated.default(coords) message with trinb

2011-07-07 Thread Edzer Pebesma
Jeff, ?tri2nb suggests you to try tri2nb(coordinates(data.sp)) On 07/07/2011 03:52 PM, Stratford, Jeffrey wrote: > Hi everyone, > > > > I'm attempting to do some spatial analyses with a relatively simple data > set and I'm running into some issues. Here's my code: > > > > sosp09 <- read.

[R-sig-Geo] Mapping contours from jpg map

2011-07-07 Thread Sarah Papworth
Dear all, I have a jpg map of altitudinal contours in my study area which I would like to map and use in R to calculate the altitude and slope of GPS coordinates. So far, I have painstakingly traced these contours in ArcMap to create a shapefile (of polylines) with altitude. I have now imported th

Re: [R-sig-Geo] Mapping contours from jpg map

2011-07-07 Thread Barry Rowlingson
On Thu, Jul 7, 2011 at 5:25 PM, Sarah Papworth wrote: > Dear all, > > I have a jpg map of altitudinal contours in my study area which I > would like to map and use in R to calculate the altitude and slope of > GPS coordinates. > So far, I have painstakingly traced these contours in ArcMap to creat

Re: [R-sig-Geo] Mapping contours from jpg map

2011-07-07 Thread Barry Rowlingson
>  I suspect there's a better way by looking at more than one nearby > contour line and doing some interpolation I've got fairly good algorithm going now, but it requires a bit of modification of some sp package code. Currently spsample can return random or regular points on a line, but do

[R-sig-Geo] Running R on a Computer Cluster in the Cloud - cloudnumbers.com

2011-07-07 Thread Markus Schmidberger
Dear analysts for geographical data and R users, cloudnumbers.com provides researchers and companies with the resources to perform high performance calculations in the cloud. As cloudnumbers.com's community manager I may invite you to register and test R on a computer cluster in the cloud for free

[R-sig-Geo] GLS model

2011-07-07 Thread Yiching Lin
Dear all, I have been using the GLS model provided by nlme package in R. I was wondering if it is possible to estimate an equivalent of R square for GLS in any R package. Nlme package did not provide such estimations. The only application that I found is Spatial Modeling Environment (SME). Tha