QGIS 3.20 | LUCCHESE, L. V.

Resolvido: Python Console do QGIS 3.20 mostra um erro ao executar a ferramenta Clip Vector by Extent

Mon, 22 Nov 2021 12:43:11 -0300

Colegas pesquisadores e entusiastas de SIG de código-livre,

Bem-vindos ao meu blog!

Gostaria de começar com um aviso - posso ser uma pesquisadora desta área, mas isso não significa que tudo o que faço ou escrevo aqui funcionará para você, em suas próprias configurações de desktop e versões de packages. Não tenho nenhuma responsabilidade se você perder dados ou estragar sua instalação. Eu também não autorizo nenhuma cópia do meu conteúdo.

Hoje, estou escrevendo sobre um erro que obtive ao usar o Console Python do QGIS 3.20. Eu estava tentando recortar uma camada vetorial preenchida com pontos que tinha acabado de criar. Estava cortando ela com base em uma extensão do tipo $xmin, xmax, ymin, ymax$, por meio da ferramenta Clip Vector by Extent, do GDAL.

O erro

Primeiro, executei a seguinte linha, para definir os parâmetros da ferramenta:

paramgdal= {'INPUT': points_layer, 'EXTENT':extent_final, 'OUTPUT':output_path}

As variáveis “extent_final” e “output_path” foram previamente setadas. “points_layer” é a camada vetorial (no caso, um shapefile) em que o universo de pontos dos quais quero recortar uma parte estão armazenados. Ele foi carregado anteriormente, usando o seguinte comando:

points_layer = QgsVectorLayer(path_to_file, 'points_loaded', 'ogr')

Até aqui, tudo bem.

Mas, depois de definir os parâmetros, executei:

processing.run("gdal:clipvectorbyextent", paramgdal)

e aí, este erro apareceu:

ERROR 1: Attempt to write non-multipoint (POINT) geometry to multipoint shapefile.
ERROR 1: Unable to write feature 0 from layer points_layer.
ERROR 1: Terminating translation prematurely after failed translation of layer points_layer (use -skipfailures to skip errors)

E o layer resultante não foi carregado (porque não foi criado corretamente).

Interpretação do problema

Relendo a mensagem de erro, percebi que o GDAL, por algum motivo, tenta salvar a camada gerada como um shapefile do tipo Multipoint, ao invés de, simplesmente, um shapefile do tipo Point. Por esse motivo, ele não consegue gravar nenhum dos pontos que estavam no layer original dentro do novo shapefile, porque os tipos não correspondem.

Solução

Verificando a documentação do QGIS, podemos notar que existem quatro possíveis parâmetros para executar gdal:clipvectorbyextent. Um deles, “OPTIONS”, é opcional, e se refere às opções de criação do GDAL. Como o algoritmo Clip Vector by Extent é baseado em um utilitário do GDAL chamado ogr2ogr, as opções disponíveis para o algoritmo original também podem ser utilizadas para construir o novo layer, no caso, no Python Console.

Uma das opções é particularmente útil, neste cenário. A opção “nlt” do ogr2ogr define a geometria do arquivo gerado, o que é exatamente o que estávamos procurando. No caso mostrado aqui, simplesmente adicionei

'OPTIONS': '-nlt point'

ao dicionário Python chamado de $gdalparams$. Isso diz ao GDAL que o vetor a ser criado é um vetor do tipo Points. E isso resolveu meu problema completamente! O fragmento de código mostrado no início desta postagem foi adaptado por mim como:

paramgdal= {'INPUT': points_layer, 'EXTENT':extent_final, 'OPTIONS': '-nlt point', 'OUTPUT':output_path}
processing.run("gdal:clipvectorbyextent", paramgdal)

E quando executei o fragmento acima, a camada resultante foi gerada corretamente, sem erros ou warnings.

Ou seja, uma solução simples!

Extras

Por que você simplesmente não usa o GDAL direto por linha de comando, ou clica na ferramenta e a usa diretamente na interface gráfica? Por que você quer usar o Console Python do QGIS?

