[R-sig-Geo] rgdal spatstat: gtiff

2012-11-06 Thread Raffaele Morelli
Hi, I am using this code snippet to produce Gtiff images from density.ppp z-density(mypattern, 6000,dimyx=200) zgrid - as.SpatialGridDataFrame.im(z) writeGDAL(zgrid, 01.tif, type=Byte, options=INTERLEAVE=PIXEL) Gtiff image is generated but gdalinfo among others reports that: Image

[R-sig-Geo] Call for Abstracts: VII Jornadas de SIG Libre, March 2013

2012-11-06 Thread Jornadas de SIG Libre
Dear, Call for Abstracts to VII Jornadas de SIG Libre (Spanish FOSS4G) is still open until November 21st. We count on your active participation! The Spanish FOSS4G, aka Jornadas de SIG Libre, is the Spanish event centered on free and open source software for Geospatial, that takes place in

[R-sig-Geo] error in Crop command from Raster package

2012-11-06 Thread kapo coulibaly
I'm running the following code to cycle through individual polygons (surface water basins) in a shapefile and averaging data from a netcdf: # Extracting basin ID values basinid-basin_latlong$value for (i in basinid) { #individual basin extraction

Re: [R-sig-Geo] raster package - zonal - standard deviations

2012-11-06 Thread Christian Levers
Thanks a lot Robert! Indeed it's (obviously) a matter of memory (although our server has approx 200GB RAM). While simulating RAM limitations, I obtain the same message: error in .local(x, z, ...) : stat can be 'sum', 'mean', 'min', or 'max'. Changing maxmemory doesn't improve things, even

Re: [R-sig-Geo] error in Crop command from Raster package

2012-11-06 Thread Robert J. Hijmans
Kapo, I am guessing that this is happening because in some cases the extent used to crop only represents a single cell (or part of a single row or column), and that ncdf returns a vector, not a matrix, but I need to look into that. It would be helpful if you could give me access to your data, or

[R-sig-Geo] how to draw a new circle/round polygon on existing plot

2012-11-06 Thread Wouter Steenbeek
I have a plot with several layers of SpatialPoints and Polygons. They are all projected and in the correct CRS, so my map looks like it should. I now want to draw a circle with a radius x on a part of this map. This should not just be a circle, but be in the correct projection (so that I can

Re: [R-sig-Geo] error in Crop command from Raster package

2012-11-06 Thread kapo coulibaly
You might be right I did just that before reading your reply (I added a percentage of xmax-xmin and ymax-ymin to the boundaries of the extent) and it worked. I'm using a loop because my netcdf has 4000 time steps so if I refine the whole thing (trmm) by 10 I run out of memory. But by doing one

Re: [R-sig-Geo] raster parallelization

2012-11-06 Thread Etienne B. Racine
When running Robert's example I get NA stripes depending on the number of cores I set. plot(stack(rf, rf2)) all.equal(rf[], rf2[]) I was wondering if it was only on my side. Looking at the source I don't see how clusterR is accounting for the buffer around the tiles. So it would be normal to have

Re: [R-sig-Geo] how to draw a new circle/round polygon on existing plot

2012-11-06 Thread Edzer Pebesma
x = 10 # center x y = 10 # center y n = 100 # nr of pts r = 5 # radius pts = seq(0, 2 * pi, length.out = n) plot(sin(pts), cos(pts), type = 'l', asp = 1) # test require(sp) xy = cbind(x + r * sin(pts), y + r * cos(pts)) sl= SpatialLines(list(Lines(list(Line(xy)), line))) plot(sl, add=FALSE,

Re: [R-sig-Geo] raster parallelization

2012-11-06 Thread Robert J. Hijmans
Etienne, Thanks, I am not sure what I was thinking. focal is an example of functions that do not work well with clusterR. Robert On Tue, Nov 6, 2012 at 2:06 PM, Etienne B. Racine etienn...@gmail.comwrote: When running Robert's example I get NA stripes depending on the number of cores I set.

[R-sig-Geo] [GSIF] what API to call for a simple reboxing test?

2012-11-06 Thread Tom Roche
summary: I'm trying to test reboxing, or 3D interpolation, using package=GSIF. The idea is to * assume a uniform input field * create a uniform input grid with box/voxel dimensions=2x2x2 with each box having value=1 * rebox to a uniform output grid with dimensions=1x1x8: if output box

Re: [R-sig-Geo] regridding netCDF?

2012-11-06 Thread Tom Roche
Tom Roche Thu, Nov 1, 2012 at 2:21 AM See code [for lon-lat to LCC regridding] @ https://github.com/TomRoche/GEIA_to_NetCDF esp https://github.com/TomRoche/GEIA_to_netCDF/blob/master/regrid.global.to.AQMEII.r ping yang Tue, 6 Nov 2012 20:36:12 -0500 (rearranged) Thank you so much for your

Re: [R-sig-Geo] [gstat] simple reboxing test fails

2012-11-06 Thread Edzer Pebesma
On 11/07/2012 04:39 AM, Tom Roche wrote: ... But all my predicted box values=NA; how to fix? details: ... # Here I fail to get the values for the output grid :-( data3D.out - krige(formula=v ~ 1, data3D.in, grid3D.out) summary(data3D.out) # gets # Object of class