Re: [R-sig-Geo] layer name in writeRaster

2012-10-26 Thread Robert J. Hijmans
Pascal, I see that names are written to file if (1) the cell values are in memory (not, as in your case on file) or (2) if you write a RasterBrick (not as in your case a RasterLayer). Bug fixed in next version. Work-arounds: 1) write a RasterBrick instead of a RasterLayer names(x) <- 'zebra' x <

Re: [R-sig-Geo] gDistance on a LARGE Spatial Object?

2012-10-26 Thread steven mosher
have a look at my blog. i solve the problem using the big memory package. there are a couple ways to speed it up but you can just download my code. On Oct 26, 2012 11:46 AM, "Wouter Steenbeek" wrote: > I am trying to use the "gDistance" from the "geos" package to calculate > the Euclidean distanc

Re: [R-sig-Geo] gDistance on a LARGE Spatial Object?

2012-10-26 Thread Wouter Steenbeek
About 15.5 seconds here. This is awesome, thanks!! From: Roger Bivand [roger.biv...@nhh.no] Sent: Friday, October 26, 2012 11:47 PM To: Wouter Steenbeek Cc: r-sig-geo@r-project.org Subject: Re: [R-sig-Geo] gDistance on a LARGE Spatial Object? On Fri, 26 O

Re: [R-sig-Geo] gDistance on a LARGE Spatial Object?

2012-10-26 Thread Roger Bivand
On Fri, 26 Oct 2012, Roger Bivand wrote: On Fri, 26 Oct 2012, Wouter Steenbeek wrote: Disregarding solution1 (;-)), let me give you a practical example where I (think I) really need the full matrix. Or more specifically, an example in which I use the function gWithinDistance (but the same mem

Re: [R-sig-Geo] gDistance on a LARGE Spatial Object?

2012-10-26 Thread Roger Bivand
On Fri, 26 Oct 2012, Wouter Steenbeek wrote: Disregarding solution1 (;-)), let me give you a practical example where I (think I) really need the full matrix. Or more specifically, an example in which I use the function gWithinDistance (but the same memory issues apply here). I have a points.

Re: [R-sig-Geo] gDistance on a LARGE Spatial Object?

2012-10-26 Thread Wouter Steenbeek
Disregarding solution1 (;-)), let me give you a practical example where I (think I) really need the full matrix. Or more specifically, an example in which I use the function gWithinDistance (but the same memory issues apply here). I have a points.sp object containing my units of analysis (e.g.,

Re: [R-sig-Geo] gDistance on a LARGE Spatial Object?

2012-10-26 Thread Edzer Pebesma
On 10/26/2012 08:46 PM, Wouter Steenbeek wrote: > I am trying to use the "gDistance" from the "geos" package to calculate the > Euclidean distance between a large number of points: to be exact, 28000 > points. > > I'm doing: > distance <- gDistance(points.sp) > > R doesn't run this on my 32-

[R-sig-Geo] gDistance on a LARGE Spatial Object?

2012-10-26 Thread Wouter Steenbeek
I am trying to use the "gDistance" from the "geos" package to calculate the Euclidean distance between a large number of points: to be exact, 28000 points. I'm doing: distance <- gDistance(points.sp) R doesn't run this on my 32-bit 3Gb RAM machine. The same function works on a subset of these

[R-sig-Geo] [raster] unexpected write behavior with writeRaster, projectRaster

2012-10-26 Thread Tom Roche
Noticing an unwelcome similarity between writeRaster and projectRaster: https://stat.ethz.ch/pipermail/r-sig-geo/2012-October/016523.html > Here I show the raster object, I then define the name, and then I > write it to disk. You can see that the layer name switches back to > "layer". > > comb >

[R-sig-Geo] layer name in writeRaster

2012-10-26 Thread Pascal Title
Hi all, I'm working with a set of rasters in R, and writing them to disk. However, when I read them back in, the layer name is back to default. How can I define the layer name and have that carry over to the raster that is saved to disk? I'm using R 2.15.1, raster package 2.0-12, Mac OSX Mountain

Re: [R-sig-Geo] RSAGA issues

2012-10-26 Thread Jean-Daniel Sylvain
Hi, For point /1 use this command: library(rgdal) image <- "E:/IRDA/Automatisation/ScriptDevelopment/data/PredictionProprietesSol/Ref_L5_2521_IBN7c.sdat" info <- try(GDALinfo(image), silent = TRUE) For point 2: You need to add extension to your 't_codes.sgrd'. rsaga.geoprocessor(lib="g

[R-sig-Geo] RSAGA issues

2012-10-26 Thread Ludwig Hilger
Dear list, I am comparatively new to RSAGA and have come across three problems I have not been able solve by searching on the internet. Your help would be very much appreciated. This is SAGA 2.0.8 and R 2.15.1. 1. Is there a way to access the header of a *.sgrd without having to load the wh