Re: [R-sig-Geo] projection of ggmap:get_map() output

2014-04-25 Thread Agustin Lobo
I've tested the effect of the projections on overlaying graphs using levelplot() + layer() In short, the user's object must be in epsg:4326, probably because layer() carries out a reprojection. The bad news is that it carries out a wrong reprojection, perhaps because the incorrect (or insufficient)

Re: [R-sig-Geo] projection of ggmap:get_map() output

2014-04-25 Thread Agustin Lobo
According to the results https://dl.dropboxusercontent.com/u/3180464/ggmap_GM_compare.odp of the code below, the matrix downloaded by get_map() is on epsg:3857, the so-called Google Pseudo Mercator, while the bounding box is given in epsg:4326 coordinates. The odp file shows the raster as created b

[R-sig-Geo] kerneloverlaphr and BRB

2014-04-25 Thread Anna-Marie Corman
Dear all, I apologize in advance for any stupid question. But I am completeIy new in this kind of analyses and confused with some of the necessary input parameters... I want to compare my BRB outputs of different individuals (gulls) from DIFFERENT COLONIES (i.e. single breeding colonies are ab

Re: [R-sig-Geo] projection of ggmap:get_map() output

2014-04-25 Thread Oscar PerpiƱan
> On 04/24/2014 03:52 PM, Barry Rowlingson wrote: >> So do you mean (and I've not tried your code yet, sorry) that the >> corners may well be at the given lat-long points but half way along any >> edge might not be at the halfway-point in lat-long? > Indeed: > > library(ggmap) > ggmap(get_map(matri

Re: [R-sig-Geo] trouble with makegrid

2014-04-25 Thread Pascal Oettli
Hello Erin, Your example is not reproducible. Regards, Pascal On Fri, Apr 25, 2014 at 12:28 PM, Hodgess, Erin wrote: > Hi everyone! > > I have been using makegrid and it worked fine. Until tonight. I upgraded to > R-3.1.0 on my Ubuntu, and am using plotKML and gstat from R-Forge. > > Here is

[R-sig-Geo] Legend of map

2014-04-25 Thread Komine
Dear all, I used this code to plot a map. >library(raster) >Biomasse<-rasterFromXYZ(tab) >PlotBio<-plot(Biomass) >PlotBio<-plot(Biomass,col=(topo.colors(100)),legend=TRUE,cex.axis=0.8,adj=1) >plot(Biomass,legend.only=TRUE,col=(topo.colors(100)),legend.args=list(text="Biomass(Kg/ha)",side=3,