Alexandre,Atualizei seu CMR, adicionando as últimas linhas para realizar o recorte com um polígono irregular.### <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 SpatialPolygonspt.borda <- SpatialPolygons(list(Polygons(list(Polygon(limite)), "1")))### Recorte do raster contido no interior da bordar.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>
_______________________________________________
R-br mailing list
[hidden email]
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.
http://r-br.2285057.n4.nabble.com/R-br-Recortar-raster-com-limite-em-data-frame-tp4661344p4661367.htmlIf you reply to this email, your message will be added to the discussion below: