首先，利用核磁共振仪器得到3D TOF MRA序列的图像，利用不同切片的二维数据进行血管的三维重建对临床医生的诊断疾病具有重要重要意义。现在简要介绍一下MIP(maximum intensity projection)的原理。

如上图所示，136张血管的切片影像，假设以slice方向从后到前射出一条光线并投影到一个二维平面，光线所经过的像素值的最大值即为该二维平面上的图像的像素值。右图即为slice方向上的投影重建图像。可以想象，如果从各个方向分别对136张血管的切片影像进行投影重建，将得到各个方向上的血管造影图像。以下仅给出三个方向上的结果，

clear, close all

%% Import sequential images into a variable
% dirpath='006_Sag_3D_BRAVO_1.0iso';%读入数据
% dirpath='006_OAx3DTOFMRA3Slab';%读入数据
dirpath='.\data\MRI_TOF';
dirinfo=dir(dirpath);
dirinfo(1:2)=[];  % remove . and ..
img=[];
for i=1:length(dirinfo)
    img(:,:,i)=dicomread([dirpath filesep dirinfo(i).name]);    
end

%% Display images slice by slice using for-loop
figure,
for i=1:size(img,3)
    imshow(img(:,:,i),[],'border','tight')
    pause(0.01)
end

%%  Perform maximum intensity projection 
MIPimg=max(img,[],3); % projection along slice direction沿着某个方向取最大值；沿着slice方向
figure,subplot(2,2,1)
imshow(MIPimg,[],'border','tight')

因为代码实现比较简单，代码这里只给出slice方向。其他的给大家练手啦。TOF数据系某高校老师免费提供，不存在侵犯个人隐私。原始数据已经上传，大家可以下载学习。