[R-br] Identificar células em um raster [RESOLVIDO]

[ASANTOS]

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:alexandresantosbr em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
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:comunello.eder em gmail.com>omunello.eder em gmail.com 
> <mailto:omunello.eder em gmail.com>>
> Dourados, MS - [22 16.5'S, 54 49'W]
>
>
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em 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 forneça código mínimo reproduzível.

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140117/93d7a5fa/attachment-0001.html>