Matlab中关于phantom生成头部仿真模体函数的使用

phantom函数使用方式

Matlab使用phantom函数会生成一个头部模体图像，具体使用方式如下：

P = phantom(def, n)
P = phantom(E,n)
[P,E] = phantom(E,n)

生成Matlab自带模体图像

P=phantom(def，n)生成头部体模的图像，该图像可用于测试radon和iradon或其他二维重建算法的数值精度。P是一个灰度强度图像，由一个大椭圆（代表大脑）组成，其中包含几个小椭圆（代表脑中的特征）。
def是一个字符串，用于指定要生成的头部幻影的类型。有效值为：
（1）“Shepp Logan”-研究人员在断层扫描中广泛使用的测试图像；
（2）“Modified Shepp-Logan”（默认）-Shepp-Logan模体的变体，其中对比度得到改善，以获得更好的视觉感知.
n是一个标量，用于指定P中的行数和列数。如果省略参数，n默认为256。
下面给出两种方式产生模体的对比，

I=phantom('Modified Shepp-Logan',256);  %生成改进的头部模型图
K = phantom('Shepp-Logan',256);         %生成头部模型图
%显示图像，选择imtool工具可有效调节窗宽窗位
imtool(I); 
imtool(K); 

在显示范围为0-1时，使用Modified Shepp-Logan参数生成的头部图像显示为：

此时使用Shepp-Logan参数生成的头部图像显示为：

修改Shepp-Logan参数生成的头部图像显示的范围为0-0.1，此时图像显示为：

可以知道使用Modified Shepp-Logan参数生成的头部图像在内部的细节中有更强的对比度，更有易于做后续的处理分析，推荐使用Modified Shepp-Logan参数。

生成自定义模体

P = phantom(E,n)
[P,E] = phantom(E,n)

这两种使用方式可以根据自己需求生成自定义模体图像，其中E为模体参数矩阵，E矩阵中每行均为六个参数，分别表示：

其中，x轴和y轴的定义域为[-1,1]，必须在此范围内指定第2列到第5列的参数值。[P,E] = phantom(E,n)还会输出参数矩阵E的值。下面简单举例如何使用：

A = [0 0.11 0.21 0.31 0.11 0.21 0.5 0 0.11 0.21 0.31 0.11 0.21]';
a = [0 0.1 0.15 0.1 0.15 0.2 1 0 0.11 0.21 0.31 0.11 0.21]';
b = [0 0.15 0.2 0.2 0.1 0.15 1 0 0.15 0.2 0.2 0.1 0.15]';
x0= [0 0.3 0.2 0.1 0.2 0.3 0 0 -0.2 -0.3 -0.4 -0.6 0.3 ]';
y0= [0 0.2 0.1 0.5 0.3 0.2 0 0 -0.2 -0.3 -0.4 -0.6 0.1]';
phi =[0 5 10 2 5 10 0 5 10 2 5 10 0]';
P = phantom([A a b x0 y0 phi],256);
%[P,E] = phantom([A a b x0 y0 phi],256);
imtool(P);

输出的自定义模体图像为：

可以根据自己的需求去修改模体参数矩阵E，达到想要的模体效果。