[R-br] Conversão de matriz em raster georeferenciado

ASANTOS alexandresantosbr em yahoo.com.br
Segunda Janeiro 18 16:44:54 BRST 2016


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 em yahoo.com.br
         alexandre.santos em 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
======================================================================



Mais detalhes sobre a lista de discussão R-br