Add open_tree function to load a ZarrDataset as a DimTree

Author: felixcremer

Project.toml

@@ -4,10 +4,24 @@ authors = ["Felix Cremer <fcremer@bgc-jena.mpg.de> and contributors"]
 version = "0.1.0-DEV"
 
 [deps]
+CommonDataModel = "1fbeeb36-5f17-413c-809b-666fb144f157"
 DimensionalData = "0703355e-b756-11e9-17c0-8b28908087d0"
+DiskArrayTools = "fcd2136c-9f69-4db6-97e5-f31981721d63"
+DiskArrays = "3c3547ce-8d99-4f5e-a174-61eb10b00ae3"
+Rasters = "a3a2b9e3-a471-40c9-b274-f788e487c689"
+YAXArrayBase = "90b8fcef-0c2d-428d-9c56-5f86629e9d14"
+YAXArrays = "c21b50f5-aa40-41ea-b809-c0f5e47bfa5c"
 Zarr = "0a941bbe-ad1d-11e8-39d9-ab76183a1d99"
+ZarrDatasets = "519a4cdf-1362-424a-9ea1-b1d782dbb24b"
 
 [compat]
+CommonDataModel = "0.3.8"
 DimensionalData = "0.29.12"
+DiskArrayTools = "0.1.12"
+DiskArrays = "0.4.11"
+Rasters = "0.14.4"
+YAXArrayBase = "0.7.5"
+YAXArrays = "0.6.1"
 Zarr = "0.9.4"
+ZarrDatasets = "0.1.3"
 julia = "1.6.7"

src/SentinelDataSource.jl

@@ -1,11 +1,38 @@
 module SentinelDataSource
-using DimensionalData
+using DimensionalData: DimTree
 using Zarr: zopen
+using ZarrDatasets: ZarrDataset
+using CommonDataModel: CommonDataModel as CDM
 
 function open_eopf(path)
-    z = zopen(path)
+    zd = ZarrDataset(path)
+    eopfstem = DimTree()
+    eopfstem.measurements = DimTree()
+end
+
+function open_tree(dataset)
+    stem = DimTree()
+    groupnames = CDM.groupnames(dataset)
+    varnames = CDM.varnames(dataset)
+    alldimnames = nesteddimnames(dataset)
+    for v in setdiff(varnames, alldimnames)
+        @show v
+        setproperty!(stem, Symbol(v), Raster(CDM.variable(dataset, v), lazy=true))
+    end
+    for g in groupnames
+        setproperty!(stem, Symbol(g), open_tree(CDM.group(dataset, g)))
+    end
+    stem
 end
 #zopen()
 # Write your package code here.
 
 end
+
+function nesteddimnames(zarrdataset)
+    alldims = []
+    for v in CDM.varnames(zarrdataset)
+        append!(alldims, CDM.dimnames(CDM.variable(zarrdataset, v)))
+    end
+    unique(alldims)
+end