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=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
_______________________________________________ R-br mailing list R-br@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@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================