
Prezados Listeiros, Gostaria de realizar a conversão de uma matriz em raster georeferenciado com duas informações, a representação de plantas em uma matriz (B) e juntar as informações geográficas conhecendo-se o espaçamento entre plantas (d1), no entanto, gostaria de saber se conforme o meu CRM abaixo essa conversão esta correta, pois não estou bem certo de que realização de um cbind() entre os objetos B e d1 correspondam a real posição desses elementos. Alguém poderia me ajudar? ##Dados artificiais injuria<-rnorm(1550, mean=12) ##Crio a matriz com os dados artificiais B = matrix( data=injuria, nrow=25, ncol=62) # Inserindo coordenadas espaciais p0<-c(680234.482625,7803361.55128)# Ponto inicial c0 <- seq(p0[1],p0[1]+182,by=2) d0 <- seq(p0[2],p0[2]+50,by=3) d1 <- expand.grid(x=c0, y=d0) ## Coordenadas do grid espaçadas 3 x 2 d1<-d1[-(1:14),] ## Removo algumas corrdenadas para ficar com o mesmo número de dados plot(d1[,1],d1[,2]) length(d1[,1]) require(raster) FORM<-cbind(d1,B) colnames(FORM)<-c('x','y','area') r <- rasterFromXYZ(as.data.frame(FORM)[, c("x", "y", "area")]) extent(r) <- c(680234.5,680416.5,7803362,7803410) projection(r) <- CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs") plot(r) # 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 OrcID: orcid.org/0000-0001-8232-6722 Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================