ocf-data-sampler 0.2.37__py3-none-any.whl → 0.3.0__py3-none-any.whl

ocf_data_sampler/load/load_dataset.py

@@ -25,7 +25,7 @@ def get_dataset_dict(
             zarr_path=input_config.gsp.zarr_path,
             boundaries_version=input_config.gsp.boundaries_version,
             public=input_config.gsp.public,
-        ).compute()
+        )
 
         if gsp_ids is None:
             # Remove national (gsp_id=0)

ocf_data_sampler/load/nwp/providers/cloudcasting.py

@@ -28,7 +28,7 @@ def open_cloudcasting(zarr_path: str | list[str]) -> xr.DataArray:
             [3] https://github.com/openclimatefix/sat_pred
     """
     # Open the data
-    ds = open_zarr_paths(zarr_path)
+    ds = open_zarr_paths(zarr_path, backend="tensorstore")
 
     # Rename
     ds = ds.rename(

ocf_data_sampler/load/nwp/providers/ecmwf.py

@@ -19,7 +19,7 @@ def open_ifs(zarr_path: str | list[str]) -> xr.DataArray:
     Returns:
         Xarray DataArray of the NWP data
     """
-    ds = open_zarr_paths(zarr_path)
+    ds = open_zarr_paths(zarr_path, backend="tensorstore")
 
     # LEGACY SUPPORT - rename variable to channel if it exists
     ds = ds.rename({"init_time": "init_time_utc", "variable": "channel"})

ocf_data_sampler/load/nwp/providers/gfs.py

@@ -23,7 +23,12 @@ def open_gfs(zarr_path: str | list[str], public: bool = False) -> xr.DataArray:
     _log.info("Loading NWP GFS data")
 
     # Open data
-    gfs: xr.Dataset = open_zarr_paths(zarr_path, time_dim="init_time_utc", public=public)
+    gfs: xr.Dataset = open_zarr_paths(
+        zarr_path,
+        time_dim="init_time_utc",
+        public=public,
+        backend="dask",
+    )
     nwp: xr.DataArray = gfs.to_array(dim="channel")
 
     del gfs

ocf_data_sampler/load/nwp/providers/icon.py

@@ -20,7 +20,7 @@ def open_icon_eu(zarr_path: str | list[str]) -> xr.DataArray:
         Xarray DataArray of the NWP data
     """
     # Open and check initially
-    ds = open_zarr_paths(zarr_path, time_dim="init_time_utc")
+    ds = open_zarr_paths(zarr_path, time_dim="init_time_utc", backend="dask")
 
     if "icon_eu_data" in ds.data_vars:
         nwp = ds["icon_eu_data"]

ocf_data_sampler/load/nwp/providers/ukv.py

@@ -19,7 +19,7 @@ def open_ukv(zarr_path: str | list[str]) -> xr.DataArray:
     Returns:
         Xarray DataArray of the NWP data
     """
