跟进:如何在R中制作日出图？

install.packages("maptools")  # 用于获取地理位置信息
install.packages("raster")    # 用于处理栅格数据
install.packages("ggplot2")   # 用于绘图
install.packages("ggthemes")  # 用于美化图表

library(maptools)
library(raster)
library(ggplot2)
library(ggthemes)

# 使用地名获取经纬度
location <- geocode("北京市")

# 或直接提供经纬度
location <- c(lon = 116.3974, lat = 39.9093)

sun <- suncalc(location, dates = Sys.Date(), keep = c("sunrise", "sunset"))
sunrise <- sun$sunrise

# 创建一个空白图表
plot <- ggplot() +
  theme_base() +
  theme(axis.text.x = element_blank(),
        axis.title.x = element_blank(),
        axis.ticks.x = element_blank())

# 添加日出时间线
plot <- plot +
  geom_vline(xintercept = sunrise, color = "red", linetype = "dashed")

# 显示图表
plot