How can I remove axis line and add labels next to grids in a projected map plot created in ggplot2? I want to produce a map similar to this one:
Here is my reproducible example code:
adf <- structure(
list(long = c(249.999938964844, 152.500122070312, 152.500122070312,
249.999938964844, 249.999938964844),
lat = c(12.5000610351562, 12.5000610351562, 80, 80, 12.5000610351562),
order = 1:5,
hole = c(FALSE,FALSE, FALSE, FALSE, FALSE),
piece = structure(c(1L, 1L, 1L, 1L, 1L), .Label = "1", class = "factor"),
id = c("0", "0", "0","0", "0"),
group = structure(c(1L, 1L, 1L, 1L, 1L), .Label = "0.1", class = "factor"),
longg = c(251, 151, 151, 251, 251),
latt = c(11, 11, 81, 81, 11)),
.Names = c("long", "lat", "order", "hole", "piece", "id", "group", "longg", "latt"),
row.names = c(NA, -5L),
class = "data.frame")
adf
projplt <- ggplot() +
geom_path(data = adf,
aes(x = longg, y = latt, group = group),
colour = "Black", alpha = 1,size = 1.5) +
scale_x_continuous(limits = c(150.5, 252.0),
breaks = c(155, 170, 185, 200, 215, 230, 245, 260, 275),
labels = c("155°E", "170°E","175°W", "160°W", "145°W",
"130°W","115°W", "100°W", "85°W"),
expand = c(0,0))+
scale_y_continuous(limits = c(10.5, 81.5),
breaks = c(10, 25, 40, 55, 70),
labels = c("10°N", "25°N", "40°N", "55°N", "70°N"),
expand = c(0, 0))+
coord_map(projection = "albers", lat0 = 45, lat1 = 50)
projplt