连续数字系列(不同长度)

时间:2011-10-21 13:49:23

标签: matlab vector run-length-encoding

如果有人向我展示了一个简单的方法,我将不胜感激。假设我在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].

任何可以帮助我的人?

3 个答案:

答案 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站数,其中连续值如上所述被标记。

相关问题