Recorte, reclassificação e vetorização de raster no gvSIG
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):
Para recortarmos a camada raster, abrimos a "Caixa ferramentas" e selecionamos o geoprocesso "Recortar camada raster com camada de polígonos".
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".
E teremos a nossa camada de estudo já inserida na área de trabalho:
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:
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":
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".
Abrirá uma nova janela, onde poderemos estabelecer os limites que desejamos aplicar à nossa 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.
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:
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:
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":
O único parâmetro que precisamos definir é a "Camada de entrada". Cuide para selecionar a camada correta.
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:
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: ). Visualizaremos a tabela de atributos da nossa camada:
Em seguida, clique no menu "Tabela > Adicionar medição > Área" (ou utilize o ícone: ). 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.
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: .
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".
Na área de trabalho poderemos visualizar as áreas selecionadas:
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: .
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:
Para calcular as áreas das outras classes devemos repetir o procedimento utilizado.
Tags: gvSIG, geoprocessos, geoprocessamento, raster