cnr-engineering / tatooinemesher Goto Github PK

View Code? Open in Web Editor NEW

11.0 6.0 6.0 22.33 MB

Interpolate and/or mesh rivers from cross-sections profiles (and constraint lines)

License: GNU General Public License v3.0

Python 90.36% Makefile 0.25% TeX 9.24% Shell 0.15%

mesher commandline-tool python3 pre-processing

tatooinemesher's Introduction

TatooineMesher

Tested versions: 3.7, 3.8, 3.9, 3.10 and 3.11.

Channel mesher and interpolator from 1D cross-sections and constraint lines

Description

Set of 4 command line scripts:

densify_cross_sections.py : interpolate initial and intermediate cross-sections
mesh_and_interpolate.py: mesh and/or interpolate from a set of cross-sections
mesh_mascaret_run.py : visualize Mascaret model geometry and results
mesh_crue10_run.py : visualize Crue10 model geometry and results

The command line scripts of TatooineMesher are located in the cli folder.

See wiki pages to learn how to use them.

Pre-requisites and installation

Python 3 with packages listed in requirements.txt. To install them use:

pip install -r requirements.txt

PyTelTools and Crue10_tools are dependencies that can be installed in another directory and linked with the PYTHONPATH environment variable. These packages are used to read/write some specific file formats of Telemac, Mascaret and Crue10 tools.

Reference

TatooineMesher: Anisotropic interpolation from 1D cross-sections and 2D channel mesher, L. Duron, F.-X. Cierco and K. Saad, XXVIth TELEMAC-MASCARET User Conference 2019, Toulouse (France). See paper and presentation.

tatooinemesher's People

Contributors

Stargazers

Watchers

Forkers

arnaud-morvan peter0016 ch-m-emery dreyerlivio acroucher nicogodet

tatooinemesher's Issues

Support superimposed constraint lines

Superimposed constraint lines are not supported yet.
Therefore meshing a 1D geometry files do not use constraint lines (minor/major/others limits).

Crossing of constraint lines should be tolerated only if the intersection corresponds to a point along a single cross-section or if constraint lines are exactly superimposed between consecutive constraint lines.

Make GDAL dependency optionnal and raise error in `mesh_crue10_run.py` if not install

GDAL can be very painfull to install (especially on windows without conda...).

As mesh_crue10_run.py is the only script that depends on GDAL when a specific argument is provided, I suggest to remove GDAL from requierements.txt and raise an error if args.infile_dem is provided and GDAL is not available.

If this is OK, I can create a PR.

Use hydraulic axis to mesh/interpolate in case of no constraint lines

Instead of having straight left and right banks, use the hydraulic axis.
The ratio between left and right part of the cross-section is assumed to be linear between consecutive cross-sections.

Allow use of Z values in axis to constraint lowest point in interpolation

Some times, lowest point in bed river is surveyed (especially when distance between sections is high). This axis can catch pit and natural weir.

Fix issues with numpy >=1.24

In numpy 1.24, np.float and np.int are deprecated and should be replace with builtin float and int.

np.float was a deprecated alias for the builtin float. To avoid this error in existing code, use float by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use np.float64 here.

https://numpy.org/devdocs/release/1.20.0-notes.html#using-the-aliases-of-builtin-types-like-np-int-is-deprecated

Remove H from a run or restart one without H

If H is present in a run, delete it if possible (I don't know if it's possible) or restart one without H.

Missing Z values above initial cross sections

Data to reproduce : tatooine.zip

Command :

cli\mesh_and_interpolate.py \
    --infile_constraint_lines constraint_lines.shp \
    --long_step 1 \
    --lat_step 1 \
    --attr_cross_sections sec_id \
    axe_cambo.shp \
    tatooine_input.shp \
    --outfile_nodes interpolate_points.shp \
    --interp_constraint_lines FINITE_DIFF

What's the issue ?

In output file, interpolated point above initial cross sections are poplated with nan values.

Which is annoying when rasterized

tatooine_input.shp was prepared by PreCourlis