PlaProduct#

class PlaProduct(product_path: Union[str, CloudPath, Path], archive_path: Optional[Union[str, CloudPath, Path]] = None, output_path: Optional[Union[str, CloudPath, Path]] = None, remove_tmp: bool = False, **kwargs)[source]#

Bases: PlanetProduct

Class of PlanetScope products. See here for more information.

The scaling factor to retrieve the calibrated radiance is 0.01.

__init__(product_path: Union[str, CloudPath, Path], archive_path: Optional[Union[str, CloudPath, Path]] = None, output_path: Optional[Union[str, CloudPath, Path]] = None, remove_tmp: bool = False, **kwargs) None#
clean_tmp()#

Clean the temporary directory of the current product

clear()#

Clear this product’s cache

get_band_paths(band_list: list, resolution: Optional[float] = None, **kwargs) dict[source]#

Return the paths of required bands.

>>> from eoreader.reader import Reader
>>> from eoreader.bands import *
>>> path = r"SENTINEL2A_20190625-105728-756_L2A_T31UEQ_C_V2-2"
>>> prod = Reader().open(path)
>>> prod.get_band_paths([GREEN, RED])
{
    <SpectralBandNames.GREEN: 'GREEN'>:
    'SENTINEL2A_20190625-105728-756_L2A_T31UEQ_C_V2-2/SENTINEL2A_20190625-105728-756_L2A_T31UEQ_C_V2-2_FRE_B3.tif',
    <SpectralBandNames.RED: 'RED'>:
    'SENTINEL2A_20190625-105728-756_L2A_T31UEQ_C_V2-2/SENTINEL2A_20190625-105728-756_L2A_T31UEQ_C_V2-2_FRE_B4.tif'
}
Parameters
  • band_list (list) – List of the wanted bands

  • resolution (float) – Band resolution

  • kwargs – Other arguments used to load bands

Returns

Dictionary containing the path of each queried band

Return type

dict

get_date(as_date: bool = False) Union[str, date]#

Get the product’s acquisition date.

>>> from eoreader.reader import Reader
>>> path = r"S2A_MSIL1C_20200824T110631_N0209_R137_T30TTK_20200824T150432.SAFE.zip"
>>> prod = Reader().open(path)
>>> prod.get_date(as_date=True)
datetime.datetime(2020, 8, 24, 0, 0)
>>> prod.get_date(as_date=False)
'20200824'
Parameters

as_date (bool) – Return the date as a datetime.date. If false, returns a string.

Returns

Its acquisition date

Return type

str

get_datetime(as_datetime: bool = False) Union[str, datetime][source]#

Get the product’s acquisition datetime, with format YYYYMMDDTHHMMSS <-> %Y%m%dT%H%M%S

>>> from eoreader.reader import Reader
>>> path = r"SENTINEL2A_20190625-105728-756_L2A_T31UEQ_C_V2-2"
>>> prod = Reader().open(path)
>>> prod.get_datetime(as_datetime=True)
datetime.datetime(2019, 6, 25, 10, 57, 28, 756000), fetched from metadata, so we have the ms
>>> prod.get_datetime(as_datetime=False)
'20190625T105728'
Parameters

as_datetime (bool) – Return the date as a datetime.datetime. If false, returns a string.

Returns

Its acquisition datetime

Return type

Union[str, datetime]

get_default_band() BandNames#

Get default band: GREEN for optical data as every optical satellite has a GREEN band.

>>> from eoreader.reader import Reader
>>> path = r"S2A_MSIL1C_20200824T110631_N0209_R137_T30TTK_20200824T150432.SAFE.zip"
>>> prod = Reader().open(path)
>>> prod.get_default_band()
<SpectralBandNames.GREEN: 'GREEN'>
Returns

Default band

Return type

str

get_default_band_path(**kwargs) Union[CloudPath, Path]#

Get default band (GREEN for optical data) path.

>>> from eoreader.reader import Reader
>>> path = r"S2A_MSIL1C_20200824T110631_N0209_R137_T30TTK_20200824T150432.SAFE.zip"
>>> prod = Reader().open(path)
>>> prod.get_default_band_path()
'zip+file://S2A_MSIL1C_20200824T110631_N0209_R137_T30TTK_20200824T150432.SAFE.zip!/S2A_MSIL1C_20200824T110631_N0209_R137_T30TTK_20200824T150432.SAFE/GRANULE/L1C_T30TTK_A027018_20200824T111345/IMG_DATA/T30TTK_20200824T110631_B03.jp2'
Parameters

kwargs – Additional arguments

Returns

Default band path

Return type

Union[CloudPath, Path]

