Perfil Longitudinal do Rio Principal Utilizando o GRASS e o R

Por definição, o perfil longitudinal de um rio apresenta a relação existente entre a variação altimétrica e o comprimento do mesmo, desde a nascente até a foz ou ponto de confluência.

Neste tutorial daremos continuidade ao post "Caracterização Morfométrica de Bacias Hidrográficas utilizando o GRASS" onde a partir do rio principal, será gerado o gráfico do perfil longitudinal do mesmo, utilizando o GRASS e o R.

Com o GRASS aberto no terminal, digite os comandos a seguir:

	# 1) Criação de um camada de pontos, contendo os vértices do rio principal
	v.to.points -i -v -t in=out_dem_filled_mainchannel out=mchannel_pts type=line dmax=1000
	# 2) Exportação dos vértices para um arquivo txt
	v.out.ascii mchannel_pts > coords.txt
	# 3) Remoção dos pipes ("|") e da coluna cat do arquivo txt
	cat coords.txt | awk -F "|" '{print $1,''$2}' > xy.txt
	# 4) Substituição dos espaços por vírgulas
	sed -i "s/ /,/g" xy.txt
	# 5) Criação do arquivo txt contendo as distâncias e as elevações do perfil
	cat xy.txt | r.profile input=dem output=profile.txt
	# 6) Início da sessão do R
	R
	# 7) Importação do arquivo txt
	profile <- read.table('profile.txt', header=F, sep=" ",col.names=c('distance','elev'))
	# 8) Plotagem do gráfico
	plot(profile$distance/1000, profile$elev, axes='true', xlab='Distância (km)',
	ylab='Elevações (m)', main='Perfil Longitudinal do Rio Principal', type='l', col='blue')
	# 9) Adição do Grid de coordenadas
	grid()

view raw gistfile1.sh hosted with ❤ by GitHub

Como resultado, teremos o gráfico abaixo:

Referências:

• Raster profile along arbitrary line segments 
• Grass Gis et R: superposition des couleurs d'un raster quelconque sur un profil topographique élaboré à partir d'un MNT