Extracting and summing data from a 3d matrix

Question

Emu 2022-9-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1801000-extracting-and-summing-data-from-a-3d-matrix

回答： Nipun 2023-9-27

Hi there,

I have some data in a 3d matrix 11*32*19. I want to a) identify all columns where there is a non-NaN and b) identify the blocks of non-NaNs (from 1 row long to 11 rowns long) and sum together the blocks in each of the 32 columns, giving a total value for each column against the column name for all the 19 'pages' combined. I have some code to do that for a 2d matrix (below), although i am just realising that it will miss the blocks that start at row 1 :/ so any suggestions on that front would be awesome too.

filename = '1001_1_2';
[locAll,typeAll]=find(isfinite(mVocs));  %extracting the locations of all the non NaNs (finite numbers)
indVocsAll=any(isfinite(mVocs));     
colindx = find(indVocsAll==1);
colres = zeros(length(colindx),1);
for i =1:length(colindx)
    tmpdat = mVocs(:,colindx(i));
    blocklength = zeros(length(tmpdat),1);
    for fi = 1:length(tmpdat)-1
        if isnan(tmpdat(fi)) && ~isnan(tmpdat(fi+1)) %finding where the column goes from NaN to 1
            count = fi+1;
            count2 = 1;
            while ~isnan(tmpdat(count)) %also count Non-Na
                count = count+1;
                count2 = count2+1;
            end
            %tmpdat2 = tmpdat(fi+1:fi+count2-1);
            blocklength(fi) = count2;
        end
    end
    blocklength = nonzeros(blocklength);
    savename = strcat(savedir, 'filename.mat');
    save(savename, 'blocklength');
    colres(i) = length(blocklength);
    colres(:,2) = colindx;
end

by counting the start of the number of blocks in each column - i am just realising that it will miss the blocks that start at row 1 :/ so any suggestions on that front would be awesome too.

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

Emu 2022-9-9

Hi there, this is the script for a 2d array. So column is the 2nd dimension (header?) (32). I want to convert this script to extract the same information across 3 dimensions (so different pages/planes yes). So identify all the values in each 2d column across all 3d pages. So if there are values in column 3 (header) in page/plane 1 and in column 3(header) in page/plane 7 - extrct that information and sum together. So using the code you gave me the other day I've extracted and just need to know how to sum together across multiple dimensions (I think).

mVocs is a 2d array with my original data (X*32). (time*variable). I have multiple of these. Within this I extracted some time blocks around events for multiple datasets - so i selected x number of rows and looked in these rows for data within the 32 header columns , and need to save this event related data seperately to compare between events in future which is why I saved it as 3d. I also need to look at the sum total for all events across each individual dataset.

So you mean I need to use something different to find - maybe ind2sub with find?

dpb 2022-9-9

If this is for a 3D array and can find what is needed on a 2D basis, then turn the 2D code into a function to which you can pass each plane of the 3D array in turn is one way.

I would still suggest that it isn't at all clear just what the end result is to be -- one can easily sum in MATLAB across any dimension by providing the optional 'DIM' parameter -- and use logical addressing across the array to select the elements needed/wanted. find itself is generally not needed for such operations; only if it is needed to return those locations themselves, not to address the elements.

But, one can generally arrange computations such that the results can be computed and then the locations of elements in the result obtained by their location rather than having to keep those first.

It would make things a whole lot easier to try to help here if you would provide a sample (small) dataset that fully illustrates the input conditions and then the desired output and describe the algorithm by which the output is obtained so that folks here can visualize the problem trying to be solved rather than having to try to figure that out from trying to debug non-working code.

Matt J 2022-9-12

编辑：Matt J 2022-9-12

Your demo.mat includes example input data, but it doesn't complete the example by showing the desired output for that data.

Emu 2022-9-13

编辑：Emu 2022-9-13

1001_1_2_preSumTypesOutput.mat

oh apologies - here is a desired output attached. So this is the summed number of events where there is at least one partial string of consecutive non-NaNs, what i want is the total number of all strings & partial strings too :) \sorry for being so unclear I am really new to this.

请先登录，再进行评论。

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

Answer 1

Nipun 2023-9-27

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1801000-extracting-and-summing-data-from-a-3d-matrix#answer_1319642

在 MATLAB Online 中打开

Hi Emu,

I understand that you have a three dimensional matrix with dimensions: 11,32,19 with few NaN values. And, you are trying to get a sum-block output against each column which has a non-NaN value. I assume that the output will be a 32 x 19 matrix with one summation for each column across all rows.

In the example code below, I have generated a random matrix with the required dimensions. Please refer to the documentation for sum function in MATLAB for any clarifications.

M = rand(11,32,19);
% Generate 100 random indices to make them nan
idx = [];
for i=1:12000
    x = randi(size(M,1),1);
    y = randi(size(M,2),1);
    z = randi(size(M,3),1);
    idx = [idx [x,y,z]];
    M(x,y,z) = NaN;
end
% sum over all rows for each column across all 'pages' (3rd dimension)
A = sum(M,1,"omitnan");
size(A)

Hope this helps.

Regards,

Nipun

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Extracting and summing data from a 3d matrix

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

回答（1 个）

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

Extracting and summing data from a 3d matrix

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

回答（1 个）

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

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

0 个评论
显示 -2更早的评论隐藏 -2更早的评论