Muito obrigado Éder, Ficou muito bom mas era quase isso, pois a coordenada do grid1 deveria marcar o centro do pixel e não os vértices -- ====================================================================== 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 ====================================================================== Em 19/01/2016 09:28, Éder Comunello escreveu:
Alexandre, bom dia!
Com seu código eu obtenho:
res(r) [1] 1.978261 2.823529
Seria esperado 2 e 3. O problema parece estar na forma de associar seu grid com os limites do raster.
Sugiro algo assim:
### ## Dados simulados dCol <- 3; nCol <- 10 dLin <- 2; nLin <- 10 injuria <- matrix(rnorm(nCol*nLin, mean=12), nLin, nCol) dim(injuria)
## Coordenadas do grid espaçadas 3 x 2 p0 <- c(680200,7803300) grid0 <- expand.grid(x=seq(p0[1], by=dCol, len=nCol), y=seq(p0[2], by=dLin, len=nLin)) grid1 <- expand.grid(x=seq(p0[1], by=dCol, len=nCol+1), y=seq(p0[2], by=dLin, len=nLin+1)) plot(grid0, asp=T, pch=21)
require(raster) r <-raster(injuria, xmn=p0[1], xmx=p0[1]+(nCol*dCol), ymn=p0[2], ymx=p0[2]+(nLin*dLin), crs=CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs"))
plot(r); points(grid0, pch=21) plot(r); points(grid1, pch=21) r ###
================================================ Éder Comunello PhD Student in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa) Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]
Em 18 de janeiro de 2016 15:44, ASANTOS <alexandresantosbr@yahoo.com.br <mailto:alexandresantosbr@yahoo.com.br>> escreveu:
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 <tel:%28%2B55%29%2065%208132-8112> (TIM) (+55) 65 9686-6970 <tel:%28%2B55%29%2065%209686-6970> (VIVO) e-mails:alexandresantosbr@yahoo.com.br <mailto:e-mails%3Aalexandresantosbr@yahoo.com.br> alexandre.santos@cas.ifmt.edu.br <mailto:alexandre.santos@cas.ifmt.edu.br> Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 <http://orcid.org/0000-0001-8232-6722> Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635 ======================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br <mailto: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 fornea cdigo mnimo reproduzvel.
_______________________________________________ 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.