一般而言,树轮数据多以折线图进行展示,清楚展示每条序列的变化,但序列一多通常难以比较。笔者尝试使用contour图进行展示,可以清楚的展示高低生长期。
使用 filled.contour(x, y, z, levels, col) 函数对一个采样点的不同树芯宽度进行可视化,效果如下:
其中横轴表示时间(年),纵轴表示不同树轮序列的宽度。
可以看出,每条序列基本都存在先高生长后低生长的趋势。另外,不同序列的敏感性是不同的,需要标准化后进行进一步分析。
全部代码如下:
library(dplR)
data.tr <- read.rwl('D:/projects/ITRDBprocess/Cleaned datasets/europe/alb001.rwl')
data.dftr <- data.frame(data.tr)
rr <- as.matrix(data.dftr)
#plot(data.tr,plot.type="spag")
colors <- c('blue','red')
RB <- colorRampPalette(colors = colors)
col <- RB(nlev)
plot_col(nlev, col)
levels <- c(0,0.5,1,3)
filled.contour(x = 1:705, y = 1:54, z = rr,levels = levels,color.palette = RB,plot.axes = { axis(1, seq(1, 705, by = 10))
axis(2, seq(1, 54, by = 1)) },)