基于自动多阈值的多层着色体绘制方法技术

本发明专利技术公开了一种基于自动多阈值的多层着色体绘制系统，该系统包括：阈值生成模块，根据输入的体素样本生成阈值序列；绘制模块，根据阈值序列对于体数据进行多层着色，生成着色后的输出图像。本发明专利技术同时还提出了一种基于自动多阈值的多层着色体绘制方法。与已有技术相比，本发明专利技术对多层结构的区分效果更加明显，可以对体数掘进行分层着色，且不需要提前分割数据，绘制过程中的内存和计算代价小，效率高。本发明专利技术在科学计算可视化和医学影像可视化领域有重要的应用价值。

【技术实现步骤摘要】
基于自动多阈值的多层着色体绘制方法
本专利技术属于计算机图形与可视化
，尤其涉及一种基于自动多阈值的多层着色体绘制方法，可用于三维标量体数据的可视化。
技术介绍
三维标量体数据的体绘制技术可以帮助人们观察、挖掘三维标量体数据(如CT切片数据)中包含的信息，在科学计算可视化和医学影像可视化领域有着重要的应用价值。三维标量体数据的体绘制有多种实现方式，从计算顺序上分，有物体顺序的体素投影法(可参考文献“Rottgeretal.2000.Hardware-AcceleratedVolumeAndIsosurfaceRenderingBasedOnCell-Projection”)和图像顺序的光线跟踪法(可参考文献“Parkeretal.2005.Interactiveraytracingforvolumevisualization”)等。在三维标量体数据的体绘制中，对不同结构的分类和着色是一个至关重要的技术。由于从三维体空间投影到二维平面会产生空间重叠，因此必须根据需要对各种结构通过透明度的设定进行取舍，并通过颜色的设定进行标记。为了区分不同的结构，现有的方法主要采用传递函数和分割两种技术。采用传递函数的方法是直接根据体数据在空间各点的取值设定透明度和颜色，或者结合局部特征或预先计算的特征信息在二维空间进行透明度和颜色的设定，可参考文献“Kindlmannetal.1998.Semi-automaticgenerationoftransferfunctionsfordirectvolumerendering”和“Haidacheretal.2010.Volumevisualizationbasedonstatisticaltransfer-functionspaces”等。采用传递函数的方法由于主要使用邻域信息，对全局结构，尤其是多层结构的区分效果不理想，有些较复杂的预处理算法耗时较长。采用分割技术的方法是通过一个额外的体数据对结构进行明确的标记，在绘制的过程中根据结构信息选择不同的传递函数，可参考文献“Hadwigeretal.2003.High-qualitytwo-levelvolumerenderingofsegmenteddatasetsonconsumergraphicshardware”和“Xiangetal.2011.Skeleton-cutsanefficientsegmentationmethodforvolumerendering”等。采用分割技术的方法需要复杂的预处理步骤，计算代价大，同时给绘制过程带来额外的空间和时间开销。
技术实现思路
基于以上现有技术中存在的问题，本专利技术提出了一种基于自动多阈值的多层着色体绘制方法，该方法根据输入体素样本自动确定多个阈值，在体绘制过程中用这些阈值将光线分为多个区段，并在每个区段应用不同的传递函数实现对体数据的多层着色显示。根据本专利技术的一方面，提出一种基于自动多阈值的多层着色体绘制系统，该系统包括：阈值生成模块和绘制模块，其中：所述阈值生成模块用于根据输入的体素样本生成阈值序列；所述绘制模块用于根据所述阈值序列对于输入的体数据进行多层着色，生成着色后的输出图像。根据本专利技术的另一方面，还提出一种基于自动多阈值的多层着色体绘制方法，该方法包括以下步骤：步骤1，根据用户输入的体素样本生成阈值序列；步骤2，根据所述阈值序列对于输入的体数据进行多层着色，生成着色后的输出图像；所述步骤1进一步包括以下步骤：步骤11，输入体素样本；步骤12，根据所述体素样本生成光线样本；步骤13，根据所述光线样本生成目标函数，并根据所述目标函数得到最优解区域；步骤14，对于所述最优解区域进行距离变换得到内部极大值点，该点的两维坐标即为所求的一组阈值；所述步骤2进一步包括以下步骤：步骤21，输入体数据；步骤22，根据所述阈值序列生成状态信息，并根据所述状态信息对于光线进行分段；步骤23，根据分段信息对于所述体数据进行多层着色，得到颜色值；步骤24，对每条光线进行遍历，得到一组样本颜色值，将其进行混合，并将混合结果写入输出图像；步骤25，为所述输出图像中的每一个像素生成一条光线，基于各条光线的遍历和混合结果，得到输出图像。本专利技术与现有技术相比所取得的有益效果为：1、本专利技术和基于传递函数的方法相比，其对全局结构，尤其是多层结构的区分效果更加明显，可以对体数据进行分层着色；2、本专利技术和基于分割技术的方法相比，虽然同样使用多个传递函数，但是不需要计算和保存标记信息，计算代价和储存代价更低。同时，由于不使用静态标记，其对动态数据具有更好的适应性。附图说明图1为本专利技术基于自动多阈值的多层着色体绘制系统的结构示意图；图2示意了由输入的体素样本生成随机光线样本的过程；图3示意了由随机光线样本生成目标函数，得到最优解区域，进而得到内部极大值点的过程；图4给出了加速算法对光线进行扫描的过程；图5给出了在体绘制过程中用阈值对光线分段的过程；图6给出了基于自动多阈值的多层着色体绘制方法的绘制效果。具体实施方式为使本专利技术的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本专利技术进一步详细说明。图1为本专利技术提出的基于自动多阈值的多层着色体绘制系统的结构示意图，如图1所示，所述系统包括阈值生成模块和绘制模块，其中：所述阈值生成模块用于根据用户输入的体素样本生成阈值序列，所述阈值序列中，每两个阈值为一组，每组阈值中含有一个高阈值和一个低阈值；所述绘制模块用于根据所述阈值序列对于输入的体数据进行多层着色，生成着色后的输出图像。进一步地，对于每组阈值，所述阈值生成模块包括：体素样本输入子模块、光线样本生成子模块、目标函数生成子模块和距离变换子模块，其中：所述体素样本输入子模块用于通过交互界面由用户输入体素样本；具体地，所述体素样本输入子模块接收用户通过鼠标、键盘、触摸屏等输入设备输入的指令，在原始三维标量体数据中标记出一系列分属于不同结构的体素样本，如图2A中由白线所围成的不同区域所示。所述光线样本生成子模块用于根据所述体素样本生成光线样本；具体地，所述光线样本生成子模块根据所述体素样本生成指定数量(比如1000)的随机光线作为光线样本，其中，每条光线穿过一个随机选定的体素样本，其方向在单位球面上随机选择，且所述光线的方向的概率密度满足均匀分布，如图2B所示。所述目标函数生成子模块用于根据所述光线样本生成目标函数，并根据所述目标函数得到最优解区域；所述目标函数是通过统计各个可选值对光线样本的正确分类数得到的，所述目标函数定义在所要求取的阈值序列中每组两个阈值的取值范围所张成的二维平面上，通过对于所述二维平面的离散化，可以确定有限个可选值(在本专利技术一实施例中，采用256x256的网格对所述二维平面进行离散化)。在不采用加速算法的情况下，可以对每个可选值依次进行测试，即基于每个可选值依次扫描每条随机光线，检测该光线到达体素样本所在位置时的状态是否正确，并统计正确率，得到的正确分类数作为所述目标函数的值，效果如图3所示。这种实现计算量较大，因此本专利技术不采用这种方法，而是采用加速算法，用较小的代价达到同样的效果。加速算法不对可选值进行测试，而是只对每条随机光线进行正向和反向两次扫描，在这两次扫本文档来自技高网...

