knitr::opts_chunk$set( collapse = TRUE, comment = "#>", fig.path = "man/figures/README-", out.width = "100%" )
This is a dead end, the WFS server stuff is bs by the looks.
The goal of thelist is to read data from the Tasmania TheLIST service, for local authority data layers.
You can install thelist from github with
remotes::install_github("mdsumner/thelist")
This is a basic example which shows how to find the matching cadastral parcel for a given address in the state of Tasmania, without downloading the data layers or reading every data element. se
library(thelist) library(lazysf) address <- lazysf(thelist_wfs(), thelist_layer("address")) library(dplyr) pid <- address %>% dplyr::filter(ST_NO_FROM == 349, STREET == "ELIZABETH") %>% select(ST_NO_FROM, STREET, PID, everything()) %>% collect() cadastre <- lazysf(thelist_wfs(), thelist_layer("cadastral")) %>% dplyr::filter(PID == !!pid$PID) %>% select(SHAPE) %>% mutate(SHAPE = sql('CAST(SHAPE AS geometry( "GEOMETRYCOLLECTION")) ')) pull(cadastre)
See that there are many layers to choose from.
writeLines(list_layers)
And we can find one of these, by pattern or at random.
(lga_layer <- thelist_layer("government_areas")) thelist_layer(random = TRUE)
We can just read direct with sf. It's pretty slow though, not sf's fault ... still exploring.
library(sf) lga <- read_sf(thelist_wfs(), lga_layer)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.