I'm not an expert for the XML but something like this:
library(xml2)
leg =
read_xml("lcv_landcover.hcl_lucas.corine.rf_p_30m_0..0cm_2000_eumap_epsg3035_v0.1.qml",
)
leg.col =
as_list(xml_children(xml_children(xml_children(leg)[[3]])[[1]])[[3]])
leg.df = do.call(rbind, lapply(leg.col,
Thanks, Tom.
How can one use the .qml file you point to in an R workflow to get
category labels and colors for the COG data you shared originally?
On 16/03/2021 12:24, Tomislav Hengl wrote:
Hi Edzer thanks for spotting this.
Cloud Optimized GeoTiffs seems to be somewhat complicated when it
Hi Tom, great work!
When reading this either with terra or stars, as in
in.tif =
"/vsicurl/http://s3.eu-central-1.wasabisys.com/eumap/lcv/lcv_landcover.hcl_lucas.corine.rf_p_30m_0..0cm_2019_eumap_epsg3035_v0.1.tif;
library(terra)
# terra version 1.1.5
plot(rast(in.tif))
#library(raster)
We have mapped land cover classes for the 2000-2019 period for
continental Europe at 30-m resolution using spatiotemporal Machine
Learning (we used R and python for modeling). Explore the dynamic EU
landscapes on your palm using the ODS-Europe viewer:
https://maps.opendatascience.eu
To