<div dir="ltr">Complementando a resposta anterior (abaixo), o arquivo topodata do CMR não corresponde à área de interesse definida. A área coincidente é dada pelo arquivo 19S435SN.zip.<div><br></div><div>### </div><div>require(raster)<br>
</div><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>dem <- raster('19S435SN.tif')</div><div>###</div><div><br></div><div>Observe as saídas transcritas abaixo:</div><div><div>> extent(SP.geo)</div><div>class       : Extent </div><div>xmin        : -42.40624 </div>
<div>xmax        : -42.37724 </div><div>ymin        : -19.51573 </div><div>ymax        : -19.48824 </div><div>> dem <- raster('19S435SN.tif')</div><div>> extent(dem)</div><div>class       : Extent </div><div>
xmin        : -43.50024 </div><div>xmax        : -42 </div><div>ymin        : -20 </div><div>ymax        : -19 </div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/7/24 Eder Comunello <span dir="ltr"><<a href="mailto:ecomunel@gmail.com" target="_blank">ecomunel@gmail.com</a>></span><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"><div dir="ltr"><div><div>Prezado Alexandre,</div><div><br>
</div><div>Quanto ao erro de spTransform(), faltou dar espaço antes do argumento '+south'.</div><div><br></div><div>### INICIO</div><div><br></div><div>### bnds fornecido no CMS<br>
</div><div><div>bnds <- structure(c(772259.1, 775259.1, 775259.1, 772259.1, 772259.1, </div><div>            7843046, 7843046, 7840046, 7840046, 7843046), .Dim = c(5L, 2L</div><div>            ), .Dimnames = list(NULL, c("x", "y")))</div>

</div><div><br></div><div>### CRS UTM</div><div>SP <- SpatialPolygons(list(Polygons(list(Polygon(bnds)), "1")))</div><div>proj4string(SP) = CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs") ## Projeção</div>

<div><br></div><div>### CRS GEO</div><div>SPgeo<- spTransform(SP, CRS("+proj=longlat +datum=WGS84")) ### aqui não é necessário +south</div><div>SPgeo@bbox</div><div><br></div><div>### FIM</div><div><br></div>

<div>Com esse código você obterá:</div><div><br></div><div>> SPgeo@bbox</div><div>        min       max</div><div>x -42.40624 -42.37724</div><div>y -19.51573 -19.48824</div></div><div><br></div><div>At.te,</div><span class=""><font color="#888888"><div>
<br>
</div><div><br></div>-- <br><div dir="ltr">Éder Comunello <<a href="mailto:comunello.eder@gmail.com" target="_blank">comunello.eder@gmail.com</a>> <br>Dourados, MS - [22 16.5'S, 54 49'W]<br></div>
</font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div dir="ltr">Éder Comunello <<a href="mailto:comunello.eder@gmail.com" target="_blank">c</a><a href="mailto:omunello.eder@gmail.com" target="_blank">omunello.eder@gmail.com</a>> <br>
Dourados, MS - [22 16.5'S, 54 49'W]<br></div>
</div></div></div></div>