如果有人向我展示了一个简单的方法,我将不胜感激。假设我在MATLAB中有一个向量,如
d = [3 2 4 2 2 2 3 5 1 1 2 1 2 2 2 2 2 9 2]
我想找到一系列连续数字“twos”以及这些系列的长度。
x=find(d==2)
可以轻松找到数字二。但我想要的是得到一个包含所有连续数字序列的长度的向量,这意味着我在这种情况下的结果将是这样的向量:
[1 3 1 5 1].
任何可以帮助我的人?
答案 0 :(得分:9)
这似乎有效:
q = diff([0 d 0] == 2);
v = find(q == -1) - find(q == 1);
给出
v =
1 3 1 5 1
对我来说
答案 1 :(得分:6)
这称为行程长度编码。在http://www.mathworks.com/matlabcentral/fileexchange/4955-rle-deencoding有一个很好的m文件可用。此方法通常比先前发布的diff / find方式更快。</ p>
tic
d_rle = rle(d==2);
d_rle{2}(d_rle{1}==1);
toc
经过的时间是0.002632秒。
tic
q = [0 diff([0 d 0] == 2)];
find(q == -1) - find(q == 1);
toc
经过的时间是0.003061秒。
答案 2 :(得分:-1)
如果我们想要连续值所在的原始矩阵的索引怎么办?此外,如果我们想要一个与原始矩阵大小相同的矩阵,其中连续值的数量存储在连续值的索引中,该怎么办?例如:
original_matrix = [1 1 1;2 2 3; 1 2 3];
output_matrix = [3 3 3;2 2 0;0 0 0];
此问题与气象数据质量控制有关。例如,如果我有一个来自多个传感器的温度数据矩阵,并且我想知道哪些天具有恒定的连续值,以及多少天是恒定的,那么我可以将数据标记为可能有故障。
温度矩阵是天数x站数,我想要一个输出矩阵,它也是天数x站数,其中连续值如上所述被标记。