mais où vont mes rivières ?

[rammus-lol]

Bonjour, je viens de terminer ma couche rivière à partir de bd topage dont voici le csv.

topage_01_07_2025.csv

j'en ai fait l'importation pour la dernière fois avec cette commande précise openig@openig-georiviere:~/georiviere/var$ docker-compose run --rm web ./manage.py load_rivers /opt/georiviere-admin/var/topage_01_07_2025.gpkg --flush --name-attribute Nom_Riv

j'obtiens ce message :

WARNING: The UID variable is not set. Defaulting to a blank string.
WARNING: The GID variable is not set. Defaulting to a blank string.
WARNING: Found orphan containers (georiviere_nginx_1) for this project. If you removed or renamed this service in your compose file, you can run this command with the --remove-orphans flag to clean it up.
Starting georiviere-memcached-1 ... done
Starting georiviere-screamshotter-1 ... done
Starting georiviere-postgres-1      ... done
Starting georiviere-convertit-1     ... done
Waiting for postgres...
PostgreSQL started
Load rivers: 2625 features to import
Delete streams.....done!
Successfully import 2625 rivers and associated morphologies / status

cependant mon site georiviere-admin n'affiche rien malgré une importation visiblement considéré comme réussi je me suis donc attelé à chercher d'abord dans le conteneur postgres sans succès je vais regarder dans le conteneur web mais ça serait bien de dire où ça va dans la sortie bash

[babastienne]

Bonjour,

Bon je suppose vu que tu importes un gpkg que les géométries sont bien dans le fichier (par rapport au csv qui ne contient que les attributs).

Est-ce que la projection(SRID) des cours d'eau dans le gpkg est bien la même que le SRID configuré dans ton fichier de configuration custom.py ?

Si oui peux-tu vérifier que les cours d'eau sont bien existants en base de données ? Tu peux te connecter à la base en faisant docker compose run --rm web dbshell puis parcourir les tables et en particulier regarder le nombre d'éléments dans ta table river_stream ou quelque chose du style (tu peux lister les tables en faisant \dt il me semble).

[rammus-lol]

Oui mon geopackage est parfait aucune erreur de topologie
Pour le srid je viens de vérifier les deux sont identiques.
pour les tables dans le postgre il y a ces 4 là

river_classificationwaterpolicy                        
river_stream
river_stream_portals  
river_topology 

après inspection j'ai ça georivieredb=# select * from river_classificationwaterpolicy limit 10;

 id |        label         | structure_id
----+----------------------+--------------
  1 | Class├® cours dÔÇÖeau   |
  2 | Non cours dÔÇÖeau      |
  3 | Reclassement demand├® |
  4 | ├Ç d├®terminer         |
(4 rows)

georivieredb=# select * from river_stream limit 10;
 id |     name      | geom | ascent | descent |

georivieredb=# select * from river_stream_portals limit 10;
 id | stream_id | portal_id
----+-----------+-----------
(0 rows)

georivieredb=# select * from river_topology limit 10;
 id | start_position | end_position | qualified | stream_id
----+----------------+--------------+-----------+-----------
 11 |              0 |            1 | f         |         6
 12 |              0 |            1 | f         |         6
 13 |              0 |            1 | f         |         7
 14 |              0 |            1 | f         |         7
(4 rows)

river_stream a ce formalisme dans l'invite de commande
ensuite je me suis rendu dans le conteneur web dans le dossier /opt/georiviere-admin/georiviere/river/sql et j'ai ce contenu

-rw-r--r-- 1 django django 2538 Apr  7 09:25 post_10_triggers.sql
-rw-r--r-- 1 django django  162 Apr  7 09:25 pre_10_cleanup.sql

[babastienne]

Qu'est-ce que renvoi la commande select count(*) from river_stream; ?

[babastienne]

Peux-tu partager ton gpkg ? Est-il topologiquement correct ? (géométries simples, valides, les intersections coupent les cours d'eau et pas d'overlap ?)

[rammus-lol]

georivieredb=# select count(*) from river_stream;
count

 2

