[发明专利]一种三维地震体矢量量化压缩体绘制方法

权利要求书

1.一种三维地震体矢量量化压缩体绘制方法，其特征在于，该方法为：

1)将原始地震数据分割成大小为N*N*N大小的数据块；

2)计算每一个数据块的梯度值，根据给定的阈值，将梯度值低于所述阈值的数据块标记为零块，并计算所述零块的均值；将梯度值不低于所述阈值的数据块标记为非零块；所述阈值的绝对值小于100；

3)对每一个数据块中相邻的（N/2）*（N/2）*（N/2）的数据求平均值，获得大小为（N/2）*（N/2）*（N/2）的小数据块；计算每一个小数据块的均值；

4)对所述数据块与所述小数据块进行差分处理，得到最高层数据，即大小为N*N*N的块差分数据；对所述小数据块与所述小数据块的均值进行差分处理，得到次高层数据，即大小为（N/2）*（N/2）*（N/2）的差分矢量；将所述最高层数据经矢量量化后，得到码书一；将所述次高层数据经矢量量化后，得到码书二；

5)将上述步骤4）处理后得到的N*N*N的块差分数据或（N/2）*(N/2)*(N/2)的差分矢量作为输入矢量I，将所述输入矢量I的集合作为一个量化单元V₁，计算量化单元V₁的失真度D₁，进入步骤6）；V₁的形心矢量Y₁就是块差分数据对应的码书一的一个码字或差分矢量对应的码书二的一个码字；

6)根据上述量化单元V₁构造双向链表，所述双向链表中每一个分组对应步骤5）的一个码字，遍历所述双向链表，将所述双向链表分裂为若干个分组；

7)确定所述双向链表每一个分组的形心，所述形心即为其对应的分组的码字；根据码字的个数，将所有码字重新排列为一个体数据codebook；

8)为所述双向链表每一个分组中的输入矢量赋给代表所述输入矢量的码字在码书一或码书二中的索引，计算所述分组对应码字在codebook中的索引，并建立索引体数据，其中码书一对应索引数据体一，码书二对应索引数据体二，两个索引数据体中均包含输入矢量对应的码字在码书一或码书二中的索引信息，索引数据体二中还包括了每一个输入矢量的对应的数据块或小数据块的类别信息，即零块或非零块；

9)确定所述最高层数据或次高层数据中每一个数据块或小数据块在所述码书一或码书二中的索引值信息，根据索引数据体二中的数据块或小数据块的类别信息分别处理非零块和零块：对于零块，直接用步骤2）中计算的均值代替该数据块的标量值；对于非零块，根据索引信息从码书一中查找最高层数据对应的码字，再从码书二中查找次高层数据对应的码字，最后将上述两个码字与该非零块对应的均值相加，即为该非零块的标量值；所述标量值的取值范围是0～255；根据标量值的范围，对每一个标量值值赋给设定颜色信息，然后用这组颜色信息生成地震数据体色标纹理；

10)根据上述步骤9）计算的标量值，查找原始地震数据的色表纹理，获取对应的颜色信息，完成绘制。

2.根据权利要求1所述的三维地震体矢量量化压缩体绘制方法，其特征在于，所述步骤1）中，所述数据块的大小为4*4*4。

3.根据权利要求1所述的三维地震体矢量量化压缩体绘制方法，其特征在于，所述步骤2）中，所述给定的阈值的取值范围为0～20。