@@ -22,115 +22,15 @@ function __init__()
2222
2323 backendregex = Pair[]
2424
25- @require DimensionalData= " 0703355e-b756-11e9-17c0-8b28908087d0" begin
26- using . DimensionalData: DimensionalArray, DimensionalData, data, Dim, metadata
27-
28- _dname (:: DimensionalData.Dim{N} ) where N = N
29- dimname (x:: DimensionalArray , i) = _dname (DimensionalData. dims (x)[i])
30-
31-
32- dimvals (x:: DimensionalArray ,i) = DimensionalData. dims (x)[i]. val
33-
34- getdata (x:: DimensionalArray ) = data (x)
35-
36- getattributes (x:: DimensionalArray ) = metadata (x)
37-
38- function yaxcreate (:: Type{<:DimensionalArray} ,data,dnames,dvals,atts)
39- d = ntuple (ndims (data)) do i
40- Dim {Symbol(dnames[i])} (dvals[i])
41- end
42- DimensionalArray (data,d,metadata = atts)
43- end
44- end
45-
46- @require AxisArrays= " 39de3d68-74b9-583c-8d2d-e117c070f3a9" begin
47- using . AxisArrays: AxisArrays, AxisArray
48-
49- dimname (a:: AxisArray , i) = AxisArrays. axisnames (a)[i]
50- dimnames (a:: AxisArray ) = AxisArrays. axisnames (a)
51- dimvals (a:: AxisArray , i) = AxisArrays. axisvalues (a)[i]
52- iscontdim (a:: AxisArray , i) = AxisArrays. axistrait (AxisArrays. axes (a,i)) <: AxisArrays.Dimensional
53- getdata (a:: AxisArray ) = parent (a)
54- function yaxcreate (:: Type{<:AxisArray} , data, dnames, dvals, atts)
55- d = ntuple (ndims (data)) do i
56- dnames[i] => dvals[i]
57- end
58- AxisArray (data; d... )
59- end
60- end
61-
62- @require AxisIndices= " f52c9ee2-1b1c-4fd8-8546-6350938c7f11" begin
63- using . AxisIndices: AbstractAxis,AxisIndices
64-
65- valfromaxis (ax:: AbstractAxis ) = keys (ax)
66-
67- getdata (a:: AxisIndices.AxisArray ) = parent (a)
68-
69- yaxcreate (:: Type{<:AxisIndices.AxisArray} , data, dnames, dvals, atts) =
70- AxisIndices. AxisArray (data, map (i-> dvals[i], 1 : ndims (data))... )
71- end
72-
73- @require ArchGDAL= " c9ce4bd3-c3d5-55b8-8973-c0e20141b8c3" begin
74- import . ArchGDAL: RasterDataset, AbstractRasterBand,
75- getgeotransform, width, height, getname, getcolorinterp,
76- getband, nraster, getdataset
77- function dimname (a:: RasterDataset , i)
78- if i == 1
79- return :Y
80- elseif i == 2
81- return :X
82- elseif i == 3
83- return :Band
84- else
85- error (" RasterDataset only has 3 dimensions" )
86- end
87- end
88- function dimvals (a:: RasterDataset , i)
89- if i == 1
90- geo= getgeotransform (a)
91- latr = range (geo[1 ],length= width (a), step= geo[2 ])
92- elseif i == 2
93- geo= getgeotransform (a)
94- range (geo[4 ],length= height (a), step= geo[6 ])
95- elseif i == 3
96- colnames = map (ib -> getname (getcolorinterp (getband (a,ib))),1 : nraster (a))
97- if ! allunique (colnames)
98- colnames = string .(" Band_" ,1 : nraster (a))
99- end
100- colnames
101- else
102- error (" RasterDataset only has 3 dimensions" )
103- end
104- end
105- iscontdim (a:: RasterDataset , i) = i < 3 ? true : nraster (a)< 8
106- getattributes (a:: RasterDataset ) =
107- Dict {String,Any} (" projection" => ArchGDAL. toPROJ4 (ArchGDAL. newspatialref (ArchGDAL. getproj (a))))
108-
109-
110- function dimname (:: AbstractRasterBand , i)
111- if i == 1
112- return :Y
113- elseif i == 2
114- return :X
115- else
116- error (" RasterDataset only has 3 dimiensions" )
117- end
118- end
119- function dimvals (b:: AbstractRasterBand , i)
120- geo = getgeotransform (getdataset (b))
121- if i == 1
122- range (geo[1 ],length= width (b), step= geo[2 ])
123- elseif i == 2
124- range (geo[4 ],length= height (b), step= geo[6 ])
125- else
126- error (" RasterDataset only has 3 dimiensions" )
127- end
128- end
129- iscontdim (a:: AbstractRasterBand , i) = true
130- getattributes (a:: AbstractRasterBand ) =
131- getattributes (ArchGDAL. RasterDataset (ArchGDAL. getdataset (a)))
132-
133- end
25+ @require NamedDims= " 356022a1-0364-5f58-8944-0da4b18d706f" include (" axisarrays/nameddims.jl" )
26+
27+ @require DimensionalData= " 0703355e-b756-11e9-17c0-8b28908087d0" include (" axisarrays/dimensionaldata.jl" )
28+
29+ @require AxisArrays= " 39de3d68-74b9-583c-8d2d-e117c070f3a9" include (" axisarrays/axisarrays.jl" )
30+
31+ @require AxisIndices= " f52c9ee2-1b1c-4fd8-8546-6350938c7f11" include (" axisarrays/axisindices.jl" )
32+
33+ @require ArchGDAL= " c9ce4bd3-c3d5-55b8-8973-c0e20141b8c3" include (" axisarrays/archgdal.jl" )
13434
13535 @require Zarr= " 0a941bbe-ad1d-11e8-39d9-ab76183a1d99" include (" datasets/zarr.jl" )
13636
0 commit comments