Tools
vtu2geo
The conversion of the vtu format to arbitrary GIS formats is provided by
vtu2geo
located in tools/vtu2geo/main.py
. Requires vtk and gdal
+ python bindings for each.
This tool produces an internal shp file that corresponds to the triangulation and uses GDAL to rasterize this to an output geotiff.
The vtu files contain multiple variables. Therefore, each output geotiff
is a 1-band file corresponding to the selected output. The variables of
interest as set in the variables
list.
For vtu variables that are parameters (and therefore constant with
time), only 1 output file is needed. These are defined in the
parameters
list. Only 1 geotiff will be produced from these.
input_path
should be either a single .vtu file or a the .pvd file.
If a pvd file is given it will produce a tiff for each vtu.
output_path
is the output path. If output into the current folder is
wanted, use '.'
.
pixel_size
is the size of the raster cells in m^2.
var_resample_method
and param_resample_method
determine what
resampling method to use when calculating on a clipped raster.
# Configuration file for vtu2geo tool
# Input path to where output vtu files are located
input_path = 'meshes/SC1506124800.vtu'
output_path = 'meshes/'
# Output variables
variables = ['snowdepthavg'] #set to None to dump all variables
var_resample_method = {'snowdepthavg':'average'} # Methods to use when calculating clipped raster
# Output parameters
parameters = ['Elevation'] # paramters are one offs we want to extract from the vtu files
param_resample_method = {'Elevation':'average'} # Methods to use when calculating clipped raster
# Output pixel size that the mesh is interpolated to (?)
pixel_size = 30 # (m)