医学三维重建MATLAB体绘制算法：多层面重建MPR

文章：（医学三维重建）MATLAB体绘制算法：多层面重建（MPR） 的测试代码

% （医学三维重建）MATLAB体绘制算法：多层面重建（MPR）
%  演示代码
%
% by HPC_ZY 20190716

clear; close all; clc

load model % 请使用自己的体数据（三维矩阵）

tic

%% 模型参数
model = im2double（model）; % 模型
[mRowsmColsmDims] = size（model）; % 模型尺寸


%% 设置采样网格参数
disp（‘设置网格参数‘）
% 切面尺寸（自行设置）
M = max（[mRowsmColsmDims]）;
N = max（[mRowsmColsmDims]）;
% 切面分辨率（自行设置）
ms = 400;
ns = 400;
% 采样进度
precison = （[MN]-1）./（[msns]-1）;
% 切面中心（默认值）（自行设置）
mCenter = （[mRowsmColsmDims]+1）/2;
% 切面旋转角
alpha = 0;
beta = 0;
gamma = 0;


%% 建立采样网格
disp（‘建立采样网格‘）
% 转为弧度
alpha = alpha/180*pi;
beta = beta/180*pi;
gamma = gamma/180*pi;
% 初始化网格
mesh.x = zeros（MN）;
mesh.y = zeros（MN）;
mesh.z = zeros（MN）;
% 计算网格
for r = 1:ms
    for c = 1:ns
        % 赋初值
        x = 1+（r-1）*precison（1）;
        y = 1+（c-1）*precison（2）;
        z = 0;
        % {中心旋转}
        % 移至原点
        x = x-mCenter（1）;
        y = y-mCenter（2）;
        % x轴逆时针旋转
        tmp = [xyz]; % 避免旋转中覆盖原值
        y = tmp（2）*cos（alpha）-tmp（3）*sin（alpha）;
        z = tmp（2）*sin（alpha）+tmp（3）*cos（alpha）;
        % y轴逆时针旋转
        tmp = [xyz]; % 避免旋转中覆盖原值
        x = tmp（1）*cos（beta）+tmp（3）*sin（beta）;
        z = -tmp（1）*sin（beta）+tmp（3）*cos（beta）;
        % z轴逆时针旋转
        tmp = [xyz]; % 避免旋转中覆盖原值
        x = tmp（1）*cos（gamma）-tmp（2）*sin（gamma）;
        y = tmp（1）*sin（gamma）+tmp（2）*cos（gamma）;
        % 移回中心
        mesh.x（rc） = x+mCenter（1）;
        mesh.y（rc） = y+mCenter（2）;
        mesh.z（rc） = z+mCenter（3）;
    end
end


%% 数据采样
disp（‘进行数据采样‘）
im = zeros（msns）;
% 注：为演示原理没有使用interp3（）函数
for r = 1:ms
    for c = 1:ns
        % 此处使用最近邻插值
        x = round（mesh.x（rc））;
        y = round（mesh.y（rc））;
        z = round（mesh.z（rc））;
        if x>=1&&x<=mRows && y>=1&&y<=mCols && z>=1&&z<=mDims
            im（rc） = model（xyz）;
        end
    end
end


%% MPR体绘制
disp（‘MPR体绘制‘）
toc

imshow（im）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2159  2019-07-16 16:50  demo1.m

     文件       1493  2019-07-17 09:17  demo2.m

----------- ---------  ---------- -----  ----

                 3652                    2