[R-br] Problema para armazenar vértices de parcelas

ASANTOS alexandresantosbr em yahoo.com.br
Sexta Abril 13 17:42:26 BRT 2012


Boa tarde pessoal,

      Tenho os centróides em UTM de 11 parcelas quadradas de 40x40m e 
estou tentando calcular os vértices, inclusive com a repetição do 
primeiro vértice para fechar o polígono. Ao final, quero armazenar os 
resultados. Para tanto fiz:

xcc<-c(672036,672332,672674,672079,671386,671364,670753,671854,669986,670860,671366)##latitude 
dos centróides
ycc<-c(6783650,6784105,6783470,6782354,6782416,6782418,6782199,6783073,6781853,6781791,6781911)##longitude 
dos centróides

Fazendo o cálculo dos quatro vértices para cada par de coordenada:

  res=NULL
   for(i in 1:(length(xcc))){
     for(j in 1:(length(ycc))){

xa<-xcc[i]-20
ya<-ycc[j]+20
v1<-cbind(xa,ya)
xb<-xcc[i]+20
yb<-ycc[j]+20
v2<-cbind(xb,yb)
xc<-xcc[i]+20
yc<-ycc[j]-20
v3<-cbind(xc,yc)
xd<-xcc[i]-20
yd<-ycc[j]-20
v4<-cbind(xd,yd)
lim.p<-rbind(v1,v2,v3,v4,v1)
write.table(lim.p, file ="parc.mod", col.names=T, row.names=F,append=T)
result=read.table("parc.mod",header=F)
}}

O arquivo result esta saindo todo bagunçado e com mais de 60 linhas que 
deveriam corresponder a cabeçalho e cinco pares de coordenadas para 
todos os centróides, porém quando inspeciono o objeto lim.p, os cálculos 
estão corretos para o último par de coordenadas do centróide. Algum 
poderia me dizer onde posso estar errando.

Obrigado,

Alexandre



-- 
Alexandre dos Santos
Engenheiro Florestal, Dr.
Universidade Federal de Lavras
Departamento de Entomologia
Laboratório de Entomologia Florestal
Caixa Postal 3037
37200-000 - Lavras/MG
Fone: +55 (35) 9223-0304



Mais detalhes sobre a lista de discussão R-br