Add leaves with setindex

Author: felixcremer

src/SentinelDataSource.jl

@@ -3,6 +3,7 @@ using DimensionalData: DimTree
 using Zarr: zopen
 using ZarrDatasets: ZarrDataset
 using CommonDataModel: CommonDataModel as CDM
+using Rasters:Raster
 
 function open_eopf(path)
     zd = ZarrDataset(path)
@@ -17,22 +18,24 @@ function open_tree(dataset)
     alldimnames = nesteddimnames(dataset)
     for v in setdiff(varnames, alldimnames)
         @show v
-        setproperty!(stem, Symbol(v), Raster(CDM.variable(dataset, v), lazy=true))
+        setindex!(stem, Raster(CDM.variable(dataset, v), lazy=true),Symbol(v))
     end
     for g in groupnames
-        setproperty!(stem, Symbol(g), open_tree(CDM.group(dataset, g)))
+        setindex!(stem,  open_tree(CDM.group(dataset, g)),Symbol(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
+
 end

test/runtests.jl

@@ -2,6 +2,7 @@ using SentinelDataSource
 using Test
 using Aqua
 using Zarr
+using ZarrDatasets
 
 @testset "SentinelDataSource.jl" begin
     @testset "Code quality (Aqua.jl)" begin
@@ -11,8 +12,8 @@ using Zarr
 
     @testset "S1 L1 GRD" begin
         s1grdpath = "https://objectstore.eodc.eu:2222/e05ab01a9d56408d82ac32d69a5aae2a:sample-data/tutorial_data/cpm_v253/S1A_IW_GRDH_1SDV_20240201T164915_20240201T164940_052368_065517_750E.zarr"
-        s1grd = zopen(s1grdpath)
-        @test s1grd isa ZGroup
+        s1grd = ZarrDataset(s1grdpath)
+        s1grdtree = SentinelDataSource.open_tree(s1grd)
     end
     @testset "S1 L1 SLC" begin
         s1slcpath = "https://objectstore.eodc.eu:2222/e05ab01a9d56408d82ac32d69a5aae2a:sample-data/tutorial_data/cpm_v253/S1A_IW_SLC__1SDV_20231119T170635_20231119T170702_051289_063021_178F.zarr"