I imported a Shapefile holding a quite detailed world map from: http://www.diva-gis.org
library(sp)
library(rgdal)
world = readOGR('./countries_shp/', 'countries')
Using R with the sp-Packages, plotting is really easy:
spplot(world, 'UNREG2')
plot(world)
Now I want to
- draw a customizable grid
- add points derived from coordinates
For example importing coordinates of Vienna (48°12′32″N 16°22′21″E) can be done using char2dms of the sp-package:
vienna.N = char2dms("48d12'32\"N")
vienna.E = char2dms("16d22'21\"E")
But how can a draw a point at this exact location on my map?
Also adding a grid is straight forward:
plot(gridlines(world), add=TRUE)
assuming, that the last plot was created with "plot(world)".
How can this be achieved with the spplot output?
Best Answer
I figured out how to print a point for Vienna on the map.
Importing the coordinates from WikiPedia into R:
Building a sp-object
The projection is taken from http://toolserver.org/~geohack/geohack.php?pagename=Vienna¶ms=48_12_32_N_16_22_21_E_type:city_region:AT where the WikiPedia article links to:
Printing this object gives the numeric representation which equals exactly the data from mentioned website:
The resulting
can be used with base- and sp-graphics. First we transform the projection to the one used in the Diva worldmap:
Then a plot for Europe can be produced: