在3D数组中按维度进行逻辑索引

Question

假设我有一个尺寸为NxMxR的数组stack。我希望在给定维度上使用逻辑索引来更新其值，而不使用临时变量。

使用临时变量的例子:

for k=1:R
    temp=stack(:,:,k);
    temp(temp<max(max(temp))/2)=NaN;
    cropped(:,:,k)=temp;
end

请注意，如果整个数组的条件相同，那么它将是微不足道的，但是它沿着R维度变化。如果我只是想在语法中减去阈值，那么bsxfun将是合适的:

cropped = bsxfun(@minus, stack, max(max(stack, [], 3))/2);

我发现临时变量的使用笨拙，无法找到为bsxfun定义函数的适当方法来做循环所做的事情。

有人有什么建议吗?

Answer 1

stack(bsxfun(@minus, stack, max(max(stack))/2) < 0) = NaN;

Answer 2

另一个bsxfun实现-

t1 = reshape(stack,size(stack,1)*size(stack,2),[])
t1(bsxfun(@lt,t1,max(t1)/2))=NaN
cropped = reshape(t1,size(stack))

t1 = max(reshape(stack,size(stack,1)*size(stack,2),[]))
stack(bsxfun(@lt,stack,permute(t1/2,[3 1 2])))=NaN

更接近路易斯的答案-

stack(bsxfun(@lt,stack,max(max(stack)/2)))=NaN