Neste tutorial vamos exemplificar o procedimento de Interpolação pela Ponderação do Inverso da Distância (IDW) com o software ArcGIS. Vamos fornecer os dados usados no passo a passo para que você possa reproduzir a atividade.
INTERPOLAÇÃO IDW NO ARCGIS
Neste exercício nós utilizamos como dados de entrada um Esri Shapefile de geometria pontual com 223 feições que entre outros dados possui os valores da população de cada um dos municípios da Paraíba em um determinado ano. Além disso, temos o arquivo com a divisão municipal do Estado.
Para fazer o download dos dois arquivos, clique aqui.
Para começar, vamos definir nossa área máxima a ser processada, ou seja vamos garantir que o raster gerado não ultrapasse a extensão geográfica da camada com a divisão municipal da Paraíba.
Acesse o menu Geoprocessing → Environments...
Nesta janela vamos informar duas vezes qual camada de dados é a base de nossa extensão geográfica máxima, que não deve ser extrapolada.
Primeiro, indique em Processing Extent, no parâmetro Extent a opção Same as layer pb (neste exemplo, nossa camada se chama pb). Ainda não clique ainda no OK. Baixe a barra de rolagem até o ponto mostrado na imagem a seguir.
Quando encontrar o parâmetro Raster Analysis indique novamente a mesma camada (pb, neste caso). Depois, clique em OK.
Agora, acesse o ArcToolBox. Uma das opções para abrir este conjunto de ferramentas é através do caminho mostrado abaixo ou do ícone vermelho que aparece nela abaixo do menu Help.
No ArcToolBox navegue pelo caminho Spatial Analyst Tools → Interpolation → IDW. Dê um duplo clique sobre esta ferramenta.
Na janela IDW temos três parâmetros principais para informar. Em Input point features indicamos a camada de pontos que possuem os dados para interpolação. Neste caso, a camada a ser escolhida se chama pb_pop.
Informamos em Z value field o campo da tabela de atributos que será usado como base para a interpolação. Neste exemplo indicamos o campo POPTOTAL.
Por fim, em Output raster indicamos o local onde será gravado o arquivo, que será um dado raster. Aqui escolhemos o formato TIF. Depois clicamos em OK.
A imagem abaixo ilustra o resultado da interpolação com o raster sob as camadas da divisão municipal (com transparência) e de pontos.
Note que já podemos ver de forma clara que no Leste do Estado a capital João Pessoa se destaca totalmente das demais áreas da camada, sendo Campina Grande o município que mais se aproxima deste resultado.
As cores que estão no mapa abaixo representam uma distribuição por valores iguais, não sendo a mais adequada para este tipo de espacialização, visto que os demais municípios do Estado tem um número bem inferior de habitantes do que os destacados no parágrafo acima.
Na imagem a seguir aplicamos um zoom na área de João Pessoa, onde estão localizados a maior parte da população da Paraíba. Vamos alterar as propriedades visuais da camada para uma melhor visualização da distribuição da população ao longo do estado.
Acesse as propriedades da camada raster. Vá em Symbology → Classified e depois clique em Classify….
Onde se pede o método (Method) de classificação altere para Desvio Padrão (Standard Deviation). Depois no tamanho do intervalo (Interval Size) escolha a opção mostrada na imagem a seguir (lembre que isso é apenas uma sugestão de valores). Depois, clique em OK.
Logicamente, é possível alterar a rampa de cores à vontade. Veja no exemplo a seguir a visualização do resultado das configurações ilustradas nos passos acima, inicialmente com a camada de pontos ativada.
Desabilitando a camada pontual, temos a seguinte visualização do resultado obtido.
Em conclusão, veja o resultado visual próximo da região metropolitana de João Pessoa/PB com um certo nível de zoom aplicado.
O que achou deste tutorial? Logo iremos produzir outros conteúdos relacionados com interpolação espacial e suas aplicações, tanto usando o ArcGIS como o QGIS. Aguardem!
Recomendamos que leia também a seguinte lista de dicas que preparamos para vocês:
- Análise Geoestatística de Casos de Dengue
- Como fazer Mapas de Kernel no ArcGIS
- Autocorrelação Espacial: Mapas de Cluster no ArcGIS
Por favor, deixem seus comentários.
10 respostas
Olá, Primeiramente queria agradecer e parabenizar você pelo tutorial sobre a Interpolação pela Ponderação do Inverso da Distância (IDW). Como graduando de engenharia agrícola estou desenvolvendo um trabalho onde expressamos as linhas de produção de forragem em mapa interpolados dos estados de Minas Gerais e Ceará. Eu gostaria de saber se tem como padronizar uma unica escala para diversos mapas interpolados cada um, com um campo diferente da tabela de atributos.
Obrigado!!! <3
Parabéns pelo tutorial Anderson!
Sou acadêmico de oceanografia e uso seus tutoriais frequentemente!
Abraço
Olá Anderson, muito bacana esse tutorial! Gostaria de saber quais procedimentos preciso realizar para transformar este tif em um vetor do tipo polígono.
Muito obrigada.
No ArcGIS, na ferramenta Search digite “raster to polygon”
Olá
Quero usar o método de IDW ´pelo arcgis no meu trabalho.
você teria algum material para estudo, para que eu tivesse uma melhor compreensão ?
Se sim, poderia mandar para o email amandasuarez00@gmail.com ?
Agradeço o contato.
Aguardo retorno.
Veja os “help” do ArcGIS. São fantásticos!
Obrigado. Fiquei procurando a solução para a minha interpolação IDW em sites gringos, mas a resposta que queria achei no seu tutorial.
Este tutorial é muito bom.
Clayton, tudo bem?
Agradeço pela sua visita e feedback.
Ficamos felizes em ver que o tutorial foi útil para você.
Abraço!