Map Australian cities - R spatial

Question

I want to draw a map of Australia and represent each city as a dot. Then highlight the cities with a high population (>1M)

library(sp)
library(maps)
data(canada.cities)
head(canada.cities)

I have checked the sp package where this can be done for Canada and some other countries. But Australia details are not there. Is there a special way to get the data for a country we like (name of cities, long, lat, pop)?

Yeah I have tried that but it doesn't have the cities or any population details — Jay
world.cities[world.cities$country.etc == "Australia",] and map.cities(country = "Australia", minpop=1e6) — thelatemail
Thank you. I got the necessary data using world.cities. Does map.cities() plot a map. I did following to mark the cities Aus_cities= world.cities[world.cities$country.etc == "Australia",] class(Aus_cities) plot(Aus_cities[Aus_cities@data$pop>1000000,],pch=16,col="red",cex=3,add=TRUE) — Jay

SymbolixAU SymbolixAU · Accepted Answer · 2016-06-23T01:52:04

Now you have the data using world.cities, you can plot them a few ways

library(maps)
df <- world.cities[world.cities$country.etc == "Australia",]

Basic plot of points

plot(df[, c("long", "lat")])

on a ggmap

library(ggmap)

myMap <- get_map(location = "Australia", zoom = 4)

ggmap(myMap) +
geom_point(data = df[, c("long","lat", "pop")], aes(x=long, y = lat, colour = pop > 1000000))

On a leaflet map

library(leaflet)

## define a palette for hte colour
pal <- colorNumeric(palette = "YlOrRd",
                    domain = df$pop)

leaflet(data = df) %>%
    addTiles() %>%
    addCircleMarkers(lat = ~lat, lng = ~long, popup = ~name, 
                     color = ~pal(pop), stroke = FALSE, fillOpacity = 0.6) %>%
    addLegend(position = "bottomleft", pal = pal, values = ~pop)

Map Australian cities - R spatial

1 Answers