Obrigado Helio e Éder, Problema resolvido, Em 28/12/2013 15:32, Hélio Gallo Rocha escreveu:
### <code r> require(sp); require(rgdal); require(raster)
## Raster----------------------------------------------------------- set.seed(20131226) r <- raster(nc=100, nr=100) r <- setValues(r, round(runif(ncell(r))* 255)) image(r,col=gray(seq(0,1,0.025)))
## Borda ----------------------------------------------------------- x<-c(-125,-100,0,100,-125) y<-c(-49,25,0,-50,-49) limite<-cbind(x,y) lines(limite, col="yellow") ##
## Converte borda para SpatialPolygon spt.borda <- SpatialPolygons(list(Polygons(list(Polygon(limite)), "1")))
### Recorte do raster contido no interior da borda r.crop <- crop(r, spt.borda, snap='out') image(r.crop, asp=1)
### Recorte com polígono irregular! r.masked <- mask(r.crop, spt.borda) image(r.masked, asp=1) ### </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@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================