如何解决如何旋转矢量/时间序列
例如我有这样一个向量 y
y <- c(1,7,4,8,6,1)
x <- c(1,2,3,5,7) # indexes
我想将向量旋转任意度数。
输入向量 --> 旋转函数 --> 旋转向量。
我找到了一个similar question,但我不擅长公式,有人可以举个例子吗
解决方法
您可以使用 sf
:
- 创建一条曲线为
linestring
- 创建一个rotation matrix
- 将曲线乘以旋转矩阵
library(sf)
x <- c(1,7,4,8,6,1)
curve <- st_linestring(cbind(1:length(x),x))
theta <- pi / 3
rotation <- matrix(c(cos(theta),sin(theta),-sin(theta),cos(theta)),ncol =2)
plot(curve)
plot(curve*rotation)
使用 {tidyverse} 和 {spdep}。
library(tidyverse)
library(spdep)
y <- c(1,1)
x <- c(1,2,3,5,7)
base <- data.frame(x,y)
rotated <-
data.frame(x = x,y = y) %>%
spdep::Rotation(angle = 0.5) %>%
data.frame() %>% rename(x = 1,y = 2)
bind_rows(base %>% mutate(id = "base"),rotated %>% mutate(id = "rotated")) %>%
ggplot(aes(x = x,y = y,colour = id)) +
geom_path()
由 reprex 包 (v2.0.0) 于 2021 年 4 月 23 日创建
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。