Matlab中的矢量化范围检查

Question

在尝试将算法从 C# 移植到 Matlab 时，我发现 Matlab 在运行 for 循环时效率低下。因此，我想对算法进行矢量化。

我有以下输入：

低档：

[ 00 10 20 30 40 50 ... ]

远距离：

[ 10 20 30 40 50 60 ... ]

这些数组长度相等。

我现在有第三个数组Values（可以是任何长度），对于这个数组，我想计算Values 元素在lowerange(i) 和highrange(i) 之间的出现次数（你可以看到我来了来自 for 循环）。

输出应该是一个长度为低范围/高范围的数组。

所以用上面的数组输入LineData:

[ 1 2 3 4 6 11 12 16 31 34 45 ]

我希望得到：

[ 05 03 00 02 01 00 ...  ]

我尝试了（对我而言）显而易见的事情：

LineData(LineData < PixelEnd & LineData > PixelStart) 

但这不起作用，因为它只是逐个元素地检查元素上的 LineData。它不会尝试对LineData 中的所有值进行比较。

不幸的是，由于我还不习惯以 Matlab 的“向量”方式思考，所以我无法想出其他任何东西，更不用说从记忆中知道所有适用的指令了。

Answer 1

当您正在寻找具有给定边缘的基本直方图时，您可以使用 Matlabs 内置函数 histc：

values = [ 1 2 3 4 6 11 12 16 31 34 45 ];
edges = 0:10:60;
histc(values, edges)

ans =

 5     3     0     2     1     0     0

Answer 2

对于具有相同间隔并从 0 开始的范围，这是基于 bsxfun 的计数方法 -

LineData = [ 1 2 3 4 6 11 12 16 31 34 45 ] %// Input

interval = 10; %// interval width
num_itervals = 6; %// number of intervals

%// Get matches for each interval and sum them within each interval for the counts
out = sum(bsxfun(@eq,ceil(LineData(:)/interval),1:num_itervals)) 

输出 -

LineData =
     1     2     3     4     6    11    12    16    31    34    45
out =
     5     3     0     2     1     0

---

假设最后一个区间是输入数据最大值的区间，您也可以尝试基于 diff + indexing 的方法 -

LineData = [ 1 2 3 4 6 11 12 16 31 34 45 ] %// Input
interval = 10; %// interval width

labels = ceil(LineData(:)/interval); %// set labels to each input entry
df_labels = diff(labels)~=0; %// mark the change of labels
df_labels_pos = find([df_labels; 1]); %// get the positions of label change
intv_pos= labels([true;df_labels]);%// position of each interval with nonzero counts

%// get counts from interval between label position change and put at right places
out(intv_pos) =  [ df_labels_pos(1) ; diff(df_labels_pos)];