EGS中ctcreate的说明

zhuer91

希望下面的内容对大家有所帮助，需要提醒的一点是，注意执行文件的目录。
说明书Pirs794讲了ctcreate， 它可以读取CT片子，生成egsphant文件给dosxyz_nrc用作模体文件。所谓模体phantom，可以这样理解，一个立方体分成很多个小格子，每个小格子里面填充一种物质。

ctcreate的执行顺序为：
1.    读入CT文件
2.    读入CT header参数
3.    读入二进制CT数据
4.    选择CT数据的子集
5.    重新安排CT数据，得到满足要求的相应的需要记录剂量的区域
6.    将CT数据转换成物质和密度，对应每个小体元
7.    生成egsphant文件，给dosxyz¬_nrc用。


Ctcreate的功能就是：
1.    选择原CT片子的一个子集作为输出文件的整体轮廓。
2.    重新细分网格小体元，重新给小体元分配物质（密度），用作后续记录剂量的小区域。

以dicom格式的CT片子为例，
例题中的slice_names包含了10个CT文件的路径名，将里面的$OMEGA_HOME改为实际的目录，这里为/home/egs/HEN_HOUSE/omega。
CT_create_DICOM.inp为ctcreate的输入文件，里面记录了需要的数据。
第一行，CT片子类型，这里为DICOM。
第二行，CT文件的路径存放文件（slices_names）。
第三行，新的模体文件的xyz三个方向的边界。
第四行，重组模体的小体元的尺寸大小
第五行，新模体的物质密度选项，默认为0。

Ctcreate的运行有两种方式，读入输入文件×××.inp或是交互式输入。
1.    终端中输入ctcreate ×××.inp，回车。执行结束，生成了×××.egsphant文件，给dosxyz_nrc用。
2.    输入ctcreate，回车。
（1）    选择CT类型，1 pinnacle 2 CADPLAN 3 DICOM
（2）    输入DICOM文件名，这里为slice_names
（3）    给定新模体xyz三个方向的边界，单位cm，中间用逗号分隔。
（4）    输入新的小体元尺寸，单位cm。
（5）    输入media号，最大为7，如果使用hard-wired抽样函数，则输入0，一般输入0.
（6）    运行结束，生成了×××.egsphant文件，给dosxyz_nrc用。


转自 52MC论坛