Apply a mask to 4D medical images in a vectorized way

2 次查看（过去 30 天）
AP
AP 2013-2-11
Dear All,
I have a 4D array V, with the size [144 144 12 25]. In fact, it is a 4D medical image with 12 stacked slices in 25 different time. I want to apply a mask, named M, to each image in the array in a vectorized way. M is a 144x144 array containing logical values. Could someone help me how to do it?
Thanks,
Ahmad
  1 个评论
José-Luis
José-Luis 2013-2-11
I notice that you don't seem accept many answers after asking questions. Some of the answers seemed valid to me. Please remember to accept answers if they help you, since it is the only "payment" the volunteers in this forum receive.

请先登录，再进行评论。

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

采纳的回答

José-Luis
José-Luis 2013-2-11
编辑：José-Luis 2013-2-12
data = rand(144,144,12,25);
your_mask = rand(144,144) > 0.5;
Assuming you want to multiply by the mask:
your_result = bsxfun(@times,data,mask);
Note that you can do operations other than multiplication. Please have a look at
doc bsxfun
  4 个评论
显示 2更早的评论
AP
AP 2013-2-11
I want the operation on each image. I thought bsxfun is only for each pixel and can't be applied to each image in a 4D array. Please correct me if I am wrong: so when we call bsxfun(@times,data,mask), bsxfun applies the mask to each 144x144. Is that right?
José-Luis
José-Luis 2013-2-12
Yes. That's what the "singleton expansion" in the documentation means.

请先登录，再进行评论。

更多回答（1 个）

Conrad
Conrad 2013-2-11
Something like this:
fV = V.*M(:,:,ones(1,size(V,3)),ones(1,size(V,4)));
Conrad
  2 个评论
AP
AP 2013-2-11
How can I convert M, with size [144 144] to M with [144 144 12 25]?
AP
AP 2013-2-11
编辑：AP 2013-2-11
I think, I have to use
repmat(M, [1 1 size(V,3) size(V,4)])
to duplicate M in other dimensions and make it [144 144 12 25].

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Matrix Indexing 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by