<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Boa tarde Éder,<br>
<br>
Funcionou perfeitamente e deixou os argumentos muito mais
lógicos, comparado ao que vinha fazendo. Meu único problema agora
são os vértices irregulares, pois quando tento recortar o mesmo
contorno utilizado (SP.geo) em uma imagem de satélite landsat 5 em
utm, por exemplo, o tamanho da minha imagem nunca fica o mesmo que o
do objeto alt.utm, sendo que se não tenho o mesmo número de nrow e
ncol não posso comparar diretamente as duas imagens com a abordagem
que emprego, de qualquer maneira muito obrigado, ajudou bastante,<br>
<br>
Abraço,<br>
<br>
<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>
<br>
<br>
<br>
<div class="moz-cite-prefix">Em 27/07/2013 09:07, Eder Comunello
escreveu:<br>
</div>
<blockquote
cite="mid:CABmC8g=kKnZ0t1DAg4CoNVmGHqJM+oZmyFBysznnm2CFR64bMA@mail.gmail.com"
type="cite">
<div dir="ltr">
<div>Alexandre e colegas, bom dia!</div>
<div><br>
</div>
<div>Meu script está um pouco diferente, mas o que vai valer
mesmo é o último bloco (reprojetar RasteLayer). Espero que
ajude.</div>
<div><br>
</div>
<div>
Além disso, tenho um bookmark pra indicar: </div>
<div>Spatial data in R: using R as a GIS (old version) <<a
moz-do-not-send="true"
href="https://sites.google.com/site/rodriguezsanchezf/news/usingrasagis">https://sites.google.com/site/rodriguezsanchezf/news/usingrasagis</a>></div>
<div><br>
</div>
<div>Vi que saiu uma nova versão em <<a
moz-do-not-send="true"
href="https://sites.google.com/site/rodriguezsanchezf/resources/rgis">https://sites.google.com/site/rodriguezsanchezf/resources/rgis</a>></div>
<div><br>
</div>
<div>
Outro dia, alguém perguntava pela geração de mapas temáticos
no R! Esses links podem ajudar bastante!</div>
<div><br>
</div>
<div>At.te,</div>
<div><br>
</div>
<div>### <BEGIN> ###</div>
<div>require(sp); require(rgdal); require(raster)</div>
<div><br>
</div>
<div>### Região de interesse</div>
<div>pc <- cbind(X1=773759.1,Y1=7841546)</div>
<div>dx <- 1500*c(-1, 1, 1,-1,-1); dx</div>
<div>dy <- 1500*c( 1, 1,-1,-1, 1); dy</div>
<div>reg <- NULL; reg <-
rbind(cbind(pc[,1]+dx,pc[,2]+dy)); reg</div>
<div>plot(reg, asp=1); points(pc, col=2)</div>
<div><br>
</div>
<div>### CRS UTM</div>
<div>SP.utm <-
SpatialPolygons(list(Polygons(list(Polygon(reg)), "1")))</div>
<div>proj4string(SP.utm) = CRS("+proj=utm +zone=23 +south
+datum=WGS84 +units=m +no_defs")</div>
<div>SP.utm@bbox; SP.utm@proj4string</div>
<div>plot(SP.utm, col=3, axes=T)</div>
<div><br>
</div>
<div>### CRS GEO (GeoTiff)</div>
<div>SP.geo<- spTransform(SP.utm, CRS("+proj=longlat
+datum=WGS84"))</div>
<div>SP.geo@bbox; SP.geo@proj4string</div>
<div>plot(SP.geo, col=5, axes=T)</div>
<div><br>
</div>
<div>### TOPODATA</div>
<div>#getwd()</div>
<div>#url=("<a moz-do-not-send="true"
href="http://www.dsr.inpe.br/topodata/data/geotiff/19S435SN.zip">http://www.dsr.inpe.br/topodata/data/geotiff/19S435SN.zip</a>")</div>
<div>#download.file(url, destfile = "19S435SN.zip")</div>
<div>#unzip("19S435SN.zip")</div>
<div>dem <- raster('19S435SN.tif') ### Abre o slope TOPODATA</div>
<div><br>
</div>
<div>### Recorte TOPODATA</div>
<div>dem.crop <- crop(dem, extent(SP.geo), snap='out')</div>
<div>image(dem.crop, asp=1)</div>
<div><br>
</div>
<div># Reprojetar RasterLayer (S4)</div>
<div>projection(dem.crop) <- CRS("+proj=longlat
+datum=WGS84")</div>
<div>alt.utm <- projectRaster(dem.crop, crs="+proj=utm
+zone=23 +south +datum=WGS84 +units=m +no_defs")</div>
<div>image(alt.utm, asp=1)</div>
<div><br>
</div>
<div>### <END> ###</div>
<div class="gmail_extra">
<br>
<br>
<div class="gmail_quote">Em 25 de julho de 2013 22:32, ASANTOS
<span dir="ltr"><<a moz-do-not-send="true"
href="mailto:alexandresantosbr@yahoo.com.br"
target="_blank">alexandresantosbr@yahoo.com.br</a>></span>
escreveu:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px
0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Boa
noite Pessoal,<br>
<br>
Consegui resolver um problema e consegui outro,
transformei as coordenadas latlong de um geotiff do
Topodata em utm usando spTransform(). Como não consegui
transformar diretamente um SpatialPixelsDataFrame do
topodata apenas colocando o novo CRS e em vários posts não
vi nenhuma solução para isso, tentei fazer da seguinte
maneira, transformei apenas o @coords do topodata de geo
para utm e depois tentei inserir com replace() os novos
valores em utm dentro do SpatialPixelsDataFrame e ocorreu
o seguinte erro:<br>
<br>
Erro em replace(dem.sp@coords[1:10026, 1:2],
dem.sp@coords, Utm@coords[1:10026, :<br>
valores negativos não são permitidos na subscrição de
matriz<br>
<br>
</blockquote>
</div>
</div>
</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">
</pre>
</body>
</html>