Vamos começar a criar um novo scripts em Python para ler uma coluna da tabela de atributo e realizar alguns cálculos estatísticos.
Abra
o Qgis e Clique no ícone do Python como mostra a figura abaixo.
Após clique no ícone do editor do Python e na área marcada
em vermelho digite as linhas de código que será mencionada abaixo.
Digite o código da imagem abaixo.
É importado a biblioteca do numpy para auxiliar na
realização dos cálculos.
A variável “layer” recebe a camada ativa na tela do qgis.
A variável “lista” recebe uma lista vazia que será
preenchida com os dados de uma dada coluna.
O loop for vai ler os dados da coluna desejada e irá inserir
os dados na “lista”.
Para realizar cada cálculo é utilizado a função referente,
sendo apontado a lista como a fonte dos dados a serem calculados e as variáveis
recebem o resultado desses cálculos.
A variável “output_file” recebe a abertura de um arquivo txt
no modo escrita (“w”). Uma observação importante é que se caso o arquivo não
existir ele será criado automaticamente.
As linhas seguintes
são para transformar os resultados numéricos em string que será concatenado os
nomes dos resultados.
A função print irá retornar os resultados dos dados.
As duas últimas linhas do script são para escrever os
resultados no arquivo de texto e fecha-lo.
*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 dos Loops for.
Coloque os espaços corretos como na figura utilizando a tecla Tab.
Observe na imagem abaixo uma parte do código dentro do editor.
Após esse procedimento clique no ícone rodar script e visualize os resultados na tela.
Também é gerado um arquivo de texto na pasta informada na
variável “output_file” com o mesmo resultado mostrado na tela.
*Caso apresente algum erro revise o código.
Abaixo se encontra o vídeo de demonstração do script sendo executado.
Abaixo se encontra o vídeo de demonstração do script sendo executado.
Até a próxima!
Quer aprender como abrir shapefile com PyQGIS? Clique aqui