[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