We use leaflet for R to locate the three national parks in Washington State, USA: Olympic NP, Mount Rainier NP and North Cascades NP. Let’s build a data frame with the relevant information.

np <- data.frame(longitude = c(-123.6044,-121.7269,-121.2985), 
                 latitude = c(47.8021,46.8800,48.7718), 
                 area = c(1442,369.3,789) )

We also want the National Park Service logo as a marker.

library(leaflet)
icon <- makeIcon(iconUrl = "https://upload.wikimedia.org/wikipedia/commons/1/1d/US-NationalParkService-Logo.svg",
                 iconWidth = 40*720/553, iconHeight = 40)

We then define the corresponding string. When clicking on the marker, the name of the park will popup as an hyperlink.

sites <- c("<a href = 'https://www.nps.gov/olym/index.htm'> Olympic NP </a>",
           "<a href = 'https://www.nps.gov/mora/index.htm'> Mount Rainier NP </a>",
           "<a href = 'https://www.nps.gov/noca/index.htm'> North Cascades NP </a>")

Here is the map ! The blue circles are proportional to area of the park.

np %>% leaflet() %>% 
    addTiles() %>% 
    addMarkers(popup = sites, icon = icon) %>% 
    addCircles(weight = 1, radius = np$area * 25)