
Boa Tarde Pessoal, Estou reorganizando minhas coordenadas geográficas (pontos) para calcular todas os arranjos possíveis de pontos em diferentes distâncias (código abaixo), até aí tudo bem, mas o problema é quando utilizo o comando apply(), dentro do loop para remover as linhas que contem NAs e ficar somente com as coordenadas em X com correspondencia em Y (xyc) e posteriormente adicionar os resultados em maille. Neste objeto maille que guarda o resultado de cada loop se eu adiciono c(xyc,n.pontos,Bk[m]) o resultado fica todo bagunçado, com n.pontos na mesma coluna que xyc , mais se deixo só xyc o resultado fica correto, tem alguma forma de fazercom que xyc,n.pontos,Bk[m] apareça individualmente em cada maille[i,], Obrigado, ############################################################################ ################### pt0 <- c(725990,7912746)#Ponto inicial h.b.a.f.p1<- expand.grid(latitude=seq(pt0[1], by=25, length.out=12),longitude=seq(pt0[2], by=25, length.out=12)) ##Malha x<-c(20,100,110,144) h.b.a.f.p1<-h.b.a.f.p1[-(x),]##Removendo alguns pontos para criar distâncias inexistentes Bk=seq(0,50,25)###Diferentes distâncias entre pontos amostrais maille<-NULL for(m in 1:length(Bk)){ Px<-c(h.b.a.f.p1[,1])###Coordenadas X Py<-c(h.b.a.f.p1[,2])###Coordenadas Y Pxy<-cbind(Px,Py) plot(Pxy) sqx<-seq(min(Px)+Bk[m],max(Px),25)###Sequência de distancias em X sqy<-seq(min(Py)+Bk[m],max(Py),25)###Sequência de distancias em Y for(j in 1:length(Px)){ for(k in 1:length(Py)){ n.pontos<-1:length(Px) xx=sqx[j] ###Coordenada X presente ou não nesta posição yy=sqy[k] ###Coordenada Y presente ou não nesta posição xyc=cbind(xx,yy) xyc=xyc[apply(xyc, 1, function(x)!any(is.na(x))), drop=F] ###Remove linhas se NA maille=rbind(maille,c(xyc,n.pontos,Bk[m]))###Resultado }}} # Alexandre dos Santos Ingenieur forestier, Msc. INRA- Biostatistique et Processus Spatiaux (BioSP) Domaine Saint-Paul Site Agroparc 84914 - Avignon - France Tél. : +33 (0)6 87 95 16 29
participantes (1)
-
Alexandre dos Santos