MDE : Anudem - Recentes avanços na modelagem do relevo

O Anudem é um interpolador presente no ArcGis com o nome TopotoRaster. A seguir é apresentado um pequeno trecho traduzido do artigo uqe Hutchinson apresentou no Geomorphometry 2011 sobre as mais recentes evoluções do Anudem, que normalmente são incorporadas nas próximas versões do ArcGis.

"O relevo desempenha um papel fundamental na modulação da superfície da terra e processos atmosféricos em uma ampla gama de escalas espaciais (Hutchinson 2008). Assim, modelos de elevação digital (DEM) têm desempenhado um papel fundamental no apoio representações de mesoescala do clima, bem como no apoio a representações mais detalhadas da hidrologia de superfície. O ANUDEM(Hutchinson 1989, 2007) é comumente usada para calcular esses modelos de elevação em forma de grade regular.Principais características do método incluem a sua eficiência computacional, permitindo-lhe ser aplicada a conjuntos de dados muito grandes, e uma gama de características localmente adaptativas, incluindo um algoritmo de imposição da drenagem que  mantem a estrutura de drenagem conectada no DEM interpolados, e algoritmos para incorporar dados de  lagos e falésias. Este artigo descreve o progresso atual do processo ANUDEM para melhor representar lagos e efetivamente mehorar a interpolação de dados ruidosos provenientes de fontes de alta resolução. Esses dados estão se tornando cada vez mais comum. O procedimento de interpolação de base multi-grade permanece eficaz em representar efectivamente lagos e falésias e estavelmente interpolando dados de elevação de alta resolução. Erros correlacionados em dados de elevação também podem ser especificamente acomodados.”

Leia o artigo completo em ingles formato PDF

Interpolando o SRTM de pixels com 90m para 10m

 

 

 

 

 

 

 

Com o SRTM em mão, dá-se início ao seguinte processo:

Será feito o seguinte: o Arquivo Raster será transformado para pontos, que será interpolado, originando um Modelo digital de elevação mais refinado.

Caso você não deseje trabalhar com toda a cena SRTM, você pode recorta-la para diminuir o tamanho do arquivo, e assim processa-lo mais rapidamente.

Iremos recorta-lo da seguinte maneira, já que apenas uma pequena área nos interessa, caso a você seja interessante toda a cena, apenas desconsidere esta parte. Vamos dar um zoom com a lupa na parte da cena que nos interessa.

Após ter dado o zoom clique com o boto direito do mouse no arquivo do SRTM, clique em DATA e posteriormente em EXPORT DATA, e uma janela se abrirá.

Nesta tela em EXTEND selecione “Data Frame” e o formato do arquivo escolha TIFF. Após isto escolha local onde o arquivo será gravado e clique em OK.

Agora com o arquivo menor faça os seguintes passos para interpolar o SRTM.

 

Em Spatial Analyst Clique em CONVERT > RASTER TO FEATURES

 

Em RASTER TO FEATURES escolha qual arquivo será convertido. Em OUTPUT GEOMETRY TYPE escolha POINT. Selecione o local onde o arquivo será gravado e clique em OK, e um arquivo de pontos será criado contendo a informação altimétrica.

 

Agora com este arquivo de pontos gerado, iremos interpola-lo para uma resolução menor.

Clique em ArcTool Box > SPATIAL ANALYST TOOLS > INTERPOLATION > SPLINE. Uma outra tela se abrirá, nela basta escolher qual o arquivo de pontos que será utilizado. O Campo na tabela que os dados altimétricos, geralmente seu nome é GRID_CODE, em OUTPUT CELL SIZE, escolha o tamanho para o qual você deseja interpolar o SRTM. Basta então clicar em OK

 

ArcGis : Como são feitas operações entre grids de diferentes tamanho de célula ?

Muita gente que trabalha com grid no ArcGis nem se da conta, mas muitas vezes estamos realizando operações entre grids com diferentes tamanhos de células. Numa operação entre dois grids de tamanho de células diferentes, qual será o tamanho de célula do grid resultante ? Depende da opção que foi escolhida na aba Cell Size do menu options do spatial analyst. Existem três opções : A primeira e Maximum of Inputs, neste caso é o tamanho da célula de saída é o do grid de maior célula, a segunda opção é Minimum of inputs, aqui é escolhido o menor tamanho de célula, e a ultima opção é As specified below onde o próprio usuário escolhe o tamanho de célula do grid de saída. Outra coisa importante é que o ArcGis interpola os grids de entrada para a mesma resolução do grid de saída usando o método do vizinho mais próximo.