文章介绍了如何在Matlab中对3D矩阵进行最大值投影，以及如何不通过遍历方式来修改矩阵元素。通过`max()`函数实现了3D矩阵的最大值投影，展示了`dim=1,2,3`的不同效果。同时，提出了两种不使用遍历的矩阵元素替换方法，包括使用`find()`函数和直接矩阵操作，并给出了示例。
Matlab 3D MIP and change matrix values without Traversing arrays------Matlab 3维矩阵（3个维度的矩阵）的最大值投影，以及不用遍历数组实现替换满足条件的数组元素
1.MIP 3个维度的矩阵的最大值投影
我们假设有一副需要最大值投影的rawdata，400x400x192,我们希望他以最大值投影得到400x400大小的图像（以最大值投影）；
A: 400x400x192
data = A;
MIPimg = max(data,[],3);
这样就得到了沿着深度方向上投影的图像了，关于详细的max(A,[],dim)可以参照Matlab 的help document。
算了，为了方便大家理解，我这里直接解释一下dim=1,=2,=3的情况，如图：
图中的箭头代表着最大值投影方向，其中dim=1，2投影出来的矩阵维度都是400x192,只有dim=3，投影出来的为400x400.
2.在不使用遍历的情况下实现具体的矩阵元素的替换修改
在很对计算成像实验中我们涉及到的数据量是非常大的，用matlab对矩阵进行遍历实现元素筛查的话会使得运行时间非常长，增加运行成本，而且有时候会因为遍历计算计算量太大导致matlab内存不够用，从而无法正常运行程序。
因此采用矩阵运算的方法会使得计算变得非常简单 和快速
比如我有一个矩阵100x100的元素，我要让其中小于5的元素变为0，那么如果用for循环遍历是绝对可以实现的，但是你需要遍历10000次，这里不用遍历有两种方法可以在matlab中实现：
方法1：使用find函数
让矩阵A中小于5的值全部变为0，
A=[1 2 3; 4 5 6; 7 8 9];
[r,c] = find(A<5);
A(r,c)=0
结果： 这里还是有一点问题的，它应该只改变前4个元素，但是somehow改了6个，应该有点问题。
`A =
0 0 0
0 0 0
7 8 9`
方法2：直接矩阵改变
让矩阵A中小于5的值全部变为0，
A=[1 2 3; 4 5 6; 7 8 9];
A(A<5) = 0
结果：
A =
0 0 0
0 5 6
7 8 9
completely correct，而且更简单！
扫一扫
2819
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
