I'm not too familiar with GIS terminology, but would you agree that
sp::over does (or allows for) a spatial table join?
On 07/06/2011 06:10 PM, Agustin Lobo wrote:
> Yes!
>
> I'm doing:
>> delme2 <- over(q1km, rndp,fn = mean)
>
> And now not only the row.names are correct and the order of the 2
Yes!
I'm doing:
> delme2 <- over(q1km, rndp,fn = mean)
And now not only the row.names are correct and the order of the 2
arguments makes more sense, but the dimensions of the result are also
the same
that that of the first argument
(which implies rows of NA for those polygons with no points: the
Agus, sp::overlay will be deprecated at some stage in favour of
sp::over, which is a more consistent and complete approach to the same
problem. For instance, overlay(x,y) would do the same as overlay(y,x),
which is not good, if you think about it a bit longer (my own mistake,
long time ago).
Could
I'm using sp::overlay between a SpPolDF and a SpPointsDF:
> delme <- overlay(rndp,q1km, fn = mean)
> class(q1km)
[1] "SpatialPolygonsDataFrame"
attr(,"package")
[1] "sp"
> class(rndp)
[1] "SpatialPointsDataFrame"
attr(,"package")
[1] "sp"
> class(delme)
[1] "data.frame"
but I'm confused by the fa
Horacio
See this thread from December:
http://www.mail-archive.com/r-sig-geo@r-project.org/msg00154.html
I haven't made much progress or further tested this approach since then. It
seems to work for me and my needs.
The example with the North Carolina map relies on pretty() to determine where
Hi,
Just to add one more solution. It uses the rasterVis package (at
R-Forge now and should be at CRAN today).
library(raster)
library(rasterVis)
levelplot(s, par.settings=RdBuTheme)
More examples here: http://rastervis.r-forge.r-project.org/
Best,
Oscar.
---
Oscar Perpiñán
I summarize here the different solutions I got to my question regarding
how to plot several raster layers on the same page and with the same
color scale.
Sorry it took a long time, I was interrupted by something else for weeks.
Hope this is useful.
Many thanks to all.
require(raster)
require(ggpl
Dear Robert,
Thank you for sorting this out. And thanks to all the others for their
valuable comments.
Christian
--
View this message in context:
http://r-sig-geo.2731867.n2.nabble.com/Subsetting-RasterBrick-very-slow-tp6537890p6554009.html
Sent from the R-sig-geo mailing list archive at Nabble
Solved: it was caused by having upgraded to raster 1.8.40 but using
raster objects created under 1.8.38
Agus
2011/7/6 Agustin Lobo :
> A function of mine outputs several raster objects in a list "l":
>> class(l$lisaobsv)
> [1] "RasterLayer"
> attr(,"package")
> [1] "raster"
>
> But by some reason
A function of mine outputs several raster objects in a list "l":
> class(l$lisaobsv)
[1] "RasterLayer"
attr(,"package")
[1] "raster"
But by some reason I get this weird error:
> a = l$lisaobsv
Error in slot(x, slotname) :
no slot of name "z" for this object of class "RasterLayer"
Calls: -> lapp
So, it seems you have one raster for your dependent variable: species
richness and one for each explanatory variable.
I have never accomplished this but I think you should make a raster stack of
all these variables using raster package: the dependent one + the
explanatory ones then fit a model, fo
11 matches
Mail list logo