get_existing_band_paths() dict#

Return the existing band paths (orthorectified if needed).

>>> from eoreader.reader import Reader
>>> path = r"S2A_MSIL1C_20200824T110631_N0209_R137_T30TTK_20200824T150432.SAFE.zip"
>>> prod = Reader().open(path)
>>> prod.get_existing_band_paths()
{
    <SpectralBandNames.CA: 'COASTAL_AEROSOL'>: 'zip+file://S2A_MSIL1C_20200824T110631_N0209_R137_T30TTK_20200824T150432.SAFE.zip!/S2A_MSIL1C_20200824T110631_N0209_R137_T30TTK_20200824T150432.SAFE/GRANULE/L1C_T30TTK_A027018_20200824T111345/IMG_DATA/T30TTK_20200824T110631_B01.jp2',
    ...,
    <SpectralBandNames.SWIR_2: 'SWIR_2'>: 'zip+file://S2A_MSIL1C_20200824T110631_N0209_R137_T30TTK_20200824T150432.SAFE.zip!/S2A_MSIL1C_20200824T110631_N0209_R137_T30TTK_20200824T150432.SAFE/GRANULE/L1C_T30TTK_A027018_20200824T111345/IMG_DATA/T30TTK_20200824T110631_B12.jp2'
}
Returns

Dictionary containing the path of each queried band

Return type

dict

get_existing_bands() list#

Return the existing band.

>>> from eoreader.reader import Reader
>>> path = r"S2A_MSIL1C_20200824T110631_N0209_R137_T30TTK_20200824T150432.SAFE.zip"
>>> prod = Reader().open(path)
>>> prod.get_existing_bands()
[<SpectralBandNames.CA: 'COASTAL_AEROSOL'>,
<SpectralBandNames.BLUE: 'BLUE'>,
<SpectralBandNames.GREEN: 'GREEN'>,
<SpectralBandNames.RED: 'RED'>,
<SpectralBandNames.VRE_1: 'VEGETATION_RED_EDGE_1'>,
<SpectralBandNames.VRE_2: 'VEGETATION_RED_EDGE_2'>,
<SpectralBandNames.VRE_3: 'VEGETATION_RED_EDGE_3'>,
<SpectralBandNames.NIR: 'NIR'>,
<SpectralBandNames.NNIR: 'NARROW_NIR'>,
<SpectralBandNames.WV: 'WATER_VAPOUR'>,
<SpectralBandNames.CIRRUS: 'CIRRUS'>,
<SpectralBandNames.SWIR_1: 'SWIR_1'>,
<SpectralBandNames.SWIR_2: 'SWIR_2'>]
Returns

List of existing bands in the products

Return type

list

get_quicklook_path() Union[None, str]#

Get quicklook path if existing (no such thing for Sentinel-2)

Returns

Quicklook path

Return type

str

get_raw_band_paths(**kwargs) dict#

Return the raw band paths.

Parameters

kwargs – Additional arguments

Returns

Dictionary containing the path of each queried band

Return type

dict

has_band(band: Union[BandNames, Callable]) bool#

Does this product has the specified band ?

By band, we mean:

  • satellite band

  • index

  • DEM band

  • cloud band

>>> from eoreader.reader import Reader
>>> from eoreader.bands import *
>>> path = r"S2A_MSIL1C_20200824T110631_N0209_R137_T30TTK_20200824T150432.SAFE.zip"
>>> prod = Reader().open(path)
>>> prod.has_band(GREEN)
True
>>> prod.has_band(TIR_2)
False
>>> prod.has_band(NDVI)
True
>>> prod.has_band(SHADOWS)
False
>>> prod.has_band(HILLSHADE)
True
Parameters

band (Union[BandNames, Callable]) – EOReader band (optical, SAR, clouds, DEM)

Returns

True if the products has the specified band

Return type

bool

has_bands(bands: Union[list, BandNames, Callable]) bool#

Does this product has the specified bands ?

By band, we mean:

  • satellite band

  • index

  • DEM band

  • cloud band

See has_bands for a code example.

Parameters

bands (Union[list, BandNames, Callable]) – EOReader bands (optical, SAR, clouds, DEM)

Returns

True if the products has the specified band

Return type

bool

load(bands: Union[list, BandNames, Callable], resolution: Optional[float] = None, size: Optional[Union[list, tuple]] = None, **kwargs) dict#

Open the bands and compute the wanted index.

  • For Optical data:

    The bands will be purged of nodata and invalid pixels (if specified with the CLEAN_OPTICAL keyword), the nodata will be set to -9999 and the bands will be DataArrays in float32.

  • For SAR data:

    The bands will be purged of nodata (not over the sea), the nodata will be set to 0 to respect SNAP’s behavior and the bands will be DataArray in float32.

