[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