1.3个维度的矩阵投影，2.不用遍历而改变满足条件的矩阵元素

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，而且更简单！

完结撒花！！