元胞数组中的操作

Question

下午好，

我有一个包含许多多波段图像的单元阵列（每个图像 6 个波段）。 我的多维数组称为“D”，有两列：第一列包含图像名称，第二列包含图像文件（包含行、带...）。首先，我想使用以下代码将其转换为双精度格式（现在是 uint16），但我收到以下错误“??? Error using ==> double Too many input arguments.”

D= double(C{:,1});

然后，我有这样的价值观：

a1=1957 
a2= 1900 
a3= 1554 
a4=1036 
a5=215 
a6=80 

我想要做的是创建一个循环（我正在简化问题），因为我需要对图像的每个波段进行相同的操作。代码如下所示：

for k = 1:band(D{:,1})
    band1= D(:,:,1);
    band2= D(:,:,2);
    ...
    band1*a1=output1
    band2*a2=output2
end 

但我真的不知道如何表达。真的提前谢谢了

---

在元胞阵列的每一行，都有一个多波段图像。这是我创建的方式：

fileNames = dir(fullfile(myPath, '*.tif'));
C = cell(length(fileNames), 2);
for k = 1:length(fileNames)
    filename = fileNames(k).name;
    C{k, 1} = filename;
    C{k, 2} = imread(filename);
end

Answer 1

我相信你可以像这样实现你所需要的：

% Put your a in an array:
a = [1957 1900 1554 1036 215 80];

fileNames = dir(fullfile(myPath, '*.tif'));
C = cell(length(fileNames), 2);
% Allocate output as well
output = cell(length(fileNames), 1);
for k = 1:length(fileNames)
    filename = fileNames(k).name;
    C{k, 1} = filename;
    imData = imread(filename);
    C{k, 2} = imData;

    % Loop over the bands and multiply
    curOutput = zeros(size(imData));
    for nBand = 1:size(imData, 3)
        curOutput(:, :, nBand) = a(nBand) * double(imData(:, :, nBand));
    end

    % Store output
    output{k} = curOutput;

    % Note: The above loop over the bands can be replaced by
    % output{k} = bsxfun(@times, double(imData), permute(a, [3 1 2]));
end

您也可以使用以下代码在一行中执行相同的操作

% After your loop...
Cdata = C(:, 2);

% OR:
%
% Generate some random data like yours (I think)
% Cdata = arrayfun(@(x){uint16((2 ^ 16 - 1) * rand(10, 10, 6))}, zeros(10, 1));

% Then do:
output = cellfun(@(x)bsxfun(@times, double(x), permute(a, [3 1 2])), Cdata, ...
    'UniformOutput', false)

这使用匿名函数一步完成您想要的转换和乘法。第三个参数周围的bsxfun 和permute 是沿元胞数组Cdata 的每个元素的第三维应用向量a 的乘法。结果是一个与Cdata 大小相同的元胞数组。