基于Matlab的MRI图像最大密度投影（MIP）

实现过程:

1）读取磁共振动态增强序列(DCE-MRI)一个序列中的多张slice
2）把多张slice的数据保存在一个三维数组中
3）返回Z轴（横轴位）方向像素最大值
4）显示MIP图像

代码

dirpath = 'H:\DCE+T2+ADC\00262266\MR\20130618\DCE00003';
dirinfo=dir(dirpath);
dirinfo(1:2)=[];
img=[];
%% 把每一张dicom文件读入img数组中
for i=1:length(dirinfo)
    img(:,:,i) = dicomread([dirpath filesep dirinfo(i).name]);
end
%% 展示每一张slice
figure,
for i =1:size(img,3)
    imshow(img(:,:,i),[],'border','tight');
    pause(0.02);
end
%% 最大密度投影
MIPimg = max(img,[],3);%沿着z轴方向（横轴位）取最大值；沿着slice方向
figure,subplot(1,1,1)
title('原图最大密度投影');
imshow(MIPimg,[],'border','tight')

实现效果：从图中可以看出病灶部分信号值明显高于其他地方

工具：Matlab R2018b

参考：
1.（医学三维重建）MATLAB体绘制算法：最大密度投影（MIP）
2.MPR、MIP、MinP定义
3.MATLAB和Python实现三维图像的最大强度投影

说明：本文为原创文章，转载或引用请注明网址和标题；有不正确的地方欢迎批评指正nick.yu.jd@qq.com