<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 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 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 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 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>