[R-sig-Geo] autoKrige: singular matrix in function LDLfactor

2013-08-05 Thread ElenaC
Dear List, I have a dataset (with an irregular grid) which I am trying to interpolate to a regular grid with much higher spatial resolution using the autoKrige function from the package 'automap'. But I get the error: > SST_data <- autoKrige(SST~1, HG1990, center_points2, nmax=50) [using ordinary

Re: [R-sig-Geo] Joining polygon and raster attributes

2013-08-05 Thread Tom Philippi
something like: test1 <- extract(dryas,teow,na.rm=TRUE,sp=TRUE,fun=function(x) sum(x>0)) might give you what you ask: the area (count of raster cells, so in units of cell areas) in each polygon of teow with raster values > 0. If dryas has small cells, relative to the size of the ecoregion polyg

[R-sig-Geo] Joining polygon and raster attributes

2013-08-05 Thread Empty Empty
[I apologize for reposting this – it showed up in the archive under the wrong subject, which would make it less useful as an archive.] Hi. I’m trying to develop a repeatable process in R that I can use with a variety of data. Basically I'm trying to join/merge/combine/intersect (I’m not sure

Re: [R-sig-Geo] ggmap route planning

2013-08-05 Thread Raphael Saldanha
Hi Oscar, thanks! --- Regards, Raphael Saldanha rfsalda...@outlook.com Date: Mon, 5 Aug 2013 15:16:06 +0200 Subject: Re: [R-sig-Geo] ggmap route planning From: oscar.perpi...@gmail.com To: saldanha.plan...@gmail.com CC: r-sig-...@stat.math.ethz.ch Hello,Try the osmar package: http://osmar.r-fo

[R-sig-Geo] Joining polygon and raster attributes

2013-08-05 Thread Empty Empty
Hi. I’m trying to develop a repeatable process in R that I can use with a variety of data. Basically I'm trying to join/merge/combine/intersect (I’m not sure technically which word is right) a polygon and a raster so that I can do summaries on attributes from both. I have A) a polygon shape file

[R-sig-Geo] raster/rgdal- problem: Too many open files (Linux)

2013-08-05 Thread Jon Olav Skoien
Dear list, We have a problem which appears to be a bug in either rgdal or raster, although it could also be a bug in base R or in our understanding of how to deal with connections. We have a process which is writing a rather large (~10-20.000) number of geoTiffs via writeRaster. However, the p

Re: [R-sig-Geo] ggmap route planning

2013-08-05 Thread Oscar Perpiñan
Hello, Try the osmar package: http://osmar.r-forge.r-project.org/. The section "R as navigator" of their paper includes an useful example: http://journal.r-project.org/archive/2013-1/eugster-schlesinger.pdf Best, Oscar. Oscar Perpiñán Lamigueiro Grupo de Sistemas Fotovoltaicos (IES-UPM) Dpto.