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

[Thiago V. dos Santos]

Alexandre, 

Não sou um especialista em cartografia ou geodésia, então posso estar enganado, mas acho que o que o Elias quis dizer é que você está criando um raster com coordenadas lat/lon e atribuindo a esse raster uma projeção utm.

O meu palpilte é de que isso está causando o erro na sua reprojeção. Veja as propriedades do objeto stackIm1:

> stackIm1
class       : RasterStack 
dimensions  : 500, 250, 125000, 1  (nrow, ncol, ncell, nlayers)
resolution  : 0.36, 0.18  (x, y)
extent      : 0, 90, 0, 90  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 
names       : layer 
min values  :     0 
max values  :   255 

As coordenadas que eu esperaria ver em projeção UTM são do tipo mapa, por exemplo 213539.91 e 6566533.33 (para uma ponto no RS - zona 22S).

Nos seus dados reais as coordenadas são de fato UTM? Se você não conseguir postar um CMR com o seu dado real, tente criar uma amostra com coordenadas UTM que talvez funcione.

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: Thursday, July 18, 2013 10:39 AM
Subject: Re: [R-br] Problema para mudar projeção em um raster
 

Elias,

     Obrigado, mas cortei um pedaço das bordas e deixei os valores 
positivos e mesmo assim:

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


Segue CRM:

require(raster)
require(rgdal)
require(sp)
#Criando raster
r <- raster(nc=1000, nr=1000)
r[] <- 1:ncell(r)
e <- extent(0, 90, 0, 90)
stackIm1.crop <- crop(r, e)
set.seed(20130622)
  stackIm1 <- stack(lapply(1, function(x) setValues(stackIm1.crop,
+ round(runif(ncell(stackIm1.crop))* 255))))##
#
# 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)

Alexandre


Em 18/07/2013 09:15, Elias Krainski escreveu:
> Quando vc faz
>    r <- raster(nc=1000, nr=1000)
> cria-se um raster com bordas:
>    bbox(r)
> totalmente incompatíveis com "+proj=utm +zone=23 +south +datum=WGS84 
> +units=m"
> _______________________________________________
> 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
======================================================================

_______________________________________________
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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130718/9d98dc02/attachment.html>