我有一个数据集df:
library(gsheet)
df <- gsheet2tbl('docs.google.com/spreadsheets/d/1uSYy4gg4QTFoI81g1_WxI1YA0d6HfWIvv-RvlvMBLR8/edit?usp=sharing')
> head(df)
date elevation DO_mg_L
1 2016-06-07 10:00:00 268.8816 7.60
2 2016-06-07 10:00:00 272.1625 8.17
3 2016-06-07 10:00:00 275.3449 8.42
4 2016-06-07 10:00:00 265.3383 7.15
5 2016-06-07 10:00:00 262.2215 6.93
6 2016-06-07 10:00:00 258.9079 6.95
该数据包含一系列在水体中进行的深度剖面。因此,文件包含数据,高程和浓度。我想使用这些数据来创建时间序列等值线图,其显示y轴上的时间,x轴上的高度以及颜色代表氧浓度的等高线图。
然而,数据包含缺失数据的时段(我想留空),高程不一致,水位波动。因此,我在创建矩阵时遇到了麻烦。
理想情况下,我的情节看起来类似于this。
我试过了:
do.dep.tim <- with(df, interp(x = date, y = elevation, z = DO_mg_L))
filled.contour(x = df$x,
y = df$y,
z = df$z,
color.palette =
colorRampPalette(c("white", "blue")),
xlab = "Time",
ylab = "Elevation (Feet ASL)",
main = "Dissolved Oxygen",
key.title = title(main = "DO (mg L)", cex.main = 1))
我感谢您的帮助和建议。