GeoContext¶
Datasets in the Descartes Labs catalog have many different resolutions and
projections. In two different images, even covering the same place on Earth,
the pixels (i, j)
usually correspond to two different points on the ground.
GeoContexts are a way to ensure multiple images from different sources are spatially compatible—that is, they all have the same shape (same width and height, in pixels), and the same pixel in each image corresponds to the same area on Earth.
They do this by simply capturing all the spatial parameters that affect how imagery is rasterized—namely output resolution, coordinate reference system, and bounding box—in one object that can be passed into different method calls. In typical use, these contexts are created for you with reasonable defaults, so you only need to understand the different parameters when you need more control.
The different subclasses of GeoContext
implement different
functionality.
AOI
clips to arbitrary geometry, and lets you specify any output resolution and projection.DLTile
helps you split large regions up into a grid of any spacing and resolution, and represents a single tile in that grid, in UTM projection.

class
AOI
(geometry=None, resolution=None, crs=None, align_pixels=True, bounds=None, bounds_crs='EPSG:4326', shape=None)[source]¶ A GeoContext that clips scenes to a geometry, and/or to square bounds, with any output resolution and CRS.
Examples
cutline_aoi = dl.scenes.AOI(my_geometry, resolution=40) aoi_with_cutline_disabled = cutline_aoi.assign(geometry=None) no_cutline_aoi = dl.scenes.AOI(geometry=None, resolution=15, bounds=(40, 35, 39, 36)) aoi_without_auto_bounds = dl.scenes.AOI(geometry=my_geometry, resolution=15, bounds=(40, 35, 39, 36)) aoi_with_specific_pixel_dimensions = dl.scenes.AOI(geometry=my_geometry, shape=(200, 400))
 Parameters
geometry (GeoJSONlike dict, object with
__geo_interface__
; optional) – Clip scenes to this geometry. Coordinates must be WGS84 (latlon). IfNone
, scenes will just be clipped tobounds
.resolution (float, optional) – Distance, in units of the CRS, that the edge of each pixel represents on the ground. Can only specify one of
resolution
andshape
.crs (str, optional) – Coordinate Reference System into which scenes will be projected, expressed as an EPSG code (like
EPSG:4326
), a PROJ.4 definition, or an OGC CRS WellKnown Text stringalign_pixels (bool, optional, default True) –
If True, this ensures that, in different Scenes rasterized with this same AOI GeoContext, pixels
(i, j)
correspond to the same area in space. This is accomplished by snapping the coordinates of the origin (topleft corner of topleft pixel) to a nonfractional interval ofresolution
.If
align_pixels
is False, when using scenes with different native resolutions and/or projections, pixels at the same indicies can be misaligned by a fraction ofresolution
(i.e. correspond to slighly different coordinates in space).However, this requires warping of the original image, which can be undesireable when you want to work with the original data in its native resolution and projection.
bounds (4tuple, optional) – Clip scenes to these
(min_x, min_y, max_x, max_y)
bounds, expressed inbounds_crs
(which defaults to WGS84 latlon).bounds
are automatically computed fromgeometry
if not specified. Otherwise,bounds
are required.bounds_crs (str, optional, default "EPSG:4326") – The Coordinate Reference System of the
bounds
, given as an EPSG code (likeEPSG:4326
), a PROJ.4 definition, or an OGC CRS WellKnown Text string.shape (2tuple, optional) –
(rows, columns)
, in pixels, the output raster should fit within; the longer side of the raster will be min(shape). Can only specify one ofresolution
andshape
.

assign
(geometry='unchanged', resolution='unchanged', crs='unchanged', align_pixels='unchanged', bounds='unchanged', bounds_crs='unchanged', shape='unchanged')[source]¶ Return a copy of the AOI with the given values assigned.
Note
If you are assigning a new geometry and want bounds to updated as well, use
bounds="update"
. This will also changebounds_crs
toEPSG:4326
, since the geometry’s coordinates are in WGS84 decimal degrees, so the new bounds determined from those coordinates must be in that CRS as well.If you assign
geometry
without changingbounds
, the new AOI GeoContext will produce rasters with the same shape and covering the same spatial area as the old one, just with pixels masked out that fall outside your new geometry. Returns
new
 Return type

