小编num*_*ers的帖子

我有20套以下数据集，并希望将两个图与两个y轴合为一个，以正确反映变量值。

hour=c(0:23)
conc=c(20.7, 19.4, 15.6, 11.7, 10.3, 9.1, 9.7, 10.6, 12.7, 12.6, 11, 9.8,
       9.3, 9.3, 10.3, 12.9, 17.1, 22.5, 22, 22.3, 23, 24.8, 24, 26.8)
ws=c(0.75, 0.7, 0.68, 0.64, 0.61, 0.64, 0.57, 0.58, .62, .65, .85, 1.12, 
     1.21, 1.29, 1.2, 1.18, 1.0, 0.84, .87, .75, .69, .69, .77, .74)
wd=c(295, 299, 288, 290, 292, 296, 306, 303, 300, 293, 259, 231, 94, 119, 
     95, 47, 75, 306, 302, 304, 319, 309, 290, 298)
mydf=data.frame(hour, conc, ws, wd)

我可以使用以下代码将它们分别绘制为1-by-2图：

conchrly …

我为特定一天的每个小时创建了风动画（单击播放动画视频）。我不想显示 19 个点，而是想创建一个在整个区域每小时使用这 19 个点进行插值/外推的等高线图 - 就像使用 ArcGIS 和样条插值工具生成的等高线图一样。

以下代码显示了我用来创建每小时风动画的 ggplot 和 gganimate。我只创建了一个小数据框作为完整 24 小时数据的子样本，因为我不熟悉将 csv 附加到此论坛中。有没有办法包含覆盖该区域的轮廓而不是geom_point？

library(ggplot2)
library(ggmap)
library(gganimate)

site <- c(1:18, 1:18)    
date <- data.frame(date=c(rep(as.POSIXct("2018-06-07 00:00:00"),18),rep(as.POSIXct("2018-06-07 01:00:00"),18)))    
long <- c(171.2496,171.1985, 171.2076, 171.2236,171.2165,171.2473,171.2448,171.2416,171.2243,171.2282,171.2344,171.2153,171.2532,171.2444,171.2443,171.2330,171.2356,171.2243)   
lati <- c(-44.40450,-44.38520,-44.38530,-44.38750,-44.39195,-44.41436,-44.38798,-44.38934,-44.37958,-44.37836,-44.37336,-44.37909,-44.40801, -44.40472,-44.39558,-44.40971,-44.39577,-44.39780)    
PM <- c(57,33,25,48,34,31,52,48,31,51,44,21,61,53,49,34,60,18,41,26,28,26,26,18,32,28,27,29,22,16,34,42,37,28,33,9)    
ws <- c(0.8, 0.1, 0.4, 0.4, 0.2, 0.1, 0.4, 0.2, 0.3, 0.3, 0.2, 0.7, NaN, 0.4, 0.3, 0.4, 0.3, 0.3, 0.8, 0.2, 0.4, 0.4, 0.1, 0.5, 0.5, 0.2, 0.3, 0.3, 0.3, 0.4, NaN, 0.5, 0.5, …