0

下午好，

我有一个包含很多多波段图像的单元阵列（每个图像有 6 个波段）。我的多维数组称为“D”，有两列：第一列包含图像名称，第二列包含图像文件（包含行、带...）。首先，我想使用以下代码将其转换为 double 格式（现在是 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
4

1 回答 1

3

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

% 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的乘法。结果是一个与 大小相同的元胞数组。permuteaCdataCdata

于 2013-04-28T17:30:05.987 回答