1import zarr
2
3root = zarr.open_group(
4  "s3://us-west-2.opendata.source.coop/tge-labs/spheer-mosaic/",
5  mode="r",
6  storage_options={"anon": True}
7)
8embeddings = root["embeddings"]
9
10# Access a spatial subset (year 2020 = index 3, all bands, y slice, x slice)
11data = embeddings[3, :, 6100:6356, 25000:25256]   # float32, shape (100, 256, 256)

1import zarr
2
3root = zarr.open_group(
4  "s3://us-west-2.opendata.source.coop/tge-labs/spheer-mosaic/",
5  mode="r",
6  storage_options={"anon": True}
7)
8embeddings = root["embeddings"]
9
10# Access a spatial subset (year 2020 = index 3, all bands, y slice, x slice)
11data = embeddings[3, :, 6100:6356, 25000:25256]   # float32, shape (100, 256, 256)

1import xarray as xr
2
3
4def fetch_point(ds, x, y, time=3):
5    """Fetch Spheer embedding for a single EPSG:32631 point (metres). Returns shape (100,)."""
6    return ds.isel(time=time).sel(x=x, y=y, method="nearest")
7
8
9def fetch_bbox(ds, bbox, time=3):
10    """Fetch Spheer embeddings for a bounding box (EPSG:32631 metres). Returns shape (100, y, x)."""
11    minx, miny, maxx, maxy = bbox
12    return ds.isel(time=time).sel(x=slice(minx, maxx), y=slice(maxy, miny))
13
14
15def main():
16    ds = xr.open_zarr(
17        "s3://us-west-2.opendata.source.coop/tge-labs/spheer-mosaic/",
18        storage_options={"anon": True}
19    )
20
21    # Fetch embedding for a single point (EPSG:32631 metres)
22    x, y = 754864.0, 5945103.0
23    point_data = fetch_point(ds, x, y).compute()
24    print(f"Point shape: {point_data.embeddings.shape}")
25
26    # Fetch embeddings for a bounding box (EPSG:32631 metres)
27    bbox = [754000.0, 5944000.0, 756000.0, 5946000.0]
28    bbox_data = fetch_bbox(ds, bbox).compute()
29    print(f"Bbox shape: {bbox_data.embeddings.shape}")
30
31
32if __name__ == "__main__":
33    main()

1import xarray as xr
2
3
4def fetch_point(ds, x, y, time=3):
5    """Fetch Spheer embedding for a single EPSG:32631 point (metres). Returns shape (100,)."""
6    return ds.isel(time=time).sel(x=x, y=y, method="nearest")
7
8
9def fetch_bbox(ds, bbox, time=3):
10    """Fetch Spheer embeddings for a bounding box (EPSG:32631 metres). Returns shape (100, y, x)."""
11    minx, miny, maxx, maxy = bbox
12    return ds.isel(time=time).sel(x=slice(minx, maxx), y=slice(maxy, miny))
13
14
15def main():
16    ds = xr.open_zarr(
17        "s3://us-west-2.opendata.source.coop/tge-labs/spheer-mosaic/",
18        storage_options={"anon": True}
19    )
20
21    # Fetch embedding for a single point (EPSG:32631 metres)
22    x, y = 754864.0, 5945103.0
23    point_data = fetch_point(ds, x, y).compute()
24    print(f"Point shape: {point_data.embeddings.shape}")
25
26    # Fetch embeddings for a bounding box (EPSG:32631 metres)
27    bbox = [754000.0, 5944000.0, 756000.0, 5946000.0]
28    bbox_data = fetch_bbox(ds, bbox).compute()
29    print(f"Bbox shape: {bbox_data.embeddings.shape}")
30
31
32if __name__ == "__main__":
33    main()

