Le format S57 et les SIG<\/h2>\n\n\n\n
Si vous vous lancez \u00e0 travailler avec des fichiers S57 dans un SIG, il y a quelques \u00e9l\u00e9ments que vous devez conna\u00eetre pour naviguer sans \u00e9cueils.<\/p>\n\n\n\n
Tout d\u2019abord, la structure des fichiers S57 ne r\u00e9pond pas aux structures adopt\u00e9es dans les SIG.<\/p>\n\n\n\n
Dans un SIG vous avez un objet g\u00e9ographique qui est repr\u00e9sent\u00e9 par une table avec deux types d\u2019informations : la g\u00e9om\u00e9trie des entit\u00e9s de l\u2019objet et les attributs de ces entit\u00e9s.<\/p>\n\n\n\n
Si vous avez d\u2019autres objets qui ont des g\u00e9om\u00e9tries identiques, l\u2019information g\u00e9om\u00e9trique est dupliqu\u00e9e, une fois dans chaque table.<\/p>\n\n\n\n
Pour le format S57 l\u2019objectif principal est d\u2019optimiser le stockage des informations et par cons\u00e9quent de ne pas dupliquer les informations. Si un objet a une entit\u00e9 ponctuelle, un point sera cr\u00e9\u00e9. Du moment ou d\u2019autres objets ont des entit\u00e9s situ\u00e9es sur ce point, on utilisera la r\u00e9f\u00e9rence du point qui a d\u00e9j\u00e0 \u00e9t\u00e9 cr\u00e9\u00e9. De cette mani\u00e8re un point n\u2019est d\u00e9crit qu\u2019une seule fois dans le fichier. Il en va de m\u00eame pour les polylignes et les surfaces. Un fichier S57 aura donc une s\u00e9rie de tables contenant l\u2019information g\u00e9om\u00e9trique, d\u00e9nomm\u00e9es \u00ab primitives \u00bb :<\/p>\n\n\n\n
\n
IsolatedNode (points)<\/li>\n\n\n\n
ConnectedNode (points)<\/li>\n\n\n\n
Edge (polylignes)<\/li>\n\n\n\n
Face (polygones)<\/li>\n<\/ul>\n\n\n\n
La table attributaire des diff\u00e9rents objets S57 contient seulement les attributs des objets.<\/p>\n\n\n\n
Ce qui complique la t\u00e2che c\u2019est qu\u2019il y a deux attributs qui se r\u00e9f\u00e8rent aux g\u00e9om\u00e9tries : posacc <\/em><\/strong>(l\u2019estimation de la pr\u00e9cision de la position, valeur quantitative) et quapos <\/em><\/strong>(qualit\u00e9 de la position, variable qualitative).<\/p>\n\n\n\n
Ces deux attributs se trouvent dans les tables des primitives.<\/p>\n\n\n\n
*Pour passer de la structure S57 \u00e0 une structure SIG (**shapefile<\/a>, geopackage, postgis<\/a>) nous utilisons la biblioth\u00e8que GDAL et sa commande ogr2ogr.<\/p>\n\n\n\n*
Cette commande va cr\u00e9er des tables SIG \u00e0 partir de la structure S57, en cr\u00e9ant une table par objet S57, en affectant \u00e0 chaque entit\u00e9 les g\u00e9om\u00e9tries correspondantes tir\u00e9es des tables primitives et en ajoutant \u00e0 chaque entit\u00e9 les attributs de l\u2019objet S57. La trace des primitives utilis\u00e9es pour la g\u00e9om\u00e9trie de chaque entit\u00e9 se trouve dans le champ NAME_RCID <\/em><\/strong>des tables SIG, \u00e0 condition d\u2019avoir ajout\u00e9 les options -oo \u00ab\u00a0RETURN_LINKAGES=ON\u00a0\u00bb \u2013oo \u00ab\u00a0LNAM_REFS=ON\u00a0\u00bb<\/em><\/strong> \u00e0 la ligne de commande ogr2ogr.<\/p>\n\n\n\n
La figure suivante montre une couche de type point. La valeur indiqu\u00e9e dans le champ NAME_RCID <\/em><\/strong>est celle RCID <\/em><\/strong>du point utilis\u00e9 de la table IsolatedNode.<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
La figure suivante montre un exemple de couche de type lin\u00e9aire. Les valeurs indiqu\u00e9es dans le champ NAME_RCID sont celles des RCID <\/em><\/strong>des polylignes utilis\u00e9es de la table Edge.<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
La figure suivante montre un exemple de couche de type polygone. Les valeurs indiqu\u00e9es dans le champ NAME_RCID sont celles des polylignes utilis\u00e9es de la table Edge.<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
Pour pouvoir r\u00e9cup\u00e9rer les attributs QUAPOS <\/strong>et POSACC <\/strong>de chaque entit\u00e9 des tables de type point nous devons r\u00e9cup\u00e9rer les valeurs du point IsolatedNode <\/strong>pour les affecter aux tables des diff\u00e9rents objets ENC. <\/p>\n\n\n\n
Si les identifiants \u00e9taient directement les RCID dans les tables ENC, on pourrait faire une jointure<\/a> entre chaque table (soundg, obstrn,…) et IsolatedNode. Mais comme vous pouvez le voir dans les images pr\u00e9c\u00e9dentes, l’attribut NAME_RCID est de type stringlist <\/em>ce qui bloque cette solution. Nous avons donc d\u00e9velopp\u00e9 un script python qui fait ce travail lors du chargement des donn\u00e9es S57->geopackage.<\/p>\n\n\n\n
T\u00e9l\u00e9chargement des fichiers V3<\/h2>\n\n\n\nVous pouvez t\u00e9l\u00e9charger l’ensemble des fichiers .bat et .py cit\u00e9s dans cet article \u00e0 partir de ce lien.<\/a><\/p>\n\n\n\n
Modification des commandes ogr2ogr d\u2019import des fichiers S57<\/h2>\n\n\n\n
Import des couches \u00ab points \u00bb<\/h3>\n\n\n\n
Pour importer les primitives des fichiers S57 il est n\u00e9cessaire d\u2019ajouter l\u2019option –oo RETURN_PRIMITIVES = ON <\/em><\/strong> aux commandes d\u2019import. Dans le cas des points il est aussi n\u00e9cessaire d\u2019ajouter l\u2019option -nlt MULTIPOINT<\/em><\/strong> car si non la table IsolatedNode <\/em><\/strong>sera cr\u00e9\u00e9e comme de type point et les couches multipoint ne pourront pas \u00eatre prises en charge.<\/p>\n\n\n\n
Le fichier .bat d\u2019import des couches de type point sera alors :<\/p>\n\n\n\n
<\/div>
mass_load_s57_points_V3.bat<\/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 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 Compter le nombre total de fichiers \u00e0 traiter set \/a total_files=0 for \/r \"%directory%\" %%a in (.000) do ( set \/a total_files+=1 )<\/code><\/pre>\n\n\n\nREM Initialiser le compteur de fichiers trait\u00e9s\nset \/a processed_files=0\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 \/a processed_files+=1\necho Traitement en cours : !processed_files! sur !total_files!\nset \"file=%%~ni\"\nogr2ogr -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 -oo RETURN_PRIMITIVES=ON -nlt MULTIPOINT -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:\/testgpkgV3\/update_geopackage_dsid_V3.py \"%output_geopackage%\" \"!file!\"\n)\necho Import termin\u00e9. Ajout de posass et quapos.\npython c:\/testgpkgV3\/add_posacc_quapos_to_pointsV3.py \"%output_geopackage%\"\necho Traitement termin\u00e9.\npause<\/code><\/pre>\n\n\n\n<\/div><\/div>\n\n\n\n Pour l’ex\u00e9cuter, ouvrez un fen\u00eatre OSGeo4W Shell, allez dans le r\u00e9pertoire contenant les fichiers .bat et rentrez la commande:<\/p>\n\n\n\n .\\mass_load_s57_points_V3.bat r\u00e9pertoire_fichiers_000<\/em><\/strong> chemin\/pointsENC.gpkg<\/em><\/strong><\/p>\n\n\n\n Au pr\u00e9alable, \u00e9ditez le fichier pour modifier les lignes:<\/p>\n\n\n\n python c:\/testgpkgV3\/<\/strong>update_geopackage_dsid_prp_prim.py \u00ab\u00a0%output_geopackage%\u00a0\u00bb \u00ab\u00a0!file!\u00a0\u00bb<\/em><\/p>\n\n\n\n et<\/p>\n\n\n\n python c:\/testgpkgV3\/<\/strong>add_posacc_quapos_to_pointsV3.py<\/em><\/p>\n\n\n\n Pour rentrer les chemins que vous aurez choisis pour stocker les scripts python.<\/p>\n\n\n\n <\/a><\/figure>\n\n\n\nLe premier script, update_geopackage_dsid_V3.py<\/em><\/strong>, a \u00e9t\u00e9 modifi\u00e9 pour ajouter en plus de enc_chart, scale <\/em>et purpose<\/em>, les attributs POSACC <\/em>et QUAPOS <\/em>\u00e0 toutes les tables point, au moment du passage S57->geopackage.<\/p>\n\n\n\n Le deuxi\u00e8me script, add_posacc_quapos_to_pointsV3.py<\/em><\/strong>, une fois termin\u00e9 le chargement de tous les fichiers S57 du r\u00e9pertoire ENC, met \u00e0 jour les valeurs de ces attributs pour chaque entit\u00e9 des tables points.<\/p>\n\n\n\n<\/div>add_posacc_quapos_to_pointsV3.py<\/div><\/div><\/div><\/p>\n\n\n\nimport sys\nfrom osgeo import ogr\n\ndef add_posacc_quapos_to_pointsV3 (geopackage_path):\n # Chemin vers le GeoPackage\n\n # Liste des tables \u00e0 exclure\n tables_exclues = [\"IsolatedNode\", \"ConnectedNode\",\"DSID\",\"C_AGGR\",\"C_ASSO\",\"layer_styles\"] \n\n # Ouvrir le GeoPackage\n driver = ogr.GetDriverByName(\"GPKG\")\n geopackage = driver.Open(geopackage_path, 1)\n\n # R\u00e9cup\u00e9rer la table IsolatedNode\n isolated_node_table = geopackage.GetLayerByName(\"IsolatedNode\")\n\n # Parcourir toutes les tables du GeoPackage\n for i in range(geopackage.GetLayerCount()):\n table = geopackage.GetLayerByIndex(i)\n table_name = table.GetName()\n \n # V\u00e9rifier si la table n'est pas dans la liste des tables exclues et si elle n'est pas IsolatedNode\n if table_name not in tables_exclues and table_name != \"IsolatedNode\" :\n print(f\"Traitement table {table_name}\")\n # R\u00e9cup\u00e9rer les noms des champs des tables\n table_defn = table.GetLayerDefn()\n field_names = [table_defn.GetFieldDefn(j).GetName() for j in range(table_defn.GetFieldCount())]\n \n # Index des champs rcid et enc_chart dans la table IsolatedNode\n isolated_node_defn = isolated_node_table.GetLayerDefn()\n rcid_index_isolated_node = isolated_node_defn.GetFieldIndex(\"RCID\")\n enc_chart_index_isolated_node = isolated_node_defn.GetFieldIndex(\"enc_chart\")\n \n # Parcourir les enregistrements de la table\n table.ResetReading()\n for feature in table:\n rcid_full = feature.GetField(\"NAME_RCID\")\n rcid = rcid_full.split(\":\")[1]\n rcid = rcid[:-1]\n enc_chart = feature.GetField(\"enc_chart\")\n \n # Rechercher le rcid et enc_chart dans la table IsolatedNode\n isolated_node_table.SetAttributeFilter(\"RCID = '{}' AND enc_chart = '{}'\".format(rcid, enc_chart))\n isolated_node_feature = isolated_node_table.GetNextFeature()\n \n if isolated_node_feature:\n # Mettre \u00e0 jour les attributs POSACC et QUAPOS de la table en cours avec ceux d'IsolatedNode\n posacc_index = table_defn.GetFieldIndex(\"POSACC\")\n quapos_index = table_defn.GetFieldIndex(\"QUAPOS\")\n if posacc_index >= 0:\n posacc = isolated_node_feature.GetField(\"POSACC\")\n feature.SetField(posacc_index, posacc)\n if quapos_index >= 0:\n quapos = isolated_node_feature.GetField(\"QUAPOS\")\n feature.SetField(quapos_index, quapos)\n \n \n # Sauvegarder les modifications\n table.SetFeature(feature)\n\n # Fermer le GeoPackage\n geopackage = None\n\nif name == \"main\":\n if len(sys.argv) != 2:\n print(\"Usage: python script.py geopackage_path \")\n sys.exit(1)\n \n geopackage_path = sys.argv[1]\n \n\n add_posacc_quapos_to_pointsV3 (geopackage_path)\n\n <\/code><\/pre>\n\n\n\n<\/div><\/div>\n\n\n\n Dans l’exemple qui suit, nous avons charg\u00e9 7 cartes ENC dans le geopackage d’import pointsENC.gpkg.<\/p>\n\n\n\n<\/a><\/figure>\n\n\n\nVous pouvez constater que les tables du geopackage contiennent les attributs POSACC et QUAPOS et que le dernier script a bien renseign\u00e9 les valeurs correspondantes de ces champs.<\/p>\n\n\n\n Dans la deuxi\u00e8me partie de cette version finale consacr\u00e9e \u00e0 la symbologie vous verrez la nouvelle gestion de la couche \u00ab\u00a0nature des fonds\u00a0\u00bb. Cette symbologie est bas\u00e9e sur la cr\u00e9ation d’un atribut \u00ab\u00a0Label\u00a0\u00bb dans la couche SBDARE et son remplissage avec un script python Label_sbdare.py<\/em><\/strong> . Cette op\u00e9ration peut se faire sur le base de donn\u00e9es compl\u00e8te ou ici, sur le geopackage pointsENC.gpkg. Comme elle prend un certain temps de calcul, nous vous conseillons de le faire \u00e0 ce stade. R\u00e9f\u00e9rez-vous \u00e0 l’article Cartes ENC Geopackage dans QGis Version finale: deuxi\u00e8me partie<\/em><\/strong> pour le mode op\u00e9ratoire.<\/p>\n\n\n\nImport des couches \u00ab lignes \u00bb<\/h3>\n\n\n\nPour importer les primitives des fichiers S57 il est n\u00e9cessaire d\u2019ajouter l\u2019option –oo RETURN_PRIMITIVES = ON<\/em><\/strong> aux commandes d\u2019import. <\/p>\n\n\n\n Le fichier .bat d\u2019import des couches de type ligne sera alors :<\/p>\n\n\n\n<\/div>mass_load_s57_lignes_V3.bat<\/div><\/div><\/div><\/p>\n\n\n\n@echo off\nsetlocal enabledelayedexpansion\n\nREM Assurez-vous qu'au moins 2 arguments ont \u00e9t\u00e9 fournis\nif \"%~2\"==\"\" (\n echo Usage: %0 directory000 output_geopackage \n exit \/b 1\n)\n\nREM R\u00e9cup\u00e9rer les arguments de la ligne de commande\nset \"directory=%~1\"\nset \"output_geopackage=%~2\"\n\nREM Compter le nombre total de fichiers \u00e0 traiter\nset \/a total_files=0\nfor \/r \"%directory%\" %%a in (.000) do (\n set \/a total_files+=1\n)\n\nREM Initialiser le compteur de fichiers trait\u00e9s\nset \/a processed_files=0\nREM It\u00e9rez sur tous les fichiers .000 dans le r\u00e9pertoire\nfor \/r \"%directory%\" %%i in (.000) do (\n echo Traitement du fichier: %%~ni\n\tset \/a processed_files+=1\n echo Traitement en cours : !processed_files! sur !total_files!\n\tset \"file=%%~ni\"\n ogr2ogr -f GPKG -skipfailures -append -update -where \"OGR_GEOMETRY='LINESTRING' or OGR_GEOMETRY='MULTILINESTRING'\" -oo RETURN_PRIMITIVES=ON -oo RETURN_LINKAGES=ON -oo LNAM_REFS=ON -mapFieldType StringList=String,IntegerList=String \"%output_geopackage%\" \"%%i\"\n ogr2ogr -f GPKG -skipfailures -append -update \"%output_geopackage%\" \"%%i\" \"DSID\"\n\togr2ogr -f GPKG -skipfailures -append -update \"%output_geopackage%\" \"%%i\" \"C_AGGR\"\n\tpython c:\/testgpkgV3\/update_geopackage_dsid_V3.py \"%output_geopackage%\" \"!file!\" \n)\necho Import termin\u00e9. Ajout de posacc et quapos.\npython c:\/testgpkgV3\/add_posacc_quapos_to_linesV3.py \"%output_geopackage%\"\necho Traitement termin\u00e9.\npause\n<\/code><\/pre>\n\n\n\n<\/div><\/div>\n\n\n\n Pour l’ex\u00e9cuter, ouvrez un fen\u00eatre OSGeo4W Shell, allez dans le r\u00e9pertoire contenant les fichiers .bat et rentrez la commande:<\/p>\n\n\n\n .\\mass_load_s57_lines_V3.bat r\u00e9pertoire_fichiers_000<\/em><\/strong> chemin\/pointsENC.gpkg<\/em><\/strong><\/p>\n\n\n\n Au pr\u00e9alable, \u00e9ditez le fichier pour modifier les lignes:<\/p>\n\n\n\n python c:\/testgpkgV3\/<\/strong>update_geopackage_dsid_prp_prim.py \u00ab\u00a0%output_geopackage%\u00a0\u00bb \u00ab\u00a0!file!\u00a0\u00bb<\/em><\/p>\n\n\n\n et<\/p>\n\n\n\n python c:\/testgpkgV3\/<\/strong>add_posacc_quapos_to_linesV3.py<\/em><\/p>\n\n\n\n pour rentrer les chemins que vous aurez choisis pour stocker les scripts python.<\/p>\n\n\n\n Le deuxi\u00e8me script, add_posacc_quapos_to_linesV3.py<\/em><\/strong>, une fois termin\u00e9 le chargement de tous les fichiers S57 du r\u00e9pertoire ENC, met \u00e0 jour les valeurs de ces attributs pour chaque entit\u00e9 des tables lignes.<\/p>\n\n\n\n <\/div>add_posacc_quapos_to_linesV3.py<\/div><\/div><\/div><\/p>\n\n\n\nimport sys\nfrom osgeo import ogr\n\ndef add_posacc_quapos_to_linesV3 (geopackage_path):\n\n # Liste des tables \u00e0 exclure\n tables_exclues = [\"Edge\", \"ConnectedNode\",\"DSID\",\"C_AGGR\",\"C_ASSO\",\"Edge\",\"FACE\"] # Ajoutez les noms des tables que vous souhaitez exclure\n\n # Ouvrir le GeoPackage\n driver = ogr.GetDriverByName(\"GPKG\")\n geopackage = driver.Open(geopackage_path, 1)\n\n # R\u00e9cup\u00e9rer la table Edge\n isolated_node_table = geopackage.GetLayerByName(\"Edge\")\n\n # Parcourir toutes les tables du GeoPackage\n for i in range(geopackage.GetLayerCount()):\n table = geopackage.GetLayerByIndex(i)\n table_name = table.GetName()\n \n # V\u00e9rifier si la table n'est pas dans la liste des tables exclues et si elle n'est pas Edge\n if table_name not in tables_exclues and table_name != \"Edge\":\n print(f\"Traitement table {table.name()}\")\n # R\u00e9cup\u00e9rer les noms des champs des tables\n table_defn = table.GetLayerDefn()\n field_names = [table_defn.GetFieldDefn(j).GetName() for j in range(table_defn.GetFieldCount())]\n \n # Index des champs rcid et enc_chart dans la table Edge\n isolated_node_defn = isolated_node_table.GetLayerDefn()\n rcid_index_isolated_node = isolated_node_defn.GetFieldIndex(\"RCID\")\n enc_chart_index_isolated_node = isolated_node_defn.GetFieldIndex(\"enc_chart\")\n \n # Parcourir les enregistrements de la table\n table.ResetReading()\n for feature in table:\n rcid_full = feature.GetField(\"NAME_RCID\")\n rcid = rcid_full.split(\":\")[1]\n rcid = rcid[:-1]\n enc_chart = feature.GetField(\"enc_chart\")\n \n # Rechercher le rcid et enc_chart dans la table Edge\n isolated_node_table.SetAttributeFilter(\"RCID = '{}' AND enc_chart = '{}'\".format(rcid, enc_chart))\n isolated_node_feature = isolated_node_table.GetNextFeature()\n \n if isolated_node_feature:\n # Mettre \u00e0 jour les attributs POSACC et QUAPOS de la table en cours avec ceux d'Edge\n posacc_index = table_defn.GetFieldIndex(\"POSACC\")\n quapos_index = table_defn.GetFieldIndex(\"QUAPOS\")\n if posacc_index >= 0:\n posacc = isolated_node_feature.GetField(\"POSACC\")\n feature.SetField(posacc_index, posacc)\n if quapos_index >= 0:\n quapos = isolated_node_feature.GetField(\"QUAPOS\")\n feature.SetField(quapos_index, quapos)\n \n \n # Sauvegarder les modifications\n table.SetFeature(feature)\n\n # Fermer le GeoPackage\n geopackage = None\n \n\nif name == \"main\":\n if len(sys.argv) != 2:\n print(\"Usage: python script.py geopackage_path \")\n sys.exit(1)\n \n geopackage_path = sys.argv[1]\n \n\n add_posacc_quapos_to_linesV3 (geopackage_path)\n<\/code><\/pre>\n\n\n\n<\/div><\/div>\n\n\n\n Dans le cas des tables de type point, le lien entre le NAME_RCID <\/em><\/strong>des tables et le RCID <\/em><\/strong>de la table IsolatedNode <\/em><\/strong>ne pose pas de question. Par contre dans le cas des polylignes, pour une entit\u00e9 g\u00e9ographique, la liste des polylignes concerne plusieurs entit\u00e9s de la couche Edge<\/em><\/strong>. Nous nous sommes pos\u00e9 la question de savoir si pour une entit\u00e9 , par exemple de DEPCNT, il \u00e9tait possible que les polylignes qui la composent aient diff\u00e9rentes valeurs de QUAPOS. <\/p>\n\n\n\n Nous avons effectu\u00e9 le test sur une base contenant 650 cartes ENC \u00e0 toutes les \u00e9chelles. La couche DEPCNT contenait 162 585 entit\u00e9s. Seules 38 731 avaient dans la liste des NAME_RCID des RCID de la table Edges avec une valeur non nulle de QUAPOS. Nous avons trouv\u00e9 193 entit\u00e9s de DEPCNT qui avaient des valeurs diff\u00e9rentes de QUAPOS dans leur liste de NAME_RCID. Etant donn\u00e9 que cela repr\u00e9sente 0,1% nous avons d\u00e9cid\u00e9 de ne pas prendre en compte cette possibilit\u00e9, car la complexit\u00e9 de d\u00e9construction des polylignes est extr\u00eame. Il s’av\u00e8re que la table Edge r\u00e9sultante du driver ogr2ogr ne pose pas de probl\u00e8me au niveau des attributs (la jointure entre NAME_RCID et RCID est bonne) mais que les g\u00e9om\u00e9tries de la table comportent presque la moiti\u00e9 de g\u00e9om\u00e9tries non valides (des lignes avec un seul point ou des g\u00e9om\u00e9tries nulles). Pourtant les polylignes des entit\u00e9s S57 n’ont aucune erreur.<\/p>\n\n\n\n Il faudrait donc, r\u00e9parer et reconstruire toutes les g\u00e9om\u00e9tries de Edge pour pouvoir d\u00e9construire les g\u00e9om\u00e9tries de DEPCNT pour ne pas avoir de QUAPOS diff\u00e9rents dans aucun cas.<\/p>\n\n\n\n Dans l’exemple suivant nous avons charg\u00e9 7 cartes ENC dans le geopackage d’import linesENC.gpkg<\/p>\n\n\n\n <\/a><\/figure>\n\n\n\nVous pouvez constater que les tables du geopackage contiennent les attributs POSACC et QUAPOS et que le dernier script a bien renseign\u00e9 les valeurs correspondantes de ces champs.<\/p>\n\n\n\n Import des couches \u00ab polygones \u00bb<\/h2>\n\n\n\nM\u00eame si le format S57 pr\u00e9voit la possibilit\u00e9 d’avoir des polygones en tant que primitives, cette option n’est pas utilis\u00e9e. Il n’est donc pas n\u00e9cessaire de modifier l’import des donn\u00e9es de type polygone. Les primitives (Face) sont demand\u00e9es dans la ligne de commande ogr2ogr, les champs POSACC et QUAPOS sont bien ajout\u00e9s aux tables, pour un \u00e9ventuelle utilisation ult\u00e9rieure, mais il n’est pas n\u00e9cessaire d’ajouter un script pour remplir les valeurs.<\/p>\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 Par rapport \u00e0 la V2 les modifications apport\u00e9es au script sont les suivantes:<\/p>\n\n\n\n\nRenseignement du SRC<\/a> des tables du geopackage lors de la cr\u00e9ation. Le code EPSG:4326<\/a> est affect\u00e9 aux g\u00e9om\u00e9tries.<\/li>\n<\/ul>\n\n\n\n# Cr\u00e9er une nouvelle couche de destination si elle n'existe pas\n output_srs = osr.SpatialReference()\n output_srs.ImportFromEPSG(4326) # D\u00e9finir EPSG 4326\n output_table = output_gpkg.CreateLayer(output_table_name, geom_type=input_table.GetGeomType(), options=[\"OVERWRITE=YES\"], srs=output_srs)<\/em>\n<\/code><\/pre>\n\n\n\n\nFor\u00e7age de la cr\u00e9ation de l’index spatial sur toutes les tables du geopackage<\/li>\n<\/ul>\n\n\n\n# Cr\u00e9er l'index spatial\n output_gpkg.ExecuteSQL(f\"CREATE SPATIAL INDEX ON {output_table_name}\")<\/em><\/code><\/pre>\n\n\n\nLe clonage est r\u00e9alis\u00e9 avec le code python : clone_or_append_tables_with_prefixV3.py<\/a> Veillez \u00e0 modifier la ligne qui indique les noms des geopackages que vous avez cr\u00e9\u00e9s avec les fichiers .bat pr\u00e9c\u00e9dents :<\/p>\n\n\n\n input_geopackages = [\u00ab c:\/testgpkg\/pointsENC.gpkg \u00bb, \u00ab c:\/testgpkg\/LinesENC.gpkg \u00bb, \u00ab c:\/testgpkg\/PolysENC.gpkg \u00bb]<\/em><\/strong><\/p>\n\n\n\n avec les geopackages d\u2019import que vous avez choisi, ainsi que la ligne<\/p>\n\n\n\n output_gpkg = ogr.Open(\u00ab c:\/testgpkg\/ENC.gpkg \u00bb, 1)<\/em><\/strong><\/p>\n\n\n\n avec le geopackage de base de donn\u00e9es que vous aurez cr\u00e9\u00e9 en suivant les indications de l’article Cr\u00e9ation d\u2019une base de donn\u00e9es Geopackage pour les cartes ENC (1\u00e8re partie: constituer la base de donn\u00e9es)<\/a>.<\/p>\n\n\n\n Une fois ce script ex\u00e9cut\u00e9, vous avez votre base de donn\u00e9es \u00e0 jour des nouveaux fichiers S57.<\/p>\n\n\n\n Mise \u00e0 jour d’une base geopackage V1 ou V2 vers la version 3<\/h2>\n\n\n\nLa proc\u00e9dure que nous venons de d\u00e9crire permet de cr\u00e9er une nouvelle base de donn\u00e9es geopackage. Mais si vous avez d\u00e9j\u00e0 une base qui a \u00e9t\u00e9 cr\u00e9\u00e9e il n’est pas n\u00e9cessaire de refaire tout le travail. Nous vous donnons ici la proc\u00e9dure et les outils n\u00e9cessaires pour mettre \u00e0 jour sa structure (POSAC et QUAPOSQ) et de renseigner les valeurs de ces attributs pour les tables d\u00e9j\u00e0 existantes.<\/p>\n\n\n\nT\u00e9l\u00e9chargez les fichiers pour la proc\u00e9dure ici.<\/a><\/p>\n\n\n\n Les diff\u00e9rentes \u00e9tapes de cet \u00ab\u00a0upgrade\u00a0\u00bb sont:<\/p>\n\n\n\n 1- A partir d’un r\u00e9pertoire contenant tous les fichiers .000 qui ont d\u00e9j\u00e0 \u00e9t\u00e9 charg\u00e9s dans la base, ex\u00e9cution d’un fichier.bat , upgrade_V2_to_V3.ba<\/em><\/strong>t, qui cr\u00e9e un geopackage update.gpkg<\/em><\/strong> contenant les primitives (IsolatedNode, ConnectedNode,Edge) des fichiers S57.<\/p>\n\n\n\n .\\ upgrade_V2_to_V3.ba r\u00e9pertoire_fichiers_000<\/em><\/strong> chemin\/update.gpkg<\/em><\/strong><\/p>\n\n\n\n 2- Dans la console Python de QGis, ex\u00e9cution du script append_tables_upgrade.py<\/em><\/strong>. Veillez \u00e0 modifier le chemin pour le fichier update.gpkg<\/p>\n\n\n\n input_geopackages = [\u00ab\u00a0c:\/testgpkgV2\/update.gpkg<\/em><\/strong>\u00ab\u00a0]<\/p>\n\n\n\n et pour votre base de donn\u00e9es ENC<\/p>\n\n\n\n output_gpkg = ogr.Open(\u00ab\u00a0c:\/enc_db\/ENC4.gpkg<\/strong><\/em>\u00ab\u00a0, 1)<\/p>\n\n\n\n Ce script ajoute les primitives \u00e0 la base de donn\u00e9es existante.<\/p>\n\n\n\n 3- Dans la console Python de QGis, ex\u00e9cution du script add_fields_upgrade.py<\/em><\/strong>. Veillez \u00e0 modifier le chemin pour le fichier .gpkg de votre base de donn\u00e9es<\/p>\n\n\n\n geopackage = ogr.Open(\u00ab\u00a0c:\/enc_db\/ENC4.gpkg<\/em><\/strong>\u00ab\u00a0, 1)<\/p>\n\n\n\n Ce script cr\u00e9e les attributs POSACC et QUAPOS dans toutes les tables de la base de donn\u00e9es.<\/p>\n\n\n\n 4-Dans la console Python de QGis, ex\u00e9cution du script add_posacc_quapos_upgrade_points.py<\/em><\/strong>. Veillez \u00e0 modifier le chemin pour le fichier .gpkg de votre base de donn\u00e9es et \u00e0 ajouter des tables \u00e0 exclure si vous en avez cr\u00e9\u00e9 d’autres tables dans la base de donn\u00e9es.<\/p>\n\n\n\n #Chemin vers le GeoPackage geopackage_path = \u00ab\u00a0c:\/enc_db\/ENC4.gpkg\u00a0\u00bb #Liste des tables \u00e0 exclure tables_exclues = [\u00ab\u00a0IsolatedNode\u00a0\u00bb, \u00ab\u00a0ConnectedNode\u00a0\u00bb,\u00a0\u00bbDSID\u00a0\u00bb,\u00a0\u00bbC_AGGR\u00a0\u00bb,\u00a0\u00bbC_ASSO\u00a0\u00bb,\u00a0\u00bblayer_styles\u00a0\u00bb] # Ajoutez les noms des tables que vous souhaitez exclure<\/em><\/strong><\/p>\n\n\n\n Ce script va renseigner les champs POSACC et QUAPOS dans toutes les tables de type Point avec les valeurs contenues dans la primitive IsolatedNode.<\/p>\n\n\n\n 5-Dans la console Python de QGis, ex\u00e9cution du script add_posacc_quapos_upgrade_lines.py<\/em><\/strong>. Veillez \u00e0 modifier le chemin pour le fichier .gpkg de votre base de donn\u00e9es et \u00e0 ajouter des tables \u00e0 exclure si vous en avez cr\u00e9\u00e9 d’autres tables dans la base de donn\u00e9es.<\/p>\n\n\n\n #Chemin vers le GeoPackage geopackage_path = \u00ab\u00a0c:\/enc_db\/ENC4.gpkg\u00a0\u00bb #Liste des tables \u00e0 exclure tables_exclues = [\u00ab\u00a0IsolatedNode\u00a0\u00bb, \u00ab\u00a0ConnectedNode\u00a0\u00bb,\u00a0\u00bbDSID\u00a0\u00bb,\u00a0\u00bbC_AGGR\u00a0\u00bb,\u00a0\u00bbC_ASSO\u00a0\u00bb,\u00a0\u00bblayer_styles\u00a0\u00bb] # Ajoutez les noms des tables que vous souhaitez exclure<\/em><\/strong><\/p>\n\n\n\n Ce script va renseigner les champs POSACC et QUAPOS dans toutes les tables de type Ligne avec les valeurs contenues dans la primitive Edge.<\/p>\n\n\n\n <\/p>\n","protected":false},"excerpt":{"rendered":" Le projet d’int\u00e9grations des cartes ENC dans QGis sous la forme de Geopackage est maintenant termin\u00e9e. Cet article compl\u00e8te les deux articles d\u00e9j\u00e0 publi\u00e9s concernant la mise en place et la gestion de la base de…<\/p>\n","protected":false},"author":1,"featured_media":12765,"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_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":"","jetpack_post_was_ever_published":false},"categories":[62],"tags":[58,85,41,192],"class_list":["post-12754","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-qgis-2","tag-qgis","tag-s57","tag-svg","tag-symbologie"],"aioseo_notices":[],"jetpack_featured_media_url":"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2024\/04\/aert5.png?fit=3734%2C1957&ssl=1","jetpack_shortlink":"https:\/\/wp.me\/p6XU0A-3jI","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/12754","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=12754"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/12754\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media\/12765"}],"wp:attachment":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media?parent=12754"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=12754"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=12754"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}