similar to: GADM Data Download

Dear All, Please see the short script at the end of the email, which I assembled looking for bits and pieces on the web. It essentially does what I need: it plots several countries as a color-coded map. I just would like to fine-tune a bit the final image, in particular (1) Select my own colors for "painting" the countries (i.e. associate manually a color to every level) (2) Be

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 only the part of my interest but the programme still plot the whole country. Is

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 region

Hi Team, Need your guidance in building SPDF objects from Latitude, Longitude Information available. I am using package "plotGoogleMaps" which is really awesome but it requires SPDF objects to build the maps. I have a data frame which have Latitude and Longitude information and wanted to convert it in SPDF for making maps. Please help me to find suitable package or function to convert

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

Hi, The following script (which I did not develop) is used to calculate and plot a skewed normal curve. The script currently requires the user to input six parameters, rather than reading these directly from a file. I've been spinning wheels here, trying to figure out how to modify the script to automate it. I have four data sets, each in excess of 300 records that I need to process. My

Hello, I would 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 <-

I have csv data that extend beyond the area I want for an existing map. I want using the boundaries of the polygon shape file as a cookie cutter so that I can overlay the csv data on map without including anything outside the map boundaries and create a dbf file or shapefile of the clipped data .  The reproducible example: ############################################### library(RColorBrewer)

Dear All, I need to do something relatively simple: generate a map of Europe and paint the various states with different colors (only 3-4 are needed) according to a rule. I would like to keep it as simple as possible and the script resorting to the sp package that I found at http://bit.ly/Oc71ub is exactly what I am looking for (I need to repeat the exercise with Europe instead of

I have a SpatialPolygonsDataFrame object (sp package). It 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")

Dear List, I am analyzing the home range area of fish and seem to have lost the individuals ID names during my manipulations, and can't find out how to rename them. I calculated the MCP of the fish using mcp() in Adehabitat. MCP's were converted to spatial points data frame and exported to qGIS for manipulations. At this point the ID names were lost. I brought the manipulated

Hello, I would like to generate a small map (say 10cm x 10cm) of France showing cumulative numbers by Distribution Center stored in a database: DISTRIBUTION_CENTER COUNT Paris 122 Paris 3 Paris 21 Lyon 12 Lyon 4444 Lyon

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

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

Dear R list, I'm trying to create an object of 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)

Hi all, ?I'm trying to iterate in a SpatialPolygonDataFrame thatcontains 110 features. I tried to use the following code? iterate.spdf =function(x){ for (i in 1:nrow(x)){ p = x[i, ]} return(p)} but with no success as it returns with aSpatialPolygonDataFrame that contains onlyone feature. The code seems to overwrite the polygons to end up with only onefeature. ?Any suggestions?Ashraf, cheers

Salutations, Long story short: I need to create an SPDF to superimpose on a Google map. I have .shp file which has a long list of locations listed in (what I guess must be) UTM. I need to transform these into lon, lat. My use of spTransform, however, keeps giving me the wrong results. For what it's worth, it only seems to be messing up the latitude. My final line of code looks like this:

Hola, # Hacemos el KED. Ver función "krige()": KED.rad <- krige( formula=pluvPcp~layer, # covariable -> radar locations=lluvia.rad.pluv.spdf, newdata=radarGrid, # podría ser cualquier objeto Spatial model=v.fit, # modelo de semivariograma. maxdist=Inf

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

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,