\n Documentos e ferramentas do Geopackage ENC<\/a>\n <\/div>\n <\/div>\n <\/div>\n<\/div>\n\n<\/div>.<\/p>\n\n\n\n

Ele lhe dar\u00e1 uma vis\u00e3o geral completa e coerente de todas as etapas envolvidas na cria\u00e7\u00e3o e no carregamento de um banco de dados ENC no QGis em formato de geopackage e de todos os elementos necess\u00e1rios (.bat, .py, .svg) para configurar uma simbologia autom\u00e1tica apropriada.<\/p>\n\n\n\n

<\/div><\/div>\n\n\n\n

O projeto para integrar os mapas ENC ao QGis na forma de Geopackage est\u00e1 conclu\u00eddo. Este artigo completa os dois artigos publicados anteriormente sobre a configura\u00e7\u00e3o e o gerenciamento do banco de dados do Geopackage.<\/p>\n\n\n\n

No primeiro artigo: Cria\u00e7\u00e3o de um banco de dados Geopackage para mapas ENC (Parte 1: cria\u00e7\u00e3o do banco de dados)<\/a>, voc\u00ea encontrar\u00e1 os fundamentos do trabalho com geopacotes para integrar dados de arquivos S57 (ENC). Voc\u00ea deve come\u00e7ar por este artigo para entender o procedimento de importa\u00e7\u00e3o de arquivos S57 com o ogr2ogr, bem como as v\u00e1rias etapas envolvidas na formata\u00e7\u00e3o dos dados. Voc\u00ea tamb\u00e9m poder\u00e1 fazer download de todas as ferramentas (arquivos de comando .bat, scripts python) necess\u00e1rias para concluir essas etapas.<\/p>\n\n\n\n

Ap\u00f3s a publica\u00e7\u00e3o desse artigo, continuamos a trabalhar em um procedimento equivalente, mas usando um banco de dados PostgreSQL\/postgis. Isso nos permitiu encontrar v\u00e1rias melhorias no fluxo de trabalho inicial.<\/p>\n\n\n\n

Um segundo artigo, Update(1): banco de dados ENC com Geopackage no QGis<\/a>, atualiza o procedimento de importa\u00e7\u00e3o (arquivos .bat com comandos ogr2ogr) para facilitar o gerenciamento de informa\u00e7\u00f5es duplicadas. Essas duplica\u00e7\u00f5es ocorrem quando mapas em escalas diferentes s\u00e3o integrados ao mesmo banco de dados. A adi\u00e7\u00e3o de um atributo de \u00ab\u00a0finalidade\u00a0\u00bb ao mapa ENC significa que somente informa\u00e7\u00f5es exclusivas podem ser exibidas para cada finalidade do mapa.<\/p>\n\n\n\n

Como o objetivo do projeto era produzir uma simbologia no QGis que fosse equivalente \u00e0 dos mapas em papel, ainda havia um ponto a ser resolvido: levar em conta o par\u00e2metro de \u00ab\u00a0qualidade\u00a0\u00bb da posi\u00e7\u00e3o dos objetos geogr\u00e1ficos. Esse \u00e9 o principal ponto abordado neste artigo.<\/p>\n\n\n\n

Os tr\u00eas artigos s\u00e3o complementares em termos de conte\u00fado, mas para os arquivos disponibilizados, baixe apenas os arquivos sugeridos neste artigo. Eles terminam com ‘V3’ e certifique-se de que est\u00e1 usando a vers\u00e3o mais recente.<\/em><\/strong><\/p>\n\n\n\n

