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

[ASANTOS]

Obrigado Thiago,
mas nenhumas das duas sugestões funcionou não, 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
#pelo project raster
 > newproj <- "+proj=longlat +datum=WGS84"
  #Nova projeção
stackIm1ll <- projectRaster(stackIm1, crs=newproj)
Erro em if (value[1] != nrow(x) | value[2] != ncol(x)) { :
   valor ausente onde TRUE/FALSE necessário
Além disso: Mensagens de aviso perdidas:
In `dim<-`(`*tmp*`, value = c(nr, nc)) : NAs introduzidos por coerção


Li alguns posts, mas sempre o caminho é bem árduo, busco uma solução 
simples, de qualquer maneira obrigado,



Alexandre





Em 17/07/2013 16:04, Thiago V. dos Santos escreveu:
> Alexandre,
>
> Parece ser um typo. Experimente adicionar um espaço entre proj e 
> datum: '+proj=longlat +datum=WGS84'.
>
> Uma abordagem alternativa é usar projectRaster, do próprio pacote raster.
> Saudações,
> --
> Thiago V. dos Santos
> PhD student
> Land and Atmospheric Science
> University of Minnesota
> http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm
> Phone: (612) 323 9898
> ------------------------------------------------------------------------
> *From:* ASANTOS <alexandresantosbr em yahoo.com.br>
> *To:* r-br em listas.c3sl.ufpr.br
> *Sent:* Wednesday, July 17, 2013 1:00 PM
> *Subject:* [R-br] Problema para mudar projeção em um raster
>
> 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 
> <mailto:alexandresantosbr em yahoo.com.br>
> alexandre.santos em cas.ifmt.edu.br <mailto:alexandre.santos em cas.ifmt.edu.br>
> Lattes: http://lattes.cnpq.br/1360403201088680
> ======================================================================
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br <mailto:R-br em 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.
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em 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 em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130718/cc12d880/attachment.html>