Alexandre, O que a área do raster tem a ver com o tipo de função usada para carregar o shapefile?
Você não consegui recortar o raster com o shapefile usando os comandos que eu eu sugeri? O que exatamente você quis dizer com criar uma máscara na pergunta inicial? Saudações, -- Thiago V. dos Santos PhD student Land and Atmospheric Science University of Minnesota http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm Phone: (612) 323 9898 ----- Original Message ----- From: ASANTOS <[email protected]> To: [email protected] Cc: Sent: Thursday, April 4, 2013 5:14 PM Subject: Re: [R-br] Criando uma mascara com shapefile sobre raster Thiago, Não deu certo não, transformei em SpatialPolygons(list(Polygons(list(Polygon(contorno2)),"contorno2"))), prefiro continuar tentando com readShapeLines("Catas_Altas.shp", proj4string=CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs")), porque consigo visualizar a área sobre o raster, achei alguns posts mais usando o GRASS, mas gostaria de fazer tudo só no R. Obrigado, Alexandre Em 04/04/2013 12:18, Thiago V. dos Santos escreveu: > 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/index.htm > Phone: (612) 323 9898 > > > ----- Original Message ----- > From: ASANTOS<[email protected]> > To:[email protected] > 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:[email protected] > [email protected] > ====================================================================== > > _______________________________________________ > R-br mailing list > [email protected] > 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. > > -- ====================================================================== 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:[email protected] [email protected] ====================================================================== _______________________________________________ R-br mailing list [email protected] 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. _______________________________________________ R-br mailing list [email protected] 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.
