3-Dimensional Matrix and Standard Deviation

2021 5 23
1 个回答
回答已采纳
更新时间：2021 5 25
7 次查看（30 天）

0 个投票

Hi there,
I have a 3-Dimensional matrix (sizeX, sizeY, 700 Frames). I would like to group the "frames" in groups of 10, so that I can take the standard deviation of each of these groups of frames. E.g. -- Since I have 700 frames (z values), I would like to take the standard deviation of 1:10, then 11:20, then 21:30. For mean, I could just bin them in groups of ten (using for loop and mean function) and get my 70 values. But for standard deviation, I'm not sure how to do this.

2 个评论
显示 无 隐藏 无

Sami Case
Sami Case 2021-5-23
**Similar to the Grouped-Z Project function in ImageJ where you can do Standard Deviation of your data with a group size of 10.
Sami Case
Sami Case 2021-5-25
See DGM’s comment below for answer help.

请先登录，再进行评论。

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

 采纳的回答

0 个投票

Hi,
Here is a relatively simple solution for the standard deviation calculation:
D = DATA; % DATA of size: X-by-Y-Zs, where Zs = 1:700;
IND = 1:10:700;
for ii=2:numel(IND)
S_D(ii-1) = std2(D(:,:,IND(ii-1):IND(ii))); % Standard deviation
M_D(ii-1)=mean2(D(:,:, IND(ii-1):IND(ii))); % Mean values
end
Good luck.

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

Sami Case
Sami Case 2021-5-23
编辑：Sami Case 2021-5-23
This is returning a 1x70 vector, I was wanting to have a standard deviation for each pixel (160 x 160) based on 10 frame bins. So the dimensions of the output would be 160x160x70. Any help?
Sulaymon Eshkabilov
编辑：Sulaymon Eshkabilov 2021-5-23
It is quite simple:
...
IND = 0:10:700;
IND(1)=1;
for ii = ... % the rest is the same
...
Good luck.
Sami Case
Sami Case 2021-5-23
This is still showing a 1x70 vector. My goal it to have a standard deviation from 1-70 for each cell in the matrix. So my matrix is 160-by-160 and the z is 700. I would like to have a standard deviation of the cell in column 1, row 1 for every z. And for column 2, row 2 for every z. The output should be 160 by 160 by 70 matrix.
DGM
DGM 2021-5-24
编辑：DGM 2021-5-24
Try this:
D = rand(10,10,100); % random sample data
blocksize = 10; % how many pages to collapse?
npages = size(D,3);
IND = 1:blocksize:npages;
stdpict = zeros(size(D,1),size(D,2),numel(IND));
meanpict = zeros(size(D,1),size(D,2),numel(IND));
for ii=1:numel(IND)
idxrange = IND(ii):(IND(ii)+blocksize-1);
% Standard deviation of block along dim 3
stdpict(:,:,ii) = std(D(:,:,idxrange),0,3);
% Mean of block along dim 3
meanpict(:,:,ii) = sum(D(:,:,idxrange),3)/blocksize;
end
% the indexing will break if npages is not integer-divisible by blocksize
This does operations along dim3 of each block of pages. That sounds like what you're after.
Prior code was also indexing 1:11, 11:21, 21:31, etc, and dropping the last sample block.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Logical 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by