3

estou fazendo um mapa com o ggplot e ggmap. Porém, as escalas dos eixos estão em graus decimais, mas eu preciso delas em graus, minutos e segundos. Agradeço pela colaboração!

1 Resposta 1

3

Eu vi o seu post na versão em inglês e consegui algum avanço. Existe a função dms do pacote GEOmap que transforma decimal em graus. Segue código da minha tentativa.

library(ggplot2)
library(ggmap)
library(GEOmap)

#get my map
city<- get_map(location = c(lon= -54.847, lat= -22.25),
               maptype = "satellite",zoom = 11,color="bw")

map<-ggmap(city,extent="normal")+
  xlab("Longitude")+ ylab("Latitude")

scale_x_longitude <- function(xmin=-180, xmax=180, step=1, ...) {
  xbreaks <- seq(xmin,xmax,step)
  xlabels <- unlist(
    lapply(xbreaks, function(x){
      ifelse(x < 0, parse(text=paste0(paste0(abs(dms(x)$d),"^{o}*"),
                               paste0(abs(dms(x)$m),expression(~minute)), "*W")), 
             ifelse(x > 0, parse(text=paste0(paste0(abs(dms(x)$d),"^{o}*"),
                                             paste0(abs(dms(x)$m),expression(~minute)),"*E")),
                    abs(dms(x))))}))
  return(scale_x_continuous("Longitude", breaks = xbreaks, labels = xlabels, expand = c(0, 0), ...))
}

scale_y_latitude <- function(ymin=-90, ymax=90, step=0.5, ...) {
  ybreaks <- seq(ymin,ymax,step)
  ylabels <- unlist(
    lapply(ybreaks, function(x){
      ifelse(x < 0, parse(text=paste0(paste0(abs(dms(x)$d),"^{o}*"),
                                      paste0(abs(dms(x)$m),expression(~minute)), "*S")), 
             ifelse(x > 0, parse(text=paste0(paste0(abs(dms(x)$d),"^{o}*"),
                                             paste0(abs(dms(x)$m),expression(~minute)),"*N")),
                    abs(dms(x))))}))
  return(scale_y_continuous("Latitude", breaks = ybreaks, labels = ylabels, expand = c(0, 0), ...))
} 

map +
  scale_x_longitude(-55.0,-54.7,.1) +
  scale_y_latitude(-22.4,-22.1,.1)

Imgur

3
  • Muito obrigado, Rafael! Iria comentar que não havia conseguido adicionar os valores dos minutos... Vou continuar tentando... qualquer coisa subo a atualização dos códigos. Valeu! Commented 16/08/2017 às 21:53
  • Consegui adicionar o símbolo ' de minutos, mas não ficou 100% porque ele tá um pouco distante do número. Utilizei a função expression(~minute). O código e a figura estão atualizadas Commented 17/08/2017 às 12:04
  • Excelente, Rafael! Commented 17/08/2017 às 17:32

Você deve fazer log-in para responder a esta pergunta.

Esta não é a resposta que você está procurando? Pesquise outras perguntas com a tag .