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.
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