Elias e Alexandre, fiz uma adaptação das duas rotinas e o resultado foi o mapa plotado em um arquivo: dados.pgn ### INICIO library(ggmap) library(RgoogleMaps) library(rgdal) ender <- rbind(c("AVENIDA BRIGADEIRO LUIS ANTONIO, 453", "01317000", "BELA VISTA", "SÃO PAULO", "SP"), c("RUA DR. COSTA JÚNIOR, 509", "05002000", "AGUA BRANCA", "SÃO PAULO", "SP"), c("AV. DUQUE DE CAXIAS, 193", "01214100", "SANTA IFIGENIA", "SÃO PAULO", "SP"), c("RUA MADRE DE DEUS, 427", "03119000", "MOOCA", "SÃO PAULO", "SP"), c("RUA CLODOMIRO AMAZONAS, 41", "04537010", "ITAIM BIBI", "SÃO PAULO", "SP"), c("RUA LOEFGREEN, 2007", "04040033", "VILA CLEMENTINO", "SÃO PAULO", "SP"), c("AV. CELIDONIO NETO, N. 165", "17120000", "CENTRO", "AGUDOS", "SP"), c("RUA CAPITÃO MIRANDA 57", "13900030", "CENTRO", "AMPARO", "SP")) colnames(ender) <- c("endereco", "cep", "bairro", "municipio", "UF") xy <- apply(ender, 1, function(x) geocode(paste0(x, collapse = ", "))) coord <- do.call("rbind", xy) ender2 <- cbind(ender, coord) ender2 options(OutDec=".") cad <- subset(ender2,select=c("lat","lon")) attach(cad,warn.conflicts=F) cad cad$size <- "tiny" cad$col <- "red" cad$char <- "" #get the bounding box: bb <- qbbox(lat = cad[,"lat"], lon = cad[,"lon"]) zoomx <- min(MaxZoom(latrange=bb$latR,lonrange=bb$lonR)); centro = c(mean(bb$latR), mean(bb$lonR)); #download the map: mapa <- GetMap.bbox(bb$lonR, bb$latR, destfile = "dados.png", GRAYSCALE =F, markers = cad, center=centro, maptype = "hybrid",zoom=zoomx); #determine the max zoom, so that all points fit on the plot plot.new() PlotOnStaticMap(mapa,lat = cad[,"lat"], lon = cad[,"lon"], cex=1.5,pch=20,col=c('blue', 'green', 'red'), add=F); ####### FIM Mauricio UFBA Em 26 de setembro de 2012 18:34, Alexandre Santos < alexandresantosbr@yahoo.com.br> escreveu:
Obrigado Elias,
É exatamente este o local, muito didático este recurso.
Redobrados agradecimentos,
Alexandre
------------------------------ *De:* Elias T. Krainski <eliaskrainski@yahoo.com.br> *Para:* "r-br@listas.c3sl.ufpr.br" <r-br@listas.c3sl.ufpr.br> *Enviadas:* Quarta-feira, 26 de Setembro de 2012 13:15
*Assunto:* Re: [R-br] RgoogleMaps e rgdal
Vc não precisa transformar para longlat. Basta informar corretamente qual a projeção das suas coordenadas e a função plotGoogleMaps() trata a projeção internamente para colocar na mesma do googleMaps. No seu script, faltou " " (espaço) na string CRS (antes de cada "+"). Para mim não ocorreu erro e visualizei um talhão ao sul do município de Cáceres - MT, a oeste da BR-070 (confere?). Veja o script:
x <- c(426588,426593,426825,426829,426588) y <- c(8216643,8216681,8216691,8216646,8216643)
### 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=utm +zone=21 +south +datum=WGS84 +units=m +no_defs")
### sobrepoe ao google Maps e visualiza require(plotGoogleMaps) plotGoogleMaps(Pontos)
Att.
Elias T. Krainski
------------------------------ *De:* Alexandre Santos <alexandresantosbr@yahoo.com.br> *Para:* "r-br@listas.c3sl.ufpr.br" <r-br@listas.c3sl.ufpr.br>; Elias T. Krainski <eliaskrainski@yahoo.com.br> *Enviadas:* Quarta-feira, 26 de Setembro de 2012 12:57
*Assunto:* Re: [R-br] RgoogleMaps e rgdal
Obrigado pelo help Elias,
Aquelas coordenadas na Bahia não eram minhas não, apenas peguei o post original e estava tentando adaptar para minhas coordenadas que são outras, atentei para a projeção, medida (m) e zona (21K) e fiz uma transformação para latlong e mesmo assim não deu certo, segue CRM abaixo:
### coordenadas dos pontos x <- c(426588,426593,426825,426829,426588) y <- c(8216643,8216681,8216691,8216646,8216643) # # ### 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)))
proj4string(Pontos) <- CRS('+proj=utm +zone=21+datum=WGS84+units=m+south') #### Pontos originalmente em UTM
t_latlong<- spTransform(Pontos, CRS("+proj=longlat +datum=WGS84")) ## Transformando para longlat
require(plotGoogleMaps) plotGoogleMaps(t_latlong) #
------------------------------ *De:* Elias T. Krainski <eliaskrainski@yahoo.com.br> *Para:* "r-br@listas.c3sl.ufpr.br" <r-br@listas.c3sl.ufpr.br> *Enviadas:* Quarta-feira, 26 de Setembro de 2012 10:43
*Assunto:* Re: [R-br] RgoogleMaps e rgdal
Isso porque suas coordenadas estão fora do range válido. Para projeção UTM, longitude varia entre 174W a 180E e latitude entre 80S a 84N. Seu primeiro ponto, por exemplo, (524227, 7978963) sugere isso e sugere também que suas coordenadas estão em alguma projeção onde a unidade de distancia é em metros, não em graus. Acho melhor você verificar qual é a projeção dos seus pontos.
No primeiro post seu sobre este assunto você estava com coordenadas em graus e, usando "+proj=longlat +datum=WGS84" não deu esse problema, pois visualizei os cinco pontos numa região noroeste da Bahia, a uns 30km da divisa com Piaui. Veja o script novamente:
### 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(Pontos)
### melhores informacoes ver vignette("plotGoogleMaps-intro")
Elias T. Krainski
------------------------------ *De:* Alexandre Santos <alexandresantosbr@yahoo.com.br> *Para:* "r-br@listas.c3sl.ufpr.br" <r-br@listas.c3sl.ufpr.br>; Alexandre Santos <alexandresantosbr@yahoo.com.br> *Enviadas:* Quarta-feira, 26 de Setembro de 2012 11:14
*Assunto:* Re: [R-br] RgoogleMaps e rgdal
Elias,
O único problema em trabalhar com UTM é que a função plotGoogleMaps() não me permite selecionar a zona/setor, então qualquer coordenada que coloco o navegador retorna: Sentimos muito, mas não temos imagens daqui. Gostaria de encontrar um meio de não ter que transformar as coordenadas para latlong,
Obrigado,
Alexandre
------------------------------ *De:* Alexandre Santos <alexandresantosbr@yahoo.com.br> *Para:* "r-br@listas.c3sl.ufpr.br" <r-br@listas.c3sl.ufpr.br>; Elias T. Krainski <eliaskrainski@yahoo.com.br> *Enviadas:* Quarta-feira, 26 de Setembro de 2012 9:54 *Assunto:* Re: [R-br] RgoogleMaps e rgdal [RESOLVIDO]
Obrigado Elias,
Passou despercebido, corrigido: ### 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)
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=utm +datum=WGS84') require(plotGoogleMaps) plotGoogleMaps(Pontos)
------------------------------ *De:* Elias T. Krainski <eliaskrainski@yahoo.com.br> *Para:* "r-br@listas.c3sl.ufpr.br" <r-br@listas.c3sl.ufpr.br> *Enviadas:* Quarta-feira, 26 de Setembro de 2012 8:26 *Assunto:* Re: [R-br] RgoogleMaps e rgdal
Essas coordenadas não estão em latlong. Coordenadas em latlong são em graus e vão de -180 a 180 (longitude) e de -90 a 90 (latitude).
Elias T. Krainski
------------------------------ *De:* ASANTOS <alexandresantosbr@yahoo.com.br> *Para:* r-br@listas.c3sl.ufpr.br *Enviadas:* Terça-feira, 25 de Setembro de 2012 22:25 *Assunto:* Re: [R-br] RgoogleMaps e rgdal
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@uol.com.br> <sznelwar@uol.com.br> *Para:* Elias T. Krainski <eliaskrainski@yahoo.com.br><eliaskrainski@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 listR-br@listas.c3sl.ufpr.brhttps://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 ======================================================================
_______________________________________________ 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.
_______________________________________________ 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.
_______________________________________________ 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.
_______________________________________________ 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.
_______________________________________________ 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.
_______________________________________________ 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.
_______________________________________________ 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.
-- Mauricio Cardeal UFBA
