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:[email protected]
[email protected]
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================
_______________________________________________
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.