Hello Edzer,
I was curious to run the example you provided (using the latest CRAN
versions of all packages), but as written the line
> idw.spdf = as(idw.out, "SpatialPolygonsDataFrame")
produces
Error in as(idw.out, "SpatialPolygonsDataFrame") : no method or default for
coercing "SpatialPixels
Ok, I'm nearly there using only RGDAL and R-base commands, if I can get
a bit of feedback I might have a decent tiled (row-by-row) processing
structure worked up. Couple of things -- first, I was mistakenly
thinking all raster formats can even really support line-by-line
writing, which is not
On Thu, 30 Oct 2008, Jean-Paul Kibambe Lubamba wrote:
Hi,
I wonder if anyone know whether there is a R package for implementing a
geographically weighted poisson regression model.
For what it's worth, see the draft ggwr() function in the spgwr package.
Roger
Thanks in advance.
Jean-Pau
Hi,
I wonder if anyone know whether there is a R package for implementing a
geographically weighted poisson regression model.
Thanks in advance.
Jean-Paul
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r
Hi,
maybe the RSAGA package has the solution to your problem; there are
actually two ways of applying functions to grids in RSAGA, (1) by
row-by-row processing (special form of tiles), or (2) using the SAGA
binaries.
1) local.function
local.function and focal.function are very flexible tool
Hi Jonathan,
I use GDAL to read subsets of raster and image files, though usually I employ
the python bindings. I don't actually have answers, but maybe ideas about how
to get there...
> 1) The code doesn't actually seem to work -- I get a tiny, unreadable
> TIFF out of the back of this algorit
Dear all,
In continuation of this thread, I've spent some time looking at kriging
on the sphere, corrected some bugs, and need further help.
First of all, distances for covariances on the sphere were computed
incorrectly as well in gstat, so I hope not too many people have been
relying on th
Hi All,
I need a function which checks whether SpatialPolygons object includes
SpatialLines object.
If the polygon is convex, it is possible to use point.in.polygon(),
however, in non-convex case, it is not sufficient.
Does anyone have a good idea?
thanks,
Satoshi Koide
Alessandro,
The rsaga.geoprocessor needs to get all SAGA parameters via the argument
"param", e.g.:
> rsaga.geoprocessor(lib="geostatistics_grid", module=6,
> param=list(ZONES="plots.sgrd",
CATLIST=0, STATLIST="DCM_1_power10.sgrd", OUTTAB="stat_DCM_1_power10"))
You seem to often forget that s