Boa noite pessoal,
Refazendo o exemplo do Elias em coordenadas UTM, não tive
sucesso (Erro: Geographical CRS given to non-conformant data). Olhei
o datum e estava OK e também no manual e não me parece que há nada
de errado, inclusive o exemplo clássico do meuse do pacote gstat
utilizado no exemplo do pacote plotGoogleMaps esta em UTM, alguém
poderia dar um help, segue CRM abaixo
#
### coordenadas dos pontos
x <- c(524227,524268,524268,524227,524227)
y <- c(7978963,7978963,7978846,7978846,7978963)
#
#
### monta SpatialPointsDataFrame, como nao ha
### dados, criado um data.frame com a coluna ID
require(sp)
Carregando pacotes exigidos: sp
Pontos <- SpatialPointsDataFrame(SpatialPoints(cbind(x,y)),
+ data=data.frame(ID=1:length(x)))
### assumindo uma particular projecao, e'
### necessario verificar se e' mesmo esta
proj4string(Pontos) <- CRS('+proj=longlat +datum=WGS84')
Erro em `proj4string<-`(`*tmp*`, value = <S4 object of class
"CRS">) :
Geographical CRS given to non-conformant data: 524268 7978963
require(plotGoogleMaps)
plotGoogleMaps(Pontos)
Erro em spTransform(xSP, CRSobj, ...) :
No transformation possible from NA reference system
Em 05/09/2012 07:44, Elias T. Krainski
escreveu:
Lamento postar codigo sem teste... Este agora foi
testado e funciona.
###
coordenadas dos pontos
x <- c(-43.565,-43.564,-43.563,-43.565,-43.564)
y <- c(-10.518,-10.466,-10.412,-10.516,-10.474)
### monta SpatialPointsDataFrame, como nao ha
### dados, criado um data.frame com a coluna ID
require(sp)
Pontos <-
SpatialPointsDataFrame(SpatialPoints(cbind(x,y)),
data=data.frame(ID=1:length(x)))
### assumindo uma particular projecao, e'
### necessario verificar se e' mesmo esta
proj4string(Pontos) <- CRS('+proj=longlat
+datum=WGS84')
### sobrepoe ao google Maps e visualiza
require(plotGoogleMaps)
plotGoogleMaps(xyd)
### melhores informacoes ver
vignette("plotGoogleMaps-intro")
Elias T. Krainski
Continua não dando certo:
Creito que tem um erro no xy <-
SpatialPoints(cbind(xy)) que deveria ser
xy
<- SpatialPoints(cbind(x,y))
E
mesmo corrigindo veja o acentece:
>
plotGoogleMaps(x,y)
Erro em function (classes, fdef,
mtable) :
unable to find an inherited method
for function "spTransform", for
signature "numeric", "CRS"
>
Troquei
os pacotes. É necessário o pacote
plotGoogleMaps. Corrigindo:
y <-
c(-43.565,-43.564,-43.563,-43.565,-43.564)
x
<- c(-10.518,-10.466,-10.412,-10.516,-10.474)
require(sp)
xy
<- SpatialPoints(cbind(xy))
require(plotGoogleMaps)
plotGoogleMaps(xy)
_______________________________________________
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
Avenida dos Ramires, s/n
Caixa Postal 744
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
======================================================================