平滑R中的方向(角度)数据

2019年8月5日 281次阅读

我试图在导出数据后处理一些运动分析软件跟踪错误.对于某些帧,方向从“真实”方向旋转180度.

我想平滑数据集,以便当方向在单个帧中改变~180时,它被转换以反映实际角度.

是否有人知道使用R语言中的任何循环统计软件包(如CircStats)解决此问题的方法？或者,我可以设想一个脚本来检查帧到帧的变化是否接近180度,如果这是真的则减去180,然后移动到下一帧.这听起来像是一种合理的方法,是否可以在R中轻松实现？

我担心我没有代表上传描述问题的图形(很容易看到),但here是一个示例数据集.

谢谢您的帮助.我一直是堆栈溢出的长期用户,但在之前需要询问之前从未找到我的答案.

大卫

编辑 – 附图

最佳答案 这是一个有趣的问题需要解决！它需要迭代,因为每当一个值改变时,它可以解决一个问题,但创建另一个…让我知道它是否成功.

threshold <- 90
correction <- 180

dat <- read.table("angle_data.txt", header=TRUE)
dat <- ts(dat)

repeat {
    diffs <- dat - lag(dat, k = 1)
    probl <- which(abs(diffs[,2]) > threshold)

    if(length(probl)==0)
        break

    obs.1 <- dat[probl[1], 2]
    obs.2 <- dat[probl[1] + 1, 2]

    dat[probl[1] + 1, 2] <- obs.2 + sign(obs.1 - obs.2) * 180
}
点赞