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)

###

================================================
É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> 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 (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
======================================================================

_______________________________________________
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 fornea cdigo mnimo reproduzvel.