[R-br] Cálculo de vértices à partir do centróide

[ASANTOS]

Boa tarde pessoal,

       Bom descobri onde estou errando, quando faço:

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

     O loop vai pegar a primeira coordenada de xcc e fazer todas as 
combinações com ycc. Na verdade, eu queria uma ajuda para o loop  e 
fazer com que ele pegue apenas o primeiro valor de xcc e comparar apenas 
com o primeiro valor de ycc, seguido de xcc[2] com ycc[2] ate xcc[11] 
com ycc[11], meu script completo é:

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

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

Obrigado

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