[R-sig-Geo] can't load rgdal on mac osx [SEC=UNCLASSIFIED]

2012-04-12 Thread Hafi, Ahmed
Dear list members, I am running R on a macbook pro with mac osx lion and trying to do some spatial modelling with R. I have run into a problem in loading rgdal. I understand that the installation of rgdal in mac is different from other packages and followed the instructions given in CRAN rgdal

[R-sig-Geo] gain/offset for stack?

2012-04-12 Thread Jonathan Greenberg
r-sig-geo'ers: Does anyone know if I can set a per-layer gain and offset for a raster brick or stack (vs. a raster)? --j -- Jonathan A. Greenberg, PhD Assistant Professor Department of Geography and Geographic Information Science University of Illinois at Urbana-Champaign 607 South Mathews Aven

Re: [R-sig-Geo] spgrass6 with GRASS 7

2012-04-12 Thread Roger Bivand
On Fri, 13 Apr 2012, Paulo van Breugel wrote: Hi, I am working on Ubuntu amd64, R 2.15.0, GRASS 7.0 I am trying to import a GRASS raster layer in R, but I am getting the following error message: Which version of spgrass6? CRAN version 0.7-10 released March 18, 2012 should have resolved this

[R-sig-Geo] spgrass6 with GRASS 7

2012-04-12 Thread Paulo van Breugel
Hi, I am working on Ubuntu amd64, R 2.15.0, GRASS 7.0 I am trying to import a GRASS raster layer in R, but I am getting the following error message: > mydata <- readRAST6(vname=c("map1","map2"), useGDAL=F)Command: r.info Description: Outputs basic information about a raster map. Keywords: raster

Re: [R-sig-Geo] Overlaying SpatialGridDataFrames over SpatialGridDataFrames? (Edzer Pebesma)

2012-04-12 Thread Jooil Kim
> does > > sPDF_grid <- overlay(as(sGDF_in, "SpatialPointsDataFrame"), sGDF_ref) > > work? Until a couple of months ago, this coercion would be automatic. No, it now returns an error as such: "Error in function (classes, fdef, mtable) : unable to find an inherited method for function

Re: [R-sig-Geo] Problems with NA values and lm() in calc() in raster package

2012-04-12 Thread Lars Dalby
Hi Robert Thanks for your reply! Your suggestion solved the problem, because in my case I do have cells where there is NA in all layers. But, I'm still a little puzzled about why it is not running in the case where there is NA in a cell in only one of the layers. lm() normally works in a case li

Re: [R-sig-Geo] Points to raster

2012-04-12 Thread Raphael Saldanha
Thanks Michael! This is my first approach for dealing with rasters inside R. On Thu, Apr 12, 2012 at 11:37 AM, Michael Sumner wrote: > The sp and raster packages both have pathways for this. > > library(sp) > coordinates(elev) <- c("x", "y") > gridded(elev) <- TRUE > > See vignette("sp") > > Or,

Re: [R-sig-Geo] Points to raster

2012-04-12 Thread Michael Sumner
The sp and raster packages both have pathways for this. library(sp) coordinates(elev) <- c("x", "y") gridded(elev) <- TRUE See vignette("sp") Or, using raster package: library(raster) r.elev <- rasterFromXYZ(elev) In either case, if your coordinates aren't quite regular you'll need to explore

[R-sig-Geo] Points to raster

2012-04-12 Thread Raphael Saldanha
Hi! Which is the best way to convert a list of points to a raster? Bellow I have a excerpt from my 'elev' object. x yelev 239937.9 7658988 407.079 239937.9 7658968 411.561 239937.9 7658948 415.924 239937.9 7658928 415.480 239937.9 7658908 414.896 239937.9 765 408.602 2399

Re: [R-sig-Geo] Overlaying SpatialGridDataFrames over SpatialGridDataFrames?

2012-04-12 Thread Edzer Pebesma
does sPDF_grid <- overlay(as(sGDF_in, "SpatialPointsDataFrame"), sGDF_ref) work? Until a couple of months ago, this coercion would be automatic. On 04/12/2012 03:00 AM, Jooil Kim wrote: > Hello list, > > I've run into a baffling problem, and I hope someone can help. > > The last few months I