<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif">Alexandre, bom dia!</div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">Com seu código eu obtenho:</div><div class="gmail_default"><font face="monospace, monospace">> res(r)<br></font></div><div class="gmail_default"><div class="gmail_default"><font face="monospace, monospace">[1] 1.978261 2.823529</font></div></div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">Seria esperado 2 e 3. O problema parece estar na forma de associar seu grid com os limites do raster.</div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">Sugiro algo assim:</div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default"><font face="monospace, monospace">###</font></div><div class="gmail_default" style=""><div class="gmail_default"><font face="monospace, monospace">## Dados simulados</font></div><div class="gmail_default"><font face="monospace, monospace">dCol <- 3; nCol <- 10</font></div><div class="gmail_default"><font face="monospace, monospace">dLin <- 2; nLin <- 10</font></div><div class="gmail_default"><font face="monospace, monospace">injuria <- matrix(rnorm(nCol*nLin, mean=12), nLin, nCol)</font></div><div class="gmail_default"><font face="monospace, monospace">dim(injuria) </font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">## Coordenadas do grid espaçadas 3 x 2</font></div><div class="gmail_default"><font face="monospace, monospace">p0 <- c(680200,7803300)</font></div><div class="gmail_default"><font face="monospace, monospace">grid0 <- expand.grid(x=seq(p0[1], by=dCol, len=nCol), y=seq(p0[2], by=dLin, len=nLin))</font></div><div class="gmail_default"><font face="monospace, monospace">grid1 <- expand.grid(x=seq(p0[1], by=dCol, len=nCol+1), y=seq(p0[2], by=dLin, len=nLin+1))</font></div><div class="gmail_default"><font face="monospace, monospace">plot(grid0, asp=T, pch=21)</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">require(raster) </font></div><div class="gmail_default"><font face="monospace, monospace">r <-raster(injuria, xmn=p0[1], xmx=p0[1]+(nCol*dCol),</font></div><div class="gmail_default"><font face="monospace, monospace"> ymn=p0[2], ymx=p0[2]+(nLin*dLin), </font></div><div class="gmail_default"><font face="monospace, monospace"> crs=CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs"))</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">plot(r); points(grid0, pch=21) </font></div><div class="gmail_default"><font face="monospace, monospace">plot(r); points(grid1, pch=21) </font></div><div class="gmail_default"><font face="monospace, monospace">r</font></div></div><div class="gmail_default"><font face="monospace, monospace">###</font></div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><font face="arial, helvetica, sans-serif" style="font-size:small"><div style="font-family:'trebuchet ms',sans-serif;display:inline"></div>================================================<br>Éder Comunello</font><div style="font-size:small"><span style="font-family:arial,helvetica,sans-serif">PhD Student in Agricultural Systems Engineering (USP/Esalq)</span><br></div><div><span style="font-size:small">Brazilian Agricultural Research Corporation (</span><font face="arial, helvetica, sans-serif" style="font-size:small">Embrapa)</font><div style="font-size:small"><font face="arial, helvetica, sans-serif">Dourados, MS, Brazil [</font>22 16.5'S, 54 49.0'W<span style="font-family:arial,helvetica,sans-serif">]</span></div><div><div><br></div><div><br></div></div><div style="font-size:small"><br></div></div></div></div></div></div></div></div></div>
<br><div class="gmail_quote">Em 18 de janeiro de 2016 15:44, ASANTOS <span dir="ltr"><<a href="mailto:alexandresantosbr@yahoo.com.br" target="_blank">alexandresantosbr@yahoo.com.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Prezados Listeiros,<br>
<br>
<br>
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?<br>
<br>
<br>
##Dados artificiais<br>
injuria<-rnorm(1550, mean=12)<br>
<br>
##Crio a matriz com os dados artificiais<br>
B = matrix(<br>
data=injuria,<br>
nrow=25,<br>
ncol=62)<br>
<br>
# Inserindo coordenadas espaciais<br>
p0<-c(680234.482625,7803361.55128)# Ponto inicial<br>
c0 <- seq(p0[1],p0[1]+182,by=2)<br>
d0 <- seq(p0[2],p0[2]+50,by=3)<br>
d1 <- expand.grid(x=c0, y=d0) ## Coordenadas do grid espaçadas 3 x 2<br>
d1<-d1[-(1:14),] ## Removo algumas corrdenadas para ficar com o mesmo número de dados<br>
plot(d1[,1],d1[,2])<br>
length(d1[,1])<br>
<br>
require(raster)<br>
FORM<-cbind(d1,B)<br>
colnames(FORM)<-c('x','y','area')<br>
r <- rasterFromXYZ(as.data.frame(FORM)[, c("x", "y", "area")])<br>
extent(r) <- c(680234.5,680416.5,7803362,7803410)<br>
projection(r) <- CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs")<br>
plot(r)<br>
#<br>
<br>
<br>
<br>
Obrigado,<br>
<br>
-- <br>
======================================================================<br>
Alexandre dos Santos<br>
Proteção Florestal<br>
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso<br>
Campus Cáceres<br>
Caixa Postal 244<br>
Avenida dos Ramires, s/n<br>
Bairro: Distrito Industrial<br>
Cáceres - MT CEP: 78.200-000<br>
Fone: <a href="tel:%28%2B55%29%2065%208132-8112" value="+556581328112" target="_blank">(+55) 65 8132-8112</a> (TIM) <a href="tel:%28%2B55%29%2065%209686-6970" value="+556596866970" target="_blank">(+55) 65 9686-6970</a> (VIVO)<br>
<a href="mailto:e-mails%3Aalexandresantosbr@yahoo.com.br" target="_blank">e-mails:alexandresantosbr@yahoo.com.br</a><br>
<a href="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank">alexandre.santos@cas.ifmt.edu.br</a><br>
Lattes: <a href="http://lattes.cnpq.br/1360403201088680" rel="noreferrer" target="_blank">http://lattes.cnpq.br/1360403201088680</a><br>
OrcID: <a href="http://orcid.org/0000-0001-8232-6722" rel="noreferrer" target="_blank">orcid.org/0000-0001-8232-6722</a><br>
Researchgate: <a href="https://www.researchgate.net/profile/Alexandre_Santos10" rel="noreferrer" target="_blank">https://www.researchgate.net/profile/Alexandre_Santos10</a><br>
LinkedIn: <a href="https://br.linkedin.com/in/alexandre-dos-santos-87961635" rel="noreferrer" target="_blank">https://br.linkedin.com/in/alexandre-dos-santos-87961635</a><br>
======================================================================<br>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e fornea cdigo mnimo reproduzvel.</blockquote></div><br></div>