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

[Elias Krainski]

Caro Alexandre,

Vc esta usando aquele geotif do outro e-mail?

### obtem o arquivo
url=("http://www.webmapit.com.br/inpe/topodata/download/?ms=137398652719842&fbn=21S465&fsx=3") 

download.file(url, destfile = "21S465HN.zip")
### descompacta
system("unzip 21S465HN.zip")

### le o raster
require(rgdal)
r <- readGDAL("21S465SN.tif")

### atribui a projeção longlat  e datum WGS83 ao raster.
proj4string(r) <- CRS("+proj=longlat +datum=WGS84")

### transformar com é muito ineficiente
### r2 <- spTransform(r, CRS("+proj=utm +zone=23 +south"))

### com projectRaster(), se vc tivesse usando
### r <- raster("21S465SN.tiff")
### r2 <- projectRaster(r, crs=CRS("+proj=utm +zone=23 +south"))
### também é muito lento

### a lentidão é devida ao fato de que extrai-se as
### coordenadas de cada pixel e converte-se cada uma...

Att.
Elias
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130718/462c177e/attachment.html>