<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Obrigado Mauro,<br>
<br>
Resolvi o erro na rotina e mesmo assim o spTransform() não esta
convertendo de utm na zona 23K para latlong no hemisfério sul, segue
CRM completo ate o recorte da área de interesse no geoTiff do
topodata:<br>
<br>
require(sp)<br>
require(rgdal)<br>
require(raster)<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<-NULL<br>
coordV
<-rbind(coordV,cbind(p.central[,1]+c(-1500,1500,1500,-1500,-1500),p.central[,2]+c(1500,1500,-1500,-1500,1500)))<br>
plot(coordV[,1],coordV[,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(coordV[,1]), y=c(coordV[,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 +south
+datum=WGS84"))<br>
#<br>
#<br>
#Obtenção do objeto declividade do TOPODATA<br>
url=(<a class="moz-txt-link-rfc2396E" href="http://www.dsr.inpe.br/topodata/data/geotiff/19S42_SN.zip">"http://www.dsr.inpe.br/topodata/data/geotiff/19S42_SN.zip"</a>)<br>
download.file(url, destfile = "19S42_SN.zip")<br>
### descompacta<br>
system("unzip 19S42_SN.zip")<br>
demlatlog<-raster('19S42_SN.tif')### Abre o slope TOPODATA<br>
#<br>
#Cortar a região de interesse<br>
demlatlog.crop <- crop(demlatlog, extent(SPlatlong), snap='out')<br>
#<br>
<br>
Obrigado,<br>
<br>
Alexandre<br>
<br>
<br>
<br>
<br>
<br>
<div class="moz-cite-prefix"> Em 22/07/2013 22:09, Mauro Sznelwar
escreveu:<br>
</div>
<blockquote cite="mid:BB73DDED0A3340EBBF6F755FFF17601C@ACER4520"
type="cite">
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
<meta name="GENERATOR" content="MSHTML 8.00.6001.23507">
<style></style>
<div><font color="#0000ff" face="Arial">Tentei rodar o seu CRM e
deu problema aqui</font></div>
<div><font size="2" face="Arial"><font size="3">> bnds <-
cbind(x=c(coordV2[21:25,1]), y=c(coordV2[21:25,2]))<br>
Error in coordV2[21:25, 1] : subscript out of bounds<br>
> <br>
<br>
</font>Segue CRM:<br>
<br>
require(sp)<br>
require(rgdal)<br>
# Criar uma área de menor dimensão que a imagem inteira que
abarque a <br>
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 <- <br>
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>
#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 <br>
+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 +south
+datum=WGS84"))<br>
#<br>
<br>
<br>
Obrigado,<br>
<br>
Alexandre<br>
<br>
<br>
<br>
<br>
</font></div>
</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>