Classificação de pontos que pertencem a uma condição dada (x,y)

Boa tarde Pessoal, Fiz a representação de 48 parcelas com plantas no espaçamento de 2x4m usando a função expand.grid(), sendo as parcelas espaçadas em 7m (grid1), após criei pontos que presentam o centro das parcelas (grid2) e sobre esses pontos eu adicionei os códigos das parcelas (parc.names). Gostaria de saber se existe uma forma de criar uma classificação onde quero saber a qual código da parcela (parc.names) pertence cada um dos pontos que representam as plantas (grid1). Segue CRM: ##Plantas ---------------------------------------------------------------------- ## Coordenadas do grid no espaçamento 2x4 ## distâncias verticais e horizontais v <- 2; h <- 4 ## número de pontos vertical e horizontal npv <- 4; nph <- 5 ## distância entre vãos vertical e hotizontal V <- 7; H <- 7 ## número de linhas e colunas ncol <- 4; nrow <- 12 grid0 <- expand.grid(L=seq(2,length.out=nrow,by=npv*v+V), C=seq(2,length.out=ncol,by=nph*h+H), x=seq(2,length.out=npv,by=v), y=seq(2,length.out=nph,by=h)) grid1 <- with(grid0, data.frame(x=L+x, y=C+y)) plot(y~x, grid1, asp=1) ## Centro das parcelas ------------------------------------------------------ grid2 <- expand.grid(x1=seq(7,length.out=12,by=15), y1=seq(12,length.out=4,by=27)) points(grid2[,1],grid2[,2],col="red") #Código das parcelas parc.names<-c("2211","1221","3221","1122","2122","3222","1113","2113","3223","1124","2214","3114", "2111","1121","3211","1112","2212","3212","1213","2123","3213","1214","2124","3224", "2221","1211","3111","1212","2222","3112","1223","2213","3113","1224","2224","3214", "2121","1111","3121","1222","2112","3122","1123","2223","3123","1114","2114","3124") text(grid2$x1,grid2$y1,labels=parc.names,cex=1,col=2) # Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

Basta você aproveitar o que tem em grid0, que tem as colunas e linhas de blocos de árvores, combiná-los para fazer um fator. Confira no CMR. str(grid0) unique(grid0$L) unique(grid0$C) grid1$parc <- with(grid0, interaction(L, C)) nlevels(grid1$parc) levels(grid1$parc) grid1$parc <- with(grid0, interaction(as.integer(as.factor(L)), as.integer(as.factor(C)), sep="-")) nlevels(grid1$parc) levels(grid1$parc) grid2 <- aggregate(cbind(x, y)~parc, data=grid1, FUN=mean) str(grid2) plot(y~x, grid1, col=as.integer(grid1$parc)) with(grid2, text(x, y, labels=as.character(parc))) subset(grid1, parc=="1-1") À disposição. Walmes.

Genial Walmes!! Muito Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== On 09/09/2014 20:06, walmes . wrote:
Basta você aproveitar o que tem em grid0, que tem as colunas e linhas de blocos de árvores, combiná-los para fazer um fator. Confira no CMR.
str(grid0) unique(grid0$L) unique(grid0$C)
grid1$parc <- with(grid0, interaction(L, C)) nlevels(grid1$parc) levels(grid1$parc)
grid1$parc <- with(grid0, interaction(as.integer(as.factor(L)), as.integer(as.factor(C)), sep="-")) nlevels(grid1$parc) levels(grid1$parc)
grid2 <- aggregate(cbind(x, y)~parc, data=grid1, FUN=mean) str(grid2)
plot(y~x, grid1, col=as.integer(grid1$parc)) with(grid2, text(x, y, labels=as.character(parc)))
subset(grid1, parc=="1-1")
À disposição. Walmes.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
participantes (2)
-
ASANTOS
-
walmes .