%?（醫學三維重建）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