Devido ao volume do resultado, estamos publicando-o em duas partes: a primeira diz respeito \u00e0 cria\u00e7\u00e3o e ao gerenciamento do banco de dados, a segunda \u00e0 simbologia S57 no QGis.<\/em><\/strong><\/p>\n\n\n\n O formato S57 \u00e9 complexo e n\u00e3o est\u00e1 totalmente em conformidade com os padr\u00f5es de formato GIS. Sua convers\u00e3o com o GDAL \u00e9, portanto, complexa ou incompleta. Neste artigo, entraremos em detalhes sobre a transforma\u00e7\u00e3o completa do S57 em tabelas de arquivos do Geopackage<\/a>.<\/p>\n\n\n\n As \u00ab\u00a0camadas\u00a0\u00bb do S57 s\u00e3o classes de objetos. Por exemplo, as diferentes \u00e1reas de terra s\u00e3o codificadas na classe de objeto LNDARE.<\/p>\n\n\n\n A defini\u00e7\u00e3o dessa classe de objeto \u00e9:<\/p>\n\n\n\n Geo object: Land area (LNDARE) (P,L,A) Embora todos os objetos LNDARE tenham os mesmos atributos, o mesmo n\u00e3o pode ser dito sobre o tipo de geometria. A informa\u00e7\u00e3o (P,L,A) indica que pontos, linhas e pol\u00edgonos podem ser encontrados nessa classe de objeto. Ao contr\u00e1rio dos padr\u00f5es do GIS, os tr\u00eas tipos de geometria coexistem na mesma classe de objeto.<\/p>\n\n\n\n Seja qual for o formato que escolhermos para integrar ao QGis, precisaremos criar uma camada para cada tipo de geometria. Se n\u00e3o fizermos isso, o GDAL criar\u00e1 o tipo de camada com base na primeira entidade encontrada durante a convers\u00e3o. Se ela for do tipo ponto, a camada criada ser\u00e1 do tipo ponto e as entidades linhas e pol\u00edgonos ser\u00e3o ignoradas. Da mesma forma, se a primeira entidade for do tipo linha, os pontos e pol\u00edgonos ser\u00e3o ignorados.<\/p>\n\n\n\n Tamb\u00e9m deve ser observado que o formato S57 n\u00e3o tem restri\u00e7\u00f5es sobre classes de objetos: se n\u00e3o houver nenhuma classe de objeto na \u00e1rea coberta pelo arquivo S57, n\u00e3o haver\u00e1 nada sobre ela no arquivo (nenhuma camada vazia). Da mesma forma, se houver apenas um tipo de geometria presente, mesmo que todos os tr\u00eas tipos sejam poss\u00edveis, n\u00e3o haver\u00e1 nenhum vest\u00edgio dos outros tipos de geometria.<\/p>\n\n\n\n Portanto, o processamento de um arquivo S57 com ogr2ogr deve ser dividido em tr\u00eas est\u00e1gios, um para cada tipo de geometria. As op\u00e7\u00f5es a seguir permitem que voc\u00ea processe cada classe de objeto S57 selecionando apenas um tipo de geometria:<\/p>\n\n\n\n -where \u00ab OGR_GEOMETRY=\u2019POINT\u2019 or OGR_GEOMETRY=\u2019MULTIPOINT\u2019 \u00bb<\/strong><\/p>\n\n\n\n -where \u00ab OGR_GEOMETRY=\u2019LINESTRING\u2019 or OGR_GEOMETRY=\u2019MULTILINESTRING\u2019 \u00bb<\/strong><\/p>\n\n\n\n -where \u00ab OGR_GEOMETRY=\u2019POLYGON\u2019 or OGR_GEOMETRY=\u2019MULTIPOLYGON\u2019 \u00bb<\/strong><\/p>\n\n\n\n Para determinados tipos de driver, o GDAL permite que voc\u00ea crie prefixos nas tabelas de sa\u00edda. Nesse caso, voc\u00ea poderia criar todas as tabelas (pontos, linhas, pol\u00edgonos) em um \u00fanico geopackage, prefixando-as com pt_, li_ e pl_, por exemplo. O problema \u00e9 que o driver S57 do GDAL n\u00e3o permite essa op\u00e7\u00e3o. Portanto, precisamos criar tr\u00eas geopacotes separados, nos quais as tabelas ser\u00e3o criadas de acordo com o tipo de geometria. Cada geopacote conter\u00e1 uma tabela com o mesmo nome, mas com uma geometria diferente. Embora essa seja a solu\u00e7\u00e3o mais simples, n\u00e3o \u00e9 a mais elegante nem a mais pr\u00e1tica. Aqui, usaremos tr\u00eas geopacotes de importa\u00e7\u00e3o para os comandos ogr2ogr nos quais importaremos tabelas de um arquivo S57. N\u00f3s os chamaremos de PointsENC, LinesENC e PolysENC. Tamb\u00e9m criaremos um \u00fanico geopackage chamado ENC para todo o banco de dados. Para cada mapa ENC, importaremos seu conte\u00fado para os tr\u00eas geopacotes de importa\u00e7\u00e3o usando ogr2ogr e, em seguida, executaremos scripts Python para atualizar o banco de dados do geopacote ENC com as novas tabelas importadas. Temos de usar scripts Python em vez de consultas e fun\u00e7\u00f5es SQL porque a vers\u00e3o SQL do SQLite, usada para geopacotes, \u00e9 muito limitada e n\u00e3o permite a cria\u00e7\u00e3o de procedimentos ou fun\u00e7\u00f5es.<\/p>\n\n\n\n Se voc\u00ea tiver um \u00fanico arquivo S57 para importar, o processo \u00e9 o seguinte:<\/p>\n\n\n\n Se voc\u00ea tiver um lote de arquivos S57 para carregar simultaneamente, o processo \u00e9 o seguinte:<\/p>\n\n\n\n As sondas batim\u00e9tricas apresentam v\u00e1rios problemas na convers\u00e3o de formatos. O primeiro \u00e9 que o valor da profundidade n\u00e3o est\u00e1 contido em um atributo, mas como Z na geometria (XYZ). O segundo \u00e9 que as sondas n\u00e3o s\u00e3o do tipo Ponto, mas do tipo Multiponto. Para obter o valor das sondas diretamente em um campo de atributo, \u00e9 necess\u00e1rio adicionar dois par\u00e2metros \u00e0 linha de comando do ogr2ogr:<\/p>\n\n\n\n -oo SPLIT_MULTIPOINT=ON -oo ADD_SOUNDG_DEPTH=ON<\/em><\/strong><\/p>\n\n\n\n <\/div><\/div>\n\n\n\n Alguns objetos podem ser agrupados por um objeto pai n\u00e3o geogr\u00e1fico. Por exemplo, os setores de inc\u00eandio s\u00e3o codificados com um setor por registro na tabela \u00ab\u00a0LIGHTS\u00a0\u00bb. Os diferentes setores do mesmo inc\u00eandio n\u00e3o t\u00eam nenhum atributo espec\u00edfico que possa indicar que eles correspondem \u00e0 mesma entidade. Essa informa\u00e7\u00e3o est\u00e1 contida em um registro \u00ab\u00a0pai\u00a0\u00bb. Para recuperar esse identificador como um atributo da tabela LIGTHS, \u00e9 necess\u00e1rio adicionar uma op\u00e7\u00e3o \u00e0 linha de comando:<\/p>\n\n\n\n -oo RETURN_LINKAGES=O<\/strong><\/p>\n\n\n\n <\/div><\/div>\n\n\n\n Essa op\u00e7\u00e3o cria um atributo name_rcid que \u00e9 comum a todos os setores no mesmo sem\u00e1foro, mas tamb\u00e9m cria uma s\u00e9rie de campos (name_rcnm,ORNT,USAG,MASK).<\/p>\n\n\n\n Al\u00e9m dos tipos de campo tradicionais (inteiro, real, texto), o formato S57 usa um tipo especial: listas de cadeias de caracteres ou inteiros. Esses tipos de campo s\u00e3o encontrados no PostgreSQL, mas n\u00e3o nos shapefiles e no geopackage.<\/p>\n\n\n\n Na aus\u00eancia de op\u00e7\u00f5es espec\u00edficas, haver\u00e1 mensagens de aviso indicando que o formato n\u00e3o \u00e9 compat\u00edvel e o conte\u00fado desses campos ser\u00e1 ignorado. Para evitar a perda do conte\u00fado, voc\u00ea precisa adicionar a extens\u00e3o :<\/p>\n\n\n\n -mapFieldType StringList=String,IntegerList=String<\/em><\/strong><\/p>\n\n\n\n <\/div><\/div>\n\n\n\n Essa op\u00e7\u00e3o usa uma lista {\u2026,\u2026} e gera uma string de texto com um primeiro n\u00famero indicando o n\u00famero de elementos na lista, dois pontos e, em seguida, os elementos da lista separados por v\u00edrgulas:<\/p>\n\n\n\n {3} -> ‘(1:3)’<\/p>\n\n\n\n {3,1} -> ‘(2:3,1)’<\/p>\n\n\n\n Isso inevitavelmente complica um pouco o processamento desses campos, mas as informa\u00e7\u00f5es n\u00e3o s\u00e3o perdidas. Por outro lado, o driver GDAL apresenta um problema para gerenciar o formato String resultante.<\/p>\n\n\n\n Em primeiro lugar, n\u00e3o h\u00e1 problema com os atributos dos objetos S57. O campo na tabela de geopacotes \u00e9 do tipo \u00ab\u00a0Text\u00a0\u00bb, que permite a cria\u00e7\u00e3o de uma cadeia de caracteres de qualquer comprimento. Em seguida, voc\u00ea ver\u00e1 mensagens do 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 A solu\u00e7\u00e3o n\u00e3o \u00e9 complicada, mas envolve uma boa quantidade de trabalho. Comece com a premissa de que, para a exibi\u00e7\u00e3o cl\u00e1ssica de cart\u00f5es ENC, isso n\u00e3o representa um problema e voc\u00ea pode simplesmente ignorar essas mensagens de aviso.<\/p>\n\n\n\n Entretanto, se voc\u00ea precisar de um processamento mais avan\u00e7ado que inclua essas informa\u00e7\u00f5es, essa \u00e9 a \u00fanica maneira que encontrei at\u00e9 o momento:<\/p>\n\n\n\n <\/div><\/div>\n\n\n\n Vamos criar um geopackage vazio que receber\u00e1 as tabelas importadas e, em seguida, atualiza\u00e7\u00f5es sucessivas quando outros arquivos S57 forem carregados.<\/p>\n\n\n\n No painel QGis Explorer, clique com o bot\u00e3o direito do mouse em Geopackages <\/em><\/strong>e selecione Create a database<\/em><\/strong> (Criar um banco de dados)<\/p>\n\n\n\n Na janela que se abre, digite o nome do geopacote, deixe a op\u00e7\u00e3o padr\u00e3o para o nome da tabela e selecione No geometry<\/em><\/strong> (Sem geometria) em Geometry type<\/em><\/strong> (Tipo de geometria).<\/p>\n\n\n\n Haga clic en Aceptar 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, abra el gestor de bases de datos de QGis desde el men\u00fa principal, sit\u00faese sobre la tabla ENC del geopackage ENC.gpkg y haga clic con el bot\u00f3n derecho del rat\u00f3n para seleccionar Borrar\u2026<\/p>\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><\/div>\n\n\n\n Para executar essas linhas de comando, basta abrir a janela do shell do OSGeo4W:<\/p>\n\n\n\n A janela Shell \u00e9 aberta<\/p>\n\n\n\n Digite as linhas de comando, tomando cuidado para inserir o caminho e o nome do arquivo com a extens\u00e3o .000.<\/p>\n\n\n\n O resultado ser\u00e1 uma s\u00e9rie de tabelas criadas em cada geopackage de importa\u00e7\u00e3o. Entretanto, algumas tabelas ficar\u00e3o completamente vazias. Se o tipo de geometria for poss\u00edvel para uma classe de objeto, a tabela ser\u00e1 criada, mas se n\u00e3o houver ocorr\u00eancias no arquivo S57 processado, a tabela n\u00e3o ter\u00e1 registros.<\/p>\n\n\n\n Usaremos o console QGis Python para todas as opera\u00e7\u00f5es de gerenciamento do banco de dados.<\/p>\n\n\n\n Para abri-lo, acesse o menu Extens\u00f5es -> Console Python. Os pain\u00e9is do console Python ser\u00e3o abertos:<\/p>\n\n\n\n Para executar um script Python, carregue o c\u00f3digo no painel Editor. Voc\u00ea pode copiar e colar ou usar o \u00edcone Diret\u00f3rio para apontar para o arquivo .py. A indenta\u00e7\u00e3o \u00e9 crucial para scripts Python, portanto, escolha a segunda op\u00e7\u00e3o.<\/p>\n\n\n\n Como vimos anteriormente, ao executar comandos ogr2ogr, s\u00e3o criadas tabelas vazias para camadas em que n\u00e3o h\u00e1 geometria solicitada na linha de comando. Portanto, removeremos essas tabelas vazias e aproveitaremos o script Python para adicionar dois atributos a todas as tabelas: o nome do mapa ENC e sua escala.<\/p>\n\n\n\n No que diz respeito ao nome, pode ser \u00fatil na manuten\u00e7\u00e3o futura do banco de dados poder selecionar as entidades correspondentes a um arquivo de origem S57. Essa informa\u00e7\u00e3o n\u00e3o est\u00e1 contida em nenhum atributo do arquivo S57. Portanto, temos de inseri-la manualmente.<\/p>\n\n\n\n O mesmo se aplica \u00e0 escala do mapa. Embora um atributo seja fornecido no S57 para inserir a escala do mapa \u00ab\u00a0em papel\u00a0\u00bb, o atributo CSCALE na tabela M_CSCL, \u00e9 muito raro que essa tabela e seu atributo estejam presentes nos arquivos S57.<\/p>\n\n\n\n Embora o nome do mapa seja de uso relativo, o mesmo n\u00e3o pode ser dito da escala, pois esse \u00e9 um crit\u00e9rio essencial na busca por duplicatas. De fato, quando entidades de v\u00e1rios arquivos S57 s\u00e3o misturadas, os dados de diferentes escalas coexistir\u00e3o de forma mais ou menos feliz.<\/p>\n\n\n\n O c\u00f3digo Python a seguir exclui tabelas vazias de um esquema de importa\u00e7\u00e3o e adiciona dois atributos (enc_chart e scale) a todas as tabelas:<\/p>\n\n\n\n <\/div><\/div>\n\n\n\n Tenha o cuidado de respeitar a indenta\u00e7\u00e3o do c\u00f3digo Python. Para simplificar a tarefa, baixe o arquivo .py diretamente deste link: delete_empty_tables_in_geopackages.py<\/a>.<\/p>\n\n\n\n Para execut\u00e1-lo, certifique-se de alterar os par\u00e2metros de entrada do script:<\/p>\n\n\n\n Voc\u00ea s\u00f3 precisar\u00e1 inserir os geopackage_paths na primeira vez que executar o script. Por outro lado, o nome do mapa ENC e a escala precisar\u00e3o ser alterados quando cada novo arquivo S57 for importado. As linhas de comando nos par\u00e1grafos anteriores pressup\u00f5em que voc\u00ea tenha um \u00fanico arquivo S57 para carregar. Mas voc\u00ea pode carregar v\u00e1rios arquivos simultaneamente. Ao usar as op\u00e7\u00f5es -append -update, voc\u00ea pode processar todos os arquivos .000 em um diret\u00f3rio em um loop. Em seguida, voc\u00ea ter\u00e1 todos os arquivos carregados nos geopacotes de importa\u00e7\u00e3o.<\/p>\n\n\n\n Diferentemente do carregamento de um \u00fanico arquivo, ao carregar um lote de arquivos, precisamos incluir a adi\u00e7\u00e3o do nome do arquivo e da escala no loop de transcodifica\u00e7\u00e3o S57->geopackage. Para fazer isso, usaremos um script python que adiciona os atributos enc_chart e scale \u00e0s tabelas criadas pelo comando ogr2ogr e as preenche com o nome do arquivo correspondente e uma escala definida pelo produtor. O script tamb\u00e9m exclui todas as tabelas vazias no geopackage de importa\u00e7\u00e3o.<\/p>\n\n\n\n Agradecemos a Janez AVSEC pela ideia de usar a tabela DSID para recuperar a escala de compila\u00e7\u00e3o de dados.<\/p>\n\n\n\n Uma atualiza\u00e7\u00e3o sobre esse procedimento est\u00e1 dispon\u00edvel no artigo aqui.<\/a><\/p>\n\n\n\n <\/div><\/div>\n\n\n\n A chamada para cada arquivo .bat tem o seguinte formato:<\/p>\n\n\n\n .\\fichier.bat \u00ab\u00a0diret\u00f3rio que cont\u00e9m os arquivos S57 a serem carregados\u00a0\u00bb \u00ab\u00a0nome do geopacote de importa\u00e7\u00e3o\u00a0\u00bb <\/em><\/strong><\/p>\n\n\n\n Por exemplo:<\/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 <\/p>\n\n\n\n Os arquivos .bat para essa tarefa s\u00e3o os seguintes:<\/p>\n\n\n\n Para importar tabelas de tipos de pontos:<\/p>\n\n\n\n <\/div><\/div>\n\n\n\n Voc\u00ea perceber\u00e1 que, al\u00e9m do comando ogr2ogr para importar camadas de pontos, importamos duas tabelas sem geometria: DSID e C_AGGR. A tabela DSID cont\u00e9m dois atributos com o nome do arquivo S57 e a escala de compila\u00e7\u00e3o de dados. Essa tabela ser\u00e1 usada pelo script Python (update_geopackage_dsid.py) para inserir o arquivo S57 original e a escala de compila\u00e7\u00e3o em cada tabela do geopackage.<\/p>\n\n\n\n Para importar tabelas do tipo linha:<\/p>\n\n\n\n <\/div><\/div>\n\n\n\n Para importar tabelas de pol\u00edgonos:<\/p>\n\n\n\n <\/div><\/div>\n\n\n\n Os arquivos .bat cont\u00eam o comando ogr2ogr e uma linha para execu\u00e7\u00e3o de um script Python, update_geopackage_dsid.py. Certifique-se de alterar o caminho desse script no arquivo .bat para o caminho de sua escolha. Esse script adiciona os dois atributos enc_chart e scale a todas as tabelas do geopackage.<\/p>\n\n\n\n O c\u00f3digo python usado \u00e9 o seguinte:<\/p>\n\n\n\n import sys def get_default_scale(geopackage_path, default_enc): def delete_empty_tables_in_geopackages_dsid(geopackage_path, default_enc): if name<\/strong> == \u00ab\u00a0main<\/strong>\u00ab\u00a0: <\/div><\/div>\n\n\n\n <\/p>\n\n\n\n Voc\u00ea poderia criar um \u00fanico arquivo .bat para todos os tr\u00eas processos, mas precisaria verificar os resultados para cada tipo de geometria.<\/p>\n\n\n\n Voc\u00ea pode baixar os arquivos .bat e o c\u00f3digo Python clicando aqui.<\/a><\/p>\n\n\n\n <\/p>\n\n\n\n Para o primeiro carregamento do S57, o script Python a seguir cria todas as tabelas presentes 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 A clonagem \u00e9 realizada usando o seguinte c\u00f3digo:<\/p>\n\n\n\n <\/div><\/div>\n\n\n\n <\/p>\n\n\n\n Voc\u00ea pode fazer o download do arquivo .py aqui: clone_or_append_tables_with_prefix.py.<\/a><\/p>\n\n\n\n Quando dois mapas ENC s\u00e3o sobrepostos, as entidades presentes na zona de sobreposi\u00e7\u00e3o ser\u00e3o inseridas tantas vezes quantas forem as sobreposi\u00e7\u00f5es. Para evitar duplicatas e mais, o banco de dados deve ser limpo ap\u00f3s cada importa\u00e7\u00e3o de um arquivo S57.<\/p>\n\n\n\n Mas h\u00e1 dois tipos diferentes de duplica\u00e7\u00e3o. Por um lado, as \u00e1reas em que dois mapas de escala semelhante se sobrep\u00f5em produzir\u00e3o duplicatas \u00ab\u00a0reais\u00a0\u00bb que precisam ser limpas. Por outro lado, a sobreposi\u00e7\u00e3o de mapas de escalas muito diferentes precisa ser tratada de forma diferente. O mapa menos detalhado (escala pequena) ter\u00e1 uma sele\u00e7\u00e3o de dados do mapa mais detalhado (escala grande). Nesse segundo caso, n\u00e3o \u00e9 aconselh\u00e1vel remover essas informa\u00e7\u00f5es.<\/p>\n\n\n\n O c\u00f3digo Python abaixo remover\u00e1 as duplicatas \u00ab\u00a0reais\u00a0\u00bb. Para outras duplicatas, recomendamos o uso de filtros QGis. Dependendo de suas necessidades, voc\u00ea pode usar o atributo \u00ab\u00a0scale\u00a0\u00bb (escala) que adicionamos ao excluir as tabelas vazias para definir quais dados s\u00e3o exibidos no seu mapa e quais n\u00e3o s\u00e3o.<\/p>\n\n\n\n Nesse exemplo, somente os recursos correspondentes a uma escala entre 12500 e 45000 ser\u00e3o exibidos em seu mapa.<\/p>\n\n\n\n Essa opera\u00e7\u00e3o corresponde \u00e0 filtragem de uma camada, mas\u2026 o que voc\u00ea faz quando tem 130 camadas exibidas? Para evitar que voc\u00ea tenha de executar essa opera\u00e7\u00e3o 130 vezes, aqui est\u00e1 um c\u00f3digo Python que permite filtrar todas as camadas do projeto:<\/p>\n\n\n\n <\/div><\/div>\n\n\n\n E outro para remover o filtro de todas as camadas exibidas:<\/p>\n\n\n\n <\/div><\/div>\n\n\n\n <\/p>\n\n\n\nTipos de geometria<\/h2>\n\n\n\n
Attributes: CONDTN OBJNAM NOBJNM STATUS INFORM NINFOM<\/strong><\/p>\n\n\n\n\n
\n
Processamento de sondagens batim\u00e9tricas<\/h2>\n\n\n\n
<\/div>
Recupera\u00e7\u00e3o de identificadores \u00ab\u00a0pai\u00a0\u00bb<\/h2>\n\n\n\n
<\/div>
Processamento do tipo \u00ab\u00a0Lista\u00a0\u00bb<\/h2>\n\n\n\n
<\/div>
No entanto, para \u00ab\u00a0parent identifiers\u00a0\u00bb, a defini\u00e7\u00e3o usada pelo driver \u00e9 \u00ab\u00a0Text()\u00a0\u00bb, ou seja, uma cadeia de caracteres com um comprimento m\u00e1ximo. Esse comprimento \u00e9 determinado a priori pelo conte\u00fado do primeiro registro processado. Em outras palavras, \u00e9 apenas por acaso que voc\u00ea tem um campo com o comprimento necess\u00e1rio para incluir todos os valores no arquivo.<\/p>\n\n\n\n<\/div>
\n
Fluxo de trabalho para fazer upload de um \u00fanico arquivo S57<\/h1>\n\n\n\n
Cria\u00e7\u00e3o do geopacote ENC<\/h2>\n\n\n\n
![\"\"](\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2024\/01\/creation-geopackage.png?resize=315%2C335&ssl=1\")
<\/a><\/figure>\n\n\n\n![\"\"](\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2024\/01\/creation-geopackage-ENC.png?resize=550%2C717&ssl=1\")
<\/a><\/figure>\n\n\n\ncomandos ogr2ogr para crear geopaquetes de importaci\u00f3n<\/h2>\n\n\n\n
<\/div>
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![\"Osgeo4W](\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2023\/12\/shell_osgeo4w.png?resize=640%2C734&ssl=1\")
<\/a><\/figure>\n\n\n\n![\"Osgeo4W](\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2023\/12\/osgeo4w_shell.jpg?resize=640%2C366&ssl=1\")
<\/a><\/figure>\n\n\n\nUso de scripts Python para gerenciar o banco de dados<\/h2>\n\n\n\n
![\"\"](\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2024\/01\/console-python-qgis-1.png?resize=640%2C400&ssl=1\")
<\/a><\/figure>\n\n\n\nRemova as tabelas vazias e adicione o nome e a escala do mapa<\/h2>\n\n\n\n
<\/div>
from 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![\"\"](\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2024\/01\/lignes_a_editer_python.png?resize=640%2C61&ssl=1\")
<\/a><\/figure>\n\n\n\n
No final dessa execu\u00e7\u00e3o, voc\u00ea ter\u00e1 os novos atributos em cada tabela<\/p>\n\n\n\n![\"\"](\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2024\/01\/Capture-decran-2024-01-28-144751.png?resize=640%2C236&ssl=1\")
<\/a><\/figure>\n\n\n\nFluxo de trabalho para carregar v\u00e1rios arquivos S57 ao mesmo tempo<\/h1>\n\n\n\n
<\/div>
<\/div>
@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>
@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>
@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>
from osgeo import ogr<\/p>\n\n\n\n
# 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\nif 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\n
# 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\n
if len(sys.argv) != 3:
print(\u00ab\u00a0Usage: python script.py geopackage_path default_enc\u00a0\u00bb)
sys.exit(1)<\/p>\n\n\n\ngeopackage_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\nContinua, independentemente do n\u00famero de arquivos carregados<\/h1>\n\n\n\n
Clonagem e\/ou adi\u00e7\u00e3o de tabelas importadas ao geopacote ENC<\/h2>\n\n\n\n
<\/div>
from 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\nGerenciamento de duplicatas<\/h2>\n\n\n\n
![\"\"](\"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2024\/01\/filtre_scale.png?resize=640%2C453&ssl=1\")
<\/a><\/figure>\n\n\n\n<\/div>
#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>
# 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