Handling grids of data — gridData

Overview

This module contains classes that allow importing and exporting of simple gridded data, A grid is an N-dimensional array that represents a discrete mesh over a region of space. The array axes are taken to be parallel to the cartesian axes of this space. Together with this array we also store the edges, which are are (essentially) the cartesian coordinates of the intersections of the grid (mesh) lines on the axes. In this way the grid is anchored in space.

The Grid object can be resampled at arbitrary resolution (by interpolating the data). Standard algebraic operations are defined for grids on a point-wise basis (same as for numpy.ndarray).

Description

The package reads grid data from files, makes them available as a Grid object, and allows one to write out the data again.

A Grid consists of a rectangular, regular, N-dimensional array of data. It contains

  1. The position of the array cell edges.

  2. The array data itself.

This is equivalent to knowing

  1. The origin of the coordinate system (i.e. which data cell corresponds to (0,0,…,0)

  2. The spacing of the grid in each dimension.

  3. The data on a grid.

Grid objects have some convenient properties:

  • The data is represented as a numpy.ndarray in Grid.grid and thus can be directly manipulated with all the tools available in NumPy.

  • Grid instances can be manipulated arithmetically, e.g. one can simply add or subtract two of them and get another one, or multiply by a constant. Note that all operations are defined point-wise (see the numpy documentation for details) and that only grids defined on the same cell edges can be combined.

  • A Grid object can also be created from within python code e.g. from the output of the numpy.histogramdd() function.

  • The representation of the data is abstracted from the format that the files are saved in. This makes it straightforward to add additional readers for new formats.

  • The data can be written out again in formats that are understood by other programs such as VMD, ChimeraX or PyMOL.

Reading grid data files

Some Formats can be read directly from a file on disk:

g = Grid(filename)

filename could be, for instance, “density.dx”.

Constructing a Grid

Data from an n-dimensional array can be packaged as a Grid for convenient handling (especially export to other formats). The Grid class acts as a universal constructor:

g = Grid(ndarray, edges=edges)                 # from histogramdd
g = Grid(ndarray, origin=origin, delta=delta)  # from arbitrary data

g.export(filename, format)   # export to the desire format

See the doc string for Grid for details.

Formats

For the available file formats see Supported file formats.