[R-br] Criando uma mascara com shapefile sobre raster

ASANTOS alexandresantosbr em yahoo.com.br
Quinta Abril 4 12:31:55 BRT 2013


Boa tarde pessoal,

         Estou com um problema e não consigo criar uma mascara com um 
arquivo no formato shapefile sobre raster usando a função crop() do 
pacote raster, por algum motivo meu objeto contorno que é um 
SpatialLinesDataFrame não esta servido para oferecer o polígono 
limítrofe da área, alguém poderia me dar um help, segue CRM:

require("raster")
A <- stack(c("LANDSAT_5_TM_20100506_217_074_L2_BAND7.tif"))## Imagem
plotRGB(A) ## Plota a imagem
contorno<- readShapeLines("Catas_Altas.shp")  ## Contorno da área
lines(contorno, col="red")Plota o contorno
C <- crop(A,contorno) ## Recorta o que esta contido no contorno na imagem
Erro em .local(x, y, ...) :
   nenhum slot de nome "legend" para esse objeto de classe "RasterStack"


-- 
======================================================================
Alexandre dos Santos
Proteção Florestal
Coordenador do curso Técnico em Florestas
Vice Coordenador do curso de Engenharia 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
======================================================================



Mais detalhes sobre a lista de discussão R-br