similar to: Change panel background color in spplot()

Hi! Is there a way to manually costumize the color legend in an spplot() - especially where to draw ticks and labels for the ticks? The reason I'm asking: Usually spplot() automatically divides the data into fitting slices and makes a color legend (also automatically). I want to assign the slices myself and have a fixed scale instead of an automatic/dynamic scale. I think what I want

Hi, I'm using R 2.5.1 and gstat 0.9-39 library and I'm working with the jura data set provided by gstat library. I tried to plot a graph of metal concentrations (let's say Cd) with the command spplot, but I realized that the default lags are equally distributed between Min and Max. I did: library(gstat) data(jura) jura.pred.xy <- jura.pred names(jura.pred.xy)[1:2] <-

Hi! Is there a way to manually costumize the color legend in an spplot() - especially where to draw ticks and labels for the ticks? The reason I'm asking: Usually spplot() automatically divides the data into fitting slices and makes a color legend (also automatically). I want to assign the slices myself and have a fixed scale instead of an automatic/dynamic scale. I think what I want gets

hi all, I have a plot to which i would like to add text labels. And i cant find a way...here is the code : enaD2<-idw(D2~1, loca=dva, newdata=grd) pts = list("sp.points", dva, pch = 20, cex=1.5, col = "darkred spplot(enaD2, "var1.pred",sp.layout=pts, main = "globina 60 cm", sub="D2",

Hi, Is there a simple way to suppress the plotting of polygon boundaries with spplot() ? # simple list of 12 colors cols <- brewer.pal(12, "Paired") # plot pile of polygons, with 12 classes: spplot(x, zcol='class2', col.regions=cols, scales=list(draw=T), xlab="Easting (m)", ylab="Northing (m)") ... seems to work well. However the polygon boundaries

Hello: I would like to superimpose vector polygons (state outlines) from a Shape file on top of a satellite image, imported into a SpatialGridDataFrame from GEOTIFF via gdal_translate and readGDAL. When I plot polygon and point shape files in R, into SpatialPointDataFrame and SpatialPolygonDataFrame, the two feature sets line up geographically, so it seems logical that a SpatialGridDataFrame

Hello! I have 2 maps - both created in ssplot and both identical in terms of outline. Is there any way to superimpose Map1 (which has black borders between Canadian provinces) onto Map2 (which is also a map of Canada)? Thanks a lot for your hints! Dimitri ### A. Reading in Canada data at the province and then at the county level: library(raster) getData('ISO3') # Canada's code is

Hello, I want to create side-by-side maps of similar attribute data in two different cities using a single legend. To simply display side-by-side census block group boundary (non-thematic) maps for Minneapolis & Cleveland I do the following: library(rgdal) library(sp) Minneapolis=readOGR("../Minneapolis/Census/2010/Census_BlockGroup_GEO/","tl_2010_27053_bg10")

Hola, Lo hice de la forma en que propones y no me sale, no sé si algo esté mal con R, porque no me respeta ni el ?key.space?, me lo sigue poniendo en la parte izquierda del gráfico Saludos El 06/05/2014, a las 02:21, Olivier Nuñez <onunez en unex.es> escribió: > El segundo argumento en spplot ha de definir el número/nombre del intervalo al que > pertenece el area a pintar. >

I have some geospatial data where two layers are thematic and the third is a percentage, so the maps need to have different themes. thumbDf <- as( stack( thumb), "SpatialGridDataFrame") names(thumbDf at data) <- c("pri", "sec", "pct") thumbDf at data$pri <- factor(thumbDf at data$pri, levels=c(0:8), labels= names( peelClasses)) thumbDf at

Hola, Intento representar en un mapa participaciones porcentuales de los sectores económicos y no logro hacerlo con spplot. He intentado con spplot(zm, c(?part88?, ?part93?), cuts=4, col.regions=brewer.pal(4, ?Set3?)) Lo pude hacer utilizando el base graphics de R, definiendo: > plotvar88 <- zm$part88 > nclr <- 8 > plotclr <- brewer.pal(nclr, "PuOr") > plotclr

Hola, El problema con la propuesta de Olivier es que los intervalos son diferentes para cada variable. La forma sencilla y rápida es: spplot(zm["part88"], col.regions=plotclr, at=class$brks) Pero para que quede más elegante hay que dar algunos pasos más: ## Intervalos en forma character op <- options(digits=4) tab <- print(class) options(op) intChar <- names(tab) ## Indice

Buenas tardes, Estoy probando a leer archivos netcdf (como raster stack) y crear animaciones de datos de modelos climáticos y estaba utilizando para generar las imágenes spplot(). Estás imágenes tienes superpuestas varias capas de imagenes, el mdt de base, nubosidad, precipitación y nieve. Mi pregunta es la siguiente,¿es posible visualizar la leyenda de cada una de las capas? Actualmente solo

Hola, no me funciona: Te contesto lo mismo que le contesté a otro compañero de la lista, espero ser más claro. Les comparto la carpeta con la capa y los datos que estoy usando. https://www.dropbox.com/s/umy0evt3qm1wr4d/dissolve.zip Esto es lo que estoy haciendo: library(rgdal) library(maptools) library(sp) library(classInt) library(RColorBrewer) zm <- readOGR(?.?, ?zmdis?) data <-

Hello- After looking through ?spplot, I would expect that I could specify the values of the cuts: "...?cuts? number of cuts or the actual cuts to use..." So in the following command, spplot(lzm.krige.dir["var1.pred"], scales=list(draw=TRUE), xlab="Easting",ylab="Northing",

Hi, I am very new to R and I am trying to write KML files. I am running R-2.13.1 and I have the following packages: - AKIMA - RGDAL - RSAGA - SP I was trying to run the simple example that I found at the following web page: http://spatial-analyst.net/wiki/index.php?title=Export_maps_to_GE The lines that I can run are: data(meuse.grid)

Hello All, I have a discontinuous dataset and I used spplot to plot the points. I was wondering whether it is possible to show the legend of the plot as color ramp. I saw that on continuous data. If anyone could tell me that would be great. Another thing I would like to know is how can we create the horizontal legend. By default, the legend is vertical. I was able to shift the legend inside

Soon on CRAN a new version of package gstat will emerge, which has a few minor changes and possible incompatibilities w.r.t. the previous version(s). The new gstat (0.9-23) now: + depends on sp, and uses internally with Spatial* classes from sp if data are provided in the old-fashoned way (as data.frame) + has a vignette to get you started with the classes in sp + defines krige as a generic;

Soon on CRAN a new version of package gstat will emerge, which has a few minor changes and possible incompatibilities w.r.t. the previous version(s). The new gstat (0.9-23) now: + depends on sp, and uses internally with Spatial* classes from sp if data are provided in the old-fashoned way (as data.frame) + has a vignette to get you started with the classes in sp + defines krige as a generic;

Dear Helper, I have a spatial lines data frame object 'spRiverDf'. The data frame consists of numbers {0,1,...,5}. And I have a vector 'colorS' of length 6 with different colours. If I make a plot with spplot I get a plot of the lines - colours depending on there number in the data frame column: spplot(spRiverDf['data.col.1'], zcol=..., names.attr=...,