[R-br] ?Problema para mudar projeção em um raster
Eder Comunello
ecomunel em gmail.com
Quarta Julho 24 17:29:26 BRT 2013
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 UTM
SP <- SpatialPolygons(list(Polygons(list(Polygon(bnds)), "1")))
proj4string(SP) = CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m
+no_defs") ## Projeção
### CRS GEO
SPgeo<- spTransform(SP, CRS("+proj=longlat +datum=WGS84")) ### aqui não é
necessário +south
SPgeo em bbox
### FIM
Com esse código você obterá:
> SPgeo em bbox
min max
x -42.40624 -42.37724
y -19.51573 -19.48824
At.te,
--
Éder Comunello <comunello.eder em gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130724/2771d4fd/attachment.html>
Mais detalhes sobre a lista de discussão R-br