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; }

quarta-feira, 24 de agosto de 2011

Programa para planejamento de voo em aerolevantamento

Uma das matérias ensinadas nas disciplinas fotointerpretação ou fotogametria é o planejamento de voo. Este post mostrará como usar um simples programa escrito em python para efetuar esses cálculos. Os cálculos disponíveis são:
- Escolha da escala da foto, levando em consideração a acuidade visual( 0.2 mm)
- Altura do voo
- Área útil do aerofotograma
- Número de fotografia por faixas
- Número de fotografias da cobertura aérea
- Tempo entre fotografias
- Tempo da cobertura aérea (Tempo de manobra, tempo de fotografia e o tempo total)
Os valores de segurança para que o voo exceda toda área, devem serem atribuídos pelo usuário. Por exemplo fotos por faixa.

Para o programa rodar deve ser baixado o Python.
Esse programa é leve e de fácil manuseio.

Depois do programa instalado deve-se abrir e copiar o código fonte

Pressione botao direito para abrir o arquivo.
O separador decimal é o ponto e o separador entre números é o ponto-vírgula

Qualquer dúvida deixem um comentário.!!

domingo, 14 de agosto de 2011

Iluminação de superfícies com SAGA-GIS

O método de simulação do sol em modelo digital de elevação (MDE) pode ser gerado com a ferramenta Analytical Hillshading do SAGA-GIS. Caso esteja trabalhando com MDE em coordenadas geográficas deve-se colocar na opção exageration o fator corresponde a latitude em que se encontra a área de estudo.


Os dados sobre o azimute e inclinação do sol de uma determinada data e área pode ser obtido no site:Cálculos. Deve-se ter cuidado com a latitude, longitude e fuso horário em que se encontram a área de estudo.


A figura abaixo representa o sombreamento do município de Erechim/RS, para a data 1 de setembro ás 15:30 hrs. Qualquer dúvida deixe um comentário!!! Bom Trabalho!!!!!






quarta-feira, 8 de junho de 2011

Configurando o GPS de navegação

Neste link abaixo você encontrará passos para configurar seu GPS de navegação, para entender como se procede a configuração do GPS de navegação, pressione as imagens das telas do GPS que aparecem no site.

http://www.quoos.com.br/Geografia/gps-garmin-etrex.html


até mais!

terça-feira, 31 de maio de 2011

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.

quinta-feira, 26 de maio de 2011

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

domingo, 24 de abril de 2011

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

sexta-feira, 8 de abril de 2011

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