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/inde... Phone: (612) 323 9898 ------------------------------------------------------------------------ *From:* ASANTOS <alexandresantosbr@yahoo.com.br> *To:* r-br@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@yahoo.com.br <mailto:alexandresantosbr@yahoo.com.br> alexandre.santos@cas.ifmt.edu.br <mailto:alexandre.santos@cas.ifmt.edu.br> Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br <mailto:R-br@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@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@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================