[R-sig-Geo] I just signed “Shri. Manmohan Singh, Prime Minister - Government of India : Save RTI Act from Amendments ”

2013-08-06 Thread Kanwal Nayan Singh
Hi, I've just signed the following petition "Shri. Manmohan Singh, Prime Minister - Government of India : Save RTI Act from Amendments " and wanted to see if you could help by adding your name. Our goal is to reach 35,000 signatures and we need more support. You can read more and sign the peti

Re: [R-sig-Geo] R2Winbugs - problem with workign directory

2013-08-06 Thread James Rooney
Hi Virgilio, Many thanks. As it happens today I switched to OpenBugs and don't have this particular problem anymore. But thanks again, James From: Virgilio Gomez Rubio [virgilio.go...@uclm.es] Sent: 06 August 2013 17:21 To: James Rooney Cc: r-sig-geo@r-p

Re: [R-sig-Geo] can gCentroid preserve polygons data frame?

2013-08-06 Thread Tom Philippi
Joseph-- if you want a SpatialPointsDataFrame, you can use PP <- gCentroid(SP, byid=TRUE, id=XX) then glue your attribute table from SP to the spatial points via SpatialPointsDataFrame(PP,SP@data). If you don't have a suitable unique id variable in SP, you can create it first. Tom Tom On Tue,

[R-sig-Geo] can gCentroid preserve polygons data frame?

2013-08-06 Thread Joseph Miller
Hello- In using gCentroid with a shapefile (SpatialPolygonsDataFrame) output is a SpatialPoints object stripped of all attribute data. Is there a way to preserve some or all of the original data frame attributes when using gCentroid? If not, another solution to creating center points from polygo

Re: [R-sig-Geo] plotRGB - negative values

2013-08-06 Thread Etienne B. Racine
Amom, you could probably offset your values and make sure you set scale to the right value. If your lowest value is 1 and the highest value then becomes 2, then something like plotRGB(all.pca + 1, scale = 2) I don't think you can use negative values in plotRGB. Etienne 2013/8/6 amluiz > Hello

[R-sig-Geo] plotRGB - negative values

2013-08-06 Thread amluiz
Hello, I'm trying to use plotRGB. But I've negative values in my rasters cells and I get the following error: plotRGB(all.pca) Error in rgb(RGB[, 1], RGB[, 2], RGB[, 3], alpha = alpha, max = scale) : color intensity -1, not in 0:255 Is there some way to solve this? Thanks in advance, Amom

[R-sig-Geo] gstat- Interpolate geostatistical data with kriging

2013-08-06 Thread Moshood Agba Bakare
Dear All, I have carried out variogram analysis on geostatistical data using gstat package of R. after rescaling my spatial coordinates, it ranges as follows easting: min=0 and max=794.256 meters northing: min= and max=621.8851 meters the maximum distance between pairs of observation is 1008.754 m

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

2013-08-06 Thread Roger Bivand
On Tue, 6 Aug 2013, Mauricio Zambrano-Bigiarini wrote: On 08/05/2013 04:37 PM, Jon Olav Skoien wrote: 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

Re: [R-sig-Geo] plotKML coerce levels ?

2013-08-06 Thread Tom Philippi
Florence-- I might misunderstand your question, but I think that your problem is the type in your SpatialPixelsDataFrame, not your call to plotKML. From what you describe, I assume that the values in your SpatialPixelsDataFrame are numeric. If you want a specific color legend, you need to use a f

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

2013-08-06 Thread Mauricio Zambrano-Bigiarini
On 08/05/2013 04:37 PM, Jon Olav Skoien wrote: 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) numb

Re: [R-sig-Geo] Changing values of raster pixels to NA

2013-08-06 Thread Zia Uddin Ahmed
You can try this: files <- list.files(pattern='.tif$') s <- stack(files) names(s) nl <- nlayers(s) # Create Raster: values <1: rc1<-function(x){ ifelse(x>1.000,0, ifelse(x<=1.000,1,NA))} r1.rc1<-calc(s,fun=rc1) plot(r1.rc1) r1<-r1.rc1*s plot(r1) # Create raster: values >=-1 and <=1 rc2<-f

Re: [R-sig-Geo] R2Winbugs - problem with workign directory

2013-08-06 Thread Virgilio Gomez Rubio
Hi James, Try to use the full path (e.g., /home/james/mycode/...) instead of just ".". Also, from the manual page: If 'useWINE=TRUE' is used, all paths (such as 'working.directory' and 'model.file', must be given in native (Unix) style, but 'bugs.directory' can be given in Windows pat

[R-sig-Geo] plotKML coerce levels ?

2013-08-06 Thread TROUDE Florence
Dear list, I use the plotKML library to plot maps of simulation results (pollen maps). It's quite great but I'd like to have some colors at some defined levels such as : 2 -> white 4 -> green 10 -> yellow 30 -> orange 60 -> red 100 -> black my object is a SpatialPixelsdataframe or a SpatialGridd

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

2013-08-06 Thread Jesse Berman
Elena, Is this projected data? I previously noted that zerodist() fails to work under certain projections (http://r-sig-geo.2731867.n2.nabble.com/zerodist-fails-with-Projected-Data-td7584211.html). Dr. Pebesma has noted this bug. Jesse - ---

Re: [R-sig-Geo] install rgdal to query postgis database

2013-08-06 Thread Roger Bivand
On Sun, 4 Aug 2013, Duncan McPherson wrote: An update that may shed some light on my issue: Mysteriously it started working, I was able to import some spatial tables form PostGIS to R, and plot them, and I was about to thank the list and close the discussion, when it broke again. This was jus

Re: [R-sig-Geo] Changing values of raster pixels to NA

2013-08-06 Thread Roman Luštrik
You should use `stackApply`. library(raster) fn <- system.file("external/test.grd", package="raster") s <- stack(fn, fn) stackApply(s, indices = 1:nlayers(s), fun = mean) Cheers, Roman On Tue, Aug 6, 2013 at 10:53 AM, Carlos Vázquez wrote: > Hi to all, > I have a problem using raster stacks

[R-sig-Geo] Changing values of raster pixels to NA

2013-08-06 Thread Carlos Vázquez
Hi to all, I have a problem using raster stacks in R. I have a set of raster layers with pixel reliability information for NDVI with pixel values of 0,1,2, and 3 Since i am not interested in values 2 and 3, I want to change all 2 and 3 values in the raster layers to NA's I am using the following