O formato S57 e o GIS<\/h2>\n\n\n\n
Se voc\u00ea vai come\u00e7ar a trabalhar com arquivos S57 em um GIS, h\u00e1 alguns aspectos que voc\u00ea precisa conhecer para navegar sem problemas.<\/p>\n\n\n\n
Em primeiro lugar, a estrutura dos arquivos S57 n\u00e3o corresponde \u00e0s estruturas adotadas no GIS.<\/p>\n\n\n\n
Em um GIS, voc\u00ea tem um objeto geogr\u00e1fico que \u00e9 representado por uma tabela com dois tipos de informa\u00e7\u00f5es: a geometria das entidades no objeto e os atributos dessas entidades.<\/p>\n\n\n\n
Se voc\u00ea tiver outros objetos com geometrias id\u00eanticas, as informa\u00e7\u00f5es geom\u00e9tricas ser\u00e3o duplicadas, uma vez em cada tabela.<\/p>\n\n\n\n
No formato S57, o principal objetivo \u00e9 otimizar o armazenamento de informa\u00e7\u00f5es e, portanto, n\u00e3o duplicar informa\u00e7\u00f5es. Se um objeto tiver uma entidade de ponto, ser\u00e1 criado um ponto. Se outros objetos tiverem entidades localizadas nesse ponto, ser\u00e1 usada a refer\u00eancia do ponto que j\u00e1 foi criado. Dessa forma, um ponto \u00e9 descrito apenas uma vez no arquivo. O mesmo se aplica a polilinhas e superf\u00edcies. Portanto, um arquivo S57 ter\u00e1 uma s\u00e9rie de tabelas contendo informa\u00e7\u00f5es geom\u00e9tricas, conhecidas como \u00ab\u00a0primitivas\u00a0\u00bb:<\/p>\n\n\n\n
\n
IsolatedNode (pontos)<\/li>\n\n\n\n
ConnectedNode (pontos)<\/li>\n\n\n\n
Edge (polilinhas)<\/li>\n\n\n\n
Face (pol\u00edgonos)<\/li>\n<\/ul>\n\n\n\n
A tabela de atributos para os v\u00e1rios objetos S57 cont\u00e9m apenas os atributos dos objetos.<\/p>\n\n\n\n
O que complica a tarefa \u00e9 que h\u00e1 dois atributos que se referem \u00e0s geometrias: posacc <\/em><\/strong>(a precis\u00e3o estimada da posi\u00e7\u00e3o, um valor quantitativo) e quapos <\/em><\/strong>(a qualidade da posi\u00e7\u00e3o, uma vari\u00e1vel qualitativa).<\/p>\n\n\n\n
Esses dois atributos podem ser encontrados nas tabelas primitivas.<\/p>\n\n\n\n
*Para mudar da estrutura S57 para uma estrutura GIS (**shapefile<\/a>, geopackage, postgis), usamos a biblioteca GDAL e seu comando ogr2ogr.<\/p>\n\n\n\n*
Esse comando criar\u00e1 tabelas GIS a partir da estrutura S57, criando uma tabela por objeto S57, atribuindo a geometria correspondente das tabelas primitivas a cada entidade e adicionando os atributos do objeto S57 a cada entidade. O rastreamento das primitivas usadas para a geometria de cada entidade pode ser encontrado no campo NAME_RCID <\/em><\/strong>das tabelas GIS, desde que as op\u00e7\u00f5es -oo \u00ab\u00a0RETURN_LINKAGES=ON\u00a0\u00bb -oo \u00ab\u00a0LNAM_REFS=ON\u00a0\u00bb<\/em><\/strong> tenham sido adicionadas \u00e0 linha de comando do ogr2ogr.<\/p>\n\n\n\n
A figura a seguir mostra uma camada do tipo ponto. O valor indicado no campo NAME_RCID <\/em><\/strong>\u00e9 o RCID <\/em><\/strong>do ponto usado na tabela IsolatedNode.<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
A figura a seguir mostra um exemplo de uma camada do tipo linear. Os valores indicados no campo NAME_RCID s\u00e3o os dos RCIDs das polilinhas usadas na tabela Edge.<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
A figura a seguir mostra um exemplo de uma camada do tipo pol\u00edgono. Os valores indicados no campo NAME_RCID s\u00e3o os das polilinhas usadas na tabela Edge.<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
Para recuperar os atributos QUAPOS <\/em><\/strong>e POSACC <\/strong><\/em>de cada entidade nas tabelas de tipo de ponto, precisamos recuperar os valores do ponto IsolatedNode e atribu\u00ed-los \u00e0s tabelas dos v\u00e1rios objetos ENC.<\/p>\n\n\n\n
Se os identificadores fossem diretamente os RCIDs nas tabelas ENC, poder\u00edamos fazer uma jun\u00e7\u00e3o entre cada tabela (soundg, obstrn,\u2026) e IsolatedNode. Mas, como voc\u00ea pode ver nas imagens anteriores, o atributo NAME_RCID <\/em>\u00e9 do tipo stringlist, o que bloqueia essa solu\u00e7\u00e3o. Portanto, desenvolvemos um script python que faz o trabalho quando os dados do S57->geopackage s\u00e3o carregados.<\/p>\n\n\n\n
Download de arquivos V3<\/h2>\n\n\n\nVoc\u00ea pode baixar todos os arquivos .bat e .py mencionados neste artigo a partir deste link.<\/a><\/p>\n\n\n\n
Modifica\u00e7\u00e3o dos comandos do ogr2ogr para importa\u00e7\u00e3o de arquivos S57<\/h2>\n\n\n\n
Importa\u00e7\u00e3o de camadas de \u00ab\u00a0pontos<\/h3>\n\n\n\n
Para importar primitivas de arquivos S57, \u00e9 necess\u00e1rio adicionar a op\u00e7\u00e3o -oo RETURN_PRIMITIVES = ON<\/em><\/strong> aos comandos de importa\u00e7\u00e3o. No caso de pontos, voc\u00ea tamb\u00e9m precisa adicionar a op\u00e7\u00e3o -nlt MULTIPOINT<\/em><\/strong>, caso contr\u00e1rio, a tabela IsolatedNode ser\u00e1 criada como um ponto e n\u00e3o haver\u00e1 suporte para camadas multiponto.<\/p>\n\n\n\n
O arquivo .bat para importa\u00e7\u00e3o de camadas do tipo ponto ser\u00e1 ent\u00e3o :<\/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 execut\u00e1-lo, abra uma janela do OSGeo4W Shell, v\u00e1 para o diret\u00f3rio que cont\u00e9m os arquivos .bat e digite o comando:<\/p>\n\n\n\n .\\mass_load_s57_points_V3.bat directory_files_000 path\/pointsENC.gpkg<\/em><\/strong><\/p>\n\n\n\n Primeiro, edite o arquivo para alterar as linhas:<\/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 e<\/p>\n\n\n\n python c:\/testgpkgV3\/add_posacc_quapos_to_pointsV3.py<\/em><\/strong><\/p>\n\n\n\n Para inserir os caminhos que voc\u00ea escolheu para armazenar os scripts python.<\/p>\n\n\n\n <\/a><\/figure>\n\n\n\nO primeiro script, update_geopackage_dsid_V3.py<\/em><\/strong>, foi modificado para adicionar os atributos POSACC e QUAPOS a todas as tabelas de pontos, al\u00e9m de enc_chart, escala e finalidade, quando a troca S57->geopackage \u00e9 feita.<\/p>\n\n\n\n O segundo script, add_posacc_quapos_to_pointsV3.py<\/em><\/strong>, uma vez que todos os arquivos S57 no diret\u00f3rio ENC tenham sido carregados, atualiza os valores desses atributos para cada entidade nas tabelas de pontos.<\/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 No exemplo a seguir, carregamos 7 mapas ENC no geopackage de importa\u00e7\u00e3o pointsENC.gpkg.<\/p>\n\n\n\n<\/a><\/figure>\n\n\n\nVoc\u00ea pode ver que as tabelas do geopacote cont\u00eam os atributos POSACC e QUAPOS e que o \u00faltimo script preencheu os valores correspondentes a esses campos.<\/p>\n\n\n\n Na segunda parte dessa vers\u00e3o final, dedicada \u00e0 simbologia, voc\u00ea ver\u00e1 o novo gerenciamento da camada \u00ab\u00a0nature of seabed\u00a0\u00bb (natureza do fundo do mar). Essa simbologia se baseia na cria\u00e7\u00e3o de um atributo \u00ab\u00a0Label\u00a0\u00bb na camada SBDARE e no seu preenchimento com um script python Label_sbdare.py<\/em><\/strong> . Essa opera\u00e7\u00e3o pode ser executada no banco de dados completo ou aqui, no geopackage pointsENC.gpkg. Como o c\u00e1lculo leva algum tempo, recomendamos que voc\u00ea o fa\u00e7a nesse est\u00e1gio. Consulte o artigo Mapas do geopacote ENC na vers\u00e3o final do QGis: parte 2<\/em><\/strong> para obter instru\u00e7\u00f5es.<\/p>\n\n\n\nImporta\u00e7\u00e3o de camadas de \u00ab\u00a0linha<\/h2>\n\n\n\nPara importar primitivas de arquivos S57, \u00e9 necess\u00e1rio adicionar a op\u00e7\u00e3o -oo RETURN_PRIMITIVES = ON<\/em><\/strong> aos comandos de importa\u00e7\u00e3o.<\/p>\n\n\n\n O arquivo .bat de importa\u00e7\u00e3o de camada de linha ser\u00e1 ent\u00e3o :<\/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 execut\u00e1-lo, abra uma janela do OSGeo4W Shell, v\u00e1 para o diret\u00f3rio que cont\u00e9m os arquivos .bat e digite o comando:<\/p>\n\n\n\n .\\mass_load_s57_lines_V3.bat directory_files_000 path\/pointsENC.gpkg<\/em><\/strong><\/p>\n\n\n\n Primeiro, edite o arquivo para alterar as linhas:<\/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 e<\/p>\n\n\n\n python c:\/testgpkgV3<\/em><\/strong>\/add_posacc_quapos_to_linesV3.py<\/p>\n\n\n\n para inserir os caminhos que voc\u00ea escolheu para armazenar os scripts python.<\/p>\n\n\n\n O segundo script, add_posacc_quapos_to_linesV3.py<\/em><\/strong>, depois que todos os arquivos S57 no diret\u00f3rio ENC tiverem sido carregados, atualiza os valores desses atributos para cada entidade nas tabelas de linhas.<\/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 No caso de tabelas de tipo de ponto, o v\u00ednculo entre o NAME_RCID das tabelas e o RCID da tabela IsolatedNode n\u00e3o levanta d\u00favidas. No entanto, no caso de polilinhas, para uma entidade geogr\u00e1fica, a lista de polilinhas diz respeito a v\u00e1rias entidades na camada Edge. N\u00f3s nos perguntamos se era poss\u00edvel que as polilinhas que comp\u00f5em uma entidade, por exemplo, DEPCNT, tivessem valores QUAPOS diferentes.<\/p>\n\n\n\n Realizamos o teste em um banco de dados com 650 mapas ENC em todas as escalas. A camada DEPCNT continha 162.585 recursos. Apenas 38.731 tinham RCIDs na lista de NAME_RCIDs na tabela Edges com um valor QUAPOS diferente de zero. Encontramos 193 entidades DEPCNT que tinham valores QUAPOS diferentes em sua lista NAME_RCID. Como isso representa 0,1%, decidimos n\u00e3o levar em conta essa possibilidade, pois a complexidade da desconstru\u00e7\u00e3o de polilinhas \u00e9 extrema. Acontece que a tabela Edge resultante do driver ogr2ogr n\u00e3o apresenta problemas no n\u00edvel do atributo (a uni\u00e3o entre NAME_RCID e RCID \u00e9 boa), mas as geometrias na tabela cont\u00eam quase metade das geometrias inv\u00e1lidas (linhas com um \u00fanico ponto ou geometrias nulas). No entanto, as polilinhas das entidades S57 n\u00e3o apresentam erros.<\/p>\n\n\n\n Portanto, seria necess\u00e1rio reparar e reconstruir todas as geometrias Edge para poder desconstruir as geometrias DEPCNT de modo a n\u00e3o ter QUAPOS diferentes em cada caso.<\/p>\n\n\n\n No exemplo a seguir, carregamos 7 mapas ENC no geopackage de importa\u00e7\u00e3o linesENC.gpkg<\/p>\n\n\n\n <\/a><\/figure>\n\n\n\nImporta\u00e7\u00e3o de camadas de pol\u00edgonos<\/h2>\n\n\n\nEmbora o formato S57 permita que os pol\u00edgonos sejam usados como primitivos, essa op\u00e7\u00e3o n\u00e3o \u00e9 usada. Portanto, n\u00e3o \u00e9 necess\u00e1rio modificar a importa\u00e7\u00e3o de dados de pol\u00edgonos. As primitivas (Face) s\u00e3o solicitadas na linha de comando ogr2ogr, e os campos POSACC e QUAPOS s\u00e3o adicionados \u00e0s tabelas para poss\u00edvel uso futuro, mas n\u00e3o \u00e9 necess\u00e1rio adicionar um script para preencher os valores.<\/p>\n\n\n\nClonagem e\/ou adi\u00e7\u00e3o de tabelas importadas ao geopacote ENC<\/h2>\n\n\n\nPara o primeiro carregamento do S57, o script Python a seguir cria todas as tabelas nos tr\u00eas geopacotes de importa\u00e7\u00e3o. Para cargas subsequentes, se a tabela j\u00e1 existir no geopacote ENC, os registros da tabela no geopacote de importa\u00e7\u00e3o ser\u00e3o adicionados aos registros j\u00e1 presentes no ENC. Se, por outro lado, a tabela ainda n\u00e3o existir, ela ser\u00e1 criada e preenchida com registros do geopacote de importa\u00e7\u00e3o. No arquivo do geopacote ENC, as tabelas s\u00e3o prefixadas com pt_, li_ e pl_ para indicar seu tipo de geometria.<\/p>\n\n\n\n As seguintes altera\u00e7\u00f5es foram feitas no script em compara\u00e7\u00e3o com a V2:<\/p>\n\n\n\n\nO SRC<\/a> das tabelas do geopacote \u00e9 preenchido durante a cria\u00e7\u00e3o. O c\u00f3digo EPSG:4326<\/a> \u00e9 atribu\u00eddo \u00e0s geometrias.<\/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\u00e7ar a cria\u00e7\u00e3o do \u00edndice espacial em todas as tabelas do 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\nA clonagem \u00e9 realizada usando o c\u00f3digo python: clone_or_append_tables_with_prefixV3.py<\/a> Certifique-se de alterar a linha que indica os nomes dos geopacotes que voc\u00ea criou com os arquivos .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 com os geopacotes de importa\u00e7\u00e3o que voc\u00ea escolheu e a linha<\/p>\n\n\n\n output_gpkg = ogr.Open(\u00ab\u00a0c:\/testgpkg\/ENC.gpkg\u00a0\u00bb, 1)<\/p>\n\n\n\n com o geopacote de banco de dados que voc\u00ea criou seguindo as instru\u00e7\u00f5es em Cria\u00e7\u00e3o de um banco de dados do Geopackage para mapas ENC (parte 1: cria\u00e7\u00e3o do banco de dados)<\/a>.<\/p>\n\n\n\n Depois que esse script for executado, seu banco de dados ser\u00e1 atualizado com os novos arquivos S57.<\/p>\n\n\n\n Atualiza\u00e7\u00e3o de um banco de dados do geopacote V1 ou V2 para a vers\u00e3o 3<\/h2>\n\n\n\nO procedimento que acabamos de descrever pode ser usado para criar um novo banco de dados de geopacote. Por\u00e9m, se voc\u00ea j\u00e1 tiver um banco de dados criado, n\u00e3o ser\u00e1 necess\u00e1rio refazer todo o trabalho. Aqui fornecemos o procedimento e as ferramentas necess\u00e1rias para atualizar sua estrutura (POSAC e QUAPOSQ) e preencher os valores desses atributos para as tabelas existentes.<\/p>\n\n\n\nFa\u00e7a o download dos arquivos para o procedimento aqui.<\/a><\/p>\n\n\n\n Os v\u00e1rios est\u00e1gios dessa atualiza\u00e7\u00e3o s\u00e3o os seguintes:<\/p>\n\n\n\n 1- Em um diret\u00f3rio que contenha todos os arquivos .000 que j\u00e1 foram carregados no banco de dados, execute um arquivo .bat, upgrade_V2_to_V3.bat, que cria um geopackage update.gpkg contendo as primitivas (IsolatedNode, ConnectedNode,Edge) dos arquivos S57.<\/p>\n\n\n\n .\\ upgrade_V2_to_V3.ba directory_files_000 path\/update.gpkg<\/em><\/strong><\/p>\n\n\n\n 2- Execute o script append_tables_upgrade.py no console do QGis Python. Certifique-se de alterar o caminho para o arquivo update.gpkg<\/p>\n\n\n\n input_geopackages = [\u00ab\u00a0c:\/testgpkgV2\/update.gpkg\u00a0\u00bb]<\/em><\/strong><\/p>\n\n\n\n e para o seu banco de dados 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 Esse script adiciona as primitivas ao banco de dados existente.<\/p>\n\n\n\n 3- No console do QGis Python, execute o script add_fields_upgrade.py. Certifique-se de alterar o caminho para o arquivo .gpkg do seu banco de dados<\/p>\n\n\n\n geopackage = ogr.Open(\u00ab\u00a0c:\/enc_db\/ENC4.gpkg\u00a0\u00bb, 1)<\/em><\/strong><\/p>\n\n\n\n Esse script cria os atributos POSACC <\/strong>e QUAPOS <\/strong>em todas as tabelas do banco de dados.<\/p>\n\n\n\n 4-No console Python do QGis, execute o script add_posacc_quapos_upgrade_points.py<\/em><\/strong>. Certifique-se de alterar o caminho para o arquivo .gpkg do seu banco de dados e adicione tabelas para excluir se tiver criado outras tabelas no banco de dados.<\/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 Esse script preencher\u00e1 os campos POSACC <\/strong>e QUAPOS <\/strong>em todas as tabelas de tipo de ponto com os valores contidos na primitiva IsolatedNode.<\/p>\n\n\n\n 5-No console do QGis Python, execute o script add_posacc_quapos_upgrade_lines.py. Certifique-se de alterar o caminho para o arquivo .gpkg do seu banco de dados e adicione tabelas para excluir se tiver criado outras tabelas no banco de dados.<\/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 Esse script preencher\u00e1 os campos POSACC <\/strong>e QUAPOS <\/strong>em todas as tabelas de tipo de linha com os valores contidos na primitiva Edge.<\/p>\n","protected":false},"excerpt":{"rendered":" O projeto para integrar os mapas ENC ao QGis na forma de Geopackage est\u00e1 conclu\u00eddo. Este artigo completa os dois artigos publicados anteriormente sobre a configura\u00e7\u00e3o e o gerenciamento do banco de dados do Geopackage. No…<\/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":[2134],"tags":[2160,2731,2733,2725],"class_list":["post-12826","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-qgis","tag-qgis-pt","tag-s57-pt","tag-svg-pt","tag-symbologia"],"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-3kS","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/12826","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=12826"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/12826\/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=12826"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=12826"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=12826"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}