# -*- coding: utf-8 -*-
from warnings import warn
import numpy as np
from numcodecs.registry import codec_registry
from zarr.core import Array
from zarr.errors import (err_array_not_found, err_contains_array,
err_contains_group)
from zarr.n5 import N5Store
from zarr.storage import (DirectoryStore, ZipStore, contains_array,
contains_group, default_compressor, init_array,
normalize_storage_path)
[docs]def create(shape, chunks=True, dtype=None, compressor='default',
fill_value=0, order='C', store=None, synchronizer=None,
overwrite=False, path=None, chunk_store=None, filters=None,
cache_metadata=True, cache_attrs=True, read_only=False,
object_codec=None, **kwargs):
"""Create an array.
Parameters
----------
shape : int or tuple of ints
Array shape.
chunks : int or tuple of ints, optional
Chunk shape. If True, will be guessed from `shape` and `dtype`. If
False, will be set to `shape`, i.e., single chunk for the whole array.
If an int, the chunk size in each dimension will be given by the value
of `chunks`. Default is True.
dtype : string or dtype, optional
NumPy dtype.
compressor : Codec, optional
Primary compressor.
fill_value : object
Default value to use for uninitialized portions of the array.
order : {'C', 'F'}, optional
Memory layout to be used within each chunk.
store : MutableMapping or string
Store or path to directory in file system or name of zip file.
synchronizer : object, optional
Array synchronizer.
overwrite : bool, optional
If True, delete all pre-existing data in `store` at `path` before
creating the array.
path : string, optional
Path under which array is stored.
chunk_store : MutableMapping, optional
Separate storage for chunks. If not provided, `store` will be used
for storage of both chunks and metadata.
filters : sequence of Codecs, optional
Sequence of filters to use to encode chunk data prior to compression.
cache_metadata : bool, optional
If True, array configuration metadata will be cached for the
lifetime of the object. If False, array metadata will be reloaded
prior to all data access and modification operations (may incur
overhead depending on storage and data access pattern).
cache_attrs : bool, optional
If True (default), user attributes will be cached for attribute read
operations. If False, user attributes are reloaded from the store prior
to all attribute read operations.
read_only : bool, optional
True if array should be protected against modification.
object_codec : Codec, optional
A codec to encode object arrays, only needed if dtype=object.
Returns
-------
z : zarr.core.Array
Examples
--------
Create an array with default settings::
>>> import zarr
>>> z = zarr.create((10000, 10000), chunks=(1000, 1000))
>>> z
<zarr.core.Array (10000, 10000) float64>
Create an array with different some different configuration options::
>>> from numcodecs import Blosc
>>> compressor = Blosc(cname='zstd', clevel=1, shuffle=Blosc.BITSHUFFLE)
>>> z = zarr.create((10000, 10000), chunks=(1000, 1000), dtype='i1', order='F',
... compressor=compressor)
>>> z
<zarr.core.Array (10000, 10000) int8>
To create an array with object dtype requires a filter that can handle Python object
encoding, e.g., `MsgPack` or `Pickle` from `numcodecs`::
>>> from numcodecs import MsgPack
>>> z = zarr.create((10000, 10000), chunks=(1000, 1000), dtype=object,
... object_codec=MsgPack())
>>> z
<zarr.core.Array (10000, 10000) object>
Example with some filters, and also storing chunks separately from metadata::
>>> from numcodecs import Quantize, Adler32
>>> store, chunk_store = dict(), dict()
>>> z = zarr.create((10000, 10000), chunks=(1000, 1000), dtype='f8',
... filters=[Quantize(digits=2, dtype='f8'), Adler32()],
... store=store, chunk_store=chunk_store)
>>> z
<zarr.core.Array (10000, 10000) float64>
"""
# handle polymorphic store arg
store = normalize_store_arg(store)
# API compatibility with h5py
compressor, fill_value = _kwargs_compat(compressor, fill_value, kwargs)
# initialize array metadata
init_array(store, shape=shape, chunks=chunks, dtype=dtype, compressor=compressor,
fill_value=fill_value, order=order, overwrite=overwrite, path=path,
chunk_store=chunk_store, filters=filters, object_codec=object_codec)
# instantiate array
z = Array(store, path=path, chunk_store=chunk_store, synchronizer=synchronizer,
cache_metadata=cache_metadata, cache_attrs=cache_attrs, read_only=read_only)
return z
def normalize_store_arg(store, clobber=False, default=dict):
if store is None:
return default()
elif isinstance(store, str):
if store.endswith('.zip'):
mode = 'w' if clobber else 'r'
return ZipStore(store, mode=mode)
elif store.endswith('.n5'):
return N5Store(store)
else:
return DirectoryStore(store)
else:
return store
def _kwargs_compat(compressor, fill_value, kwargs):
# to be compatible with h5py, as well as backwards-compatible with Zarr
# 1.x, accept 'compression' and 'compression_opts' keyword arguments
if compressor != 'default':
# 'compressor' overrides 'compression'
if 'compression' in kwargs:
warn("'compression' keyword argument overridden by 'compressor'")
del kwargs['compression']
if 'compression_opts' in kwargs:
warn("'compression_opts' keyword argument overridden by 'compressor'")
del kwargs['compression_opts']
elif 'compression' in kwargs:
compression = kwargs.pop('compression')
compression_opts = kwargs.pop('compression_opts', None)
if compression is None or compression == 'none':
compressor = None
elif compression == 'default':
compressor = default_compressor
elif isinstance(compression, str):
codec_cls = codec_registry[compression]
# handle compression_opts
if isinstance(compression_opts, dict):
compressor = codec_cls(**compression_opts)
elif isinstance(compression_opts, (list, tuple)):
compressor = codec_cls(*compression_opts)
elif compression_opts is None:
compressor = codec_cls()
else:
# assume single argument, e.g., int
compressor = codec_cls(compression_opts)
# be lenient here if user gives compressor as 'compression'
elif hasattr(compression, 'get_config'):
compressor = compression
else:
raise ValueError('bad value for compression: %r' % compression)
# handle 'fillvalue'
if 'fillvalue' in kwargs:
# to be compatible with h5py, accept 'fillvalue' instead of
# 'fill_value'
fill_value = kwargs.pop('fillvalue')
# ignore other keyword arguments
for k in kwargs:
warn('ignoring keyword argument %r' % k)
return compressor, fill_value
[docs]def empty(shape, **kwargs):
"""Create an empty array.
For parameter definitions see :func:`zarr.creation.create`.
Notes
-----
The contents of an empty Zarr array are not defined. On attempting to
retrieve data from an empty Zarr array, any values may be returned,
and these are not guaranteed to be stable from one access to the next.
"""
return create(shape=shape, fill_value=None, **kwargs)
[docs]def zeros(shape, **kwargs):
"""Create an array, with zero being used as the default value for
uninitialized portions of the array.
For parameter definitions see :func:`zarr.creation.create`.
Examples
--------
>>> import zarr
>>> z = zarr.zeros((10000, 10000), chunks=(1000, 1000))
>>> z
<zarr.core.Array (10000, 10000) float64>
>>> z[:2, :2]
array([[0., 0.],
[0., 0.]])
"""
return create(shape=shape, fill_value=0, **kwargs)
[docs]def ones(shape, **kwargs):
"""Create an array, with one being used as the default value for
uninitialized portions of the array.
For parameter definitions see :func:`zarr.creation.create`.
Examples
--------
>>> import zarr
>>> z = zarr.ones((10000, 10000), chunks=(1000, 1000))
>>> z
<zarr.core.Array (10000, 10000) float64>
>>> z[:2, :2]
array([[1., 1.],
[1., 1.]])
"""
return create(shape=shape, fill_value=1, **kwargs)
[docs]def full(shape, fill_value, **kwargs):
"""Create an array, with `fill_value` being used as the default value for
uninitialized portions of the array.
For parameter definitions see :func:`zarr.creation.create`.
Examples
--------
>>> import zarr
>>> z = zarr.full((10000, 10000), chunks=(1000, 1000), fill_value=42)
>>> z
<zarr.core.Array (10000, 10000) float64>
>>> z[:2, :2]
array([[42., 42.],
[42., 42.]])
"""
return create(shape=shape, fill_value=fill_value, **kwargs)
def _get_shape_chunks(a):
shape = None
chunks = None
if hasattr(a, 'shape') and \
isinstance(a.shape, tuple):
shape = a.shape
if hasattr(a, 'chunks') and \
isinstance(a.chunks, tuple) and \
(len(a.chunks) == len(a.shape)):
chunks = a.chunks
elif hasattr(a, 'chunklen'):
# bcolz carray
chunks = (a.chunklen,) + a.shape[1:]
return shape, chunks
[docs]def array(data, **kwargs):
"""Create an array filled with `data`.
The `data` argument should be a NumPy array or array-like object. For
other parameter definitions see :func:`zarr.creation.create`.
Examples
--------
>>> import numpy as np
>>> import zarr
>>> a = np.arange(100000000).reshape(10000, 10000)
>>> z = zarr.array(a, chunks=(1000, 1000))
>>> z
<zarr.core.Array (10000, 10000) int64>
"""
# ensure data is array-like
if not hasattr(data, 'shape') or not hasattr(data, 'dtype'):
data = np.asanyarray(data)
# setup dtype
kw_dtype = kwargs.get('dtype', None)
if kw_dtype is None:
kwargs['dtype'] = data.dtype
else:
kwargs['dtype'] = kw_dtype
# setup shape and chunks
data_shape, data_chunks = _get_shape_chunks(data)
kwargs['shape'] = data_shape
kw_chunks = kwargs.get('chunks', None)
if kw_chunks is None:
kwargs['chunks'] = data_chunks
else:
kwargs['chunks'] = kw_chunks
# pop read-only to apply after storing the data
read_only = kwargs.pop('read_only', False)
# instantiate array
z = create(**kwargs)
# fill with data
z[...] = data
# set read_only property afterwards
z.read_only = read_only
return z
[docs]def open_array(store=None, mode='a', shape=None, chunks=True, dtype=None,
compressor='default', fill_value=0, order='C', synchronizer=None,
filters=None, cache_metadata=True, cache_attrs=True, path=None,
object_codec=None, chunk_store=None, **kwargs):
"""Open an array using file-mode-like semantics.
Parameters
----------
store : MutableMapping or string, optional
Store or path to directory in file system or name of zip file.
mode : {'r', 'r+', 'a', 'w', 'w-'}, optional
Persistence mode: 'r' means read only (must exist); 'r+' means
read/write (must exist); 'a' means read/write (create if doesn't
exist); 'w' means create (overwrite if exists); 'w-' means create
(fail if exists).
shape : int or tuple of ints, optional
Array shape.
chunks : int or tuple of ints, optional
Chunk shape. If True, will be guessed from `shape` and `dtype`. If
False, will be set to `shape`, i.e., single chunk for the whole array.
If an int, the chunk size in each dimension will be given by the value
of `chunks`. Default is True.
dtype : string or dtype, optional
NumPy dtype.
compressor : Codec, optional
Primary compressor.
fill_value : object, optional
Default value to use for uninitialized portions of the array.
order : {'C', 'F'}, optional
Memory layout to be used within each chunk.
synchronizer : object, optional
Array synchronizer.
filters : sequence, optional
Sequence of filters to use to encode chunk data prior to compression.
cache_metadata : bool, optional
If True, array configuration metadata will be cached for the
lifetime of the object. If False, array metadata will be reloaded
prior to all data access and modification operations (may incur
overhead depending on storage and data access pattern).
cache_attrs : bool, optional
If True (default), user attributes will be cached for attribute read
operations. If False, user attributes are reloaded from the store prior
to all attribute read operations.
path : string, optional
Array path within store.
object_codec : Codec, optional
A codec to encode object arrays, only needed if dtype=object.
chunk_store : MutableMapping or string, optional
Store or path to directory in file system or name of zip file.
Returns
-------
z : zarr.core.Array
Examples
--------
>>> import numpy as np
>>> import zarr
>>> z1 = zarr.open_array('data/example.zarr', mode='w', shape=(10000, 10000),
... chunks=(1000, 1000), fill_value=0)
>>> z1[:] = np.arange(100000000).reshape(10000, 10000)
>>> z1
<zarr.core.Array (10000, 10000) float64>
>>> z2 = zarr.open_array('data/example.zarr', mode='r')
>>> z2
<zarr.core.Array (10000, 10000) float64 read-only>
>>> np.all(z1[:] == z2[:])
True
Notes
-----
There is no need to close an array. Data are automatically flushed to the
file system.
"""
# use same mode semantics as h5py
# r : read only, must exist
# r+ : read/write, must exist
# w : create, delete if exists
# w- or x : create, fail if exists
# a : read/write if exists, create otherwise (default)
# handle polymorphic store arg
clobber = mode == 'w'
store = normalize_store_arg(store, clobber=clobber)
if chunk_store is not None:
chunk_store = normalize_store_arg(chunk_store, clobber=clobber)
path = normalize_storage_path(path)
# API compatibility with h5py
compressor, fill_value = _kwargs_compat(compressor, fill_value, kwargs)
# ensure fill_value of correct type
if fill_value is not None:
fill_value = np.array(fill_value, dtype=dtype)[()]
# ensure store is initialized
if mode in ['r', 'r+']:
if contains_group(store, path=path):
err_contains_group(path)
elif not contains_array(store, path=path):
err_array_not_found(path)
elif mode == 'w':
init_array(store, shape=shape, chunks=chunks, dtype=dtype,
compressor=compressor, fill_value=fill_value,
order=order, filters=filters, overwrite=True, path=path,
object_codec=object_codec, chunk_store=chunk_store)
elif mode == 'a':
if contains_group(store, path=path):
err_contains_group(path)
elif not contains_array(store, path=path):
init_array(store, shape=shape, chunks=chunks, dtype=dtype,
compressor=compressor, fill_value=fill_value,
order=order, filters=filters, path=path,
object_codec=object_codec, chunk_store=chunk_store)
elif mode in ['w-', 'x']:
if contains_group(store, path=path):
err_contains_group(path)
elif contains_array(store, path=path):
err_contains_array(path)
else:
init_array(store, shape=shape, chunks=chunks, dtype=dtype,
compressor=compressor, fill_value=fill_value,
order=order, filters=filters, path=path,
object_codec=object_codec, chunk_store=chunk_store)
# determine read only status
read_only = mode == 'r'
# instantiate array
z = Array(store, read_only=read_only, synchronizer=synchronizer,
cache_metadata=cache_metadata, cache_attrs=cache_attrs, path=path,
chunk_store=chunk_store)
return z
def _like_args(a, kwargs):
shape, chunks = _get_shape_chunks(a)
if shape is not None:
kwargs.setdefault('shape', shape)
if chunks is not None:
kwargs.setdefault('chunks', chunks)
if hasattr(a, 'dtype'):
kwargs.setdefault('dtype', a.dtype)
if isinstance(a, Array):
kwargs.setdefault('compressor', a.compressor)
kwargs.setdefault('order', a.order)
kwargs.setdefault('filters', a.filters)
else:
kwargs.setdefault('compressor', 'default')
kwargs.setdefault('order', 'C')
[docs]def empty_like(a, **kwargs):
"""Create an empty array like `a`."""
_like_args(a, kwargs)
return empty(**kwargs)
[docs]def zeros_like(a, **kwargs):
"""Create an array of zeros like `a`."""
_like_args(a, kwargs)
return zeros(**kwargs)
[docs]def ones_like(a, **kwargs):
"""Create an array of ones like `a`."""
_like_args(a, kwargs)
return ones(**kwargs)
[docs]def full_like(a, **kwargs):
"""Create a filled array like `a`."""
_like_args(a, kwargs)
if isinstance(a, Array):
kwargs.setdefault('fill_value', a.fill_value)
return full(**kwargs)
[docs]def open_like(a, path, **kwargs):
"""Open a persistent array like `a`."""
_like_args(a, kwargs)
if isinstance(a, Array):
kwargs.setdefault('fill_value', a.fill_value)
return open_array(path, **kwargs)