JuliaDataCubes · Feb 21, 2023 · Apr 21, 2023 · Jul 31, 2024 · Aug 1, 2024 · Aug 1, 2024
diff --git a/.github/dependabot.yml b/.github/dependabot.yml
@@ -0,0 +1,7 @@
+# https://docs.github.com/github/administering-a-repository/configuration-options-for-dependency-updates
+version: 2
+updates:
+  - package-ecosystem: "github-actions"
+    directory: "/" # Location of package manifests
+    schedule:
+      interval: "weekly"
diff --git a/.github/workflows/CI.yml b/.github/workflows/CI.yml
@@ -13,21 +13,22 @@ jobs:
       fail-fast: false
       matrix:
         version:
-          - '1.6'
+          - '1.10'
           - '1'
+          - 'nightly'
         os:
           - ubuntu-latest
           - macOS-latest
           - windows-latest
         arch:
           - x64
     steps:
-      - uses: actions/checkout@v2
-      - uses: julia-actions/setup-julia@v1
+      - uses: actions/checkout@v4
+      - uses: julia-actions/setup-julia@v2
         with:
           version: ${{ matrix.version }}
           arch: ${{ matrix.arch }}
-      - uses: actions/cache@v1
+      - uses: actions/cache@v4
         env:
           cache-name: cache-artifacts
         with:
@@ -40,6 +41,6 @@ jobs:
       - uses: julia-actions/julia-buildpkg@v1
       - uses: julia-actions/julia-runtest@v1
       - uses: julia-actions/julia-processcoverage@v1
-      - uses: codecov/codecov-action@v1
+      - uses: codecov/codecov-action@v5
         with:
           file: lcov.info
diff --git a/.gitignore b/.gitignore
@@ -1,5 +1,7 @@
+.vscode
 *.jl.cov
 *.jl.*.cov
 *.jl.mem
 /deps/deps.jl
 /docs/build
+Manifest.toml
diff --git a/Manifest.toml b/Manifest.toml
diff --git a/Project.toml b/Project.toml
@@ -1,16 +1,38 @@
 name = "YAXArrayBase"
 uuid = "90b8fcef-0c2d-428d-9c56-5f86629e9d14"
 authors = ["Fabian Gans <fgans@bgc-jena.mpg.de>"]
-version = "0.6.1"
+version = "0.7.6"
 
 [deps]
 DataStructures = "864edb3b-99cc-5e75-8d2d-829cb0a9cfe8"
 Dates = "ade2ca70-3891-5945-98fb-dc099432e06a"
 Downloads = "f43a241f-c20a-4ad4-852c-f6b1247861c6"
-Requires = "ae029012-a4dd-5104-9daa-d747884805df"
+
 
 [compat]
-DataStructures = "0.17, 0.18"
-Requires = "1"
-julia = "1.6"
+DataStructures = "0.17,0.18"
+julia = "1.9"
+ArchGDAL = "0.10"
+AxisArrays = "0.4"
+AxisKeys = "0.2"
+DimensionalData = "0.27, 0.28, 0.29"
+NetCDF = "0.11, 0.12"
+Zarr = "0.8, 0.9"
+
+[extensions]
+ArchGDALExt = "ArchGDAL"
+AxisArraysExt = "AxisArrays"
+AxisKeysExt = "AxisKeys"
+DimensionalDataExt = "DimensionalData"
+NamedDimsExt = "NamedDims"
+NetCDFExt = "NetCDF"
+ZarrExt = "Zarr"
 
+[weakdeps]
+ArchGDAL = "c9ce4bd3-c3d5-55b8-8973-c0e20141b8c3"
+AxisArrays = "39de3d68-74b9-583c-8d2d-e117c070f3a9"
+AxisKeys = "94b1ba4f-4ee9-5380-92f1-94cde586c3c5"
+DimensionalData = "0703355e-b756-11e9-17c0-8b28908087d0"
+NamedDims = "356022a1-0364-5f58-8944-0da4b18d706f"
+NetCDF = "30363a11-5582-574a-97bb-aa9a979735b9"
+Zarr = "0a941bbe-ad1d-11e8-39d9-ab76183a1d99"
diff --git a/README.md b/README.md
@@ -1,12 +1,10 @@
 # YAXArrayBase.jl
 
