diff --git a/.DS_Store b/.DS_Store new file mode 100644 index 00000000..30cea4d6 Binary files /dev/null and b/.DS_Store differ diff --git a/0.1/index.bs b/0.1/index.bs index b60204f1..52a33115 100644 --- a/0.1/index.bs +++ b/0.1/index.bs @@ -501,8 +501,7 @@ Projects which support reading and/or writing OME-NGFF data include: -Diagram of related projects +Diagram of related projects All implementations prevent an equivalent representation of a dataset which can be downloaded or uploaded freely. An interactive version of this diagram is available from the [OME2020 Workshop](https://downloads.openmicroscopy.org/presentations/2020/Dundee/Workshops/NGFF/zarr_diagram/). diff --git a/0.2/index.bs b/0.2/index.bs index 4e25ad2c..487f91a3 100644 --- a/0.2/index.bs +++ b/0.2/index.bs @@ -559,8 +559,7 @@ Projects which support reading and/or writing OME-NGFF data include: -Diagram of related projects +Diagram of related projects All implementations prevent an equivalent representation of a dataset which can be downloaded or uploaded freely. An interactive version of this diagram is available from the [OME2020 Workshop](https://downloads.openmicroscopy.org/presentations/2020/Dundee/Workshops/NGFF/zarr_diagram/). diff --git a/0.3/index.bs b/0.3/index.bs index 84192b65..47bc2dcf 100644 --- a/0.3/index.bs +++ b/0.3/index.bs @@ -569,8 +569,7 @@ Projects which support reading and/or writing OME-NGFF data include: -Diagram of related projects +Diagram of related projects All implementations prevent an equivalent representation of a dataset which can be downloaded or uploaded freely. An interactive version of this diagram is available from the [OME2020 Workshop](https://downloads.openmicroscopy.org/presentations/2020/Dundee/Workshops/NGFF/zarr_diagram/). diff --git a/0.4/index.bs b/0.4/index.bs index a0c3f83d..54786e9f 100644 --- a/0.4/index.bs +++ b/0.4/index.bs @@ -578,8 +578,7 @@ Projects which support reading and/or writing OME-NGFF data include: -Diagram of related projects +Diagram of related projects All implementations prevent an equivalent representation of a dataset which can be downloaded or uploaded freely. An interactive version of this diagram is available from the [OME2020 Workshop](https://downloads.openmicroscopy.org/presentations/2020/Dundee/Workshops/NGFF/zarr_diagram/). diff --git a/index.html b/index.html index 74f803c4..e69de29b 100644 --- a/index.html +++ b/index.html @@ -1,6 +0,0 @@ - - - - -

If you are not redirected in five seconds, click here.

- diff --git a/latest/generate_anndata_example.py b/latest/generate_anndata_example.py new file mode 100644 index 00000000..9b5b9be8 --- /dev/null +++ b/latest/generate_anndata_example.py @@ -0,0 +1,149 @@ +import numpy as np +import pandas as pd +import anndata as ad +import zarr +import numcodecs +from scipy.sparse import csr_matrix, csc_matrix + +def generate_example(n_obs, n_var): + # X and layers + adata = ad.AnnData(csr_matrix(np.random.poisson(1, size=(n_obs, n_var)), dtype=np.float32)) + adata.layers["log_transformed"] = np.log1p(adata.X) + adata.layers["other_data"] = np.random.poisson(1, size=(n_obs, n_var)) + 1.0 + + # obs and var names + adata.obs_names = [f"Cell_{i:d}" for i in range(adata.n_obs)] + adata.var_names = [f"Gene_{i:d}" for i in range(adata.n_vars)] + + # annotations + # the tutorial mentions that string arrays are automatically converted to categoricals if convenient + adata.obs["cell_type"] = pd.Categorical(np.random.choice(["B", "T", "Monocyte"], size=(adata.n_obs,))) + adata.obsm["X_umap"] = np.random.normal(0, 1, size=(adata.n_obs, 2)) + adata.varm["gene_stuff"] = np.random.normal(0, 1, size=(adata.n_vars, 5)) + adata.obsp["pairwise_data"] = csc_matrix(np.random.poisson(1, size=(n_obs, n_obs)), dtype=np.int32) + + # uns + adata.uns["random"] = [1, 2, 3] + + return adata + + +def write_anndata(adata, filename, chunks): + # root and tables + root = zarr.open(filename, mode="w") + root.array("some_image", np.array([0]), chunks=(1,)) + tables = root.create_group("tables") + tables.attrs["tables"] = ["/anndata/obs", "/anndata/var", "/anndata/obsm", "/anndata/varm", "/anndata/obsp", "/anndata/varp"] + adgroup = tables.create_group("anndata") + adgroup.attrs["anndata"] = "0.9.1" + adgroup.attrs["other-metadata"] = "metadata describing how the anndata annotates some_image" + + # X/layers + X = adgroup.array('X', adata.X.todense(), chunks=chunks) + layers = adgroup.create_group("layers") + layers.array("log_transformed", adata.layers["log_transformed"].todense(), chunks=chunks) + layers.array("other_data", adata.layers["other_data"], chunks=chunks) + + # obs + localChunks = (chunks[0],) + obs = adgroup.create_group("obs") + obs.create_dataset("row_names", data=np.array(adata.obs_names), dtype=object, object_codec=numcodecs.VLenUTF8()) + obs.create_dataset("cell_type", data=np.array(adata.obs["cell_type"]), chunks=localChunks, object_codec=numcodecs.VLenUTF8()) + obs.attrs["annotated-data"] = get_annotated_data_map(dimension=0) + obs.attrs["column-order"] = ["row_names", "cell_type"] + obs.attrs["row-names"] = "row_names" + + # obsm + obsm = adgroup.create_group("obsm") + obsm.create_dataset("X_umap", data=adata.obsm["X_umap"], chunks=(chunks[0], 2)) + obsm.attrs["annotated-data"] = get_annotated_data_map(dimension=0) + obsm.attrs["column-order"] = ["X_umap"] + + # obsp + obsp = adgroup.create_group("obsp") + obsp.create_dataset("pairwise_data", data=adata.obsp["pairwise_data"].todense(), chunks=(chunks[0], chunks[0])) + obsp.attrs["annotated-data"] = get_annotated_data_map(dimension=0) + obsp.attrs["column-order"] = ["pairwise_data"] + + # var + var = adgroup.create_group("var") + var.create_dataset("row_names", data=np.array(adata.var_names), chunks=(chunks[1],), object_codec=numcodecs.VLenUTF8()) + var.attrs["annotated-data"] = get_annotated_data_map(dimension=1) + var.attrs["column-order"] = ["row_names"] + var.attrs["row-names"] = "row_names" + + # varm + varm = adgroup.create_group("varm") + varm.create_dataset("gene_stuff", data=adata.varm["gene_stuff"], chunks=(chunks[1], 5)) + varm.attrs["annotated-data"] = get_annotated_data_map(dimension=1) + varm.attrs["column-order"] = ["gene_stuff"] + + # varm + varp = adgroup.create_group("varp") + varp.attrs["annotated-data"] = get_annotated_data_map(dimension=1) + varp.attrs["column-order"] = [] + + # uns + uns = adgroup.create_group("uns") + uns.create_dataset("random", data=adata.uns["random"], chunks=(3,)) + + +def get_annotated_data_map(*, dimension): + return [{"array": "/tables/anndata/X", "dimension": str(dimension)}, + {"array": "/tables/anndata/layers/log_transformed", "dimension": str(dimension)}, + {"array": "/tables/anndata/layers/other_data", "dimension": str(dimension)}] + + +def write_anndata_suggestion(adata, filename, chunks): + # root and tables + root = zarr.open(filename, mode="w") + root.array("some_image", np.array([0]), chunks=(1,)) + tables = root.create_group("tables") + tables.attrs["tables"] = ["/anndata/obs", "/anndata/var"] + adgroup = tables.create_group("anndata") + adgroup.attrs["anndata"] = "0.9.1" + adgroup.attrs["other-metadata"] = "metadata describing how the anndata annotates some_image" + + # X and layers are combined into one array, a table is used to name the layers + all_together = np.stack([np.array(adata.X.todense()), np.array(adata.layers["log_transformed"].todense()), adata.layers["other_data"]], axis=2) + X = adgroup.array('X', all_together, chunks=(*chunks,1)) + layers = adgroup.create_group("layers") + row_names = np.array(["X", "log_transformed", "other_data"]) + layers.create_dataset("row_names", data=row_names, dtype=object, object_codec=numcodecs.VLenUTF8()) + layers.attrs["annotated-data"] = [{"array": "/tables/anndata/X", "dimension": "2"}] + layers.attrs["row-names"] = "row_names" + + # obs (combines obs, obsm, obsp) + localChunks = (chunks[0],) + obs = adgroup.create_group("obs") + obs.create_dataset("row_names", data=np.array(adata.obs_names), dtype=object, object_codec=numcodecs.VLenUTF8()) + obs.create_dataset("cell_type", data=np.array(adata.obs["cell_type"]), chunks=localChunks, object_codec=numcodecs.VLenUTF8()) + obs.create_dataset("X_umap", data=adata.obsm["X_umap"], chunks=(chunks[0], 2)) + obs.create_dataset("pairwise_data", data=adata.obsp["pairwise_data"].todense(), chunks=(chunks[0], chunks[0])) + obs.attrs["annotated-data"] = [{"array": "/tables/anndata/X", "dimension": "0"}] + obs.attrs["column-order"] = ["row_names", "cell_type", "X_umap", "pairwise_data"] + obs.attrs["row-names"] = "row_names" + + # var (combines var, varm, varp) + var = adgroup.create_group("var") + var.create_dataset("row_names", data=np.array(adata.var_names), chunks=(chunks[1],), object_codec=numcodecs.VLenUTF8()) + var.create_dataset("gene_stuff", data=adata.varm["gene_stuff"], chunks=(chunks[1], 5)) + var.attrs["annotated-data"] = [{"array": "/tables/anndata/X", "dimension": "1"}] + var.attrs["column-order"] = ["row_names", "gene_stuff"] + var.attrs["row-names"] = "row_names" + + # uns + uns = adgroup.create_group("uns") + uns.create_dataset("random", data=adata.uns["random"], chunks=(3,)) + + +# generate example and store as zarr using the minimal table spec proposal +n_obs = 10 +n_var = 200 +chunks = (10, 40) + +adata = generate_example(n_obs, n_var) +write_anndata(adata, "example.zarr", chunks) + +# store example in an alternative way, exploiting the properties of the suggested minimal table spec a bit more +write_anndata_suggestion(adata, "example_suggestion.zarr", chunks) diff --git a/latest/index.bs b/latest/index.bs index ed8b6aa4..996a01bf 100644 --- a/latest/index.bs +++ b/latest/index.bs @@ -177,308 +177,112 @@ For this example we assume an image with 5 dimensions and axes called `t,c,z,y,x Tables {#table-layout} ---------------------- -The following describes the expected layout for tabular data. -OME-NGFF tables are compatible with the [AnnData model](https://github.com/scverse/anndata). +A table is a 2-dimensional data structure consisting of rows and columns. +Tables are an intuitive way of organizing data or metadata consisting of variables +and records, which are often organized as columns and rows, respectively. + +In OME-NGFF, a table is a Zarr group containing zero, one, or more Zarr arrays, where each +array represents one column of the table. Columns are ordered, and each column in a +table MUST have the same number of rows. While the table itself MUST be 2-dimensional, +the columns need not be. 1-dimensional columns will be typical. Nevertheless, there +may be use cases in which it is sensible to conceptualize a 2D, 3D, or higher-dimensional +array as a single column. A row of the table is the set of arrays obtained by slicing all +columns at the same index in the 0-th dimension. + +Tables are located in the `tables/` directory, in the root of an image's +Zarr group, alongside the `labels/` directory if one is present. The .zattrs file +of `tables` MUST contain the "tables" property, which indicates the names of the +tables as a JSON array. + +The .zattrs file in each table of `tables/` MAY contain the property "row-names", +and it MUST contain the two properties "column-order", and "annotated-data". + +"row-names" is a string indicating which column, if any, contains the table's row names. +An array of row names is a special column containing no duplicate or missing values. +For example, in a table in which the rows represent genes, it may be sensible to +designate the column called "geneID" as the "row-names" column. The "row-names" column +may be either an array of strings or an array of integers. The "row-names" column need not appear +first in the ordered array of column names. Either a missing "row-names" property or an empty +string value (i.e. "row-names": "") indicates that the table does not have row names. + +"column-order" is a JSON array containing all of the column names, in order. Column names listed +in "column-order" MUST be strings, even if they are string representations of a digit, e.g. "1". +If the table contains no columns, then "column-order" MUST be an empty array. + +In some cases, the table provides information about--in other words, annotates--a +particular dimension of an array. The "annotated-data" property describes these cases. +Tables MAY annotate arrays, and can only annotate arrays. If a table is annotating an +array, then the 0-th dimension of the column(s) (i.e. the number of rows in the table) +MUST correspond to the length of the dimension that is being annotated in the source +array. The required "annotated-data" property is a JSON array of Zarr array paths and +dimensions (0-based indexing), as shown below: + +```json +"annotated-data": [ // A JSON array containing all array / dimension pairs this table is annotating + { + "array": "/path/to/arrayA", // This table annotates the 0-th dimension of arrayA ... + "dimension": 0, + }, + { + "array": "/path/to/arrayA", // ... and the first dimension of arrayA ... + "dimension": 1, + }, + { + "array": "/path/to/arrayB", // ... and the first dimension of arrayB. + "dimension": 1, + } +], +``` + +The "annotated-data" property MUST be present, but it MAY be an empty JSON array. If +`annotated-data` contains exactly one entry, then the 0-th dimension of all arrays +(columns) in the table group SHOULD be chunked in the same way as the corresponding +dimension of the source array. + +Note: The AnnData data model is based on a main array with additional "annotation" +tables that have the same number of rows, columns, or both, as the main array. AnnData +objects can be stored in the `tables/` directory of an OME-NGFF file. See the AnnData +documentation for detailed recommendations on formatting AnnData data structures within Zarr. + +There MAY be one or more intermediate directories between `tables/` and a particular table. +These SHOULD NOT contain metadata, unless the intermediate directory represents the parent +directory for an AnnData object. The names of directories beneath the `tables/` +directory are arbitrary, except in the AnnData case. + + +The on-disk format of a table group looks like this:
-.                             # Root folder, potentially in S3,
-│                             # with a flat list of images.
-│
-└── 123.zarr
+    .                            # Root folder, potentially in S3.
     |
-    ├── .zgroup
-    |
-    ├── .zattrs
-    |
-    └── tables                # The tables group is a container which holds one or multiple tables that are compatible with AnnData.
+    └── 123.zarr
         |
-        │                     # The tables group MAY be in the root of the zarr file.
-        ├── .zgroup           # The tables group MAY be in root or in another group.
+        ├── .zgroup
         |
-        ├── .zattrs           # `.zattrs` MUST contain "tables", which lists the keys of the subgroups that are tables. In this case, the only table is "my_table".
-                              # hence `.zattrs` should be equal to `{ "tables": [ "my_table" ] }`.
-        |
-        └── my_table
-            │                     # The table group MAY be in the root of the zarr file.
-            ├── .zgroup           # The table group MAY be in root or in another group.
-            |
-            ├── .zattrs           # `.zattrs` MUST contain "type", which is set to `"ngff:region_table"`
-            |                     # `.zattrs` MUST contain "region", which is the path to the data the table is annotating.
-            |                     # "region" MUST be a single path (single region) or an array of paths (multiple regions).
-            |                     # "region" paths MUST be objects with a key "path" and the path value MUST be a string.
-            |                     # `.zattrs` MUST contain "region_key" if "region" is an array. "region_key" is the key in `obs` denoting which region a given row corresponds to.
-            |                     # `.zattrs` MAY contain "instance_key", which is the key in `obs` that denotes which instance in "region" the row corresponds to. If "instance_key" is not provided, the values from the `obs` `.zattrs` "_index" key is used.
-            │
-            ├── X                 # You MAY add an zarr array `X`.
-            │   │                 # `X` MUST not be a complex type (i.e., MUST be a single type)
-            │   │                 # `X` MAY be chunked as the user desires.
-            │   ├── .zarray
-            │   ├── 0.0
-            │   │   ...
-            │   └── n.m
-            |
-            ├── layers            # You MAY add a `layers` group, which contains dense matrices with the same shape as X.
-            │   │
-            │   ├── .zgroup
-            │   ├── .zattrs       # `.zattrs` MUST contain `"keys"`, which is an array of the names of the subgroups containing a `layer`.
-            │   │
-            │   └── layer_0       # You MAY add a zarr array for each layer
-            |       |             # Each layer array MUST have the same shape as X
-            |       |             # Each layer array SHOULD be chunked the same as X
-            |       ├── .zarray
-            |       |
-            |       ├── 0.0
-            │       │   ...
-            │       └── n.m
-            │        
-            ├── obs               # You MUST add an obs group container. The obs group holds a table of annotations on the rows in X.
-            │   │                 # The rows in obs MUST be index-matched to the rows in X.
-            │   ├── .zgroup
-            │   │                     
-            │   ├── .zattrs       # `.zattrs` MUST contain `"_index"`, which is the name of the column in obs to be used as the index.           
-            │   │                 # `.zattrs` MUST contain `"column-order"`, which is a list of the order of the non-_index columns.
-            │   │                 # `.zattrs` MUST contain `"encoding-type"`, which is set to `"dataframe"` by AnnData.
-            │   │                 # `.zattrs` MUST contain `"encoding-version"`, which is set to `"0.2.0"` by AnnData.
-            │   │      
-            │   └── col_0         # Each column in the obs table is a 1D zarr array. The rows can be chunked as the user desires.
-            │       ├── .zarray   # However, the obs columns SHOULD be chunked in the same way as the rows in X (if present).
-            │       │
-            │       └─ 0
-            ├── var               # You MAY add a var group container. The var group holds a table of annotations on the columns in X.
-            |   │                 # The rows in var MUST be index-matched to the columns in X (if present). 
-            |   |
-            |   ├── .zattrs       # `.zattrs` MUST contain `"_index"`, which is the name of the column in obs to be used as the index.           
-            |   │                 # `.zattrs` MUST contain `"column-order"`, which is a list of the order of the non-_index columns.
-            |   │                 # `.zattrs` MUST contain `"encoding-type"`, which is set to `"dataframe"` by AnnData.
-            |   │                 # `.zattrs` MUST contain `"encoding-version"`, which is set to `"0.2.0"` by AnnData.
-            |   │
-            |   ├── array_col     # Columns in the var table MAY be a 1D zarr array. The rows can be chunked as the user desires.
-            |   |   ├── .zarray   # However, the var columns SHOULD be chunked in the same way as the columns in X.
-            |   |   │
-            |   |   └─ 0
-            |   |
-            |   └── cat_col       # Columns in the var table MAY be categorical
-            |       ├── .zattrs.  # `.zattrs` MUST contain `"encoding-type"`, which is set to `"categorical"` by AnnData.
-            |       |             # `.zattrs` MUST contain `"encoding-version"`, which is set to `"0.2.0"` by AnnData.
-            |       |
-            |       ├── categories
-            |       |  ├── .zarray # categories MUST be a 1D zarr array. The rows can be chunked as the user desires.
-            |       |  |
-            |       |  └─ 0
-            |       ├── codes
-            |       |  ├── .zarray # codes MUST be a 1D zarr array. The rows can be chunked as the user desires.
-            |       |  |
-            |       |  └─ 0
-            |       |
-            |       ├── null_col  # Columns in the var table MAY nullable integer
-            |       ├── .zattrs.  # `.zattrs` MUST contain `"encoding-type"`, which is set to `"nullable-integer"` by AnnData.
-            |       |             # `.zattrs` MUST contain `"encoding-version"`, which is set to `"0.1.0"` by AnnData.
-            |       |
-            |       ├── mask
-            |       |  ├── .zarray # categories MUST be a 1D zarr array. The rows can be chunked as the user desires.
-            |       |  |
-            |       |  └─ 0
-            |       └── values
-            |          ├── .zarray # codes MUST be a 1D zarr array. The rows can be chunked as the user desires.
-            |          |
-            |          └─ 0
-            |
-            ├── obsm              # You MAY add a obsm group comtainer. The obsm group contains arrays that annotate the rows in X.
-            |   │                 # The rows in each array MUST be index-matched to the rows in X (if present). 
-            |   |
-            │   ├── .zgroup
-            |   |
-            |   ├── .zattrs       # `.zattrs` MUST contain `"encoding-type"`, which is set to `"dict"` by AnnData.           
-            |   │                 # `.zattrs` MUST contain `"encoding-version"`, which is set to `"0.1.0"` by AnnData.
-            |   |                 # `.zattrs` MUST contain `"keys"`, which is an array of the names of the subgroups containing `obsm` arrays.
-            |   │
-            │   └── obsm_0        # You MAY add a zarr array for each obsm matrix.
-            |       |             # Each obsm array MUST have the same number of rows as X.
-            |       |             # The rows in each obsm array SHOULD be chunked the same as the rows in X.
-            |       ├── .zarray
-            |       |
-            |       ├── 0.0
-            │       │   ...
-            │       └── n.m
-            |
-            ├── varm              # You MAY add a varm group comtainer. The varm group contains arrays that annotate the columns in X.
-            |   │                 # The rows in each array MUST be index-matched to the columns in X (if present). 
-            |   |
-            │   ├── .zgroup
-            |   |
-            |   ├── .zattrs       # `.zattrs` MUST contain `"encoding-type"`, which is set to `"dict"` by AnnData.           
-            |   │                 # `.zattrs` MUST contain `"encoding-version"`, which is set to `"0.1.0"` by AnnData.
-            |   |                 # `.zattrs` MUST contain `"keys"`, which is an array of the names of the subgroups containing `varm` arrays.
-            |   │
-            │   └── varm_0        # You MAY add a zarr array for each varm matrix.
-            |       |             # Each varm array MUST have the same number of rows as columns in X.
-            |       |             # The rows in each obsm array SHOULD be chunked the same as the columns in X.
-            |       ├── .zarray
-            |       ├── 0.0
-            │       │   ...
-            │       └── n.m
+        └──tables
             |
-            ├── obsp              # You MAY add a obsp group comtainer. The obsp group contains sparse arrays that annotate the rows in X.
-            |   │                 # The rows in each array MUST be index-matched to the columns in X (if present). 
-            |   |
-            │   ├── .zgroup
-            |   |
-            |   ├── .zattrs       # `.zattrs` MUST contain `"encoding-type"`, which is set to `"dict"` by AnnData.           
-            |   │                 # `.zattrs` MUST contain `"encoding-version"`, which is set to `"0.1.0"` by AnnData.
-            |   |                 # `.zattrs` MUST contain `"keys"`, which is an array of the names of the subgroups containing `obsp` arrays.
-            |   │
-            │   └── obsp_0        # You MAY add a zarr group for each obsp array.
-            |       |             # Each obsp array MUST have the same number of rows as rows in X.
-            |       |
-            │       ├── .zgroup
-            |       |
-            |       ├── .zattrs   # `.zattrs` MUST contain `"encoding-type"`, which is set to `"csr_matrix"` or `"csc_matrix"` for compressed sparse row and compressed sparse column, respectively.          
-            |       │             # `.zattrs` MUST contain `"encoding-version"`, which is set to `"0.1.0"` by AnnData.
-            |       |             # `.zattrs` MUST contain `"shape"` which is an array giving the shape of the densified array.
-            |       |
-            |       ├── data      # You MUST add a one-dimensional zarr array named "data". 
-            |       |   |         # `data` MAY be chunked as the user desires.
-            |       |   ├── .zarray
-            |       |   |
-            |       |   ├── 0
-            │       │   |   ...
-            │       |   └── n
-            |       |
-            |       ├── indices   # You MUST add a one-dimensional zarr array named "indices".
-            |       |   |         # `indices` MAY be chunked as the user desires.
-            |       |   ├── .zarray  # `indices` MUST be an `int` dtype.
-            |       |   |
-            |       |   ├── 0
-            │       │   |   ...
-            │       |   └── n
-            |       |
-            |       └── indptr    # You MUST add a one-dimensional zarr array named "indptr".
-            |           |         # `indptr` MAY be chunked as the user desires.
-            |           ├── .zarray  # `indptr` MUST be an `int` dtype.
-            |           |
-            |           ├── 0
-            │           |   ...
-            │           └── n
+            ├── .zgroup
+            ├── .zattrs          # A table group's `.zattrs` file MUST contain an attribute "tables" that lists all
+            |                    # tables in the group with relative paths, e.g., `{"tables": ["table_1", "some/nested/table_2"]}` 
             |
-            ├── varp              # You MAY add a varp group comtainer. The varp group contains sparse arrays that annotate the columns in X.
-            |   │                 # The rows in each array MUST be index-matched to the columns in X (if present). 
-            |   |
-            │   ├── .zgroup
-            |   |
-            |   ├── .zattrs       # `.zattrs` MUST contain `"encoding-type"`, which is set to `"dict"` by AnnData.           
-            |   │                 # `.zattrs` MUST contain `"encoding-version"`, which is set to `"0.1.0"` by AnnData.
-            |   |                 # `.zattrs` MUST contain `"keys"`, which is an array of the names of the subgroups containing `varp` arrays.
-            |   │
-            │   └── varp_0        # You MAY add a zarr group for each varp array.
-            |       |             # Each varp array MUST have the same number of rows as columns in X.
-            |       |
-            │       ├── .zgroup
-            |       |
-            |       ├── .zattrs   # `.zattrs` MUST contain `"encoding-type"`, which is set to `"csr_matrix"` or `"csc_matrix"` for compressed sparse row and compressed sparse column, respectively.           
-            |       │             # `.zattrs` MUST contain `"encoding-version"`, which is set to `"0.1.0"` by AnnData.
-            |       |             # `.zattrs` MUST contain `"shape"` which is an array giving the shape of the densified array.
-            |       |
-            |       ├── data      # You MUST add a one-dimensional zarr array named "data". 
-            |       |   |         # `data` MAY be chunked as the user desires.
-            |       |   ├── .zarray
-            |       |   |
-            |       |   ├── 0
-            │       │   |   ...
-            │       |   └── n
-            |       |
-            |       ├── indices   # You MUST add a one-dimensional zarr array named "indices".
-            |       |   |         # `indices` MAY be chunked as the user desires.
-            |       |   ├── .zarray  # `indices` MUST be an `int` dtype.
-            |       |   |
-            |       |   ├── 0
-            │       │   |   ...
-            │       |   └── n
-            |       |
-            |       └── indptr    # You MUST add a one-dimensional zarr array named "indptr".
-            |           |         # `indptr` MAY be chunked as the user desires.
-            |           ├── .zarray  # `indptr` MUST be an `int` dtype.
-            |           |
-            |           ├── 0
-            │           |   ...
-            │           └── n
-            |
-            └── uns               # You MAY add a uns containter to store unstructured data.
-                |
+            └── table_1          # A table group holds arrays (called columns) with one or more dimensions, 
+                |                # where the number of rows n MUST be the same across all columns.
                 ├── .zgroup
-                |
-                ├── .zattrs       # `.zattrs` MUST contain `"encoding-type"`, which is set to `"dict"` by AnnData.           
-                │                 # `.zattrs` MUST contain `"encoding-version"`, which is set to `"0.1.0"` by AnnData.
-                │
-                ├── group         # You MAY add zarr groups.
-                |   |             # `uns` groups MAY contain groups, dataframes, dense arrays, and sparse arrays.
-                |   |
-                |   ├── .zgroup
-                |   |
-                |   ├── .zattrs   # `.zattrs` MUST contain `"encoding-type"`, which is set to `"csr_matrix"` by AnnData.           
-                |   │             # `.zattrs` MUST contain `"encoding-version"`, which is set to `"0.1.0"` by AnnData.
-                |   ...
-                |
-                ├── dataframe_0   # You MAY add dataframe group containers.
-                |   |                 # dataframes MAY be in the `uns` group or in a subgroup.
-                |   │
-                |   ├── .zgroup
-                |   │                     
-                |   ├── .zattrs       # `.zattrs` MUST contain `"_index"`, which is the name of the column in obs to be used as the index.           
-                |   │                 # `.zattrs` MUST contain `"column-order"`, which is a list of the order of the non-_index columns.
-                |   │                 # `.zattrs` MUST contain `"encoding-type"`, which is set to `"dataframe"` by AnnData.
-                |   │                 # `.zattrs` MUST contain `"encoding-version"`, which is set to `"0.2.0"` by AnnData.
-                |   │      
-                |   └── col_0         # Each column in the obs table is a 1D zarr array.
-                |       ├── .zarray   # Each columns MUST be chunked the same, but the chunking may be chosen by the user.
-                |       │
-                |       └─ 0
-                |
-                ├── dense_array       # You MAY dense arrays as n n-dimensional zarr arrays.
-                |   │                 # `dense_array` MUST not be a complex type (i.e., MUST be a single type)
-                |   │                 # `dense_array` MAY be chunked as the user desires.
-                |   |                 # `dense array` MAY be in the `uns` group or in a subgroup.
-                |   |
-                |   ├── .zarray
-                |   ├── 0.0
-                |   │   ...
-                |   └── n.m
-                |
-                └── sparse_array  # You MAY add sparse arrays as a zarr group for each sparse array.
-                    |             # sparse arrays MAY be in the `uns` group or in a subgroup.
-                    |
-                    ├── .zgroup
-                    |
-                    ├── .zattrs      # `.zattrs` MUST contain `"encoding-type"`, which is set to `"csr_matrix"` or `"csc_matrix"` for compressed sparse row and compressed sparse column, respectively.           
-                    │                # `.zattrs` MUST contain `"encoding-version"`, which is set to `"0.1.0"` by AnnData.
-                    |                # `.zattrs` MUST contain `"shape"` which is an array giving the shape of the densified array.
+                ├── .zattrs      # `.zattrs` MAY contain "row-names", which states which column contains row names.
+                |                # `.zattrs` MUST contain "column-order", which is a list of the order of the columns.
+                |                # `.zattrs` MUST contain "annotated-data", which is specified above.
+                |                
+                |     
+                └── column_1     # The table group MAY hold zero, one, or an arbitrary number of columns.
                     |
-                    ├── data         # You MUST add a one-dimensional zarr array named "data". 
-                    |   |            # `data` MAY be chunked as the user desires.
-                    |   ├── .zarray
-                    |   |
-                    |   ├── 0
-                    │   |   ...
-                    |   └── n
+                    ├── .zarray
                     |
-                    ├── indices      # You MUST add a one-dimensional zarr array named "indices".
-                    |   |            # `indices` MAY be chunked as the user desires.
-                    |   ├── .zarray  # `indices` MUST be an `int` dtype.
-                    |   |
-                    |   ├── 0
-                    │   |   ...
-                    |   └── n
-                    |
-                    └── indptr       # You MUST add a one-dimensional zarr array named "indptr".
-                        |            # `indptr` MAY be chunked as the user desires.
-                        ├── .zarray  # `indptr` MUST be an `int` dtype.
-                        |
-                        ├── 0
-                        |   ...
-                        └── n
-        
-
+                    ├── 0.0
+                    |   ...
+                    └── n.m
 
+ High-content screening {#hcs-layout} ------------------------------------ @@ -882,8 +686,7 @@ Projects which support reading and/or writing OME-NGFF data include: -Diagram of related projects +Diagram of related projects All implementations prevent an equivalent representation of a dataset which can be downloaded or uploaded freely. An interactive version of this diagram is available from the [OME2020 Workshop](https://downloads.openmicroscopy.org/presentations/2020/Dundee/Workshops/NGFF/zarr_diagram/). @@ -896,7 +699,7 @@ Citing {#citing} [Next-generation file format (NGFF) specifications for storing bioimaging data in the cloud.](https://ngff.openmicroscopy.org/0.4) J. Moore, *et al*. Editors. Open Microscopy Environment Consortium, 8 February 2022. -This edition of the specification is [https://ngff.openmicroscopy.org/0.4/](https://ngff.openmicroscopy.org/0.4/]). +This edition of the specification is [https://ngff.openmicroscopy.org/0.4/](https://ngff.openmicroscopy.org/0.4/). The latest edition is available at [https://ngff.openmicroscopy.org/latest/](https://ngff.openmicroscopy.org/latest/). [(doi:10.5281/zenodo.4282107)](https://doi.org/10.5281/zenodo.4282107) diff --git a/latest/index.html b/latest/index.html new file mode 100644 index 00000000..5eea766b --- /dev/null +++ b/latest/index.html @@ -0,0 +1,2992 @@ + + + + + Next-generation file formats (NGFF) + + + + + + + + + + + + + +
+ OME logo (6 circles in a hexagon) +

Next-generation file formats (NGFF)

+

Editor’s Draft,

+
+
+
This version: +
https://ngff.openmicroscopy.org/latest/ +
Issue Tracking: +
Forums +
GitHub +
Editors: +
Josh Moore (University of Dundee (UoD)) +
Sébastien Besson (University of Dundee (UoD)) +
Constantin Pape (European Molecular Biology Laboratory (EMBL)) +
+
+
+ +
+
+
+

Abstract

+

This document contains next-generation file format (NGFF) + +specifications for storing bioimaging data in the cloud. +All specifications are submitted to the https://image.sc community for review.

+
+

Status of this document

+
+

+

The current released version of this specification is 0.4. Migration scripts +will be provided between numbered versions. Data written with these latest changes +(an "editor’s draft") will not necessarily be supported.

+
+
+ +
+

1. Introduction

+

Bioimaging science is at a crossroads. Currently, the drive to acquire more, +larger, preciser spatial measurements is unfortunately at odds with our ability +to structure and share those measurements with others. During a global pandemic +more than ever, we believe fervently that global, collaborative discovery as +opposed to the post-publication, "data-on-request" mode of operation is the +path forward. Bioimaging data should be shareable via open and commercial cloud +resources without the need to download entire datasets.

+

At the moment, that is not the norm. The plethora of data formats produced by +imaging systems are ill-suited to remote sharing. Individual scientists +typically lack the infrastructure they need to host these data themselves. When +they acquire images from elsewhere, time-consuming translations and data +cleaning are needed to interpret findings. Those same costs are multiplied when +gathering data into online repositories where curator time can be the limiting +factor before publication is possible. Without a common effort, each lab or +resource is left building the tools they need and maintaining that +infrastructure often without dedicated funding.

+

This document defines a specification for bioimaging data to make it possible +to enable the conversion of proprietary formats into a common, cloud-ready one. +Such next-generation file formats layout data so that individual portions, or +"chunks", of large data are reference-able eliminating the need to download +entire datasets.

+

1.1. Why "NGFF"?

+

A short description of what is needed for an imaging format is "a hierarchy +of n-dimensional (dense) arrays with metadata". This combination of features +is certainly provided by HDF5 from the HDF Group, which a number of +bioimaging formats do use. HDF5 and other larger binary structures, however, +are ill-suited for storage in the cloud where accessing individual chunks +of data by name rather than seeking through a large file is at the heart of +parallelization.

+

As a result, a number of formats have been developed more recently which provide +the basic data structure of an HDF5 file, but do so in a more cloud-friendly way. +In the PyData community, the Zarr [zarr] format was developed +for easily storing collections of NumPy arrays. In the ImageJ community, N5 [n5] was developed to work around +the limitations of HDF5 ("N5" was originally short for "Not-HDF5"). +Both of these formats permit storing individual chunks of data either locally in +separate files or in cloud-based object stores as separate keys.

+

A current effort is underway to unify the two similar specifications to provide a single binary +specification. The editor’s draft will soon be entering a request for comments (RFC) phase with the goal of having a first version early in 2021. As that +process comes to an end, this document will be updated.

+

1.2. OME-NGFF

+

The conventions and specifications defined in this document are designed to +enable next-generation file formats to represent the same bioimaging data +that can be represented in OME-TIFF and beyond. However, the conventions will also be usable by HDF5 and other sufficiently advanced +binary containers. Eventually, we hope, the moniker "next-generation" will no longer be +applicable, and this will simply be the most efficient, common, and useful representation +of bioimaging data, whether during acquisition or sharing in the cloud.

+

Note: The following text makes use of OME-Zarr [ome-zarr-py], the current prototype implementation, +for all examples.

+

1.3. Document conventions

+

The key words “MUST”, “MUST NOT”, “REQUIRED”, “SHALL”, “SHALL NOT”, “SHOULD”, “SHOULD NOT”, +“RECOMMENDED”, “MAY”, and “OPTIONAL” are to be interpreted as described in RFC 2119.

+

Some of the JSON examples in this document include commments. However, these are only for +clarity purposes and comments MUST NOT be included in JSON objects.

+

2. On-disk (or in-cloud) layout

+

An overview of the layout of an OME-Zarr fileset should make +understanding the following metadata sections easier. The hierarchy +is represented here as it would appear locally but could equally +be stored on a web server to be accessed via HTTP or in object storage +like S3 or GCS.

+

OME-Zarr is an implementation of the OME-NGFF specification using the Zarr +format. Arrays MUST be defined and stored in a hierarchical organization as +defined by the version 2 of the Zarr specification . +OME-NGFF metadata MUST be stored as attributes in the corresponding Zarr +groups.

+

2.1. Images

+

The following layout describes the expected Zarr hierarchy for images with +multiple levels of resolutions and optionally associated labels. +Note that the number of dimensions is variable between 2 and 5 and that axis names are arbitrary, see § 3.3 "multiscales" metadata for details. +For this example we assume an image with 5 dimensions and axes called t,c,z,y,x.

+
.                             # Root folder, potentially in S3,
+│                             # with a flat list of images by image ID.
+│
+├── 123.zarr                  # One image (id=123) converted to Zarr.
+│
+└── 456.zarr                  # Another image (id=456) converted to Zarr.
+    │
+    ├── .zgroup               # Each image is a Zarr group, or a folder, of other groups and arrays.
+    ├── .zattrs               # Group level attributes are stored in the .zattrs file and include
+    │                         # "multiscales" and "omero" (see below). In addition, the group level attributes
+    │                         # must also contain "_ARRAY_DIMENSIONS" if this group directly contains multi-scale arrays.
+    │
+    ├── 0                     # Each multiscale level is stored as a separate Zarr array,
+    │   ...                   # which is a folder containing chunk files which compose the array.
+    ├── n                     # The name of the array is arbitrary with the ordering defined by
+    │   │                     # by the "multiscales" metadata, but is often a sequence starting at 0.
+    │   │
+    │   ├── .zarray           # All image arrays must be up to 5-dimensional
+    │   │                     # with the axis of type time before type channel, before spatial axes.
+    │   │
+    │   └─ t                  # Chunks are stored with the nested directory layout.
+    │      └─ c               # All but the last chunk element are stored as directories.
+    │         └─ z            # The terminal chunk is a file. Together the directory and file names
+    │            └─ y         # provide the "chunk coordinate" (t, c, z, y, x), where the maximum coordinate
+    │               └─ x      # will be dimension_size / chunk_size.
+    │
+    └── labels
+        │
+        ├── .zgroup           # The labels group is a container which holds a list of labels to make the objects easily discoverable
+        │
+        ├── .zattrs           # All labels will be listed in .zattrs e.g. { "labels": [ "original/0" ] }
+        │                     # Each dimension of the label (t, c, z, y, x) should be either the same as the
+        │                     # corresponding dimension of the image, or 1 if that dimension of the label
+        │                     # is irrelevant.
+        │
+        └── original          # Intermediate folders are permitted but not necessary and currently contain no extra metadata.
+            │
+            └── 0             # Multiscale, labeled image. The name is unimportant but is registered in the "labels" group above.
+                ├── .zgroup   # Zarr Group which is both a multiscaled image as well as a labeled image.
+                ├── .zattrs   # Metadata of the related image and as well as display information under the "image-label" key.
+                │
+                ├── 0         # Each multiscale level is stored as a separate Zarr array, as above, but only integer values
+                │   ...       # are supported.
+                └── n
+
+

2.2. Tables

+ A table is a 2-dimensional data structure consisting of rows and columns. +Tables are an intuitive way of organizing data or metadata consisting of variables +and records, which are often organized as columns and rows, respectively. +

In OME-NGFF, a table is a Zarr group containing zero, one, or more Zarr arrays, where each +array represents one column of the table. Columns are ordered, and each column in a +table MUST have the same number of rows. While the table itself MUST be 2-dimensional, +the columns need not be. 1-dimensional columns will be typical. Nevertheless, there +may be use cases in which it is sensible to conceptualize a 2D, 3D, or higher-dimensional +array as a single column. A row of the table is the set of arrays obtained by slicing all +columns at the same index in the 0-th dimension.

+

Tables are located in the tables/ directory, in the root of an image’s +Zarr group, alongside the labels/ directory if one is present. The .zattrs file +of tables MUST contain the "tables" property, which indicates the names of the +tables as a JSON array.

+

The .zattrs file in each table of tables/ MAY contain the property "row-names", +and it MUST contain the two properties "column-order", and "annotated-data".

+

"row-names" is a string indicating which column, if any, contains the table’s row names. +An array of row names is a special column containing no duplicate or missing values. +For example, in a table in which the rows represent genes, it may be sensible to +designate the column called "geneID" as the "row-names" column. The "row-names" column +may be either an array of strings or an array of integers. The "row-names" column need not appear +first in the ordered array of column names. Either a missing "row-names" property or an empty +string value (i.e. "row-names": "") indicates that the table does not have row names.

+

"column-order" is a JSON array containing all of the column names, in order. Column names listed +in "column-order" MUST be strings, even if they are string representations of a digit, e.g. "1". +If the table contains no columns, then "column-order" MUST be an empty array.

+

In some cases, the table provides information about--in other words, annotates--a +particular dimension of an array. The "annotated-data" property describes these cases. +Tables MAY annotate arrays, and can only annotate arrays. If a table is annotating an +array, then the 0-th dimension of the column(s) (i.e. the number of rows in the table) +MUST correspond to the length of the dimension that is being annotated in the source +array. The required "annotated-data" property is a JSON array of Zarr array paths and +dimensions (0-based indexing), as shown below:

+
"annotated-data": [                  // A JSON array containing all array / dimension pairs this table is annotating
+    {
+        "array": "/path/to/arrayA",  // This table annotates the 0-th dimension of arrayA ...
+        "dimension": 0,
+    },
+    {
+        "array": "/path/to/arrayA",  // ... and the first dimension of arrayA ...
+        "dimension": 1,
+    },
+    {
+        "array": "/path/to/arrayB",  // ... and the first dimension of arrayB.
+        "dimension": 1,
+    }
+],
+
+

The "annotated-data" property MUST be present, but it MAY be an empty JSON array. If annotated-data contains exactly one entry, then the 0-th dimension of all arrays +(columns) in the table group SHOULD be chunked in the same way as the corresponding +dimension of the source array.

+

Note: The AnnData data model is based on a main array with additional "annotation" +tables that have the same number of rows, columns, or both, as the main array. AnnData +objects can be stored in the tables/ directory of an OME-NGFF file. See the AnnData +documentation for detailed recommendations on formatting AnnData data structures within Zarr.

+

There MAY be one or more intermediate directories between tables/ and a particular table. +These SHOULD NOT contain metadata, unless the intermediate directory represents the parent +directory for an AnnData object. The names of directories beneath the tables/ directory are arbitrary, except in the AnnData case.

+

The on-disk format of a table group looks like this:

+
.                            # Root folder, potentially in S3.
+|
+└── 123.zarr
+    |
+    ├── .zgroup
+    |
+    └──tables
+        |
+        ├── .zgroup
+        ├── .zattrs          # A table group’s .zattrs file MUST contain an attribute "tables" that lists all
+        |                    # tables in the group with relative paths, e.g., {"tables": ["table_1", "some/nested/table_2"]} 
+        |
+        └── table_1          # A table group holds arrays (called columns) with one or more dimensions, 
+            |                # where the number of rows n MUST be the same across all columns.
+            ├── .zgroup
+            ├── .zattrs      # .zattrs MAY contain "row-names", which states which column contains row names.
+            |                # .zattrs MUST contain "column-order", which is a list of the order of the columns.
+            |                # .zattrs MUST contain "annotated-data", which is specified above.
+            |                
+            |     
+            └── column_1     # The table group MAY hold zero, one, or an arbitrary number of columns.
+                |
+                ├── .zarray
+                |
+                ├── 0.0
+                |   ...
+                └── n.m
+
+

2.3. High-content screening

+

The following specification defines the hierarchy for a high-content screening +dataset. Three groups MUST be defined above the images:

+ +

A well row group SHOULD NOT be present if there are no images in the well row. +A well group SHOULD NOT be present if there are no images in the well.

+
.                             # Root folder, potentially in S3,
+│
+└── 5966.zarr                 # One plate (id=5966) converted to Zarr
+    ├── .zgroup
+    ├── .zattrs               # Implements "plate" specification
+    ├── A                     # First row of the plate
+    │   ├── .zgroup
+    │   │
+    │   ├── 1                 # First column of row A
+    │   │   ├── .zgroup
+    │   │   ├── .zattrs       # Implements "well" specification
+    │   │   │
+    │   │   ├── 0             # First field of view of well A1
+    │   │   │   │
+    │   │   │   ├── .zgroup
+    │   │   │   ├── .zattrs   # Implements "multiscales", "omero"
+    │   │   │   ├── 0
+    │   │   │   │   ...       # Resolution levels
+    │   │   │   ├── n
+    │   │   │   └── labels    # Labels (optional)
+    │   │   ├── ...           # Fields of view
+    │   │   └── m
+    │   ├── ...               # Columns
+    │   └── 12
+    ├── ...                   # Rows
+    └── H
+
+

3. Metadata

+

The various .zattrs files throughout the above array hierarchy may contain metadata +keys as specified below for discovering certain types of data, especially images.

+

3.1. "axes" metadata

+

"axes" describes the dimensions of a physical coordinate space. It is a list of dictionaries, where each dictionary describes a dimension (axis) and:

+ +

If part of § 3.3 "multiscales" metadata, the length of "axes" MUST be equal to the number of dimensions of the arrays that contain the image data.

+

3.2. "coordinateTransformations" metadata

+

"coordinateTransformations" describe a series of transformations that map between two coordinate spaces (defined by "axes"). +For example, to map a discrete data space of an array to the corresponding physical space. +It is a list of dictionaries. Each entry describes a single transformation and MUST contain the field "type". +The value of "type" MUST be one of the elements of the type column in the table below. +Additional fields for the entry depend on "type" and are defined by the column fields.

+ + + + + + + +
identity + + identity transformation, is the default transformation and is typically not explicitly defined +
translation + one of: "translation":List[float], "path":str + translation vector, stored either as a list of floats ("translation") or as binary data at a location in this container (path). The length of vector defines number of dimensions. | +
scale + one of: "scale":List[float], "path":str + scale vector, stored either as a list of floats (scale) or as binary data at a location in this container (path). The length of vector defines number of dimensions. | +
type + fields + description +
+

The transformations in the list are applied sequentially and in order.

+

3.3. "multiscales" metadata

+

Metadata about an image can be found under the "multiscales" key in the group-level metadata. Here, image refers to 2 to 5 dimensional data representing image or volumetric data with optional time or channel axes. It is stored in a multiple resolution representation.

+

"multiscales" contains a list of dictionaries where each entry describes a multiscale image.

+

Each "multiscales" dictionary MUST contain the field "axes", see § 3.1 "axes" metadata. +The length of "axes" must be between 2 and 5 and MUST be equal to the dimensionality of the zarr arrays storing the image data (see "datasets:path"). +The "axes" MUST contain 2 or 3 entries of "type:space" and MAY contain one additional entry of "type:time" and MAY contain one additional entry of "type:channel" or a null / custom type. +The order of the entries MUST correspond to the order of dimensions of the zarr arrays. In addition, the entries MUST be ordered by "type" where the "time" axis must come first (if present), followed by the "channel" or custom axis (if present) and the axes of type "space". +If there are three spatial axes where two correspond to the image plane ("yx") and images are stacked along the other (anisotropic) axis ("z"), the spatial axes SHOULD be ordered as "zyx".

+

Each "multiscales" dictionary MUST contain the field "datasets", which is a list of dictionaries describing the arrays storing the individual resolution levels. +Each dictionary in "datasets" MUST contain the field "path", whose value contains the path to the array for this resolution relative +to the current zarr group. The "path"s MUST be ordered from largest (i.e. highest resolution) to smallest.

+

Each "datasets" dictionary MUST have the same number of dimensions and MUST NOT have more than 5 dimensions. The number of dimensions and order MUST correspond to number and order of "axes". +Each dictionary in "datasets" MUST contain the field "coordinateTransformations", which contains a list of transformations that map the data coordinates to the physical coordinates (as specified by "axes") for this resolution level. +The transformations are defined according to § 3.2 "coordinateTransformations" metadata. The transformation MUST only be of type translation or scale. +They MUST contain exactly one scale transformation that specifies the pixel size in physical units or time duration. If scaling information is not available or applicable for one of the axes, the value MUST express the scaling factor between the current resolution and the first resolution for the given axis, defaulting to 1.0 if there is no downsampling along the axis. +It MAY contain exactly one translation that specifies the offset from the origin in physical units. If translation is given it MUST be listed after scale to ensure that it is given in physical coordinates. +The length of the scale and translation array MUST be the same as the length of "axes". +The requirements (only scale and translation, restrictions on order) are in place to provide a simple mapping from data coordinates to physical coordinates while being compatible with the general transformation spec.

+

Each "multiscales" dictionary MAY contain the field "coordinateTransformations", describing transformations that are applied to all resolution levels in the same manner. +The transformations MUST follow the same rules about allowed types, order, etc. as in "datasets:coordinateTransformations" and are applied after them. +They can for example be used to specify the scale for a dimension that is the same for all resolutions.

+

Each "multiscales" dictionary SHOULD contain the field "name". It SHOULD contain the field "version", which indicates the version of the multiscale metadata of this image (current version is 0.5-dev).

+

Each "multiscales" dictionary SHOULD contain the field "type", which gives the type of downscaling method used to generate the multiscale image pyramid. +It SHOULD contain the field "metadata", which contains a dictionary with additional information about the downscaling method.

+
{
+    "multiscales": [
+        {
+            "version": "0.5-dev",
+            "name": "example",
+            "axes": [
+                {"name": "t", "type": "time", "unit": "millisecond"},
+                {"name": "c", "type": "channel"},
+                {"name": "z", "type": "space", "unit": "micrometer"},
+                {"name": "y", "type": "space", "unit": "micrometer"},
+                {"name": "x", "type": "space", "unit": "micrometer"}
+            ],
+            "datasets": [
+                {
+                    "path": "0",
+                    "coordinateTransformations": [{
+                        // the voxel size for the first scale level (0.5 micrometer)
+                        "type": "scale",
+                        "scale": [1.0, 1.0, 0.5, 0.5, 0.5]
+                    }]
+                },
+                {
+                    "path": "1",
+                    "coordinateTransformations": [{
+                        // the voxel size for the second scale level (downscaled by a factor of 2 -> 1 micrometer)
+                        "type": "scale",
+                        "scale": [1.0, 1.0, 1.0, 1.0, 1.0]
+                    }]
+                },
+                {
+                    "path": "2",
+                    "coordinateTransformations": [{
+                        // the voxel size for the third scale level (downscaled by a factor of 4 -> 2 micrometer)
+                        "type": "scale",
+                        "scale": [1.0, 1.0, 2.0, 2.0, 2.0]
+                    }]
+                }
+            ],
+            "coordinateTransformations": [{
+                // the time unit (0.1 milliseconds), which is the same for each scale level
+                "type": "scale",
+                "scale": [0.1, 1.0, 1.0, 1.0, 1.0]
+            }],
+            "type": "gaussian",
+            "metadata": {
+                "description": "the fields in metadata depend on the downscaling implementation. Here, the parameters passed to the skimage function are given",
+                "method": "skimage.transform.pyramid_gaussian",
+                "version": "0.16.1",
+                "args": "[true]",
+                "kwargs": {"multichannel": true}
+            }
+        }
+    ]
+}
+

If only one multiscale is provided, use it. Otherwise, the user can choose by +name, using the first multiscale as a fallback:

+
datasets = []
+for named in multiscales:
+    if named["name"] == "3D":
+        datasets = [x["path"] for x in named["datasets"]]
+        break
+if not datasets:
+    # Use the first by default. Or perhaps choose based on chunk size.
+    datasets = [x["path"] for x in multiscales[0]["datasets"]]
+
+

3.4. "omero" metadata

+

Information specific to the channels of an image and how to render it +can be found under the "omero" key in the group-level metadata:

+
"id": 1,                              # ID in OMERO
+"name": "example.tif",                # Name as shown in the UI
+"version": "0.5-dev",                 # Current version
+"channels": [                         # Array matching the c dimension size
+    {
+        "active": true,
+        "coefficient": 1,
+        "color": "0000FF",
+        "family": "linear",
+        "inverted": false,
+        "label": "LaminB1",
+        "window": {
+            "end": 1500,
+            "max": 65535,
+            "min": 0,
+            "start": 0
+        }
+    }
+],
+"rdefs": {
+    "defaultT": 0,                    # First timepoint to show the user
+    "defaultZ": 118,                  # First Z section to show the user
+    "model": "color"                  # "color" or "greyscale"
+}
+
+

See https://docs.openmicroscopy.org/omero/5.6.1/developers/Web/WebGateway.html#imgdata +for more information.

+

3.5. "labels" metadata

+

The special group "labels" found under an image Zarr contains the key labels containing +the paths to label objects which can be found underneath the group:

+
{
+  "labels": [
+    "orphaned/0"
+  ]
+}
+
+

Unlisted groups MAY be labels.

+

3.6. "image-label" metadata

+

Groups containing the image-label dictionary represent an image segmentation +in which each unique pixel value represents a separate segmented object. image-label groups MUST also contain multiscales metadata and the two +"datasets" series MUST have the same number of entries.

+

The colors key defines a list of JSON objects describing the unique label +values. Each entry in the list MUST contain the key "label-value" with the +pixel value for that label. Additionally, the "rgba" key MAY be present, the +value for which is an RGBA unsigned-int 4-tuple: [uint8, uint8, uint8, uint8] All label-values must be unique. Clients who choose to not throw an error +should ignore all except the _last_ entry.

+

Some implementations may represent overlapping labels by using a specially assigned +value, for example the highest integer available in the pixel range.

+

The properties key defines a list of JSON objects which also describes the unique +label values. Each entry in the list MUST contain the key "label-value" with the +pixel value for that label. Additionally, an arbitrary number of key-value pairs +MAY be present for each label value denoting associated metadata. Not all label +values must share the same key-value pairs within the properties list.

+

The source key is an optional dictionary which contains information on the +image the label is associated with. If included it MAY include a key image whose value is the relative path to a Zarr image group. The default value is +"../../" since most labels are stored under a subgroup named "labels/" (see +above).

+
"image-label":
+  {
+    "version": "0.5-dev",
+    "colors": [
+      {
+        "label-value": 1,
+        "rgba": [255, 255, 255, 0]
+      },
+      {
+        "label-value": 4,
+        "rgba": [0, 255, 255, 128]
+      },
+      ...
+      ],
+    "properties": [
+      {
+        "label-value": 1,
+        "area (pixels)": 1200,
+        "class": "foo"
+
+      },
+      {
+        "label-value": 4,
+        "area (pixels)": 1650
+      },
+      ...
+      ]
+  },
+  "source": {
+    "image": "../../"
+  }
+]
+
+

3.7. "plate" metadata

+

For high-content screening datasets, the plate layout can be found under the +custom attributes of the plate group under the plate key in the group-level metadata.

+

The plate dictionary MAY contain an acquisitions key whose value MUST be a list of +JSON objects defining the acquisitions for a given plate to which wells can refer to. Each +acquisition object MUST contain an id key whose value MUST be an unique integer identifier +greater than or equal to 0 within the context of the plate to which fields of view can refer +to (see #well-md). +Each acquisition object SHOULD contain a name key whose value MUST be a string identifying +the name of the acquisition. Each acquisition object SHOULD contain a maximumfieldcount key whose value MUST be a positive integer indicating the maximum number of fields of view for the +acquisition. Each acquisition object MAY contain a description key whose value MUST be a +string specifying a description for the acquisition. Each acquisition object MAY contain +a starttime and/or endtime key whose values MUST be integer epoch timestamps specifying +the start and/or end timestamp of the acquisition.

+

The plate dictionary MUST contain a columns key whose value MUST be a list of JSON objects +defining the columns of the plate. Each column object defines the properties of +the column at the index of the object in the list. Each column in the physical plate +MUST be defined, even if no wells in the column are defined. Each column object MUST +contain a name key whose value is a string specifying the column name. The name MUST +contain only alphanumeric characters, MUST be case-sensitive, and MUST NOT be a duplicate of any +other name in the columns list. Care SHOULD be taken to avoid collisions on +case-insensitive filesystems (e.g. avoid using both Aa and aA).

+

The plate dictionary SHOULD contain a field_count key whose value MUST be a positive integer +defining the maximum number of fields per view across all wells.

+

The plate dictionary SHOULD contain a name key whose value MUST be a string defining the +name of the plate.

+

The plate dictionary MUST contain a rows key whose value MUST be a list of JSON objects +defining the rows of the plate. Each row object defines the properties of +the row at the index of the object in the list. Each row in the physical plate +MUST be defined, even if no wells in the row are defined. Each defined row MUST +contain a name key whose value MUST be a string defining the row name. The name MUST +contain only alphanumeric characters, MUST be case-sensitive, and MUST NOT be a duplicate of any +other name in the rows list. Care SHOULD be taken to avoid collisions on +case-insensitive filesystems (e.g. avoid using both Aa and aA).

+

The plate dictionary SHOULD contain a version key whose value MUST be a string specifying the +version of the plate specification.

+

The plate dictionary MUST contain a wells key whose value MUST be a list of JSON objects +defining the wells of the plate. Each well object MUST contain a path key whose value MUST +be a string specifying the path to the well subgroup. The path MUST consist of a name in +the rows list, a file separator (/), and a name from the columns list, in that order. +The path MUST NOT contain additional leading or trailing directories. +Each well object MUST contain both a rowIndex key whose value MUST be an integer identifying +the index into the rows list and a columnIndex key whose value MUST be an integer indentifying +the index into the columns list. rowIndex and columnIndex MUST be 0-based. The rowIndex, columnIndex, and path MUST all refer to the same row/column pair.

+

For example the following JSON object defines a plate with two acquisitions and +6 wells (2 rows and 3 columns), containing up to 2 fields of view per acquisition.

+
{
+    "plate": {
+        "acquisitions": [
+            {
+                "id": 1,
+                "maximumfieldcount": 2,
+                "name": "Meas_01(2012-07-31_10-41-12)",
+                "starttime": 1343731272000
+            },
+            {
+                "id": 2,
+                "maximumfieldcount": 2,
+                "name": "Meas_02(201207-31_11-56-41)",
+                "starttime": 1343735801000
+            }
+        ],
+        "columns": [
+            {
+                "name": "1"
+            },
+            {
+                "name": "2"
+            },
+            {
+                "name": "3"
+            }
+        ],
+        "field_count": 4,
+        "name": "test",
+        "rows": [
+            {
+                "name": "A"
+            },
+            {
+                "name": "B"
+            }
+        ],
+        "version": "0.5-dev",
+        "wells": [
+            {
+                "path": "A/1",
+                "rowIndex": 0,
+                "columnIndex": 0
+            },
+            {
+                "path": "A/2",
+                "rowIndex": 0,
+                "columnIndex": 1
+            },
+            {
+                "path": "A/3",
+                "rowIndex": 0,
+                "columnIndex": 2
+            },
+            {
+                "path": "B/1",
+                "rowIndex": 1,
+                "columnIndex": 0
+            },
+            {
+                "path": "B/2",
+                "rowIndex": 1,
+                "columnIndex": 1
+            },
+            {
+                "path": "B/3",
+                "rowIndex": 1,
+                "columnIndex": 2
+            }
+        ]
+    }
+}
+
+

The following JSON object defines a sparse plate with one acquisition and +2 wells in a 96 well plate, containing one field of view per acquisition.

+
{
+    "plate": {
+        "acquisitions": [
+            {
+                "id": 1,
+                "maximumfieldcount": 1,
+                "name": "single acquisition",
+                "starttime": 1343731272000
+            }
+        ],
+        "columns": [
+            {
+                "name": "1"
+            },
+            {
+                "name": "2"
+            },
+            {
+                "name": "3"
+            },
+            {
+                "name": "4"
+            },
+            {
+                "name": "5"
+            },
+            {
+                "name": "6"
+            },
+            {
+                "name": "7"
+            },
+            {
+                "name": "8"
+            },
+            {
+                "name": "9"
+            },
+            {
+                "name": "10"
+            },
+            {
+                "name": "11"
+            },
+            {
+                "name": "12"
+            }
+        ],
+        "field_count": 1,
+        "name": "sparse test",
+        "rows": [
+            {
+                "name": "A"
+            },
+            {
+                "name": "B"
+            },
+            {
+                "name": "C"
+            },
+            {
+                "name": "D"
+            },
+            {
+                "name": "E"
+            },
+            {
+                "name": "F"
+            },
+            {
+                "name": "G"
+            },
+            {
+                "name": "H"
+            }
+        ],
+        "version": "0.5-dev",
+        "wells": [
+            {
+                "path": "C/5",
+                "rowIndex": 2,
+                "columnIndex": 4
+            },
+            {
+                "path": "D/7",
+                "rowIndex": 3,
+                "columnIndex": 6
+            }
+        ]
+    }
+}
+
+

3.8. "well" metadata

+

For high-content screening datasets, the metadata about all fields of views +under a given well can be found under the "well" key in the attributes of the +well group.

+

The well dictionary MUST contain an images key whose value MUST be a list of JSON objects +specifying all fields of views for a given well. Each image object MUST contain a path key whose value MUST be a string specifying the path to the field of view. The path MUST contain only alphanumeric characters, MUST be case-sensitive, and MUST NOT be a duplicate +of any other path in the images list. If multiple acquisitions were performed in the plate, +it MUST contain an acquisition key whose value MUST be an integer identifying the acquisition +which MUST match one of the acquisition JSON objects defined in the plate metadata (see #plate-md).

+

The well dictionary SHOULD contain a version key whose value MUST be a string specifying the +version of the well specification.

+

For example the following JSON object defines a well with four fields of +view. The first two fields of view were part of the first acquisition while +the last two fields of view were part of the second acquisition.

+
{
+    "well": {
+        "images": [
+            {
+                "acquisition": 1,
+                "path": "0"
+            },
+            {
+                "acquisition": 1,
+                "path": "1"
+            },
+            {
+                "acquisition": 2,
+                "path": "2"
+            },
+            {
+                "acquisition": 2,
+                "path": "3"
+            }
+        ],
+        "version": "0.5-dev"
+    }
+}
+

The following JSON object defines a well with two fields of view in a plate with +four acquisitions. The first field is part of the first acquisition, and the second +field is part of the last acquisition.

+
{
+    "well": {
+        "images": [
+            {
+                "acquisition": 0,
+                "path": "0"
+            },
+            {
+                "acquisition": 3,
+                "path": "1"
+            }
+        ],
+        "version": "0.5-dev"
+    }
+}
+

4. Specification naming style

+

Multi-word keys in this specification should use the camelCase style. +NB: some parts of the specification don’t obey this convention as they +were added before this was adopted, but they should be updated in due course.

+

5. Implementations

+

Projects which support reading and/or writing OME-NGFF data include:

+
+
bigdataviewer-ome-zarr +
Fiji-plugin for reading OME-Zarr. +
bioformats2raw +
A performant, Bio-Formats image file format converter. +
omero-ms-zarr +
A microservice for OMERO.server that converts images stored in OMERO to OME-Zarr files on the fly, served via a web API. +
idr-zarr-tools +
A full workflow demonstrating the conversion of IDR images to OME-Zarr images on S3. +
OMERO CLI Zarr plugin +
An OMERO CLI plugin that converts images stored in OMERO.server into a local Zarr file. +
ome-zarr-py +
A napari plugin for reading ome-zarr files. +
vizarr +
A minimal, purely client-side program for viewing Zarr-based images with Viv & ImJoy. +
+

Diagram of related projects

+

All implementations prevent an equivalent representation of a dataset which can be downloaded or uploaded freely. An interactive +version of this diagram is available from the OME2020 Workshop. +Mouseover the blackboxes representing the implementations above to get a quick tip on how to use them.

+

Note: If you would like to see your project listed, please open an issue or PR on the ome/ngff repository.

+

6. Citing

+

Next-generation file format (NGFF) specifications for storing bioimaging data in the cloud. J. Moore, et al. Editors. Open Microscopy Environment Consortium, 8 February 2022. +This edition of the specification is https://ngff.openmicroscopy.org/0.4/. +The latest edition is available at https://ngff.openmicroscopy.org/latest/. (doi:10.5281/zenodo.4282107)

+

7. Version History

+ + + + + + + + + + + + + +
Revision + Date + Description +
0.4.0 + 2022-02-08 + multiscales: add axes type, units and coordinateTransformations +
0.4.0 + 2022-02-08 + plate: add rowIndex/columnIndex +
0.3.0 + 2021-08-24 + Add axes field to multiscale metadata +
0.2.0 + 2021-03-29 + Change chunk dimension separator to "/" +
0.1.4 + 2020-11-26 + Add HCS specification +
0.1.3 + 2020-09-14 + Add labels specification +
0.1.2 + 2020-05-07 + Add description of "omero" metadata +
0.1.1 + 2020-05-06 + Add info on the ordering of resolutions +
0.1.0 + 2020-04-20 + First version for internal demo +
+
+
+

Conformance

+

Document conventions

+

Conformance requirements are expressed + with a combination of descriptive assertions + and RFC 2119 terminology. + The key words “MUST”, “MUST NOT”, “REQUIRED”, “SHALL”, “SHALL NOT”, “SHOULD”, “SHOULD NOT”, “RECOMMENDED”, “MAY”, and “OPTIONAL” + in the normative parts of this document + are to be interpreted as described in RFC 2119. + However, for readability, + these words do not appear in all uppercase letters in this specification.

+

All of the text of this specification is normative + except sections explicitly marked as non-normative, examples, and notes. [RFC2119]

+

Examples in this specification are introduced with the words “for example” + or are set apart from the normative text + with class="example", + like this:

+
+ +

This is an example of an informative example.

+
+

Informative notes begin with the word “Note” + and are set apart from the normative text + with class="note", + like this:

+

Note, this is an informative note.

+
+

Conformant Algorithms

+

Requirements phrased in the imperative as part of algorithms + (such as "strip any leading space characters" + or "return false and abort these steps") + are to be interpreted with the meaning of the key word + ("must", "should", "may", etc) + used in introducing the algorithm.

+

Conformance requirements phrased as algorithms or specific steps + can be implemented in any manner, + so long as the end result is equivalent. + In particular, the algorithms defined in this specification + are intended to be easy to understand + and are not intended to be performant. + Implementers are encouraged to optimize.

+
+
+ +

Index

+

Terms defined by this specification

+ +

References

+

Normative References

+
+
[RFC2119] +
S. Bradner. Key words for use in RFCs to Indicate Requirement Levels. March 1997. Best Current Practice. URL: https://datatracker.ietf.org/doc/html/rfc2119 +
+

Informative References

+
+
[N5] +
John A. Bogovic; et al. N5---a scalable Java API for hierarchies of chunked n-dimensional tensors and structured meta-data. 2020. Informational. URL: https://github.com/saalfeldlab/n5/issues/62 +
[OME-ZARR-PY] +
OME; et al. ome-zarr-py: Experimental implementation of next-generation file format (NGFF) specifications for storing bioimaging data in the cloud.. 06 October 2020. Informational. URL: https://doi.org/10.5281/zenodo.4113931 +
[ZARR] +
Alistair Miles; et al. Zarr: An implementation of chunked, compressed, N-dimensional arrays for Python.. 06 October 2020. Informational. URL: https://doi.org/10.5281/zenodo.4069231 +
\ No newline at end of file