The goal of neogeonames is to provide a useful subset of the GeoNames Gazetteer Data with functions to infer ISO3166 codes for place name queries in a hierarchical manner without a REST API. This package also includes coordinates and shape data for plotting maps, language codes, and timezone data.
You can install the development version from GitHub with:
# install.packages("devtools")
devtools::install_github("dnanto/neogeonames")This is a basic example which shows you how to standardize potentially misspelled place name into a set of ISO3166 codes.
library(neogeonames)
geo <- adminify_delim("USA: Fairfax County, Virginia", delim = "[:,]")
geo
#> $id
#> ac0 ac1 ac2 ac3 ac4
#> 6252001 6254928 4758041 NA NA
#>
#> $ac
#> ac0 ac1 ac2 ac3 ac4
#> "US" "VA" "059" NA NA
paste(Filter(Negate(is.na), geo$ac), collapse = ".")
#> [1] "US.VA.059"Here’s another example with misspelled name…
geo <- adminify_delim("USA: Furfax County, Virginia", delim = "[:,]")
geo
#> $id
#> ac0 ac1 ac2 ac3 ac4
#> 6252001 6254928 4758041 NA NA
#>
#> $ac
#> ac0 ac1 ac2 ac3 ac4
#> "US" "VA" "059" NA NAUse the geonameid to get the coordinates.
# get the id that occurs before the first NA value
idx <- which(is.na(c(geo$id, NA)))[[1]] - 1
with(geoname, geoname[geonameid == geo$id[idx], c("longitude", "latitude")])
#> longitude latitude
#> 403895 -77.27622 38.83469Here’s another example using regular expressions…
adminify_regex(
"USA: Furfax County, Virginia",
list(pattern = "(.+):\\s*(.+)\\s*,\\s*(.+)", names = c("ac0", "ac2", "ac1"))
)
#> $id
#> ac0 ac1 ac2 ac3 ac4
#> 6252001 6254928 4758041 NA NA
#>
#> $ac
#> ac0 ac1 ac2 ac3 ac4
#> "US" "VA" "059" NA NAPlot all feature codes in the US state of Virginia.
library(ggplot2)
df.virginia <- with(geoname, geoname[which(country_code == "US" & admin1_code == "VA"), ])
ggplot(df.virginia, aes(longitude, latitude)) +
geom_point(aes(fill = feature_code), pch = 21, size = 2, alpha = 0.75) +
guides(fill = guide_legend(nrow = 1)) +
coord_map() +
theme_minimal() +
theme(legend.position = "bottom")
Plot all world capitals.
df.capital <- with(geoname, geoname[which(feature_code == "PPLC"), ])
ggplot() +
geom_polygon(
data = shape, color = "black", fill = "white",
aes(long, lat, group = group)
) +
geom_point(
data = df.capital, fill = "blue", pch = 21,
aes(longitude, latitude)
) +
theme_minimal()
Also, check out the vignette("neogeonames").
- Data updated: Mon Sep 14 14:45:03 2020
- Data license: CC BY 4.0