Porque, às vezes, o console Python é a maneira mais prática de realizar um cálculo. Este trecho de código mostrado aqui é uma pequena parte de um código muito maior. Por esse motivo, não vale a pena mudar toda a interface ou linguagem na qual estamos desenvolvendo o código.

Carreguei meu arquivo vetorial no Python Console, mas ele não apareceu na lista de camadas do QGIS! O que aconteceu?

Isto não é um problema. O trecho de código usado para carregar a camada vetorial não lista a camada no QGIS, nem a mostra a visualização da camada no mapa. Para fazer isso, você deve executar

QgsProject.instance().addMapLayers([points_layer])

Não sei o que devo importar antes de usar a ferramenta “gdal:clipvectorbyextent” no Console Python do QGIS.

Você deve importar os algoritmos do QGIS e a ferramenta de processamento. Veja o fragmento de código abaixo:

from qgis.core import *
import processing

Amostragem de rasters em pontos pré-determinados no QGIS

Mon, 04 Oct 2021 10:10:08 -0300

Colegas pesquisadores e entusiastas de SIG de código-livre,

Bem-vindos ao meu blog!

Hoje, vou escrever sobre como fazer amostragem de valores de rasters em pontos pré-determinados. Isso pode ser útil para uma variedade de aplicações, em especial, análises estatísticas.

Neste post, vou mostrar como realizar este processo usando duas ferramentas diferentes, uma delas é nativa do QGIS e a outra é um plugin disponível no repositório de plugins. Em seguida, vou comentar sobre como amostrar um raster usando interpolação do tipo cúbica, spline cúbica, bilinear, entre outras.

Para o exemplo mostrado aqui, estou usando o mesmo conjunto de dados que usei para a palestra no QGIS Open Day Brasil Setembro/2021, disponível em https://github.com/luisalucchese/QGISday_092021.

Os rasters são baseados em um Modelo Digital de Elevação (MDE) “imaginário”, criado usando o código dem-utils, desenvolvido por ArMoraer e disponível no GitHub. As projeções foram atribuídas às camadas apenas para fins de exemplo, pois o MDE não retrata uma parte real do globo.

Abrindo os rasters no QGIS 3.20 “Odense”, vejo isto:

Então, para este exemplo, eu crio uma camada vetorial de pontos e adiciono três pontos a ela.

Existem duas maneiras principais de extrair o valor dos rasters nos pontos fornecidos.

Usando a ferramenta nativa “Sample raster values”

Encontre a ferramenta “Sample raster values” ou “native: rastersampling” na sua caixa de ferramentas de processamento.

Usando este algoritmo, você deve amostrar cada camada raster individualmente.

Mudei o prefixo da coluna de saída para “DEM_” porque tenho vários rasters para amostrar e quero que as colunas resultantes fiquem com nomes diferentes.

Ao clicar em Run, o algoritmo foi executado sem problemas.

Abrindo a tabela de atributos da minha nova camada vetorial, encontro:

A coluna recebeu o sufixo “1”, provavelmente por ser a primeira e única banda do raster.

Dando zoom no raster na região em que cada ponto está localizado e colhendo os valores das células onde os pontos se localizam, encontro os valores:

Os valores das células do raster em que cada ponto está contido são iguais aos extraídos pelo algoritmo. Isso significa que o valor é simplesmente extraído do raster e nenhuma interpolação é realizada, o que é o mesmo que dizer que a ferramenta executou uma interpolação do tipo “Vizinho Mais Próximo”.

Para extrair os valores de todos os quatro rasters, execute a ferramenta mais três vezes. Em cada rodada, selecione o vetor que foi gerado na rodada anterior, mas escolha um raster diferente dos anteriores e ajuste o nome da coluna de acordo. Aqui, quando o processo foi concluído, abri a tabela de atributos da camada de pontos resultante:

Esse processo se torna tedioso se tiver muitas camadas raster para extrair as informações. E o processamento em lote disponível na própria interface do QGIS não é uma opção neste caso, já que você precisa da saída de uma rodada para usar como entrada na outra. Porém, uma opção para lidar com isso seria construir um raster multibanda a partir dos rasters que você tem em layers separados e, em seguida, executar o algoritmo neste raster multibanda.

Usando o plugin “Point Sampling Tool”

O plugin “Point Sampling Tool” pode ser encontrado no gerenciador de plugin do QGIS. Eu estou usando a versão 0.5.3 deste plugin.

Quando abro o plug-in, vejo a seguinte tela:

Selecione mais de um raster segurando Ctrl enquanto clica nos nomes dos layers.

Neste plugin, não é possível salvar um arquivo como temporário.

Seleciono a pasta do arquivo a ser gerado no botão “Browse” e, em seguida, clico em OK.

O plugin me mostrou este aviso “All layers must have the same coordinate refere system. The <layer name> layer seems to have different CRS id (0) than the point layer (32722) […]”

Nesse caso, estou usando MDEs criados apenas para fins de exemplo e tenho certeza de que todos eles “descendem” do mesmo raster, então cliquei em “Yes”. No entanto, você deve analisar caso a caso. Normalmente, é melhor reprojetar a camada raster ou a camada vetorial, para que estejam na mesma projeção.

Abrindo a tabela de atributos do vetor resultante:

Os valores na tabela de atributos ficaram semelhantes aos selecionados pela ferramenta nativa do QGIS, porém com menor número de algarismos significativos.

Eles também são amostrados sem interpolação, ou seja, pelo método “Vizinho mais próximo” (Nearest Neighbor).

Os nomes das colunas são baseados nos nomes com que as camadas estão em seu projeto do QGIS.

Extras

- Qual dos métodos devo usar?

Ambas as ferramentas funcionam muito bem. Na minha opinião, depende do número de rasters de que você quer extrair informações. Se forem apenas um ou dois rasters, eu usaria a ferramenta nativa, mas se forem mais do que isso, acho que vale mais a pena usar o plugin.

- Você mencionou que ambas as ferramentas usam a interpolação do tipo Vizinho Mais Próximo. Como faço para usar outro tipo de interpolação?

Bom, por meio das ferramentas em si, não é possível. Mas tem uma solução alternativa para isso. Usando a ferramenta “Warp (reproject)” (“gdal: warpreproject”), você pode fazer resample dos seus rasters usando um tamanho de pixel menor e um tipo diferente de interpolação e, em seguida, executar qualquer uma das ferramentas de amostragem que eu mostrei aqui.

Exemplo:

Raster original:

Escolhendo a interpolação “Cubic” e escolhendo um tamanho de pixel menor para o raster de saída, este é o raster resultante:

Agora, se eu extrair o valor deste ponto, ele será 61,9293, em vez de 61,7967. Neste caso, em que os pixels originais não eram muito grandes nem diferiam muito em valor entre si, a diferença no valor extraído foi da ordem de $10^{-1}$. Mas, dependendo do caso, pode ser bem maior. É sempre uma boa ideia examinar o seu raster e o uso pretendido dos valores extraídos, analisando a situação para ver se você deve reamostrar seu raster antes de iniciar a extração dos valores nos pontos.

Testando algumas das minhas novidades favoritas do QGIS 3.20 Odense

Mon, 21 Jun 2021 16:57:19 -0300

Colegas pesquisadores e entusiastas de SIG de código-livre,

Bem-vindos ao meu blog!

Hoje foi lançado o QGIS 3.20 “Odense”!

Na postagem de hoje, testarei alguns dos meus recursos favoritos da nova versão.

Para ver todos os novos recursos e os bugs corrigidos, verifique o Log de alterações do QGIS 3.20.

Abrindo a tabela de atributos apenas com os vetores selecionados

Especialmente útil ao lidar com shapefiles com muitos polígonos, linhas ou pontos.