property
align_pixels
¶ If True, this ensures that, in different
Scenes
rasterized with this same AOI GeoContext, pixels(i, j)
correspond to the same area in space. This is accomplished by snapping the coordinates of the origin (topleft corner of topleft pixel) to a nonfractional interval ofresolution
.If
align_pixels
is False, when using scenes with different native resolutions and/or projections, pixels at the same indicies can be misaligned by a fraction ofresolution
(i.e. correspond to slighly different coordinates in space).However, this requires warping of the original image, which can be undesireable when you want to work with the original data in its native resolution and projection.
 Type
bool

property
bounds
¶ Clip scenes to these
(min_x, min_y, max_x, max_y)
bounds, expressed in the coordinate reference system inbounds_crs
. Type
tuple

property
bounds_crs
¶ The coordinate reference system of the
bounds
, given as an EPSG code (likeEPSG:4326
), a PROJ.4 definition, or an OGC CRS WellKnown Text string. Type
str

property
crs
¶ Coordinate reference system into which scenes will be projected, expressed as an EPSG code (like
EPSG:4326
), a PROJ.4 definition, or an OGC CRS WellKnown Text string. Type
str

property
geometry
¶ Clip scenes to this geometry Coordinates must be WGS84 (latlon). If
None
, scenes will just be clipped tobounds
. Type
shapely geometry

property
raster_params
¶ The properties of this
AOI
, as keyword arguments to use forndarray
orraster
.Raises ValueError if
bounds
,crs
,bounds_crs
,resolution
, oralign_pixels
isNone
. Type
dict

property
resolution
¶ Distance, in units of the CRS, that the edge of each pixel represents on the ground.
 Type
float

property
shape
¶ (rows, columns)
, in pixels, the output raster should fit within; the longer side of the raster will be min(shape). Type
tuple

class
DLTile
(dltile_dict)[source]¶ A GeoContext that clips and projects
Scenes
to a single DLTile.DLTiles allow you to define a grid of arbitrary spacing, resolution, and overlap that can cover the globe. DLTiles are always in a UTM projection.
__init__
instantiates a DLTile from a dict returned byRaster.dltile
.It’s preferred to use the
DLTile.from_latlon
,DLTile.from_shape
, orDLTile.from_key
class methods to construct a DLTile GeoContext.
assign
(pad, raster_client=None)[source]¶ Return a copy of the DLTile with the pad value modified.
 Returns
new
 Return type

classmethod
from_key
(dltile_key, raster_client=None)[source]¶ Return a DLTile GeoContext from a DLTile key.

classmethod
from_latlon
(lat, lon, resolution, tilesize, pad, raster_client=None)[source]¶ Return a DLTile GeoContext that covers a latitude/longitude.
Where the point falls within the tile will vary, depending on the point and tiling parameters.
 Parameters
lat (float) – Latitude (WGS84)
lon (float) – Longitude (WGS84)
resolution (float) – Distance, in meters, that the edge of each pixel represents on the ground
tilesize (int) – Length of each side of the tile, in pixels
pad (int) – Number of extra pixels by which each side of the tile is buffered. This determines the number of pixels by which two tiles overlap.
raster_client (descarteslabs.client.services.Raster, optional, default None) – Unneeded in general use; lets you use a specific client instance with nondefault auth and parameters.
 Returns
tile
 Return type

classmethod
from_shape
(shape, resolution, tilesize, pad, raster_client=None)[source]¶ Return a list of DLTiles that intersect the given geometry.
 Parameters
shape (GeoJSONlike) – A GeoJSON dict, or object with a
__geo_interface__
. Must be inEPSG:4326
(WGS84 latlon) projection.resolution (float) – Distance, in meters, that the edge of each pixel represents on the ground.
tilesize (int) – Length of each side of the tile, in pixels.
pad (int) – Number of extra pixels by which each side of the tile is buffered. This determines the number of pixels by which two tiles overlap.
raster_client (
Raster
, optional, defaultNone
.) – Unneeded in general use; lets you use a specific client instance with nondefault auth and parameters.
 Returns
tiles
 Return type

property
bounds
¶ The
(min_x, min_y, max_x, max_y)
of the area covered by this DLTile, in the UTM coordinate reference system given inbounds_crs
. Type
tuple

property
bounds_crs
¶ The coordinate reference system of the
bounds
, given as an EPSG code (likeEPSG:32615
). A DLTile’s CRS is always UTM. Type
str

