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

% （医学三维重建）MATLAB体绘制算法：多层面重建（MPR） % 演示代码 % % by HPC_ZY 20190716 clear; close all; clc load model % 请使用自己的体数据（三维矩阵） tic %% 模型参数 model = im2double(model); % 模型 [mRows,mCols,mDims] = size(model); % 模型尺寸 %% 设置采样网格参数 disp('设置网格参数') % 切面尺寸（自行设置） M = max([mRows,mCols,mDims]); N = max([mRows,mCols,mDims]); % 切面分辨率（自行设置） ms = 400; ns = 400; % 采样进度 precison = ([M,N]-1)./([ms,ns]-1); % 切面中心(默认值）（自行设置） mCenter = ([mRows,mCols,mDims]+1)/2; % 切面旋转角 alpha = 0; beta = 0; gamma = 0; %% 建立采样网格 disp('建立采样网格') % 转为弧度 alpha = alpha/180*pi; beta = beta/180*pi; gamma = gamma/180*pi; % 初始化网格 mesh.x = zeros(M,N); mesh.y = zeros(M,N); mesh.z = zeros(M,N); % 计算网格 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 = [x,y,z]; % 避免旋转中覆盖原值 y = tmp(2)*cos(alpha)-tmp(3)*sin(alpha); z = tmp(2)*sin(alpha)+tmp(3)*cos(alpha); % y轴逆时针旋转 tmp = [x,y,z]; % 避免旋转中覆盖原值 x = tmp(1)*cos(beta)+tmp(3)*sin(beta); z = -tmp(1)*sin(beta)+tmp(3)*cos(beta); % z轴逆时针旋转 tmp = [x,y,z]; % 避免旋转中覆盖原值 x = tmp(1)*cos(gamma)-tmp(2)*sin(gamma); y = tmp(1)*sin(gamma)+tmp(2)*cos(gamma); % 移回中心 mesh.x(r,c) = x+mCenter(1); mesh.y(r,c) = y+mCenter(2); mesh.z(r,c) = z+mCenter(3); end end %% 数据采样 disp('进行数据采样') im = zeros(ms,ns); % 注：为演示原理没有使用interp3()函数 for r = 1:ms for c = 1:ns % 此处使用最近邻插值 x = round(mesh.x(r,c)); y = round(mesh.y(r,c)); z = round(mesh.z(r,c)); if x>=1&&x<=mRows && y>=1&&y<=mCols && z>=1&&z<=mDims im(r,c) = model(x,y,z); end end end %% MPR体绘制 disp('MPR体绘制') toc imshow(im)