Adicionar o R no SAGA,GRASS, QGIS ou ILWIS

Em algumas pesquisas pela internet descobri esse site que mostra como adicionar o programa R de estatística em softwares livres de SIG. Tenho o plugin do R no QGIS, mas ainda está acontecendo um erro quando tento abrí-lo.

attt!!

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.

Apresentando os dados através de gráficos no R, histograma e BoxPlot.

Para criarmos os gráficos primeiramente vamos abrir uma tabela que contén dados de altitude do relevo.
alt <-read.table(“C:/dados.txt”, h=T) Essas palavras chaves já foram vista em outra postagem. Como queremos colocar o histograma e o boxplot em uma mesma área gráfica, temos que usar os comando “par()” e “mfrow=”. O código usado nesse exemplo foi.

>par(mfrow=c(1,2))
> hist(alt$ALTITUDE, main="Dados ltimétricos", xlab="Altitudes", ylab="Frequência")
“main=” O título do histograma
“ylab=” Coluna das frequências
“xlab=”linha dos valores
> boxplot(summary(alt$ALTITUDE))
O gráfico final ilustrado abaixo.


> par(mfrow=c(1,2))
“mfrow=c(1,2)” : Cria um vetor que separa a área gráfica,o primeiro número separa informa a quantidade de divisoes horizontais, nesse caso foi 1. O segundo número informa a quantidade de divisoes na vertical, nesse caso 2, um é o histograma e outro o boxplot.
A etapa posterior é informar quais representacoes gráficas desejamos.
> hist(alt$ALTITUDE)
> boxplot(summary(alt$ALTITUDE))

O gráfico obtido está abaixo.

Outra questão interessante é colocar um título e melhorar a apresentação do histograma.

>par(mfrow=c(1,2))
> hist(alt$ALTITUDE, main="Dados ltimétricos", xlab="Altitudes", ylab="Frequência")
“main=” O título do histograma
“ylab=” Coluna das frequências
“xlab=”linha dos valores
> boxplot(summary(alt$ALTITUDE))
O gráfico final ilustrado abaixo.

Abrindo um arquivo texto no R e efetuando cálculos estatísticos básicos.

O arquivo está na repartição “C:/” com dados da declividade de um terreno. Para abrir o arquivo usamos o comando:

Declividade <- read.table(“C:/dados.txt”, h=T)

“/” = deve-se usar esse sentido da barra


“h=T” = Significa que a primeira linha será o nome da coluna.


sep = “ ; ”= caso tenha um separador entre colunas que não seja tabulações.

Faremos alguns cálculos estatísticos com o comando “ summary( )”

est <- summary(Declividade$DECLIVIDADE)

Tabela com nome Declividade(primeira letra maiúscula) e campo da tabela com nome DECLIVIDADE(todas letras maiúscula).


“$”= permite acessar individualmente o campo da tabela

Criando um boxplot com os dados.”boxplot( )”

boxplot(est)

Bom trabalho a todos!!