property
crs
¶ Coordinate reference system into which scenes will be projected. For DLTiles, this is always a UTM projection, given as an EPSG code.
 Type
str

property
geometry
¶ The polygon covered by this DLTile in WGS84 (latlon) coordinates
 Type
shapely.geometry.Polygon

property
geotrans
¶ The 6tuple GDAL geotrans for this DLTile in the shape
(a, b, c, d, e, f)
wherea is the top left pixel’s xcoordinateb is the westeast pixel resolutionc is the row rotation, always 0 for DLTilesd is the top left pixel’s ycoordinatee is the column rotation, always 0 for DLTilesf is the northsouth pixel resolution, always a negative value Type
tuple

property
key
¶ The DLTile’s key, which encodes the tiling parameters, and which number in the grid this tile is.
 Type
str

property
pad
¶ Number of extra pixels by which each side of the tile is buffered. This determines the number of pixels by which two tiles overlap.
 Type
int

property
proj4
¶ PROJ.4 definition for this DLTile’s coordinate reference system
 Type
str

property
raster_params
¶ The properties of this DLTile, as keyword arguments to use for
Raster.ndarray
orRaster.raster
. Type
dict

property
resolution
¶ Distance, in meters, that the edge of each pixel represents on the ground
 Type
float

property
ti
¶ The yindex of this tile in its grid
 Type
int

property
tilesize
¶ Length of each side of the tile, in pixels. Note that the total number of pixels along each side of an image is
tile_size + 2*padding
 Type
int

property
tj
¶ The xindex of this tile in its grid
 Type
int

property
wkt
¶ OGC WellKnown Text definition for this DLTile’s coordinate reference system
 Type
str

property
zone
¶ The UTM zone of this tile
 Type
int


class
GeoContext
[source]¶ Specifies spatial parameters to use when loading a raster from the Descartes Labs catalog.
Two Scenes loaded with the same GeoContext will result in images with the same shape (in pixels), covering the same spatial extent, regardless of the dimensions or projection of the original data.
Specifically, a fullydefined GeoContext specifies:
geometry to use as a cutline (WGS84), and/or bounds
resolution (m)
EPSG code of the output coordinate reference system
whether to align pixels to the output CRS (see docstring for
AOI.align_pixels
for more information)
GeoContexts are immutable.

property
raster_params
¶ The properties of this GeoContext, as keyword arguments to use for
Raster.ndarray
orRaster.raster
. Type
dict

class
XYZTile
(x, y, z)[source]¶ A GeoContext for XYZ tiles, such as those used in web maps.
The tiles are always 256x256 pixels, in the spherical Mercator or “Web Mercator” coordinate reference system (
EPSG:3857
).Requires the optional
mercantile
package. Parameters
x (int) – Xindex of the tile (increases going east)
y (int) – Yindex of the tile (increases going south)
z (int) – Zoom level of the tile

property
bounds
¶ The
(min_x, min_y, max_x, max_y)
of the area covered by this XYZTile, in spherical Mercator coordinates (EPSG:3857). Type
tuple

property
bounds_crs
¶ The coordinate reference system of the
bounds
. AlwaysEPSG:3857
(spherical Mercator, aka “Web Mercator”) Type
str

property
crs
¶ Coordinate reference system into which scenes will be projected. Always
EPSG:3857
(spherical Mercator, aka “Web Mercator”) Type
str

property
geometry
¶ The polygon covered by this XYZTile in
WGS84
(latlon) coordinates Type
shapely.geometry.Polygon

property
raster_params
¶ The properties of this XYZTile, as keyword arguments to use for
Raster.ndarray
orRaster.raster
. Type
dict

property
resolution
¶ Distance, in meters, that the edge of each pixel represents in the spherical Mercator (“Web Mercator”, EPSG:3857) projection.
 Type
float

property
tilesize
¶ Length of each side of the tile, in pixels. Always 256.
 Type
int

property
x
¶ Xindex of the tile (increases going east)
 Type
int

property
y
¶ Yindex of the tile (increases going south)
 Type
int

property
z
¶ Zoom level of the tile
 Type
int

EARTH_CIRCUMFERENCE_WGS84
= 40075016.68557849¶ Circumference of the earth, in meters, on the WGS84 ellipsoid