【技术保护点】
一种基于自动多阈值的多层着色体绘制系统，其特征在于，该系统包括：阈值生成模块和绘制模块，其中：所述阈值生成模块用于根据输入的体素样本生成阈值序列；所述绘制模块用于根据所述阈值序列对于输入的体数据进行多层着色，生成着色后的输出图像。

【技术特征摘要】
1.一种基于自动多阈值的多层着色体绘制系统，其特征在于，该系统包括：阈值生成模块和绘制模块，其中：所述阈值生成模块用于根据输入的体素样本生成阈值序列；所述绘制模块用于根据所述阈值序列对于输入的体数据进行多层着色，生成着色后的输出图像；所述阈值生成模块包括：体素样本输入子模块、光线样本生成子模块、目标函数生成子模块和距离变换子模块，其中：所述体素样本输入子模块用于通过交互界面由用户输入体素样本；所述光线样本生成子模块用于根据所述体素样本生成光线样本；所述目标函数生成子模块用于根据所述光线样本生成目标函数，并根据所述目标函数得到最优解区域；所述距离变换子模块用于对于所述最优解区域进行距离变换得到内部极大值点，该点的两维坐标即为所求的一组阈值；其中，所述目标函数是通过统计各个可选值对光线样本的正确分类数得到的，所述目标函数定义在所要求取的阈值序列中每组两个阈值的取值范围所形成的二维平面上，通过对于所述二维平面的离散化，确定有限个可选值，采用加速算法实现；加速算法不对可选值进行测试，而是只对每条随机光线进行正向和反向两次扫描，在这两次扫描中确定能够对该光线进行正确分类的可选值区域，在前向扫描时计算最大值函数vmax(t)，在反向扫描时计算最小值函数vmin(t)，在函数值f(t)的每个下降区间的起点位置处取最大值vmax和最小vmin值，则所形成的三角形区域，所围成的区域内的可选值会对此光线做出“可以通过”的分类判断，如果这种分类判断是正确的，则对目标函数在此区域内进行加一操作，否则进行减一操作，最后得到的目标函数值也代表正确率。2.根据权利要求1所述的系统，其特征在于，所述阈值序列中，每两个阈值为一组，每组阈值中含有一个高阈值和一个低阈值。3.根据权利要求1所述的系统，其特征在于，所述光线样本生成子模块根据所述体素样本生成指定数量的随机光线作为光线样本，其中，每条光线穿过一个随机选定的体素样本，其方向在单位球面上随机选择，且所述光线的方向的概率密度满足均匀分布。4.根据权利要求1所述的系统，其特征在于，所述绘制模块包括：体数据输入子模块、分类器子模块、状态控制子模块、光线投射子模块、光线生成子模块，其中：所述体数据输入子模块用于通过交互界面由用户输入体数据；所述状态控制子模块用于根据所述阈值序列生成状态信息，并根据所述状态信息对于光线进行分段；所述分类器子模块用于根据所述分段信息对于所述体数据进行多层着色，并将得到的颜色值发送给所述光线投射子模块；所述光线投射子模块用于基于所述分类器子模块的输出，对每条光线进行遍历，得到一组样本颜...

【专利技术属性】
技术研发人员：田捷，杨飞，李秀丽，杨凤，杨彩云，
申请(专利权)人：中国科学院自动化研究所，
类型：发明
国别省市：