我正在尝试在栅格地图上绘制纬度/经度网格。我用这样的shapefile完成了它:
ggplot(a, aes(long,lat, group=group,fill=group)) + #' a is the fortify object od my shp
geom_polygon(alpha=0.5,color='black',fill="lavenderblush3") +
coord_map("azequalarea") +
hrbrthemes::theme_ipsum_rc() +theme(legend.position="none")+ggsn::scalebar(a, dist = 100, st.size=3, height=0.01, dd2km = TRUE, model = 'WGS84')
但我不知道如何使用光栅文件执行此操作。有任何想法吗?
ggplot() + geom_raster(data = z, aes(x=x, y = y, fill=zabsm)) + #the z is the csv of my ratser
coord_fixed(ratio = 1) +
scale_fill_viridis(direction = -1) +
theme_bw()
我见过this example,但我希望将我的栅格地图放在albert / equal area投影中。
有什么想法吗?