MATLAB:如何在连续 block 中按列计算数字频率

Question

我在 Matlab 中有一个矩阵 a，如下所示:

a = zeros(10,3); 
a(3:6,1)=2; a(5:9,3)=1; a(5:7,2)=3; a(8:10,1)=2;

a =
 0     0     0
 0     0     0
 2     0     0
 2     0     0
 2     3     1
 2     3     1
 0     3     1
 2     0     1
 2     0     1
 2     0     0

我想获取一个元胞数组，其中包含每个数字在一列中出现的次数。此外，无论列号如何，都应根据元素值对其进行排序。在上面的示例中，我想获取单元格:

b = {[5],[4,3],[3]}

因为数字 1 在 5 次中出现一次，数字 2 在 4 和 3 的 block 中出现两次，数字 3 在 3 次中出现一次。如您所见，循环是根据元素值而不是元素出现的列数排序的。

Answer 1

由于您不关心列，您可以将所有列串成一个列向量，在两端填充零以防止列开始和结束处的跨度一起运行:

v = reshape(padarray(a, [1 0]), [], 1);
% Or if you don't have the Image Processing Toolbox function padarray...
v = reshape([zeros(1, size(a, 2)); a; zeros(1, size(a, 2))], [], 1);

现在，假设跨度总是由 1 个或多个零分隔，您可以按如下方式找到每个跨度的长度:

endPoints = find(diff(v) ~= 0);  % Find where transitions to or from 0 occur
spans = endPoints(2:2:end)-endPoints(1:2:end); % Index of transitions to 0 minus
                                               % index of transitions from 0

最后你可以 accumulate基于这些跨度中存在的值的跨度:

b = accumarray(v(endPoints(1:2:end)+1), spans, [], @(v) {v(:).'}).';

以你的例子为例:

b =

  1×3 cell array

    [5]    [1×2 double]    [3]

注意事项:

生成的元胞数组中值的顺序不能保证与 spans 中的顺序匹配(即上面的 b{2} 是 [3 4] 而不是 [4 3])。如果顺序很重要，您需要对下标进行排序 as per this section of the documentation .以下是如何更改 b 的计算:

[vals, index] = sort(v(endPoints(1:2:end)+1));
b = accumarray(vals, spans(index), [], @(v) {v(:).'}).';