[Donnees] Données (ouvertes) des cartes géologiques du BRGM

Message posté par : jfmoyen ---------------------------------------- Alors, en tant que géologue (principalement) et pas du tout gémoaticien, à part comme utilisateur, je peux commenter un peu sur la construction de la BD-Charm du point de vue géologique. Si mes étudiants de licence passent par là ils peuvent lire aussi, c'était le cours de la semaine dernière et ce sera à l'exam :-) En plusieurs parties parce que c'est .. long. I. Des cartes géologiques à la BD-Charm Les cartes géologiques (en général, mais en France en particulier pour ce qui nous intéresse) ont été levées feuille par feuille. Il y a 1000 et quelques coupures au 1/50 000 en France métropolitaine, et elles ont été levées depuis les années 60 jusqu'aux années 2020, les dernières feuilles datent de un an ou deux. Le travail a donc été fait à des époques différentes, par des personnes différentes. Malgré les efforts de standardisation, ça implique des différences: - une carte géologique est une interprétation du sous-sol (si vous regardez n'importe quel paysage, vous verrez surtout des champs et des arbres, pas tellement de roches, il faut donc décider ce qu'il y a sous le champ !). Les interprétations ne sont pas faites au hasard, elles sont ancrées sur la façon dont on comprend la géologie au moment où on faut la feuille, on va donc dessiner les choses de façons différentes à différents moments (Le point de rencontre entre 745 St Etienne, années 60; 768 Monistrol et 769, Annonay, années 80; et 744 Firminy, années 2010, est assez mythique...). - c'est en particulier vrai dans les terrains métamorphiques, où depuis les années 2000 on a fait des progrès énormes en géochoronologie (datation) et quantification des conditions de température/pression, et ça a très largement renouvellé notre façon de comprendre, et partant de dessiner ces choses. Voyez par exemple la lmite entre 742 (St Germain Lembron, 2005) et ses voisines (années 80). - des géologues différents ont des expertises différentes, et sont plus ou moins bons dans différents types de terrain. En principe le BRGM essaye de mettre des gens pertinents, mais des fois il y a une partie d'une feuille qui se retrouve confiée à quelqu'un dont ce n'est pas la principale compétence, et c'est souvent moins bon. Comparez la représentation, d'une part des granites, d'autre part des volcans entre Mende (862) et St Geniez d'Olt (861). - il y a aussi un choix fondamental sur la façon de représenter les "formations superficielles" (éboulis, alluvions de rivières....). Est-ce qu'on les ignore pour se consacrer à la "vraie" géologie, ou est-ce qu'on les dessine parce que c'est capital par exemple pour de la géotechnie ? Le BRGM n'a pas eu une politique très cohérente sur ce point et les choix ont flucutés de carte en carte, parfois à l'intérieur même d'une carte. (Voyez la limite entre Mende 862 et St Chély 838; sur 918 La Javie comparez le secteur de Barles et la vallée de la Bléone près de Prads). D'où l'idée de l'harmonisation. Elle a été faite par département, le processus est pas mal décrit dans certains rapports publics du BRGM, celui-ci si je me souviens bien (http://infoterre.brgm.fr/rapports/RP-57097-FR.pdf) donne des explications assez détaillées par exemple. L'harmonisation ce fait à deux niveaux: 1) Harmonisation des légendes - on épluche les légendes de toutes les feuilles couvrant le département, et on fait une liste des "cartons" (les entrées de légende); - quand deux cartons se retrouvent sur différentes cartes, éventuellement sous des noms différents, on les regroupe; - dans certains cas on est obligé de regrouper deux formations pour harmoniser, par exemple si une feuille a du j1 et et du j2-3 et sa voisine du j1-2 et du j3, on peut harmoniser en j1-3, au prix d'une perte de détails - le géologique qui fait l'harmonisation décide alors d'une légende unifiée pour son département, à la fois en ce qui concerne les notations et les descriptions (issues des légendes des différentes feuilles) 2) Harmonisation des contours - quand on a des contours qui ne se rejoignent pas, il faut les redessiner près de la jonction (l'exemple de Firminy plus haut, regardez autour du point quadruple comment le cartographe a "bricolé" de différentes façons, en faisant s'amortir des formations); - idéalement, il faudrait apporter le même niveau d'information sur les deux cartes (formations superficielles), mais ça c'est pas encore fait... Comme l'harmonisation est faite par département, on a maintenant les mêmes problèmes aux frontières entre département ! Voyez la limite entre 42 et 43 dans mon exemple... Le prochain travail sera donc une harmonisation à l'échelle nationale (et même Européenne....). Le BRGM s'est lancé depuis une quinzaine d'années dans un projet ambitieux, le RGF (Référentiel Géologique de France, https://rgf.brgm.fr ), qui vise à faire ce travail. Il se fait d'abord par entité géologique (les Vosges, les Pyrénées, le Bassin Parisien...) mais il a vocation à couvrir tout le territoire. II. Les fichiers de la BD-Charm Malheureusement, autant je suis admiratif du travail de géologie, autant je trouve que le travail de géomatique a été un peu mal fait ... en fait, il est sans doute fait pour très bien fonctionner en interne, sur un poste BRGM avec accès avec toutes les configurations de la maison, mais la diffusion donne l'impression d'avoir été faite "parce qu'on nous le demande" et pas "parce qu'on veut que les données soient vraiment utilisables". Voyez mes messages de Mars 2022 pour les détails. Concrètement, ça veut dire qu'un utilisateur lambda dispose des utilisations suivantes: 1) Pour ce qui est de la géologie: - la notation (dans le champ éponyme), sous forme d'une combinaison de lettres et de chiffres plus ou moins conventionelles, harmonisée à l'échelle du département. La convention générale est que les formations superficielles ont une lettre majuscule qui indique leur nature, et une minuscule leur âge (coucou les L1 :-) ), par exemple Gz, Fx etc. Les roches sédimentaires une lettre minuscule et un chiffre (j1, c4, d2...). Les roches magmatiques et métamorphiques une combinaison de lettres grecques et romaines, et là c'est le drame : dans le fichier du BRGM ces lettres ne sont codées dans AUCUN encoding identifiable. Par exemple, un "béta" (en unicode, U+03B2) est codé avec U+00E2, qui devrait correspondre à un a circonflexe). Je n'ai trouvé aucun encoding dans lequel U+00E2 soit un beta... Pour afficher correctement ces notations il faut utiliser la police "BRGM_NOT" dans laquelle le code U+00E2 est représenté par un glyphe qui a la forme d'un béta (mais est, en fait, un â). On ne peut donc pas représenter les étiquettes avec une autre police (ni les lire correctement dans un autre logiciel d'ailleurs). Une partie de mon boulot de l'an passé avait été de convertir ces caractères foireux en leur équivalent unicode, ce qui permet de représenter la légende dans la police qu'on veut. - une DESCRiption, assez brève (par exemple "Eboulis récents à actuels, localement anciens. Parfois éboulis remaniés par ruisselement, éboulis à gros blocs, éboulis glissé, dépôts soliflués"), reprise des cartons des cartes d'origine. Il n'y a ni decription longue (pour celà il faut lire les notices, où on peut facilement avoir une page par carton dans les cartes récentes), ni organisation de la donnée, en tout cas dans la version 2019 de la BD-Charm (la version 2015 avait des champs AGE, ERE, NATURE qui permettait par exemple de savoir que m2a-b, "Burdigalien inférieur et moyen. Molasses indifférenciées et marnes", est un âge Burdigalien, une ère Cénozoïque, une nature sédimentaire.... cette info a disparu dans la version 2019). Et c'est tout. Vous n'avez rien de plus, débrouillez vous avec ça. 2) Pour ce qui est de la symbologie (qui n'est pas totalement normalisée, même si le BRGM essaye d'avoir des couleurs à peu près conventionelles) : elle se compose d'une couleur de fond, et éventuellement d'une "surcharge" avec des symboles superposés. Dans la BD-Charm elle est décidée par le cartographe, qui essaye normalement de maintenair une cohérence avec ses voisins et avec les cartes d'origines (mais en fait c'est un travail impossible, il y a je ne sais plus combien de milliers d'entrées différentes dans la légende..). - C_FOND, M_FOND etc. sont les codes CMJK des couleurs de fond, qu'on peut facilement visualiser dans QGIS avec color_cmyk("C_FOND","M_FOND","J_FOND","B_FOND"). - C_SYMB, M_SYMB etc. sont les couleurs des symboles de la surcharge. - NOM_SYMB est un code (comme S02_007) qui fait référence à un symbole svg mis à disposition par le BRGM dans son fichier "outils" (https://infoterre.brgm.fr/actualites/nouvelle-version-cartes-geologiques-ha… , le fichier est http://infoterre.brgm.fr/telechargements/BDCharm50/Outils.zip). Dans ce zip on a un dossier "Polygon_Patterns" qui contient une centaine de fichiers svg avec le bon nom, par exemple S02_007.svg Vous pouvez donc reconstruire à la main la symbologie en utilisant (dans QGIS, désolé pour les utilisateurs d'ARcGIS) en utilisant largement les valeurs définies par les epxressions. Il suffit de donner à chaque polygone: - une couleur de fond définie par l'expression color_cmyk("C_FOND","M_FOND","J_FOND","B_FOND") - une couche de symboles svg, où le fichier svg a un chemin défini par une expression bien sentie (un truc du genre 'chemin/vers/le/dossier/outils' || "NOM_SYMB" ) et une couleur définie par color_cmyk("C_SYMB","M_SYMB","J_SYMB","B_SYMB") - une étiquette basée sur NOTATION et formattée en utilisant la police BRGM_NOT, à l'exclusion de toute autre. L'ensemble de cette symbologie est aussi fournie, je crois dans l'archive de chaque département, sous forme d'un fichier qml. Je ne l'utilise plus, donc j'ai un peu oublié comment il est construit, mais - il est sous la forme d'un "ensemble de règles" ("si NOTATION = 'm2a-b' ALORS utilisez telle couleur et tels symboles), donc PAS une symbologie catégorisée qui permettrait de voir apparaître la légende dans la liste des couches; - il doit faire référence aux symboles svg et s'attendre à ce qu'ils soient dans un dossier particulier, j'ai oublié lequel. III. Mes fichiers Pour ces raisons, j'ai donc fait ma propre compilation de la BD-Charm, sous forme d'un geopackage qui comporte la BD complète (tout les départements), requalifiés de mon mieux (les attributs manquants), et avec les notations reconverties. Et aussi avec une symbologie pas tout à fait complète. Concrètement mon geopackage inclut - Tout les polygones - Le plus d'attributs possibles, c'est à dire ceux de la version 2019 (cf plus haut) mais aussi, quand ils existent, ceux de 2015 un peu plus détaillés (il y a en particulier des ERE, AGE, etc) - Les NOTATIONS reconverties en utf-8 - Les svg des symboles tous regroupés dans une table (on peut donc assez facilement les retrouver avec quelque chose comme 'base64:'|| to_base64(attribute(get_feature( 'svg','name',"NOM_SYMB"),'image')) ) - Je crois aussi un projet QGIS prêt à l'usage (j'en suis certain pour la carte au 1/1 000 000, moins pour les 1/50 000). J'ai mis ce fichier à disposition, je crois qu'il a été archivé chez C. Quest (http://data.cquest.org/ ) mais je ne le retrouve plus, le site a été réorganisé. Si le fichier est quelque part, le signaler sur ce fil. Si vous en avez besoin, me contacter par PM (je ne passe pas très souvent sur le forum mais normalement il y a des notifications), ou par mail (avec mon nom et quelques mots clés bien choisis je ne suis pas très dur à trouver). Dernière remarque, enfin : la carte au 1/50 000 est très (trop ?) détaillée et probablement plus que de détails que nécessaires pour la plupart des applications. En pratique, la carte au 1/ 1 000 000 (qui est d'ailleurs mieux qualifiée !) a des contours qui restent utilisables jusqu'à des échelles entre le 250 et le 500k. Il est donc possible, selon vos applications (un département entier) que vous soyez mieux servi par la 1M. Il existe aussi quelques feuiilles papier au 250k, mais je ne crois pas qu'elles soient disponibles en format vecteur, et quoi qu'il en soit la couverture n'est pas complète... ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=363324#p363324 Pour y répondre : donnees(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net