Fala Éder...
muito bom

Abraço

Hélio


2013/12/28 Eder Comunello [via R-br] <ml-node+s2285057n4661367h45@n4.nabble.com>
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 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>



_______________________________________________
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.


If you reply to this email, your message will be added to the discussion below:
http://r-br.2285057.n4.nabble.com/R-br-Recortar-raster-com-limite-em-data-frame-tp4661344p4661367.html
To unsubscribe from R-br, click here.
NAML



--
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho