<\/div><\/div>

<\/p>\n\n\n\n

-oo SPLIT_MULTIPOINT=ON -oo ADD_SOUNDG_DEPTH=ON<\/em><\/strong><\/p>\n\n\n\n

<\/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>
<\/p>\n\n\n\n
-oo RETURN_LINKAGES=O<\/strong><\/p>\n\n\n\n
<\/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>
<\/p>\n\n\n\n
-mapFieldType StringList=String,IntegerList=String<\/em><\/strong><\/p>\n\n\n\n
<\/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>
<\/p>\n\n\n\n
\n
Ex\u00e9cutez la commande ogr2ogr pour cr\u00e9er le geopackage d’import<\/li>\n\n\n\n
Rep\u00e9rez la valeur maximale indiqu\u00e9e pour chaque type de champ<\/li>\n\n\n\n
Pour la table ou les tables en sortie qui auront \u00e9t\u00e9 cr\u00e9\u00e9es et dont vous souhaitez r\u00e9cup\u00e9rer toute l’information, ouvrez le gestionnaire de base de donn\u00e9es de QGis<\/li>\n\n\n\n
Effacez tous les enregistrements de la table (DELETE FROM nom_de la table)<\/li>\n\n\n\n
Allez dans le menu Table->Modifiez une table…<\/em><\/strong>
S\u00e9lectionnez dans la liste des champs le champ qui vous concerne, puis cliquez sur Editer une colonne<\/em><\/strong><\/li>\n\n\n\n
Rentrez la valeur souhait\u00e9e dans le champ Longueur, puis sur OK<\/li>\n\n\n\n
Faites la m\u00eame op\u00e9ration pour tous les autres champs qui vous int\u00e9ressent<\/li>\n\n\n\n
Effacez du geopackage d’import toutes les autres tables qui ne posent pas de probl\u00e8me <\/li>\n\n\n\n
R\u00e9-ex\u00e9cutez la commande ogr2ogr d’import.<\/li>\n<\/ol>\n\n\n\n
<\/div><\/div>\n\n\n\n
Workflow pour le chargement d’un seul fichier S57<\/h1>\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
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
<\/div>
Commandes ogr2ogr<\/div>
<\/div><\/div>
<\/p>\n\n\n\n
ogr2ogr -f GPKG -skipfailures -append -update -where \"OGR_GEOMETRY='POINT' or OGR_GEOMETRY='MULTIPOINT'\" -oo SPLIT_MULTIPOINT=ON -oo RETURN_LINKAGES=ON -oo LNAM_REFS=ON -oo ADD_SOUNDG_DEPTH=ON -mapFieldType StringList=String,IntegerList=String C:\/testgpkg\/pointsENC.gpkg C:\\testgpkg\\US1GC09M\\US1GC09M.000\n\nogr2ogr -f GPKG -skipfailures -append -update -where \"OGR_GEOMETRY='LINESTRING' or OGR_GEOMETRY='MULTILINESTRING'\" -oo RETURN_LINKAGES=ON -oo LNAM_REFS=ON -mapFieldType StringList=String,IntegerList=String C:\/testgpkg\/LinesENC.gpkg C:\\testgpkg\\US1GC09M\\US1GC09M.000\n\nogr2ogr -f GPKG -skipfailures -append -update -where \"OGR_GEOMETRY='POLYGON' or OGR_GEOMETRY='MULTIPOLYGON'\" -oo RETURN_LINKAGES=ON -oo LNAM_REFS=ON -mapFieldType StringList=String,IntegerList=String C:\/testgpkg\/PolysENC.gpkg C:\\testgpkg\\US1GC09M\\US1GC09M.000<\/code><\/pre>\n\n\n\n <\/div><\/div>\n\n\n\n Vous devrez modifier le texte C:\/testgpkg\/pointsENC.gpkg C:\\testgpkg\\US1GC09M\\US1GC09M.000 <\/p>\n\n\n\n par l’emplacement des geopackages d’import que vous souhaitez et par le chemin et nom de fichier du fichier S547 \u00e0 charger.<\/p>\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 <\/a><\/figure>\n\n\n\nLa fen\u00eatre de Shell s\u2019ouvre<\/p>\n\n\n\n<\/a><\/figure>\n\n\n\nRentrez 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\nUtilisation des scripts Python pour la gestion de la base de donn\u00e9es<\/h2>\n\n\n\nPour 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<\/em><\/strong>. Les panneaux de la console Python s\u2019ouvrent:<\/p>\n\n\n\n<\/a><\/figure>\n\n\n\nPour 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\nSuppression des tables vides et ajout du nom de la carte et de l’\u00e9chelle<\/h2>\n\n\n\nComme nous l’avons vu plus haut, 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> <\/p>\n\n\n\nfrom osgeo import ogr<\/code><\/pre>\n\n\n\ndef delete_empty_tables_in_geopackages(geopackage_paths, valeur_par_defaut_enc, valeur_par_defaut_scale): for geopackage_path in geopackage_paths: delete_empty_tables(geopackage_path, valeur_par_defaut_enc, valeur_par_defaut_scale) print(f\"Fin du traitement\") def delete_empty_tables(geopackage_path, valeur_par_defaut_enc, valeur_par_defaut_scale): # Ouvrir le GeoPackage en mode \u00e9dition geopackage = ogr.Open(geopackage_path, 1)<\/code><\/pre>\n\n\n\nif geopackage is not None:\n # R\u00e9cup\u00e9rer le nombre de tables dans le GeoPackage\n num_tables = geopackage.GetLayerCount()\n\n # Liste des tables \u00e0 supprimer\n tables_to_delete = []\n tables_to_update = []\n\n # Identifier les tables \u00e0 supprimer\n for i in range(num_tables):\n table = geopackage.GetLayerByIndex(i)\n\n # V\u00e9rifier si la table est vide (aucun enregistrement)\n if table.GetFeatureCount() == 0:\n tables_to_delete.append(table.GetName())\n else:\n tables_to_update.append(table.GetName())\n\n # Supprimer les tables\n for table_name in tables_to_delete:\n geopackage.DeleteLayer(table_name)\n print(f\"Table supprim\u00e9e dans {geopackage_path}: {table_name}\")\n\n # Mettre \u00e0 jour les tables restantes\n for table_name in tables_to_update:\n table = geopackage.GetLayerByName(table_name)\n\n # Ajouter le premier champ texte\n champ1 = ogr.FieldDefn('enc_chart', ogr.OFTString)\n champ1.SetWidth(50)\n champ1.SetNullable(True)\n champ1.SetDefault(valeur_par_defaut_enc)\n table.CreateField(champ1)\n\n # Ajouter le deuxi\u00e8me champ texte\n champ2 = ogr.FieldDefn('scale', ogr.OFTReal)\n champ2.SetWidth(10)\n champ2.SetPrecision(2)\n champ2.SetDefault(valeur_par_defaut_scale)\n table.CreateField(champ2)\n\n # Mettre \u00e0 jour les valeurs dans les champs apr\u00e8s leur cr\u00e9ation\n for feature in table:\n feature.SetField('enc_chart', valeur_par_defaut_enc)\n feature.SetField('scale', valeur_par_defaut_scale)\n table.SetFeature(feature)\n print(f\"ENC & scale ajout\u00e9s \u00e0 {geopackage_path}: {table_name}\")\n\n # Fermer le GeoPackage\n geopackage = None\n\nelse:\n print(f\"Impossible d'ouvrir le GeoPackage {geopackage_path} en mode \u00e9dition.\")<\/code><\/pre>\n\n\n\n#Lignes \u00e0 \u00e9diter avant l'ex\u00e9cution du script<\/code><\/pre>\n\n\n\ngeopackage_paths = [\"c:\/testgpkg\/pointsENC.gpkg\", \"c:\/testgpkg\/LinesENC.gpkg\", \"c:\/testgpkg\/PolysENC.gpkg\"] valeur_par_defaut_enc = 'US1GC09M' valeur_par_defaut_scale = '2160000' delete_empty_tables_in_geopackages(geopackage_paths, valeur_par_defaut_enc, valeur_par_defaut_scale)<\/code><\/pre>\n\n\n\n <\/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.<\/a><\/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\nLes 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\nWorkflow pour le chargement de plusieurs fichiers S57 en m\u00eame temps<\/h1>\n\n\n\nLes lignes de commandes des paragraphes pr\u00e9c\u00e9dents 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.<\/p>\n\n\n\n Contrairement au chargement d’un seul fichier, lors du chargement d’un lot de fichiers il faut int\u00e9grer l’ajout du nom du fichier et de l’\u00e9chelle \u00e0 la boucle de transcodage S57->geopackage. Pour cela, nous utiliserons un script python qui ajoute les attributs enc_chart et scale aux tables cr\u00e9\u00e9es par la commande ogr2ogr et qui les renseigne avec le nom du fichier correspondant et une \u00e9chelle d\u00e9finie par le producteur. Le script efface aussi les tables vides pr\u00e9sentes dans le geopackage d’import.<\/p>\n\n\n\n Merci \u00e0 Janez AVSEC pour l’id\u00e9e d’utiliser la table DSID pour r\u00e9cup\u00e9rer l’\u00e9chelle de compilation des donn\u00e9es.<\/p>\n\n\n\n<\/div>Mise \u00e0 jour disponible<\/div><\/div><\/div><\/p>\n\n\n\nUne mise \u00e0 jour de cette proc\u00e9dure est disponible dans l’article ici.<\/a><\/p>\n\n\n\n <\/div><\/div>\n\n\n\n Les fichiers .bat doivent \u00eatre ex\u00e9cut\u00e9s \u00e0 partir d’une fen\u00eatre OSGeo4W pour pouvoir b\u00e9n\u00e9ficier de l’environnement Python de QGis.<\/p>\n\n\n\n L’appel de chaque fichier .bat est de la forme:<\/p>\n\n\n\n .\\fichier.bat \u00ab\u00a0r\u00e9pertoire contenant les fichiers S57 \u00e0 charger\u00a0\u00bb \u00ab\u00a0nom du geopackage d’import\u00a0\u00bb <\/em><\/strong><\/p>\n\n\n\n Par exemple:<\/p>\n\n\n\n .\\mass_load_s57_polys.bat \u00ab\u00a0c:\/enc_db\/enc000\u00a0\u00bb \u00ab\u00a0c:\/enc_db\/polys.gpkg\u00a0\u00bb <\/em><\/strong> Les fichiers .bat pour r\u00e9aliser cette t\u00e2che sont les suivants:<\/p>\n\n\n\n Pour importer les tables de type point:<\/p>\n\n\n\n <\/div>Fichier .bat pour charger les g\u00e9om\u00e9tries point<\/div><\/div><\/div> <\/p>\n\n\n\n@echo off setlocal enabledelayedexpansion<\/code><\/pre>\n\n\n\nREM Assurez-vous qu'au moins 2 arguments ont \u00e9t\u00e9 fournis if \"%~2\"==\"\" ( echo Usage: %0 directory000 output_geopackage scale exit \/b 1 )<\/code><\/pre>\n\n\n\nREM R\u00e9cup\u00e9rer les arguments de la ligne de commande set \"directory=%~1\" set \"output_geopackage=%~2\"<\/code><\/pre>\n\n\n\nREM It\u00e9rez sur tous les fichiers .000 dans le r\u00e9pertoire for \/r \"%directory%\" %%i in (.000) do ( echo Traitement du fichier: %%~ni set \"file=%%~ni\" ogr2ogr -f GPKG -skipfailures -append -update -where \"OGR_GEOMETRY='POINT' or OGR_GEOMETRY='MULTIPOINT'\" -oo SPLIT_MULTIPOINT=ON -oo RETURN_LINKAGES=ON -oo LNAM_REFS=ON -oo ADD_SOUNDG_DEPTH=ON -mapFieldType StringList=String,IntegerList=String \"%output_geopackage%\" \"%%i\" ogr2ogr -f GPKG -skipfailures -append -update \"%output_geopackage%\" \"%%i\" \"DSID\" ogr2ogr -f GPKG -skipfailures -append -update \"%output_geopackage%\" \"%%i\" \"C_AGGR\" python c:\/testgpkg\/update_geopackage_dsid.py \"%output_geopackage%\" \"!file!\" )<\/code><\/pre>\n\n\n\necho Traitement termin\u00e9. pause<\/code><\/pre>\n\n\n\n <\/div><\/div>\n\n\n\n Vous remarquerez que, en plus de la commande ogr2ogr pour importer les couches de points, on importe deux tables sans g\u00e9om\u00e9trie: DSID et C_AGGR. <\/p>\n\n\n\n Pour importer les tables de type ligne:<\/p>\n\n\n\n <\/div>Fichier .bat pour charger les g\u00e9om\u00e9tries ligne<\/div><\/div><\/div> <\/p>\n\n\n\n@echo off setlocal enabledelayedexpansion<\/code><\/pre>\n\n\n\nREM Assurez-vous qu'au moins 2 arguments ont \u00e9t\u00e9 fournis\nif \"%~2\"==\"\" (\necho Usage: %0 directory000 output_geopackage scale\nexit \/b 1\n)<\/code><\/pre>\n\n\n\nREM R\u00e9cup\u00e9rer les arguments de la ligne de commande set \"directory=%~1\" set \"output_geopackage=%~2\" set \"scale=%~3\"<\/code><\/pre>\n\n\n\nREM It\u00e9rez sur tous les fichiers .000 dans le r\u00e9pertoire\nfor \/r \"%directory%\" %%i in (.000) do (\necho Traitement du fichier: %%~ni\nset \"file=%%~ni\"\nogr2ogr -f GPKG -skipfailures -append -update -where \"OGR_GEOMETRY='LINESTRING' or OGR_GEOMETRY='MULTILINESTRING'\" -oo RETURN_LINKAGES=ON -oo LNAM_REFS=ON -mapFieldType StringList=String,IntegerList=String \"%output_geopackage%\" \"%%i\"\nogr2ogr -f GPKG -skipfailures -append -update \"%output_geopackage%\" \"%%i\" \"DSID\"\nogr2ogr -f GPKG -skipfailures -append -update \"%output_geopackage%\" \"%%i\" \"C_AGGR\"\npython c:\/testgpkg\/update_geopackage_dsid.py \"%output_geopackage%\" \"!file!\" \n)<\/code><\/pre>\n\n\n\necho Traitement termin\u00e9. pause<\/code><\/pre>\n\n\n\n<\/div><\/div>\n\n\n\n Pour importer les tables de type polygone:<\/p>\n\n\n\n <\/div>Fichier .bat pour charger les g\u00e9om\u00e9tries polygone<\/div><\/div><\/div> <\/p>\n\n\n\n@echo off setlocal enabledelayedexpansion<\/code><\/pre>\n\n\n\nREM Assurez-vous qu'au moins 2 arguments ont \u00e9t\u00e9 fournis\nif \"%~2\"==\"\" (\necho Usage: %0 directory000 output_geopackage scale\nexit \/b 1\n)<\/code><\/pre>\n\n\n\nREM R\u00e9cup\u00e9rer les arguments de la ligne de commande set \"directory=%~1\" set \"output_geopackage=%~2\" set \"scale=%~3\"<\/code><\/pre>\n\n\n\nREM It\u00e9rez sur tous les fichiers .000 dans le r\u00e9pertoire\nfor \/r \"%directory%\" %%i in (.000) do (\necho Traitement du fichier: %%~ni\nset \"file=%%~ni\"\nogr2ogr -f GPKG -skipfailures -append -update -where \"OGR_GEOMETRY='POLYGON' or OGR_GEOMETRY='MULTIPOLYGON'\" -oo RETURN_LINKAGES=ON -oo LNAM_REFS=ON -mapFieldType StringList=String,IntegerList=String \"%output_geopackage%\" \"%%i\"\nogr2ogr -f GPKG -skipfailures -append -update \"%output_geopackage%\" \"%%i\" \"DSID\"\nogr2ogr -f GPKG -skipfailures -append -update \"%output_geopackage%\" \"%%i\" \"C_AGGR\"\npython c:\/testgpkg\/update_geopackage_dsid.py \"%output_geopackage%\" \"!file!\" \n)<\/code><\/pre>\n\n\n\necho Traitement termin\u00e9. pause<\/code><\/pre>\n\n\n\n <\/div><\/div>\n\n\n\n Les fichiers .bat contiennent la commande ogr2ogr ainsi qu’une ligne pour ex\u00e9cuter un script Python, update_geopackage_dsid.py. Faites attention \u00e0 modifier le chemin de ce script dans le fichier .bat pour qu’il corresponde \u00e0 celui de votre choix. Ce script ajoute \u00e0 toutes les tables du geopackage les deux attributs enc_chart <\/strong><\/em>et scale<\/em><\/strong>.<\/p>\n\n\n\n Le code python utilis\u00e9 est le suivant:<\/p>\n\n\n\n <\/div>update_geopackage_dsid.py<\/div><\/div><\/div> <\/p>\n\n\n\nimport sys from osgeo import ogr<\/p>\n\n\n\n def get_default_scale(geopackage_path, default_enc): # Ouvrir le GeoPackage en mode lecture seule geopackage = ogr.Open(geopackage_path, 0) default_scale = ‘0’ default_enc = default_enc +’.000′<\/p>\n\n\n\n if geopackage is not None:\n # R\u00e9cup\u00e9rer la couche DSID\n dsid_layer = geopackage.GetLayerByName('DSID')\n\n # V\u00e9rifier si la couche DSID existe\n if dsid_layer:\n # Rechercher la valeur de DSID_CSCL correspondant \u00e0 DSID_DSNM\n dsid_layer.SetAttributeFilter(f\"DSID_DSNM = '{default_enc}'\")\n feature = dsid_layer.GetNextFeature()\n if feature:\n default_scale = feature.GetField('DSPM_CSCL')\n\nreturn default_scale<\/code><\/pre>\n\n\n\ndef delete_empty_tables_in_geopackages_dsid(geopackage_path, default_enc): # R\u00e9cup\u00e9rer la valeur de scale \u00e0 partir de la table DSID default_scale = get_default_scale(geopackage_path, default_enc)<\/p>\n\n\n\n # Continuer avec le reste du code\nif default_scale is None:\n print(\"La valeur par d\u00e9faut de l'\u00e9chelle n'a pas \u00e9t\u00e9 trouv\u00e9e dans la table DSID.\")\n return\ndefault_scale_str = str(default_scale)\n# Ouvrir le GeoPackage en mode \u00e9dition\ngeopackage = ogr.Open(geopackage_path, 1)\n\nif geopackage is not None:\n # R\u00e9cup\u00e9rer le nombre de tables dans le GeoPackage\n num_tables = geopackage.GetLayerCount()\n\n # Liste des tables \u00e0 supprimer\n tables_to_delete = []\n tables_to_update = []\n\n # Identifier les tables \u00e0 supprimer\n for i in range(num_tables):\n table = geopackage.GetLayerByIndex(i)\n\n # V\u00e9rifier si la table est vide (aucun enregistrement)\n if table.GetFeatureCount() == 0:\n tables_to_delete.append(table.GetName())\n else:\n tables_to_update.append(table.GetName())\n\n # Supprimer les tables\n for table_name in tables_to_delete:\n geopackage.DeleteLayer(table_name)\n print(f\"Table supprim\u00e9e dans {geopackage_path}: {table_name}\")\n\n # Mettre \u00e0 jour les tables restantes\n for table_name in tables_to_update:\n table = geopackage.GetLayerByName(table_name)\n\n # V\u00e9rifier si le champ 'enc_chart' existe\n enc_chart_index = table.FindFieldIndex('enc_chart', 1)\n if enc_chart_index < 0:\n # Ajouter le champ 'enc_chart' s'il n'existe pas\n champ1 = ogr.FieldDefn('enc_chart', ogr.OFTString)\n champ1.SetWidth(50)\n champ1.SetNullable(True)\n champ1.SetDefault(default_enc)\n table.CreateField(champ1)\n\n # V\u00e9rifier si le champ 'scale' existe\n scale_index = table.FindFieldIndex('scale', 1)\n if scale_index < 0:\n # Ajouter le champ 'scale' s'il n'existe pas\n champ2 = ogr.FieldDefn('scale', ogr.OFTReal)\n champ2.SetWidth(10)\n champ2.SetPrecision(2)\n champ2.SetDefault(default_scale_str) # Convertir en nombre flottant\n table.CreateField(champ2)\n # Mettre \u00e0 jour les valeurs dans les champs seulement si les champs n'ont pas d\u00e9j\u00e0 de valeurs\n for feature in table:\n enc_chart_value = feature.GetField('enc_chart')\n scale_value = feature.GetField('scale')\n\n # V\u00e9rifier si les champs n'ont pas d\u00e9j\u00e0 de valeurs d\u00e9finies\n if enc_chart_value is None or scale_value is None:\n feature.SetField('enc_chart', default_enc)\n feature.SetField('scale', default_scale)\n table.SetFeature(feature)\n\n\n # Fermer le GeoPackage\n geopackage = None\n\nelse:\n print(f\"Impossible d'ouvrir le GeoPackage {geopackage_path} en mode \u00e9dition.\")<\/code><\/pre>\n\n\n\nif name<\/strong> == \u00ab\u00a0main<\/strong>\u00ab\u00a0: if len(sys.argv) != 3: print(\u00ab\u00a0Usage: python script.py geopackage_path default_enc\u00a0\u00bb) sys.exit(1)<\/p>\n\n\n\n geopackage_path = sys.argv[1]\ndefault_enc = sys.argv[2]\n\ndelete_empty_tables_in_geopackages_dsid(geopackage_path, default_enc)<\/code><\/pre>\n\n\n\n<\/div><\/div>\n\n\n\n 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 Vous pouvez t\u00e9l\u00e9charger les fichiers .bat et le code Python en cliquant ici.<\/a><\/p>\n\n\n\n Suite quel qu’il soit le nombre de fichiers charg\u00e9s<\/h1>\n\n\n\nClonage et\/ou ajout des tables import\u00e9es dans le geopackage ENC<\/h2>\n\n\n\nDans 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> <\/p>\n\n\n\nfrom osgeo import ogr<\/code><\/pre>\n\n\n\ndef clone_or_append_tables_with_prefix(): input_geopackages = [\"c:\/testgpkg\/pointsENC.gpkg\", \"c:\/testgpkg\/LinesENC.gpkg\", \"c:\/testgpkg\/PolysENC.gpkg\"] prefixes = [\"pt_\", \"li_\", \"pl_\"]<\/code><\/pre>\n\n\n\nfor i in range(len(input_geopackages)):\n # Ouvrir le GeoPackage source en mode lecture\n input_gpkg = ogr.Open(input_geopackages[i], 0)\n\n if input_gpkg is not None:\n # Ouvrir le GeoPackage de destination en mode \u00e9criture\n output_gpkg = ogr.Open(\"c:\/testgpkg\/ENC.gpkg\", 1)\n prefix = prefixes[i]\n\n if output_gpkg is not None:\n # R\u00e9cup\u00e9rer le nombre de tables dans le GeoPackage source\n num_tables = input_gpkg.GetLayerCount()\n\n # Ajouter le contenu des tables avec le pr\u00e9fixe appropri\u00e9\n for j in range(num_tables):\n input_table = input_gpkg.GetLayerByIndex(j)\n output_table_name = f\"{prefix}{input_table.GetName()}\"\n input_table_name = input_table.GetName()\n\n # V\u00e9rifier si la table de destination existe d\u00e9j\u00e0\n output_table = output_gpkg.GetLayerByName(output_table_name)\n if output_table is None:\n # Cr\u00e9er une nouvelle couche de destination si elle n'existe pas\n output_table = output_gpkg.CreateLayer(output_table_name, geom_type=input_table.GetGeomType(), options=[\"OVERWRITE=YES\"])\n\n # Copiez la structure de la table d'origine vers la table temporaire\n input_layer = input_gpkg.GetLayerByName(input_table_name)\n input_layer_defn = input_layer.GetLayerDefn()\n\n for k in range(input_layer_defn.GetFieldCount()):\n field_defn = input_layer_defn.GetFieldDefn(k)\n output_table.CreateField(field_defn)\n\n # Copier les entit\u00e9s de la table source vers la table de destination\n for feature in input_table:\n output_feature = ogr.Feature(output_table.GetLayerDefn())\n output_feature.SetGeometry(feature.GetGeometryRef())\n\n # Copier les valeurs des champs\n for field_index in range(feature.GetFieldCount()):\n output_feature.SetField(field_index, feature.GetField(field_index))\n\n output_table.CreateFeature(output_feature)\n output_feature = None # Lib\u00e9rer la m\u00e9moire\n\n print(f\"Contenu ajout\u00e9 de {input_geopackages[i]}.{input_table.GetName()} vers ENC.gpkg.{output_table_name}\")\n\n # Fermer les GeoPackages\n input_gpkg = None\n output_gpkg = None\n else:\n print(f\"Impossible d'ouvrir le GeoPackage {output_geopackage} en mode \u00e9criture.\")\n else:\n print(f\"Impossible d'ouvrir le GeoPackage {input_geopackages[i]} en mode lecture.\")<\/code><\/pre>\n\n\n\n#Appel de la fonction append_tables_with_prefix<\/code><\/pre>\n\n\n\nclone_or_append_tables_with_prefix()<\/code><\/pre>\n\n\n\n <\/div><\/div>\n\n\n\n Veillez \u00e0 modifier la ligne <\/p>\n\n\n\n input_geopackages = [\u00ab\u00a0c:\/testgpkg\/pointsENC.gpkg\u00a0\u00bb, \u00ab\u00a0c:\/testgpkg\/LinesENC.gpkg\u00a0\u00bb, \u00ab\u00a0c:\/testgpkg\/PolysENC.gpkg\u00a0\u00bb]<\/em><\/strong><\/p>\n\n\n\n avec les geopackages d’import que vous avez choisi, ainsi que la ligne<\/p>\n\n\n\n output_gpkg = ogr.Open(\u00ab\u00a0c:\/testgpkg\/ENC.gpkg\u00a0\u00bb, 1)<\/em><\/strong><\/p>\n\n\n\n avec le geopackage de base de donn\u00e9es de votre choix.<\/p>\n\n\n\n Vous pouvez t\u00e9l\u00e9charger le fichier .py avec ce lien: clone_or_append_tables_with_prefix.py.<\/a><\/p>\n\n\n\n Gestion des doublons<\/h2>\n\n\n\nQuand 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\nDans 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 Cette op\u00e9ration correspond au filtrage d’une couche, mais… que faire quand vous en avez 130 couches affich\u00e9es? Pour vous \u00e9viter de faire 130 fois cette op\u00e9ration, voici un code Python qui vous permet de filtrer toutes les couches du projet:<\/p>\n\n\n\n<\/div>filter<br>.py<\/div><\/div><\/div> <\/p>\n\n\n\n#Importer le module QGIS<\/code><\/pre>\n\n\n\nfrom qgis.core import QgsProject<\/code><\/pre>\n\n\n\n#D\u00e9finir la valeur de l'attribut \"scale\" \u00e0 filtrer<\/code><\/pre>\n\n\n\nvaleur_echelle = 50000<\/code><\/pre>\n\n\n\n#Obtenir le projet actif<\/code><\/pre>\n\n\n\nprojet = QgsProject.instance()<\/code><\/pre>\n\n\n\n#Obtenir la liste des couches charg\u00e9es dans le projet<\/code><\/pre>\n\n\n\ncouches = projet.mapLayers().values()<\/code><\/pre>\n\n\n\n#Parcourir toutes les couches<\/code><\/pre>\n\n\n\nfor couche in couches: # V\u00e9rifier si la couche a un champ nomm\u00e9 \"scale\" if couche.fields().indexFromName('scale') != -1: # D\u00e9finir le filtre sur l'attribut \"scale\" \u00e9gal \u00e0 50000 filtre = f'\"scale\" = {valeur_echelle}' couche.setSubsetString(filtre) print(f\"Filre d\u00e9fini pour la couche {couche.name()}\") else: print(f\"La couche {couche.name()} n'a pas d'attribut 'scale'\") <\/code><\/pre>\n\n\n\n<\/div><\/div>\n\n\n\n Et un autre pour enlever le filtre \u00e0 toutes les couches affich\u00e9es:<\/p>\n\n\n\n <\/div>unfilter<br>.py<\/div><\/div><\/div> <\/p>\n\n\n\n# Importer le module QGIS\nfrom qgis.core import QgsProject\n\n# Obtenir le projet actif\nprojet = QgsProject.instance()\n\n# Obtenir la liste des couches charg\u00e9es dans le projet\ncouches = projet.mapLayers().values()\n\n# Parcourir toutes les couches\nfor couche in couches:\n # V\u00e9rifier si la couche est une couche vectorielle\n if couche.type() == QgsMapLayer.VectorLayer:\n # Supprimer le filtre de la couche\n couche.setSubsetString('')\n print(f\"Filtre supprim\u00e9 pour la couche {couche.name()}\")\n<\/code><\/pre>\n\n\n\n<\/div><\/div>\n\n\n\n T\u00e9l\u00e9charger les codes python ici.<\/a><\/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 LNAM. 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:<\/p>\n\n\n\n \naux tables des feux (LIGHTS) et des sondes bathym\u00e9triques (SOUNDG) o\u00f9 l’attribut LNAM n’est pas utile. En effet, chaque secteur de feu a un m\u00eame LNAM et name_rcid<\/code> mais correspond \u00e0 une couleur est un secteur g\u00e9ographique diff\u00e9rents. De m\u00eame, le LNAM et 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.<\/li>\n\n\n\n aux tables SEAARE, LNDRGN et ADMARE qui couvrent des vastes zones et qui sont d\u00e9coup\u00e9es selon l’emprise de la carte, tout en gardant le m\u00eame LNAM. Dans ce cas m\u00eame si le LNAM est le m\u00eame, les polygones n’ont pas la m\u00eame forme et surface.<\/li>\n\n\n\nles tables DSID et C_AGGR qui n’ont pas de LNAM et qui, par prinipe, ne peuvent pas avoir des doublons.<\/li>\n<\/ul>\n\n\n\n<\/div>effacer doublons .py<\/div><\/div><\/div> <\/p>\n\n\n\nfrom osgeo import ogr<\/p>\n\n\n\n def delete_doublons(): # Ouvre le GeoPackage en mode mise \u00e0 jour geopackage_path = \u00ab\u00a0c:\/enc_db\/enc2.gpkg\u00a0\u00bb geopackage = ogr.Open(geopackage_path, update=True)<\/p>\n\n\n\nif geopackage is not None:\n # R\u00e9cup\u00e8re la liste de noms de toutes les tables du GeoPackage\n table_names = [geopackage.GetLayerByIndex(i).GetName() for i in range(geopackage.GetLayerCount())]\n\n # Boucle sur les tables et les traite individuellement\n for table_name in table_names:\n if table_name not in [\"pl_SEAARE\",\"pl_LNDRGN\",\"pl_ADMARE\",\"layer_styles\",\"natsurf\",\"DSID\",\"C_AGGR\"]:\n # Cr\u00e9e une table temporaire en ajoutant un suffixe \"_temp\"\n temp_table_name = table_name + \"_temp\"\n temp_table = geopackage.CreateLayer(temp_table_name)\n # Copie la structure de la table d'origine vers la table temporaire\n input_layer = geopackage.GetLayerByName(table_name)\n input_layer_defn = input_layer.GetLayerDefn()\n\n for i in range(input_layer_defn.GetFieldCount()):\n field_defn = input_layer_defn.GetFieldDefn(i)\n temp_table.CreateField(field_defn)\n if table_name not in [\"pt_LIGHTS\", \"pt_SOUNDG\"]:\n attribut = \"LNAM\"\n geopackage.ExecuteSQL(f\"INSERT INTO {temp_table_name} SELECT FROM {table_name} GROUP BY {attribut};\")\n elif table_name == \"pt_SOUNDG\":\n attribut = \"name_rcid,depth, ST_X(geom),ST_Y(geom)\"\n geopackage.ExecuteSQL(f\"INSERT INTO {temp_table_name} SELECT * FROM {table_name} GROUP BY {attribut};\")\n elif table_name == \"pt_LIGHTS\":\n attribut = \"name_rcid, colour, sectr1, sectr2\"\n geopackage.ExecuteSQL(f\"INSERT INTO {temp_table_name} SELECT * FROM {table_name} GROUP BY {attribut};\")\n geopackage.ExecuteSQL(f\"DROP TABLE {table_name}\")\n\n # Renomme la table temporaire pour qu'elle ait le m\u00eame nom que la table d'origine\n geopackage.ExecuteSQL(f\"ALTER TABLE {temp_table_name} RENAME TO {table_name}\")\n print(f\"alter table {table_name}\")\n\n # Ferme le GeoPackage\n geopackage = None\nelse:\n print(f\"Impossible d'ouvrir le GeoPackage {geopackage} en mode \u00e9criture.\")<\/code><\/pre>\n\n\n\n#Appel de la fonction delete_doublons<\/code><\/pre>\n\n\n\ndelete_doublons()<\/code><\/pre>\n\n\n\n <\/div><\/div>\n\n\n\n Veillez \u00e0 modifier la ligne<\/p>\n\n\n\n geopackage_path = \u00ab\u00a0c:\/testgpkg\/ENC.gpkg\u00a0\u00bb<\/em><\/strong><\/p>\n\n\n\n avec votre geopackage de base de donn\u00e9es.<\/p>\n\n\n\n Vous pouvez t\u00e9l\u00e9charger le fichier .py \u00e0 partir de ce lien : supprime_doublons.py.<\/a> 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\u00a0pl_SEAARE\u00a0\u00bb,\u00a0\u00bbpl_LNDRGN\u00a0\u00bb,\u00a0\u00bbpl_ADMARE\u00a0\u00bb,\u00a0\u00bblayer_styles\u00a0\u00bb,\u00a0\u00bbnatsurf\u00a0\u00bb,\u00a0\u00bbDSID\u00a0\u00bb,\u00a0\u00bbC_AGGR\u00a0\u00bb]:<\/em><\/strong><\/p>\n\n\n\n En principe, vous aurez que les tables pr\u00e9fix\u00e9es et les tables layer_styles et natsurf. Cette ligne permet de ne pas traiter pl_SEAARE, pl_LNDRGN et pl_ADMARE, mais surtout, elle emp\u00eache de vider les tables layer_styles<\/em><\/strong> et natsurf <\/em><\/strong>que nous utilisons pour les symbologies. 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 \u00e0 cette ligne d’exceptions.<\/p>\n\n\n\n <\/p>\n","protected":false},"excerpt":{"rendered":" La premi\u00e8re partie du projet Financement Collaboratif pour l’Int\u00e9gration de Donn\u00e9es Marines dans QGIS est arriv\u00e9e \u00e0 son terme gr\u00e2ce aux contributions d’ORANGE Marine, Geoconceptos Uruguay, Janez Avzec et Gaetan Mas. Nous les remercions vivement. Nous…<\/p>\n","protected":false},"author":1,"featured_media":12508,"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":[222,62],"tags":[1478,1634,131,85,223,192,180],"class_list":["post-12480","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-python","category-qgis-2","tag-geopackage","tag-ogr2ogr","tag-python","tag-s57","tag-script","tag-symbologie","tag-tutoriel"],"aioseo_notices":[],"jetpack_featured_media_url":"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2024\/03\/resultat-1.png?fit=1402%2C898&ssl=1","jetpack_shortlink":"https:\/\/wp.me\/p6XU0A-3fi","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/12480","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=12480"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/12480\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media\/12508"}],"wp:attachment":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media?parent=12480"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=12480"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=12480"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}