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 <
alexandresantosbr@yahoo.com.br>
> To:
r-br@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@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.
>>
>>
--
======================================================================
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
======================================================================