Elias,
Na verdade o problema esta na função spTransform() que não esta
OK, porque o meu bbox em UTM (Zona 23K) é:
> bbox(SP)
min max
x 772259.1 775259.1
y 7840046.3 7843046.3
e após spTransform():
> bbox(SPlatlong)
min max
x -37.67016 -37.58042
y 70.51439 70.54432
O que esta incorreto pois 772259.1 em UTM (Zona 23K) no bbox(SP) deveria
ser -42.4058070334 em graus e não -37.67016 no bbox(SPlatlong) e para os
outros valores a mesma coisa, o que pode estar acontecendo?
Segue CRM:
require(sp)
require(rgdal)
# 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]),]
#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 +south +datum=WGS84"))
#
Obrigado,
Alexandre
Em 19/07/2013 17:33, Elias Krainski escreveu:
isso significa que o quadrado que você quer 'recortar'
bbox(SPlatlong)
não tem intersecção com o do raster
bbox(demlatlong)
_______________________________________________
R-br mailing list
[email protected]
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:[email protected]
[email protected]
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================
_______________________________________________
R-br mailing list
[email protected]
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.
