Obrigado Éder,
Queria também extrair o valor dos pixels e adicionei a rotina,
vou deixar o scrip para depósito na lista,
### <code r>
require(raster); require(sp)
## Criando um raster
r <- raster(nc=10, nr=10)
r <- setValues(r, round(runif(ncell(r))* 255))
##Pontos selecionados
x <- c(44.25,33.06,63.22,71.35)
y <- c(22.00,81.90,15.18,71.03)
pontos <- SpatialPoints(cbind(x,y))
cells <- cellFromXY(r, pontos); cells
coordinates(r)[cells,] ## alternativa 1
p.cen <- xyFromCell(r, cells); p.cen ## alternativa2
plot(r)
points(pontos, col="red")
points(p.cen)
values <- extract(r,p.cen) ##Extrai as coordendas
dados<-cbind(p.cen,cells,values)
dados
### </code>
--
======================================================================
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
======================================================================
Em 16/01/2014 21:45, Éder Comunello escreveu:
Alexandre,
Pelo que entendi você quer a coordenada do centro do pixel/célula
identificado. Se for isso mesmo posso sugerir duas opções, conforme
código abaixo:
### <code r>
require(raster); require(sp)
## Criando um raster
r <- raster(nc=10, nr=10)
r <- setValues(r, round(runif(ncell(r))* 255))
##Pontos selecionados
x <- c(44.25,33.06,63.22,71.35)
y <- c(22.00,81.90,15.18,71.03)
pontos <- SpatialPoints(cbind(x,y))
cells <- cellFromXY(r, pontos); cells
coordinates(r)[cells,] ## alternativa 1
p.cen <- xyFromCell(r, cells); p.cen ## alternativa2
plot(r)
points(pontos, col="red")
points(p.cen)
### </code>
Éder Comunello <c
<mailto:[email protected]>[email protected]
<mailto:[email protected]>>
Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________
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.
