similar to: lattice fun: multiple themes in panels with spplot()

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. >

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

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! How does one change the background color of the map-panel in spplot()? Example: library(sp) data(meuse.grid) gridded(meuse.grid) = ~x+y spplot(meuse.grid, "part.a") How would I get another background-color for the map-panel (but not for the whole plot) here? Thank you! Marcel

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

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",

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

Hola He hecho el típico mapa con gadm y spplot: library(sp) ES <- local(get(load(url("http://biogeo.ucdavis.edu/data/gadm2/R/ESP_adm2.RData ")))) acerco las canarias: A <- which(ES en data$NAME_1 == "Islas Canarias") L <- length(ES en polygons[[A]]@Polygons) L <- length(ES en polygons[[A]]@Polygons) for (i in 1:L){ ES en polygons[[A]]@Polygons[[i]]@coords

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

Hola, Yo lo resuelvo usando grid.rect. Tienes una posible solución aquí: https://github.com/oscarperpinan/spacetime-vis/blob/master/choropleth.R#L216 (que es una versión actualizada de este artículo https://procomun.wordpress.com/2012/02/18/maps_with_r_1/) Saludos. Oscar. ----------------------------------------------------------------- Oscar Perpiñán Lamigueiro Dpto. Ing. Eléctrica,

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, I would like to do coloring of map regions based on the region values "weight". The approach I am taking is first to break regions into equal intervals, classIntervals(spdf$weight,4)$brks #4 intervals in this case and coloring all regions within the interval with the same color col = brewer.pal(4,"RdYlGn")) The max "weight" is as well the boundary of the

http://www.nabble.com/file/p25336596/Conductivity1.jpeg I need a little help making modifications to the image included with this post. First, rather than using a linear color legend to display the output I would like to use a log-scale legend. Thus, the legend on the right would go from 1 to 1000, for example, with a classic log-scale gradation. What I hope to avoid is taking the log of the

Hi everyone. I posted this on R-sig-geo but got no response. Can select.spatial() be used in an existing spplot? I have tried selecting points (eq) from a plot generated from sp. However, when I invoke select.spatial(eq). It generates only the points without the background containing the faults. I need the background to select which earthquakes coalesce on which fault. Is there an alternative?

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=...,

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