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:[email protected]
[email protected]
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
<[email protected]
<mailto:[email protected]>> 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:[email protected]
<mailto:e-mails%[email protected]>
[email protected]
<mailto:[email protected]>
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
[email protected] <mailto:[email protected]>
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
[email protected]
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.
_______________________________________________
R-br mailing list
[email protected]
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.
