gstools.field.Field
- class gstools.field.Field(model=None, value_type='scalar', mean=None, normalizer=None, trend=None, dim=None)[source]
Bases:
object
A base class for random fields, kriging fields, etc.
- Parameters
model (
CovModel
, optional) – Covariance Model related to the field.value_type (
str
, optional) – Value type of the field. Either “scalar” or “vector”. The default is “scalar”.mean (
None
orfloat
orcallable
, optional) – Mean of the field if wanted. Could also be a callable. The default is None.normalizer (
None
orNormalizer
, optional) – Normalizer to be applied to the field. The default is None.trend (
None
orfloat
orcallable
, optional) – Trend of the denormalized fields. If no normalizer is applied, this behaves equal to ‘mean’. The default is None.dim (
None
orint
, optional) – Dimension of the field if no model is given.
- Attributes
all_fields
list
: All fields as stacked list.dim
int
: Dimension of the field.field_names
list
: Names of present fields.field_shape
tuple
: The shape of the field.latlon
bool
: Whether the field depends on geographical coords.mean
mesh_type
str
: The mesh type of the field.model
CovModel
: The covariance model of the field.name
str
: The name of the class.normalizer
Normalizer
: Normalizer of the field.pos
tuple
: The position tuple of the field.temporal
bool
: Whether the field depends on time.trend
value_type
str
: Type of the field values (scalar, vector).
Methods
__call__
([pos, field, mesh_type, ...])Generate the field.
delete_fields
([select])Delete selected fields.
get_store_config
(store[, default, fld_cnt])Get storage configuration from given selection.
mesh
(mesh[, points, direction, name])Generate a field on a given meshio, ogs5py or PyVista mesh.
plot
([field, fig, ax])Plot the spatial random field.
post_field
(field[, name, process, save])Postprocessing field values.
pre_pos
([pos, mesh_type, info])Preprocessing positions and mesh_type.
set_pos
(pos[, mesh_type, info])Set positions and mesh_type.
structured
(*args, **kwargs)Generate a field on a structured mesh.
to_pyvista
([field_select, fieldname])Create a VTK/PyVista grid of the stored field.
transform
(method[, field, store, process])Apply field transformation.
unstructured
(*args, **kwargs)Generate a field on an unstructured mesh.
vtk_export
(filename[, field_select, fieldname])Export the stored field to vtk.
- __call__(pos=None, field=None, mesh_type='unstructured', post_process=True, store=True)[source]
Generate the field.
- Parameters
pos (
list
, optional) – the position tuple, containing main direction and transversal directionsfield (
numpy.ndarray
orNone
, optional) – the field values. Will be all zeros ifNone
is given.mesh_type (
str
, optional) – ‘structured’ / ‘unstructured’. Default: ‘unstructured’post_process (
bool
, optional) – Whether to apply mean, normalizer and trend to the field. Default: Truestore (
str
orbool
, optional) – Whether to store field (True/False) with default name or with specified name. The default isTrue
for default name “field”.
- Returns
field – the field values.
- Return type
- get_store_config(store, default=None, fld_cnt=None)[source]
Get storage configuration from given selection.
- Parameters
store (
str
orbool
orlist
, optional) – Whether to store fields (True/False) with default names or with specified names. The default isTrue
for default names.default (
str
orlist
, optional) – Default field names. The default is “field”.fld_cnt (
None
orint
, optional) – Number of fields when using lists. The default is None.
- Returns
- mesh(mesh, points='centroids', direction='all', name='field', **kwargs)[source]
Generate a field on a given meshio, ogs5py or PyVista mesh.
- Parameters
mesh (meshio.Mesh or ogs5py.MSH or PyVista mesh) – The given mesh
points (
str
, optional) – The points to evaluate the field at. Either the “centroids” of the mesh cells (calculated as mean of the cell vertices) or the “points” of the given mesh. Default: “centroids”direction (
str
orlist
, optional) – Here you can state which direction should be chosen for lower dimension. For example, if you got a 2D mesh in xz direction, you have to pass “xz”. By default, all directions are used. One can also pass a list of indices. Default: “all”name (
str
orlist
ofstr
, optional) – Name(s) to store the field(s) in the given mesh as point_data or cell_data. If to few names are given, digits will be appended. Default: “field”**kwargs – Keyword arguments forwarded to
__call__
.
Notes
This will store the field in the given mesh under the given name, if a meshio or PyVista mesh was given.
- See:
- plot(field='field', fig=None, ax=None, **kwargs)[source]
Plot the spatial random field.
- Parameters
field (
str
, optional) – Field that should be plotted. Default: “field”fig (
Figure
orNone
) – Figure to plot the axes on. If None, a new one will be created. Default: Noneax (
Axes
orNone
) – Axes to plot on. If None, a new one will be added to the figure. Default: None**kwargs – Forwarded to the plotting routine.
- post_field(field, name='field', process=True, save=True)[source]
Postprocessing field values.
- Parameters
field (
numpy.ndarray
) – Field values.name (
str
, optional) – Name. to store the field. The default is “field”.process (
bool
, optional) – Whether to process field to apply mean, normalizer and trend. The default is True.save (
bool
, optional) – Whether to store the field under the given name. The default is True.
- Returns
field – Processed field values.
- Return type
- pre_pos(pos=None, mesh_type='unstructured', info=False)[source]
Preprocessing positions and mesh_type.
- Parameters
- Returns
iso_pos ((d, n),
numpy.ndarray
) – Isometrized position tuple.shape (
tuple
) – Shape of the resulting field.info (
dict
, optional) – Information about settings.
Warning
When setting a new position tuple that differs from the present one, all stored fields will be deleted.
- set_pos(pos, mesh_type='unstructured', info=False)[source]
Set positions and mesh_type.
- Parameters
- Returns
info – Information about settings.
- Return type
dict
, optional
Warning
When setting a new position tuple that differs from the present one, all stored fields will be deleted.
- to_pyvista(field_select='field', fieldname='field')[source]
Create a VTK/PyVista grid of the stored field.
- transform(method, field='field', store=True, process=False, **kwargs)[source]
Apply field transformation.
- Parameters
method (
str
) – Method to use. See gstools.transform for available transformations.field (
str
, optional) – Name of field to be transformed. The default is “field”.store (
str
orbool
, optional) – Whether to store field inplace (True/False) or under a given name. The default is True.process (
bool
, optional) – Whether to process in/out fields with trend, normalizer and mean of given Field instance. The default is False.**kwargs – Keyword arguments forwarded to selected method.
- Raises
ValueError – When method is unknown.
- Returns
Transformed field.
- Return type
- vtk_export(filename, field_select='field', fieldname='field')[source]
Export the stored field to vtk.
- Parameters
filename (
str
) – Filename of the file to be saved, including the path. Note that an ending (.vtr or .vtu) will be added to the name.field_select (
str
, optional) – Field that should be stored. Can be: “field”, “raw_field”, “krige_field”, “err_field” or “krige_var”. Default: “field”fieldname (
str
, optional) – Name of the field in the VTK file. Default: “field”
- property normalizer
Normalizer of the field.
- Type