search for: gadm

Dear All, Please consider the snippet at the end of the email. I often download some maps (in the R format) from http://www.gadm.org/ However, when I run (typically more than once) a variation of the script below (based on http://bit.ly/1b3W0Aa ), I often get Error in load(url(paste("http://gadm.org/data/rda/", fileName, "_adm", : cannot open the connection In addition: Warning message: In load(u...

...lor to every level) (2) Be able to control the position of the legend and the size of the character used in the legend itself. Any suggestion is welcome. Cheers Lorenzo ################################################? ## you will need the sp-package library('sp') ## load a file from GADM (you just have to specify the countries "special part" of the file name, like "ARG" for Argentina. Optionally you can specify which level you want to have loadGADM <- function (fileName, level = 0, ...) { load(url(paste("http://gadm.org/data/rda/", fileName...

...at the end of the email), well, nothing happens and I certainly do not see a map of Switzerland. I am not familiar with the sp package, but am I missing the obvious? Any help is appreciated. Best Regards Lorenzo ########################################? library(sp) con <- url("http://gadm.org/data/rda/CHE_adm1.RData") print(load(con)) close(con) language <- c("german", "german", "german","german", "german","german","french", "french", "german","german","french&quot...

Dear all, in GADM map there are three levels (nation, province and precinct) for each country of the world but for all of them you are never able to plot only one part of a chosen country. To be sure, I am trying to plot only one region of “Italy” and colour the different precincts in it. So far I am able to colour...

...brary(RColorBrewer) library(sp) require(foreign) require(latticeExtra) reguire(lattice) #Example data lat <- seq(-12, -1, 0.1) lon <- seq(28, 45, 0.1) data <- expand.grid(lat=lat, lon=lon) n <- nrow(data) data$X <- runif(n, 0, 250) #Read map boundaries TZA_adm2 <- url("http://gadm.org/data/rda/TZA_adm2.RData") print(load(TZA_adm2)) close(TZA_adm2) col = rainbow(length(levels(gadm$NAME_2))) spplot(gadm, "NAME_2", col.regions=col, auto.key = FALSE, main="TZA Regions") #How do I clipp or overaly X in data on TZA_adm2 polygon? #Also how to remove the leg...

I am having trouble getting data labels to display over the provinces in a GADM map of Canada. Specifically, I need the variable "Number" from the data set "by_province", grouped by "region", to appear on the corresponding regions of the map. The data set "by_province" looks like this: long lat order hole piece...

...ould like to place region names (or abbreviations if space is too limited) on a map. I have tried the following, copying and adapting various examples I have seen but without success. The first example is self-contained, but the second one requires the .shp file, which can retrieved at http://www.gadm.org/country. ----- Example 1 ----- library(sp) con <- url("http://gadm.org/data/rda/ITA_adm1.RData") print(load(con)) close(con) regionnames <- c("Piemonte", "Valle d'Aosta", "Liguria", "Lombardia", "Trentino-Alto Adige/S?d...

Dear All, I like very much figure 2.a) and 2.b) of this paper http://www.nature.com/srep/2013/130410/srep01640/pdf/srep01640.pdf and I probably need a similar visualization. Is anything like that doable in R? I have some experience with R and gadm (gadm.org), but that is only to produce maps colored according to a scalar. I do not know how to visualize maps with a "perspective", let alone adding the histograms in 3D. Sorry for not providing an example script, but if I had one I would not be posting at all. Any suggestion is appreci...

...999 Grenoble 99 ... I have looked at different graphing packages such as 'sp', 'raster' and 'ggplot2', but am not sure which one is most appropriate for my purposes. Furthermore, I don't understand how you retrieve city coordinates from shp, gadm, etc., maps generated, or downloaded from the web and then processed, with those packages (SHP, GADM, etc.), and plot data at those coordinates. Any help would be much appreciated. Many thanks. Philippe

