<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;"><span>Alexandre, </span></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 16px; color: rgb(0, 0, 0); background-color: transparent; font-style: normal;"><span><br></span></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 16px; color: rgb(0, 0, 0); background-color: transparent; font-style: normal;"><span>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 </span><span style="background-color: transparent;">e atribuindo a esse raster uma projeção utm.</span></div><div style="font-family: 'times new roman', 'new york', times, serif;
font-size: 16px; color: rgb(0, 0, 0); background-color: transparent; font-style: normal;"><span style="background-color: transparent;"><br></span></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 16px; color: rgb(0, 0, 0); background-color: transparent; font-style: normal;">O meu palpilte é de que isso está causando o erro na sua reprojeção. Veja as propriedades do objeto stackIm1:</div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 16px; color: rgb(0, 0, 0); background-color: transparent; font-style: normal;"><br></div><div style="background-color: transparent;">> stackIm1</div><div style="background-color: transparent;">class : RasterStack </div><div style="background-color: transparent;">dimensions : 500, 250, 125000, 1 (nrow, ncol, ncell, nlayers)</div><div style="background-color: transparent;">resolution : 0.36, 0.18 (x,
y)</div><div style="background-color: transparent;">extent : 0, 90, 0, 90 (xmin, xmax, ymin, ymax)</div><div style="background-color: transparent;">coord. ref. : +proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 </div><div style="background-color: transparent;">names : layer </div><div style="background-color: transparent;">min values : 0 </div><div style="background-color: transparent;">max values : 255 </div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;"><br></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;"></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;">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).</div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;"><br></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;">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.</div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;"><br></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;">Saudações,<br>--<br>Thiago V. dos Santos<br>PhD student<br>Land and Atmospheric Science<br>University of Minnesota<br>http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm<br>Phone: (612) 323 9898<br></div> <div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;"> <div style="font-family: 'times new roman', 'new york', times,
serif; font-size: 12pt;"> <div dir="ltr"> <hr size="1"> <font size="2" face="Arial"> <b><span style="font-weight:bold;">From:</span></b> ASANTOS <alexandresantosbr@yahoo.com.br><br> <b><span style="font-weight: bold;">To:</span></b> r-br@listas.c3sl.ufpr.br <br> <b><span style="font-weight: bold;">Sent:</span></b> Thursday, July 18, 2013 10:39 AM<br> <b><span style="font-weight: bold;">Subject:</span></b> Re: [R-br]
Problema para mudar projeção em um raster<br> </font> </div> <div class="y_msg_container"><br>Elias,<br><br> Obrigado, mas cortei um pedaço das bordas e deixei os valores <br>positivos e mesmo assim:<br><br>Erro em if (value[1] != nrow(x) | value[2] != ncol(x)) { :<br> valor ausente onde TRUE/FALSE necessário<br>Além disso: Mensagens de aviso perdidas:<br>In `dim<-`(`*tmp*`, value = c(nr, nc)) : NAs introduzidos por coerção<br><br><br>Segue CRM:<br><br>require(raster)<br>require(rgdal)<br>require(sp)<br>#Criando raster<br>r <- raster(nc=1000, nr=1000)<br>r[] <- 1:ncell(r)<br>e <- extent(0, 90, 0, 90)<br>stackIm1.crop <- crop(r, e)<br>set.seed(20130622)<br> stackIm1 <- stack(lapply(1, function(x) setValues(stackIm1.crop,<br>+ round(runif(ncell(stackIm1.crop))* 255))))##<br>#<br># Definindo sistema de projeção<br>r.geo <- CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m<br>+
+no_defs") # geographical datum WGS84<br> proj4string(stackIm1) <- r.geo<br>#pelo project raster<br>newproj <- "+proj=longlat +datum=WGS84"<br> #Nova projeção<br>stackIm1ll <- projectRaster(stackIm1, crs=newproj)<br><br>Alexandre<br><br><br>Em 18/07/2013 09:15, Elias Krainski escreveu:<br>> Quando vc faz<br>> r <- raster(nc=1000, nr=1000)<br>> cria-se um raster com bordas:<br>> bbox(r)<br>> totalmente incompatíveis com "+proj=utm +zone=23 +south +datum=WGS84 <br>> +units=m"<br>> _______________________________________________<br>> R-br mailing list<br>> <a ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>> Leia o guia de postagem (<a
href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça <br>> código mínimo reproduzível.<br>><br><br>-- <br>======================================================================<br>Alexandre dos Santos<br>Proteção Florestal<br>IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso<br>Campus Cáceres<br>Caixa Postal 244<br>Avenida dos Ramires, s/n<br>Bairro: Distrito Industrial<br>Cáceres - MT CEP: 78.200-000<br>Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)<br>e-mails:<a ymailto="mailto:alexandresantosbr@yahoo.com.br" href="mailto:alexandresantosbr@yahoo.com.br">alexandresantosbr@yahoo.com.br</a><br> <a ymailto="mailto:alexandre.santos@cas.ifmt.edu.br" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a><br>Lattes: <a
href="http://lattes.cnpq.br/1360403201088680" target="_blank">http://lattes.cnpq.br/1360403201088680</a><br>======================================================================<br><br>_______________________________________________<br>R-br mailing list<br><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br><br></div> </div> </div> </div></body></html>