时间戳

首页 > 解决方案 > 从元胞数组中提取内容

matlab - 从元胞数组中提取内容

问题描述

我有一系列图像，存储在数组 A 中。所以 A 的每个条目都包含一个图像（矩阵）。所有矩阵的大小都相同。

现在我想提取特定位置（像素）的值，但我目前的方法似乎很慢，我认为可能有更好的方法来做到这一点。

% Create data that resembles my problem
N = 5
for i = 1:N
A{i} = rand(5,5);
end

% my current approach
I = size(A{1},1);
J = size(A{1},2);
val = zeros(N,1);
for i = 1:I
    for j = 1:J
        for k = 1:N
            B(k) = A{k}(i,j);
        end
        % do further operations on B for current i,j, don't save B
    end
end

我在想应该有一些方法，A{:}(i,j)或者vertcat(A{:}(i,j))两者都导致

??? Bad cell reference operation.

我正在使用 Matlab2008b。

有关更多信息，我之后在 B 上使用 fft。

以下是 Cris 的回答结果

|     Code     | # images | Extracting Values |   FFT    |  Overall  |
|--------------|----------|-------------------|----------|-----------|
| Original     | 16       | 12.809 s          | 19.728 s | 62.884 s  |
| Original     | 128      | 105.974 s         | 23.242 s | 177.280 s |
| ------------ | -------- | ----------------- | -------  | --------- |
| Answer       | 16       | 42.122 s          | 27.382 s | 104.565 s |
| Answer       | 128      | 36.807 s          | 26.623 s | 102.601 s |
| ------------ | -------- | ----------------- | -------  | --------- |
| Answer (mod) | 16       | 14.772 s          | 27.797 s | 77.784 s  |
| Answer (mod) | 128      | 13.637 s          | 28.095 s | 83.839 s  |

答案代码被修改为double(squeeze(A(i,j,:)));因为没有双倍的 FFT 需要更长的时间。

答案（mod）使用double(A(i,j,:));

因此，改进似乎确实适用于更大的图像集，但我目前计划每次运行处理约 500 张图像。

更新

用函数测量profile，使用/省略挤压的结果

|              Code              | # Calls |   Time   |
|--------------------------------|---------|----------|
| B = double(squeeze(A(i,j,:))); | 1431040 | 36.325 s |
| B= double(A(i,j,:));           | 1431040 | 14.289 s |

标签： matlabcell-array

解决方案

A{:}(i,j)不起作用，因为A{:}它是一个以逗号分隔的元素列表，相当于A{1},A{2},A{3},...A{end}. 索引到这样的数组是没有意义的。

为了加快您的操作，我建议您从数据中创建一个 3D 矩阵，如下所示：

A3 = cat(3,A{:});

当然，这仅在所有元素A具有相同大小的情况下才有效（如问题中最初指定的那样）。

现在您可以像这样快速访问数据：

for i = 1:I
    for j = 1:J
        B = squeeze(A3(i,j,:));
        % do further operations on B for current i,j, don't save B
    end
end

根据您应用于每个 的操作B，您也可以对这些操作进行矢量化。

编辑：由于您适用fft于每个 B，因此您也可以在不循环的情况下获得它：

B_fft = fft(A3,[],3); % 3 is the dimension along which to apply the FFT

推荐阅读