Question

我是这里的新手。

1，目的是绘制有关我的研究站点（名为RB1）的时间段（从2019-05到2019-10的8个日期）中的平均NDVI值的图表。并绘制垂直线以显示割草事件的日期。

2，现在我已经计算出这8个选定日期的NDVI值，并制作了一个CSV文件。（PS。“切割”是指研究地点的草地已被切割，因此应使用geom_vline将相应的日期显示为垂直线）

infor <- read_csv("plotting information.csv")
infor
# A tibble: 142 x 3
   date         NDVI cutting
   <date>      <dbl> <lgl>  
 1 2019-05-12 NA     NA     
 2 2019-05-13 NA     NA     
 3 2019-05-14 NA     NA     
 4 2019-05-15 NA     NA     
 5 2019-05-16 NA     NA     
 6 2019-05-17  0.787 TRUE      
# ... with 132 more rows

3，问题是，当我执行ggplot时，首先我想将x轴保留为整个时间段（2019-05至2019-10），但是当然不显示其间的所有日期，否则x轴上会显示太多日期）。因此，我执行scale_x_discrte(breaks=, labels=)以显示具有NDVI值的特定日期。

第二，我还想显示割草geom_vline的日期。

但是，看来scale_x_discrte的前提是factor我的日期，而geom_vline的前提是将日期保持为nummeric。这两个电话似乎是矛盾的。

y1 <- ggplot(infor, aes(factor(date), NDVI, group = 1)) +
  geom_point() +
  geom_line(data=infor[!is.na(infor$NDVI),]) + 
  scale_x_discrete(breaks = c("2019-05-17", "2019-06-18", "2019-06-26", "2019-06-28","2019-07-23","2019-07-28", "2019-08-27","2019-08-30", "2019-09-21"), 
                   labels = c("0517","0618","0626","0628","0723","0728", "0827","0830","0921"))) 


y2 <- ggplot(infor, aes(date, NDVI, group = 1)) +
  geom_point() +
  geom_line(data=infor[!is.na(infor$NDVI),]))

当我在y1中添加geom_vline时，我的绘图上不显示垂直线： y1 + geom_vline

当我将其添加到y2中时，显示了垂直线，但是日期（x轴）很奇怪（不显示为y1，因为我们不在此处运行scale_x_） y2 + geom_vline

   y1 + 
      geom_vline(data=filter(infor,cutting == "TRUE"), aes(xintercept = as.numeric(date)), color = "red", linetype ="dashed")

geom_vline和scale_x_discrete不兼容

0 个答案: