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=21...") 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