R语言树轮数据可视化二——以contour图展示

最新推荐文章于 2023-04-04 21:35:38 发布

niewzh

最新推荐文章于 2023-04-04 21:35:38 发布

阅读量2.2k

点赞数 1

分类专栏：统计学算法与技术 R语言

本文链接：https://blog.csdn.net/ScapeD/article/details/117215342

版权

R语言同时被 2 个专栏收录

19 篇文章 6 订阅

订阅专栏

统计学算法与技术

17 篇文章 5 订阅

订阅专栏

本文探讨了如何使用contour图展示复杂的树轮数据，通过实例展示了不同树芯宽度的时间序列变化，并强调了高低生长期的识别和序列敏感性的标准化处理。通过代码展示了如何利用filled.contour函数实现这一过程。

摘要由CSDN通过智能技术生成

一般而言，树轮数据多以折线图进行展示，清楚展示每条序列的变化，但序列一多通常难以比较。笔者尝试使用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)) },)