[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