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
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)
-
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çãoexpression(~minute)
. O código e a figura estão atualizadas Commented 17/08/2017 às 12:04 -