[R-sig-Geo] Legend unit for raster layer and create a color map exported to tiff

Hallo, i have a problem to understand how to set unit value, for the parameter, of a raster layer to the bar legend. The same problem occur for SpatialGridDataFrame. Another problem is export a tiff image with writeRaster with a colours "map". Thanks in advance [[alternative HTML version

Re: [R-sig-Geo] gwr.morantest in spgwr package.

On Wed, 25 Jul 2012, Lilis Husna wrote: I want to ask about gwr.morantest in spgwr package. What is the purpose of doing this test? Which variabel is tested? residual from global model (OLS), residual from local model (GWR), or the dependent variabel (same as testing spatial autocorrelatio

Re: [R-sig-Geo] Donut holes

rgeos::gDifference should do what you are after. HTH, Simon Simon O'Hanlon, BSc MSc Department of Infectious Disease Epidemiology Imperial College London St. Mary's Hospital London W2 1PG From: r-sig-geo-boun...@r-project.

[R-sig-Geo] Donut holes

Can somebody suggest a method for subtracting one SpatialPolygon* object from another, so that I'm left with a polygon with a donut hole? I can't seem to find one. Thanks. -- View this message in context: http://r-sig-geo.2731867.n2.nabble.com/Donut-holes-tp7580607.html Sent from the R-sig-geo

Re: [R-sig-Geo] How can I plot a shapefile with geom_map from ggplot2

Hi Barry, Just a small addition based on your example. ## fortify(columbus,region="id") ## did not work for me # but this works: library(spdep) example(columbus) columbus@data$AREA fco <- fortify(columbus, region= "AREA") ggplot(columbus@data) + geom_map(map=fco, aes(map_id=AREA,fill=AREA)) + ex

Re: [R-sig-Geo] Google Mercator EPSG:900913, in which init file to add it?

On Wed, Jul 25, 2012 at 4:21 PM, Mathieu Rajerison wrote: > Hi, > > I'm using Google Mercator (EPSG:900913) for one purpose but it doesn't seem > recognized by rgdal.. > > Where to add it please? 900913 was a 'cute' number used unofficially because it looks like 'google'. Eventually after a coup

[R-sig-Geo] Google Mercator EPSG:900913, in which init file to add it?

Hi, I'm using Google Mercator (EPSG:900913) for one purpose but it doesn't seem recognized by rgdal.. Where to add it please? Thanks! Mathieu [[alternative HTML version deleted]] ___ R-sig-Geo mailing list R-sig-Geo@r-project.org https://sta

Re: [R-sig-Geo] how to hange the size of raster maps

Dear Agus, I'd suggest using the pdf() or png() functions in your script, rather than relying on RStudio's exporting. This way you can specify the size. best, Adam On Wed, Jul 25, 2012 at 8:36 AM, Agus Camacho wrote: > Dear all, > > I am trying to change the size to raster maps generated in R

[R-sig-Geo] how to hange the size of raster maps

Dear all, I am trying to change the size to raster maps generated in R. However, the differently layers separate when i do this. Any body knows if its possible and how to avoid this? Thanks There it goes an reproducible example alt <- getData('alt', country='CHE') slope <- terrain(alt, opt='slope

Re: [R-sig-Geo] problem extracting Z value from line shapefile

On Wed, Jul 25, 2012 at 8:22 AM, zbynek.jano...@centrum.cz wrote: > Warning message: > In readOGR("C:/Users/Janoska/Documents/GPS", "tesikov") : > Z-dimension discarded > > Clearly, Z dimension is stored in the data, but readOGR does not read it. > Hence, I followed this post: > https://stat.et

[R-sig-Geo] problem extracting Z value from line shapefile

Hello, I have following problem: I wish to extract z coordinates from a polyline shapefile. Using readOGR() i get following: > data <- readOGR(path,'tesikov') OGR data source with driver: ESRI Shapefile Source: "C:/Users/Janoska/Documents/GPS", layer: "tesikov" with 1 features and 1 fields Featu