<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Elias,<br>
<br>
Estou tentando selecionar uma pequena área de interesse na
imagem para não demorar muito, mas agora quando tento o crop()
tenho: Erro em function (classes, fdef, mtable) : <br>
unable to find an inherited method for function "crop", for
signature "SpatialGridDataFrame"<br>
<br>
Na verdade criei um contorno SpatialPolygons e não
SpatialGridDataFrame, mas a função não aceita segue CRM:<br>
<br>
# Criar uma área de menor dimensão que a imagem inteira que abarque
a região de interesse<br>
xcc<-773759.1<br>
ycc<-7841546<br>
p.central<-cbind(xcc,ycc)<br>
<br>
###Criando os vértices da área<br>
coordV <-
rbind(coordV,cbind(p.central[,1]+c(-1500,1500,1500,-1500,-1500),p.central[,2]+c(1500,1500,-1500,-1500,1500)))<br>
coordV<br>
coordV2=coordV[!is.na(coordV[,1]),]<br>
coordV2<br>
plot(coordV2[21:25,1],coordV2[21:25,2])<br>
points(p.central[,1],p.central[,2], col="red")<br>
#<br>
<br>
# Cria um polígono com o contorno definido<br>
bnds <- cbind(x=c(coordV2[21:25,1]), y=c(coordV2[21:25,2]))<br>
<br>
# CRS UTM<br>
SP <- SpatialPolygons(list(Polygons(list(Polygon(bnds)), "1")))<br>
proj4string(SP) = CRS("+proj=utm +zone=23+south +datum=WGS84
+units=m +no_defs") ## Projeção<br>
<br>
# CRS em lat long para o polígomno ficar na mesma projeção que o
GeoTif<br>
SPlatlong<- spTransform(SP, CRS("+proj=longlat +datum=WGS84"))<br>
#<br>
<br>
demlatlog<-readGDAL('19S435HN.tif')### Abre o slope TOPODATA - <font
face="Courier New, Courier, monospace">geotif do outro e-mail</font><br>
<br>
### atribui a projeção longlat e datum WGS83 ao raster.<br>
proj4string(demlatlog) <- CRS("+proj=longlat +datum=WGS84")<br>
<br>
#Cortar uma área menor de interesse<br>
demlatlog.crop <- crop(demlatlog, extent(SPlatlong), snap='out')<br>
#<br>
<br>
Obrigado,<br>
<br>
Alexandre<br>
<br>
<div class="moz-cite-prefix">Em 18/07/2013 19:01, Elias Krainski
escreveu:<br>
</div>
<blockquote cite="mid:51E873C7.6020402@yahoo.com.br" type="cite">
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
<div><font face="Courier New, Courier, monospace">Caro Alexandre,
<br>
<br>
Vc esta usando aquele geotif do outro e-mail? <br>
<br>
### obtem o arquivo <br>
url=("<a moz-do-not-send="true"
href="http://www.webmapit.com.br/inpe/topodata/download/?ms=137398652719842&fbn=21S465&fsx=3">http://www.webmapit.com.br/inpe/topodata/download/?ms=137398652719842&fbn=21S465&fsx=3</a>")
<br>
download.file(url, destfile = "21S465HN.zip") <br>
### descompacta<br>
system("unzip 21S465HN.zip") <br>
</font><br>
<font face="Courier New, Courier, monospace">### le o raster<br>
require(rgdal) <br>
r <- readGDAL("21S465SN.tif")<br>
<br>
### atribui a projeção longlat e datum WGS83 ao raster. <br>
proj4string(r) <- CRS("+proj=longlat +datum=WGS84")<br>
</font><font face="Courier New, Courier, monospace"><br>
### transformar com é muito ineficiente <br>
### r2 <- spTransform(r, CRS("+proj=utm +zone=23 +south"))<br>
</font><br>
<font face="Courier New, Courier, monospace">### com
projectRaster(), se vc tivesse usando <br>
### r <- raster("21S465SN.tiff")<br>
### r2 <- projectRaster(r, crs=</font><font face="Courier
New, Courier, monospace"><font face="Courier New, Courier,
monospace">CRS("+proj=utm +zone=23 +south"))</font> <br>
### também é muito lento <br>
<br>
### a lentidão é devida ao fato de que extrai-se as <br>
### </font><font face="Courier New, Courier, monospace"><font
face="Courier New, Courier, monospace">coordenadas </font>de
cada pixel e converte-se cada uma...</font><font face="Courier
New, Courier, monospace"> </font><font face="Courier New,
Courier, monospace"><font face="Courier New, Courier,
monospace"><font face="Courier New, Courier, monospace"><font
face="Courier New, Courier, monospace"><font
face="Courier New, Courier, monospace"><font
face="Courier New, Courier, monospace"><font
face="Courier New, Courier, monospace"><font
face="Courier New, Courier, monospace"><font
face="Courier New, Courier, monospace"><font
face="Courier New, Courier, monospace"><font
face="Courier New, Courier, monospace"><font
face="Courier New, Courier, monospace"><br>
<br>
</font></font></font></font></font></font></font></font></font></font></font></font>Att.<br>
Elias<br>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">--
======================================================================
Alexandre dos Santos
Proteção 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)
<a class="moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a>
<a class="moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>
Lattes: <a class="moz-txt-link-freetext" href="http://lattes.cnpq.br/1360403201088680">http://lattes.cnpq.br/1360403201088680</a>
======================================================================
</pre>
</body>
</html>