get_tilecodes: Intersecting Albers Australian Tile Codes

In sustainablefarms/sflddata: Data Extraction Methods for the Linking Data Project

Intersecting Albers Australian Tile Codes

Description

Given a spatial object returns the tile codes for Albers Tiles used by Geoscience Australia and others

Usage

get_tilecodes(spobj)

tiles_as_sf(tilecodes)

Arguments

spobj	an sp object
tilecodes	A list of tile codes.

Value

a named vector of tile codes

Functions

• tiles_as_sf: Make sf polygons showing the tiles

Examples

tilecodes <- unlist(read.csv(system.file("austilecodes.txt", package = "sflddata")))
names(tilecodes) <- NULL
tiles <- tiles_as_sf(tilecodes)
ausstates <- readRDS("./private/data/basemaps/ausstates.rds") 
sf::st_crs(ausstates) <- sf::st_crs(ausstates)
ausstates <- sf::st_transform(ausstates, sf::st_crs(tiles))
# check that tiles cover all land mass
library(ggplot2)
ggplot() +
  geom_sf(aes(), fill = "red", data = ausstates) +
  geom_sf(aes(), data = tiles)

sustainablefarms/sflddata documentation built on April 19, 2022, 11:19 a.m.