%?（醫(yī)學(xué)三維重建）MATLAB體繪制算法：多層面重建（MPR）
%??演示代碼
%
%?by?HPC_ZY?20190716

clear;?close?all;?clc

load?model?%?請使用自己的體數(shù)據(jù)（三維矩陣）

tic

%%?模型參數(shù)
model?=?im2double（model）;?%?模型
[mRowsmColsmDims]?=?size（model）;?%?模型尺寸


%%?設(shè)置采樣網(wǎng)格參數(shù)
disp（‘設(shè)置網(wǎng)格參數(shù)‘）
%?切面尺寸（自行設(shè)置）
M?=?max（[mRowsmColsmDims]）;
N?=?max（[mRowsmColsmDims]）;
%?切面分辨率（自行設(shè)置）
ms?=?400;
ns?=?400;
%?采樣進(jìn)度
precison?=?（[MN]-1）./（[msns]-1）;
%?切面中心（默認(rèn)值）（自行設(shè)置）
mCenter?=?（[mRowsmColsmDims]+1）/2;
%?切面旋轉(zhuǎn)角
alpha?=?0;
beta?=?0;
gamma?=?0;


%%?建立采樣網(wǎng)格
disp（‘建立采樣網(wǎng)格‘）
%?轉(zhuǎn)為弧度
alpha?=?alpha/180*pi;
beta?=?beta/180*pi;
gamma?=?gamma/180*pi;
%?初始化網(wǎng)格
mesh.x?=?zeros（MN）;
mesh.y?=?zeros（MN）;
mesh.z?=?zeros（MN）;
%?計(jì)算網(wǎng)格
for?r?=?1:ms
????for?c?=?1:ns
????????%?賦初值
????????x?=?1+（r-1）*precison（1）;
????????y?=?1+（c-1）*precison（2）;
????????z?=?0;
????????%?{中心旋轉(zhuǎn)}
????????%?移至原點(diǎn)
????????x?=?x-mCenter（1）;
????????y?=?y-mCenter（2）;
????????%?x軸逆時(shí)針旋轉(zhuǎn)
????????tmp?=?[xyz];?%?避免旋轉(zhuǎn)中覆蓋原值
????????y?=?tmp（2）*cos（alpha）-tmp（3）*sin（alpha）;
????????z?=?tmp（2）*sin（alpha）+tmp（3）*cos（alpha）;
????????%?y軸逆時(shí)針旋轉(zhuǎn)
????????tmp?=?[xyz];?%?避免旋轉(zhuǎn)中覆蓋原值
????????x?=?tmp（1）*cos（beta）+tmp（3）*sin（beta）;
????????z?=?-tmp（1）*sin（beta）+tmp（3）*cos（beta）;
????????%?z軸逆時(shí)針旋轉(zhuǎn)
????????tmp?=?[xyz];?%?避免旋轉(zhuǎn)中覆蓋原值
????????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


%%?數(shù)據(jù)采樣
disp（‘進(jìn)行數(shù)據(jù)采樣‘）
im?=?zeros（msns）;
%?注：為演示原理沒有使用interp3（）函數(shù)
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）

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

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

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

-----------?---------??----------?-----??----

?????????????????3652????????????????????2