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

Jônatan jdtatsch em gmail.com
Quinta Julho 18 17:45:15 BRT 2013


 r <- raster(nc=1000, nr=1000)
 r[] <- 1:ncell(r)

## usando uma extent com dominio espacial em torno de -23, aparentemente
funciona
e <- extent(-28, -15, -50, -45)
 stackIm1.crop <- crop(r, e)
 plot(stackIm1.crop)
  stackIm1 <- raster(stackIm1.crop)
  stackIm1[] <- round(runif(ncell(stackIm1.crop))* 255)
 x11();plot(stackIm1)
 r.geo <- "+proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs"
 stackIm2 <- projectRaster(stackIm1, crs=r.geo,method="ngb")
 x11(); plot(stackIm2)



2013/7/18 ASANTOS <alexandresantosbr em yahoo.com.br>

> Não Elias,
>
>      Na realidade quero fazer uma correção topográfica em uma banda 4 do
> landsat em utm, para isso preciso converter uma imagem geoTiff de elevação
> do Topodata do Inpe que esta em lat long para utm e não estou conseguindo,
>
> Obrigado,
>
> Alexandre
>
>
> Em 18/07/2013 11:51, Elias Krainski escreveu:
>
>  Note que o sistema de projeção default já é o que você quer:
>> > r <- raster(nc=3, nr=5)
>> > r
>> class       : RasterLayer
>> dimensions  : 5, 3, 15  (nrow, ncol, ncell)
>> resolution  : 120, 36  (x, y)
>> extent      : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
>> coord. ref. : +proj=longlat +datum=WGS84
>>
>> Curiosidade: Você quer criar um raster para cobrir uma particular região
>> do globo?
>> ______________________________**_________________
>> R-br mailing list
>> R-br em listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-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<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<e-mails%3Aalexandresantosbr em yahoo.com.br>
>         alexandre.santos em cas.ifmt.edu.**br<alexandre.santos em cas.ifmt.edu.br>
> Lattes: http://lattes.cnpq.br/**1360403201088680<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<https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br>
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-**guia<http://www.leg.ufpr.br/r-br-guia>)
> e forneça código mínimo reproduzível.
>



-- 
###############################################################
##  Jônatan Dupont Tatsch
##  Professor do Departamento de Física
##  Centro de Ciências Exatas e Naturais (CCNE)
##  Universidade Federal de Santa Maria
##  Faixa de Camobi, Prédio 13 - Campus UFSM - Santa Maria, RS, Brasil -
97105-900
##  Telefone: +55(55)33012084
##  www.ufsm.br/meteorologia
###############################################################
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130718/89b9bbad/attachment.html>


Mais detalhes sobre a lista de discussão R-br