El formato S57 y el SIG<\/h2>\n\n\n\n
Si va a empezar a trabajar con ficheros S57 en un SIG, hay algunas cosas que debe tener en cuenta para navegar sin problemas.<\/p>\n\n\n\n
En primer lugar, la estructura de los ficheros S57 no se corresponde con las estructuras adoptadas en los SIG.<\/p>\n\n\n\n
En un SIG tienes un objeto geogr\u00e1fico que est\u00e1 representado por una tabla con dos tipos de informaci\u00f3n: la geometr\u00eda de las entidades del objeto y los atributos de estas entidades.<\/p>\n\n\n\n
Si tiene otros objetos con geometr\u00edas id\u00e9nticas, la informaci\u00f3n geom\u00e9trica se duplica, una vez en cada tabla.<\/p>\n\n\n\n
En el formato S57, el objetivo principal es optimizar el almacenamiento de la informaci\u00f3n y, por tanto, no duplicarla. Si un objeto tiene una entidad puntual, se crear\u00e1 un punto. Si otros objetos tienen entidades situadas en este punto, se utilizar\u00e1 la referencia del punto ya creado. De este modo, un punto s\u00f3lo se describe una vez en el fichero. Lo mismo ocurre con las polil\u00edneas y las superficies. Por tanto, un fichero S57 tendr\u00e1 una serie de tablas que contienen informaci\u00f3n geom\u00e9trica, conocidas como \u00ab\u00a0primitivas\u00a0\u00bb:<\/p>\n\n\n\n
\n
IsolatedNode (puntos)<\/li>\n\n\n\n
ConnectedNode (puntos)<\/li>\n\n\n\n
Edge (polil\u00edneas)<\/li>\n\n\n\n
Face (pol\u00edgonos)<\/li>\n<\/ul>\n\n\n\n
La tabla de atributos de los distintos objetos S57 s\u00f3lo contiene los atributos de los objetos.<\/p>\n\n\n\n
Lo que complica la tarea es que hay dos atributos que se refieren a las geometr\u00edas: posacc <\/em><\/strong>(la precisi\u00f3n estimada de la posici\u00f3n, un valor cuantitativo) y quapos <\/em><\/strong>(la calidad de la posici\u00f3n, una variable cualitativa).<\/p>\n\n\n\n
Estos dos atributos se encuentran en las tablas primitivas.<\/p>\n\n\n\n
*Para pasar de la estructura S57 a una estructura SIG (**shapefile<\/a>, geopackage, postgis) utilizamos la biblioteca GDAL y su comando ogr2ogr.<\/p>\n\n\n\n*
Este comando crear\u00e1 tablas SIG a partir de la estructura S57, creando una tabla por cada objeto S57, asignando la geometr\u00eda correspondiente de las tablas de primitivas a cada entidad y a\u00f1adiendo los atributos del objeto S57 a cada entidad. La traza de las primitivas utilizadas para la geometr\u00eda de cada entidad puede encontrarse en el campo NAME_RCID <\/em><\/strong>de las tablas SIG, siempre que se hayan a\u00f1adido las opciones -oo \u00ab\u00a0RETURN_LINKAGES=ON\u00a0\u00bb -oo \u00ab\u00a0LNAM_REFS=ON\u00a0\u00bb<\/em><\/strong> a la l\u00ednea de comandos de ogr2ogr.<\/p>\n\n\n\n
La siguiente figura muestra una capa de tipo punto. El valor indicado en el campo NAME_RCID<\/em><\/strong> es el RCID <\/em><\/strong>del punto utilizado en la tabla IsolatedNode.<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
La siguiente figura muestra un ejemplo de capa de tipo lineal. Los valores indicados en el campo NAME_RCID son los de los RCID de las polil\u00edneas utilizadas en la tabla Edge.<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
La siguiente figura muestra un ejemplo de capa de tipo poligonal. Los valores indicados en el campo NAME_RCID son los de las polil\u00edneas utilizadas en la tabla Edge.<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
Para recuperar los atributos QUAPOS <\/em><\/strong>y POSACC <\/strong><\/em>de cada entidad en las tablas de tipo punto, necesitamos recuperar los valores del punto IsolatedNode y asignarlos a las tablas de los distintos objetos ENC.<\/p>\n\n\n\n
Si los identificadores fueran directamente los RCID de las tablas ENC, podr\u00edamos hacer un join entre cada tabla (soundg, obstrn,\u2026) y IsolatedNode. Pero como se puede ver en las im\u00e1genes anteriores, el atributo NAME_RCID <\/em><\/strong>es de tipo stringlist, lo que bloquea esta soluci\u00f3n. Por ello hemos desarrollado un script en python que hace el trabajo cuando se cargan los datos de S57->geopackage.<\/p>\n\n\n\n
Descarga de archivos V3<\/h2>\n\n\n\nPuede descargar todos los archivos .bat y .py mencionados en este art\u00edculo desde este enlace.<\/a><\/p>\n\n\n\n
Modificaci\u00f3n de los comandos de ogr2ogr para importar archivos S57<\/h2>\n\n\n\n
Importaci\u00f3n de capas \u00ab\u00a0puntuales<\/h3>\n\n\n\n
Para importar primitivas desde ficheros S57, es necesario a\u00f1adir la opci\u00f3n -oo RETURN_PRIMITIVES = ON<\/em><\/strong> a los comandos de importaci\u00f3n. En el caso de los puntos, tambi\u00e9n es necesario a\u00f1adir la opci\u00f3n -nlt MULTIPOINT<\/em><\/strong>, ya que de lo contrario la tabla IsolatedNode se crear\u00e1 como un punto y las capas multipunto no ser\u00e1n compatibles.<\/p>\n\n\n\n
El archivo .bat para importar capas de tipo punto ser\u00e1 entonces :<\/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 Para ejecutarlo, abra una ventana OSGeo4W Shell, vaya al directorio que contiene los archivos .bat e introduzca el comando:<\/p>\n\n\n\n .\\mass_load_s57_points_V3.bat directorio_archivos_000 ruta\/pointsENC.gpkg<\/em><\/strong><\/p>\n\n\n\n En primer lugar, editar el archivo para cambiar las l\u00edneas:<\/p>\n\n\n\n python c:\/testgpkgV3\/update_geopackage_dsid_prp_prim.py \u00ab\u00a0%output_geopackage%\u00a0\u00bb \u00ab\u00a0!file!\u00a0\u00bb<\/em><\/strong><\/p>\n\n\n\n y<\/p>\n\n\n\n python c:\/testgpkgV3\/add_posacc_quapos_to_pointsV3.py<\/em><\/strong><\/p>\n\n\n\n Para introducir las rutas que ha elegido para almacenar los scripts de python.<\/p>\n\n\n\n <\/a><\/figure>\n\n\n\nEl primer script, update_geopackage_dsid_V3.py,<\/em><\/strong> se ha modificado para a\u00f1adir los atributos POSACC y QUAPOS a todas las tablas de puntos, adem\u00e1s de enc_chart, scale y purpose, cuando se realiza el cambio S57->geopackage.<\/p>\n\n\n\n El segundo script, add_posacc_quapos_to_pointsV3.py<\/em><\/strong>, una vez cargados todos los ficheros S57 en el directorio ENC, actualiza los valores de estos atributos para cada entidad en las tablas de puntos.<\/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 En el siguiente ejemplo, hemos cargado 7 mapas ENC en el geopackage import pointsENC.gpkg.<\/p>\n\n\n\n<\/a><\/figure>\n\n\n\nPuedes ver que las tablas del geopackage contienen los atributos POSACC y QUAPOS y que el \u00faltimo script ha rellenado los valores correspondientes para estos campos.<\/p>\n\n\n\n En la segunda parte de esta versi\u00f3n final, dedicada a la simbolog\u00eda, ver\u00e1s la nueva gesti\u00f3n de la capa \u00ab\u00a0naturaleza del fondo marino\u00a0\u00bb. Esta simbolog\u00eda se basa en la creaci\u00f3n de un atributo \u00ab\u00a0Label\u00a0\u00bb en la capa SBDARE y su relleno con un script python Label_sbdare.py<\/em><\/strong> . Esta operaci\u00f3n puede realizarse sobre la base de datos completa o aqu\u00ed, sobre el geopackage pointsENC.gpkg. Como se tarda alg\u00fan tiempo en calcular, le aconsejamos que lo haga en esta fase. Consulte el art\u00edculo Mapas del geopackage ENC en QGis : versi\u00f3n final parte 2<\/em><\/strong> para obtener instrucciones.<\/p>\n\n\n\nImportar capas de ‘l\u00ednea<\/h2>\n\n\n\nPara importar primitivas desde archivos S57, necesita a\u00f1adir la opci\u00f3n -oo RETURN_PRIMITIVES = ON<\/em><\/strong> a los comandos de importaci\u00f3n.<\/p>\n\n\n\n El archivo .bat de importaci\u00f3n de capas de l\u00ednea ser\u00e1 entonces :<\/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 Para ejecutarlo, abra una ventana OSGeo4W Shell, vaya al directorio que contiene los archivos .bat e introduzca el comando:<\/p>\n\n\n\n .\\mass_load_s57_lines_V3.bat directorio_archivos_000 ruta\/puntosENC.gpkg<\/em><\/strong><\/p>\n\n\n\n En primer lugar, editar el archivo para cambiar las l\u00edneas:<\/p>\n\n\n\n python c:\/testgpkgV3<\/em><\/strong>\/update_geopackage_dsid_prp_prim.py \u00ab\u00a0%output_geopackage%\u00a0\u00bb \u00ab\u00a0!file!\u00a0\u00bb<\/p>\n\n\n\n y<\/p>\n\n\n\n python c:\/testgpkgV3<\/em><\/strong>\/add_posacc_quapos_to_linesV3.py<\/p>\n\n\n\n para introducir las rutas que ha elegido para almacenar los scripts python.<\/p>\n\n\n\n El segundo script, add_posacc_quapos_to_linesV3.py<\/em><\/strong>, una vez cargados todos los ficheros S57 en el directorio ENC, actualiza los valores de estos atributos para cada entidad en las tablas de l\u00edneas.<\/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 En el caso de las tablas de tipo punto, el v\u00ednculo entre el NAME_RCID de las tablas y el RCID de la tabla IsolatedNode no plantea ninguna duda. Sin embargo, en el caso de las polil\u00edneas, para una entidad geogr\u00e1fica, la lista de polil\u00edneas concierne a varias entidades de la capa Edge. Nos preguntamos si para una caracter\u00edstica, por ejemplo una caracter\u00edstica DEPCNT, era posible que las polil\u00edneas que la componen tuvieran valores de QUAPOS diferentes.<\/p>\n\n\n\n Realizamos la prueba en una base de datos que conten\u00eda 650 mapas ENC a todas las escalas. La capa DEPCNT conten\u00eda 162.585 elementos. S\u00f3lo 38.731 ten\u00edan RCID en la lista NAME_RCID de la tabla Edges con un valor QUAPOS distinto de cero. Encontramos 193 entidades DEPCNT que ten\u00edan diferentes valores QUAPOS en su lista NAME_RCID. Dado que esto representa el 0,1%, decidimos no tener en cuenta esta posibilidad, ya que la complejidad de la deconstrucci\u00f3n de polil\u00edneas es extrema. Resulta que la tabla Edge resultante del controlador ogr2ogr no plantea ning\u00fan problema a nivel de atributos (la uni\u00f3n entre NAME_RCID y RCID es buena), pero que las geometr\u00edas de la tabla contienen casi la mitad de geometr\u00edas no v\u00e1lidas (l\u00edneas con un solo punto o geometr\u00edas nulas). Sin embargo, las polil\u00edneas de las entidades S57 no presentan errores.<\/p>\n\n\n\n Por tanto, ser\u00eda necesario reparar y reconstruir todas las geometr\u00edas Edge para poder deconstruir las geometr\u00edas DEPCNT y no tener QUAPOS diferentes en cada caso.<\/p>\n\n\n\n En el siguiente ejemplo, hemos cargado 7 mapas ENC en el geopackage de importaci\u00f3n linesENC.gpkg<\/p>\n\n\n\n <\/a><\/figure>\n\n\n\nImportaci\u00f3n de capas poligonales<\/h2>\n\n\n\nAunque el formato S57 permite utilizar pol\u00edgonos como primitivas, esta opci\u00f3n no se utiliza. Por lo tanto, no es necesario modificar la importaci\u00f3n de datos poligonales. Las primitivas (Face) se solicitan en la l\u00ednea de comandos ogr2ogr, y los campos POSACC y QUAPOS se a\u00f1aden a las tablas para su posible uso futuro, pero no es necesario a\u00f1adir un script para rellenar los valores.<\/p>\n\n\n\nClonaci\u00f3n y\/o adici\u00f3n de tablas importadas al geopackage ENC<\/h2>\n\n\n\nPara la primera carga del S57, el siguiente script de Python crea todas las tablas en los tres geopacks de importaci\u00f3n. Para las cargas posteriores, si la tabla ya existe en el geopackage ENC, los registros de la tabla en el geopackage de importaci\u00f3n se a\u00f1aden a los registros ya presentes en ENC. Si, por el contrario, la tabla a\u00fan no existe, se crea y se rellena con los registros del geopackage de importaci\u00f3n. En el archivo del geopackage ENC, las tablas llevan los prefijos pt_, li_ y pl_ para indicar su tipo de geometr\u00eda.<\/p>\n\n\n\n Se han introducido los siguientes cambios en el script en comparaci\u00f3n con la V2:<\/p>\n\n\n\n\nEl SRC<\/a> de las tablas del geopackage se rellena durante la creaci\u00f3n. Se asigna el c\u00f3digo EPSG:4326<\/a> a las geometr\u00edas.<\/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\nForzar la creaci\u00f3n del \u00edndice espacial en todas las tablas de paquetes geogr\u00e1ficos<\/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\nLa clonaci\u00f3n se realiza mediante el c\u00f3digo python: clone_or_append_tables_with_prefixV3.py<\/a> Aseg\u00farate de cambiar la l\u00ednea que indica los nombres de los geopackages que creaste con los archivos .bat anteriores:<\/p>\n\n\n\n input_geopackages = [\u00a0\u00bb c:\/testgpkg\/pointsENC.gpkg \u00ab\u00a0, \u00a0\u00bb c:\/testgpkg\/LinesENC.gpkg \u00ab\u00a0, \u00a0\u00bb c:\/testgpkg\/PolysENC.gpkg \u00ab\u00a0].<\/p>\n\n\n\n con los geopackages de importaci\u00f3n que hayas elegido, y la l\u00ednea<\/p>\n\n\n\n output_gpkg = ogr.Open(\u00ab\u00a0c:\/testgpkg\/ENC.gpkg\u00a0\u00bb, 1)<\/p>\n\n\n\n con el geopackage de base de datos que ha creado siguiendo las instrucciones en Creaci\u00f3n de una base de datos Geopackage para mapas ENC (parte 1: creaci\u00f3n de la base de datos).<\/a><\/p>\n\n\n\n Una vez ejecutado este script, su base de datos se actualizar\u00e1 con los nuevos archivos S57.<\/p>\n\n\n\n Actualizaci\u00f3n de una base de datos de geopacks V1 o V2 a la versi\u00f3n 3<\/h2>\n\n\n\nEl procedimiento que acabamos de describir puede utilizarse para crear una nueva base de datos de geopacks. Pero si ya tienes una base de datos creada, no hace falta que vuelvas a hacer todo el trabajo. Aqu\u00ed le damos el procedimiento y las herramientas que necesita para actualizar su estructura (POSAC y QUAPOSQ) y rellenar los valores de estos atributos para las tablas existentes.<\/p>\n\n\n\nDescargue aqu\u00ed los archivos del procedimiento.<\/a><\/p>\n\n\n\n Las distintas etapas de esta actualizaci\u00f3n son las siguientes:<\/p>\n\n\n\n 1- Desde un directorio que contenga todos los archivos .000 que ya se han cargado en la base de datos, ejecute un archivo .bat, upgrade_V2_to_V3.bat, que crea un geopackage update.gpkg que contiene las primitivas (IsolatedNode, ConnectedNode,Edge) de los archivos S57.<\/p>\n\n\n\n .\\ upgrade_V2_to_V3.ba directorio_archivos_000 ruta\/update.gpkg<\/em><\/strong><\/p>\n\n\n\n 2- Ejecuta el script append_tables_upgrade.py en la consola Python de QGis. Aseg\u00farese de cambiar la ruta para el archivo update.gpkg<\/p>\n\n\n\n input_geopackages = [\u00ab\u00a0c:\/testgpkgV2\/update.gpkg\u00a0\u00bb]<\/em><\/strong><\/p>\n\n\n\n y para su base de datos ENC<\/p>\n\n\n\n output_gpkg = ogr.Open(\u00ab\u00a0c:\/enc_db\/ENC4.gpkg\u00a0\u00bb, 1)<\/em><\/strong><\/p>\n\n\n\n Este script a\u00f1ade las primitivas a la base de datos existente.<\/p>\n\n\n\n 3- En la consola Python de QGis, ejecuta el script add_fields_upgrade.py. Aseg\u00farate de cambiar la ruta al archivo .gpkg de tu base de datos<\/p>\n\n\n\n geopackage = ogr.Open(\u00ab\u00a0c:\/enc_db\/ENC4.gpkg\u00a0\u00bb, 1)<\/em><\/strong><\/p>\n\n\n\n Este script crea los atributos POSACC <\/em><\/strong>y QUAPOS <\/em><\/strong>en todas las tablas de la base de datos.<\/p>\n\n\n\n 4-En la consola Python de QGis, ejecute el script add_posacc_quapos_upgrade_points.py<\/em><\/strong>. Aseg\u00farese de cambiar la ruta al archivo .gpkg de su base de datos y a\u00f1ada tablas para excluir si ha creado otras tablas en la base de datos.<\/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 Este script rellenar\u00e1 los campos POSACC <\/em><\/strong>y QUAPOS <\/em><\/strong>en todas las tablas de tipo Punto con los valores contenidos en la primitiva IsolatedNode.<\/p>\n\n\n\n 5-En la consola Python de QGis, ejecute el script add_posacc_quapos_upgrade_lines.py. Aseg\u00farese de cambiar la ruta al archivo .gpkg de su base de datos y a\u00f1ada tablas a excluir si ha creado otras tablas en la base de datos.<\/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 Este script rellenar\u00e1 los campos POSACC y QUAPOS en todas las tablas de tipo L\u00ednea con los valores contenidos en la primitiva Edge.<\/p>\n","protected":false},"excerpt":{"rendered":" El proyecto de integraci\u00f3n de los mapas ENC en QGis en forma de Geopackage ha finalizado. Este art\u00edculo completa los dos art\u00edculos publicados anteriormente sobre la configuraci\u00f3n y gesti\u00f3n de la base de datos Geopackage. En…<\/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_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":[1357],"tags":[1369,2739,2741,2743],"class_list":["post-12825","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-qgis-es","tag-qgis-es","tag-s57-es","tag-svg-es","tag-symbologia-es"],"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-3kR","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/12825","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=12825"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/12825\/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=12825"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=12825"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=12825"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}