[R-br] Seleção de pixels em raster em função de coordenadas não inteiras

ASANTOS alexandresantosbr em yahoo.com.br
Terça Janeiro 14 10:00:10 BRST 2014


Caros Membros

          Gostaria de fazer a seleção e representar graficamente os 
pixels de um um raster (r) que contem pontos de interesse (pontos). No 
entanto, o raster possui coordenadas com valores inteiros (Ex. 6, 8, 9 
..) e as coordenadas dos pontos possuem subunidades(6.03, 8.22, 9.65), 
gostaria de saber como posso fazer para extrair os pixels que contem os 
pontos selecionados e sua representação gráfica, segue CRM que não 
funcionou,


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(4.25,3.06,6.22,7.35)
y<-c(7.00,8.90,1.18,7.03)
pontos <- SpatialPoints(cbind(x,y))#cbind-junta duas colunas em uma, no 
caso as coordenadas


### Recorte dos pontos de interesse do raster
r.crop <- crop(r, pontos, snap='out')
image(r.crop, asp=1)

### Recorte dos pixels nas coordenadas de interesse
r.masked <- mask(r.crop, pontos)
image(r.masked, asp=1)
###


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 (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
======================================================================



Mais detalhes sobre a lista de discussão R-br