Interpolação com Python no QGIS





Nesse tutorial vamos criar um script para interpolar dados de um shapefile de pontos.

A imagem abaixo mostra o editor do python aberto no QGIS 3.2.



Com o editor do python aberto digite o seguinte script.



A primeira variável recebe o caminho do shapefiles de ponto que será a base para realizar a interpolação. Uma observação importante é que depois do “.shp” tem uma pequena configuração que aponta se existe ou não a dimensão “Z” que nesse exemplo não existe, o número 1 aponta a coluna onde estão presente os dados de peso e o último 0 aponta o tipo que nesse caso é de pontos.

A variável “saída” vai receber o caminho onde será salvo o arquivo gerado pela interpolação.

A variável “caminhoArea” vai receber o caminho da área limite da interpolação.

A variável seguinte recebe o nome da camada de ponto.

A função “abrir” é a junção das três próximas linhas e é responsável por abrir um determinado shapefile quando chamada abaixo e é necessário passar dois parâmetros, sendo eles o caminho dos shapefile e o nome da camada.

As variáveis “camadaAtiva” e “camadaAtiva_01” recebem as camadas ativas, que são a camada de ponto e a área limite da interpolação.

A função processing vai ser responsável por realizar a interpolação e nele você terá que passar os parâmetros do tipo de interpolação, qual dado que é a base da interpolação, o coeficiente, a camada limite, a camada de saída e outros parâmetros.

*Obs: O nome do usuário no script está em "*****" e você deverá inserir o caminho do arquivo do seu computador.

*Obs: Todas as “ASPAS” fazem parte do código e se caso não coloca-las o Python irá retornar um Erro. Outra observação importante é a correta indentação dentro da função. Coloque os espaços corretos como na figura utilizando a tecla Tab.

A imagem a seguir mostra uma parte do script dentro do editor.




Após esse procedimento clique no ícone rodar script para ver o resultado na tela.



Na imagem a seguir foi alterado a rampa de cor para uma melhor visualização.



*Caso apresente algum erro revise o código.

Abaixo se encontra o vídeo de demonstração do script sendo executado.



Até a próxima!


Será que você precisa de um mapa? Entre aqui se você precisa de uma consultoria.

Quer aprender mais sobre Python no QGIS? Clique aqui