基于目标cf的矢量量化压缩体绘制方法
【专利摘要】本发明提供一种基于目标CF的矢量量化压缩体绘制方法，包括：根据预先定义的可信度计算准则，对输入矢量的每一个点计算可信度，并统计对应各CF区间的比重；CF值高对应的优先级也高，CF值低对应的优先级也低；查看第一优先级的比重，若该优先级的比重低于第一阈值，继续，否则给该矢量最高优先级的CF值，结束；当第一优先级的比重过低，比较第二优先级的比重，若该比重低于第二阈值，继续，否则给该矢量第二优先级的CF值，结束；当第二优先级的比重仍不满足要求，则计算第三优先级，若比重低于第三阈值，那么赋给第四优先级的CF值，否则赋给第三优先级的CF值，利用得到的CF值，去影响失真测度计算，码书分裂偏向可信度高的区域。
【专利说明】基于目标CF的矢量量化压缩体绘制方法
【技术领域】
[0001]本发明属于体绘制领域，具体涉及一种CF矢量量化压缩体绘制。
【背景技术】
[0002]体绘制是一种直接由三维数据场产生二维图像的技术。由于该技术可以高效的呈现出数据的整体信息和内部细节，在众多领域得到广泛的应用。随着图形硬件的发展，基于GPU加速的体绘制方法成为目前体绘制领域的主流方法。压缩体绘制(compressed volumerendering,CVR)同时涉及体压缩和体绘制两种技术,是在二者紧密结合之下产生出的一种新的体绘制方法。压缩体绘制主要解决海量数据体绘制时面临主存和显存容量限制的问题。现有压缩体绘制方法大都采用一种均匀压缩或基于数值分布的非均匀的方式，但是在处理海量地震数据时，上述方式无法平衡失真度与压缩比之间的矛盾，要求高压缩比就导致高失真，重要信息损失严重。
[0003]施耐德等第一次使用GPU实现了分层矢量量化压缩体渲染。矢量量化第一次被引入体绘制，分层矢量量化(hierarchical vector quantizat1n HVQ)是一种改进的VQ策略。该方法通过引入对原始数据的多分辨率的协方差分析，来获得一种能有效处理多为数据的VQ方法。赵利平等针对地震数据领域具有的情境，引入了分类的思想对提出分层分类矢量量化压缩体绘制方法。将矢量量化压缩引入地震数据处理领域，获得较优秀的成果。
[0004]在实际地震数据应用中，存在另一种应用情景。在该应用场景下，目标不能通过一个区域来描述，需要给出更明确的描述方式。例如地震数据中的某些微构造，对于这种存在于数据体中的重要信息，其往往具有一定的特征，但是对这种特征的描述需要引入行业内的经验判断。由于引入了经验判断，这种描述就存在一定的不确定性，基于此本发明提出基于目标可信度的矢量量化压缩体绘制方法。通过可信度的概念来对无法精确描述的目标进行建模。
[0005]矢量量化在压缩性能上要优于标量量化，已被广泛的应用在数据压缩领域，特别是三维体数据压缩方面。对于三维体数据，矢量量化中输入矢量来自于对原始数据的分块，对于三维体数据，划分出来的小块的数量通常非常巨大。应用传统方法，由于巨大的输入矢量集合，使得码书的训练生成，码字的搜索匹配变得异常的困难，因此现有的对体数据进行矢量量化的方法都是在传统方法基础之上进行的改进。
[0006]LBG方法是最经典的矢量量化方法，该方法通过给定一个起始的码书，然后通过不断迭代的方式，直到满足用户设定的失真测度。流程图如图1所示。
[0007]传统意义的LBG方法主要有以下几个缺点:
[0008]在根据现有码本对输入矢量进行最佳划分时，所需要的计算量太大，时间和空间复杂度都很高，导致方法的实用性过低。
[0009]该方法在每次改变形心后，进行最佳划分时，矢量的变化不会很大，通常只会在小范围内变化，导致最后的量化结果满足局部最优。因而对初始码本的选定较为依赖。
[0010]因为LBG方法对初始码本的依赖较为严重，学者们研究提出了几种初始码本的选取方法:
[0011 ] 随机选取法，该方法具有一定的随机性，来构成初始码书
[0012]^^^'^^…，^，此方式特点在于^间复杂度低^且没有空胞腔问题。
[0013]分裂法，从初始只有一个码字，是通过计算所有输入矢量的形心来计算得至IJ。然后用一个适当的系数A，乘以码字Vt"，形成第二个码字Ctn。以为简单的初始码书，运用LBG方法设计码书，然后根据乘以合适的系数B，来确定4个码字
广吨⑼，然后以这四个去构建码书，依次继续进行，直到码书个数达到要求。
【发明内容】
[0014]为解决现有技术存在的上述问题，本发明提出一种基于目标CF的矢量量化压缩体绘制方法。
[0015]本发明具体的技术方案为:基于目标CF的矢量量化压缩体绘制方法，CF计算的具体包括以下步骤:
[0016]S1:根据预先定义的可信度计算准则，对输入矢量的每一个点计算其可信度，并统计对应各个CF区间的比重；
[0017]S2:CF的优先级主要根据专家对可信度准则的定义来给出，也即CF值高的其对应的优先级同样高，CF值低的其对应的优先级也低；
[0018]S3:查看第一优先级的比重，如果该优先级的比重低于第一阈值，那么继续，否则给该矢量最高优先级的CF值，结束；
[0019]S4:当第一优先级的比重过低时，比较第二优先级的比重，如果该比重低于第二阈值时，继续，否则给该矢量第二优先级的CF值，结束；
[0020]S5:当第二优先级的比重依然不能满足要求时，则计算第三优先级，如果比重低于第三阈值，那么赋给第四优先级的CF值，否则赋给第三优先级的CF值；
[0021]S6:利用计算出来的CF值，作为可信度的因子，乘以失真测度，得到可信度的值，进行码书分裂时，选取可信度高的块进行分裂；
[0022]进一步地，所述第一阈值为10% ；
[0023]进一步地，所述第二阈值为20% ;
[0024]更进一步地，所述第三阈值为30%。
[0025]所述码书设计具体包括:
[0026]对于一个输入的单一量化单元V1，该量化单元包含了所有的输入矢量集合I，设Y1为该集合的形心，也就是对应的码书，计算当前量化单元的失真测度A = Σλ%.，
[0027]其中，λ为计算出来的CF值，=为均方误差欧式距离；
k i=i
[0028]计算得到失真测度后，接着构造一个双向链表，将
[0029](Di Y飞=(/_e ^F1))插入到链表中，然后对链表中的元素从大到小进行排序，每一次对失真测度I最大的量化单元进行分裂，每次分裂的步骤如下:
[0030]S61.在链表中选择包含最大失真测度Dj的组；
[0031]S62.计算自相关协方差矩阵
【权利要求】
1.基于目标CF的矢量量化压缩体绘制方法,其特征在于,包括: 51:根据预先定义的可信度计算准则，对输入矢量的每一个点计算其可信度，并统计对应各个CF区间的比重； 52:CF的优先级主要根据专家对可信度准则的定义来给出，也即CF值高的其对应的优先级同样高，CF值低的其对应的优先级也低； S3:查看第一优先级的比重，如果该优先级的比重低于第一阈值，那么继续，否则给该矢量最高优先级的CF值，结束； 54:当第一优先级的比重过低时，比较第二优先级的比重，如果该比重低于第二阈值时，继续，否则给该矢量第二优先级的CF值，结束； 55:当第二优先级的比重依然不能满足要求时，则计算第三优先级，如果比重低于第三阈值，那么赋给第四优先级的CF值，否则赋给第三优先级的CF值。 S6:利用计算出来的CF值，作为可信度的因子，乘以失真测度，得到可信度的值，进行码书分裂时，选取可信度高的块进行分裂。
2.根据权利要求1基于目标CF的矢量量化压缩体绘制方法，其特征在于，CF通过如下形式给出:
其中，CF表示当数据值X在(an，bn]范围内的时候，判定其为目标的可信程度λη。
3.根据权利要求1基于目标CF的矢量量化压缩体绘制方法，其特征在于，所述第一阈值为10%。
4.根据权利要求1基于目标CF的矢量量化压缩体绘制方法，其特征在于，所述第二阈值为20%。
5.根据权利要求1基于目标CF的矢量量化压缩体绘制方法，其特征在于，所述第三阈值为30%。
6.根据权利要求1基于目标CF的矢量量化压缩体绘制方法，其特征在于，步骤S6所述码书的设计具体包括:对于一个输入的单一量化单元V1，该量化单元包含了所有的输入矢量集合1，设Y1为该集合的形心，所述形心即对应的码书，计算当前量化单元的失真测度
.其中，λ为计算出来的CF值，
力均方误差欧式距离； 计算得到失真测度后，接着构造一个双向链表，将 (D1，3, ) = (/￡ 3 ^插入到链表中，然后对链表中的元素从大到小进行排序，每一次对失真测度h最大的量化单元进行分裂； 其中，$表示一个分组里面量化单元的编号。
7.根据权利要求6基于目标CF的矢量量化压缩体绘制方法，其特征在于，所述分裂的具体步骤如下: S61:在链表中选择包含最大失真测度Dj的组； S62:计算自相关协方差矩阵:
S63:计算M的最大特征向量和与其对应的最大特征值λ max ； S64:将该组分为 Ileft 和 3 ,两个组，
S65:对于新分出来的组3〃,和，分别计算他们的形心Ylrft和YHght，以及各自新的失真测度Dleft和Dright ； S66:将得到的两个组插入到链表中； S67:如果组数的数目等于2%也就是码书的数目等于Ar自己给出的一个值，则停止，否则跳转到步骤S61。
【文档编号】G01V1/32GK104200510SQ201410416240
【公开日】2014年12月10日 申请日期:2014年8月22日 优先权日:2014年8月22日
【发明者】鲁才, 唐志梁, 胡光岷 申请人:电子科技大学