Source code for torchgeo.datasets.spacenet
# Copyright (c) Microsoft Corporation. All rights reserved.
# Licensed under the MIT License.
"""SpaceNet datasets."""
import abc
import copy
import glob
import math
import os
import re
from typing import Any, Callable, Optional
import fiona
import matplotlib.pyplot as plt
import numpy as np
import rasterio as rio
import torch
from fiona.errors import FionaValueError
from fiona.transform import transform_geom
from matplotlib.figure import Figure
from rasterio.crs import CRS
from rasterio.features import rasterize
from rasterio.transform import Affine
from torch import Tensor
from .geo import NonGeoDataset
from .utils import (
check_integrity,
download_radiant_mlhub_collection,
download_radiant_mlhub_dataset,
extract_archive,
percentile_normalization,
)
[docs]class SpaceNet(NonGeoDataset, abc.ABC):
"""Abstract base class for the SpaceNet datasets.
The `SpaceNet <https://spacenet.ai/datasets/>`__ datasets are a set of
datasets that all together contain >11M building footprints and ~20,000 km
of road labels mapped over high-resolution satellite imagery obtained from
a variety of sensors such as Worldview-2, Worldview-3 and Dove.
.. note::
The SpaceNet datasets require the following additional library to be installed:
* `radiant-mlhub <https://pypi.org/project/radiant-mlhub/>`_ to download the
imagery and labels from the Radiant Earth MLHub
"""
@property
@abc.abstractmethod
def dataset_id(self) -> str:
"""Dataset ID."""
@property
@abc.abstractmethod
def imagery(self) -> dict[str, str]:
"""Mapping of image identifier and filename."""
@property
@abc.abstractmethod
def label_glob(self) -> str:
"""Label filename."""
@property
@abc.abstractmethod
def collection_md5_dict(self) -> dict[str, str]:
"""Mapping of collection id and md5 checksum."""
@property
@abc.abstractmethod
def chip_size(self) -> dict[str, tuple[int, int]]:
"""Mapping of images and their chip size."""
[docs] def __init__(
self,
root: str,
image: str,
collections: list[str] = [],
transforms: Optional[Callable[[dict[str, Any]], dict[str, Any]]] = None,
download: bool = False,
api_key: Optional[str] = None,
checksum: bool = False,
) -> None:
"""Initialize a new SpaceNet Dataset instance.
Args:
root: root directory where dataset can be found
image: image selection
collections: collection selection
transforms: a function/transform that takes input sample and its target as
entry and returns a transformed version.
download: if True, download dataset and store it in the root directory.
api_key: a RadiantEarth MLHub API key to use for downloading the dataset
checksum: if True, check the MD5 of the downloaded files (may be slow)
Raises:
RuntimeError: if ``download=False`` but dataset is missing
"""
self.root = root
self.image = image # For testing
if collections:
for collection in collections:
assert collection in self.collection_md5_dict
self.collections = collections or list(self.collection_md5_dict.keys())
self.filename = self.imagery[image]
self.transforms = transforms
self.checksum = checksum
to_be_downloaded = self._check_integrity()
if to_be_downloaded:
if not download:
raise RuntimeError(
f"Dataset not found in `root={self.root}` and `download=False`, "
"either specify a different `root` directory or use "
"`download=True` to automatically download the dataset."
)
else:
self._download(to_be_downloaded, api_key)
self.files = self._load_files(root)
def _load_files(self, root: str) -> list[dict[str, str]]:
"""Return the paths of the files in the dataset.
Args:
root: root dir of dataset
Returns:
list of dicts containing paths for each pair of image and label
"""
files = []
for collection in self.collections:
images = glob.glob(os.path.join(root, collection, "*", self.filename))
images = sorted(images)
for imgpath in images:
lbl_path = os.path.join(
f"{os.path.dirname(imgpath)}-labels", self.label_glob
)
files.append({"image_path": imgpath, "label_path": lbl_path})
return files
def _load_image(self, path: str) -> tuple[Tensor, Affine, CRS]:
"""Load a single image.
Args:
path: path to the image
Returns:
the image
"""
filename = os.path.join(path)
with rio.open(filename) as img:
array = img.read().astype(np.int32)
tensor = torch.from_numpy(array).float()
return tensor, img.transform, img.crs
def _load_mask(
self, path: str, tfm: Affine, raster_crs: CRS, shape: tuple[int, int]
) -> Tensor:
"""Rasterizes the dataset's labels (in geojson format).
Args:
path: path to the label
tfm: transform of corresponding image
shape: shape of corresponding image
Returns:
Tensor: label tensor
"""
try:
with fiona.open(path) as src:
vector_crs = CRS(src.crs)
if raster_crs == vector_crs:
labels = [feature["geometry"] for feature in src]
else:
labels = [
transform_geom(
vector_crs.to_string(),
raster_crs.to_string(),
feature["geometry"],
)
for feature in src
]
except FionaValueError:
labels = []
if not labels:
mask_data = np.zeros(shape=shape)
else:
mask_data = rasterize(
labels,
out_shape=shape,
fill=0, # nodata value
transform=tfm,
all_touched=False,
dtype=np.uint8,
)
mask = torch.from_numpy(mask_data).long()
return mask
[docs] def __len__(self) -> int:
"""Return the number of samples in the dataset.
Returns:
length of the dataset
"""
return len(self.files)
[docs] def __getitem__(self, index: int) -> dict[str, Tensor]:
"""Return an index within the dataset.
Args:
index: index to return
Returns:
data and label at that index
"""
files = self.files[index]
img, tfm, raster_crs = self._load_image(files["image_path"])
h, w = img.shape[1:]
mask = self._load_mask(files["label_path"], tfm, raster_crs, (h, w))
ch, cw = self.chip_size[self.image]
sample = {"image": img[:, :ch, :cw], "mask": mask[:ch, :cw]}
if self.transforms is not None:
sample = self.transforms(sample)
return sample
def _check_integrity(self) -> list[str]:
"""Checks the integrity of the dataset structure.
Returns:
List of collections to be downloaded
"""
# Check if collections exist
missing_collections = []
for collection in self.collections:
stacpath = os.path.join(self.root, collection, "collection.json")
if not os.path.exists(stacpath):
missing_collections.append(collection)
if not missing_collections:
return []
to_be_downloaded = []
for collection in missing_collections:
archive_path = os.path.join(self.root, f"{collection}.tar.gz")
if os.path.exists(archive_path):
print(f"Found {collection} archive")
if (
self.checksum
and check_integrity(
archive_path, self.collection_md5_dict[collection]
)
or not self.checksum
):
print("Extracting...")
extract_archive(archive_path)
else:
print(f"Collection {collection} is corrupted")
to_be_downloaded.append(collection)
else:
print(f"{collection} not found")
to_be_downloaded.append(collection)
return to_be_downloaded
def _download(self, collections: list[str], api_key: Optional[str] = None) -> None:
"""Download the dataset and extract it.
Args:
collections: Collections to be downloaded
api_key: a RadiantEarth MLHub API key to use for downloading the dataset
Raises:
RuntimeError: if download doesn't work correctly or checksums don't match
"""
for collection in collections:
download_radiant_mlhub_collection(collection, self.root, api_key)
archive_path = os.path.join(self.root, f"{collection}.tar.gz")
if (
not self.checksum
or not check_integrity(
archive_path, self.collection_md5_dict[collection]
)
) and self.checksum:
raise RuntimeError(f"Collection {collection} corrupted")
print("Extracting...")
extract_archive(archive_path)
[docs] def plot(
self,
sample: dict[str, Tensor],
show_titles: bool = True,
suptitle: Optional[str] = None,
) -> Figure:
"""Plot a sample from the dataset.
Args:
sample: a sample returned by :meth:`__getitem__`
show_titles: flag indicating whether to show titles above each panel
suptitle: optional string to use as a suptitle
Returns:
a matplotlib Figure with the rendered sample
.. versionadded:: 0.2
"""
# image can be 1 channel or >3 channels
if sample["image"].shape[0] == 1:
image = np.rollaxis(sample["image"].numpy(), 0, 3)
else:
image = np.rollaxis(sample["image"][:3].numpy(), 0, 3)
image = percentile_normalization(image, axis=(0, 1))
ncols = 1
show_mask = "mask" in sample
show_predictions = "prediction" in sample
if show_mask:
mask = sample["mask"].numpy()
ncols += 1
if show_predictions:
prediction = sample["prediction"].numpy()
ncols += 1
fig, axs = plt.subplots(ncols=ncols, figsize=(ncols * 8, 8))
if not isinstance(axs, np.ndarray):
axs = [axs]
axs[0].imshow(image)
axs[0].axis("off")
if show_titles:
axs[0].set_title("Image")
if show_mask:
axs[1].imshow(mask, interpolation="none")
axs[1].axis("off")
if show_titles:
axs[1].set_title("Label")
if show_predictions:
axs[2].imshow(prediction, interpolation="none")
axs[2].axis("off")
if show_titles:
axs[2].set_title("Prediction")
if suptitle is not None:
plt.suptitle(suptitle)
return fig
[docs]class SpaceNet1(SpaceNet):
"""SpaceNet 1: Building Detection v1 Dataset.
`SpaceNet 1 <https://spacenet.ai/spacenet-buildings-dataset-v1/>`_
is a dataset of building footprints over the city of Rio de Janeiro.
Dataset features:
* No. of images: 6940 (8 Band) + 6940 (RGB)
* No. of polygons: 382,534 building labels
* Area Coverage: 2544 sq km
* GSD: 1 m (8 band), 50 cm (rgb)
* Chip size: 101 x 110 (8 band), 406 x 438 (rgb)
Dataset format:
* Imagery - Worldview-2 GeoTIFFs
* 8Band.tif (Multispectral)
* RGB.tif (Pansharpened RGB)
* Labels - GeoJSON
* labels.geojson
If you use this dataset in your research, please cite the following paper:
* https://arxiv.org/abs/1807.01232
"""
dataset_id = "spacenet1"
imagery = {"rgb": "RGB.tif", "8band": "8Band.tif"}
chip_size = {"rgb": (406, 438), "8band": (101, 110)}
label_glob = "labels.geojson"
collection_md5_dict = {"sn1_AOI_1_RIO": "e6ea35331636fa0c036c04b3d1cbf226"}
[docs] def __init__(
self,
root: str = "data",
image: str = "rgb",
transforms: Optional[Callable[[dict[str, Any]], dict[str, Any]]] = None,
download: bool = False,
api_key: Optional[str] = None,
checksum: bool = False,
) -> None:
"""Initialize a new SpaceNet 1 Dataset instance.
Args:
root: root directory where dataset can be found
image: image selection which must be "rgb" or "8band"
transforms: a function/transform that takes input sample and its target as
entry and returns a transformed version.
download: if True, download dataset and store it in the root directory.
api_key: a RadiantEarth MLHub API key to use for downloading the dataset
checksum: if True, check the MD5 of the downloaded files (may be slow)
Raises:
RuntimeError: if ``download=False`` but dataset is missing
"""
collections = ["sn1_AOI_1_RIO"]
assert image in {"rgb", "8band"}
super().__init__(
root, image, collections, transforms, download, api_key, checksum
)
[docs]class SpaceNet2(SpaceNet):
r"""SpaceNet 2: Building Detection v2 Dataset.
`SpaceNet 2 <https://spacenet.ai/spacenet-buildings-dataset-v2/>`_
is a dataset of building footprints over the cities of Las Vegas,
Paris, Shanghai and Khartoum.
Collection features:
+------------+---------------------+------------+------------+
| AOI | Area (km\ :sup:`2`\)| # Images | # Buildings|
+============+=====================+============+============+
| Las Vegas | 216 | 3850 | 151,367 |
+------------+---------------------+------------+------------+
| Paris | 1030 | 1148 | 23,816 |
+------------+---------------------+------------+------------+
| Shanghai | 1000 | 4582 | 92,015 |
+------------+---------------------+------------+------------+
| Khartoum | 765 | 1012 | 35,503 |
+------------+---------------------+------------+------------+
Imagery features:
.. list-table::
:widths: 10 10 10 10 10
:header-rows: 1
:stub-columns: 1
* -
- PAN
- MS
- PS-MS
- PS-RGB
* - GSD (m)
- 0.31
- 1.24
- 0.30
- 0.30
* - Chip size (px)
- 650 x 650
- 162 x 162
- 650 x 650
- 650 x 650
Dataset format:
* Imagery - Worldview-3 GeoTIFFs
* PAN.tif (Panchromatic)
* MS.tif (Multispectral)
* PS-MS (Pansharpened Multispectral)
* PS-RGB (Pansharpened RGB)
* Labels - GeoJSON
* label.geojson
If you use this dataset in your research, please cite the following paper:
* https://arxiv.org/abs/1807.01232
"""
dataset_id = "spacenet2"
collection_md5_dict = {
"sn2_AOI_2_Vegas": "a5a8de355290783b88ac4d69c7ef0694",
"sn2_AOI_3_Paris": "8299186b7bbfb9a256d515bad1b7f146",
"sn2_AOI_4_Shanghai": "4e3e80f2f437faca10ca2e6e6df0ef99",
"sn2_AOI_5_Khartoum": "8070ff9050f94cd9f0efe9417205d7c3",
}
imagery = {
"MS": "MS.tif",
"PAN": "PAN.tif",
"PS-MS": "PS-MS.tif",
"PS-RGB": "PS-RGB.tif",
}
chip_size = {
"MS": (162, 162),
"PAN": (650, 650),
"PS-MS": (650, 650),
"PS-RGB": (650, 650),
}
label_glob = "label.geojson"
[docs] def __init__(
self,
root: str = "data",
image: str = "PS-RGB",
collections: list[str] = [],
transforms: Optional[Callable[[dict[str, Any]], dict[str, Any]]] = None,
download: bool = False,
api_key: Optional[str] = None,
checksum: bool = False,
) -> None:
"""Initialize a new SpaceNet 2 Dataset instance.
Args:
root: root directory where dataset can be found
image: image selection which must be in ["MS", "PAN", "PS-MS", "PS-RGB"]
collections: collection selection which must be a subset of:
[sn2_AOI_2_Vegas, sn2_AOI_3_Paris, sn2_AOI_4_Shanghai,
sn2_AOI_5_Khartoum]. If unspecified, all collections will be
used.
transforms: a function/transform that takes input sample and its target as
entry and returns a transformed version
download: if True, download dataset and store it in the root directory.
api_key: a RadiantEarth MLHub API key to use for downloading the dataset
checksum: if True, check the MD5 of the downloaded files (may be slow)
Raises:
RuntimeError: if ``download=False`` but dataset is missing
"""
assert image in {"MS", "PAN", "PS-MS", "PS-RGB"}
super().__init__(
root, image, collections, transforms, download, api_key, checksum
)
[docs]class SpaceNet3(SpaceNet):
r"""SpaceNet 3: Road Network Detection.
`SpaceNet 3 <https://spacenet.ai/spacenet-roads-dataset/>`_
is a dataset of road networks over the cities of Las Vegas, Paris, Shanghai,
and Khartoum.
Collection features:
+------------+---------------------+------------+---------------------------+
| AOI | Area (km\ :sup:`2`\)| # Images | # Road Network Labels (km)|
+============+=====================+============+===========================+
| Vegas | 216 | 854 | 3685 |
+------------+---------------------+------------+---------------------------+
| Paris | 1030 | 257 | 425 |
+------------+---------------------+------------+---------------------------+
| Shanghai | 1000 | 1028 | 3537 |
+------------+---------------------+------------+---------------------------+
| Khartoum | 765 | 283 | 1030 |
+------------+---------------------+------------+---------------------------+
Imagery features:
.. list-table::
:widths: 10 10 10 10 10
:header-rows: 1
:stub-columns: 1
* -
- PAN
- MS
- PS-MS
- PS-RGB
* - GSD (m)
- 0.31
- 1.24
- 0.30
- 0.30
* - Chip size (px)
- 1300 x 1300
- 325 x 325
- 1300 x 1300
- 1300 x 1300
Dataset format:
* Imagery - Worldview-3 GeoTIFFs
* PAN.tif (Panchromatic)
* MS.tif (Multispectral)
* PS-MS (Pansharpened Multispectral)
* PS-RGB (Pansharpened RGB)
* Labels - GeoJSON
* labels.geojson
If you use this dataset in your research, please cite the following paper:
* https://arxiv.org/abs/1807.01232
.. versionadded:: 0.3
"""
dataset_id = "spacenet3"
collection_md5_dict = {
"sn3_AOI_2_Vegas": "8ce7e6abffb8849eb88885035f061ee8",
"sn3_AOI_3_Paris": "90b9ebd64cd83dc8d3d4773f45050d8f",
"sn3_AOI_4_Shanghai": "3ea291df34548962dfba8b5ed37d700c",
"sn3_AOI_5_Khartoum": "b8d549ac9a6d7456c0f7a8e6de23d9f9",
}
imagery = {
"MS": "MS.tif",
"PAN": "PAN.tif",
"PS-MS": "PS-MS.tif",
"PS-RGB": "PS-RGB.tif",
}
chip_size = {
"MS": (325, 325),
"PAN": (1300, 1300),
"PS-MS": (1300, 1300),
"PS-RGB": (1300, 1300),
}
label_glob = "labels.geojson"
[docs] def __init__(
self,
root: str = "data",
image: str = "PS-RGB",
speed_mask: Optional[bool] = False,
collections: list[str] = [],
transforms: Optional[Callable[[dict[str, Any]], dict[str, Any]]] = None,
download: bool = False,
api_key: Optional[str] = None,
checksum: bool = False,
) -> None:
"""Initialize a new SpaceNet 3 Dataset instance.
Args:
root: root directory where dataset can be found
image: image selection which must be in ["MS", "PAN", "PS-MS", "PS-RGB"]
speed_mask: use multi-class speed mask (created by binning roads at
10 mph increments) as label if true, else use binary mask
collections: collection selection which must be a subset of:
[sn3_AOI_2_Vegas, sn3_AOI_3_Paris, sn3_AOI_4_Shanghai,
sn3_AOI_5_Khartoum]. If unspecified, all collections will be
used.
transforms: a function/transform that takes input sample and its target as
entry and returns a transformed version
download: if True, download dataset and store it in the root directory.
api_key: a RadiantEarth MLHub API key to use for downloading the dataset
checksum: if True, check the MD5 of the downloaded files (may be slow)
Raises:
RuntimeError: if ``download=False`` but dataset is missing
"""
assert image in {"MS", "PAN", "PS-MS", "PS-RGB"}
self.speed_mask = speed_mask
super().__init__(
root, image, collections, transforms, download, api_key, checksum
)
def _load_mask(
self, path: str, tfm: Affine, raster_crs: CRS, shape: tuple[int, int]
) -> Tensor:
"""Rasterizes the dataset's labels (in geojson format).
Args:
path: path to the label
tfm: transform of corresponding image
shape: shape of corresponding image
Returns:
Tensor: label tensor
"""
min_speed_bin = 1
max_speed_bin = 65
speed_arr_bin = np.arange(min_speed_bin, max_speed_bin + 1)
bin_size_mph = 10.0
speed_cls_arr: "np.typing.NDArray[np.int_]" = np.array(
[math.ceil(s / bin_size_mph) for s in speed_arr_bin]
)
try:
with fiona.open(path) as src:
vector_crs = CRS(src.crs)
labels = []
for feature in src:
if raster_crs != vector_crs:
geom = transform_geom(
vector_crs.to_string(),
raster_crs.to_string(),
feature["geometry"],
)
else:
geom = feature["geometry"]
if self.speed_mask:
val = speed_cls_arr[
int(feature["properties"]["inferred_speed_mph"]) - 1
]
else:
val = 1
labels.append((geom, val))
except FionaValueError:
labels = []
if not labels:
mask_data = np.zeros(shape=shape)
else:
mask_data = rasterize(
labels,
out_shape=shape,
fill=0, # nodata value
transform=tfm,
all_touched=False,
dtype=np.uint8,
)
mask = torch.from_numpy(mask_data).long()
return mask
[docs] def plot(
self,
sample: dict[str, Tensor],
show_titles: bool = True,
suptitle: Optional[str] = None,
) -> Figure:
"""Plot a sample from the dataset.
Args:
sample: a sample returned by :meth:`SpaceNet.__getitem__`
show_titles: flag indicating whether to show titles above each panel
suptitle: optional string to use as a suptitle
Returns:
a matplotlib Figure with the rendered sample
"""
# image can be 1 channel or >3 channels
if sample["image"].shape[0] == 1:
image = np.rollaxis(sample["image"].numpy(), 0, 3)
else:
image = np.rollaxis(sample["image"][:3].numpy(), 0, 3)
image = percentile_normalization(image, axis=(0, 1))
ncols = 1
show_mask = "mask" in sample
show_predictions = "prediction" in sample
if show_mask:
mask = sample["mask"].numpy()
ncols += 1
if show_predictions:
prediction = sample["prediction"].numpy()
ncols += 1
fig, axs = plt.subplots(ncols=ncols, figsize=(ncols * 8, 8))
if not isinstance(axs, np.ndarray):
axs = [axs]
axs[0].imshow(image)
axs[0].axis("off")
if show_titles:
axs[0].set_title("Image")
if show_mask:
if self.speed_mask:
cmap = copy.copy(plt.get_cmap("autumn_r"))
cmap.set_under(color="black")
axs[1].imshow(mask, vmin=0.1, vmax=7, cmap=cmap, interpolation="none")
else:
axs[1].imshow(mask, cmap="Greys_r", interpolation="none")
axs[1].axis("off")
if show_titles:
axs[1].set_title("Label")
if show_predictions:
if self.speed_mask:
cmap = copy.copy(plt.get_cmap("autumn_r"))
cmap.set_under(color="black")
axs[2].imshow(
prediction, vmin=0.1, vmax=7, cmap=cmap, interpolation="none"
)
else:
axs[2].imshow(prediction, cmap="Greys_r", interpolation="none")
axs[2].axis("off")
if show_titles:
axs[2].set_title("Prediction")
if suptitle is not None:
plt.suptitle(suptitle)
return fig
[docs]class SpaceNet4(SpaceNet):
"""SpaceNet 4: Off-Nadir Buildings Dataset.
`SpaceNet 4 <https://spacenet.ai/off-nadir-building-detection/>`_ is a
dataset of 27 WV-2 imagery captured at varying off-nadir angles and
associated building footprints over the city of Atlanta. The off-nadir angle
ranges from 7 degrees to 54 degrees.
Dataset features:
* No. of chipped images: 28,728 (PAN/MS/PS-RGBNIR)
* No. of label files: 1064
* No. of building footprints: >120,000
* Area Coverage: 665 sq km
* Chip size: 225 x 225 (MS), 900 x 900 (PAN/PS-RGBNIR)
Dataset format:
* Imagery - Worldview-2 GeoTIFFs
* PAN.tif (Panchromatic)
* MS.tif (Multispectral)
* PS-RGBNIR (Pansharpened RGBNIR)
* Labels - GeoJSON
* labels.geojson
If you use this dataset in your research, please cite the following paper:
* https://arxiv.org/abs/1903.12239
"""
dataset_id = "spacenet4"
collection_md5_dict = {"sn4_AOI_6_Atlanta": "c597d639cba5257927a97e3eff07b753"}
imagery = {"MS": "MS.tif", "PAN": "PAN.tif", "PS-RGBNIR": "PS-RGBNIR.tif"}
chip_size = {"MS": (225, 225), "PAN": (900, 900), "PS-RGBNIR": (900, 900)}
label_glob = "labels.geojson"
angle_catalog_map = {
"nadir": [
"1030010003D22F00",
"10300100023BC100",
"1030010003993E00",
"1030010003CAF100",
"1030010002B7D800",
"10300100039AB000",
"1030010002649200",
"1030010003C92000",
"1030010003127500",
"103001000352C200",
"103001000307D800",
],
"off-nadir": [
"1030010003472200",
"1030010003315300",
"10300100036D5200",
"103001000392F600",
"1030010003697400",
"1030010003895500",
"1030010003832800",
],
"very-off-nadir": [
"10300100035D1B00",
"1030010003CCD700",
"1030010003713C00",
"10300100033C5200",
"1030010003492700",
"10300100039E6200",
"1030010003BDDC00",
"1030010003CD4300",
"1030010003193D00",
],
}
[docs] def __init__(
self,
root: str = "data",
image: str = "PS-RGBNIR",
angles: list[str] = [],
transforms: Optional[Callable[[dict[str, Any]], dict[str, Any]]] = None,
download: bool = False,
api_key: Optional[str] = None,
checksum: bool = False,
) -> None:
"""Initialize a new SpaceNet 4 Dataset instance.
Args:
root: root directory where dataset can be found
image: image selection which must be in ["MS", "PAN", "PS-RGBNIR"]
angles: angle selection which must be in ["nadir", "off-nadir",
"very-off-nadir"]
transforms: a function/transform that takes input sample and its target as
entry and returns a transformed version
download: if True, download dataset and store it in the root directory.
api_key: a RadiantEarth MLHub API key to use for downloading the dataset
checksum: if True, check the MD5 of the downloaded files (may be slow)
Raises:
RuntimeError: if ``download=False`` but dataset is missing
"""
collections = ["sn4_AOI_6_Atlanta"]
assert image in {"MS", "PAN", "PS-RGBNIR"}
self.angles = angles
if self.angles:
for angle in self.angles:
assert angle in self.angle_catalog_map.keys()
super().__init__(
root, image, collections, transforms, download, api_key, checksum
)
def _load_files(self, root: str) -> list[dict[str, str]]:
"""Return the paths of the files in the dataset.
Args:
root: root dir of dataset
Returns:
list of dicts containing paths for each pair of image and label
"""
files = []
nadir = []
offnadir = []
veryoffnadir = []
images = glob.glob(os.path.join(root, self.collections[0], "*", self.filename))
images = sorted(images)
catalog_id_pattern = re.compile(r"(_[A-Z0-9])\w+$")
for imgpath in images:
imgdir = os.path.basename(os.path.dirname(imgpath))
match = catalog_id_pattern.search(imgdir)
assert match is not None, "Invalid image directory"
catalog_id = match.group()[1:]
lbl_dir = os.path.dirname(imgpath).split("-nadir")[0]
lbl_path = os.path.join(f"{lbl_dir}-labels", self.label_glob)
assert os.path.exists(lbl_path)
_file = {"image_path": imgpath, "label_path": lbl_path}
if catalog_id in self.angle_catalog_map["very-off-nadir"]:
veryoffnadir.append(_file)
elif catalog_id in self.angle_catalog_map["off-nadir"]:
offnadir.append(_file)
elif catalog_id in self.angle_catalog_map["nadir"]:
nadir.append(_file)
angle_file_map = {
"nadir": nadir,
"off-nadir": offnadir,
"very-off-nadir": veryoffnadir,
}
if not self.angles:
files.extend(nadir + offnadir + veryoffnadir)
else:
for angle in self.angles:
files.extend(angle_file_map[angle])
return files
[docs]class SpaceNet5(SpaceNet3):
r"""SpaceNet 5: Automated Road Network Extraction and Route Travel Time Estimation.
`SpaceNet 5 <https://spacenet.ai/sn5-challenge/>`_
is a dataset of road networks over the cities of Moscow, Mumbai and San
Juan (unavailable).
Collection features:
+------------+---------------------+------------+---------------------------+
| AOI | Area (km\ :sup:`2`\)| # Images | # Road Network Labels (km)|
+============+=====================+============+===========================+
| Moscow | 1353 | 1353 | 3066 |
+------------+---------------------+------------+---------------------------+
| Mumbai | 1021 | 1016 | 1951 |
+------------+---------------------+------------+---------------------------+
Imagery features:
.. list-table::
:widths: 10 10 10 10 10
:header-rows: 1
:stub-columns: 1
* -
- PAN
- MS
- PS-MS
- PS-RGB
* - GSD (m)
- 0.31
- 1.24
- 0.30
- 0.30
* - Chip size (px)
- 1300 x 1300
- 325 x 325
- 1300 x 1300
- 1300 x 1300
Dataset format:
* Imagery - Worldview-3 GeoTIFFs
* PAN.tif (Panchromatic)
* MS.tif (Multispectral)
* PS-MS (Pansharpened Multispectral)
* PS-RGB (Pansharpened RGB)
* Labels - GeoJSON
* labels.geojson
If you use this dataset in your research, please use the following citation:
* The SpaceNet Partners, “SpaceNet5: Automated Road Network Extraction and
Route Travel Time Estimation from Satellite Imagery”,
https://spacenet.ai/sn5-challenge/
.. versionadded:: 0.2
"""
dataset_id = "spacenet5"
collection_md5_dict = {
"sn5_AOI_7_Moscow": "b18107f878152fe7e75444373c320cba",
"sn5_AOI_8_Mumbai": "1f1e2b3c26fbd15bfbcdbb6b02ae051c",
}
imagery = {
"MS": "MS.tif",
"PAN": "PAN.tif",
"PS-MS": "PS-MS.tif",
"PS-RGB": "PS-RGB.tif",
}
chip_size = {
"MS": (325, 325),
"PAN": (1300, 1300),
"PS-MS": (1300, 1300),
"PS-RGB": (1300, 1300),
}
label_glob = "labels.geojson"
[docs] def __init__(
self,
root: str = "data",
image: str = "PS-RGB",
speed_mask: Optional[bool] = False,
collections: list[str] = [],
transforms: Optional[Callable[[dict[str, Any]], dict[str, Any]]] = None,
download: bool = False,
api_key: Optional[str] = None,
checksum: bool = False,
) -> None:
"""Initialize a new SpaceNet 5 Dataset instance.
Args:
root: root directory where dataset can be found
image: image selection which must be in ["MS", "PAN", "PS-MS", "PS-RGB"]
speed_mask: use multi-class speed mask (created by binning roads at
10 mph increments) as label if true, else use binary mask
collections: collection selection which must be a subset of:
[sn5_AOI_7_Moscow, sn5_AOI_8_Mumbai]. If unspecified, all
collections will be used.
transforms: a function/transform that takes input sample and its target as
entry and returns a transformed version
download: if True, download dataset and store it in the root directory.
api_key: a RadiantEarth MLHub API key to use for downloading the dataset
checksum: if True, check the MD5 of the downloaded files (may be slow)
Raises:
RuntimeError: if ``download=False`` but dataset is missing
"""
super().__init__(
root,
image,
speed_mask,
collections,
transforms,
download,
api_key,
checksum,
)
[docs]class SpaceNet6(SpaceNet):
r"""SpaceNet 6: Multi-Sensor All-Weather Mapping.
`SpaceNet 6 <https://spacenet.ai/sn6-challenge/>`_ is a dataset
of optical and SAR imagery over the city of Rotterdam.
Collection features:
+------------+---------------------+------------+-----------------------------+
| AOI | Area (km\ :sup:`2`\)| # Images | # Building Footprint Labels |
+============+=====================+============+=============================+
| Rotterdam | 120 | 3401 | 48000 |
+------------+---------------------+------------+-----------------------------+
Imagery features:
.. list-table::
:widths: 10 10 10 10 10 10
:header-rows: 1
:stub-columns: 1
* -
- PAN
- RGBNIR
- PS-RGB
- PS-RGBNIR
- SAR-Intensity
* - GSD (m)
- 0.5
- 2.0
- 0.5
- 0.5
- 0.5
* - Chip size (px)
- 900 x 900
- 450 x 450
- 900 x 900
- 900 x 900
- 900 x 900
Dataset format:
* Imagery - GeoTIFFs from Worldview-2 (optical) and Capella Space (SAR)
* PAN.tif (Panchromatic)
* RGBNIR.tif (Multispectral)
* PS-RGB (Pansharpened RGB)
* PS-RGBNIR (Pansharpened RGBNIR)
* SAR-Intensity (SAR Intensity)
* Labels - GeoJSON
* labels.geojson
If you use this dataset in your research, please cite the following paper:
* https://arxiv.org/abs/2004.06500
.. note::
This dataset requires the following additional library to be installed:
* `radiant-mlhub <https://pypi.org/project/radiant-mlhub/>`_ to download the
imagery and labels from the Radiant Earth MLHub
.. versionadded:: 0.4
"""
dataset_id = "spacenet6"
collections = ["sn6_AOI_11_Rotterdam"]
# This is actually the metadata hash
collection_md5_dict = {"sn6_AOI_11_Rotterdam": "66f7312218fec67a1e0b3b02b22c95cc"}
imagery = {
"PAN": "PAN.tif",
"RGBNIR": "RGBNIR.tif",
"PS-RGB": "PS-RGB.tif",
"PS-RGBNIR": "PS-RGBNIR.tif",
"SAR-Intensity": "SAR-Intensity.tif",
}
chip_size = {
"PAN": (900, 900),
"RGBNIR": (450, 450),
"PS-RGB": (900, 900),
"PS-RGBNIR": (900, 900),
"SAR-Intensity": (900, 900),
}
label_glob = "labels.geojson"
[docs] def __init__(
self,
root: str = "data",
image: str = "PS-RGB",
transforms: Optional[Callable[[dict[str, Any]], dict[str, Any]]] = None,
download: bool = False,
api_key: Optional[str] = None,
) -> None:
"""Initialize a new SpaceNet 6 Dataset instance.
Args:
root: root directory where dataset can be found
image: image selection which must be in ["PAN", "RGBNIR",
"PS-RGB", "PS-RGBNIR", "SAR-Intensity"]
transforms: a function/transform that takes input sample and its target as
entry and returns a transformed version
download: if True, download dataset and store it in the root directory.
api_key: a RadiantEarth MLHub API key to use for downloading the dataset
Raises:
RuntimeError: if ``download=False`` but dataset is missing
"""
self.root = root
self.image = image # For testing
self.filename = self.imagery[image]
self.transforms = transforms
if download:
self.__download(api_key)
self.files = self._load_files(os.path.join(root, self.dataset_id))
def __download(self, api_key: Optional[str] = None) -> None:
"""Download the dataset and extract it.
Args:
api_key: a RadiantEarth MLHub API key to use for downloading the dataset
Raises:
RuntimeError: if download doesn't work correctly or checksums don't match
"""
if os.path.exists(
os.path.join(
self.root, self.dataset_id, self.collections[0], "collection.json"
)
):
print("Files already downloaded and verified")
return
download_radiant_mlhub_dataset(self.dataset_id, self.root, api_key)
[docs]class SpaceNet7(SpaceNet):
"""SpaceNet 7: Multi-Temporal Urban Development Challenge.
`SpaceNet 7 <https://spacenet.ai/sn7-challenge/>`_ is a dataset which
consist of medium resolution (4.0m) satellite imagery mosaics acquired from
Planet Labs’ Dove constellation between 2017 and 2020. It includes ≈ 24
images (one per month) covering > 100 unique geographies, and comprises >
40,000 km2 of imagery and exhaustive polygon labels of building footprints
therein, totaling over 11M individual annotations.
Dataset features:
* No. of train samples: 1423
* No. of test samples: 466
* No. of building footprints: 11,080,000
* Area Coverage: 41,000 sq km
* Chip size: 1023 x 1023
* GSD: ~4m
Dataset format:
* Imagery - Planet Dove GeoTIFF
* mosaic.tif
* Labels - GeoJSON
* labels.geojson
If you use this dataset in your research, please cite the following paper:
* https://arxiv.org/abs/2102.04420
.. versionadded:: 0.2
"""
dataset_id = "spacenet7"
collection_md5_dict = {
"sn7_train_source": "9f8cc109d744537d087bd6ff33132340",
"sn7_train_labels": "16f873e3f0f914d95a916fb39b5111b5",
"sn7_test_source": "e97914f58e962bba3e898f08a14f83b2",
}
imagery = {"img": "mosaic.tif"}
chip_size = {"img": (1023, 1023)}
label_glob = "labels.geojson"
[docs] def __init__(
self,
root: str = "data",
split: str = "train",
transforms: Optional[Callable[[dict[str, Any]], dict[str, Any]]] = None,
download: bool = False,
api_key: Optional[str] = None,
checksum: bool = False,
) -> None:
"""Initialize a new SpaceNet 7 Dataset instance.
Args:
root: root directory where dataset can be found
split: split selection which must be in ["train", "test"]
transforms: a function/transform that takes input sample and its target as
entry and returns a transformed version
download: if True, download dataset and store it in the root directory.
api_key: a RadiantEarth MLHub API key to use for downloading the dataset
checksum: if True, check the MD5 of the downloaded files (may be slow)
Raises:
RuntimeError: if ``download=False`` but dataset is missing
"""
self.root = root
self.split = split
self.filename = self.imagery["img"]
self.transforms = transforms
self.checksum = checksum
assert split in {"train", "test"}, "Invalid split"
if split == "test":
self.collections = ["sn7_test_source"]
else:
self.collections = ["sn7_train_source", "sn7_train_labels"]
to_be_downloaded = self._check_integrity()
if to_be_downloaded:
if not download:
raise RuntimeError(
f"Dataset not found in `root={self.root}` and `download=False`, "
"either specify a different `root` directory or use "
"`download=True` to automatically download the dataset."
)
else:
self._download(to_be_downloaded, api_key)
self.files = self._load_files(root)
def _load_files(self, root: str) -> list[dict[str, str]]:
"""Return the paths of the files in the dataset.
Args:
root: root dir of dataset
Returns:
list of dicts containing paths for images and labels (if train split)
"""
files = []
if self.split == "train":
imgs = sorted(
glob.glob(os.path.join(root, "sn7_train_source", "*", self.filename))
)
lbls = sorted(
glob.glob(os.path.join(root, "sn7_train_labels", "*", self.label_glob))
)
for img, lbl in zip(imgs, lbls):
files.append({"image_path": img, "label_path": lbl})
else:
imgs = sorted(
glob.glob(os.path.join(root, "sn7_test_source", "*", self.filename))
)
for img in imgs:
files.append({"image_path": img})
return files
[docs] def __getitem__(self, index: int) -> dict[str, Tensor]:
"""Return an index within the dataset.
Args:
index: index to return
Returns:
data at that index
"""
files = self.files[index]
img, tfm, raster_crs = self._load_image(files["image_path"])
h, w = img.shape[1:]
ch, cw = self.chip_size["img"]
sample = {"image": img[:, :ch, :cw]}
if self.split == "train":
mask = self._load_mask(files["label_path"], tfm, raster_crs, (h, w))
sample["mask"] = mask[:ch, :cw]
if self.transforms is not None:
sample = self.transforms(sample)
return sample
