1.一种血管管腔结构重建方法,其特征在于,包括有如下步骤:
步骤S1,通过三维影像数据采集方式获取血管的体数据文件;
步骤S2,从体数据文件中提取血管文件,所述血管文件包括中心线点坐标和血管半径;
步骤S3,将所述血管文件中的数据由长度单位转换为体素单位;
步骤S4,构建一个呈三维矩阵的采样空间,所述采样空间中每个采样点的数据包括以体素为单位的中心线点坐标和血管半径;
步骤S5,计算采样空间范围内每一采样点的灰度值grey_scale_value:其中,center_x、center_y和center_z表示当前采样点的坐标,radius表示当前采样点的血管半径,乘数1000表示将计算结果放大1000倍;
步骤S6,利用VTK系统将采样空间中的数据在该vtk系统的渲染窗口中进行可视化体绘制。
2.如权利要求1所述的血管管腔结构重建方法,其特征在于,所述步骤S1中的体数据文件包括有mhd格式文件和raw格式文件,所述步骤S2中的血管文件是txt文本格式文件。
3.如权利要求2所述的血管管腔结构重建方法,其特征在于,所述血管文件中包括五列数据,其中第一至第三列数据表示中心线点的X轴、Y轴和Z轴坐标,第四列为中心线点的半径,第五列为误差。
4.如权利要求2所述的血管管腔结构重建方法,其特征在于,所述mhd格式文件中包括有转换系数ElementSpacing,所述转换系数ElementSpacing与中心线点一一对应,所述步骤S3中由长度单位转换为体素单位的过程包括:将当前中心线点的坐标值(x、y、z)除以对应的转换系数ElementSpacing。
5.如权利要求1所述的血管管腔结构重建方法,其特征在于,所述步骤S3与步骤S4之间还包括有线性插值计算步骤:根据预设数量的体素单位内中心线点的坐标,计算得出插值点坐标coordinate:其中,start_coordinate表示开始点的坐标,end_coordinate表示结束点的坐标,accumulate_dis表示由开始点至结束点的长度累加。
6.如权利要求5所述的血管管腔结构重建方法,其特征在于,预设数量的体素单位为1个体素单位。
7.如权利要求2所述的血管管腔结构重建方法,其特征在于,所述mhd格式文件中包括有三维空间的尺寸数据,所述步骤S4中,按照该三维空间的尺寸数据构建采样空间。
8.如权利要求1所述的血管管腔结构重建方法,其特征在于,所述步骤S4中,在所述采样空间中建立子采样空间,所述子采样空间以当前中心线点坐标为中心,并利用步骤S5中的计算方法,得出子采样空间中各采样点的灰度值。
9.如权利要求8所述的血管管腔结构重建方法,其特征在于,所述步骤S4中,当所述子采样空间到达或者超过采样空间的边界时,终止步骤S5,结束对采样点灰度值的计算。
10.如权利要求8所述的血管管腔结构重建方法,其特征在于,所述子采样空间中,中心线点坐标至子采样空间边缘的距离大于所述血管半径的最大值。
11.如权利要求1所述的血管管腔结构重建方法,其特征在于,所述步骤S6中,进行可视化体绘制之前包括灰度值过滤步骤:以灰度值45作为一个临界值进行过滤,将灰度值小于
45的采样点设置为透明。
12.如权利要求11所述的血管管腔结构重建方法,其特征在于,所述灰度值大于45的采样点,其透明度为:0<透明度≤0.8。
13.如权利要求1所述的血管管腔结构重建方法,其特征在于,所述步骤S6中还包括有上色过程:利用VTK系统中的vtkColorTransferFunction函数对体数据进行上色。