iDea Plus Geo

Geo ideias em formato aberto

Recorte, reclassificação e vetorização de raster no gvSIG

August 28, 2019 — Eliazer Kosciuk

Buenas, pessoal! Fazendo uma pequena parada antes de estudarmos o próximo índice radiométrico da série, vamos acatar a sugestão do nosso amigo Gilberto Cugler de aproveitarmos o resultado do raster NDVI que geramos no último artigo e estudarmos algumas das ferramentas que o gvSIG disponibiliza para lidarmos com arquivos raster.

Recortando a camada raster

Para começarmos nossos trabalhos, vamos recortar do raster a nossa área de interesse, que nesse exemplo será a área do município que estamos estudando. Os arquivos para esse exercício (camada raster do NDVI e o shape do município) podem ser baixados aqui. Criamos uma Vista no gvSIG com EPSG 23722 e inserimos nela os arquivos que baixamos. Devemos estar com uma área de trabalho semelhante a essa, conforme visto anteriormente (apliquei a tabela de cores ao NDVI):

Área de trabalho inicial

Para recortarmos a camada raster, abrimos a "Caixa ferramentas" e selecionamos o geoprocesso "Recortar camada raster com camada de polígonos".

Caixa de ferramentas - Reclassificação

Em "Camada raster > Camada a cortar" selecionamos a camada raster que queremos cortar (NDVI, nesse caso), em "Camada vetorial > Polígonos" selecionamos a camada contendo o limite do município, e em seguida clicamos em "Aceitar".

Janela - Recortar camada raster com camada de polígonos

E teremos a nossa camada de estudo já inserida na área de trabalho:

Área de trabalho - imagem raster recortada

Podemos excluir as camadas que usamos para gerar essa última, e passaremos para o próximo passo.

Reclassificação de arquivos raster

Para fins de visualização, apliquei novamente a tabela de cores NDVI no raster gerado:

NDVI com tabela de cores

Observe que, muito embora estamos visualizando o raster dividido em 10 classes, cada pixel dele continua com valores entre -1,0 e 1,0, o que pode ser comprovado se utilizarmos a ferramenta de "Informação por ponto", conforme já vimos no artigo anterior.

Se quisermos calcular qual a área de vegetação apresenta um desenvolvimento vigoroso, por exemplo, precisaremos aplicar uma reclassificação à imagem raster que temos, de modo que cada pixel da imagem tenha um valor apenas para cada faixa de valores que definirmos.

Para aplicar uma reclassificação, abrimos novamente a "Caixa de ferramentas" e selecionamos o geoprocesso "Reclassificação":

Caixa de ferrmentas - Reclassificação

Na opção "Camadas raster > Camada a reclassificar" selecionamos a camada que queremos reclassificar. Em "Opções", "Método", selecionamos a opção "Min < x < Máx"; em "Tabela de correspondências", vamos clicar no botão "..." ao lado da opção "Tabela fixa 1x3".

Janela geoprocesso Reclassificação

Abrirá uma nova janela, onde poderemos estabelecer os limites que desejamos aplicar à nossa reclassificação:

Janela de entrada de dados da tabela de reclassificação Vamos aplicar a seguinte tabela para reclassificar nosso raster, lembrando sempre que os valores esperados estão entre -1,0 e 1,0: | Faixas | Interpretação | Novo valor | |--|--|--| | -1,0 a 0,0 | Plantas mortas ou objeto inanimado (solo nú, cursos de água, construções, ...) | 1 | | 0,0 a 0,33 | Plantas não saudáveis | 2 | | 0,33 a 0,66 | Plantas moderadamente saudáveis | 3 | | 0,66 a 1,0 | Plantas saudáveis | 4 |

Para isso, vamos inserindo os valores "mínimo" e "máximo" e "novo valor" nas células correspondentes da tabela (lembre de utilizar o ponto (".") como separador de decimal), clicando no botão "Adicionar linha" para adicionar as linhas que necessitamos.

Tabela de reclassificação a ser aplicada

