{"id":12338,"date":"2024-01-30T09:00:00","date_gmt":"2024-01-30T08:00:00","guid":{"rendered":"https:\/\/www.sigterritoires.fr\/?p=12338"},"modified":"2025-05-08T09:16:44","modified_gmt":"2025-05-08T07:16:44","slug":"creation-dune-base-de-donnees-geopackage-pour-les-cartes-enc-actualisation-du-projet-en-cours","status":"publish","type":"post","link":"https:\/\/www.sigterritoires.fr\/index.php\/creation-dune-base-de-donnees-geopackage-pour-les-cartes-enc-actualisation-du-projet-en-cours\/","title":{"rendered":"Cr\u00e9ation d\u2019une base de donn\u00e9es Geopackage pour les cartes ENC (Actualisation du projet en cours)"},"content":{"rendered":"\n

Cet article sous sa forme compl\u00e8te (codes et liens) sera un des r\u00e9sultats du travail r\u00e9sultant du financement collaboratif que nous avons mis en place: https:\/\/www.sigterritoires.fr\/index.php\/financement-collaboratif-pour-lintegration-de-donnees-marines-dans-qgis\/<\/a>. D\u00e8s que le montant vis\u00e9 sera atteint, l’article sera mis \u00e0 jours dans sa version compl\u00e8te.<\/em><\/strong><\/p>\n\n\n\n

Dans le cas o\u00f9 nous n’atteindrions pas le montant vis\u00e9, nous enverrons la version compl\u00e8te \u00e0 tous les contributeurs.<\/em><\/strong><\/p>\n\n\n\n

Le contenu de l’article indique, au jour d’aujourd’hui, le contenu du projet pour les fichiers Geopackage<\/a>, \u00e0 l’issue des multiples \u00e9changes avec nos contributeurs. <\/em><\/strong><\/p>\n\n\n\n

…………………………………………………………………………………………………..<\/p>\n\n\n\n

Le format S57<\/a> est complexe et ne respecte pas totalement les standards des formats SIG. Sa conversion avec GDAL est donc soit complexe, soit incompl\u00e8te. Dans cet article nous allons rentrer dans le d\u00e9tail de la transformation compl\u00e8te du S57 en tables de fichier Geopackage.<\/p>\n\n\n\n

Les types de g\u00e9om\u00e9trie<\/h2>\n\n\n\n
Les \u00ab couches \u00bb S57 sont des classes d\u2019objets. Par exemple, les diff\u00e9rentes zones terrestres sont cod\u00e9es dans la classe d\u2019objets LNDARE.<\/p>\n\n\n\n
La d\u00e9finition de cette classe d\u2019objets est:<\/p>\n\n\n\n
Geo object: Land area (LNDARE) (P,L,A)
Attributes: CONDTN OBJNAM NOBJNM STATUS INFORM NINFOM<\/strong><\/p>\n\n\n\n
Si tous les objets LNDARE poss\u00e8dent les m\u00eames attributs, il n\u2019en est pas de m\u00eame du type de g\u00e9om\u00e9trie. L\u2019information (P,L,A) indique que, dans cette classe d\u2019objets on peut trouver des points, des lignes et des polygones. Contrairement aux standards SIG, les trois types de g\u00e9om\u00e9trie coexistent au sein de la m\u00eame classe d\u2019objets.<\/p>\n\n\n\n
Quel qu\u2019il soit le format choisi pour l\u2019int\u00e9grer dans QGis, nous aurons \u00e0 cr\u00e9er une couche par type de g\u00e9om\u00e9trie. Si on ne le fait pas, GDAL va cr\u00e9er le type de couche en fonction de la premi\u00e8re entit\u00e9 trouv\u00e9e lors de la conversion. Si elle est de type point, la couche cr\u00e9\u00e9e sera de type point et les entit\u00e9s lignes et polygones seront ignor\u00e9es. De m\u00eame, si la premi\u00e8re est de type ligne, ce seront les points et les polygones qui seront ignor\u00e9s.<\/p>\n\n\n\n
Il faut aussi noter que le format S57 n\u2019a aucune contrainte sur les classes d\u2019objets: si dans la zone couverte par le fichier S57 il n\u2019y a pas une classe d\u2019objet, il n\u2019y aura rien la concernant dans le fichier (pas de couche vide). De m\u00eame, s\u2019il n\u2019y a qu\u2019un type de g\u00e9om\u00e9trie pr\u00e9sente, alors que les trois types sont possibles, il n\u2019y aura pas de trace des autres types de g\u00e9om\u00e9trie.<\/p>\n\n\n\n
Il faudra alors d\u00e9composer le traitement d\u2019un fichier S57 avec ogr2ogr en trois \u00e9tapes, une pour chaque type de g\u00e9om\u00e9trie. Les options suivantes permettent de traiter chaque classe d\u2019objet S57 en ne s\u00e9lectionnant qu\u2019un type de g\u00e9om\u00e9trie :<\/p>\n\n\n\n
-where \u00ab OGR_GEOMETRY=\u2019POINT\u2019 or OGR_GEOMETRY=\u2019MULTIPOINT\u2019 \u00bb<\/strong><\/p>\n\n\n\n
-where \u00ab OGR_GEOMETRY=\u2019LINESTRING\u2019 or OGR_GEOMETRY=\u2019MULTILINESTRING\u2019 \u00bb<\/strong><\/p>\n\n\n\n
-where \u00ab OGR_GEOMETRY=\u2019POLYGON\u2019 or OGR_GEOMETRY=\u2019MULTIPOLYGON\u2019 \u00bb<\/strong><\/p>\n\n\n\n
GDAL permet pour certains types de drivers de cr\u00e9er des pr\u00e9fixes dans les tables en sortie. Dans ce cas on pourrait cr\u00e9er toutes les tables (points, lignes, polygones) dans un seul geopackage<\/a>, en les pr\u00e9fixant par exemple avec pt_, li_ et pl_. Le probl\u00e8me est que le driver S57 de GDAL ne permet pas cette option. On doit alors cr\u00e9er trois geopackages distincts, dans lesquels les tables seront cr\u00e9\u00e9es selon le type de g\u00e9om\u00e9trie. Dans chaque geopackage il y aura alors une table de m\u00eame nom mais de g\u00e9om\u00e9trie diff\u00e9rente. Si cette solution est la plus simple, elle n\u2019est pas la plus \u00e9l\u00e9gante et pratique. Nous allons utiliser ici trois geopackages d\u2019import pour les commandes ogr2ogr dans lesquels nous allons importer les tables d\u2019un fichier S57. Nous les appellerons pointsENC, LinesENC et PolysENC. D\u2019autre part nous allons cr\u00e9er un geopackage unique appel\u00e9 ENC<\/a> pour la base de donn\u00e9es compl\u00e8te. Pour chaque carte ENC nous allons importer son contenu dans les trois geopackages d\u2019import avec ogr2ogr puis nous allons ex\u00e9cuter des scripts Python pour mettre \u00e0 jour la base de donn\u00e9es du geopackage ENC avec les nouvelles tables import\u00e9es. Nous sommes oblig\u00e9s d\u2019utiliser les scripts Python au lieu de requ\u00eates et fonctions SQL car la version SQL de SQLite, utilis\u00e9e pour les geopackages, est tr\u00e8s limit\u00e9e.<\/p>\n\n\n\n
Pour le premier fichier S57 \u00e0 importer le processus est le suivant:<\/p>\n\n\n\n
\n
Cr\u00e9ation du geopackage ENC, vide.<\/li>\n\n\n\n
Chargement des classes d\u2019objets Point dans le geopackage pointsENC avec ogr2ogr<\/li>\n\n\n\n
Chargement des classes d\u2019objets Lignes dans le geopackage LinesENC avec ogr2ogr<\/li>\n\n\n\n
Chargement des classes d\u2019objets Polygone dans le geopackage PolysENC avec ogr2ogr<\/li>\n\n\n\n
Suppression des tables vides des trois geopackages.<\/li>\n\n\n\n
Clonage des tables des trois geopackages d\u2019import dans le geopackage ENC, en pr\u00e9fixant les tables de pointsENC par pt_, de LinesENC par li_ et celles de PolysENC par pl_.<\/li>\n<\/ol>\n\n\n\n
Pour les fichiers S57 suivants, le processus est le suivant:<\/p>\n\n\n\n
\n
Chargement des classes d\u2019objets Point dans le geopackage pointsENC avec ogr2ogr<\/li>\n\n\n\n
Chargement des classes d\u2019objets Lignes dans le geopackage LinesENC avec ogr2ogr<\/li>\n\n\n\n
Chargement des classes d\u2019objets Polygone dans le geopackage PolysENC avec ogr2ogr<\/li>\n\n\n\n
Suppression des tables vides des trois geopackages.<\/li>\n\n\n\n
Mise \u00e0 jour des tables existantes dans le geopackage ENC \u00e0 partir des trois geopackages d\u2019import ET<\/li>\n\n\n\n
Clonage des tables des trois geopackages d\u2019import qui n\u2019\u00e9taient pas pr\u00e9sentes dans le geopackage ENC,<\/li>\n\n\n\n
D\u00e9tection et effacement des doublons \u00e9ventuels r\u00e9sultant de la mise \u00e0 jour des tables du geopackage ENC<\/li>\n<\/ol>\n\n\n\n
Le traitement des sondes bathym\u00e9triques<\/h2>\n\n\n\n
Les sondes bathym\u00e9triques posent plusieurs probl\u00e8mes lors de la conversion de format. Le premier est que la valeur de profondeur n\u2019est pas une contenue dans un attribut, mais en tant que Z dans la g\u00e9om\u00e9trie (XYZ). Le deuxi\u00e8me c\u2019est que les sondes ne sont pas du type Point <\/em><\/strong>mais Multipoint<\/em><\/strong>. Pour avoir directement la valeur des sondes dans un champ attributaire, il faut donc ajouter deux param\u00e8tres \u00e0 la ligne de commande ogr2ogr:<\/p>\n\n\n\n
<\/div>
sondes bathym\u00e9triques<\/div>
<\/div><\/div>
<\/div><\/div>\n\n\n\n
R\u00e9cup\u00e9ration des identifiants \u00ab parents \u00bb<\/h2>\n\n\n\n
Certains objets peuvent \u00eatre regroup\u00e9s par un objet non g\u00e9ographique parent. Par exemple, les secteurs de feu sont cod\u00e9s avec un secteur par enregistrement de la table \u00ab LIGHTS \u00bb. Les diff\u00e9rents secteurs d\u2019un m\u00eame feu n\u2019ont aucun attribut particulier qui puisse indiquer qu\u2019ils correspondant \u00e0 la m\u00eame entit\u00e9. Cette information se trouve dans un enregistrement \u00ab parent \u00bb. Pour r\u00e9cup\u00e9rer cet identifiant en tant qu\u2019attribut de la table LIGTHS, il faut ajouter une option \u00e0 la ligne de commande:<\/p>\n\n\n\n
<\/div>
identifiants parents<\/div>
<\/div><\/div>
<\/div><\/div>\n\n\n\n
Cette option cr\u00e9e un attribut name_rcid <\/em><\/strong>qui est commun \u00e0 tous les secteurs d\u2019un m\u00eame feu, mais cr\u00e9\u00e9 aussi une s\u00e9rie de champs (name_rcnm,ORNT,USAG,MASK<\/em><\/strong>).<\/p>\n\n\n\n
Le traitement des type \u00ab Liste \u00bb<\/h2>\n\n\n\n
Le format S57 utilise, en plus des types de champ classiques (entier, r\u00e9el, texte) un type particulier: les listes, de cha\u00eenes de caract\u00e8res ou d\u2019entiers. Si on retrouve ces types de champ dans PostgreSQL, il n\u2019en va pas de m\u00eame dans les shapefiles et le geopackage.<\/p>\n\n\n\n
En absence des options sp\u00e9cifiques, il y aura des messages d\u2019avertissement indiquant que le format n\u2019est pas support\u00e9 et le contenu de ces champs sera ignor\u00e9. Pour ne pas perdre le contenu, il faut ajouter dans la ligne de commande ogr2ogr l\u2019option :<\/p>\n\n\n\n
<\/div>
listes<\/div>
<\/div><\/div>
<\/div><\/div>\n\n\n\n
Cette option prend une liste {..,\u2026} et g\u00e9n\u00e8re une cha\u00eene de texte avec un premier chiffre indiquant le nombre d\u2019\u00e9l\u00e9ments de la liste, deux points, puis les \u00e9l\u00e9ments de la liste s\u00e9par\u00e9s par des virgules:<\/p>\n\n\n\n
{3} -> \u2018(1:3)\u2019<\/p>\n\n\n\n
{3,1} -> \u2018(2:3,1)\u2019<\/p>\n\n\n\n
Forc\u00e9ment, ceci complique un peu le traitement de ces champs mais l\u2019information n\u2019est pas perdue. Par contre, le driver GDAL pr\u00e9sente un probl\u00e8me pour la gestion du format String r\u00e9sultant.<\/p>\n\n\n\n
Tout d\u2019abord, pour ce qui est des attributs des objets S57 il n\u2019y a aucun probl\u00e8me. Le champ de la table geopackage est de type \u2018Text\u2019, tout court, ce qui permet la cr\u00e9ation d\u2019une cha\u00eene de caract\u00e8res de n\u2019importe quelle longueur.
Par contre, pour les \u00ab identifiants parents \u00bb, la d\u00e9finition utilis\u00e9e par le driver est \u2018Text()\u2019, c\u2019est \u00e0 dire une cha\u00eene de caract\u00e8res avec une longueur maximum. Cette longueur est, a priori, d\u00e9termin\u00e9e par le contenu du premier enregistrement trait\u00e9. C\u2019est dire que ce n\u2019est que du hasard si vous avez un champ avec la longueur n\u00e9cessaire pour int\u00e9grer toutes les valeurs du fichier.<\/p>\n\n\n\n
Vous verrez alors des messages du type:<\/p>\n\n\n\n
WARNING : Value of field ‘NAME_RCID’ has 371 characters, whereas maximum allowed is 30<\/em><\/strong><\/p>\n\n\n\n
La solution n\u2019est pas compliqu\u00e9e mais implique pas mal de temps de travail. Partez de la base que pour l\u2019affichage classique des cartes ENC, ceci ne repr\u00e9sente pas un probl\u00e8me et vous pouvez tout simplement ignorer ces messages d\u2019avertissement.<\/p>\n\n\n\n
Par contre si vous avez besoin de traitements plus pouss\u00e9s incluant ces informations, voici le seul moyen que j\u2019ai trouv\u00e9, pour le moment:<\/p>\n\n\n\n
<\/div>
r\u00e9cup\u00e9ration des values<\/div>
<\/div><\/div>
<\/div><\/div>\n\n\n\n
<\/p>\n\n\n\n
Commandes ogr2ogr pour cr\u00e9er les geopackages d\u2019import<\/h2>\n\n\n\n
Fort de tous ces \u00e9l\u00e9ments, voici donc les trois lignes de commandes ogr2ogr pour cr\u00e9er les tables dans les geopackages d\u2019import:<\/p>\n\n\n\n
Commandes ogr2ogr<\/p>\n\n\n\n
<\/div>
Commandes ogr2ogr<\/div>
<\/div><\/div>
<\/div><\/div>\n\n\n\n
Pour ex\u00e9cuter ces lignes de commande il suffit d\u2019ouvrir la fen\u00eatre de shell d\u2019OSGeo4W:<\/p>\n\n\n\n
$\"Osgeo4W$ <\/a><\/figure>\n\n\n\n
La fen\u00eatre de Shell s\u2019ouvre<\/p>\n\n\n\n
$\"Osgeo4W$ <\/a><\/figure>\n\n\n\n
Rentrez alors les lignes de commande, en faisant attention \u00e0 bien rentrer le chemin et nom de fichier avec l\u2019extension .000 .<\/p>\n\n\n\n
On aura alors comme r\u00e9sultat une s\u00e9rie de tables cr\u00e9\u00e9es dans chaque geopackage d\u2019import. Par contre il y aura des tables compl\u00e8tement vides. En effet, si le type de g\u00e9om\u00e9trie est possible pour une classe d\u2019objet, la table sera cr\u00e9\u00e9e mais s\u2019il n\u2019y a aucune occurrence dans le fichier S57 trait\u00e9, la table n\u2019aura aucun enregistrement.<\/p>\n\n\n\n
Charger plusieurs fichiers S57 en m\u00eame temps<\/h2>\n\n\n\n
Les lignes de commandes du paragraphe pr\u00e9c\u00e9dent assument que vous avez un seul fichier S57 \u00e0 charger. Mais vous pouvez charger plusieurs fichiers simultan\u00e9ment. L\u2019utilisation des options -append -update font que vous pouvez traiter en boucle tous les fichiers .000 d\u2019un r\u00e9pertoire. Vous aurez ainsi tous les fichiers charg\u00e9s dans les geopackages d\u2019import.
Les fichiers .bat pour r\u00e9aliser cette t\u00e2che sont les suivants:<\/p>\n\n\n\n
Fichier .bat pour charger les g\u00e9om\u00e9tries point<\/p>\n\n\n\n
Fichier .bat pour charger les g\u00e9om\u00e9tries ligne<\/p>\n\n\n\n
Fichier .bat pour charger les g\u00e9om\u00e9tries polygone<\/p>\n\n\n\n
Faites attention \u00e0 modifier les chemins pour qu\u2019ils correspondent \u00e0 vos geopackages d\u2019import et au r\u00e9pertoire \u00e0 charger. On pourrait faire un seul fichier.bat pour les trois traitements, mais il convient de bien v\u00e9rifier le r\u00e9sultat de chacun des types de g\u00e9om\u00e9trie.<\/p>\n\n\n\n
Cr\u00e9ation du geopackage ENC<\/h2>\n\n\n\n
Nous allons cr\u00e9er un geopackage vide qui recevra les tables import\u00e9es, puis les mises \u00e0 jour successives lors des chargements d\u2019autres fichiers S57.<\/p>\n\n\n\n
Dans le panneau Explorateur <\/em><\/strong>de QGis, faites un clic droit sur Geopackages <\/em><\/strong>et s\u00e9lectionnez Cr\u00e9er une base de donn\u00e9es<\/em><\/strong><\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
Dans la fen\u00eatre qui s\u2019ouvre, renseignez le nom du geopackage, laissez l\u2019option par d\u00e9faut pour le nom de la table et s\u00e9lectionnez Pas de g\u00e9om\u00e9trie<\/em> dans Type de g\u00e9om\u00e9trie<\/em><\/strong>.<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
Cliquez sur OK pour avoir le nouveau geopackage vide ajout\u00e9 aux connexions geopackage de QGis. Il est conseill\u00e9 d\u2019effacer la table vide ENC qui est cr\u00e9\u00e9e lors de la cr\u00e9ation du geopackage. Pour cela ouvrez le gestionnaire de base de donn\u00e9es de QGis \u00e0 partir du menu principal, positionnez-vous sur la table ENC du geopackage ENC.gpkg et avec un clic droit s\u00e9lectionnez Effacer\u2026<\/em><\/strong><\/p>\n\n\n\n
Utilisation des scripts Python pour la gestion de la base de donn\u00e9es<\/h2>\n\n\n\n
Pour toutes les op\u00e9rations de gestion de la base de donn\u00e9es nous allons utiliser la console Python de QGis.<\/p>\n\n\n\n
Pour l\u2019ouvrir, allez dans le menu Extensions -> Console Python. Les panneaux de la console Python s\u2019ouvrent:<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
Pour ex\u00e9cuter un script Python, chargez le code dans le panneau \u00c9diteur<\/strong>. Vous pouvez faire un copier-coller ou bien utiliser l\u2019ic\u00f4ne R\u00e9pertoire <\/strong>pour pointer sur le fichier .py. L\u2019indentation \u00e9tant cruciale pour les scripts Python, privil\u00e9giez cette deuxi\u00e8me option.<\/p>\n\n\n\n
Suppression des tables vides et ajout du nom de la carte et de l’\u00e9chelle<\/h2>\n\n\n\n
Comme nous l’avons vu plus haute, lors de l’ex\u00e9cution des commandes ogr2ogr, des tables vides sont cr\u00e9\u00e9es pour les couches o\u00f9 il n’y a pas la g\u00e9om\u00e9trie demand\u00e9e dans la ligne de commande. Nous allons donc supprimer ces tables vides et profiter du script Python pour ajouter deux attributs \u00e0 toutes les tables: le nom de la carte ENC et l’\u00e9chelle de celle-ci.<\/p>\n\n\n\n
Pour ce qui est du nom, il peut \u00eatre utile dans la maintenance ult\u00e9rieure de la base de donn\u00e9es de pouvoir s\u00e9lectionner les entit\u00e9s correspondantes \u00e0 un fichier source S57. Cette information n’est contenue dans aucun attribut du fichier S57. Nous devons donc la renseigner manuellement.<\/p>\n\n\n\n
Il en va de m\u00eame pour l’\u00e9chelle de la carte. Si un attribut est bien pr\u00e9vu dans le S57 pour rentrer l’\u00e9chelle de la carte \u00ab\u00a0papier\u00a0\u00bb, l\u2019attribut CSCALE dans la table (M_CSCL, il est bien rare que cette table et son attribut soient pr\u00e9sents dans les fichiers S57.<\/p>\n\n\n\n
Si le nom de la carte est d’une utilit\u00e9 relative, il n’en va pas de m\u00eame de l’\u00e9chelle, car c’est un crit\u00e8re indispensable dans la recherche de doublons. En effet, lors du m\u00e9lange des entit\u00e9s en provenance de plusieurs fichiers S57 il y aura une coexistence plus ou moins heureuse de donn\u00e9es issues de diff\u00e9rentes \u00e9chelles.<\/p>\n\n\n\n
Le code Python suivant permet de supprimer les tables vides d\u2019un sch\u00e9ma d\u2019import et d’ajouter deux attributs (enc_chart et scale) \u00e0 toutes les tables:<\/p>\n\n\n\n
<\/div>
tables vides .py<\/div>
<\/div><\/div>
<\/div><\/div>\n\n\n\n
Attention \u00e0 respecter l\u2019indentation du code Python. Pour simplifier la t\u00e2che, t\u00e9l\u00e9chargez directement le fichier .py avec ce lien: delete_empty_tables_in_geopackages.py.<\/p>\n\n\n\n
Pour l’ex\u00e9cution, veillez \u00e0 modifier les param\u00e8tres en entr\u00e9e du script:<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
Les geopackage_paths vous n’aurez qu’\u00e0 les rentrer que la premi\u00e8re fois que vous ex\u00e9cuterez le script. Par contre le nom de la carte ENC et l’\u00e9chelle seront \u00e0 modifier lors de l’import de chaque nouveau fichier S57.
A l’issue de cette ex\u00e9cution vous aurez dans chaque table les nouveaux attributs:<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
Clonage et\/ou ajout des tables import\u00e9es dans le geopackage ENC<\/h2>\n\n\n\n
Dans le cas du premier chargement S57, le script Python suivant cr\u00e9\u00e9 toutes les tables pr\u00e9sentes dans les trois geopackages d\u2019import. Pour les chargements ult\u00e9rieurs, si la table existe d\u00e9j\u00e0 dans le geopackage ENC, les enregistrements de la table du geopackage d\u2019import sont ajout\u00e9s aux enregistrements d\u00e9j\u00e0 pr\u00e9sents dans ENC. Si, par contre, la table n\u2019existe pas encore, elle est cr\u00e9\u00e9e et aliment\u00e9e par les enregistrements du geopackage d\u2019import. Dans le fichier geopackage ENC les tables sont pr\u00e9fix\u00e9es par pt_, li_ et pl_ pour indiquer leur type de g\u00e9om\u00e9trie.<\/p>\n\n\n\n
Le clonage est r\u00e9alis\u00e9 avec le code suivant:<\/p>\n\n\n\n
<\/div>
Clonage ou mise \u00e0 jour .py<\/div>
<\/div><\/div>
<\/div><\/div>\n\n\n\n
Vous pouvez t\u00e9l\u00e9charger le fichier .py avec ce lien: clone_or_append_tables_with_prefix.py.<\/p>\n\n\n\n
Ajout de la symbologie au Geopackage<\/h2>\n\n\n\n
La symbologie par d\u00e9faut des tables d’un geopackage est contenue dans une table nomm\u00e9e layer_styles qui n’est cr\u00e9\u00e9e que si vous utilisez l’option enregistrer par d\u00e9faut du bouton style.<\/p>\n\n\n\n
Apr\u00e8s le clonage des tables dans le fichier ENC, la table n’existe pas. Nous allons donc importer la table contenant toute la symbologie des cartes ENC. Pour cela; t\u00e9l\u00e9chargez le geopackage layer_styles en cliquant ici.<\/p>\n\n\n\n
Ce geopackage ne contient qu’une table, nomm\u00e9e layer_styles, avec les symbologies de toutes les tables S57.<\/p>\n\n\n\n
Cette symbologie utilise des symboles svg que vous devez t\u00e9l\u00e9charger sur votre machine. En cliquant ici vous t\u00e9l\u00e9chargerez un r\u00e9pertoire ‘nautical’ contenant tous les symboles svg n\u00e9cessaires.<\/p>\n\n\n\n
Par d\u00e9faut, les r\u00e9f\u00e9rences aux symboles svg sont faites vers un r\u00e9pertoire ‘C:\/nautical’. Vous avez plusieurs possibilit\u00e9s:<\/p>\n\n\n\n
1- D\u00e9compresser le fichier nautical dans un r\u00e9pertoire c:\/nautical de votre machine<\/h3>\n\n\n\n
Dans ce cas vous n’avez plus rien \u00e0 faire. Les layer_styles trouveront les symboles svg sans probl\u00e8me.<\/p>\n\n\n\n
2- Vous souhaitez enregistrer le r\u00e9pertoire nautical dans un r\u00e9pertoire de votre choix<\/h3>\n\n\n\n
Dans ce cas, vous devez modifier la table layer_styles pour qu’elle trouve les symboles.<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
Ouvrez le gestionnaire de base de donn\u00e9es de QGis.<\/p>\n\n\n\n
\n
S\u00e9lectionnez la table layer_styles de votre geopackage<\/li>\n\n\n\n
Ouvrez une fen\u00eatre SQL<\/li>\n\n\n\n
Entrez la requ\u00eate:
UPDATE layer_styles
SET styleQML = REPLACE(styleQML, ‘C:\/nautical’<\/em>, ‘Votre chemin’<\/em>);<\/li>\n\n\n\n
Cliquez sur le bouton Ex\u00e9cuter<\/li>\n<\/ol>\n\n\n\n
3- Vous souhaitez enregistrer le r\u00e9pertoire nautical <\/em><\/strong>dans un le r\u00e9pertoire de votre geopackage ENC<\/h3>\n\n\n\n
Dans ce cas, les r\u00e9f\u00e9rences aux symboles svg seront relatives au Geopackage et vous pourrez l’utiliser \u00e0 partir de plusieurs machines. On fera alors la m\u00eame d\u00e9marche que pr\u00e9c\u00e9demment mais la requ\u00eate SQL sera:<\/p>\n\n\n\n
UPDATE layer_styles
SET styleQML = REPLACE(styleQML, ‘C:\/nautical’<\/em>, ‘<\/em>nautical‘<\/em>);<\/p>\n\n\n\n
Bien s\u00fbr, il faut que le projet soit configur\u00e9 pour stocker les chemins relatifs (Propri\u00e9t\u00e9s du projet -> G\u00e9n\u00e9ral->Enregistrer les chemins->Relatif<\/em><\/strong>)<\/p>\n\n\n\n
Gestion des doublons<\/h2>\n\n\n\n
Quand deux cartes ENC se superposent, les entit\u00e9s pr\u00e9sentes dans la zone de superposition vont \u00eatre ins\u00e9r\u00e9es autant de fois que de superpositions il y en a. Il faut donc, pour \u00e9viter des doublons et plus, nettoyer la base de donn\u00e9es apr\u00e8s chaque import de fichier S57.<\/p>\n\n\n\n
Mais il y a deux types de doublons diff\u00e9rents. D’une part, les zones de recouvrement de deux cartes d’\u00e9chelle similaire vont produire de \u00ab\u00a0vrais\u00a0\u00bb doublons qu’il faudra nettoyer. Mais la superposition de cartes d’\u00e9chelle tr\u00e8s diff\u00e9rente est \u00e0\u00e7 traiter diff\u00e9remment. La carte moins d\u00e9taill\u00e9 (petite \u00e9chelle) aura une s\u00e9lection des donn\u00e9es de la carte plus d\u00e9taill\u00e9 (grande \u00e9chelle). Dans ce deuxi\u00e8me cas il n’est pas conseill\u00e9 de supprimer cette information.<\/p>\n\n\n\n
Le code Python ci-apr\u00e8s va supprimer les \u00ab\u00a0vrais\u00a0\u00bb doublons. Pour ce qui est des autres doublons, il est conseill\u00e9 de travailler avec des filtres QGis. Selon vos besoins, vous pourrez d\u00e9finir, gr\u00e2ce \u00e0 l’attribut \u00ab\u00a0scale\u00a0\u00bb que nous avons ajout\u00e9 lors de la suppression des tables vides, quelles sont les donn\u00e9es affich\u00e9es ou pas sur votre carte.<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
Dans cet exemple, seules les entit\u00e9s correspondantes \u00e0 une \u00e9chelle comprise entre 12500 et 45000 seront affich\u00e9es sur votre carte.<\/p>\n\n\n\n
Le code suivant cr\u00e9e, pour chaque table du geopackage ENC, une table temporaire (`temp_table<\/code>). Puis il remplit cette table avec les enregistrements distincts bas\u00e9s sur name_rcid<\/code>. Ensuite, il supprime la table d\u2019origine et renomme la table temporaire avec le nom de la table d\u2019origine.<\/p>\n\n\n\n`
Ce traitement s\u2019applique \u00e0 toutes les tables S57, sauf les tables des feux (LIGHTS) et des sondes bathym\u00e9triques (SOUNDG). En effet, chaque secteur de feu a le m\u00eame name_rcid<\/code> mais correspond \u00e0 une couleur est un secteur g\u00e9ographique diff\u00e9rents. De m\u00eame, le name_rcid<\/code> des sondes bathym\u00e9triques est le m\u00eame pour tous les points de sondes de l\u2019entit\u00e9 multipoint d\u2019origine. Dans le premier cas nous devons chercher les couleurs et secteurs en double pour un m\u00eame name_rcid<\/code>, et pour les sondes nous devons trouver les sondes ayant les m\u00eames XY dans leur g\u00e9om\u00e9trie.<\/p>\n\n\n\n
<\/div>effacer doublons .py<\/div><\/div><\/div> <\/div><\/div>\n\n\n\n<\/p>\n\n\n\n Vous pouvez t\u00e9l\u00e9charger le fichier .py \u00e0 partir de ce lien : supprime_doublons.py. N\u2019oubliez pas de modifier les chemins et noms de fichiers pour qu\u2019ils correspondent aux v\u00f4tres. Remarquez aussi la ligne :<\/p>\n\n\n\n if table_name not in [\u00ab pt_LIGHTS \u00bb, \u00ab pt_SOUNDG \u00bb, \u00bblayer_styles \u00bb]:<\/em><\/strong><\/p>\n\n\n\n En principe, vous aurez que les tables pr\u00e9fix\u00e9es et la table layer_styles. Cette ligne permet de traiter diff\u00e9remment les tables LIGTHS et SOUNDG, mais surtout, elle emp\u00eache de vider la table layer_styles<\/em><\/strong>. Si vous \u00eates amen\u00e9 \u00e0 cr\u00e9er d\u2019autres tables que celles issues du fichier s57 dans le geopackage, il faut prendre en compte les noms de ces tables dans ce script, en les ajoutant aux exceptions.<\/p>\n\n\n\n <\/p>\n","protected":false},"excerpt":{"rendered":" Cet article sous sa forme compl\u00e8te (codes et liens) sera un des r\u00e9sultats du travail r\u00e9sultant du financement collaboratif que nous avons mis en place: https:\/\/www.sigterritoires.fr\/index.php\/financement-collaboratif-pour-lintegration-de-donnees-marines-dans-qgis\/. D\u00e8s que le montant vis\u00e9 sera atteint, l’article sera mis…<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"give_campaign_id":0,"_bbp_topic_count":0,"_bbp_reply_count":0,"_bbp_total_topic_count":0,"_bbp_total_reply_count":0,"_bbp_voice_count":0,"_bbp_anonymous_reply_count":0,"_bbp_topic_count_hidden":0,"_bbp_reply_count_hidden":0,"_bbp_forum_subforum_count":0,"sfsi_plus_gutenberg_text_before_share":"","sfsi_plus_gutenberg_show_text_before_share":"","sfsi_plus_gutenberg_icon_type":"","sfsi_plus_gutenberg_icon_alignemt":"","sfsi_plus_gutenburg_max_per_row":"","_monsterinsights_skip_tracking":false,"_monsterinsights_sitenote_active":false,"_monsterinsights_sitenote_note":"","_monsterinsights_sitenote_category":0,"jetpack_post_was_ever_published":false,"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_memberships_contains_paid_content":false,"footnotes":""},"categories":[425],"tags":[],"class_list":["post-12338","post","type-post","status-publish","format-standard","hentry","category-non-classe-fr"],"aioseo_notices":[],"jetpack_featured_media_url":"","jetpack_shortlink":"https:\/\/wp.me\/p6XU0A-3d0","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/12338","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/comments?post=12338"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/12338\/revisions"}],"wp:attachment":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media?parent=12338"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=12338"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=12338"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}