(1 row)
(j'ai deux rivière de test digitalisées)

topage_01_07_2025.zip
je n'ai pas trouvé comment on partage de geopackage directement sur github donc je vous ai zippé le fichier

[rammus-lol]

je précise que je viens de retenter avec un geopackage avec des noms en minuscule et ça n'a rien changé.

[babastienne]

Hum. A première vue ça m'a l'air plutôt correct. Bon y a pas de segmentation à chaque intersection comme recommandé mais en soit GeoRivière est supposé le gérer il me semble.

Il faudrait que j'installe un GeoRivière et que j'essaye d'importer ton fichier pour pousser l'analyse et trouver le problème mais là j'ai pas vraiment le temps.

En attendant tu peux essayer une approche itérative :

• Je suppose que tu as déjà suivi la doc et vérifié que tu pouvais créer manuellement des cours d'eau dans l'interface ?
• plutôt que d'importer 2000+ cours d'eau essaye déjà d'importer un fichier avec un seul cours d'eau et regarde si ça fonctionne ?
• Sinon est-ce que les migrations en bdd ont toutes été correctement lancées sans erreur (lors de l'étape update.sh) ?

[rammus-lol]

Je n'ai aucun problème pour rentrer des cours d'eau via l'interface, cependant rentrer un seul cours d'eau n'a eu aucune différence avec les fois où j'en ai rentré 2000.
Par contre pour le troisième point je ne vois même pas de quoi vous parlez
et qu'est ce que vous appelez segmentation ?

[babastienne]

Ok j'ai trouvé le problème. Toutes les features de ta couche sont des multi-linestrings. Pour être importées il faut que ce soit des linestrings uniquement.

Pour transformer tes données et mieux comprendre comment obtenir un réseau topologiquement correct tu peux regarder :

• Un guide sur la segmentation dynamique
• La doc de Geotrek pour nettoyer les données (section requirements et clean paths).

Normalement en appliquant cette procédure tu devrai avoir un réseau importable dans GeoRiviere.

[rammus-lol]

Merci beuacoup cependant je vois que dans geotrek il est possible d'importer des données dans la base via l'utilitaire de connexion postgis de qgis c'est possible avec georivière ?

[babastienne]

Je ne crois pas que ce soit possible car Geotrek utilise des triggers pour calculer ensuite des données mais GeoRivière n'effectue pas ces calculs en sql, directement en python (et donc ça implique d'utiliser la commande d'import).

De mon côté après traitement de la couche j'ai réussi à importer les données :

[babastienne]

Procédure suivie :

• Ouverture de ta couche gpkg intiale dans QGis
• Clean geometries
  • search “v_clean” in “Processing toolbox”
  • select following options in cleaning tool: break, snap, duplicate (ou rmdup), rmline, rmdangle, chdangle, bpol, prune
  • in threshold enter 2,2,2,2,2,2,2,2 (2 meters for each option)
• delete duplicate geometries
  • search “duplicate” in “Processing toolbox”
• Export de la couche obtenue et import dans GeoRivière ok.

---

La solution de facilité ...

topage-clean.gpkg.tar.gz

[rammus-lol]

alors ça a marché merci beaucoup mais c'est largement trop segmenté, je me retrouve avec des cours d'eau de 4 km coupé en 3 tronçons sans fondement réel finalement

[babastienne]

Ok. En effet j'ai fait pris quelques minutes pour faire ce fichier modifié mais je n'ai pas la connaissance territoriale ni l'expertise adaptée. C'est à vous d'ajuster ces paramètres en fonction de vos besoins, de fusionner des segments selon les spécificités de votre territoire, etc.

L'idée principale est qu'il faut une couche uniquement avec des linestrings, on a prouvé que ça fonctionnait comme ça. Maintenant le reste est parfaitement ajustable, à vous de jouer.

[babastienne]

D'ailleurs ce qu'illustre ce ticket c'est qu'il manque une partie de documentation sur la préparation des données pour l'import. N'hésitez pas à proposer une contribution pour enrichir la documentation afin que les autres utilisateurs à l'avenir ne se heurtent pas au même problème.