Remove aggregate_diskarray way as well

Author: felixcremer

src/auxil.jl

@@ -154,7 +154,6 @@ function gdalcube(filenames::AbstractVector{<:AbstractString}, stackgroups=:lazy
 
     #@show sdates
     # Put the dates which are 200 seconds apart into groups
-    if stackgroups in [:dae, :lazyagg]
         groupinds = grouptimes(sdates, 200000)
         onefile = first(sfiles)
         gd = backendlist[:gdal]
@@ -173,13 +172,8 @@ function gdalcube(filenames::AbstractVector{<:AbstractString}, stackgroups=:lazy
 
         cubelist = CFDiskArray.(group_gdbs, (gdbattrs,))
         stackinds = stackindices(sdates)
-        aggdata = if stackgroups == :dae
-            gcube = diskstack(cubelist)
-            aggdata = DAE.aggregate_diskarray(gcube, mean ∘ skipmissing, (3 => stackinds,); strategy=:direct)
-        else
-            println("Construct lazy diskarray")
-            LazyAggDiskArray(skipmissingmean, cubelist, stackinds)
-        end
+        println("Construct lazy diskarray")
+        LazyAggDiskArray(skipmissingmean, cubelist, stackinds)
         #    data = DiskArrays.ConcatDiskArray(reshape(groupcubes, (1,1,length(groupcubes))))
         dates_grouped = [sdates[group[begin]] for group in groupinds]
 

src/main.jl

@@ -77,7 +77,7 @@ function main(;
     orbit="D",
     threshold=3.0,
     folders=["V1M0R1", "V1M1R1", "V1M1R2"],
-    stack=:dae,
+    stack=:lazyagg,
     delete_intermediate=false
 )
     outdir = Path(outdir)

test/testdata.jl

@@ -21,7 +21,8 @@
         indir=indir,
         start_date=Date("2021-01-01"),
         end_date=Date("2022-01-01"),
-        outdir=Path(outdir)
+        outdir=Path(outdir),
+        stack=:lazyagg,
     )
     a = open_dataset(joinpath(outdir, "E051N018T3_rqatrend_VH_D022_thresh_3.0_2021-01-01_2022-01-01.zarr")).layer