Thiago,
Explicando melhor teoricamente, pois um CRM implicaria em uma
imagem de 50MB, tenho:
1) Uma imagem do landsat 5 com um quadrante de interesse;
2)Um contorno em shapefile de uma área de 200 hectares.
Quero utilizar apenas os pixels da imagem contidos no interior da
área de 200 hectares, para tanto extend() não resolve pois utiliza
apenas utiliza coordenadas max e min e não considera a forma (contorno)
do meu polígono para seleção dos pixels contidos nele.
Não consegui realizar o que desejo com a rotina sugerida. Mascara,
usei como um termo que implica em área útil da imagem, não quero o
quadrante todo oferecido pela imagem, mas aquilo que esta contido no meu
shapefile.
Obrigado pela atenção,
Alexandre
Em 04/04/2013 21:03, Thiago V. dos Santos escreveu:
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.
