ggplot2를 사용하여지도에 경도와 위도를 추가하는 방법은 무엇입니까?

Question

나는 이제 ggplot2을 사용하여 캐나다지도를 계획하고 있습니다. 기본 투영 방법은 "aea"(Albers Equal Area)이므로 위도와 경도는지도의 직선입니다. 지도에서 "110W, 100W, 90W"및 "50N, 60N, 70N"의 형태로 경도와 위도를 어떻게 표시 할 수 있을지 궁금합니다. 곡선이어야합니다. 고마워.

의 ArcGIS shapfile이 같아야 https://www.arcgis.com/home/item.html?id=dcbcdf86939548af81efbd2d732336db

library(ggplot2) 
library(rgdal) 
countries<-readOGR("Canada.shp", layer="Canada") 
ggplot()+geom_polygon(data=countries,aes(x=long,y=lat,group=group),fill='white',color = "black") 

최종 결과로부터 다운로드된다.

Answer 1

당신은이 좌표 격자를 변경하는 예측을 사용 ggplot documented here

의 coord_map 인수와 함께이 작업을 수행 할 수 있습니다. 곡선 선은 등거리 투영을 포함하지만 허용되는 모든 투영 목록을 보려면 here으로 표시해야합니다. 어느 쪽을 선택 하는가는 선호하는 방식입니다.

(I이이 Azimuthal equidistant projection 생각) azequidistant를 사용하여 수동으로 레이블을 추가 :

axis_labels <- rbind(
        data.frame(long = rep(-140,5),lat = seq(40,80,10), labels = seq(40,80,10)), # x axis labels 
        data.frame(long = seq(-140,-60,40),lat = rep(85,3), labels = seq(140,60,-40)) # y axis labels 
) 

    ggplot() + 
    geom_polygon(data=countries,aes(x=long,y=lat,group=group),fill='white',color = "black") + 
    coord_map("azequidistant") + 
    scale_x_continuous(breaks = seq(-140,60, by = 20))+ 
    scale_y_continuous(breaks = seq(40,80, by = 10)) + 
    geom_text(data = axis_labels, aes(x = long, y = lat, label = labels)) + 
    theme_bw() + 
    theme(panel.grid.major = element_line(colour = "grey"), 
     panel.border = element_blank(), 
     axis.text = element_blank()) 

Answer 2

공간 데이터의 별도 계수 선 레이어를 사용할 수 있습니다.이 계수 데이터는 캐나다 레이어를 기반으로 투영됩니다.

무료 계수 선은 NaturalEarthData에서 다운로드 할 수 있습니다.

countries<-readOGR("Canada.shp", layer="Canada") 
grat <- readOGR("graticule.shp", layer="graticule") 

grat_prj <- spTransform(grat, CRS(countries)) 

ggplot() + 
geom_polygon(data=countries, aes(x=long,y=lat,group=group),fill='white',color = "black") + 
geom_path(data=grat_prj, aes(long, lat, group=group, fill=NULL), linetype="solid", color="grey50")