How to concatenate a set of rows created in a for loop?

Question

Shaun 2015-2-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/179397-how-to-concatenate-a-set-of-rows-created-in-a-for-loop

评论： Shaun 2015-2-25

I'm trying to create a 2D array by concatenating a set of 1D column arrays created in a for loop. The code I have so far.

     fileList = {};
  warning('off', 'all');
% Ask user for directory input and gather subfolders
[filePath] = uigetdir('','Select Directory of Folders to be Loaded');
fileList = dir(filePath);
fileList = fileList(arrayfun(@(x) x.name(1), fileList) ~= '.');
% Create array of files
dat = cell(size(fileList, 1), 6);
for j = 1:size(fileList, 1)
      % Open array of file and gather info
      fID = fileList(j).name;
      fID = strcat(filePath, '\',fID);
      info = dicominfo(fID);
      pID = info.PatientID;
      serNo = info.SeriesNumber;
      sliceLoc = info.SliceLocation;
      repet = info.RepetitionTime;
      echo = info.EchoTime;
      % Construct array of files with info
      dat(j,:) = {fID, pID, serNo, sliceLoc, repet, echo};
  end
  % Sort array by echo time
  dat = sortrows(dat, 6);
% Read files according to echo time
for j = 1:length(dat)
    [img] = dicomread(char(dat(j,1)));
    images(:,:,j) = img;
end
% plot T2 by matching each pixel, use T(1, 1, :) and iterate for all
% pixels.
vals = cell(size(images, 3), 2);
T2ind = cell(size(images, 1), 1);
T2 = cell(size(images, 1), size(images, 2));
% Loop through all rows
for y = 1:size(images,1)
    % Loop through all colunms
    for x = 1:size(images,2)
        % If the average value is greater than 100 for the given row and colunm number then loop through all 'pages'
        if ((sum(images(y,x,:))/size(images, 3)) > 100)
            for z = 1:size(images, 3)
                % Read echo time and pixel values and constuct array
                echo = dat(z, 6);
                imgVal = images(y,x,z);
                vals(z,:) = {echo, imgVal};
            end
            % Fit an exponential curve to pixel values and calculate T2
            f = fit(cellfun(@cell2mat, vals(:,1)), cell2mat(vals(:,2)), 'exp1');
            coeff = coeffvalues(f);
            a= coeff(1);
            b = coeff(2);
            S1 = f(10);
            S2 = S1/exp(1);
            T2ind{x} = ((log(S2/a))/b);
        end
    end
    T2(:,y) = (T2ind(:));
end
% T2(T2 < 0) = 0;
% T2(isnan(T2)) = 0;

The error I get is...

Subscripted assignment dimension mismatch.
Error in T2_mapping (line 61)
T2(:,y) = cat(2, T2ind{:});

Any help? Thanks!

6 个评论
显示 4更早的评论隐藏 4更早的评论

Shaun 2015-2-24

Any more help anyone?

Shaun 2015-2-24

在 MATLAB Online 中打开

So I now have a 256 by 256 array of cells! When I try to convert the cells I get a concatenation error....

cell2mat(T2)
Error using cat
Dimensions of matrices being concatenated are not consistent.
Error in cell2mat (line 86)
        m = cat(2,m{:});
 Any ideas on that?

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Guillaume 2015-2-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/179397-how-to-concatenate-a-set-of-rows-created-in-a-for-loop#answer_169084

编辑：Guillaume 2015-2-24

在 MATLAB Online 中打开

If I read your code right (comments needed!) T2Ind is a cell array of scalar or empties. Therefore horzcat(T2Ind{:}) is a vector of numbers. You're trying to assign that to a m*1 cell array (with m = size(images, 2)), so of course it's not going to work. I'm not sure what you're trying to do, so can't really suggest how to correct it.

Further to my comment above, it's possible you don't need all your loops. Again, I've not wrapped my head around the whole of your code, but

sum(images(y,x,:))/size(images, 3)

could be just calculated outside the loop as:

imagesmean = mean(images, 3);

Furthermore, what the heck is that supposed to do:

cell2mat(cellfun(@cell2mat, vals(:,1), 'UniformOutput', false))

5 个评论
显示 3更早的评论隐藏 3更早的评论

Shaun 2015-2-24

在 MATLAB Online 中打开

Hi,

Sorry for confusion, I've added some comments to the question. Hopefully they will clear a few things up. I'm a beginner at Matlab and so my code is messy...

To clarify, I am comparing pixel values across images. So I want to create a 256 by 256 array of numbers. The innermost loop compares the pixel values across the images for a given row and column. The outer two loops are to move along the rows and columns so that all pixels have been compared.

So the innermost loop gives me a single scalar, the second loop gives me a set of scalars (256 by 1) and I then want to outermost loop to give me an array of scalars (256 by 256)... I think....

Regarding the imagesmean comment, I want to assess whether the loop should be run on a pixel set by pixel set basis. Pixel set referring to a set of pixels with the same row and column throughout the images.

As far as I am aware....

cell2mat(cellfun(@cell2mat, vals(:,1), 'UniformOutput', false))

prepares the values for the fitting function to convert from cells..... I was given that code by someone and it seemed to work so I stuck with it. Very amateur, I know!!

Guillaume 2015-2-24

@Shaun,

Yes, I understood the purpose of your loop and I'm fairly certain they're not necessary or can be reduced to just one loop. The obscure bit is all this weird cell manipulation with echo and as mentioned what you're trying to do with your cell array assignment on the line with the error.

Shaun 2015-2-25

Thanks for your comments everyone, I think I have it working now. I re-wrote the code after actually reading up on the vectorization Stephen mentioned.... I can see why cell arrays aren't what I should have been using now! My code is much neater and nicer now and is giving me what I wanted. Thanks!

请先登录，再进行评论。