sondes bathym\u00e9triques<\/div>

<\/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><\/p>\n\n\n\n

Recuperaci\u00f3n de identificadores \u00ab\u00a0padre<\/h2>\n\n\n\n
Algunos objetos pueden agruparse por un objeto padre no geogr\u00e1fico. Por ejemplo, los sectores de incendio se codifican con un sector por registro en la tabla \u00ab\u00a0LUCES\u00a0\u00bb. Los distintos sectores de un mismo incendio no tienen ning\u00fan atributo particular que pueda indicar que corresponden a la misma entidad. Esta informaci\u00f3n figura en un registro \u00ab\u00a0padre\u00a0\u00bb. Para recuperar este identificador como atributo de la tabla \u00ab\u00a0LIGTHS\u00a0\u00bb, es necesario a\u00f1adir una opci\u00f3n a la l\u00ednea de comandos:<\/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><\/p>\n\n\n\n
Esta opci\u00f3n crea un atributo name_rcid com\u00fan a todos los sectores de un mismo sem\u00e1foro, pero tambi\u00e9n crea una serie de campos (name_rcnm,ORNT,USAG,MASK).<\/p>\n\n\n\n
Tratamiento de tipo \u00ab\u00a0lista\u00a0\u00bb<\/h2>\n\n\n\n
Adem\u00e1s de los tipos de campo tradicionales (entero, real, texto), el formato S57 utiliza un tipo especial: las listas de cadenas de caracteres o enteros. Estos tipos de campo se encuentran en PostgreSQL, pero no en shapefiles y geopackage.<\/p>\n\n\n\n
En ausencia de opciones espec\u00edficas, aparecer\u00e1n mensajes de advertencia indicando que el formato no est\u00e1 soportado y el contenido de estos campos ser\u00e1 ignorado. Para evitar perder el contenido, es necesario a\u00f1adir la extensi\u00f3n :<\/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><\/p>\n\n\n\n
Esta opci\u00f3n toma una lista {..,\u2026} y genera una cadena de texto con un primer n\u00famero que indica el n\u00famero de elementos de la lista, dos puntos y, a continuaci\u00f3n, los elementos de la lista separados por comas:<\/p>\n\n\n\n
{3} -> ‘(1:3)’<\/p>\n\n\n\n
{3,1} -> ‘(2:3,1)’<\/p>\n\n\n\n
Esto complica inevitablemente un poco el tratamiento de estos campos, pero la informaci\u00f3n no se pierde. Por otro lado, el controlador GDAL presenta un problema para gestionar el formato String resultante.<\/p>\n\n\n\n
En primer lugar, no hay ning\u00fan problema con los atributos de los objetos S57. El campo de la tabla geopackage es de tipo ‘Texto’, lo que permite crear una cadena de caracteres de cualquier longitud.
Sin embargo, para los \u00ab\u00a0identificadores padre\u00a0\u00bb, la definici\u00f3n utilizada por el controlador es \u00ab\u00a0Text()\u00a0\u00bb, es decir, una cadena de caracteres con una longitud m\u00e1xima. Esta longitud viene determinada a priori por el contenido del primer registro procesado. En otras palabras, s\u00f3lo por casualidad tendr\u00e1 un campo con la longitud necesaria para incluir todos los valores del fichero.<\/p>\n\n\n\n
Ver\u00e1 entonces mensajes del tipo<\/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 soluci\u00f3n no es complicada, pero implica bastante trabajo. Partimos de la base de que para la visualizaci\u00f3n cl\u00e1sica de las tarjetas ENC, esto no representa ning\u00fan problema y puedes simplemente ignorar estos mensajes de advertencia.<\/p>\n\n\n\n
Sin embargo, si necesitas un procesamiento m\u00e1s avanzado que incluya esta informaci\u00f3n, esta es la \u00fanica manera que he encontrado hasta ahora:<\/p>\n\n\n\n
<\/div>
r\u00e9cup\u00e9ration des values<\/div>
<\/div><\/div>
<\/p>\n\n\n\n
\n
Ejecuta el comando ogr2ogr para crear el geopackage de importaci\u00f3n<\/li>\n\n\n\n
Localiza el valor m\u00e1ximo indicado para cada tipo de campo<\/li>\n\n\n\n
Para la tabla o tablas de salida que se han creado y de las que se quiere recuperar toda la informaci\u00f3n, abre el gestor de base de datos QGis<\/li>\n\n\n\n
Elimina todos los registros de la tabla (DELETE FROM nombre_tabla)<\/li>\n\n\n\n
V\u00e9 al men\u00fa Tabla->Modificar una tabla\u2026<\/li>\n\n\n\n
Selecciona el campo que desees de la lista de campos, luego haz clic en Editar una columna<\/li>\n\n\n\n
Introduce el valor deseado en el campo Longitud, luego haz clic en OK<\/li>\n\n\n\n
Haz lo mismo para todos los dem\u00e1s campos que te interesen<\/li>\n\n\n\n
Elimina todas las dem\u00e1s tablas del geopackage de importaci\u00f3n que no planteen problemas<\/li>\n\n\n\n
Vuelve a ejecutar el comando de importaci\u00f3n ogr2ogr.<\/li>\n<\/ol>\n\n\n\n
<\/div><\/div><\/p>\n\n\n\n
Flujo de trabajo para cargar un \u00fanico archivo S57<\/h1>\n\n\n\n
Creaci\u00f3n del geopackage ENC<\/h2>\n\n\n\n
Vamos a crear un geopackage vac\u00edo que recibir\u00e1 las tablas importadas, y luego las sucesivas actualizaciones cuando se carguen otros archivos S57.<\/p>\n\n\n\n
En el panel Explorador de QGis, haz clic con el bot\u00f3n derecho del rat\u00f3n en Geopackages <\/em><\/strong>y selecciona Crear una base de datos<\/em><\/strong><\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
En la ventana que se abre, introduce el nombre del geopackage, deja la opci\u00f3n por defecto para el nombre de la tabla y selecciona Sin<\/em><\/strong> geometr\u00eda <\/em><\/strong>en Tipo de geometr\u00eda.<\/em><\/strong><\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
Haz clic en Aceptar<\/em><\/strong> para que el nuevo geopackage vac\u00edo se a\u00f1ada a las conexiones del geopackage de QGis. Es aconsejable borrar la tabla ENC vac\u00eda creada al crear el geopackage. Para ello, abre el gestor de bases de datos de QGis desde el men\u00fa principal, sit\u00faate sobre la tabla ENC del geopackage ENC.gpkg y haz clic con el bot\u00f3n derecho del rat\u00f3n para seleccionar Borrar\u2026<\/p>\n\n\n\n
comandos ogr2ogr para crear geopaquetes de importaci\u00f3n<\/h2>\n\n\n\n
Con todos estos elementos en mente, he aqu\u00ed las tres l\u00edneas de comandos ogr2ogr para crear las tablas de los geopaquetes de importaci\u00f3n:<\/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><\/p>\n\n\n\n Para ejecutar estas l\u00edneas de comando, simplemente abre la ventana shell de OSGeo4W:<\/p>\n\n\n\n <\/a><\/figure>\n\n\n\nSe abre la ventana Shell<\/p>\n\n\n\n<\/a><\/figure>\n\n\n\nIntroduce las l\u00edneas de comando, teniendo cuidado de introducir la ruta y el nombre del archivo con la extensi\u00f3n .000.<\/p>\n\n\n\n El resultado ser\u00e1 una serie de tablas creadas en cada geopackage de importaci\u00f3n. Sin embargo, algunas tablas estar\u00e1n completamente vac\u00edas. Si el tipo de geometr\u00eda es posible para una clase de objeto, la tabla se crear\u00e1, pero si no hay ocurrencias en el archivo S57 procesado, la tabla no tendr\u00e1 registros.<\/p>\n\n\n\nUso de scripts Python para gestionar la base de datos<\/h2>\n\n\n\nUtilizaremos la consola Python de QGis para todas las operaciones de gesti\u00f3n de la base de datos.<\/p>\n\n\n\n Para abrirla, vaya al men\u00fa Extensiones -> Consola Python<\/em><\/strong>. Se abrir\u00e1n los paneles de la consola Python:<\/p>\n\n\n\n<\/a><\/figure>\n\n\n\nPara ejecutar un script Python, carga el c\u00f3digo en el panel Editor. Puedes copiar y pegar o utilizar el icono Directorio para apuntar al archivo .py. La sangr\u00eda es crucial para los scripts Python, as\u00ed que elige la segunda opci\u00f3n.<\/p>\n\n\n\nElimine las tablas vac\u00edas y a\u00f1ada el nombre y la escala del mapa<\/h2>\n\n\n\nComo hemos visto antes, al ejecutar los comandos ogr2ogr, se crean tablas vac\u00edas para las capas en las que no se solicita geometr\u00eda en la l\u00ednea de comandos. Por tanto, eliminaremos estas tablas vac\u00edas y aprovecharemos el script de Python para a\u00f1adir dos atributos a todas las tablas: el nombre del mapa ENC y su escala.<\/p>\n\n\n\n En cuanto al nombre, puede ser \u00fatil en el futuro mantenimiento de la base de datos para poder seleccionar las entidades correspondientes a un fichero fuente S57. Esta informaci\u00f3n no figura en ning\u00fan atributo del fichero S57. Por lo tanto, hay que introducirla manualmente.<\/p>\n\n\n\n Lo mismo ocurre con la escala del mapa. Aunque en el S57 se proporciona un atributo para introducir la escala del mapa \u00ab\u00a0de papel\u00a0\u00bb, el atributo CSCALE de la tabla M_CSCL, es muy raro que esta tabla y su atributo est\u00e9n presentes en los ficheros S57.<\/p>\n\n\n\n Si bien el nombre del mapa tiene una utilidad relativa, no puede decirse lo mismo de la escala, ya que se trata de un criterio esencial en la b\u00fasqueda de duplicados. En efecto, cuando se mezclan entidades procedentes de varios ficheros S57, los datos de diferentes escalas coexistir\u00e1n m\u00e1s o menos felizmente.<\/p>\n\n\n\n El siguiente c\u00f3digo Python suprime las tablas vac\u00edas de un esquema de importaci\u00f3n y a\u00f1ade dos atributos (enc_chart y scale) a todas las tablas:<\/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><\/p>\n\n\n\n Ten cuidado de respetar la sangr\u00eda del c\u00f3digo Python. Para simplificar la tarea, descarga el archivo .py directamente desde este enlace: delete_empty_tables_in_geopackages.py<\/a>.<\/p>\n\n\n\n Para ejecutarlo, aseg\u00farate de cambiar los par\u00e1metros de entrada del script:<\/p>\n\n\n\n <\/a><\/figure>\n\n\n\nS\u00f3lo tendr\u00e1s que introducir las geopackage_paths la primera vez que ejecutes el script. Por otro lado, el nombre del mapa ENC y la escala tendr\u00e1n que cambiarse cuando se importe cada nuevo archivo S57. Al final de esta ejecuci\u00f3n tendr\u00e1s los nuevos atributos en cada tabla<\/p>\n\n\n\n<\/a><\/figure>\n\n\n\nFlujo de trabajo para cargar varios archivos S57 al mismo tiempo<\/h1>\n\n\n\nLas l\u00edneas de comando de los p\u00e1rrafos anteriores suponen que tienes un \u00fanico archivo S57 para cargar. Pero puedes cargar varios archivos simult\u00e1neamente. Utilizando las opciones -append -update<\/em><\/strong>, puedes procesar todos los archivos .000 de un directorio en un bucle. Entonces tendr\u00e1s todos los archivos cargados en los geopacks de importaci\u00f3n.<\/p>\n\n\n\n A diferencia de la carga de un \u00fanico archivo, cuando cargamos un lote de archivos necesitamos incluir la adici\u00f3n del nombre del archivo y la escala en el bucle de transcodificaci\u00f3n S57->geopackage. Para ello, utilizaremos un script python que a\u00f1ade los atributos enc_chart y scale a las tablas creadas por el comando ogr2ogr y las rellena con el nombre del archivo correspondiente y una escala definida por el productor. El script tambi\u00e9n elimina cualquier tabla vac\u00eda en el geopackage de importaci\u00f3n.<\/p>\n\n\n\n Gracias a Janez AVSEC por la idea de utilizar la tabla DSID para recuperar la escala de compilaci\u00f3n de datos.<\/p>\n\n\n\n <\/div>Actualizaci\u00f3n disponible<\/div><\/div><\/div><\/p>\n\n\n\nUna actualizaci\u00f3n de este procedimiento est\u00e1 disponible en el art\u00edculo aqu\u00ed.<\/a><\/p>\n\n\n\n <\/div><\/div><\/p>\n\n\n\n Los archivos .bat deben ser ejecutados desde una ventana de OSGeo4W para poder beneficiarse del entorno Python de QGis.<\/p>\n\n\n\n La llamada a cada archivo .bat es de la forma<\/p>\n\n\n\n .\\fichier.bat \u00ab\u00a0directorio que contiene los archivos S57 a cargar\u00a0\u00bb \u00ab\u00a0nombre del geopackage de importaci\u00f3n\u00a0\u00bb <\/em><\/strong><\/p>\n\n\n\n Por ejemplo:<\/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><\/p>\n\n\n\n Los archivos .bat para esta tarea son los siguientes:<\/p>\n\n\n\n Para importar tablas de tipo de punto:<\/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><\/p>\n\n\n\n Observar\u00e1 que, adem\u00e1s del comando ogr2ogr para importar capas de puntos, importamos dos tablas sin geometr\u00eda: DSID y C_AGGR. La tabla DSID contiene dos atributos con el nombre del archivo S57 y la escala de compilaci\u00f3n de datos. Esta tabla ser\u00e1 utilizada por el script de Python (update_geopackage_dsid.py) para introducir el archivo S57 original y la escala de compilaci\u00f3n en cada tabla de geopackage.<\/p>\n\n\n\n Para importar tablas de tipo l\u00ednea:<\/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><\/p>\n\n\n\n Para importar tablas de pol\u00edgonos:<\/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><\/p>\n\n\n\n Los archivos .bat contienen el comando ogr2ogr y una l\u00ednea para ejecutar un script de Python, update_geopackage_dsid.py. Aseg\u00farate de cambiar la ruta de este script en el archivo .bat por la de tu elecci\u00f3n. Este script a\u00f1ade los dos atributos enc_chart<\/em><\/strong> y scale <\/em><\/strong>a todas las tablas del geopackage.<\/p>\n\n\n\n El c\u00f3digo python utilizado es el siguiente:<\/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><\/p>\n\n\n\n <\/p>\n\n\n\n Podr\u00edas hacer un \u00fanico archivo .bat para los tres procesos, pero necesitas comprobar los resultados para cada tipo de geometr\u00eda.<\/p>\n\n\n\n Puedes descargar los archivos .bat y el c\u00f3digo Python haciendo clic aqu\u00ed.<\/a><\/p>\n\n\n\n <\/p>\n\n\n\n Contin\u00faa independientemente del n\u00famero de archivos cargados<\/h1>\n\n\n\nClonar y\/o a\u00f1adir tablas importadas al geopackage ENC<\/h2>\n\n\n\nPara la primera carga del S57, el siguiente script de Python crea todas las tablas presentes en los tres geopacks de importaci\u00f3n. Para las cargas siguientes, 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 La clonaci\u00f3n se realiza mediante el siguiente c\u00f3digo:<\/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><\/p>\n\n\n\n <\/p>\n\n\n\n Puedes descargar el archivo .py aqu\u00ed: clone_or_append_tables_with_prefix.py.<\/a><\/p>\n\n\n\n Gesti\u00f3n de duplicados<\/h2>\n\n\n\nCuando se superponen dos mapas ENC, las entidades presentes en la zona de superposici\u00f3n se insertar\u00e1n tantas veces como superposiciones haya. Para evitar duplicados y m\u00e1s, la base de datos debe limpiarse despu\u00e9s de cada importaci\u00f3n de un fichero S57.<\/p>\n\n\n\n Pero hay dos tipos diferentes de duplicaci\u00f3n. Por un lado, las zonas en las que se solapan dos mapas de escala similar producir\u00e1n duplicados \u00ab\u00a0reales\u00a0\u00bb que habr\u00e1 que limpiar. Por otro lado, los mapas superpuestos de escalas muy diferentes deben tratarse de forma diferente. El mapa menos detallado (escala peque\u00f1a) tendr\u00e1 una selecci\u00f3n de datos del mapa m\u00e1s detallado (escala grande). En este segundo caso, no es aconsejable eliminar esta informaci\u00f3n.<\/p>\n\n\n\n El c\u00f3digo Python siguiente eliminar\u00e1 los duplicados \u00ab\u00a0reales\u00a0\u00bb. Para otros duplicados, recomendamos utilizar filtros QGis. Dependiendo de sus necesidades, puedes utilizar el atributo \u00ab\u00a0escala\u00a0\u00bb que a\u00f1adimos al eliminar las tablas vac\u00edas para definir qu\u00e9 datos se muestran en su mapa y cu\u00e1les no.<\/p>\n\n\n\n<\/a><\/figure>\n\n\n\nEn este ejemplo, s\u00f3lo se mostrar\u00e1n en el mapa los elementos correspondientes a una escala comprendida entre 12500 y 45000.<\/p>\n\n\n\n Esta operaci\u00f3n corresponde a filtrar una capa, pero\u2026 \u00bfqu\u00e9 haces cuando tienes 130 capas desplegadas? Para ahorrarte tener que realizar esta operaci\u00f3n 130 veces, aqu\u00ed tienes un c\u00f3digo Python que te permite filtrar todas las capas del proyecto:<\/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><\/p>\n\n\n\n Y otra para eliminar el filtro de todas las capas mostradas:<\/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><\/p>\n\n\n\n <\/p>\n\n\n\n Descarga los c\u00f3digos python aqu\u00ed.<\/a><\/p>\n\n\n\n El siguiente c\u00f3digo crea una tabla temporal (temp_table) para cada tabla del geopackage ENC. A continuaci\u00f3n, rellena esta tabla con los registros separados basados en LNAM. A continuaci\u00f3n, borra la tabla original y renombra la tabla temporal con el nombre de la tabla original.<\/p>\n\n\n\n Este tratamiento se aplica a todas las tablas S57 excepto:<\/p>\n\n\n\n \nlas tablas de luces (LIGHTS) y de sondeos batim\u00e9tricos (SOUNDG), en las que el atributo LNAM no es \u00fatil. En efecto, cada sector de incendio tiene el mismo LNAM y name_rcid pero corresponde a un color y a un sector geogr\u00e1fico diferentes. Del mismo modo, la LNAM y el name_rcid de los sondeos batim\u00e9tricos son los mismos para todos los puntos de sondeo de la entidad multipunto original. En el primer caso tenemos que buscar colores y sectores duplicados para el mismo name_rcid, y para los sondeos tenemos que encontrar los sondeos con el mismo XY en su geometr\u00eda.<\/li>\n\n\n\n Tablas SEAARE, LNDRGN y ADMARE que cubren grandes \u00e1reas y que se dividen seg\u00fan la huella del mapa, manteniendo el mismo LNAM. En este caso, aunque la LNAM sea la misma, los pol\u00edgonos no tienen la misma forma y \u00e1rea.<\/li>\n\n\n\nlas tablas DSID y C_AGGR, que no tienen LNAM y que, por principio, no pueden tener duplicados.<\/li>\n<\/ul>\n\n\n\n<\/div>effacer doublons .py<\/div><\/div><\/div> <\/p>\n\n\n\nfrom osgeo import ogr<\/code><\/pre>\n\n\n\ndef delete_doublons(): # Ouvre le GeoPackage en mode mise \u00e0 jour geopackage_path = \"c:\/enc_db\/enc2.gpkg\" geopackage = ogr.Open(geopackage_path, update=True)<\/code><\/pre>\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><\/p>\n\n\n\n <\/p>\n\n\n\n Puedes descargar el archivo .py desde este enlace: supprime_doublons.py.<\/a> No olvides cambiar las rutas y los nombres de archivo para que coincidan con los tuyos. F\u00edjate tambi\u00e9n en la l\u00ednea :<\/p>\n\n\n\n if table_name not in [\u00ab\u00a0pl_SEAARE\u00a0\u00bb, \u00ab\u00a0pl_LNDRGN\u00a0\u00bb, \u00ab\u00a0pl_ADMARE\u00a0\u00bb, \u00ab\u00a0layer_styles\u00a0\u00bb, \u00ab\u00a0natsurf\u00a0\u00bb,\u00a0\u00bbDSID\u00a0\u00bb,\u00a0\u00bbC_AGGR\u00a0\u00bb]:<\/em><\/strong><\/p>\n\n\n\n En principio, s\u00f3lo tendr\u00e1s las tablas prefijadas y las tablas layer_styles y natsurf. Esta l\u00ednea evita que se procesen pl_SEAARE, pl_LNDRGN y pl_ADMARE, pero sobre todo evita que se vac\u00eden las tablas layer_styles y natsurf que utilizamos para las simbolog\u00edas. Si necesitas crear tablas distintas a las del fichero s57 del geopackage, debes tener en cuenta los nombres de estas tablas en este script, a\u00f1adi\u00e9ndolos a esta l\u00ednea de excepciones.<\/p>\n\n\n\n <\/p>\n","protected":false},"excerpt":{"rendered":" La primera parte del proyecto Financiaci\u00f3n Colaborativa para la Integraci\u00f3n de Datos Marinos en QGIS ha llegado a su fin, gracias a las contribuciones de ORANGE Marine, Geoconceptos Uruguay, Janez Avzec y Gaetan Mas. Nuestro m\u00e1s…<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"_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_feature_clip_id":0,"_jetpack_memberships_contains_paid_content":false,"footnotes":"","jetpack_post_was_ever_published":false},"categories":[1357],"tags":[1989,2801,2797,2739,2799,2743,2032],"class_list":["post-12511","post","type-post","status-publish","format-standard","hentry","category-qgis-es","tag-geopackage-es","tag-ogr2ogr-es","tag-python-es","tag-s57-es","tag-script-es","tag-symbologia-es","tag-tutorial"],"aioseo_notices":[],"aioseo_head":"\n\t\t\n\t\n\t\n\t\n\t\n\t\n\t\n\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t