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
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
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
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
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
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
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
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
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,
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.
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
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
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
13 matches
Mail list logo