[R-br] Problema para mudar projeção em um raster

ASANTOS alexandresantosbr em yahoo.com.br
Quarta Julho 17 15:00:58 BRT 2013


Boa tarde Pessoal,

Criei um raster e referencie em UTM sendo:

require(raster)
require(rgdal)
require(sp)
#Criando raster
r <- raster(nc=1000, nr=1000)
set.seed(20130622)
stackIm1 <- stack(lapply(1, function(x) setValues(r, 
round(runif(ncell(r))* 255))))## Simulation red band
#
# Definindo sistema de projeção
r.geo <- CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m 
+no_defs")    # geographical datum WGS84
proj4string(stackIm1) <- r.geo
#

mas quando tento converter para latlong com spTransform:

stackIm1.ll<-spTransform(stackIm1, CRS("+proj=longlat+datum=WGS84"))


Erro em spTransform(stackIm1, CRS("+proj=longlat+datum=WGS84")) :
   erro na avaliação do argumento 'CRSobj' na seleção do método para a 
função 'spTransform': 'Erro em CRS("+proj=longlat+datum=WGS84") : 
unknown projection id

Obtenho este erro e não consigo achar a solução, alguém poderia dar um help,

Obrigado

-- 
======================================================================
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 em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================



Mais detalhes sobre a lista de discussão R-br