Bands that come out this function at the same time are collocated and therefore have the same shapes. This can be broken if you load data separately. Its is best to always load DEM data with some real bands.

>>> from eoreader.reader import Reader
>>> from eoreader.bands import *
>>> path = r"S2A_MSIL1C_20200824T110631_N0209_R137_T30TTK_20200824T150432.SAFE.zip"
>>> prod = Reader().open(path)
>>> bands = prod.load([GREEN, NDVI], resolution=20)
Parameters
  • bands (Union[list, BandNames, Callable]) – Band list

  • resolution (float) – Resolution of the band, in meters

  • size (Union[tuple, list]) – Size of the array (width, height). Not used if resolution is provided.

  • kwargs – Other arguments used to load bands

Returns

{band_name, band xarray}

Return type

dict

open_mask(mask_id: str, resolution: Optional[float] = None, size: Optional[Union[list, tuple]] = None) Optional[DataArray]#

Open a Planet UDM2 (Usable Data Mask) mask, band by band, as a xarray. Returns None if the mask is not available.

Do not open cloud mask with this function. Use load instead.

See UDM2 specifications for more information.

Accepted mask IDs:

  • CLEAR: Band 1 Clear map [0, 1] 0: not clear, 1: clear

  • SNOW: Band 2 Snow map [0, 1] 0: no snow or ice, 1: snow or ice

  • SHADOW: Band 3 Shadow map [0, 1] 0: no shadow, 1: shadow

  • LIGHT_HAZE: Band 4 Light haze map [0, 1] 0: no light haze, 1: light haze

  • HEAVY_HAZE: Band 5 Heavy haze map [0, 1] 0: no heavy haze, 1: heavy haze

  • CLOUD: Band 6 Cloud map [0, 1] 0: no cloud, 1: cloud

  • CONFIDENCE: Band 7 Confidence map [0-100] %age value: per-pixel algorithmic confidence in classif

  • UNUSABLE: Band 8 Unusable pixels – Equivalent to the UDM asset

>>> from eoreader.bands import *
>>> from eoreader.reader import Reader
>>> path = r"SENTINEL2B_20190401-105726-885_L2A_T31UEQ_D_V2-0.zip"
>>> prod = Reader().open(path)
>>> prod.open_mask("EDG", GREEN)
array([[[0, ..., 0]]], dtype=uint8)
Parameters
  • mask_id – Mask ID

  • resolution (float) – Band resolution in meters

  • size (Union[tuple, list]) – Size of the array (width, height). Not used if resolution is provided.

Returns

Mask array

Return type

Union[xarray.DataArray, None]

open_mask_udm(resolution: Optional[float] = None, size: Optional[Union[list, tuple]] = None) Optional[DataArray]#

Open a Planet UDM (Unusable Data Mask) mask as a xarray. For RapidEye, the mask is subsampled to 50m, so this function will interpolate to make it to the correct resolution Returns None if the mask is not available.

Do not open cloud mask with this function. Use load instead.

See here for more information.

Parameters
  • resolution (float) – Band resolution in meters

  • size (Union[tuple, list]) – Size of the array (width, height). Not used if resolution is provided.

Returns

Mask array

Return type

Union[xarray.DataArray, None]

open_mask_udm2(mask_id: str, resolution: Optional[float] = None, size: Optional[Union[list, tuple]] = None) Optional[DataArray]#

Open a Planet UDM2 (Usable Data Mask) mask, band by band, as a xarray. Returns None if the mask is not available.

Do not open cloud mask with this function. Use load instead.

See UDM2 specifications for more information.

Accepted mask IDs:

  • CLEAR: Band 1 Clear map [0, 1] 0: not clear, 1: clear

  • SNOW: Band 2 Snow map [0, 1] 0: no snow or ice, 1: snow or ice

  • SHADOW: Band 3 Shadow map [0, 1] 0: no shadow, 1: shadow

  • LIGHT_HAZE: Band 4 Light haze map [0, 1] 0: no light haze, 1: light haze

  • HEAVY_HAZE: Band 5 Heavy haze map [0, 1] 0: no heavy haze, 1: heavy haze

  • CLOUD: Band 6 Cloud map [0, 1] 0: no cloud, 1: cloud

  • CONFIDENCE: Band 7 Confidence map [0-100] %age value: per-pixel algorithmic confidence in classif

  • UNUSABLE: Band 8 Unusable pixels – Equivalent to the UDM asset

