Working with shapefiles
[1]:
!wget https://gisco-services.ec.europa.eu/distribution/v2/nuts/shp/NUTS_RG_20M_2021_4326_LEVL_0.shp.zip
--2024-12-16 19:22:58-- https://gisco-services.ec.europa.eu/distribution/v2/nuts/shp/NUTS_RG_20M_2021_4326_LEVL_0.shp.zip
Resolving gisco-services.ec.europa.eu (gisco-services.ec.europa.eu)... 13.69.248.157
Connecting to gisco-services.ec.europa.eu (gisco-services.ec.europa.eu)|13.69.248.157|:443... connected.
HTTP request sent, awaiting response... 200 OK
Length: 93608 (91K) [application/zip]
Saving to: ‘NUTS_RG_20M_2021_4326_LEVL_0.shp.zip.2’
NUTS_RG_20M_2021_43 100%[===================>] 91.41K --.-KB/s in 0.07s
2024-12-16 19:22:59 (1.37 MB/s) - ‘NUTS_RG_20M_2021_4326_LEVL_0.shp.zip.2’ saved [93608/93608]
[2]:
import earthkit.plots
chart = earthkit.plots.Map(domain="Europe")
chart.coastlines()
chart.shapes(
"NUTS_RG_20M_2021_4326_LEVL_0.shp.zip",
labels="NUTS_ID", # Label taken directly from shapefile metadata
adjust_labels=True,
)
chart.gridlines(draw_labels=False)
chart.title("Nomenclature of Territorial Units for Statistics (NUTS) Level 0")
chart.show()
