from pystac import Extensions
from pystac.item import Item
from pystac.extensions.base import (ItemExtension, ExtensionDefinition, ExtendedObject)
[docs]class ProjectionItemExt(ItemExtension):
"""ProjectionItemExt is the extension of an Item in the Projection Extension.
The Projection extension adds projection information to STAC Items.
Args:
item (Item): The item to be extended.
Attributes:
item (Item): The Item that is being extended.
Note:
Using ProjectionItemExt to directly wrap an item will add the 'proj' extension ID to
the item's stac_extensions.
"""
def __init__(self, item):
if item.stac_extensions is None:
item.stac_extensions = [Extensions.PROJECTION]
elif Extensions.PROJECTION not in item.stac_extensions:
item.stac_extensions.append(Extensions.PROJECTION)
self.item = item
[docs] def apply(self,
epsg,
wkt2=None,
projjson=None,
geometry=None,
bbox=None,
centroid=None,
shape=None,
transform=None):
"""Applies Projection extension properties to the extended Item.
Args:
epsg (int or None): REQUIRED. EPSG code of the datasource.
wkt2 (str or None): WKT2 string representing the Coordinate Reference System (CRS) that
the ``geometry`` and ``bbox`` fields represent
projjson (dict or None): PROJJSON dict representing the
Coordinate Reference System (CRS) that the ``geometry`` and ``bbox``
fields represent
geometry (dict or None): GeoJSON Polygon dict that defines the footprint of this Item.
bbox (List[float] or None): Bounding box of the Item in the asset CRS in
2 or 3 dimensions.
centroid (dict or None): A dict with members 'lat' and 'lon' that defines
coordinates representing the centroid of the item in the asset data CRS.
Coordinates are defined in latitude and longitude, even if the data coordinate
system may not use lat/long.
shape (List[int] or None): Number of pixels in Y and X directions for the default grid.
transform (List[float] or None): The affine transformation coefficients for the
default grid
"""
self.epsg = epsg
self.wkt2 = wkt2
self.projjson = projjson
self.geometry = geometry
self.bbox = bbox
self.centroid = centroid
self.shape = shape
self.transform = transform
@property
def epsg(self):
"""Get or sets the EPSG code of the datasource.
A Coordinate Reference System (CRS) is the data reference system (sometimes called a
'projection') used by the asset data, and can usually be referenced using an
`EPSG code <http://epsg.io/>`_.
If the asset data does not have a CRS, such as in the case of non-rectified imagery with
Ground Control Points, epsg should be set to None.
It should also be set to null if a CRS exists, but for which there is no valid EPSG code.
Returns:
int
"""
return self.get_epsg()
@epsg.setter
def epsg(self, v):
self.set_epsg(v)
[docs] def get_epsg(self, asset=None):
"""Gets an Item or an Asset epsg.
If an Asset is supplied and the Item property exists on the Asset,
returns the Asset's value. Otherwise returns the Item's value
Returns:
int
"""
if asset is None or 'proj:epsg' not in asset.properties:
return self.item.properties.get('proj:epsg')
else:
return asset.properties.get('proj:epsg')
[docs] def set_epsg(self, epsg, asset=None):
"""Set an Item or an Asset epsg.
If an Asset is supplied, sets the property on the Asset.
Otherwise sets the Item's value.
"""
self._set_property('proj:epsg', epsg, asset)
@property
def wkt2(self):
"""Get or sets the WKT2 string representing the Coordinate Reference System (CRS)
that the proj:geometry and proj:bbox fields represent
This value is a `WKT2 string <http://docs.opengeospatial.org/is/12-063r5/12-063r5.html>`_.
If the data does not have a CRS, such as in the case of non-rectified imagery with Ground
Control Points, wkt2 should be set to null. It should also be set to null if a CRS exists,
but for which a WKT2 string does not exist.
Returns:
str
"""
return self.get_wkt2()
@wkt2.setter
def wkt2(self, v):
self.set_wkt2(v)
[docs] def get_wkt2(self, asset=None):
"""Gets an Item or an Asset wkt2.
If an Asset is supplied and the Item property exists on the Asset,
returns the Asset's value. Otherwise returns the Item's value
Returns:
str
"""
if asset is None or 'proj:wkt2' not in asset.properties:
return self.item.properties.get('proj:wkt2')
else:
return asset.properties.get('proj:wkt2')
[docs] def set_wkt2(self, wkt2, asset=None):
"""Set an Item or an Asset wkt2.
If an Asset is supplied, sets the property on the Asset.
Otherwise sets the Item's value.
"""
self._set_property('proj:wkt2', wkt2, asset)
@property
def projjson(self):
"""Get or sets the PROJJSON string representing the Coordinate Reference System (CRS)
that the proj:geometry and proj:bbox fields represent
This value is a `PROJJSON object <https://proj.org/specifications/projjson.html>`_.
If the data does not have a CRS, such as in the case of non-rectified imagery with Ground
Control Points, projjson should be set to null. It should also be set to null if a
CRS exists, but for which a PROJJSON string does not exist.
The schema for this object can be found
`here <https://proj.org/schemas/v0.2/projjson.schema.json>`_.
Returns:
dict
"""
return self.get_projjson()
@projjson.setter
def projjson(self, v):
self.set_projjson(v)
[docs] def get_projjson(self, asset=None):
"""Gets an Item or an Asset projjson.
If an Asset is supplied and the Item property exists on the Asset,
returns the Asset's value. Otherwise returns the Item's value
Returns:
dict
"""
if asset is None or 'proj:projjson' not in asset.properties:
return self.item.properties.get('proj:projjson')
else:
return asset.properties.get('proj:projjson')
[docs] def set_projjson(self, projjson, asset=None):
"""Set an Item or an Asset projjson.
If an Asset is supplied, sets the property on the Asset.
Otherwise sets the Item's value.
"""
self._set_property('proj:projjson', projjson, asset)
@property
def geometry(self):
"""Get or sets a Polygon GeoJSON dict representing the footprint of this item.
This dict should be formatted according the Polygon object format specified in
`RFC 7946, sections 3.1.6 <https://tools.ietf.org/html/rfc7946>`_,
except not necessarily in EPSG:4326 as required by RFC7946. Specified based on the
``epsg``, ``projjson`` or ``wkt2`` fields (not necessarily EPSG:4326).
Ideally, this will be represented by a Polygon with five coordinates, as the item in
the asset data CRS should be a square aligned to the original CRS grid.
Returns:
dict
"""
return self.get_geometry()
@geometry.setter
def geometry(self, v):
self.set_geometry(v)
[docs] def get_geometry(self, asset=None):
"""Gets an Item or an Asset projection geometry.
If an Asset is supplied and the Item property exists on the Asset,
returns the Asset's value. Otherwise returns the Item's value
Returns:
dict
"""
if asset is None or 'proj:geometry' not in asset.properties:
return self.item.properties.get('proj:geometry')
else:
return asset.properties.get('proj:geometry')
[docs] def set_geometry(self, geometry, asset=None):
"""Set an Item or an Asset projection geometry.
If an Asset is supplied, sets the property on the Asset.
Otherwise sets the Item's value.
"""
self._set_property('proj:geometry', geometry, asset)
@property
def bbox(self):
"""Get or sets the bounding box of the assets represented by this item in the asset
data CRS.
Specified as 4 or 6 coordinates based on the CRS defined in the ``epsg``, ``projjson``
or ``wkt2`` properties. First two numbers are coordinates of the lower left corner,
followed by coordinates of upper right corner, e.g.,
[west, south, east, north], [xmin, ymin, xmax, ymax], [left, down, right, up],
or [west, south, lowest, east, north, highest]. The length of the array must be 2*n
where n is the number of dimensions.
Returns:
List[float]
"""
return self.get_bbox()
@bbox.setter
def bbox(self, v):
self.set_bbox(v)
[docs] def get_bbox(self, asset=None):
"""Gets an Item or an Asset projection bbox.
If an Asset is supplied and the Item property exists on the Asset,
returns the Asset's value. Otherwise returns the Item's value
Returns:
List[float]
"""
if asset is None or 'proj:bbox' not in asset.properties:
return self.item.properties.get('proj:bbox')
else:
return asset.properties.get('proj:bbox')
[docs] def set_bbox(self, bbox, asset=None):
"""Set an Item or an Asset projection bbox.
If an Asset is supplied, sets the property on the Asset.
Otherwise sets the Item's value.
"""
self._set_property('proj:bbox', bbox, asset)
@property
def centroid(self):
"""Get or sets coordinates representing the centroid of the item in the asset data CRS.
Coordinates are defined in latitude and longitude, even if the data coordinate system
does not use lat/long.
Exmample::
item.ext.proj.centroid = { 'lat': 0.0, 'lon': 0.0 }
Returns:
dict
"""
return self.get_centroid()
@centroid.setter
def centroid(self, v):
self.set_centroid(v)
[docs] def get_centroid(self, asset=None):
"""Gets an Item or an Asset centroid.
If an Asset is supplied and the Item property exists on the Asset,
returns the Asset's value. Otherwise returns the Item's value
Returns:
dict
"""
if asset is None or 'proj:centroid' not in asset.properties:
return self.item.properties.get('proj:centroid')
else:
return asset.properties.get('proj:centroid')
[docs] def set_centroid(self, centroid, asset=None):
"""Set an Item or an Asset centroid.
If an Asset is supplied, sets the property on the Asset.
Otherwise sets the Item's value.
"""
self._set_property('proj:centroid', centroid, asset)
@property
def shape(self):
"""Get or sets the number of pixels in Y and X directions for the default grid.
The shape is an array of integers that represents the number of pixels in the most
common pixel grid used by the item's assets. The number of pixels should be specified
in Y, X order. If the shape is defined in an item's properties it is used as the default
shape for all assets that don't have an overriding shape.
Returns:
List[int]
"""
return self.get_shape()
@shape.setter
def shape(self, v):
self.set_shape(v)
[docs] def get_shape(self, asset=None):
"""Gets an Item or an Asset shape.
If an Asset is supplied and the Item property exists on the Asset,
returns the Asset's value. Otherwise returns the Item's value
Returns:
List[int]
"""
if asset is None or 'proj:shape' not in asset.properties:
return self.item.properties.get('proj:shape')
else:
return asset.properties.get('proj:shape')
[docs] def set_shape(self, shape, asset=None):
"""Set an Item or an Asset shape.
If an Asset is supplied, sets the property on the Asset.
Otherwise sets the Item's value.
"""
self._set_property('proj:shape', shape, asset)
@property
def transform(self):
"""Get or sets the the affine transformation coefficients for the default grid.
The transform is a linear mapping from pixel coordinate space (Pixel, Line) to
projection coordinate space (Xp, Yp). It is a 3x3 matrix stored as a flat array of 9
elements in row major order. Since the last row is always 0,0,1 it can be omitted, in
which case only 6 elements are recorded. This mapping can be obtained from
GDAL `GetGeoTransform <https://gdal.org/api/gdaldataset_cpp.html#_CPPv4N11GDALDataset15GetGeoTransformEPd>`_
or the
Rasterio `Transform <https://rasterio.readthedocs.io/en/stable/api/rasterio.io.html#rasterio.io.BufferedDatasetWriter.transform>`_.
Returns:
List[float]
""" # noqa: E501
return self.get_transform()
@transform.setter
def transform(self, v):
self.set_transform(v)
@classmethod
def _object_links(cls):
return []
[docs] @classmethod
def from_item(cls, item):
return cls(item)
PROJECTION_EXTENSION_DEFINITION = ExtensionDefinition(Extensions.PROJECTION,
[ExtendedObject(Item, ProjectionItemExt)])