Skip to content

Commit 05a9fbb

Browse files
lesserwhirlslesserwhirls
committed
Handle scalar xarray-Zarr arrays
Fixes #1100
1 parent 4927201 commit 05a9fbb

24 files changed

Lines changed: 411 additions & 29 deletions

File tree

‎cdm/core/src/main/java/ucar/nc2/constants/CDM.java‎

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -66,6 +66,7 @@ public class CDM {
6666
public static final String NCPROPERTIES = "_NCProperties";
6767
public static final String ISNETCDF4 = "_IsNetcdf4";
6868
public static final String SUPERBLOCKVERSION = "_SuperblockVersion";
69+
public static final String ARRAYDIMENSIONS = "_ARRAY_DIMENSIONS";
6970

7071
public static final String[] SPECIALS = {NCPROPERTIES, ISNETCDF4, SUPERBLOCKVERSION};
7172

‎cdm/zarr/src/main/java/ucar/nc2/iosp/zarr/ZarrHeader.java‎

Lines changed: 43 additions & 21 deletions
Original file line numberDiff line numberDiff line change
@@ -1,10 +1,12 @@
11
/*
2-
* Copyright (c) 2021 University Corporation for Atmospheric Research/Unidata
2+
* Copyright (c) 2021-2025 University Corporation for Atmospheric Research/Unidata
33
* See LICENSE for license information.
44
*/
55

66
package ucar.nc2.iosp.zarr;
77

8+
import static ucar.nc2.constants.CDM.ARRAYDIMENSIONS;
9+
810
import com.fasterxml.jackson.databind.ObjectMapper;
911

1012
import ucar.ma2.ArrayObject;
@@ -34,6 +36,7 @@ public class ZarrHeader {
3436
private final RandomAccessDirectory rootRaf;
3537
private final Group.Builder rootGroup;
3638
private final String rootLocation;
39+
3740
private static final ObjectMapper objectMapper = new ObjectMapper();
3841

3942
public ZarrHeader(RandomAccessDirectory raf, Group.Builder rootGroup) {
@@ -219,20 +222,25 @@ private void makeVariable(RandomAccessDirectoryItem item, long dataOffset, ZArra
219222

220223
for (Attribute attr : attrs) {
221224
final String attrName = attr.getName();
222-
if ("_ARRAY_DIMENSIONS".equals(attrName)) {
225+
if (ARRAYDIMENSIONS.equals(attrName)) {
223226
try {
224-
final ArrayObject.D1 aod1 = (ArrayObject.D1) attr.getValues();
225-
226-
// getSize returns a long
227-
final int aodSize = (int) aod1.getSize();
228-
dimNames = new String[aodSize];
229-
230-
for (int i = 0; i < aodSize; ++i) {
231-
dimNames[i] = (String) aod1.get(i);
227+
if (attr.getLength() == 1 && attr.getStringValue().equals("")) {
228+
// scalar array without a named dimension
229+
logger.debug(" {} is a scalar array without a named dimension", vname);
230+
} else {
231+
final ArrayObject.D1 aod1 = (ArrayObject.D1) attr.getValues();
232+
233+
// getSize returns a long
234+
final int aodSize = (int) aod1.getSize();
235+
dimNames = new String[aodSize];
236+
237+
for (int i = 0; i < aodSize; ++i) {
238+
dimNames[i] = (String) aod1.get(i);
239+
}
240+
hasNamedDimensions = true;
232241
}
233-
hasNamedDimensions = true;
234242
} catch (final Exception exc) {
235-
logger.debug(" Could not extract _ARRAY_DIMENSIONS for {}, {}", vname, exc.getMessage());
243+
logger.debug(" Could not extract {} for {}, {}", ARRAYDIMENSIONS, vname, exc.getMessage());
236244
}
237245
}
238246
}
@@ -327,7 +335,13 @@ private List<Attribute> makeAttributes(RandomAccessDirectoryItem item) {
327335
Attribute.Builder attr = Attribute.builder(key);
328336
Object val = attrMap.get(key);
329337
if (val instanceof Collection<?>) {
330-
attr.setValues(Arrays.asList(((Collection) val).toArray()), false);
338+
Collection<?> collection = (Collection<?>) val;
339+
if (collection.isEmpty() && key.equals(ARRAYDIMENSIONS)) {
340+
// scalar array
341+
attr.setValues(Collections.singletonList(""), false);
342+
} else {
343+
attr.setValues(Arrays.asList(collection.toArray()), false);
344+
}
331345
} else if (val instanceof Number) {
332346
attr.setNumericValue((Number) val, false);
333347
} else {
@@ -354,7 +368,9 @@ private static int getChunkIndex(RandomAccessDirectoryItem item, ZArray zarray)
354368

355369
int nDims = zarray.getShape().length;
356370
// verify is data file, else return -1
357-
String pattern = String.format("([0-9]+%c){%d}[0-9]+", zarray.getSeparator().charAt(0), nDims - 1);
371+
String pattern = String.format("([0-9]+%c){%d}[0-9]+", zarray.getSeparator().charAt(0),
372+
nDims == 0 ? 0 : nDims - 1);
373+
358374
if (!fileName.matches(pattern)) {
359375
return -1;
360376
}
@@ -363,14 +379,20 @@ private static int getChunkIndex(RandomAccessDirectoryItem item, ZArray zarray)
363379
String[] dims = fileName.split(String.format("\\%c", zarray.getSeparator().charAt(0)));
364380
int[] subs = Arrays.stream(dims).mapToInt(dim -> Integer.parseInt(dim)).toArray();
365381

366-
// get number of chunks in each dimension
367-
int[] nChunks = new int[nDims];
368-
int[] shape = zarray.getShape();
369-
int[] chunkSize = zarray.getChunks();
370-
for (int i = 0; i < nDims; i++) {
371-
nChunks[i] = (int) Math.ceil(shape[i] / chunkSize[i]);
382+
// find chunk number as a flat index
383+
if (nDims != 0) {
384+
// get number of chunks in each dimension
385+
int[] nChunks = new int[nDims];
386+
int[] shape = zarray.getShape();
387+
int[] chunkSize = zarray.getChunks();
388+
for (int i = 0; i < nDims; i++) {
389+
nChunks[i] = (int) Math.ceil(shape[i] / chunkSize[i]);
390+
}
391+
return ZarrUtils.subscriptsToIndex(subs, nChunks);
392+
} else {
393+
// scalar array
394+
return 0;
372395
}
373-
return ZarrUtils.subscriptsToIndex(subs, nChunks);
374396
}
375397

376398
/**

‎cdm/zarr/src/main/java/ucar/nc2/iosp/zarr/ZarrLayoutBB.java‎

Lines changed: 18 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -1,3 +1,8 @@
1+
/*
2+
* Copyright (c) 2021-2025 University Corporation for Atmospheric Research/Unidata
3+
* See LICENSE for license information.
4+
*/
5+
16
package ucar.nc2.iosp.zarr;
27

38
import ucar.ma2.Range;
@@ -127,14 +132,20 @@ public LayoutBBTiled.DataChunk next() {
127132

128133
private void incrementChunk() {
129134
// increment index from inner dimension outward
130-
int i = this.currChunk.length - 1;
131-
while (this.currChunk[i] + 1 >= nChunks[i] && i > 0) {
132-
this.currChunk[i] = 0;
133-
i--;
135+
if (this.currChunk.length != 0) {
136+
int i = this.currChunk.length - 1;
137+
while (this.currChunk[i] + 1 >= nChunks[i] && i > 0) {
138+
this.currChunk[i] = 0;
139+
i--;
140+
}
141+
this.currChunk[i]++;
142+
this.currOffset += initializedChunks.getOrDefault(this.chunkNum, (long) 0);
143+
this.chunkNum = ZarrUtils.subscriptsToIndex(this.currChunk, nChunks);
144+
} else {
145+
// scalar array
146+
this.currOffset = 0;
147+
this.chunkNum = 0;
134148
}
135-
this.currChunk[i]++;
136-
this.currOffset += initializedChunks.getOrDefault(this.chunkNum, (long) 0);
137-
this.chunkNum = ZarrUtils.subscriptsToIndex(this.currChunk, nChunks);
138149
}
139150
}
140151

‎cdm/zarr/src/test/data/geozarr/README‎

Lines changed: 5 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,5 @@
1+
Examples obtained from https://github.com/zarr-developers/geozarr-spec/tree/cnl-examples/ on 2025-04-28 (commit d1eabecbea63d6f300600d529d7818bbf6df5bcc).
2+
No changes were made.
3+
4+
Examples are licensed under Creative Commons Attribution 4.0 International.
5+
Full text can be found at https://raw.githubusercontent.com/zarr-developers/geozarr-spec/refs/heads/cnl-examples/LICENSE

‎cdm/zarr/src/test/data/geozarr/xyt-raster.zarr/.zattrs‎

Lines changed: 7 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,7 @@
1+
{
2+
"Conventions": "CF-1.8",
3+
"profile": [
4+
"time-series-raster",
5+
"scalar-raster"
6+
]
7+
}

‎cdm/zarr/src/test/data/geozarr/xyt-raster.zarr/.zgroup‎

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,3 @@
1+
{
2+
"zarr_format": 2
3+
}

‎cdm/zarr/src/test/data/geozarr/xyt-raster.zarr/.zmetadata‎

Lines changed: 164 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,164 @@
1+
{
2+
"metadata": {
3+
".zattrs": {
4+
"Conventions": "CF-1.8",
5+
"profile": [
6+
"time-series-raster",
7+
"scalar-raster"
8+
]
9+
},
10+
".zgroup": {
11+
"zarr_format": 2
12+
},
13+
"spatial_ref/.zarray": {
14+
"chunks": [],
15+
"compressor": null,
16+
"dtype": "<i8",
17+
"fill_value": null,
18+
"filters": null,
19+
"order": "C",
20+
"shape": [],
21+
"zarr_format": 2
22+
},
23+
"spatial_ref/.zattrs": {
24+
"GeoTransform": "319281.23 12.5 0 5639331.75 0 -12.5",
25+
"_ARRAY_DIMENSIONS": [],
26+
"crs_wkt": "PROJCS[\"WGS 84 / UTM zone 30N\",GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS 84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4326\"]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"latitude_of_origin\",0],PARAMETER[\"central_meridian\",-3],PARAMETER[\"scale_factor\",0.9996],PARAMETER[\"false_easting\",500000],PARAMETER[\"false_northing\",0],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"Easting\",EAST],AXIS[\"Northing\",NORTH],AUTHORITY[\"EPSG\",\"32630\"]]",
27+
"false_easting": 500000.0,
28+
"false_northing": 0.0,
29+
"geographic_crs_name": "WGS 84",
30+
"grid_mapping_name": "transverse_mercator",
31+
"horizontal_datum_name": "World Geodetic System 1984",
32+
"inverse_flattening": 298.257223563,
33+
"latitude_of_projection_origin": 0.0,
34+
"longitude_of_central_meridian": -3.0,
35+
"prime_meridian_name": "Greenwich",
36+
"projected_crs_name": "WGS 84 / UTM zone 30N",
37+
"reference_ellipsoid_name": "WGS 84",
38+
"scale_factor_at_central_meridian": 0.9996,
39+
"semi_major_axis": 6378137.0,
40+
"semi_minor_axis": 6356752.314245179
41+
},
42+
"temperature/.zarray": {
43+
"chunks": [
44+
2,
45+
128,
46+
128
47+
],
48+
"compressor": {
49+
"blocksize": 0,
50+
"clevel": 5,
51+
"cname": "lz4",
52+
"id": "blosc",
53+
"shuffle": 1
54+
},
55+
"dtype": "<f8",
56+
"fill_value": "NaN",
57+
"filters": null,
58+
"order": "C",
59+
"shape": [
60+
4,
61+
128,
62+
128
63+
],
64+
"zarr_format": 2
65+
},
66+
"temperature/.zattrs": {
67+
"AREA_OR_POINT": "Point",
68+
"_ARRAY_DIMENSIONS": [
69+
"time",
70+
"y",
71+
"x"
72+
],
73+
"grid_mapping": "spatial_ref",
74+
"standard_name": "air_temperature",
75+
"units": "K"
76+
},
77+
"time/.zarray": {
78+
"chunks": [
79+
4
80+
],
81+
"compressor": {
82+
"blocksize": 0,
83+
"clevel": 5,
84+
"cname": "lz4",
85+
"id": "blosc",
86+
"shuffle": 1
87+
},
88+
"dtype": "<i8",
89+
"fill_value": null,
90+
"filters": null,
91+
"order": "C",
92+
"shape": [
93+
4
94+
],
95+
"zarr_format": 2
96+
},
97+
"time/.zattrs": {
98+
"_ARRAY_DIMENSIONS": [
99+
"time"
100+
],
101+
"calendar": "proleptic_gregorian",
102+
"long_name": "time of observation",
103+
"standard_name": "time",
104+
"units": "days since 2020-01-01 00:00:00"
105+
},
106+
"x/.zarray": {
107+
"chunks": [
108+
128
109+
],
110+
"compressor": {
111+
"blocksize": 0,
112+
"clevel": 5,
113+
"cname": "lz4",
114+
"id": "blosc",
115+
"shuffle": 1
116+
},
117+
"dtype": "<f8",
118+
"fill_value": "NaN",
119+
"filters": null,
120+
"order": "C",
121+
"shape": [
122+
128
123+
],
124+
"zarr_format": 2
125+
},
126+
"x/.zattrs": {
127+
"_ARRAY_DIMENSIONS": [
128+
"x"
129+
],
130+
"long_name": "x coordinate of projection",
131+
"standard_name": "projection_x_coordinate",
132+
"units": "metre"
133+
},
134+
"y/.zarray": {
135+
"chunks": [
136+
128
137+
],
138+
"compressor": {
139+
"blocksize": 0,
140+
"clevel": 5,
141+
"cname": "lz4",
142+
"id": "blosc",
143+
"shuffle": 1
144+
},
145+
"dtype": "<f8",
146+
"fill_value": "NaN",
147+
"filters": null,
148+
"order": "C",
149+
"shape": [
150+
128
151+
],
152+
"zarr_format": 2
153+
},
154+
"y/.zattrs": {
155+
"_ARRAY_DIMENSIONS": [
156+
"y"
157+
],
158+
"long_name": "y coordinate of projection",
159+
"standard_name": "projection_y_coordinate",
160+
"units": "metre"
161+
}
162+
},
163+
"zarr_consolidated_format": 1
164+
}

‎cdm/zarr/src/test/data/geozarr/xyt-raster.zarr/spatial_ref/.zarray‎

Lines changed: 10 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,10 @@
1+
{
2+
"chunks": [],
3+
"compressor": null,
4+
"dtype": "<i8",
5+
"fill_value": null,
6+
"filters": null,
7+
"order": "C",
8+
"shape": [],
9+
"zarr_format": 2
10+
}

‎cdm/zarr/src/test/data/geozarr/xyt-raster.zarr/spatial_ref/.zattrs‎

Lines changed: 19 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,19 @@
1+
{
2+
"GeoTransform": "319281.23 12.5 0 5639331.75 0 -12.5",
3+
"_ARRAY_DIMENSIONS": [],
4+
"crs_wkt": "PROJCS[\"WGS 84 / UTM zone 30N\",GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS 84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4326\"]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"latitude_of_origin\",0],PARAMETER[\"central_meridian\",-3],PARAMETER[\"scale_factor\",0.9996],PARAMETER[\"false_easting\",500000],PARAMETER[\"false_northing\",0],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"Easting\",EAST],AXIS[\"Northing\",NORTH],AUTHORITY[\"EPSG\",\"32630\"]]",
5+
"false_easting": 500000.0,
6+
"false_northing": 0.0,
7+
"geographic_crs_name": "WGS 84",
8+
"grid_mapping_name": "transverse_mercator",
9+
"horizontal_datum_name": "World Geodetic System 1984",
10+
"inverse_flattening": 298.257223563,
11+
"latitude_of_projection_origin": 0.0,
12+
"longitude_of_central_meridian": -3.0,
13+
"prime_meridian_name": "Greenwich",
14+
"projected_crs_name": "WGS 84 / UTM zone 30N",
15+
"reference_ellipsoid_name": "WGS 84",
16+
"scale_factor_at_central_meridian": 0.9996,
17+
"semi_major_axis": 6378137.0,
18+
"semi_minor_axis": 6356752.314245179
19+
}

‎cdm/zarr/src/test/data/geozarr/xyt-raster.zarr/spatial_ref/0‎

8 Bytes
Binary file not shown.

0 commit comments

Comments
 (0)