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 <[email protected]>
*To:* [email protected]
*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:[email protected]
<mailto:[email protected]>
[email protected] <mailto:[email protected]>
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================
_______________________________________________
R-br mailing list
[email protected] <mailto:[email protected]>
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
[email protected]
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:[email protected]
[email protected]
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================
_______________________________________________
R-br mailing list
[email protected]
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.
