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