CN114189692B - 基于连续子空间图变换的点云属性编码方法及解码方法 - Google Patents
基于连续子空间图变换的点云属性编码方法及解码方法 Download PDFInfo
- Publication number
- CN114189692B CN114189692B CN202010966930.8A CN202010966930A CN114189692B CN 114189692 B CN114189692 B CN 114189692B CN 202010966930 A CN202010966930 A CN 202010966930A CN 114189692 B CN114189692 B CN 114189692B
- Authority
- CN
- China
- Prior art keywords
- attribute
- point cloud
- target point
- coefficient
- level
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 77
- 230000009466 transformation Effects 0.000 title claims abstract description 73
- 238000010586 diagram Methods 0.000 title description 19
- 239000011159 matrix material Substances 0.000 claims description 24
- 238000013139 quantization Methods 0.000 claims description 21
- 239000013598 vector Substances 0.000 claims description 20
- 230000001131 transforming effect Effects 0.000 claims description 9
- 230000000694 effects Effects 0.000 description 9
- 238000012545 processing Methods 0.000 description 4
- 238000004590 computer program Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
- H04N19/30—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using hierarchical techniques, e.g. scalability
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
- H04N19/10—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding
- H04N19/169—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding characterised by the coding unit, i.e. the structural portion or semantic portion of the video signal being the object or the subject of the adaptive coding
- H04N19/18—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using adaptive coding characterised by the coding unit, i.e. the structural portion or semantic portion of the video signal being the object or the subject of the adaptive coding the unit being a set of transform coefficients
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Compression Or Coding Systems Of Tv Signals (AREA)
- Image Processing (AREA)
Abstract
本发明公开了基于连续子空间图变换的点云属性编码方法以及解码方法,通过将目标点云划分为多层级结构,对多层级结构的第二级以及第二级以上的目标点云空间,都由该目标点云空间在下一级的子空间的第一属性系数进行图变换得到该节点的第一属性系数和第二属性系数,目标点云空间的第一属性系数在上一级继续被变换,这样进行连续变换后,最终得到所述目标点云的第一属性系数和各个目标点云空间的第二属性系数,再根据所述目标点云的第一属性系数和各个目标点云空间的第二属性系数进行点云属性编码,在编码的过程中,不涉及寻找变换点的过程,避免了寻找方向的顺序影响点云属性编码效果,保证了点云属性编码质量。
Description
技术领域
本发明涉及点云处理技术领域,特别涉及一种基于连续子空间图变换的点云属性编码方法及解码方法。
背景技术
三维点云是现实世界数字化的重要表现形式。随着三维扫描设备(如激光、雷达等)的快速发展,点云的精度和分辨率变得更高。高精度点云广泛应用于城市数字化地图的构建,在智慧城市、无人驾驶、文物保护等众多热门研究中起技术支撑作用。点云是三维扫描设备对物体表面采样所获取的,一帧点云的点数一般是百万级别,其中每个点包含几何信息和颜色、反射率等属性信息,数据量十分庞大。三维点云庞大的数据量给数据存储、传输等带来了巨大挑战,因此对点云进行压缩编码变得十分重要。
在现有技术中,点云编码常使用RAHT(基于区域自适应的分层变换)方法进行点云属性编码,但是,RAHT作为一种多层次分层变换,每次变换沿着x,y和z方向,寻找两个变换点进行变换,这就会涉及到寻找方向的顺序问题,不是最佳解的寻找方向会导致点云编码效果差。
因此,现有技术还有待改进和提高。
发明内容
针对现有技术的上述缺陷,本发明提供基于连续子空间图变换的点云属性编码方法及解码方法,旨在解决现有技术中采用RAHT方法进行点云属性编码时寻找方向的顺序影响点云编码效果的问题。
为了解决上述技术问题,本发明所采用的技术方案如下:
本发明的第一方面,提供一种基于连续子空间图变换的点云属性编码方法,包括:
根据目标点云内点的位置将所述目标点云划分为多层级结构,所述多层级结构包括M个层级,每级中包括至少两个点云空间,每一级的点云空间是对上一级的点云空间进行划分得到的,其中,M为大于1的整数;
对于所述多层级结构中的每一个目标点云空间,获取所述目标点云空间的子空间,根据所述子空间构建图模型,其中,所述目标点云空间为所述多层级结构中第一级之外的层级中的点云空间;
根据所述图模型对所述子空间的第一属性系数进行变换,获取所述目标点云空间的第一属性系数和第二属性系数;
根据所述多层级结构中的第M级中的点云空间的第一属性系数、以及所述多层级结构中各个所述目标点云空间的第二属性系数对所述目标点云进行属性编码;
其中,所述多层级结构中的第一级中的点云空间的第一属性系数为点云空间的原始属性值。
本发明的第二方面,提供一种基于连续子空间图变换的点云属性解码方法,包括:
根据目标点云内点的位置将所述目标点云划分为多层级结构,所述多层级结构包括M个层级,每级中包括至少两个点云空间,每一级的点云空间是对上一级的点云空间进行划分得到的,其中,M为大于1的整数;
对于所述多层级结构中的每一个目标点云空间,获取所述目标点云空间的子空间,根据所述子空间构建图模型,其中,所述目标点云空间为所述多层级结构中第一级之外的层级中的点云空间;
对所述目标点云的编码码流进行解码,获取所述多层级结构的第M级中的点云空间的第一属性系数和各个所述目标点云空间分别对应的第二属性系数,并根据所述图模型对所述目标点云空间的第一属性系数和第二属性系数进行变换,获取各个所述子空间的第一属性系数;
根据所述多层级结构中的第一级中各个点云空间的第一属性系数获取所述目标点云的属性值。
本发明的第三方面,提供一种终端,所述终端包括处理器、与处理器通信连接的存储介质,所述存储介质适于存储多条指令,所述处理器适于调用所述存储介质中的指令,以执行实现上述的基于连续子空间图变换的点云属性编码方法的步骤和/或上述的基于连续子空间图变换的点云属性解码方法的步骤。
本发明的第四方面,提供一种存储介质,所述存储介质存储有一个或者多个程序,所述一个或者多个程序可被一个或者多个处理器执行,以实现上述的基于连续子空间图变换的点云属性编码方法的步骤和/或上述的基于连续子空间图变换的点云属性解码方法的步骤。
与现有技术相比,本发明提供了基于连续子空间图变换的点云属性编码方法及解码方法,通过将目标点云划分为多层级结构,对多层级结构的第二级以及第二级以上的目标点云空间,都由该目标点云空间在下一级的子空间的第一属性系数进行图变换得到该节点的第一属性系数和第二属性系数,目标点云空间的第一属性系数在上一级继续被变换,这样进行连续变换后,最终得到所述目标点云的第一属性系数和各个目标点云空间的第二属性系数,再根据所述目标点云的第一属性系数和各个目标点云空间的第二属性系数进行点云属性编码,在编码的过程中,不涉及寻找变换点的过程,避免了寻找方向的顺序影响点云属性编码效果,保证了点云属性编码质量。
附图说明
图1为本发明提供的基于连续子空间图变换的点云属性编码方法的实施例的流程图;
图2为本发明提供的基于连续子空间图变换的点云属性编码方法的实施例中多层级结构的示意图;
图3为本发明提供的基于连续子空间图变换的点云属性编码方法的实施例中图模型的示意图;
图4为本发明提供的基于连续子空间图变换的点云属性编码方法的实施例的效果示意图;
图5为本发明提供的基于连续子空间图变换的点云属性解码方法的实施例的流程图;
图6为本发明提供的终端的实施例的原理示意图。
具体实施方式
为使本发明的目的、技术方案及效果更加清楚、明确,以下参照附图并举实施例对本发明进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
本发明提供的基于连续子空间图变换的点云属性编码方法及解码方法,可以是应用在终端中,终端可以通过本发明提供的基于连续子空间图变换的点云属性编码方法及解码方法进行点云属性编码和点云属性解码。终端可以但不限于是各种计算机、手机、平板电脑、车载电脑和便携式可穿戴设备。
实施例一
如图1所示,本发明提供的基于连续子空间图变换的点云属性编码方法的一个实施例中,包括步骤:
S110、根据目标点云内点的位置将所述目标点云划分为多层级结构。
所述多层级结构包括M个层级,每级中包括至少两个点云空间,每一级的点云空间是对上一级的点云空间进行划分得到的,其中,M为大于1的整数。将目标点云看做是一个根空间,可以根据空间位置将所述目标点云划分为多个点云空间,每个点云空间可以继续划分为下一级的点云空间。如图2所示,图2中是一个4级的多层级结构,G01为给多层级结构中的所述目标点云,位于多层级结构的第4级。
具体地,所述所述根据目标点云内点的位置将所述目标点云划分为多层级结构,包括:
对点云空间进行划分,得到K个子空间,每个子空间的空间大小相等或包含的数据点数量相等,每个子空间内包括至少一个所述目标点云的数据点,其中,K为大于1的整数;
当点云空间内包含的数据点的数量为1时,不再继续划分。
也就是说,在所述多层级结构中的第1级中的点云空间中包含1个数据点。在每次对点云空间进行划分,得到K个子空间时,K的数量可以自行设定。
S120、对于所述多层级结构中的每一个目标点云空间,获取所述目标点云空间的子空间,根据所述子空间构建图模型。
所述目标点云空间为所述多层级结构中第一级之外的层级中的点云空间,也就是说,位于所述多层级结构的第2级及第2级以上的层级中的点云空间均为目标点云空间。所述目标点云空间的子空间时对所述目标点云空间进行划分后得到的点云空间,对于每个目标点云空间,通过下面的步骤S300和S400进行处理,获取该目标点云空间对应的第一属性系数和第二属性系数,具体将在后文被说明。不难看出,所述目标点云也是目标点云空间。具体地,所述根据所述子空间构建图模型包括:
将所述子空间作为所述图模型中的顶点,将所述子空间两两连接的连接线作为所述图模型中的边。也就是说,当所述目标点云空间有N个所述子空间,那么,所述图模型中的顶点的个数为N,每个顶点都分别通过边与其他的各个顶点连接,值得说明的是,在所述图模型中,认为每个顶点还存在和自身连接的边,即所述图模型中共有N*(N-1)/2+N条边。
假设目标点云空间有三个被占据的子空间v1、v2和v3,那么可以构建图模型如图3所示。
S130、根据所述图模型对所述子空间的第一属性系数进行变换,获取所述目标点云空间的第一属性系数和第二属性系数。
具体地,所述步骤S130包括子步骤:
S131、根据各个所述子空间之间的距离获取所述图模型中每条边的权重。
所述根据各个所述子空间之间的距离获取所述图模型中每条边的权重包括:
根据预设公式获取所述图模型中每条边的权重;
所述预设公式为:
其中,W(i,j)表示连接顶点i和顶点j的边的权重,ai、aj分别为顶点i和顶点j对应的点云子空间内包含的点数,xi、xj分别为顶点i和顶点j对应的点云子空间的位置,α为预设参数,α为可调值,控制权重随两点距离的增加而减小的速度。
从上面说明不难看出,本发明提供的点云属性编码方法中引入了点间距离的信息,从而可以更加准确地表达点间关系,实现更好的变换效果。
S132、根据所述图模型中每条边的权重获取所述图模型对应的归一化对称图拉普拉斯矩阵。
在获取到所述图模型中每条边的权重后,可以获取所述图模型对应的归一化对称图拉普拉斯矩阵,此为现有技术,在此不再赘述。
S133、将所述归一化对称图拉普拉斯矩阵对应的特征向量组作为变换矩阵,对所述子空间的第一属性系数进行图傅里叶变换,获取所述目标点云空间的第一属性系数和第二属性系数。
在获取到所述归一化对称图拉普拉斯矩阵后,求解所述归一化对称图拉普拉斯矩阵的特征向量组,将所述特征向量组作为变换矩阵,对所述子空间的第一属性系数进行图傅里叶变换,变换后可以得到所述第一属性系数和所述第二属性系数。变换可用公式表示为:C=F×A,其中,A为变换矩阵,C为所述目标点云空间的属性系数(包括第一属性系数和第二属性系数)组成的向量,F为所述子空间的第一属性系数组成的向量。
所述第一属性系数是所述特征向量组中特征值为0的向量对应的变换结果,所述第二属性系数为所述特征向量组中特征值不为0的特征向量对应的变换结果,具体地,拉普拉斯矩阵对应的特征向量组里,一定有一个特征值为0的特征向量,所述第一属性系数的数量为1个,所述第二属性系数的数量可以有多个。将所述第一属性系数记为DC系数(直流系数),将所述第二属性系数记为AC系数(交流系数)。
从上面的说明不难看出,在从所述多层级结构的第二级开始,对于所述多层级结构中的每一级进行逐层变换的过程中,每个点云空间的第一属性系数会在上一级中继续被变换,最终会得到一个第一属性系数和多个第二属性系数,即所述多层级结构中的所述目标点云的第一属性系数,以及各个所述目标点云空间的第二属性系数。
所述多层级结构中的第一级中的点云空间的第一属性系数为点云空间中的数据点的属性值。在逐层变换后,最终获取的第一属性系数和第二属性系数能够实现保留原有属性信息的基础上尽可能地去冗余,提高了存储和传输效率。
在对所述多层级结构逐层变换完成后,本实施例提供的基于连续子空间图变换的点云属性编码方法还包括步骤:
S140、根据所述多层级结构中的第M级中的点云空间的第一属性系数、以及所述多层级结构中各个所述目标点云空间的第二属性系数对所述目标点云进行属性编码。
所述多层级结构中的第M级中的点云空间即为所述目标点云,根据所述多层级结构中的第M级中的点云空间的第一属性系数、以及所述多层级结构中各个所述目标点云空间的第二属性系数对所述目标点云进行属性编码,可以是将第一属性系数和第二属性系数作为编码系数,进行量化后再熵编码,或者直接熵编码。
在一种可能的实现方式中,为了使得在对所述目标点云进行属性编码后得到的码流进行解码后得到的属性值能够尽可能地更接近所述目标点云的属性值,所述所述根据所述多层级结构中的第M级中的点云空间的第一属性系数、以及所述多层级结构中各个所述目标点云空间的第二属性系数对所述目标点云进行属性编码之后,包括:
对所述第M级中的点云空间的第一属性系数、以及所述多层级结构中各个所述目标点云的第二属性系数进行量化得到量化属性系数,对所述量化属性系数进行逆量化得到重建属性系数,根据所述图模型对所述重建属性系数进行逆变换,获取所述多层级结构中第1级的点云空间分别对应的重建属性值;
计算所述多层级结构中第1级的点云空间分别对应的所述原始属性值与所述重建属性值的差值作为属性残差值;
对所述属性残差值进行量化后再进行熵编码,得到属性残差码流;
或者,直接对所述属性残差值进行熵编码得到属性残差码流;
或者,对所述属性残差值进行量化得到属性残差量化系数并获取量化过程得到的余数,对所述属性量化残差系数进行熵编码得到属性残差码流,对所述余数进行熵编码,得到残差量化余数码流。
不难看出,所述第M级中的点云空间的第一属性系数、以及所述多层级结构中各个所述目标点云的第二属性系数分别对应有一个重建属性系数,所述根据所述图模型对所述重建属性系数进行逆变换,是指对于每个所述目标点云空间,都通过所述归一化对称图拉普拉斯矩阵对应的特征向量组的转置对所述目标点云空间的重建第一属性系数和第二属性系数对应的重建属性系数进行逆图傅里叶变换,得到所述目标点云空间的子空间的重建第一属性系数,这样逐层逆变换,最终得到所述多层级结构中第1级的点云空间的重建第一属性系数,将所述第1级的点云空间的重建第一属性系数作为对应的重建属性值。
所述属性残差码流、所述残差量化余数码流可以与对所述目标点云进行属性编码之后得到的重建点云的码流一同被存储,在对所述重建点云进行解码时被应用。具体在实施例二中被说明。
本实施例提供的基于连续子空间图变换的点云属性编码方法的效果图可以如图4所示,图4中是采用本发明提供的基于连续子空间图变换的点云属性编码方法和RAHT方法对同一点云进行颜色属性编码的效果图,图4中的横坐标标识压缩率,纵坐标表示颜色的亮度成分的失真程度,四条曲线中最下面一条为RAHT方法对应的曲线,上面三条至上而下分别为本发明提供的方法中多层级结构的中子空间大小为8x8x8个体素、4x4x4个体素、2x2x2个体素对应的曲线,可以看到,本发明提供的编码方法相对于RAHT方法,具有更小的压缩率和更低失真程度,具有更好的压缩效果。
实施例二
基于实施例一中提供的基于连续子空间图变换的点云属性编码方法,本发明还提供了一种基于连续子空间图变换的点云属性解码方法,如图5所示,本发明提供的基于连续子空间图变换的点云属性解码方法的一个实施例中,包括步骤:
S210、根据目标点云内点的位置将所述目标点云划分为多层级结构。S220、对于所述多层级结构中的每一个目标点云空间,获取所述目标点云空间的子空间,根据所述子空间构建图模型。
所述多层级结构包括M个层级,每级中包括至少两个点云空间,每一级的点云空间是对上一级的点云空间进行划分得到的,其中,M为大于1的整数。步骤S210以及步骤S220的具体说明可以参照实施例一中步骤S110和S120的说明。值得说明的是,在本实施例提供的基于连续子空间图变换的点云属性解码方法中,将所述目标点云进行划分得到的多层级结构应与对所述目标点云进行属性编码时对所述目标点云进行划分得到的多层级结构一致,二者中的点云空间一一对应。
S230、对所述目标点云的编码码流进行解码,获取所述多层级结构的第M级中的点云空间的第一属性系数和各个所述目标点云空间分别对应的第二属性系数,并根据所述图模型对所述目标点云空间的第一属性系数和第二属性系数进行变换,获取各个所述子空间的第一属性系数。
在本实施例中,所述目标点云的码流是通过实施例一中提供的基于连续子空间图变换的点云属性编码方法编码得到的,对所述目标点云的码流进行解码,可以获取到所述多层级结构的第M级中的点云空间的第一属性系数和各个所述目标点云空间分别对应的第二属性系数。
对所述目标点云的码流进行解码的方法根据对所述目标点云进行编码生成码流的方式确定,当对所述目标点云进行属性编码时,是将属性系数进行量化后再熵编码时,那么解码时是先对点云码流进行熵解码再进行逆量化,当对所述目标点云进行属性编码时,是直接将属性系数进行熵编码,那么解码使对点云码流进行熵解码即可。
所述根据所述子空间构建图模型可以参照实施例一中对步骤S120的说明。所述根据所述图模型对所述目标点云空间的第一属性系数和第二属性系数进行变换,获取各个所述子空间的第一属性系数,包括:
S231、根据各个所述子空间之间的距离获取所述图模型中每条边的权重;
S232、根据所述图模型中每条边的权重获取所述图模型对应的归一化对称图拉普拉斯矩阵;
S233、将所述归一化对称图拉普拉斯矩阵对应的特征向量组的转置作为变换矩阵,对所述目标点云空间的第一属性系数和第二属性系数进行逆图傅里叶变换,获取各个所述子空间的第一属性系数。
所述步骤S231、步骤S232的具体说明可以参照实施例一中对步骤S131、步骤S132的说明。
在获取到所述归一化对称图拉普拉斯矩阵后,求解所述归一化图拉普拉斯矩阵的特征向量组的转置作为变换矩阵,将所述目标点云空间的第一属性系数和第二属性系数作为输入进行逆图傅里叶变换,得到各个所述子空间的第一属性系数。
具体地,逆变换可用公式表示为:F=C×AT,其中,AT为变换矩阵,是归一化对称图拉普拉斯矩阵的转置,C为所述目标点云空间的属性系数(包括第一属性系数和第二属性系数)组成的向量,F为所述子空间的第一属性系数组成的向量。
请再次参阅图5,本实施例提供的基于连续子空间图变换的编码方法中还包括步骤:
S240、根据所述多层级结构中的第一级中各个点云空间的第一属性系数获取所述目标点云的属性值。
对所述目标点云进行属性编码后得到的码流进行解码,目的是为了获取所述目标点云的属性值。从上面的说明不难看出,在经过对所述多层级结构中从所述目标点云开始到第2级的目标点云空间逐级进行变换后,最终可以得到所述多层级结构中的第一级中的点云空间的第一属性系数,在一种可能的实现方式中,可以直接将所述多层级结构中的第一级中的点云空间的属性系数作为点云空间的属性值,这样就获取到了所述目标点云中各个数据点的属性值,完成对所述目标点云对应的属性编码码流进行解码的过程。
而在本实施例中,为了使得解码后的到所述目标点云的属性值能够尽可能地接近真实值,降低编解码带来的失真程度,还通过对所述目标点云进行属性编码时生成的属性残差值等参数来获取所述目标点云的属性值,具体地,根据所述多层级结构中的第一级中各个点云空间的第一属性系数获取所述目标点云的属性值包括:
获取所述目标点云对应的属性残差码流,对所述属性残差码流进行熵解码后再进行逆量化操作,得到属性残差值,获取所述属性残差值和所述多层级结构中的第一级中各个点云空间的第一属性系数的和作为所述目标点云的属性值;
或者,获取所述目标点云对应的属性残差码流,直接对所述属性残差码流进行熵解码,得到属性残差值,获取所述属性残差值和所述多层级结构中的第一级中各个点云空间的第一属性系数的和作为所述目标点云的属性值;
或者,获取所述目标点云对应的属性残差码流以及残差量化余数码流,对所述属性残差码流和所述残差量化余数码流进行熵解码得到属性残差量化系数和量化过程得到的余数,对所述属性参数量化系数进行逆量化后与所述余数相加,得到属性残差值,获取所述属性残差值和所述多层级结构中的第一级中各个点云空间的第一属性系数的和作为所述目标点云的属性值。
从实施例一的说明中不难看出,所述属性残差值包括多个值,分别与所述多层级结构中第1级的点云空间一一对应,在获取所述属性残差值与所述多层级结构中的第一级中各个点云空间的第一属性系数的和时,是根据点云空间将属性残差值和第一属性系数对应相加,即,对于某一个点云空间来说,将该点云空间对应的属性残差值和对应的第一属性系数进行相加,得到该点云空间对应的属性值。综上所述,本发明提供了一种基于连续子空间图变换的点云属性编码方法及解码方法,通过将目标点云划分为多层级结构,对多层级结构的第二级以及第二级以上的目标点云空间,都由该目标点云空间在下一级的子空间的第一属性系数进行图变换得到该节点的第一属性系数和第二属性系数,目标点云空间的第一属性系数在上一级继续被变换,这样进行连续变换后,最终得到所述目标点云的第一属性系数和各个目标点云空间的第二属性系数,再根据所述目标点云的第一属性系数和各个目标点云空间的第二属性系数进行点云属性编码,在编码的过程中,不涉及寻找变换点的过程,避免了寻找方向的顺序影响点云属性编码效果,保证了点云属性编码质量。
应该理解的是,虽然本发明说明书附图中给出的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,流程图中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本发明所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限,RAM以多种形式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synch l i nk)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
实施例三
基于上述实施例,本发明还相应提供了一种终端,如图6所示,所述终端包括处理器10以及存储器20。可以理解的是,图6仅示出了终端的部分组件,但是应理解的是,并不要求实施所有示出的组件,可以替代的实施更多或者更少的组件。
所述存储器20在一些实施例中可以是所述终端的内部存储单元,例如终端的硬盘或内存。所述存储器20在另一些实施例中也可以是所述终端的外部存储设备,例如所述终端上配备的插接式硬盘,智能存储卡(Smart Med ia Card,SMC),安全数字(Secure Digital,SD)卡,闪存卡(F l ash Card)等。进一步地,所述存储器20还可以既包括所述终端的内部存储单元也包括外部存储设备。所述存储器20用于存储安装于所述终端的应用软件及各类数据。所述存储器20还可以用于暂时地存储已经输出或者将要输出的数据。在一些实施例中,存储器20上存储有计算机程序30,该计算机程序30可被处理器10所执行,从而实现本申请中基于连续子空间图变换的点云属性编码方法和/或基于连续子空间图变换的点云属性解码方法。
所述处理器10在一些实施例中可以是一中央处理器(Centra l Process i ngUnit,CPU),微处理器或其他芯片,用于运行所述存储器20中存储的程序代码或处理数据,例如执行所述基于连续子空间图变换的点云属性编码方法和/或基于连续子空间图变换的点云属性解码方法等。
实施例四
本发明还提供一种存储介质,其中,存储有一个或者多个程序,所述一个或者多个程序可被一个或者多个处理器执行,以实现如上所述的基于连续子空间图变换的点云属性编码方法的步骤。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (9)
1.一种基于连续子空间图变换的点云属性编码方法,其特征在于,所述方法包括:
根据目标点云内点的位置将所述目标点云划分为多层级结构,所述多层级结构包括M个层级,每级中包括至少两个点云空间,每一级的点云空间是对上一级的点云空间进行划分得到的,其中,M为大于1的整数;
对于所述多层级结构中的每一个目标点云空间,获取所述目标点云空间的子空间,根据所述子空间构建图模型,其中,所述目标点云空间为所述多层级结构中第一级之外的层级中的点云空间;
根据所述图模型对所述子空间的第一属性系数进行变换,获取所述目标点云空间的第一属性系数和第二属性系数;
所述根据所述图模型对所述子空间的第一属性系数进行变换,获取所述目标点云空间的第一属性系数和第二属性系数,包括:
根据各个所述子空间之间的距离获取所述图模型中每条边的权重;
根据所述图模型中每条边的权重获取所述图模型对应的归一化对称图拉普拉斯矩阵;
将所述归一化对称图拉普拉斯矩阵对应的特征向量组作为变换矩阵,对所述子空间的第一属性系数进行图傅里叶变换,获取所述目标点云空间的第一属性系数和第二属性系数;
其中,所述目标点云空间的第一属性系数为所述特征向量组中特征值为0的特征向量对应的变换结果,所述目标点云空间的第二属性系数为所述特征向量组中特征值不为0的特征向量对应的变换结果;
根据所述多层级结构中的第M级中的点云空间的第一属性系数、以及所述多层级结构中各个所述目标点云空间的第二属性系数对所述目标点云进行属性编码;
其中,所述多层级结构中的第一级中的点云空间的第一属性系数为点云空间的原始属性值。
2.根据权利要求1所述的基于连续子空间图变换的点云属性编码方法,其特征在于,所述根据目标点云内点的位置将所述目标点云划分为多层级结构,包括:
对点云空间进行划分,得到K个子空间,每个子空间的空间大小相等或包含的数据点数量相等,每个子空间内包括至少一个所述目标点云的数据点,其中,K为大于1的整数;
当点云空间内包含的数据点的数量为1时,不再继续划分。
3.根据权利要求1所述的基于连续子空间图变换的点云属性编码方法,其特征在于,所述根据所述子空间构建图模型包括:
将所述子空间作为所述图模型中的顶点,将所述子空间两两连接的连接线作为所述图模型中的边。
4.根据权利要求1所述的基于连续子空间图变换的点云属性编码方法,其特征在于,所述根据各个所述子空间之间的距离获取所述图模型中每条边的权重,包括:
根据预设公式获取所述图模型中每条边的权重;
所述预设公式为:
其中,W(i,j)表示连接顶点i和顶点j的边的权重,ai、aj分别为顶点i和顶点j对应的点云子空间内包含的点数,Xi、Xj分别为顶点i和顶点j对应的点云子空间的位置,a为预设参数。
5.根据权利要求1所述的基于连续子空间图变换的点云属性编码方法,其特征在于,所述根据所述多层级结构中的第M级中的点云空间的第一属性系数、以及所述多层级结构中各个所述目标点云空间的第二属性系数对所述目标点云进行属性编码之后,包括:
对所述第M级中的点云空间的第一属性系数、以及所述多层级结构中各个所述目标点云的第二属性系数进行量化得到量化属性系数,对所述量化属性系数进行逆量化得到重建属性系数,根据所述图模型对所述重建属性系数进行逆变换,获取所述多层级结构中第1级的点云空间分别对应的重建属性值;计算所述多层级结构中第1级的点云空间分别对应的所述原始属性值与所述重建属性值的差值作为属性残差值;
对所述属性残差值进行量化后再进行熵编码,得到属性残差码流;
或者,直接对所述属性残差值进行熵编码得到属性残差码流;
或者,对所述属性残差值进行量化得到属性残差量化系数并获取量化过程得到的余数,对所述属性量化残差系数进行熵编码得到属性残差码流,对所述余数进行熵编码,得到残差量化余数码流。
6.一种基于连续子空间图变换的点云属性解码方法,其特征在于,所述方法包括:
根据目标点云内点的位置将所述目标点云划分为多层级结构,所述多层级结构包括M个层级,每级中包括至少两个点云空间,每一级的点云空间是对上一级的点云空间进行划分得到的,其中,M为大于1的整数;
对于所述多层级结构中的每一个目标点云空间,获取所述目标点云空间的子空间,根据所述子空间构建图模型,其中,所述目标点云空间为所述多层级结构中第一级之外的层级中的点云空间;
对所述目标点云的编码码流进行解码,获取所述多层级结构的第M级中的点云空间的第一属性系数和各个所述目标点云空间分别对应的第二属性系数,并根据所述图模型对所述目标点云空间的第一属性系数和第二属性系数进行变换,获取各个所述子空间的第一属性系数;
根据所述多层级结构中的第一级中各个点云空间的第一属性系数获取所述目标点云的属性值;
所述根据所述图模型对所述目标点云空间的第一属性系数和第二属性系数进行变换,获取各个所述子空间的第一属性系数,包括:
根据各个所述子空间之间的距离获取所述图模型中每条边的权重;
根据所述图模型中每条边的权重获取所述图模型对应的归一化对称图拉普拉斯矩阵;
将所述归一化对称图拉普拉斯矩阵对应的特征向量组的转置作为变换矩阵,对所述目标点云空间的第一属性系数和第二属性系数进行逆图傅里叶变换,获取各个所述子空间的第一属性系数。
7.根据权利要求6所述的基于连续子空间图变换的点云属性解码方法,其特征在于,所述根据所述多层级结构中的第一级中各个点云空间的第一属性系数获取所述目标点云的属性值,包括:
获取所述目标点云对应的属性残差码流,对所述属性残差码流进行熵解码后再进行逆量化操作,得到属性残差值,获取所述属性残差值和所述多层级结构中的第一级中各个点云空间的第一属性系数的和作为所述目标点云的属性值;
或者,获取所述目标点云对应的属性残差码流,直接对所述属性残差码流进行熵解码,得到属性残差值,获取所述属性残差值和所述多层级结构中的第一级中各个点云空间的第一属性系数的和作为所述目标点云的属性值;
或者,获取所述目标点云对应的属性残差码流以及残差量化余数码流,对所述属性残差码流和所述残差量化余数码流进行熵解码得到属性残差量化系数和量化过程得到的余数,对所述属性参数量化系数进行逆量化后与所述余数相加,得到属性残差值,获取所述属性残差值和所述多层级结构中的第一级中各个点云空间的第一属性系数的和作为所述目标点云的属性值。
8.一种终端,其特征在于,所述终端包括:处理器、与处理器通信连接的存储介质,所述存储介质适于存储多条指令,所述处理器适于调用所述存储介质中的指令,以执行实现上述权利要求1-5任一项所述的基于连续子空间图变换的点云属性编码方法的步骤和/或上述权利要求6-7任一项所述的基于连续子空间图变换的点云属性解码方法的步骤。
9.一种存储介质,其特征在于,所述存储介质存储有一个或者多个程序,所述一个或者多个程序可被一个或者多个处理器执行,以实现如权利要求1-5任一项所述的基于连续子空间图变换的点云属性编码方法的步骤和/或上述权利要求6-7任一项所述的基于连续子空间图变换的点云属性解码方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010966930.8A CN114189692B (zh) | 2020-09-15 | 2020-09-15 | 基于连续子空间图变换的点云属性编码方法及解码方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010966930.8A CN114189692B (zh) | 2020-09-15 | 2020-09-15 | 基于连续子空间图变换的点云属性编码方法及解码方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114189692A CN114189692A (zh) | 2022-03-15 |
CN114189692B true CN114189692B (zh) | 2024-05-14 |
Family
ID=80601196
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010966930.8A Active CN114189692B (zh) | 2020-09-15 | 2020-09-15 | 基于连续子空间图变换的点云属性编码方法及解码方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114189692B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107403456A (zh) * | 2017-07-28 | 2017-11-28 | 北京大学深圳研究生院 | 一种基于kd树和优化图变换的点云属性压缩方法 |
CN110572655A (zh) * | 2019-09-30 | 2019-12-13 | 北京大学深圳研究生院 | 一种基于邻居权重的参数选取和传递的点云属性编码和解码的方法及设备 |
CN111145090A (zh) * | 2019-11-29 | 2020-05-12 | 鹏城实验室 | 一种点云属性编码方法、解码方法、编码设备及解码设备 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10223810B2 (en) * | 2016-05-28 | 2019-03-05 | Microsoft Technology Licensing, Llc | Region-adaptive hierarchical transform and entropy coding for point cloud compression, and corresponding decompression |
-
2020
- 2020-09-15 CN CN202010966930.8A patent/CN114189692B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107403456A (zh) * | 2017-07-28 | 2017-11-28 | 北京大学深圳研究生院 | 一种基于kd树和优化图变换的点云属性压缩方法 |
CN110572655A (zh) * | 2019-09-30 | 2019-12-13 | 北京大学深圳研究生院 | 一种基于邻居权重的参数选取和传递的点云属性编码和解码的方法及设备 |
CN111145090A (zh) * | 2019-11-29 | 2020-05-12 | 鹏城实验室 | 一种点云属性编码方法、解码方法、编码设备及解码设备 |
Also Published As
Publication number | Publication date |
---|---|
CN114189692A (zh) | 2022-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Huang et al. | 3d point cloud geometry compression on deep learning | |
JP6676193B2 (ja) | シーンを表す点群を符号化する方法、符号化器システム、及びプログラムを記憶した非一時的コンピューター可読記録媒体 | |
US10964068B2 (en) | Methods and devices for predictive point cloud attribute coding | |
US11836954B2 (en) | 3D point cloud compression system based on multi-scale structured dictionary learning | |
Shen et al. | Dct-mask: Discrete cosine transform mask representation for instance segmentation | |
JP7476631B2 (ja) | 画像コーディング方法及び装置並びに画像デコーディング方法及び装置 | |
US10599935B2 (en) | Processing artificial neural network weights | |
US9349072B2 (en) | Local feature based image compression | |
CN113766229B (zh) | 一种编码方法、解码方法、装置、设备及可读存储介质 | |
CN111641832A (zh) | 编码方法、解码方法、装置、电子设备及存储介质 | |
US11570477B2 (en) | Data preprocessing and data augmentation in frequency domain | |
CN112119408A (zh) | 获取画质增强网络的方法、图像画质增强方法、装置、可移动平台、相机及存储介质 | |
WO2023277974A1 (en) | Systems and methods of hierarchical implicit representation in octree for 3d modeling | |
CN114598891A (zh) | 点云数据编码方法、解码方法、点云数据处理方法及装置 | |
CN113438481B (zh) | 训练方法、图像编码方法、图像解码方法及装置 | |
CN115393452A (zh) | 一种基于非对称自编码器结构的点云几何压缩方法 | |
CN114189692B (zh) | 基于连续子空间图变换的点云属性编码方法及解码方法 | |
JP6748022B2 (ja) | 分割形状決定装置、学習装置、分割形状決定方法及び分割形状決定プログラム | |
Zhang et al. | Graph-based transform for 2D piecewise smooth signals with random discontinuity locations | |
TW202406344A (zh) | 一種點雲幾何資料增強、編解碼方法、裝置、碼流、編解碼器、系統和儲存媒介 | |
CN113411615B (zh) | 一种面向虚拟现实的纬度自适应的全景图像编码方法 | |
WO2023205969A1 (zh) | 点云几何信息的压缩、解压缩及点云视频编解码方法、装置 | |
CN113554719B (zh) | 一种图像编码方法、解码方法、存储介质及终端设备 | |
CN116648915A (zh) | 点云编解码方法、编码器、解码器以及计算机存储介质 | |
WO2023248486A1 (ja) | 情報処理装置および方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |