How to change map scale in ggplot2?

Asked

Viewed 423 times

3

I’m making a map with ggplot and ggmap. However, the axis scales are in decimal degrees, but I need them in degrees, minutes and seconds. Thank you for your cooperation!

1 answer

3


I saw your post in the English version and got some progress. There is the function dms package GEOmap which turns decimal into degrees. Follow my attempt code.

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

  • Thank you very much, Rafael! I would comment that I had not been able to add the values of the minutes... I’ll keep trying... anything up the code update. Thanks!

  • I managed to add the symbol ' minutes, but not 100% because it is a little far from the number. I used the function expression(~minute). The code and figure are updated

  • Excelente, Rafael!

Browser other questions tagged

You are not signed in. Login or sign up in order to post.