code{ -webkit-background-clip:initial; -webkit-background-origin:initial; background-attachment:initial; background-color:#FFF; background-image:url(http://arquivos.noticiaeblog.com/code.jpg); background-position:0 0; background-repeat:no-repeat; color:#666; display:block; font-size:13px; font:normal normal normal 1em/normal 'Courier New', Fixed, monospace; line-height:17px; overflow-x:auto; overflow-y:auto; text-align:left; border-color:#9c6; border-style:solid; border-width:1px; margin:1em 0; padding:0 20px 0 30px; }

terça-feira, 30 de agosto de 2011

Criando perfil topográfico no programa R

Abrir o modelo digital de elevação (MDE) no SAGA-GIS, de acordo com a extensão em que se encontra o arquivo. Para isso pode-se utilizar as ferramentas de import dos módulos. Posteriormente abrir a linha que servirá para amostrar os dados.


Para amostrar os pixeis do Raster será utilizado a ferramenta no Terrain Analysis – Profiles -> Profiles from lines


DEM: especificar o MDE
Values: Caso estivesse mais Grids poderia amostrá-los e isso iria criar colunas correspondentes
Lines: Qual o vetor no caso é Perfil
Name: Caso estivesse mais linhas no vetor como eu iria diferenciá-las
Profiles: A opção Create
O resultado dessa operação será um vetor de pontos que conterá na tabela colunas X,Y e Z, com a distância entre ponto de 30 metros, essa distancia é de acordo com pixel .


Após esse procedimento vamos exportar a tabela em texto com o módulo Import/Export -> Eport Text Table.E posteriormente abrir essa tabela no programa R, essa etapa de como exportar a tabela e de como abrir uma arquivo texto no R está especificado em outro post desse blog.
Com o arquivo aberto no Programa R podemos usar a funcao plot() para visualizar o perfil.
Código para efetuar a plotagem do gráfico
>plot(perfil$Z, type="l", main="Perfil topográfico", xlab="", ylab="Elevação")
Entao o nome da minha tabela no R é perfil, utilizei os dados dos campo Z, o tipo do gráfico é linha (type=”l”). A figura resultante se encontra abaixo.


Nenhum comentário:

Postar um comentário