Esplorazione dei dati

Introduzione

La funzione getCapability() permette innanzitutto di restituire le diverse coverage offerte dal servizio di Rasdaman di ARPA LOMBARDIA.

[Utilizza la funzione GetCapabilities del WCS]


Attraverso le seguenti funzioni è invece possibile interrogare i metadati (ottenere descrizione e alcune proprietà) di una specifica coverage selezionata:

[Utilizzano la funzione DescribeCoverage del servizio WCS]

Importazione delle librerie necessarie

library(myCubeR)

Utilizzo delle funzioni implementate

Richiedere l’elenco delle Coverage/Servizi WCS offerti dal server Rasdaman:

coverages<-getCapability()
coverages
#[1]  "rh_ana"    "rh_idi"    "t2m_idi"   "enne_mean" "t2m_ana"   "prec_ana"  "pwet_idi"  "bui"       "fwi"       "ffmc"      "pdry_idi" 
#[12] "dmc"       "rh_hdx"    "isi"       "dc"        "t2m_bkg"

Scegliere una coverage (in questo caso rh_ana):

coverage<-getCapability()[1]
coverage
#[1] "rh_ana"

Per la coverage selezionata (rh_ana) è possibile:

Richiedere gli “assi” che descrivono i Datacube (solitamente X,Y e come terza variabile/asse il tempo):

axes<-coverage_get_coordsys(coverage = coverage)
axes
#[1] "ansi" "E"    "N"

Richiedere il sistema di riferimento della coverage

ref_sys=coverage_get_coordinate_reference(coverage=coverage)
ref_sys
#[1] "32632"

Richiedere “l’estensione temporale” della Coverage

temp_ext=coverage_get_temporal_extent(coverage=coverage)
temp_ext
#[1] "2020-05-18T12:00:00.000Z" "2020-11-10T06:00:00.000Z"

Richiedere l’entensione (BBOX) spaziale della Coverage

bb_box=coverage_get_bounding_box(coverage=coverage)
bb_box
#[1] "436000"  "687000"  "4918000" "5167000"

Richiedere la risoluzione spaziale della Coverage

res=coverage_get_resolution(coverage=coverage)
res
#[1] 1000

Richiedere il nome delle bande (di default per Rasdaman di ARPA LOMBARDIA la variabile principale viene nominata field_1. Es per rh_ana corrisponde a……)

bands<-coverage_get_bands(coverage = coverage)
bands
#[1] "field_1"

Richiedere tutte i “timestamp”/le date per cui sono disponibili i raster che costruiscono il Datacube

time_stamp=coverage_get_timestamps(coverage=coverage)
time_stamp
#[1] "2020-05-18T12:00:00.000Z" "2020-05-18T13:00:00.000Z" "2020-05-18T14:00:00.000Z" "2020-05-18T15:00:00.000Z" "2020-05-18T16:00:00.000Z"
#[6] "2020-05-18T17:00:00.000Z" "2020-05-18T18:00:00.000Z" "2020-05-18T19:00:00.000Z" "2020-05-18T20:00:00.000Z" "2020-05-18T21:00:00.000Z"
# ................

print(paste("Un totale di",length(time_stamp), "dati"))
#[1] "Un totale di 4133 dati"