Histórico Topografia
O trabalho da topografia em campo inciou-se em novembro de 2003, coordenado pelo Alexandre Farelo, nosso primeiro estagiário. A principio contratamos um topógrafo que acompanhou os trabalhos em campo. Farelo foi treinado em BCI antes de iniciar os trabalhos em Manaus.
Foram 15 expedições a campo, de uma semana cada, sendo feitas 1395 medidas entre piquetes na parcela. Os dados originais estão no arquivo a seguir:
Ao trabalhar os dados em 2006 percebi uma redundância de informações, com três dados repetidos que exclui no cálculo das altitude. O arquivo seguinte contem os dados excluindo os tres registro. UTILIZEI ESSE ARQUIVO
{manausdeltatopo21abr06.xls|}}
Utilizando a função solve.topo, criada por Richard Condit, transformamos os dados da topografia em altitude relativa ao inicio da parcela.
atençao: retirar a extensão pdf
A partir dessa função criamos o arquivo de elevação do terreno abaixo onde todos os piquetes da parcela correspondem a uma elevação do terreno, relativo à origem da parcela (00×00).
atençao: retirar a extensão pdf
Para criar mapas com maior definição, fizemos uma krigagem dos dados para que tivessemos informação de elevação a cada 5 metros a partir dos dados coletados a cada 20 metros. Na época utilizei a função kriging.r e também utilities.r .
Há necessidade de carregar o pacote sgeostat no R.
library(sgeostat)
Abaixo os argumentos que usei para outra parcela…
kriging(elevorig=km41elev,plotdim=c(1000,1000),num.lags=3,maxdist.pair=151,type.variogram=“exponential”,maxdist.krige=61,iter=30,range=100,nugget=2,sill=20,plot.it=F,filepath=“C:/km415x5”)
NÃO LEMBRO O PORQUE, MAS LEMBRO DE TER TESTADO VÁRIOS MÉTODOS
O arquivo criado foi:
A partir dele é possível fazer os mapas. Primeiro é preciso criar a matrix
mat.elev5=matrix(elev5$elev, nrow=141)
str(mat.elev5)
Depois é só fazer o mapa
contour(seq(0,700,by=5),seq(0,700,by=5),t(mat.elev5),nlevels=35)
ou
persp(t(mat.elev5),phi=40,theta=30,d=0.8,asp=1,scale=T,expand=0.2,col=“green”, xlab=“ ”, ylab=“ ”, zlab=“alt”, axes=F)
2009 Fevereiro
Função para plotar o mapa com simples controles de cor [usa rainbow()] e do numero de niveis desejados.
Exemplo:
plot.topo.map(obj=t(mat.elev5),opacity=0.8,end.col=0.7,start.col=0,invert=T,n.intervals=10)
Função:
plot.topo.map <- function(obj=topo.map,n.intervals=10, start.col=0.1,end.col=max(1,n.intervals-1)/n.intervals,opacity=0.8,invert=T) {
n = n.intervals
col.vec = rainbow(n, start = start.col, end = end.col, gamma = 1, alpha = opacity)
if (invert==T) {cl=rev(col.vec)} else {cl=col.vec}
filled.contour(obj,nlevels=n,col=col.vec,plot.axes=F)
}