一种用于突出显示体数据中目标区域的体绘制方法
【专利摘要】本发明涉及一种用于突出显示体数据中目标区域的体绘制方法，属于计算机图形学中科学可视化直接体绘制【技术领域】。首先计算原始数据的多个特征，对特征数据进行评估选择合适的特征组成传输函数。用归一化割方法分割二维直方图图像对目标数据进行分离。如果目标分离不够精细，就继续选择合适的特征对当前结果进一步分离，从而逐步得到精细的分离结果。之后利用分离结果合成最终的传输函数。最后在体绘制时利用传输函数突出显示目标区域，得到理想的可视化效果。本发明方法强化了传输函数的特征区分能力，使用了基于二维直方图图像分割和层次聚类的交互方式，操作便捷。而且适用于多种数据。
【专利说明】—种用于突出显示体数据中目标区域的体绘制方法
【技术领域】
[0001]本发明涉及一种用于突出显示体数据中目标区域的体绘制方法，属于计算机图形学中科学可视化直接体绘制【技术领域】。
【背景技术】
[0002]直接体绘制方法是探索三维图像整体信息和内部结构的一种有效手段。直接体绘制方法通过传输函数为不同的体素设置不同的颜色和不透明度，然后采用光学积分模型对体数据采样，得到绘制结果。其中传输函数的设计是决定体绘制结果好坏的关键。通过传输函数，用户可以从原始体数据中分离出用户感兴趣的区域，并对其突出显示。按使用的特征维度分类，传输函数的几种类型及相应存在的问题如下:
[0003](I)以灰度作为特征的一维传输函数是目前应用最广的传输函数，它具有良好的交互性和较好的可视化效果。但对于具有相同灰度的不同特征，一维灰度传输函数并不能很好地区分。
[0004](2) 二维传输函数在科研领域十分流行。它凭借特征维度的优势相比一维传输函数拥有更强的特征区分能力。但是将三维以上的体素信息投射到二维传输函数上也会引起维度丢失，而且二维传输函数的交互不直观，操作不便捷。
[0005](3)三维以上的传输函数具有更好的特征维度优势，特征区分能力更突出。但是其传输函数的交互非常复杂，无论在科研领域还是工业界都没有得到广泛的应用。
[0006]除此之外，还有很多应用于某种特定体数据的传输函数设计方法，如医学中用于突出血管的传输函数设计。也有很多强化交互模式的传输函数设计方法，如直接与体绘制结果交互的所见即所得的传输函数设计。
[0007]以上这些方法都有自己独特的优势，但不能同时具有简单的交互模式，强大的特征区分能力和多种数据应用范围这些特点。如何让传输函数对多种数据同时拥有强大的特征区分能力和交互的便捷性仍然是一个亟待解决的问题。
【发明内容】
[0008]本发明的目的是提出一种用于突出显示体数据中目标区域的体绘制方法，在二维传输函数的基础上，依次使用不同的特征对原始体数据进行处理，以求利用特征维度的拓展，逐步精确目标区域的分离。同时利用面向数据分类的特征空间度量差异性的评估方法，在分离目标区时下一二维传输函数特征维度的选取，使得选用特征之间的差异性较大，从而逐步构造合适的特征序列。
[0009]本发明提出的用于突出显示体数据中目标区域的体绘制方法，包括以下步骤:
[0010](I)分别计算体数据的梯度特征、曲率特征、分型特征、光滑度、线型特征和三维不变矩，具体过程如下:
[0011](1-1)设体数据中任意一个体素的梯度特征为G，G = |g|，G为g的模，g为体数据中体素的梯度向量，
【权利要求】
1.一种用于突出显示体数据中目标区域的体绘制方法，其特征在于该方法包括以下步骤: (I)分别计算体数据的梯度特征、曲率特征、分型特征、光滑度、线型特征和三维不变矩，具体过程如下:(1-1)设体数据中任意一个体素的梯度特征为G，G = |g|，G为g的模，g为体数据中体素的梯度向量，
f为体数据灰度特征在x、y、z三个方向上的函数映射，灰度特征为体素的数值特征
的转置，则体数据的梯度特征为体数据中所有体素的梯度特征按照体素在体数据中的位置排列得到的一个三维矩阵；(1-2)设体数据中任意一个体素的曲率特征为Cm,cn=k1+k2/2,K1为体素所在位置的主曲率，K 2为体素所在位置的副曲率
为几何张量ζ的迹，F为几何张量Z的弗罗贝尼乌斯范数，几何张量
其中扎是体素的海森矩阵，Pffl是体素所在位置的等值面切平面矩阵，Pffl = Ι-ηνη；, I为单位矩阵，ην为体素所在位置的表面向量，nv = η；为ην的转置，则体数据的曲率特征为体数据中所有体素的曲率特征按照体素在体数据中的位置排列得到的三维矩阵；(1-3)设体数据中任意一个体素的分型特征为D，对体数据中任意一个体素Vi,分型特征
为体素 Vi 的邻域，r 为邻域半径，f(vj)为体素Vj的灰度特征，Vj为体素Vi的邻域U(Vi，r)中任意体素，r))表示体素Vi的以r为半径的邻域中所有体素的灰度特征的和，kd的取值为2，则体数据的分型特征为体数据中所有体素的分型特征按照体素在体数据中的位置排列得到的三维矩阵；(1-4)设体数据中任意一个体素的光滑度为S，对体数据中任意一个体素Vi,光滑度5 = 1-$，其中σ2为体素Vi的邻域U(Vi，r)中所有体素的灰度特征方差，体素Vi的邻域半径r取2，则体数据的光滑度为体数据中所有体素的光滑度按照体素在体数据中的位置排列得到的三维矩阵；(1-5)设体数据中任意一个体素的线型特征为Fflit，对体数据中任意一个体素Vi，线型特征
，其中为权值，取值为任意非负实数，r为体素Vi的邻域半径
，其中5力~均一化后的
结果，
为体素Vi的邻域内以邻域边界为两端端点且通过Vi的线段，线段总数为12r2+l，'为线段Lk上的一个体素，体素Vi的邻域半径r取1，则体数据的线型特征为体数据中所有体素的线型特征按照体素在体数据中的位置排列得到的三维矩阵； (1-6)设体数据中任意一个体素的三维不变矩为τ，三维不变矩τ = μ 200+μ 020+μ 002^其中μ 200> μ 020和μ 002分别为中心距，对任意一个体素 '，中心距计算通式为:
其中Vj是体素Vi邻域U (Vi, r)内的任意一个体素，X1, X2, X3分别为体素Vj在体数据中三个维度的坐标值，石，石，石分别为体素Vi的邻域U(Vi，r)的灰度特征重心在体数据中三个维度的坐标值，f (x1； x2, X3)为Vj的灰度特征，体素Vi的邻域半径r取4,则体数据的三维不变矩为体数据中所有体素的三维不变矩按照体素在体数据中的位置排列得到的三维矩阵； (1-7)将原始体数据进行备份，备份的原始体数据记为Vcw ； (2)分别计算上述步骤(1)中得到的六个特征的特征评估量E，E=Q1Ln-W2R^w3Cn,其中ωι，ω2，ω3为权值，取值为任意非负实数，Ln、Rn、Cn分别为对特征距离评估量L、特征相关性评估量R、特征价值评估量C进行均一化处理的结果，L、R、C分别为:
'其中为用k均值聚类算法对特征Ρι聚类后的
第I1个类的重心，P1为上一次评估时选取的特征评估量E为最大的特征，若当前评估为第一次评估，则特EP1为灰度特征，h为用k均值聚类算法对特征P1聚类后得到的类的个数，为用k均值聚类算法对特征P2聚类后的第J1个类的重心，P2为待评估的特征，k2为用k均值聚类算法对待评估的特征P2聚类后得到的类的个数； R = H(P1)+H(P2)-H(P1, P2)，其中H(P1)为特征P1的信息熵，H(P2)为待评估特征P2的信息熵，H(P1, P2)为特征P1和待评估的特征P2的联合熵；
C = H(P17P2)-H(P1)； (3)从步骤(2)的六个特征中选取特征评估量E为最大的特征，记为P，用特征P与原始体数据的灰度特征组成一个二维传输函数，用该二维传输函数对原始体数据中的目标区域进行分离，得到原始体数据中位置在目标区域中的体素，分离过程如下: (3-1)用特征P与体数据的灰度特征组成一个二维统计直方图，用归一化割方法对该直方图图像进行分割，得到多个二维统计直方图图像子区域，记该多个子区域组成的集合为Y0 ； (3-2)计算原始体数据中所有体素在二维统计直方图中的坐标，根据该坐标，从原始体数据中分别找到与每个子区域相对应的体素，利用光线投射方法，分别对与每个子区域相对应的体素进行体绘制，得到多个体绘制图像； (3-3)从子区域集合Ytl中得到最终子区域集合Ofinal，过程如下: (3-3-1)初始化时，最终子区域集合为空集，对Ytl中所有子区域进行从I到w排序，其中w是子区域集合Ytl中子区域的数量； (3-3-2)从子区域集合Ytl中取出标号最小的子区域％，若与子区域％相对应的体绘制图像中只包含目标区域，则将10放入最终子区域集合中，并将10从Ytl中删除，将删除10后的\记为Y1 ;若与子区域％相对应的体绘制图像中只包含非目标区域，则将％从\中删除，将删除％后的Ytl记为Y1，若与子区域10相对应的体绘制图像中既包含目标区域又包含非目标区域，则用归一化割方法对子区域10进行分割，得到两个新的子区域Y1和12，分别将Y1和I2加入到子区域集合Ytl中，将L从Yci中删除，得到新的子区域集合Y1, Y1中Y1标号为Nmax+1, Y2标号为N_+2，其中Nniax为子区域集合Ytl中子区域标号的最大值； 重复上述步骤，若子区域集合Y变为空集，则记最终子区域集合为Ofinal，进行步骤(3-4);若上述步骤重复次数超过30次，则将Y3tl中剩余的所有子区域加入到最终子区域集合中，并记最终子区域集合为Ofinal，进行步骤(3-4)； (3-4)根据步骤(3-2)得到原始体数据中所有体素的坐标，分别分离出坐标在最终子区域集合Ofinal中的每个子区域的所有体素，所有体素组成的集合记为I，按照体素在原始体数据中的坐标，将\中的所有体素加入到原始体数据中，将加入\中的所有体素后的原始体数据记为\，根据下式对 \的灰度特征进行处理，得到Vn的新的灰度特征如下:
其中Vh为Vn中的任意一个体素，fold (Vh)为Vn中Vh的原来的灰度特征，fnew (Vh)为\中Vh的新的灰度特征； (3-5)对上述Vn中的体素进行判断，若Vn中的任意体素Vh满足Vh e Vp则进行步骤(4)，若Vn中存在一个体素vh，满足仏t片，则根据下式对步骤⑴的6种特征数据进行处理，得到与每种特征相对应的新的特征值:
其中Vh为Vn中任意一个体素，Pf-(Vh)表示Vh的原来特征值，
Pfeatutr—new(Vh)表不Vh的新特征值，返回步骤(2); (4)利用步骤(3)得到的体数据Vn和步骤(1)中的备份原始体数据Vmw，合成一个传输函数，利用该传输函数得到体绘制图像，具体过程如下: (4-1)将步骤(3)得到的体数据Vn中的体素归为一个类，记为Hitl; (4-2)根据备份的原始体数据Vrapy，将步骤(3)得到的体数据Vn中灰度特征变为O的体素的灰度特征值恢复为原来的灰度特征值，将恢复灰度特征值的体素组成一个集合，记为V。，用V。中所有体素的梯度特征和灰度特征组成一个二维统计直方图，采用归一化割方法将此直方图分割成e个子区域，分别记为a1、afae ； (4-3)计算上述V。中所有体素在梯度特征和灰度特征组成的二维统计直方图中的坐标，根据V。中所有体素在该二维统计直方图中的坐标，将V。划分为e个类，分别记为m1、m2...mb...me,b = I, 2,..., e,分类方法为:对V0中任意一个体素Vi,若Vi在该二维统计直方图中的坐标在子区域ab内部，则将Vi归到类mb中； (4-4)根据显示需要，为体数据Vn中的体素的类Hitl设置不透明度和颜色，为类设置不透明度和互相不同的颜色，不透明度取值范围为[O，I]，颜色采用红绿蓝颜色模型，红、绿、蓝颜色分量的强度值范围为[0，255]，根据体数据中每个类的颜色和不透明度，采用光线投射算法得到最终的突出显示目标区域的体绘制图像。
【文档编号】G06T7/00GK104167013SQ201410378697
【公开日】2014年11月26日 申请日期:2014年8月4日 优先权日:2014年8月4日
【发明者】陈莉, 孙思远 申请人:清华大学