1import asyncio
2from affine import Affine
3from obstore.store import S3Store
4from zarr.storage import ObjectStore
5from zarr.api.asynchronous import open_group
6
7
8store = S3Store(
9    "us-west-2.opendata.source.coop",
10    prefix="tge-labs/spheer-mosaic/",
11    region="us-west-2",
12    skip_signature=True,
13)
14zarr_store = ObjectStore(store, read_only=True)
15
16
17def xy_to_rowcol(x, y, transform):
18    """Convert EPSG:32631 x/y (metres) to row/col using the inverse transform."""
19    inv = ~transform
20    col, row = inv * (x, y)
21    return int(row), int(col)
22
23
24def bbox_to_slices(bbox, transform):
25    """Convert a projected bounding box to row/col slices."""
26    minx, miny, maxx, maxy = bbox
27    inv = ~transform
28    col1, row1 = inv * (minx, maxy)
29    col2, row2 = inv * (maxx, miny)
30    row_start, row_stop = int(min(row1, row2)), int(max(row1, row2)) + 1
31    col_start, col_stop = int(min(col1, col2)), int(max(col1, col2)) + 1
32    return slice(row_start, row_stop), slice(col_start, col_stop)
33
34
35async def fetch_point(embeddings, transform, x, y, year=3):
36    """Fetch Spheer embedding for a single EPSG:32631 point. Returns shape (100,)."""
37    row, col = xy_to_rowcol(x, y, transform)
38    return await embeddings.getitem((year, slice(None), row, col))
39
40
41async def fetch_bbox(embeddings, transform, bbox, year=3):
42    """Fetch Spheer embeddings for a projected bounding box. Returns shape (100, rows, cols)."""
43    row_slice, col_slice = bbox_to_slices(bbox, transform)
44    return await embeddings.getitem((year, slice(None), row_slice, col_slice))
45
46
47async def main():
48    root = await open_group(store=zarr_store, mode="r")
49    embeddings = await root.get("embeddings")
50    transform = Affine(*root.metadata.attributes["spatial:transform"])
51
52    # Fetch embedding for a single point (EPSG:32631 metres)
53    x, y = 754864.0, 5945103.0
54    point_data = await fetch_point(embeddings, transform, x, y)
55    print(f"Point shape: {point_data.shape}")
56
57    # Fetch embeddings for a bounding box (EPSG:32631 metres)
58    bbox = [754000.0, 5944000.0, 756000.0, 5946000.0]
59    bbox_data = await fetch_bbox(embeddings, transform, bbox)
60    print(f"Bbox shape: {bbox_data.shape}")
61
62
63if __name__ == "__main__":
64    asyncio.run(main())

1import asyncio
2from affine import Affine
3from obstore.store import S3Store
4from zarr.storage import ObjectStore
5from zarr.api.asynchronous import open_group
6
7
8store = S3Store(
9    "us-west-2.opendata.source.coop",
10    prefix="tge-labs/spheer-mosaic/",
11    region="us-west-2",
12    skip_signature=True,
13)
14zarr_store = ObjectStore(store, read_only=True)
15
16
17def xy_to_rowcol(x, y, transform):
18    """Convert EPSG:32631 x/y (metres) to row/col using the inverse transform."""
19    inv = ~transform
20    col, row = inv * (x, y)
21    return int(row), int(col)
22
23
24def bbox_to_slices(bbox, transform):
25    """Convert a projected bounding box to row/col slices."""
26    minx, miny, maxx, maxy = bbox
27    inv = ~transform
28    col1, row1 = inv * (minx, maxy)
29    col2, row2 = inv * (maxx, miny)
30    row_start, row_stop = int(min(row1, row2)), int(max(row1, row2)) + 1
31    col_start, col_stop = int(min(col1, col2)), int(max(col1, col2)) + 1
32    return slice(row_start, row_stop), slice(col_start, col_stop)
33
34
35async def fetch_point(embeddings, transform, x, y, year=3):
36    """Fetch Spheer embedding for a single EPSG:32631 point. Returns shape (100,)."""
37    row, col = xy_to_rowcol(x, y, transform)
38    return await embeddings.getitem((year, slice(None), row, col))
39
40
41async def fetch_bbox(embeddings, transform, bbox, year=3):
42    """Fetch Spheer embeddings for a projected bounding box. Returns shape (100, rows, cols)."""
43    row_slice, col_slice = bbox_to_slices(bbox, transform)
44    return await embeddings.getitem((year, slice(None), row_slice, col_slice))
45
46
47async def main():
48    root = await open_group(store=zarr_store, mode="r")
49    embeddings = await root.get("embeddings")
50    transform = Affine(*root.metadata.attributes["spatial:transform"])
51
52    # Fetch embedding for a single point (EPSG:32631 metres)
53    x, y = 754864.0, 5945103.0
54    point_data = await fetch_point(embeddings, transform, x, y)
55    print(f"Point shape: {point_data.shape}")
56
57    # Fetch embeddings for a bounding box (EPSG:32631 metres)
58    bbox = [754000.0, 5944000.0, 756000.0, 5946000.0]
59    bbox_data = await fetch_bbox(embeddings, transform, bbox)
60    print(f"Bbox shape: {bbox_data.shape}")
61
62
63if __name__ == "__main__":
64    asyncio.run(main())

1@misc{spheer_fm_embeddings,
2  title  = {Spheer FM Embeddings},
3  author = {Spheer},
4  year   = {2026},
5  url    = {https://huggingface.co/datasets/spheer/spheer-fm-embeddings}
6}

1@misc{spheer_fm_embeddings,
2  title  = {Spheer FM Embeddings},
3  author = {Spheer},
4  year   = {2026},
5  url    = {https://huggingface.co/datasets/spheer/spheer-fm-embeddings}
6}