<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div><span>Mas você recebe algum erro ou o resultado é diferente do esperado?</span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal;"><span><br></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: 'times new roman', 'new york', times, serif; background-color: transparent; font-style: normal;"><span>O melhor mesmo seria fornecer seus dados para tentarmos reproduzir o problema. Não dá pra colocar no dropbox? 50Mb não é nada tão grande assim...</span></div><div></div><div> </div><div>Saudações,<br>--<br>Thiago V. dos Santos<br>PhD student<br>Land and Atmospheric Science<br>University of
Minnesota<br>http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm<br>Phone: (612) 323 9898<br></div> <div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;"> <div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;"> <div dir="ltr"> <font size="2" face="Arial"> <hr size="1"> <b><span style="font-weight:bold;">From:</span></b> ASANTOS <alexandresantosbr@yahoo.com.br><br> <b><span style="font-weight: bold;">To:</span></b> r-br@listas.c3sl.ufpr.br; Thiago V. dos Santos <thi_veloso@yahoo.com.br> <br> <b><span style="font-weight: bold;">Sent:</span></b> Thursday, April 4, 2013 9:23 PM<br> <b><span style="font-weight: bold;">Subject:</span></b> Re: [R-br] Criando uma mascara com shapefile sobre raster<br> </font> </div> <br>Thiago,<br><br> Explicando melhor teoricamente, pois um CRM implicaria em uma <br>imagem de 50MB, tenho:<br><br>1)
Uma imagem do landsat 5 com um quadrante de interesse;<br>2)Um contorno em shapefile de uma área de 200 hectares.<br><br> Quero utilizar apenas os pixels da imagem contidos no interior da <br>área de 200 hectares, para tanto extend() não resolve pois utiliza <br>apenas utiliza coordenadas max e min e não considera a forma (contorno) <br>do meu polígono para seleção dos pixels contidos nele.<br><br> Não consegui realizar o que desejo com a rotina sugerida. Mascara, <br>usei como um termo que implica em área útil da imagem, não quero o <br>quadrante todo oferecido pela imagem, mas aquilo que esta contido no meu <br>shapefile.<br><br>Obrigado pela atenção,<br><br>Alexandre<br><br><br><br>Em 04/04/2013 21:03, Thiago V. dos Santos escreveu:<br>> Alexandre,<br>><br>> O que a área do raster tem a ver com o tipo de função usada para carregar o shapefile?<br>><br>> 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?<br>><br>> Saudações,<br>> --<br>> Thiago V. dos Santos<br>> PhD student<br>> Land and Atmospheric Science<br>> University of Minnesota<br>> <a href="http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm" target="_blank">http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm</a><br>> Phone: (612) 323 9898<br>><br>><br>> ----- Original Message -----<br>> From: ASANTOS <<a ymailto="mailto:alexandresantosbr@yahoo.com.br" href="mailto:alexandresantosbr@yahoo.com.br">alexandresantosbr@yahoo.com.br</a>><br>> To: <a ymailto="mailto:r-br@listas.c3sl.ufpr.br" href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>> Cc:<br>> Sent: Thursday, April 4, 2013 5:14 PM<br>> Subject: Re: [R-br] Criando uma mascara
com shapefile sobre raster<br>><br>> Thiago,<br>><br>> Não deu certo não, transformei em<br>> SpatialPolygons(list(Polygons(list(Polygon(contorno2)),"contorno2"))),<br>> prefiro continuar tentando com readShapeLines("Catas_Altas.shp",<br>> proj4string=CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m<br>> +no_defs")), porque consigo visualizar a área sobre o raster, achei<br>> alguns posts mais usando o GRASS, mas gostaria de fazer tudo só no R.<br>><br>> Obrigado,<br>><br>> Alexandre<br>><br>><br>> Em 04/04/2013 12:18, Thiago V. dos Santos escreveu:<br>>> Alexandre,<br>>><br>>> a) Tente ler o seu shapefile com a funcao readShapePoly.<br>>><br>>> b) Se nao der certo, sugiro uma combinacao das funcoes crop e raster. Nesse caso, experimente rodar o seu codigo com as seguintes mudancas:<br>>><br>>>
require("raster")<br>>> A <- stack("LANDSAT_5_TM_20100506_217_074_L2_BAND7.tif")<br>>> contorno<- readShapePoly("Catas_Altas.shp")<br>>><br>>> A.crop <- crop(A, extent(contorno), snap='out')<br>>><br>>> contorno.na<-setValues(A.crop, NA)<br>>> contorno.r<-rasterize(contorno, contorno.na) ### isso pode demorar, dependendo do seu shapefile<br>>> C <- mask(x=A.crop, mask=contorno.r)<br>>><br>>><br>>> Saudações,<br>>><br>>> --<br>>> Thiago V. dos Santos<br>>> PhD student<br>>> Land and Atmospheric Science<br>>> University of Minnesota<br>>> <a href="http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm" target="_blank">http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm</a><br>>> Phone: (612) 323 9898<br>>><br>>><br>>> ----- Original Message
-----<br>>> From: ASANTOS<<a ymailto="mailto:alexandresantosbr@yahoo.com.br" href="mailto:alexandresantosbr@yahoo.com.br">alexandresantosbr@yahoo.com.br</a>><br>>> To:<a ymailto="mailto:r-br@listas.c3sl.ufpr.br" href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>>> Cc:<br>>> Sent: Thursday, April 4, 2013 10:31 AM<br>>> Subject: [R-br] Criando uma mascara com shapefile sobre raster<br>>><br>>> Boa tarde pessoal,<br>>><br>>> 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:<br>>><br>>> require("raster")<br>>> A <-
stack(c("LANDSAT_5_TM_20100506_217_074_L2_BAND7.tif"))## Imagem<br>>> plotRGB(A) ## Plota a imagem<br>>> contorno<- readShapeLines("Catas_Altas.shp") ## Contorno da área<br>>> lines(contorno, col="red")Plota o contorno<br>>> C <- crop(A,contorno) ## Recorta o que esta contido no contorno na imagem<br>>> Erro em .local(x, y, ...) :<br>>> nenhum slot de nome "legend" para esse objeto de classe "RasterStack"<br>>><br>>><br>>> -- ======================================================================<br>>> Alexandre dos Santos<br>>> Proteção Florestal<br>>> Coordenador do curso Técnico em Florestas<br>>> Vice Coordenador do curso de Engenharia Florestal<br>>> IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso<br>>> Campus Cáceres<br>>> Caixa Postal 244<br>>> Avenida dos Ramires, s/n<br>>>
Bairro: Distrito Industrial<br>>> Cáceres - MT CEP: 78.200-000<br>>> Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)<br>>> e-mails:<a ymailto="mailto:alexandresantosbr@yahoo.com.br" href="mailto:alexandresantosbr@yahoo.com.br">alexandresantosbr@yahoo.com.br</a><br>>> <a ymailto="mailto:alexandre.santos@cas.ifmt.edu.br" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a><br>>> ======================================================================<br>>><br>>> _______________________________________________<br>>> R-br mailing list<br>>> <a ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>>> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br"
target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>>> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br>>><br>>><br><br>-- <br>======================================================================<br>Alexandre dos Santos<br>Proteção Florestal<br>Coordenador do curso Técnico em Florestas<br>Vice Coordenador do curso de Engenharia Florestal<br>IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso<br>Campus Cáceres<br>Caixa Postal 244<br>Avenida dos Ramires, s/n<br>Bairro: Distrito Industrial<br>Cáceres - MT CEP: 78.200-000<br>Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)<br>e-mails:<a ymailto="mailto:alexandresantosbr@yahoo.com.br"
href="mailto:alexandresantosbr@yahoo.com.br">alexandresantosbr@yahoo.com.br</a><br> <a ymailto="mailto:alexandre.santos@cas.ifmt.edu.br" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a><br>======================================================================<br><br><br><br> </div> </div> </div></body></html>