[R-br] RgoogleMaps e rgdal

ASANTOS alexandresantosbr em yahoo.com.br
Terça Setembro 25 22:25:32 BRT 2012


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
>
>     ------------------------------------------------------------------------
>     *De:* Mauro Sznelwar <sznelwar em uol.com.br>
>     *Para:* Elias T. Krainski <eliaskrainski em yahoo.com.br>
>     *Enviadas:* Terça-feira, 4 de Setembro de 2012 23:03
>     *Assunto:* Re: [R-br] RgoogleMaps e rgdal
>
>     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)
>         Elias T. Krainski
>
>
>
>
>
> _______________________________________________
> 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
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 em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
======================================================================

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


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