-    ds = open_zarr_paths(zarr_path)
+    ds = open_zarr_paths(zarr_path, backend="tensorstore")
 
     ds = ds.rename(
         {

ocf_data_sampler/load/nwp/providers/utils.py

@@ -1,21 +1,48 @@
 """Utility functions for the NWP data processing."""
 
+from glob import glob
+
 import xarray as xr
+from xarray_tensorstore import open_zarr
+
+from ocf_data_sampler.load.open_tensorstore_zarrs import open_zarrs
 
 
 def open_zarr_paths(
-    zarr_path: str | list[str], time_dim: str = "init_time", public: bool = False,
+    zarr_path: str | list[str],
+    time_dim: str = "init_time",
+    public: bool = False,
+    backend: str = "dask",
 ) -> xr.Dataset:
     """Opens the NWP data.
 
     Args:
         zarr_path: Path to the zarr(s) to open
         time_dim: Name of the time dimension
-        public: Whether the data is public or private
+        public: Whether the data is public or private. Only available for the dask backend.
+        backend: The xarray backend to use.
 
     Returns:
         The opened Xarray Dataset
     """
+    if backend not in ["dask", "tensorstore"]:
+        raise ValueError(
+            f"Unsupported backend: {backend}. Supported backends are 'dask' and 'tensorstore'.",
+        )
+
+    if public and backend == "tensorstore":
+        raise ValueError("Public data is only supported with the 'dask' backend.")
+
+    if backend == "tensorstore":
+        ds = _tensostore_open_zarr_paths(zarr_path, time_dim)
+
+    elif backend == "dask":
+        ds = _dask_open_zarr_paths(zarr_path, time_dim, public)
+
+    return ds
+
+
+def _dask_open_zarr_paths(zarr_path: str | list[str], time_dim: str, public: bool) -> xr.Dataset:
     general_kwargs = {
         "engine": "zarr",
         "chunks": "auto",
@@ -26,7 +53,7 @@ def open_zarr_paths(
         # note this only works for s3 zarr paths at the moment
         general_kwargs["storage_options"] = {"anon": True}
 
-    if type(zarr_path) in [list, tuple] or "*" in str(zarr_path):  # Multi-file dataset
+    if isinstance(zarr_path, list | tuple) or "*" in str(zarr_path):  # Multi-file dataset
         ds = xr.open_mfdataset(
             zarr_path,
             concat_dim=time_dim,
@@ -41,3 +68,16 @@ def open_zarr_paths(
             **general_kwargs,
         )
     return ds
+
+
+def _tensostore_open_zarr_paths(zarr_path: str | list[str], time_dim: str) -> xr.Dataset:
+
+    if "*" in str(zarr_path):
+        zarr_path = sorted(glob(zarr_path))
+
+    if isinstance(zarr_path, list | tuple):
+        ds = open_zarrs(zarr_path, concat_dim=time_dim).sortby(time_dim)
+    else:
+        ds = open_zarr(zarr_path)
+    return ds
+

ocf_data_sampler/load/open_tensorstore_zarrs.py

@@ -0,0 +1,92 @@
+"""Open multiple zarrs with TensorStore.
+
+This extendds the functionality of xarray_tensorstore to open multiple zarr stores
+"""
+
+import os
+
+import tensorstore as ts
+import xarray as xr
+from xarray_tensorstore import (
+    _raise_if_mask_and_scale_used_for_data_vars,
+    _TensorStoreAdapter,
+    _zarr_spec_from_path,
+)
+
+
+def tensorstore_open_multi_zarrs(
+    paths: list[str],
+    data_vars: list[str],
+    concat_axes: list[int],
+    context: ts.Context,
+    write: bool,
+) -> dict[str, ts.TensorStore]:
+    """Open multiple zarrs with TensorStore.
+
+    Args:
+        paths: List of paths to zarr stores.
+        data_vars: List of data variable names to open.
+        concat_axes: List of axes along which to concatenate the data variables.
+        context: TensorStore context.
+        write: Whether to open the stores for writing.
+    """
+    arrays_list = []
+    for path in paths:
+        specs = {k: _zarr_spec_from_path(os.path.join(path, k)) for k in data_vars}
+        array_futures = {
+          k: ts.open(spec, read=True, write=write, context=context)
+          for k, spec in specs.items()
+        }
+        arrays_list.append({k: v.result() for k, v in array_futures.items()})
+
+    arrays = {}
+    for k, axis in zip(data_vars, concat_axes, strict=False):
+        datasets = [d[k] for d in arrays_list]
+        arrays[k] = ts.concat(datasets, axis=axis)
+
+    return arrays
+
+
+def open_zarrs(
+    paths: list[str],
+    concat_dim: str,
+    *,
+    context: ts.Context | None = None,
+    mask_and_scale: bool = True,
+    write: bool = False,
+) -> xr.Dataset:
+    """Open multiple zarrs with TensorStore.
+
+    Args:
+        paths: List of paths to zarr stores.
+        concat_dim: Dimension along which to concatenate the data variables.
+        context: TensorStore context.
+        mask_and_scale: Whether to mask and scale the data.
+        write: Whether to open the stores for writing.
+    """
+    if context is None:
+        context = ts.Context()
+
+    ds = xr.open_mfdataset(
+        paths,
+        concat_dim=concat_dim,
+        combine="nested",
+        mask_and_scale=mask_and_scale,
+        decode_timedelta=True,
+    )
+
+    if mask_and_scale:
+        # Data variables get replaced below with _TensorStoreAdapter arrays, which
+        # don't get masked or scaled. Raising an error avoids surprising users with
+        # incorrect data values.
+        _raise_if_mask_and_scale_used_for_data_vars(ds)
+
+    data_vars = list(ds.data_vars)
+
+    concat_axes = [ds[v].dims.index(concat_dim) for v in data_vars]
+
+    arrays = tensorstore_open_multi_zarrs(paths, data_vars, concat_axes, context, write)
+
+    new_data = {k: _TensorStoreAdapter(v) for k, v in arrays.items()}
+
+    return ds.copy(data=new_data)

ocf_data_sampler/load/satellite.py

@@ -1,6 +1,7 @@
 """Satellite loader."""
 import numpy as np
 import xarray as xr
+from xarray_tensorstore import open_zarr
 
 from ocf_data_sampler.load.utils import (
     check_time_unique_increasing,
@@ -8,39 +9,7 @@ from ocf_data_sampler.load.utils import (
     make_spatial_coords_increasing,
 )
 
-
-def get_single_sat_data(zarr_path: str) -> xr.Dataset:
-    """Helper function to open a zarr from either a local or GCP path.
-
-    Args:
-        zarr_path: path to a zarr file. Wildcards (*) are supported only for local paths
-                   GCS paths (gs://) do not support wildcards
-
-    Returns:
-        An xarray Dataset containing satellite data
-
-    Raises:
-        ValueError: If a wildcard (*) is used in a GCS (gs://) path
-    """
-    # Raise an error if a wildcard is used in a GCP path
-    if "gs://" in str(zarr_path) and "*" in str(zarr_path):
-        raise ValueError("Wildcard (*) paths are not supported for GCP (gs://) URLs")
-
-    # Handle multi-file dataset for local paths
-    if "*" in str(zarr_path):
-        ds = xr.open_mfdataset(
-            zarr_path,
-            engine="zarr",
-            concat_dim="time",
-            combine="nested",
-            chunks="auto",
-            join="override",
-        )
-        check_time_unique_increasing(ds.time)
-    else:
-        ds = xr.open_dataset(zarr_path, engine="zarr", chunks="auto")
-
-    return ds
+from .open_tensorstore_zarrs import open_zarrs
 
 
 def open_sat_data(zarr_path: str | list[str]) -> xr.DataArray:
@@ -52,14 +21,11 @@ def open_sat_data(zarr_path: str | list[str]) -> xr.DataArray:
     """
     # Open the data
     if isinstance(zarr_path, list | tuple):
-        ds = xr.combine_nested(
-            [get_single_sat_data(path) for path in zarr_path],
-            concat_dim="time",
-            combine_attrs="override",
-            join="override",
-        )
+        ds = open_zarrs(zarr_path, concat_dim="time")
     else:
-        ds = get_single_sat_data(zarr_path)
+        ds = open_zarr(zarr_path)
+
+    check_time_unique_increasing(ds.time)
 
     ds = ds.rename(
         {

ocf_data_sampler/load/utils.py

@@ -47,7 +47,7 @@ def get_xr_data_array_from_xr_dataset(ds: xr.Dataset) -> xr.DataArray:
     Args:
         ds: xr.Dataset to extract xr.DataArray from
     """
-    datavars = list(ds.var())
+    datavars = list(ds.data_vars)
     if len(datavars) != 1:
         raise ValueError("Cannot open as xr.DataArray: dataset contains multiple variables")
     return ds[datavars[0]]

ocf_data_sampler/select/find_contiguous_time_periods.py

@@ -242,6 +242,11 @@ def intersection_of_2_dataframes_of_periods(a: pd.DataFrame, b: pd.DataFrame) ->
     if a.empty or b.empty:
         return pd.DataFrame(columns=["start_dt", "end_dt"])
 
+    # Maybe switch these for efficiency in the next section. We will do the native python loop over
+    # the shorter dataframe
+    if len(a) > len(b):
+        a, b = b, a
+
     all_intersecting_periods = []
     for a_period in a.itertuples():
         # Five ways in which two periods may overlap:
@@ -250,12 +255,12 @@ def intersection_of_2_dataframes_of_periods(a: pd.DataFrame, b: pd.DataFrame) ->
         # In all five, `a` must always start before (or equal to) where `b` ends,
         # and `a` must always end after (or equal to) where `b` starts.
 
-        overlapping_periods = b[(a_period.start_dt <= b.end_dt) & (a_period.end_dt >= b.start_dt)]
-
         # There are two ways in which two periods may *not* overlap:
         # a: |---|        or        |---|
         # b:       |---|      |---|
-        # `overlapping` will not include periods which do *not* overlap.
+        # `overlapping_periods` will not include periods which do *not* overlap.
+
+        overlapping_periods = b[(a_period.start_dt <= b.end_dt) & (a_period.end_dt >= b.start_dt)]
 
         # Now find the intersection of each period in `overlapping_periods` with
         # the period from `a` that starts at `a_start_dt` and ends at `a_end_dt`.
@@ -269,5 +274,6 @@ def intersection_of_2_dataframes_of_periods(a: pd.DataFrame, b: pd.DataFrame) ->
 
         all_intersecting_periods.append(intersection)
 
+
     all_intersecting_periods = pd.concat(all_intersecting_periods)
     return all_intersecting_periods.sort_values(by="start_dt").reset_index(drop=True)

ocf_data_sampler/torch_datasets/datasets/pvnet_uk.py

@@ -270,6 +270,8 @@ class PVNetUKRegionalDataset(AbstractPVNetUKDataset):
     def __getitem__(self, idx: int) -> NumpySample:
         # Get the coordinates of the sample
 
+        idx = int(idx)
+
         if idx >= len(self):
             raise ValueError(f"Index {idx} out of range for dataset of length {len(self)}")
 

{ocf_data_sampler-0.2.37.dist-info → ocf_data_sampler-0.3.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: ocf-data-sampler
-Version: 0.2.37
+Version: 0.3.0
 Author: James Fulton, Peter Dudfield
 Author-email: Open Climate Fix team <info@openclimatefix.org>
 License: MIT License
@@ -44,6 +44,7 @@ Requires-Dist: pyproj
 Requires-Dist: pyaml_env
 Requires-Dist: pyresample
 Requires-Dist: h5netcdf
+Requires-Dist: xarray-tensorstore==0.1.5
 
 # ocf-data-sampler
  

{ocf_data_sampler-0.2.37.dist-info → ocf_data_sampler-0.3.0.dist-info}/RECORD

@@ -8,19 +8,20 @@ ocf_data_sampler/data/uk_gsp_locations_20220314.csv,sha256=RSh7DRh55E3n8lVAaWXGT
 ocf_data_sampler/data/uk_gsp_locations_20250109.csv,sha256=XZISFatnbpO9j8LwaxNKFzQSjs6hcHFsV8a9uDDpy2E,9055334
 ocf_data_sampler/load/__init__.py,sha256=-vQP9g0UOWdVbjEGyVX_ipa7R1btmiETIKAf6aw4d78,201
 ocf_data_sampler/load/gsp.py,sha256=IrTA6z9quN08imKGHJLf8gRktarxn1-utNMNFD0zWQs,2944
-ocf_data_sampler/load/load_dataset.py,sha256=WjB3DvHbDQQYYnPmDFOWg_TQPgARZ5pu8fiRZSGtIg0,2099
-ocf_data_sampler/load/satellite.py,sha256=Gsc3oyPydEZLy6slUDtIpBCYLxWy9P3pD1VyI4W9-2w,2944
+ocf_data_sampler/load/load_dataset.py,sha256=K8rWykjII-3g127If7WRRFivzHNx3SshCvZj4uQlf28,2089
+ocf_data_sampler/load/open_tensorstore_zarrs.py,sha256=_RHWe0GmrBSA9s1TH5I9VCMPpeZEsuRuhDt5Vyyx5Fo,2725
+ocf_data_sampler/load/satellite.py,sha256=RylkJz8avxdM5pK_liaTlD1DTboyPMgykXJ4_Ek9WBA,1840
 ocf_data_sampler/load/site.py,sha256=WtOy20VMHJIY0IwEemCdcecSDUGcVaLUown-4ixJw90,2147
-ocf_data_sampler/load/utils.py,sha256=sZ0-zzconcLkVQwAkCYrqKDo98Hrh5ChdiQJv5Bh91g,2040
+ocf_data_sampler/load/utils.py,sha256=AGL0aOOQPrgqNBTjlBtR7Qg1PyQov3DFJo-y198u8pY,2044
 ocf_data_sampler/load/nwp/__init__.py,sha256=SmcrnbygO5xtCKmGR4wtHrj-HI7nOAvnAtfuvRufBGQ,25
 ocf_data_sampler/load/nwp/nwp.py,sha256=0E9shei3Mq1N7F-fBlEKY5Hm0_kI7ysY_rffnWIshvk,3612
 ocf_data_sampler/load/nwp/providers/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-ocf_data_sampler/load/nwp/providers/cloudcasting.py,sha256=fozXpB3a2rNqQgnpRDC7xunxffh7Wwmc0kkCiYmDVJ4,1521
-ocf_data_sampler/load/nwp/providers/ecmwf.py,sha256=an-gXsZwkPQvRXeza1U_4MNU5yEnVm0_8tn03rxTudI,997
-ocf_data_sampler/load/nwp/providers/gfs.py,sha256=anNSITZ8MSCMnJhGdwUUJ1Ed-ArGKeqJ6HB6eJ814c8,1011
-ocf_data_sampler/load/nwp/providers/icon.py,sha256=BnY3vAa5pHn1cyrImj0ymaRFKHanNtfD9_JO-4p2IZY,1241
-ocf_data_sampler/load/nwp/providers/ukv.py,sha256=Ka1KFZcJYPwr5vuxo-xWGVQC0pudheqGBonUnbyJCMg,1016
-ocf_data_sampler/load/nwp/providers/utils.py,sha256=NrzE3JAtoc6oEywJHxPUdi_I4UJgJ_l5dxLZ4DLKvcg,1124
+ocf_data_sampler/load/nwp/providers/cloudcasting.py,sha256=EbhJO1QNtTtumuyh3I7HQa7xZ50WeD-SQe7AsMxZC84,1544
+ocf_data_sampler/load/nwp/providers/ecmwf.py,sha256=P7JqfssmQq8eHKKXaBexsxts325AnK-JF3bl5ymVCjY,1020
+ocf_data_sampler/load/nwp/providers/gfs.py,sha256=h6vm-Rfz1JGOE4P_fP1_XQJ3bugNbeNAIyt56N8B1Dc,1066
+ocf_data_sampler/load/nwp/providers/icon.py,sha256=iVZwLKRr_D74_kAu5MHir6pRKEfbTmIxFRZAxzmiYdI,1257
+ocf_data_sampler/load/nwp/providers/ukv.py,sha256=2i32VM9gnmWUpbL0qBSp_AKzuyKucXZPS8yklbcGlbc,1039
+ocf_data_sampler/load/nwp/providers/utils.py,sha256=cVwCiC8FqNpkZFSUGv1CRqIQlKdjx1sIsb2SIUlvWV8,2333
 ocf_data_sampler/numpy_sample/__init__.py,sha256=nY5C6CcuxiWZ_jrXRzWtN7WyKXhJImSiVTIG6Rz4B_4,401
 ocf_data_sampler/numpy_sample/collate.py,sha256=hoxIc5SoHoIs3Nx37aRZzWChpswjy9lHUgaKgHIoo80,2039
 ocf_data_sampler/numpy_sample/common_types.py,sha256=9CjYHkUTx0ObduWh43fhsybZCTXvexql7qC2ptMDoek,377
@@ -33,13 +34,13 @@ ocf_data_sampler/numpy_sample/sun_position.py,sha256=5tt-zNm6aRuZMsxZPaAxyg7Heik
 ocf_data_sampler/select/__init__.py,sha256=mK7Wu_-j9IXGTYrOuDf5yDDuU5a306b0iGKTAooNg_s,210
 ocf_data_sampler/select/dropout.py,sha256=9gPyDF7bGmvSoMjMPu1j0gTZFHNFqsT3ToIo9mFNA00,1565
 ocf_data_sampler/select/fill_time_periods.py,sha256=TlGxp1xiAqnhdWfLy0pv3FuZc00dtimjWdLzr4JoTGA,865
-ocf_data_sampler/select/find_contiguous_time_periods.py,sha256=8lkWsV5i7iLCVGqQ-PGZbvWxsz3wBvLO70GSf6WeR0k,11363
+ocf_data_sampler/select/find_contiguous_time_periods.py,sha256=etkr6LuB7zxkfzWJ6SgHiULdRuFzFlq5bOUNd257Qx4,11545
 ocf_data_sampler/select/geospatial.py,sha256=CDExkl36eZOKmdJPzUr_K0Wn3axHqv5nYo-EkSiINcc,5032
 ocf_data_sampler/select/location.py,sha256=AZvGR8y62opiW7zACGXjoOtBEWRfSLOZIA73O5Deu0c,1037
 ocf_data_sampler/select/select_spatial_slice.py,sha256=Hd4jGRUfIZRoWCirOQZeoLpaUnStB6KyFSTPX69wZLw,8790
 ocf_data_sampler/select/select_time_slice.py,sha256=HeHbwZ0CP03x0-LaJtpbSdtpLufwVTR73p6wH6O_PS8,5513
 ocf_data_sampler/torch_datasets/datasets/__init__.py,sha256=jfJSFcR0eO1AqeH7S3KnGjsBqVZT5w3oyi784PUR6Q0,146
-ocf_data_sampler/torch_datasets/datasets/pvnet_uk.py,sha256=GuwDB27Ja0o8OwLt0-7KzZBxg5AU1uIDjvKpLBUI90I,12233
+ocf_data_sampler/torch_datasets/datasets/pvnet_uk.py,sha256=v63goKEMI6UgBPnQCnIbxhFFdwuP_sxgcPYY6iNfGkc,12257
 ocf_data_sampler/torch_datasets/datasets/site.py,sha256=R9sYZz3e1zr8NAtlYQp8_DgI3wIfC7Zvaeo_73TyiA8,24936
 ocf_data_sampler/torch_datasets/sample/__init__.py,sha256=GL84vdZl_SjHDGVyh9Uekx2XhPYuZ0dnO3l6f6KXnHI,100
 ocf_data_sampler/torch_datasets/sample/base.py,sha256=cQ1oIyhdmlotejZK8B3Cw6MNvpdnBPD8G_o2h7Ye4Vc,2206
@@ -55,7 +56,7 @@ ocf_data_sampler/torch_datasets/utils/validation_utils.py,sha256=YqmT-lExWlI8_ul
 scripts/download_gsp_location_data.py,sha256=rRDXMoqX-RYY4jPdxhdlxJGhWdl6r245F5UARgKV6P4,3121
 scripts/refactor_site.py,sha256=skzvsPP0Cn9yTKndzkilyNcGz4DZ88ctvCJ0XrBdc2A,3135
 utils/compute_icon_mean_stddev.py,sha256=a1oWMRMnny39rV-dvu8rcx85sb4bXzPFrR1gkUr4Jpg,2296
-ocf_data_sampler-0.2.37.dist-info/METADATA,sha256=KejDSbMqRFuGL-XF8CVoM4QFqzgCbzzTEDrkY8C3VJc,12184
-ocf_data_sampler-0.2.37.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-ocf_data_sampler-0.2.37.dist-info/top_level.txt,sha256=LEFU4Uk-PEo72QGLAfnVZIUEm37Q8mKuMeg_Xk-p33g,31
-ocf_data_sampler-0.2.37.dist-info/RECORD,,
+ocf_data_sampler-0.3.0.dist-info/METADATA,sha256=Kq7LhwcpxOpfu4S4NOq-JHFJYI7eeeuxPleNPx6UMLE,12224
+ocf_data_sampler-0.3.0.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+ocf_data_sampler-0.3.0.dist-info/top_level.txt,sha256=LEFU4Uk-PEo72QGLAfnVZIUEm37Q8mKuMeg_Xk-p33g,31
+ocf_data_sampler-0.3.0.dist-info/RECORD,,