Vou testar isso em um shapefile com os países do mundo.

Selecione os polígonos

Eu estou selecionando todos os países dentro desse limite de latitudes.

Resultando nesta seleção:

Abra a tabela de atributos

Clique em Open Attribute Table (Selected Features)

O que abre uma tabela de atributos apenas com os países selecionados:

Para parar de ver apenas os vetores selecionados, clique em “Show Selected Features” e em seguida, clique em “Show All Features”.

Digitalização à mão livre

O QGIS agora suporta o desenho à mão livre de linhas e polígonos.

Na barra de ferramentas de digitalização avançada, selecione “Stream Digitizing”:

Após clicar no primeiro ponto, mova o mouse pela tela para desenhar a linha ou polígono.

Cursor se move automaticamente para os pontos de início e fim de linha.

Como nos softwares de Computer Aided Design (CAD), agora no QGIS você também pode escolher os pontos inicial e final de uma linha ao desenhar uma nova linha ou polígono.

Primeiro, desenhei uma linha (em verde). Depois, cliquei para começar a desenhar outra linha.

Certifique-se de que a barra de ferramentas Snapping Toolbar esteja ativada em View, Toolbars, Snapping Toolbar. Selecione a opção “Line Endpoints”.

Após selecionar esta opção, quando o cursor estiver próximo ao ponto final de uma linha, ele deverá mostrar um marcador especial sobre o ponto.

Apenas clique quando o marcador aparecer, para desenhar a linha.

Linhas de espessura variável e linhas em gradiente

~~Este recurso está listado no changelog do QGIS 3.20. Infelizmente, não consegui fazer esse tipo de simbologia funcionar no meu QGIS 3.20.0 Odense.~~

~~Vou atualizar este post caso eu consiga fazer isso funcionar.~~

Atualizado! 22/06/21

Este recurso me parece especialmente útil para representações ligadas à hidrologia. Por exemplo, configurar os rios para que tenham linhas mais finas à montante e mais grossas à jusante.

Primeiramente, abra o painel de simbologia da sua linha

Clique em Simple Line

Escolha o tipo de símbolo como Interpolated Line

Deve abrir este painel:

Troque o tipo de largura da linha

Troque de Fixed width para Varying width. Isto deve aparecer:

Se você clicar na flechinha do lado de Start value, o QGIS vai oferecer “id”. Não clique!

Na verdade, nesse campo devemos escrever um número, diferente do número do campo “end value”. Isso que vai fazer que a espessura seja variável.

Coloco 0.1 em Start Value e 10 em End Value. Clico em Apply.

Pronto! Espessura variável.

Agora, cor variável

No seletor, troque de Single color para Varying color. Este painel irá abrir:

Novamente, não selecione o id em Start Value ou em End Value. Digite números.

Digito 0.1 como Start Value e 10 como End Value, e clico em Apply.

A esse ponto, a linha está assim:

Podemos trocar a escala de cores em “Color ramp”. Eu troquei para a denominada “Blues”, que tem tons de azul, e servirá bem para representações de hidrografia.

Pronto!

Estas são as minhas novidades favoritas do QGIS 3.20, até agora!

Se você trabalha com rótulos em vetores, muitos recursos novos sobre esse assunto foram implementados, então sugiro que você dê uma olhada no log completo.

Vale a pena instalar o QGIS 3.20, sendo que o QGIS 3.22 LTS está programado para sair em outubro?

Eu não tive erros graves no QGIS 3.18 e ele inclui a implementação de legendas de raster contínuas, que eu uso com grande frequência. Eu preferi não esperar até outubro para ter legendas de raster contínuas no Layout de impressão, sabendo que esse recurso já existia. Há mais de um ano, tenho instalado versões não-LTS do QGIS e, até agora, não tive muitos problemas. Dito isso, algumas pessoas têm experiências diferentes da minha, e instalar uma versão não-LTS é considerado algo arriscado por alguns profissionais.