...lack 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 "CAN" can1<-getData('GADM', country="CAN", level=1) can2<-getData('GADM', country="CAN", level=2) ### B. Creating Map1- with borders between provinces: can1 at data[["groups"]]<-1 # want all provinces to have the same color, but ideally I'd like the whole map to be tran...

Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20120614/518ed68f/attachment.html>

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

...contains 40 polygons. Now I want to join some of the polygones. I cannot figure out how to do that? Any ideas? The code below produces a map of Germany and I need to join some counties. library(sp) library(RColorBrewer) # get spatial data for Germany on county level con <- url("http://gadm.org/data/rda/DEU_adm3.RData") print(load(con)) close(con) # plot Germany with random colors col = rainbow(length(levels(gadm$NAME_3))) spplot(gadm, "NAME_3", col.regions=col, main="German Regions", colorkey = FALSE, lwd=.4, col="white") ????????????????...

Dear R'ers, I am wondering what is the smallest geographicterritorial unit available for formatting in Canada. Provinces? I know that in the US it is the county so that I can color US counties any way I want, for example: ### Example for coloring US counties ### Creating an ARTIFICIAL criterion for coloring US counties: library(maps) allcounties<-data.frame(county=map('county',

...car Perpiñán Lamigueiro Dpto. Ing. Eléctrica, Electrónica, Automática y Física Aplicada (ETSIDI-UPM) Grupo de Sistemas Fotovoltaicos (IES-UPM) URL: http://oscarperpinan.github.io 2015-04-02 23:33 GMT+02:00 Ivan Aguilar <iaguilar en gmail.com>: > 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]]@Poly...

...co la información es muy clara ni abarcativa. ¿Alguno/a trabaja con mapas en R y me puede indicar por donde encarar? En principio mi interés es realizar mapas sencillos con información por provincias y departamentos de Argentina. (En cuanto a los shapefiles de Argentina tomé uno de aquí http://www.gadm.org/ que no incluye las Islas Malvinas y el correcto del Instituto Geográfico Nacional pero que al incluir el territorio antártico lo hace menos practico). Se agradece cualquier apunte o pista de por donde empezar a paso seguro. Saludos y buen año! Santiago [[alternative HTML version deleted]]

Hola a todos. Soy nuevo en R, estoy aprendiendo con tutoriales y lo hago sin un fin comercial y/o laboral. Necesito ayuda para mostrar un mapa y sobre ese mapa poner un punto dado por una coordenada. con esto consigo mostrar el mapa. library(raster) arg<-getData('GADM', country='ARG', level=0) plot(arg) En ese mapa quiero mostrar un punto que está en el sistema de coordenadas de google en -32.571723, -58.991303 gracias! -- Rafael Bidegain Ya que los cuerdos no hablan, hablará el loco. (The Fool, Padraic Pearse) # yo quiero educación libre. ¡YA!...

...f class "owin" (observation window) in the package spatstat from GIS mapping data. Here's an example of my problem. Everything goes well until the last line of code. I get the error message shown at the bottom: library(spatstat) library(sp) library(maptools) con <- url("http://gadm.org/data/rda/CHE_adm1.RData") class(con) print(load(con)) close(con) SP <- as(gadm, "SpatialPolygons") W <- as(SP, "owin") "Error in owin(poly = opls) : Polygon data contain overlaps between polygons" Thanks for any help, Axel. [[alternative HTML versi...

Hi R users I need the map of France?s ? communes ? (towns) to build a map Is there a way to get it? More generally: How to do to get the map of a country and its different geographical levels? Best regards -- View this message in context: http://old.nabble.com/map-of-a-country-and-its-different-geographical-levels-tp26225645p26225645.html Sent from the R help mailing list archive at

Hello, I am wondering about the X11() graphic device on Windows. I try to plot a raster image but nothing gets displayed. I found some pages where it is mentioned that x11() not always supports raster rendering. Is there any add on for x11, any update or any R-package which solves that displaying problem in Windows? What I try to test it is an example from the package {raster}: library(raster)