1majortom-core/
2├── COLLECTION.json
3├── METADATA/
4│   ├── collection.parquet
5│   ├── sample.parquet
6│   ├── sample__s2.parquet
7│   ├── sample__l8.parquet
8│   ├── sample__dem.parquet
9│   └── sample__wc.parquet
10└── DATA/
11    └── {tile_id}/
12        ├── s2/data.tif
13        ├── l8/data.tif
14        ├── dem/data.tif
15        └── wc/data.tif

1majortom-core/
2├── COLLECTION.json
3├── METADATA/
4│   ├── collection.parquet
5│   ├── sample.parquet
6│   ├── sample__s2.parquet
7│   ├── sample__l8.parquet
8│   ├── sample__dem.parquet
9│   └── sample__wc.parquet
10└── DATA/
11    └── {tile_id}/
12        ├── s2/data.tif
13        ├── l8/data.tif
14        ├── dem/data.tif
15        └── wc/data.tif

1import numpy as np
2import rasterio
3import matplotlib.pyplot as plt
4
5BASE = "https://data.source.coop/major-tom/core"
6IDX = 42  # ← change this to explore different tiles
7
8MODS = ["s2", "l8", "dem", "wc"]
9fig, axes = plt.subplots(1, 4, figsize=(20, 5))
10
11for j, mod in enumerate(MODS):
12    url = f"/vsicurl/{BASE}/DATA/{IDX}/{mod}/data.tif"
13    with rasterio.open(url) as src:
14        if mod in ("s2", "l8"):
15            rgb = src.read([4, 3, 2]).astype(np.float32)
16            axes[j].imshow(np.clip(rgb / 3000, 0, 1).transpose(1, 2, 0))
17        elif mod == "dem":
18            axes[j].imshow(src.read(1), cmap="terrain")
19        else:
20            axes[j].imshow(src.read(1))
21    axes[j].set_title(mod)
22    axes[j].axis("off")
23
24plt.tight_layout()
25plt.show()

1import numpy as np
2import rasterio
3import matplotlib.pyplot as plt
4
5BASE = "https://data.source.coop/major-tom/core"
6IDX = 42  # ← change this to explore different tiles
7
8MODS = ["s2", "l8", "dem", "wc"]
9fig, axes = plt.subplots(1, 4, figsize=(20, 5))
10
11for j, mod in enumerate(MODS):
12    url = f"/vsicurl/{BASE}/DATA/{IDX}/{mod}/data.tif"
13    with rasterio.open(url) as src:
14        if mod in ("s2", "l8"):
15            rgb = src.read([4, 3, 2]).astype(np.float32)
16            axes[j].imshow(np.clip(rgb / 3000, 0, 1).transpose(1, 2, 0))
17        elif mod == "dem":
18            axes[j].imshow(src.read(1), cmap="terrain")
19        else:
20            axes[j].imshow(src.read(1))
21    axes[j].set_title(mod)
22    axes[j].axis("off")
23
24plt.tight_layout()
25plt.show()

1@inproceedings{francis2024majortom,
2  title={Major TOM: Expandable Datasets for Earth Observation},
3  author={Francis, Alistair and Czerkawski, Mikolaj},
4  booktitle={IGARSS 2024},
5  pages={2935--2940},
6  year={2024},
7  doi={10.1109/IGARSS53475.2024.10640760}
8}

1@inproceedings{francis2024majortom,
2  title={Major TOM: Expandable Datasets for Earth Observation},
3  author={Francis, Alistair and Czerkawski, Mikolaj},
4  booktitle={IGARSS 2024},
5  pages={2935--2940},
6  year={2024},
7  doi={10.1109/IGARSS53475.2024.10640760}
8}