Kort-baserede plot med ggplot2 i R

I dette indlæg vil jeg demonstrere, hvordan rumlige placeringer kan plottes på et kort ved hjælp af ggplot2-pakken i R.

Til dette indlæser vi først de relevante pakker. Derudover definerer jeg en geokodningsfunktion ved hjælp af Open Street Map API. Geokodningsfunktionen blev kopieret fra https://datascienceplus.com/osm-nominatim-with-r-getting-locations-geo-coordinates-by-its-address/.

I det næste trin anvender jeg map_data-funktionen fra ggplot2-pakken for at indlæse et kort i form af en dataramme. Kortet . Vi opretter derefter også en liste over bynavne manuelt og konverterer disse bynavne til koordinater ved hjælp af geokodningsfunktionen.

Ved hjælp af ggplot2 visualiseres de geokodede datapunkter på et verdenskort.

# importer relevante biblioteker
library(ggplot2)
library(ggrepel)

# osm geocoder
osm_geocoder <- function(address = NULL)
{
  if(suppressWarnings(is.null(address)))
    return(data.frame())
  
  tryCatch(
    d <- jsonlite::fromJSON( 
      gsub('\\@addr\\@', gsub('\\s+', '\\%20', address), 
           'http://nominatim.openstreetmap.org/search/@addr@?format=json&addressdetails=0&limit=1')
    ), error = function(c) return(data.frame())
  )
  
  if(length(d) == 0) 
    return(data.frame())
  
  return(data.frame(lon = as.numeric(d$lon), lat = as.numeric(d$lat)))
}

# kerne del
# dataforberedelse
world <- map_data("world")

points_df <- data.frame(matrix(nrow=0,ncol=2))
colnames(points_df) <- c("long","lat")
city_list <- c("New York City, New York, USA",
               "Detroit, Michigan, USA",
               "Atlanta, Georgia, USA",
               "Birmingham, Alabama, USA",
               "Miami, Florida, USA",
               "Nashville, Tennessee, USA",
               "Knoxville, Tennessee, USA")

for(i in city_list){
  dum_coord <- osm_geocoder(i)
  points_df <- rbind(points_df,data.frame("long"=c(as.numeric(dum_coord[1])),"lat"=c(as.numeric(dum_coord[2]))))
}

gg_1 <- ggplot(world) +
  geom_polygon(mapping=aes(x=long,y=lat, group=group),fill="white",colour="black")

gg_dots <- gg_1 + geom_point(data=points_df,mapping=aes(x=long, y=lat),colour="red",size=3)

gg_text <- gg_dots + geom_text_repel(data=data.frame("long"=points_df$long,"lat"=points_df$lat,"city"=city_list),
                                                                    mapping=aes(x=long,y=lat,label=c("New York",
                                                                                                     "Detroit",
                                                                                                     "Atlanta",
                                                                                                     "Birmingham",
                                                                                                     "Miami",
                                                                                                     "Nashville",
                                                                                                     "Knoxville")),
                                     color="blue",
                                     size=6)

# output the plot
gg_text

You May Also Like

Leave a Reply

Leave a Reply

Din e-mailadresse vil ikke blive publiceret. Krævede felter er markeret med *

This site uses Akismet to reduce spam. Learn how your comment data is processed.