<div dir="ltr"><div class="gmail_default" style=""><div class="gmail_default"><font face="verdana, sans-serif">Alterando grid para representar centro dos pixels...</font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div class="gmail_default"><font face="monospace, monospace">###</font></div><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) ### Lower Left</font></div><div class="gmail_default"><font face="monospace, monospace">p1 <- p0 + c(dCol/2, dLin/2)</font></div><div class="gmail_default"><font face="monospace, monospace">grid <- expand.grid(x=seq(p1[1], by=dCol, len=nCol), </font></div><div class="gmail_default"><font face="monospace, monospace"> y=seq(p1[2], by=dLin, len=nLin))</font></div><div class="gmail_default"><font face="monospace, monospace">plot(grid, 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(grid, pch=21) </font></div><div class="gmail_default"><font face="monospace, monospace">r</font></div><div class="gmail_default"><font face="monospace, monospace">###</font></div></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 19 de janeiro de 2016 09:43, 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">
<div bgcolor="#FFFFFF" text="#000000">
Muito obrigado Éder,<br>
<br>
Ficou muito bom mas era quase isso, pois a coordenada do grid1
deveria marcar o centro do pixel e não os vértices<span class=""><br>
<pre cols="72">--
======================================================================
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: <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)
<a href="mailto:e-mails:alexandresantosbr@yahoo.com.br" target="_blank">e-mails:alexandresantosbr@yahoo.com.br</a>
<a href="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank">alexandre.santos@cas.ifmt.edu.br</a>
Lattes: <a href="http://lattes.cnpq.br/1360403201088680" target="_blank">http://lattes.cnpq.br/1360403201088680</a>
OrcID: <a href="http://orcid.org/0000-0001-8232-6722" target="_blank">orcid.org/0000-0001-8232-6722</a>
Researchgate: <a href="https://www.researchgate.net/profile/Alexandre_Santos10" target="_blank">https://www.researchgate.net/profile/Alexandre_Santos10</a>
LinkedIn: <a href="https://br.linkedin.com/in/alexandre-dos-santos-87961635" target="_blank">https://br.linkedin.com/in/alexandre-dos-santos-87961635</a>
======================================================================
</pre>
</span><div><div class="h5"><div>Em 19/01/2016 09:28, Éder Comunello
escreveu:<br>
</div>
</div></div><blockquote type="cite"><div><div class="h5">
<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">
<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>
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr"><font style="font-size:small" face="arial, helvetica, sans-serif">
<div></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 style="font-size:small" face="arial,
helvetica, sans-serif">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>
<br>
<fieldset></fieldset>
<br>
</div></div><pre><div><div class="h5">_______________________________________________
R-br mailing list
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a></div></div>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forne�a c�digo m�nimo reproduz�vel.</pre>
</blockquote>
<br>
</div>
</blockquote></div><br></div>