Observe que temos a opção de "Salvar" a tabela criada, para que possamos utilizá-la em outras situações, o que pode ser feito através da opção "Abrir".

Ao terminarmos, podemos clicar em "Aceitar", que teremos a nossa tabela de reclassificação pronta para ser aplicada sobre o raster (Note que agora alterou para "Tabela fixa 4x3). Para isso, basta clicarmos na opção "Aceitar" na janela "Reclassificação", e teremos a nossa nova camada inserida na área de trabalho:

Resultado da reclassificação

Utilizando novamente a ferramenta "Informação por ponto" podemos verificar que os pixels do raster agora apresentam os valores variando entre 1 e 4, conforme definimos na nossa tabela de reclassificação.

Se quisermos , podemos definir uma nova tabela de cores para o nosso raster, refletindo os novos valores:

Reclassificação com tabela de cor aplicada

Vetorizando

Para que possamos calcular a área de cada classe que definimos, precisamos transformar o nosso raster em um vetor, através do geoprocesso "Vetorizar camada raster (polígonos)", que pode ser encontrada na "Caixa de ferramentas":

Caixa de ferramentas - Vetorizar camada raster (polígonos)

O único parâmetro que precisamos definir é a "Camada de entrada". Cuide para selecionar a camada correta.

Janela Vetorizar camada raster (polígonos)

Clicando em "Aceitar", o geoprocesso será aplicado, e a nova camada será inserida na área de trabalho. Observe que já aplicamos uma nova simbologia, para refletir as quatro classes que definimos:

Camada vetorizada, com simbologia aplicada

Concluída esta fase podemos passar para o nosso último passo:

Calculando as áreas de cada classe de vegetação

Primeiramente, vamos calcular a área de cada polígono: com a camada vetorial selecionada no ToC, abra a sua tabela de atributos (menu "Camada > Mostrar tabela de atributos" ou através do ícone correspondente na barra de ícones: Ícone Mostrar tabela de atributos). Visualizaremos a tabela de atributos da nossa camada:

Tabela de atributos da camada

Em seguida, clique no menu "Tabela > Adicionar medição > Área" (ou utilize o ícone: Ícone: Área). Abrirá uma janela com o aviso: "Serão adicionados novos campos na tabela. Deseja continuar?". Clique em "Sim". Terminado o processamento poderemos conferir que foi acrescentado o campo "Área" na tabela, contendo o valor da área de cada polígono.

Tabela de atributos com campo Área inserido

Para calcular a área de vegetação com desenvolvimento vigoroso, por exemplo, precisamos clicar sobre a coluna "Área" na tabela de atributos e, em seguida, utilizar a ferramenta de seleção por atributos, que pode ser acessada através do menu "Seleção > Seleção por atributos" ou do ícone correspondente: Ícone Seleção por atributos.

Na janela que abre, selecionaremos o campo "Resultados" na lista de "Campos", clicamos no botão "=" e, em seguida, no campo "Valores conhecidos", selecionamos o campo "4", que corresponde ao valor referente às áreas de vegetação com desenvolvimento vigoroso. Depois, clicamos no botão "Novo conjunto" e, depois, em "Aceitar".

Janela de Seleção por atributos

Na área de trabalho poderemos visualizar as áreas selecionadas:

Camada vetorizada com seleção aplicada

Agora, para calcular a área de vegetação com desenvolvimento vigoroso, precisamos clicar sobre a coluna "Área" na tabela de atributos e, em seguida, utilizarmos a ferramenta de Estatísticas, através do menu "Tabela > Estatísticas", ou através do ícone: Ícone Estatísticas.

Janela de Estatísticas da tabela - campo Área

Na janela que abre, no primeiro campo nos é informada a soma de todos os campos "Área" das geometrias selecionadas: 1.506923E8, isto é, 150.692.300 m² ou 15.069,23 ha.

Para sabermos qual a unidade de medida da área, podemos observar a barra de status do gvSIG, no canto inferior direito: gvSIG, Barra de status

Para calcular as áreas das outras classes devemos repetir o procedimento utilizado.

Tags: gvSIG, geoprocessos, geoprocessamento, raster