need CONVERT_TO_LINEAR -nlt for geometry read

[mdsumner]

or general conversion from as yet unsupported types in wk

pl <- place::place("Hobart", wh = 1.852e3 * 9, family = "EPSG:3577")

dea <- "WFS:http://geoserver.dea.ga.gov.au/geoserver/wfs?OUTPUT_FORMAT=SHAPE-ZIP"
library(vapour)
layers <- vapour_layer_names(dea)
fields <- vapour_read_fields(dea, extent = pl$extent, 
                             sql = sprintf("SELECT * FROM \"%s\"", layers[1]))
g <- vapour_read_geometry(dea, extent = pl$extent, 
                        sql = sprintf("SELECT geom FROM \"%s\"", layers[1]))

plot(wk::as_wkb(sf::st_cast(sf::st_as_sfc(g), "MULTILINESTRING")))

[mdsumner]