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.
###require(raster)
download.file(url, destfile = "19S435SN.zip")dem <- raster('19S435SN.tif')###
Observe as saídas transcritas abaixo:> extent(SP.geo)class : Extentxmin : -42.40624xmax : -42.37724ymin : -19.51573ymax : -19.48824> dem <- raster('19S435SN.tif')> extent(dem)class : Extentxmin : -43.50024xmax : -42ymin : -20ymax : -19
2013/7/24 Eder Comunello <ecomunel@gmail.com>
Prezado Alexandre,
Quanto ao erro de spTransform(), faltou dar espaço antes do argumento '+south'.
### INICIO
### bnds fornecido no CMS
bnds <- structure(c(772259.1, 775259.1, 775259.1, 772259.1, 772259.1,7843046, 7843046, 7840046, 7840046, 7843046), .Dim = c(5L, 2L), .Dimnames = list(NULL, c("x", "y")))
### CRS UTMSP <- SpatialPolygons(list(Polygons(list(Polygon(bnds)), "1")))proj4string(SP) = CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs") ## Projeção
### CRS GEOSPgeo<- spTransform(SP, CRS("+proj=longlat +datum=WGS84")) ### aqui não é necessário +southSPgeo@bbox
### FIM
Com esse código você obterá:
> SPgeo@bboxmin maxx -42.40624 -42.37724y -19.51573 -19.48824
At.te,
--
--
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- ====================================================================== 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) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================