Alexandre, a) Tente ler o seu shapefile com a funcao readShapePoly. b) Se nao der certo, sugiro uma combinacao das funcoes crop e raster. Nesse caso, experimente rodar o seu codigo com as seguintes mudancas: require("raster") A <- stack("LANDSAT_5_TM_20100506_217_074_L2_BAND7.tif") contorno<- readShapePoly("Catas_Altas.shp") A.crop <- crop(A, extent(contorno), snap='out') contorno.na<-setValues(A.crop, NA) contorno.r<-rasterize(contorno, contorno.na) ### isso pode demorar, dependendo do seu shapefile C <- mask(x=A.crop, mask=contorno.r) Saudações, -- Thiago V. dos Santos PhD student Land and Atmospheric Science University of Minnesota http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/inde... Phone: (612) 323 9898 ----- Original Message ----- From: ASANTOS <alexandresantosbr@yahoo.com.br> To: r-br@listas.c3sl.ufpr.br Cc: Sent: Thursday, April 4, 2013 10:31 AM Subject: [R-br] Criando uma mascara com shapefile sobre raster 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@yahoo.com.br alexandre.santos@cas.ifmt.edu.br ====================================================================== _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.