limitation de la construction du maillage avec l'argument bbox ?

[CatherineFouchier]

Salut
quand on construit un maillage avec l'argument bbox=[xmin, xmax, ymin, ymax] et les attributs d'exutoires, on n'a plus les attributs mesh["flwdst"] et mesh["gauge_pos"]. J'ai ce problème sur mes données et aussi sur les données du tuto Cance (code utilisé avec un simple copié-collé).
Merci pour les infos sur cette apparente limitation de la construction du maillage avec l'argument bbox.
Catherine

[inoelloc]

Salut @CatherineFouchier, oui avec une bbox, les coordonnées des exutoires ne sont pas prises en compte dans la génération du maillage. L'idée à la base était que si on choisit une bbox, on veut faire une simulation sur tout le domaine et que quand on choisit des exutoires, on fait une simulation uniquement sur des points d’intérêt. Il y a peut etre quelque chose à faire pour permettre un maillage mix, mais je ne crois pas que ce soit trop d'actualité dans les développements à venir !

[CatherineFouchier]

Salut François
Merci pour tes retours.
Le contexte pour cette demande ce sont des petits bv soumis à aléas gravitaires dans les Alpes Maritimes et affluents de différents cours d’eau indépendants : Roya, Var, Durance, …
Il m’avait semblé plus pratique de construire un mesh avec une bbox sur la zone d’intérêt (qui n’a pas besoin d’aller jusqu’aux cours d’eau majeurs) et de mettre les exutoires de ces petits bv en arguments du mesh pour y obtenir les Qsim. L'exemple dans le tuto donne à penser que l'utilisation de la bbox est compatible avec le passage des coordonnées des exutoires :
mesh = smash.factory.generate_mesh(
flwdir_path="./Cance-dataset/France_flwdir.tif",
bbox=[800000., 850000., 6440000., 6490000.],
x=list(gauge_attributes["x"]),
y=list(gauge_attributes["y"]),
area=list(gauge_attributes["area"] * 1e6), # Convert km² to m²
code=list(gauge_attributes["code"]),
)
Pour sortir les qsim aux exutoires, Pierre m'a suggéré de faire un post-traitement sur les grilles issues de model.forward_run(return_options={"q_domain": True}) mais c'est peut-être consommer pas mal de mémoire et de puissance de calcul si on veut juste récupérer les débits aux exutoires.
Pas de souci si cette question n'est pas dans les priorités :) !
A+

[maximejay]

Bonjour Catherine,

Si cette fonctionnalité existe dans la version 1.2 de SMASH. Je l'ai ajouté cet automne via cette pull request:

#430

Kevin l'utilise depuis avec SmashBox. Quelle version de Smash utilises tu ?

[CatherineFouchier]

Salut Maxime
tu peux me passer la syntaxe python stp pour demander la version de smash ? Je l'ai oubliée. Merci !
Catherine

[nghi-truyen]

Salut Catherine,
c'est pip install hydro-smash==1.2.0 pour installer la version 1.2.

[maximejay]

Salut,
Coté Python ca fait un truc comme ca:

bbox_tmp = {"left": 925000.0, "bottom": 6230000.0, "right": 974000.0, "top": 6255000.0}
bbox = [
bbox_tmp["left"],
bbox_tmp["right"],
bbox_tmp["bottom"],
bbox_tmp["top"],
]

smash.factory.generate_mesh(
flwdir_path="",
bbox=bbox,
x=np.array(stations_calage[columns["coord_x"]][:]),
y=np.array(stations_calage[columns["coord_y"]][:]),
area=np.array(
stations_calage[columns["area"]][:] * 1e6
), # Convert km² to m²
code=np.array(stations_calage[columns["id"]][:]),
epsg="",
shp_path="",
max_depth=1,
area_error_th=0.2,
)
Nouveau aussi : area_error_th: Tolerance error during the positionning of the outlets. If the Error (Ssim-Sobs)/Sobs > area_error_th, the outlet will be excluded.

[CatherineFouchier]

Salut Catherine, c'est pip install hydro-smash==1.2.0 pour installer la version 1.2.

Merci Truyen, as-tu également stp la syntaxe pour seulement vérifier ma version actuelle de Smash, sans installer une autre version ?

[nghi-truyen]

Salut Catherine, c'est pip install hydro-smash==1.2.0 pour installer la version 1.2.

Merci Truyen, as-tu également stp la syntaxe pour seulement vérifier ma version actuelle de Smash, sans installer une autre version ?

Pour vérifier ta version smash :
python -c "import smash; print(smash.__version__)"

Après tu peux le supprimer et reinstaller s'il est pas à jour (la 1.2.1 vient d'être sorti donc tu peux installer cette version) :
pip uninstall hydro-smash
pip install hydro-smash==1.2.1

[ahmad-3100]

Hello everyone,
I have tested the mix mesh generation with bbox and outlets over cance, the problem is that when I add outlets the mesh is generated with respect to the outlets domain not bbox domain. See attached figures.

[maximejay]

Salut Ahmad,
Pour moi ton mesh semble cohérent. En effet le mesh est généré seulement pour les cellules actives. Il y a t-il moyen de désactiver les active_cell only ?

Je vais regarder cela...

[maximejay]

Salut Ahmad,

Dans ce dépôt https://github.com/maximejay/smash/tree/mesh_bbox_with_no_active_cell, il y a une version de smash qui génère des mesh avec tout les pixels. Il faut utiliser la branche mesh_bbox_with_no_active_cell.

J'ai tester la génération du maillage. Peux tu tester le code complet avec une simulation ?