Mark grid_mapping variable (crs) as a coordinate

This is related to #837
opendatacube · Apr 6, 2020 · a140ad7 · a140ad7
1 parent a31b222
commit a140ad7
Showing 1 changed file with 8 additions and 0 deletions.
diff --git a/datacube/drivers/netcdf/writer.py b/datacube/drivers/netcdf/writer.py
@@ -237,6 +237,14 @@ def create_grid_mapping_variable(nco, crs, name=DEFAULT_GRID_MAPPING):
         crs_var = _create_projected_grid_mapping_variable(nco, crs, name)
     else:
         raise ValueError('Unknown CRS')
+
+    # mark crs variable as a coordinate
+    coords = getattr(nco, 'coordinates', None)
+    coords = [] if coords is None else coords.split(',')
+    if name not in coords:
+        coords.append(name)
+    nco.coordinates = ','.join(coords)
+
     crs_var.semi_major_axis = crs.semi_major_axis
     crs_var.semi_minor_axis = crs.semi_minor_axis
     crs_var.inverse_flattening = crs.inverse_flattening