-![Lifecycle](https://img.shields.io/badge/lifecycle-maturing-blue.svg)<!--
-![Lifecycle](https://img.shields.io/badge/lifecycle-stable-green.svg)
-![Lifecycle](https://img.shields.io/badge/lifecycle-retired-orange.svg)
-![Lifecycle](https://img.shields.io/badge/lifecycle-archived-red.svg)
-![Lifecycle](https://img.shields.io/badge/lifecycle-dormant-blue.svg) -->
-[![Build Status](https://travis-ci.com/JuliaDataCubes/YAXArrayBase.jl.svg?branch=master)](https://travis-ci.com/JuliaDataCubes/YAXArrayBase.jl)
-[![codecov.io](http://codecov.io/github/JuliaDataCubes/YAXArrayBase.jl/coverage.svg?branch=master)](http://codecov.io/github/JuliaDataCubes/YAXArrayBase.jl?branch=master)
+![Lifecycle](https://img.shields.io/badge/lifecycle-maturing-blue.svg)
+[![][codecov-img]][codecov-url]
+
+[codecov-img]: https://codecov.io/gh/JuliaDataCubes/YAXArrayBase.jl/branch/master/graph/badge.svg
+[codecov-url]: https://codecov.io/gh/JuliaDataCubes/YAXArrayBase.jl
 
 # YAXArrayBase
 

diff --git a/ext/ArchGDALExt/ArchGDALExt.jl b/ext/ArchGDALExt/ArchGDALExt.jl
@@ -0,0 +1,99 @@
+module ArchGDALExt
+    import ArchGDAL: RasterDataset, AbstractRasterBand, getgeotransform, width, height
+    import ArchGDAL: getname, getcolorinterp, getband, nraster, getdataset
+    using ArchGDAL: ArchGDAL as AG
+
+    import ArchGDAL.DiskArrays: GridChunks, eachchunk
+    import ArchGDAL.DiskArrays
+
+    import YAXArrayBase: dimname, dimnames, dimvals, iscontdim, getattributes, getdata, yaxcreate
+    import YAXArrayBase: YAXArrayBase as YAB
+    using DataStructures: OrderedDict
+
+    include("archgdaldataset.jl")
+
+    function dimname(a::RasterDataset, i)
+        if i == 1
+            return :Y
+        elseif i == 2
+            return :X
+        elseif i == 3
+            return :Band
+        else
+            error("RasterDataset only has 3 dimensions")
+        end
+    end
+    function dimvals(a::RasterDataset, i)
+        if i == 1
+            geo=getgeotransform(a)
+            latr = range(geo[1],length=width(a), step=geo[2])
+        elseif i == 2
+            geo=getgeotransform(a)
+            range(geo[4],length=height(a), step=geo[6])
+        elseif i == 3
+            colnames = map(ib -> getname(getcolorinterp(getband(a,ib))),1:nraster(a))
+            if !allunique(colnames)
+                colnames = string.("Band_",1:nraster(a))
+            end
+            colnames
+        else
+            error("RasterDataset only has 3 dimensions")
+        end
+    end
+
+    iscontdim(a::RasterDataset, i) = i < 3 ? true : nraster(a)<8
+    function getattributes(a::RasterDataset)
+        globatts = Dict{String,Any}(
+            "projection_PROJ4"=>AG.toPROJ4(AG.newspatialref(AG.getproj(a))),
+            "projection_WKT"=>AG.toWKT(AG.newspatialref(AG.getproj(a))),
+        )
+        bands = (getbandattributes(AG.getband(a, i)) for i in 1:size(a, 3))
+        allbands = mergewith(bands...) do a1,a2
+            isequal(a1,a2) ? a1 : missing
+        end
+        return merge(globatts, allbands)
+    end
+
+
+    function dimname(::AbstractRasterBand, i)
+        if i == 1
+            return :Y
+        elseif i == 2
+            return :X
+        else
+            error("RasterDataset only has 3 dimensions")
+        end
+    end
+    function dimvals(b::AbstractRasterBand, i)
+        geo = getgeotransform(getdataset(b))
+        if i == 1
+            range(geo[1],length=width(b), step=geo[2])
+        elseif i == 2
+            range(geo[4],length=height(b), step=geo[6])
+        else
+            error("RasterDataset only has 3 dimensions")
+        end
+    end
+    iscontdim(a::AbstractRasterBand, i) = true
+    function getattributes(a::AbstractRasterBand)
+        atts = getattributes(AG.RasterDataset(AG.getdataset(a)))
+        bandatts = getbandattributes(a)
+        return merge(atts, bandatts)
+    end
+
+    function insertattifnot!(attrs, val, name, condition)
+        if !condition(val)
+            attrs[name] = val
+        end
+    end
+    function getbandattributes(a::AbstractRasterBand)
+        atts = Dict{String,Any}()
+        catdict = Dict((i-1)=>v for (i,v) in enumerate(AG.getcategorynames(a)))
+        insertattifnot!(atts, AG.getnodatavalue(a), "missing_value", isnothing)
+        insertattifnot!(atts, catdict, "labels", isempty)
+        insertattifnot!(atts, AG.getunittype(a), "units", isempty)
+        insertattifnot!(atts, AG.getoffset(a), "add_offset", iszero)
+        insertattifnot!(atts, AG.getscale(a), "scale_factor", x->isequal(x, one(x)))
+        return atts
+    end
+end