Quantization‎ > ‎

Vector Data


PYXIS DGGS SDK provides 2 methods to read, sample, and quantize data values from vector (feature) data files into DGGS cells:

  1. Rapid Rasterization: At any DGGS resolution a vector is represented as a rasterization of the vector with its associated attributes set as data values in the DGGS cells.
  2. Multi-Resolution GeoPacket: A vector tree is formed that represents the rasterization of the vectors at all resolutions of the DGGS while data values are represented as statistical summaries of the associated attribute values.

Sampling Vector Geometry

The first step in representing a DGGS GeoSource is to open the native vector file and extract the following from the file:.
  1. Spatial Domain Information - Information about the native coordinate system, native bbox, scale, data types, etc. 
  2. Data Definition - Describes the attribute value domain and value types return from the feature.
The second step is the rasterization of the associated vector by representing it as a collection of DGGS cells.  The subject of representing vector geometries in a regular square raster grid is well understood and widely implemented.  However, the rasterization of spherical geometries into hexagonal tessellations require some attention.  The PYXIS DGGS SDK provide optimized implementations of several techniques of vector rasterization:
  • Cell tracing of great circle arcs (O(cells x arcs) - border only)
  • Multi-resolution distance based vector rasterization ( O(log(cells) x arcs) ) - including support for buffers.
  • Multi-resolution, multi-details simplification based on tessellation error bounds ( O(log(cells) x log(arcs)) ) - including support for buffers.
  • Hybrid representation and indexing of vectors. (indexing O(log(cells) x log(arcs)), raster: O(1))
At code level, PYXVectorRegion provide the basic class for all vector based regions. and region can be rasterized on demand into DGGS set of cells using PYXVectorGeometryPYXVectorGeometry2 (optimized implementation).

Quantization of Vector Data

Once the cells of a DGGS area associated with a particular vector geometry the final step is to quantize vector attribute data into associated DGGS cells.  A bounding cell for each feature entity is calculated and the attribute is assigned to the cell.  The GDAL/OGR File Reader is used to identify all features that intersects a given cell. 

If a mutli-resolution GeoPacket is to be generated, statistical summaries are created for each numerical attribute value associated with the feature.  The rasterization of each feature and statistical summary are stored for each ISEA3H resolution to form a vector tree that enable progressive streaming of very large vectors and their attribute data.

Interactive Example

Rapid Rasterization logic as a naive implementation using PYXIS HTML5 canvas API. 

Vector rasterization

Vector Features on DGGS