Wanted: SF EEZ layers in SOmap Data

[Maschette]

It would be good to have an SF version of the EEZ layers, probably separate in the first instance or stuff may break. For both of them it would be great if we can add some level of ID to them to make it easier to filter to set EEZ's I cant seem to see one in the data.

[mdsumner]

here's a way

dsn <- "WFS:https://geo.vliz.be/geoserver/MarineRegions/wfs?Request=getCapabilities"
layers <- sf::st_layers(dsn)
layers$name      layers$geomtype  layers$driver    layers$features  layers$fields    layers$crs
grep("eez", layers$name, value = TRUE)
[1] "MarineRegions:eez_archipelagic_waters"
[2] "MarineRegions:eez_24nm"
[3] "MarineRegions:eez"
[4] "MarineRegions:eez_internal_waters"
[5] "MarineRegions:eez_land"
[6] "MarineRegions:eez_iho_union_v2"
[7] "MarineRegions:eez_boundaries"
[8] "MarineRegions:eez_12nm"
[9] "MarineRegions:eez_iho"

eez <- sf::read_sf(dsn, "MarineRegions:eez")

That's equivalent to what {mregions} would get I think.

ouch, hang on - this will present in the criminally-imagined form of a "MultiSurface", when it's nothing of the sort (a kind of esri gate-keeping). I'll extract and make it useable ...

[mdsumner]

ok it's also quite large, so I'll have to look at this in depth a bit later

notes for self:

ogr2ogr eez.fgb -f FlatGeobuf -nlt MULTIPOLYGON  "WFS:https://geo.vliz.be/geoserver/MarineRegions/wfs?Request=getCapabilities" "MarineRegions:eez"