>>> from eoreader.bands import *
>>> from eoreader.reader import Reader
>>> path = r"SENTINEL2B_20190401-105726-885_L2A_T31UEQ_D_V2-0.zip"
>>> prod = Reader().open(path)
>>> prod.open_mask("EDG", GREEN)
array([[[0, ..., 0]]], dtype=uint8)
Parameters
  • mask_id – Mask ID

  • resolution (float) – Band resolution in meters

  • size (Union[tuple, list]) – Size of the array (width, height). Not used if resolution is provided.

Returns

Mask array

Return type

Union[xarray.DataArray, None]

plot() None#

Plot the quicklook if existing

read_mtd()#

Read metadata and outputs the metadata XML root and its namespaces as a dict.

>>> from eoreader.reader import Reader
>>> path = r"S1A_IW_GRDH_1SDV_20191215T060906_20191215T060931_030355_0378F7_3696.zip"
>>> prod = Reader().open(path)
>>> prod.read_mtd()
(<Element product at 0x1832895d788>, '')
Returns

Metadata XML root and its namespace

Return type

(etree._Element, dict)

stack(bands: list, resolution: Optional[float] = None, size: Optional[Union[list, tuple]] = None, stack_path: Optional[Union[str, CloudPath, Path]] = None, save_as_int: bool = False, **kwargs) DataArray#

Stack bands and index of a products.

>>> from eoreader.reader import Reader
>>> from eoreader.bands import *
>>> path = r"S2A_MSIL1C_20200824T110631_N0209_R137_T30TTK_20200824T150432.SAFE.zip"
>>> prod = Reader().open(path)
>>> stack = prod.stack([NDVI, MNDWI, GREEN], resolution=20)  # In meters
Parameters
  • bands (list) – Bands and index combination

  • resolution (float) – Stack resolution. . If not specified, use the product resolution.

  • size (Union[tuple, list]) – Size of the array (width, height). Not used if resolution is provided.

  • stack_path (Union[str, CloudPath, Path]) – Stack path

  • save_as_int (bool) – Convert stack to uint16 to save disk space (and therefore multiply the values by 10.000)

  • **kwargs – Other arguments passed to load or rioxarray.to_raster() (such as compress)

Returns

Stack as a DataArray

Return type

xr.DataArray

archive_path#

Archive path, same as the product path if not specified. Useful when you want to know where both the extracted and archived version of your product are stored.

bands#

Band mapping between band wrapping names such as GREEN and band real number such as 03 for Sentinel-2.

condensed_name#

Condensed name, the filename with only useful data to keep the name unique (ie. 20191215T110441_S2_30TXP_L2A_122756). Used to shorten names and paths.

constellation#

Product constellation, such as Sentinel-2

constellation_id#

Constellation ID, i.e. S2 for Sentinel-2

corresponding_ref#

The corresponding reference products to the current one (if the product is not a reference but has a reference data corresponding to it). A list because of multiple ref in case of non-stackable products (S3, S1…)

crs = <methodtools._LruCacheWire object>#
date#

Acquisition date.

datetime#

Acquisition datetime.

default_transform = <methodtools._LruCacheWire object>#
extent = <methodtools._LruCacheWire object>#
filename#

Product filename

footprint = <methodtools._LruCacheWire object>#
get_cloud_cover = <methodtools._LruCacheWire object>#
get_mean_sun_angles = <methodtools._LruCacheWire object>#
get_mean_viewing_angles = <methodtools._LruCacheWire object>#
get_orbit_direction = <methodtools._LruCacheWire object>#
instrument#

Product instrument, such as MSI for Sentinel-2 data.

is_archived#

Is the archived product is processed (a products is considered as archived if its products path is a directory).

is_ortho#

True if the images are orthorectified and the footprint is retrieved easily.

is_reference#

If the product is a reference, used for algorithms that need pre and post data, such as fire detection.

name#

Product true name (as specified in the metadata)

needs_extraction#

Does this product needs to be extracted to be processed ? (True by default).

nodata#

Product nodata, set to -9999 by default

property output: Union[cloudpathlib.cloudpath.CloudPath, Path]#

Output directory of the product, to write orthorectified data for example.

path#

Usable path to the product, either extracted or archived path, according to the satellite.

product_type#

Product type, satellite-related field, such as L1C or L2A for Sentinel-2 data.

resolution#

Default resolution in meters of the current product. For SAR product, we use Ground Range resolution as we will automatically orthorectify the tiles.

sensor_type#

Sensor type, SAR or optical.

split_name#

Split name, to retrieve every information from its true name (dates, tile, product type…).

property stac: StacItem#
tile_name#

Tile if possible (for data that can be piled, for example S2 and Landsats).