[R-br] Recortar raster com limite em data.frame [RESOLVIDO]

[ASANTOS]

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 em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131229/e13e260a/attachment.html>