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

[Thiago V. dos Santos]

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 <alexandresantosbr em yahoo.com.br>
To: r-br em listas.c3sl.ufpr.br
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<alexandresantosbr em yahoo.com.br>
> To:r-br em 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 em yahoo.com.br
>          alexandre.santos em cas.ifmt.edu.br
> ======================================================================
>
> _______________________________________________
> R-br mailing list
> R-br em 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.
>
>

-- 
======================================================================
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                    
======================================================================

_______________________________________________
R-br mailing list
R-br em 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.