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

ASANTOS alexandresantosbr em yahoo.com.br
Sexta Julho 19 10:39:22 BRT 2013


Elias,

      Estou tentando selecionar uma pequena área de interesse na imagem 
para não demorar muito, mas agora quando tento o crop() tenho: Erro em 
function (classes, fdef, mtable)  :
   unable to find an inherited method for function "crop", for signature 
"SpatialGridDataFrame"

Na verdade criei um contorno SpatialPolygons e não SpatialGridDataFrame, 
mas a função não aceita segue CRM:

# Criar uma área de menor dimensão que a imagem inteira que abarque a 
região de interesse
xcc<-773759.1
ycc<-7841546
p.central<-cbind(xcc,ycc)

###Criando os vértices da área
coordV <- 
rbind(coordV,cbind(p.central[,1]+c(-1500,1500,1500,-1500,-1500),p.central[,2]+c(1500,1500,-1500,-1500,1500)))
coordV
coordV2=coordV[!is.na(coordV[,1]),]
coordV2
plot(coordV2[21:25,1],coordV2[21:25,2])
points(p.central[,1],p.central[,2], col="red")
#

# Cria um polígono com o contorno definido
bnds <- cbind(x=c(coordV2[21:25,1]), y=c(coordV2[21:25,2]))

# CRS UTM
SP <- SpatialPolygons(list(Polygons(list(Polygon(bnds)), "1")))
proj4string(SP) = CRS("+proj=utm +zone=23+south +datum=WGS84 +units=m 
+no_defs") ## Projeção

# CRS em lat long para o polígomno ficar na mesma projeção que o GeoTif
SPlatlong<- spTransform(SP, CRS("+proj=longlat +datum=WGS84"))
#

demlatlog<-readGDAL('19S435HN.tif')### Abre o slope TOPODATA - geotif do 
outro e-mail

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

#Cortar uma área menor de interesse
demlatlog.crop <- crop(demlatlog, extent(SPlatlong), snap='out')
#

Obrigado,

Alexandre

Em 18/07/2013 19:01, Elias Krainski escreveu:
> 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
>
>
> _______________________________________________
> 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
======================================================================

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130719/3fac013e/attachment.html>


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