gemgis.visualization.read_raster

gemgis.visualization.read_raster(path=<class 'str'>, nodata_val: typing.Union[float, int] = None, name: str = 'Elevation [m]') → pyvista.core.pointset.PolyData

Reading a raster and returning a mesh

Parameters

• path (str) – Path to the raster, e.g. path='raster.tif'

• nodata_val (Union[float, int]) – Nodata value of the raster, e.g. nodata_val=9999.0

• name (str) – Name of the data array, e.g. name='Elevation [m], default is 'Elevation [m]'

Returns

mesh – PyVista mesh containing the raster values

Return type

pyvista.core.pointset.PolyData

New in version 1.0.x.

Changed in version 1.1.1: Set nodata value manually if no data value is provided and raster does not contain nodata values

Example

>>> # Loading Libraries and outputting mesh
>>> import gemgis as gg
>>> polydata = gg.visualization.read_raster(path='raster.tif', nodata_val=9999.0, name='Elevation [m]')
>>> polydata
Header
StructuredGrid  Information
N Cells         5595201
N Points        5600000
X Bounds        3.236e+07, 3.250e+07
Y Bounds        5.700e+06, 5.800e+06
Z Bounds        0.000e+00, 0.000e+00
Dimensions      2000, 2800, 1
N Arrays        1
Data Arrays
Name            Field   Type    N Comp  Min         Max
Elevation [m]   Points  float32 1       0.000e+00   5.038e+02

See also

convert_to_rgb

Converting bands to RGB values for plotting

drape_array_over_dem

Draping an array of the Digital Elevation Model