The goal of ozmaps is to get maps of Australia to plot!
Install from CRAN with
install.packages("ozmaps")
The development version of ozmaps may be installed directly from github.
devtools::install_github("mdsumner/ozmaps")
The package includes some simple features data, which can be used
independently of ozmaps with the sf
package. If required, install sf
from CRAN.
NOTE: Since April 2020) we must have sf installed because of requirements of the tibble package.
If sf causes you problems or you can’t work with it get in touch and I will help you.
install.packages("sf")
Plot Australia with states.
library(ozmaps)
ozmap()
Plot Australia without states.
ozmap(x = "country")
Add to an existing plot.
plot(quakes[c("long", "lat")], xlim = c(120, 190))
ozmap(add = TRUE)
Obtain the data used in sf
form.
sf_oz <- ozmap_data("states")
tibble::as_tibble(sf_oz)
#> # A tibble: 9 x 2
#> NAME geometry
#> <chr> <MULTIPOLYGON [°]>
#> 1 New South Wales (((150.7016 -35.12286, 150.6611 -35.11782, 150.6373 -35.…
#> 2 Victoria (((146.6196 -38.70196, 146.6721 -38.70259, 146.6744 -38.…
#> 3 Queensland (((148.8473 -20.3457, 148.8722 -20.37575, 148.8515 -20.3…
#> 4 South Australia (((137.3481 -34.48242, 137.3749 -34.46885, 137.3805 -34.…
#> 5 Western Australia (((126.3868 -14.01168, 126.3625 -13.98264, 126.3765 -13.…
#> 6 Tasmania (((147.8397 -40.29844, 147.8902 -40.30258, 147.8812 -40.…
#> 7 Northern Territory (((136.3669 -13.84237, 136.3339 -13.83922, 136.3532 -13.…
#> 8 Australian Capital … (((149.2317 -35.222, 149.2346 -35.24047, 149.2716 -35.27…
#> 9 Other Territories (((167.9333 -29.05421, 167.9188 -29.0344, 167.9313 -29.0…
Plot with a custom palette.
library(sf)
if (utils::packageVersion("paletteer") < '1.0.0') {
pal <- paletteer::paletteer_d(package = "ochRe", palette = "namatjira_qual")
} else {
pal <- paletteer::paletteer_d(palette = "ochRe::namatjira_qual")
}
opal <- colorRampPalette(pal)
nmjr <- opal(nrow(sf_oz))
plot(st_geometry(sf_oz), col = nmjr)
## plot directly with ggplot2
library(ggplot2)
ggplot(sf_oz, aes(fill = NAME)) + geom_sf() + coord_sf(crs = "+proj=lcc +lon_0=135 +lat_0=-30 +lat_1=-10 +lat_2=-45 +datum=WGS84") + scale_fill_manual(values = nmjr)
Plot the ABS layers (from 2016).
ozmap("abs_ced", col = opal(nrow(abs_ced)))
ozmap("abs_ste", col = opal(nrow(abs_ste)))
These ABS layers abs_ced
, abs_lga
, and abs_ste
are derived from
the 2016 sources and simplified using rmapshaper::ms_simplify(, keep = 0.05, keep_shapes = TRUE)
so all the original polygons are there.
There is sufficient detail to map many (most?) of the regions on their
own, which was a major goal for this package.
The cache of the source data at original resolution is available in ozmaps.data.
Compare the detail of Bruny Island here in this box, compared with the very basic maps package layer.
library(dplyr)
kbor <- abs_lga %>% dplyr::filter(grepl("Kingborough", NAME))
bb <- st_bbox(kbor)
layout(matrix(c(1, 1, 1, 2, 2, 2, 2, 2, 2), nrow = 3))
plot(kbor, reset = FALSE, main = "Kingborough (TAS)")
rect(bb["xmin"], bb["ymin"], bb["xmax"], bb["ymax"])
library(mapdata)
#> Loading required package: maps
par(mar = rep(0, 4))
plot(c(145, 148.5), c(-43.6, -40.8), type = "n", asp = 1/cos(mean(bb[c(2, 4)]) * pi/180), axes = FALSE, xlab = "", ylab = "")
maps::map(database = "worldHires", regions = "australia", xlim = c(145, 148.5), ylim = c(-43.6, -40.8), add = TRUE)
rect(bb["xmin"], bb["ymin"], bb["xmax"], bb["ymax"])
Please note that the ‘ozmaps’ project is released with a Contributor Code of Conduct. By contributing to this project, you agree to abide by its terms.