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

Atribuindo valores de grids para shapefile com SAGA-GIS

Modules ->Shapes-Grid -> Grid Values to Shapes


Grid system -> sistema de grids
Grids -> pressione as reticências que estão no final da linha.


Selecione os grids que serão utilizados e pressione a seta que a ponta para direita.


Polygons -> Escolha o polígono que será utilizado como máscara para extrair os valores dos grids
Shapes -> Selecione a palabra “CREATE”
Marque a opção Exclude NoDat Cells
Type -> cells

Como usar o Grid Calculator do SAGA-GIS.

Nesse exemplo iremos calcular o Índice de úmidade topográfica, utilizando o Grid Calculator.





Janela modules -> Grid-Calculus -> Grid Calculator



Selecione o sistema de Grid
Grids (Pressione as reticências no final da linha para adicionar os grids que serão usados nos cálculos).



Selecione o grid do lado esquerdo e presione a flexa que está apontada para a direita, assim o grid será inserido no lado direio e usado na fórmula pretendida. A ordem final do lado direito terá como referencia as letra (a,b,c,d,e,f…), assim por diante. Isso significa que quando escrever a fórmula terá que usar letras para se referir as camadas grids.



Result-> create (sempre quando querer um novo resultado terá que escolher “Create”)
Formula -> ln((a/30)/(tan(b))), Logaritmo natural de ((área de captação /resolução do pixel)/(tangente da declividade(slope)))
a -> área de captação
b -> declividade
Name -> IUT

Revista Inside GNSS

Um site interessante para quem gosta de estar informado sobre novas tecnologias que envolvam posicionamento por satélite.

http://www.insidegnss.com/

Georreferenciamento de carta topográfica no Erdas 9.1

Gerando atributos do terreno no SAGA-GIS

Esse tutorial explicará como gerar atributos do terreno no SAGA-GIS, utilizando o modelo digital de elevação explicado na postagem anterior.


Modules -> Terrain Analysis – Compound Analyses


Grid system -> Grupo de camadas
Elevation -> Modelo digital de elevação (triangulação)
Todos os itens abaixo devem estar com a opção CREATE . Entre esses atributos do terreno estão alguns muito conhecidos como: curvatura planar e vertical, sombreamento onde sua configuração está sendo mostrada no quadro azul abaixo (azimute e inclinação), declividade do terreno, Aspecto do terreno, índice de umidade topográfica, delineamento de bacias hidrográficas e outros. No quadro verde será criada a rede de drenagem.