CN108022221A - 一种各项异性扩散滤波的三维散乱点云平滑去噪方法 - Google Patents
一种各项异性扩散滤波的三维散乱点云平滑去噪方法 Download PDFInfo
- Publication number
- CN108022221A CN108022221A CN201711278418.9A CN201711278418A CN108022221A CN 108022221 A CN108022221 A CN 108022221A CN 201711278418 A CN201711278418 A CN 201711278418A CN 108022221 A CN108022221 A CN 108022221A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mtd
- msubsup
- mtr
- 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.)
- Pending
Links
- 238000009792 diffusion process Methods 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000001914 filtration Methods 0.000 title claims abstract description 17
- 239000006185 dispersion Substances 0.000 title claims abstract description 11
- 238000009499 grossing Methods 0.000 title claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims abstract description 32
- 238000013461 design Methods 0.000 claims description 7
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000005070 sampling Methods 0.000 abstract description 3
- 230000007812 deficiency Effects 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000010146 3D printing Methods 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000004806 packaging method and process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10028—Range image; Depth image; 3D point clouds
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种各项异性扩散滤波的三维散乱点云平滑去噪方法,通过对采样点与其有效邻域点进行张量投票得到张量矩阵结构张量矩阵,求解特征值和特征向量并根据结构张量矩阵的特征值和特征向量来分析采样点的局部特性,根据采样点的不同几何特征信息来设计扩散张量矩阵的特征值,并根据不同的几何特征信息来设计扩散速率,使得在不同主特征方向上的扩散速率不同,重构出修正后的扩散张量矩阵,最后将重构好的扩散张量代入三维扩散异性滤波方程中进行微分求解,经过一定次数的迭代得到滤波因子,对噪声进行平滑处理。本方法对散乱点云去除噪声的同时,可以有效的保持原始模型的特征信息,避免了过光顺和局部失真。
Description
技术领域
本发明涉及计算机视觉技术与三维重建的逆向工程领域,具体是一种各项异性扩散滤波的三维散乱点云平滑去噪方法。
背景技术
物体三维点云数据的获取是逆向工程的重要环节。但由于测量设备本身的缺陷以及测量环境的影响,导致测量得到的三维点云数据不可避免的存在噪声。因此,三维点云数据的平滑去噪处理已经成为逆向工程研究领域的热点问题之一。
随着计算机与三维数字化的不断发展,三维点模型也被广泛应用于3D打印、虚拟现实和计算机建模等领域。然而通过三维扫描仪得到的点云数据中含有噪声,传统的去噪方法在去除噪声时,易导致模型特征细节的丢失,所以一种点云保特征的去噪方法显得尤为重要。
发明内容
针对现有技术的不足,本发明拟解决的技术问题是,提供一种各向异性扩散滤波的三维点云去噪方法。
本发明解决所述技术问题的技术方案是,提供一种各项异性扩散滤波的三维散乱点云平滑去噪方法,其特征在于该方法包括以下步骤:
步骤1,输入三维散乱点云数据P={p1,...,pn},pi为其中的一个采样点,确定其邻域点个数k;
步骤2,运用K-D树算法搜索邻域点:设邻域点集Ni={pi1,pi2,pi3,...,pik},i=1,2,...,n,k≤n;
步骤3,对采样点和邻域点进行张量投票,得到结构张量矩阵T,采样点pi和其邻域点集Ni进行张量投票得到的结构张量矩阵T如式(1)所示:
其中T3为3×3的位矩阵,为采样点pi与其邻域点pik连线的方向向量,σ为自定义参数;结构张量矩阵T的特征值和特征向量的表示形式如式(2)所示:
式中,λ1,λ2,λ3分别代表了采样点的结构张量矩阵T的3个不同的特征值,令λ1≥λ2≥λ3≥0,e1,e2,e3分别代表了不同特征值对应的特征向量,同时也代表了采样点的三个主特征方向;
步骤4,构造扩散张量矩阵D:扩散张量矩阵D为3×3的对称矩阵,如式(3)所示:
扩散张量矩阵D的特征值和特征向量的表示形式如式(4)所示:
扩散张量矩阵D的特征向量和结构张量矩阵T的特征向量相同都为e1、e2和e3,所以只需要根据点云的方向结构特征设计相应的特征值u1,u2,u3,并且u1,u2,u3的值都属于[0,1];
步骤5,为了区分采样点的结构张量矩阵T的特征值的大小定义了采样点的相干系数,如式(5)所示:
步骤6,根据相干系数的大小对扩散张量矩阵D的特征值进行设计,设计方法如式(6)所示,其中α是用于控制扩散速率的因子;
u1=α
步骤7,利用所设计的扩散张量矩阵D的特征值μ1、μ2和μ3与结构张量矩阵T的特征向量e1、e2和e3来重构扩散张量矩阵D,得到重构后的扩散张量矩阵D如式(7)所示:
步骤8,将重构后的扩散张量矩阵D带入式(8),求解过程如下:
式中,j1、j2和j3表示为通量,具体是
在式(8)中div表示散度,则式(8)的右侧可以表示为式(9):
式(8)的左侧利用差分代替微分来计算,展开后如式(10)所示:
其中,和分别为I(x,y,z,τ)在时刻(k+1)τ和时刻kτ的位置坐标信息;τ为时间步长因子;式(8)的离散可以表示为式(11):
步骤9,对所有采样点重复步骤2-步骤8,得到去噪后的点云模型。
与现有技术相比,本发明有益效果在于:
利用K-D树搜索采样点的邻域点,并采用张量投票的方法获取三维点云的结构张量矩阵,求解特征值和特征向量并根据结构张量矩阵的特征值和特征向量来分析采样点的局部特性,根据采样点的不同几何特征信息来设计扩散张量矩阵的特征值,并根据不同的几何特征信息来设计扩散速率,使得在不同主特征方向上的扩散速率不同,重构出修正后的扩散张量矩阵,保证在不同特征方向的扩散速率能自适应调整。最后将重构好的扩散张量代入三维扩散异性滤波方程中进行微分求解,经过一定次数的迭代得到滤波因子,对噪声进行平滑处理。
利用重构的扩散张量进行滤波处理时,滤波因子能根据采样点的几何特征自适应调整,避免了在模型尖锐导致的过光顺或平坦区域的光顺不足。
对散乱点云去除噪声的同时,可以有效的保持原始模型的特征信息,保证点云模型的高频特征不失真,避免了过光顺和局部失真。
附图说明
图1为本发明各项异性扩散滤波的三维散乱点云平滑去噪方法一种实施例的采用该方法的处理对象-Fandisk原始模型;
图2为本发明各项异性扩散滤波的三维散乱点云平滑去噪方法一种实施例的采用该方法的处理对象-Fandisk加噪模型;
图3为本发明各项异性扩散滤波的三维散乱点云平滑去噪方法一种实施例的采用该方法的处理对象-Fandisk去噪后模型;
具体实施方式
下面给出本发明的具体实施例。具体实施例仅用于进一步详细说明本发明,不限制本申请权利要求的保护范围。
本发明提供了一种各项异性扩散滤波的三维散乱点云平滑去噪方法(简称方法),其特征在于该方法包括以下步骤:
步骤1,输入三维散乱点云数据P={p1,...,pn},pi为其中的一个采样点,确定其邻域点个数k;
步骤2,运用K-D树算法搜索邻域点:设邻域点集Ni={pi1,pi2,pi3,...,pik},i=1,2,...,n,k≤n。
步骤3,对采样点和邻域点进行张量投票,得到结构张量矩阵T,采样点pi和其邻域点集Ni进行张量投票得到的结构张量矩阵T如式(1)所示:
其中T3为3×3的位矩阵,为采样点pi与其邻域点pik连线的方向向量,σ为自定义参数;结构张量矩阵T的特征值和特征向量的表示形式如式(2)所示:
式中,λ1,λ2,λ3分别代表了采样点的结构张量矩阵T的3个不同的特征值,令λ1≥λ2≥λ3≥0,e1,e2,e3分别代表了不同特征值对应的特征向量,同时也代表了采样点的三个主特征方向;
步骤4,构造了一个对点云局部特征具有鲁棒特性的扩散张量,扩散张量的构建不仅要考虑它的方向信息而且还要考虑它的梯度模值,本发明所构造的扩散张量矩阵D为3×3的对称矩阵,如式(3)所示:
扩散张量矩阵D的特征值和特征向量的表示形式如式(4)所示:
扩散张量矩阵D的特征向量和结构张量矩阵T的特征向量相同都为e1、e2和e3,所以只需要根据点云的方向结构特征设计相应的特征值u1,u2,u3,并且u1,u2,u3的值都属于[0,1];
步骤5,为了区分采样点的结构张量矩阵T的特征值的大小定义了采样点的相干系数,如式(5)所示:
步骤6,根据相干系数的大小对扩散张量矩阵D的特征值进行设计,设计方法如式(6)所示,其中α是用于控制扩散速率的因子,本实施例中取α=0.005;
u1=α
步骤7,利用所设计的扩散张量矩阵D的特征值μ1、μ2和μ3与结构张量矩阵T的特征向量e1、e2和e3来重构扩散张量矩阵D,得到重构后的扩散张量矩阵D如式(7)所示:
步骤8,通过扩散张量矩阵D的重构可以得到扩散张量矩阵D中的每个元素的具体值:将重构后的扩散张量矩阵D带入式(8),式(8)为三维扩散异性滤波方程,求解过程如下:
式中,j1、j2和j3表示为通量,具体是
在式(8)中div表示散度,则式(8)的右侧可以表示为式(9):
式(8)的左侧利用差分代替微分来计算,展开后如式(10)所示:
其中,和分别为I(x,y,z,τ)在时刻(k+1)τ和时刻kτ的位置坐标信息;τ为时间步长因子,取值为1;式(8)的离散可以表示为式(11):
步骤9,对所有采样点重复步骤2-步骤8,得到去噪后的点云模型。
实施例
本文利用Fandisk点云模型进行了实验仿真研究,图1为Fandisk的原始点云封装模型,图2为加入30dB的高斯白噪声后的封装模型,图3为利用本文算法迭代2次的去噪结果,从去噪结果看,去噪后模型未发生特征细节的丢失,模型的体积为发生变形。
本发明未述及之处适用于现有技术。
Claims (1)
1.一种各项异性扩散滤波的三维散乱点云平滑去噪方法,其特征在于该方法包括以下步骤:
步骤1,输入三维散乱点云数据P={p1,...,pn},pi为其中的一个采样点,确定其邻域点个数k;
步骤2,运用K-D树算法搜索邻域点:设邻域点集Ni={pi1,pi2,pi3,...,pik},i=1,2,...,n,k≤n;
步骤3,对采样点和邻域点进行张量投票,得到结构张量矩阵T,采样点pi和其邻域点集Ni进行张量投票得到的结构张量矩阵T如式(1)所示:
<mrow>
<mi>T</mi>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>&Element;</mo>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
</mrow>
</munder>
<msub>
<mi>&mu;</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mover>
<mi>&upsi;</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>k</mi>
</msub>
<msubsup>
<mover>
<mi>&upsi;</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>k</mi>
<mi>T</mi>
</msubsup>
</mrow>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mover>
<mi>&upsi;</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>k</mi>
</msub>
<msubsup>
<mover>
<mi>&upsi;</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>k</mi>
<mi>T</mi>
</msubsup>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中T3为3×3的位矩阵,为采样点pi与其邻域点pik连线的方向向量,σ为自定义参数;结构张量矩阵T的特征值和特征向量的表示形式如式(2)所示:
<mrow>
<mi>T</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>3</mn>
</munderover>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<msub>
<mi>e</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>e</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>e</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>e</mi>
<mn>2</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>e</mi>
<mn>3</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mn>3</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>e</mi>
<mn>1</mn>
<mi>T</mi>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>e</mi>
<mn>2</mn>
<mi>T</mi>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>e</mi>
<mn>3</mn>
<mi>T</mi>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,λ1,λ2,λ3分别代表了采样点的结构张量矩阵T的3个不同的特征值,令λ1≥λ2≥λ3≥0,e1,e2,e3分别代表了不同特征值对应的特征向量,同时也代表了采样点的三个主特征方向;
步骤4,构造扩散张量矩阵D:扩散张量矩阵D为3×3的对称矩阵,如式(3)所示:
<mrow>
<mi>D</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>D</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>D</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>D</mi>
<mn>13</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>D</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>D</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>D</mi>
<mn>23</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>D</mi>
<mn>31</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>D</mi>
<mn>32</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>D</mi>
<mn>33</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
扩散张量矩阵D的特征值和特征向量的表示形式如式(4)所示:
<mrow>
<mi>D</mi>
<mo>=</mo>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
<msub>
<mi>e</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>1</mn>
<mi>T</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
<msub>
<mi>e</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>2</mn>
<mi>T</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>u</mi>
<mn>3</mn>
</msub>
<msub>
<mi>e</mi>
<mn>3</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>3</mn>
<mi>T</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
扩散张量矩阵D的特征向量和结构张量矩阵T的特征向量相同都为e1、e2和e3,所以只需要根据点云的方向结构特征设计相应的特征值u1,u2,u3,并且u1,u2,u3的值都属于[0,1];
步骤5,为了区分采样点的结构张量矩阵T的特征值的大小定义了采样点的相干系数,如式(5)所示:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>k</mi>
<mn>12</mn>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>k</mi>
<mn>13</mn>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&lambda;</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>k</mi>
<mn>23</mn>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&lambda;</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤6,根据相干系数的大小对扩散张量矩阵D的特征值进行设计,设计方法如式(6)所示,其中α是用于控制扩散速率的因子;
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>&alpha;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mi>&alpha;</mi>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
<mi> </mi>
<msub>
<mi>k</mi>
<mn>12</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&alpha;</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>k</mi>
<mn>12</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>e</mi>
<mi>l</mi>
<mi>s</mi>
<mi>e</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mi>&alpha;</mi>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
<mi> </mi>
<msub>
<mi>k</mi>
<mn>23</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&alpha;</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>k</mi>
<mn>23</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>e</mi>
<mi>l</mi>
<mi>s</mi>
<mi>e</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤7,利用所设计的扩散张量矩阵D的特征值μ1、μ2和μ3与结构张量矩阵T的特征向量e1、e2和e3来重构扩散张量矩阵D,得到重构后的扩散张量矩阵D如式(7)所示:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mn>11</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&mu;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>11</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>21</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mn>3</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>31</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mn>22</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&mu;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>12</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>22</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mn>3</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>32</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mn>33</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&mu;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>13</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>23</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mn>3</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>33</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mn>12</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&mu;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>11</mn>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>e</mi>
<mn>12</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>21</mn>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>e</mi>
<mn>22</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mn>3</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>31</mn>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>e</mi>
<mn>32</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mn>11</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&mu;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>11</mn>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>e</mi>
<mn>13</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>21</mn>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>e</mi>
<mn>23</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mn>3</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>31</mn>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>e</mi>
<mn>33</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mn>11</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&mu;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>12</mn>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>e</mi>
<mn>13</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>22</mn>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>e</mi>
<mn>23</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&mu;</mi>
<mn>3</mn>
</msub>
<msubsup>
<mi>e</mi>
<mn>32</mn>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>e</mi>
<mn>33</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤8,将重构后的扩散张量矩阵D带入式(8),求解过程如下:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>I</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>d</mi>
<mi>i</mi>
<mi>v</mi>
<mo>&lsqb;</mo>
<mi>D</mi>
<mo>&dtri;</mo>
<mi>I</mi>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>,</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>I</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,j1、j2和j3表示为通量,具体是
在式(8)中div表示散度,则式(8)的右侧可以表示为式(9):
<mrow>
<mi>d</mi>
<mi>i</mi>
<mi>v</mi>
<mo>&lsqb;</mo>
<mi>D</mi>
<mo>&dtri;</mo>
<mi>I</mi>
<mo>&rsqb;</mo>
<mo>=</mo>
<msub>
<mo>&part;</mo>
<mi>x</mi>
</msub>
<msub>
<mi>j</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mo>&part;</mo>
<mi>y</mi>
</msub>
<msub>
<mi>j</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mo>&part;</mo>
<mi>z</mi>
</msub>
<msub>
<mi>j</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
式(8)的左侧利用差分代替微分来计算,展开后如式(10)所示:
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>I</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>I</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>I</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
</mrow>
<mi>k</mi>
</msubsup>
</mrow>
<mi>&tau;</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,和分别为I(x,y,z,τ)在时刻(k+1)τ和时刻kτ的位置坐标信息;τ为时间步长因子;式(8)的离散可以表示为式(11):
<mrow>
<msubsup>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mi>k</mi>
</msubsup>
<mo>+</mo>
<mi>&tau;</mi>
<mi>d</mi>
<mi>i</mi>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>&dtri;</mo>
<mi>I</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤9,对所有采样点重复步骤2-步骤8,得到去噪后的点云模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711278418.9A CN108022221A (zh) | 2017-12-06 | 2017-12-06 | 一种各项异性扩散滤波的三维散乱点云平滑去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711278418.9A CN108022221A (zh) | 2017-12-06 | 2017-12-06 | 一种各项异性扩散滤波的三维散乱点云平滑去噪方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108022221A true CN108022221A (zh) | 2018-05-11 |
Family
ID=62078625
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711278418.9A Pending CN108022221A (zh) | 2017-12-06 | 2017-12-06 | 一种各项异性扩散滤波的三维散乱点云平滑去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108022221A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109508489A (zh) * | 2018-11-07 | 2019-03-22 | 山东大学 | 一种各向异性多孔结构的建模方法及*** |
CN112381823A (zh) * | 2021-01-13 | 2021-02-19 | 北京智拓视界科技有限责任公司 | 用于图像的几何特征的提取方法和相关产品 |
CN112837331A (zh) * | 2021-03-08 | 2021-05-25 | 电子科技大学 | 一种基于自适应形态重建模糊三维sar图像目标提取方法 |
CN113034638A (zh) * | 2021-03-22 | 2021-06-25 | 山东大学 | 一种基于定向平滑约束反演的场地污染物刻画方法及*** |
CN113902779A (zh) * | 2021-09-28 | 2022-01-07 | 燕山大学 | 一种基于张量投票方法的点云配准方法 |
CN114428295A (zh) * | 2020-09-24 | 2022-05-03 | 中国石油化工股份有限公司 | 一种基于断层置信度参数控制的边缘保持扩散滤波方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102629367A (zh) * | 2012-01-17 | 2012-08-08 | 安徽建筑工业学院 | 基于KDTree的点云数据双边滤波去噪的方法 |
CN103530899A (zh) * | 2013-10-10 | 2014-01-22 | 浙江万里学院 | 基于几何特征性的点云简化方法 |
CN106530238A (zh) * | 2016-09-21 | 2017-03-22 | 燕山大学 | 一种散乱点云的保特征滤波方法 |
US20170193692A1 (en) * | 2015-12-30 | 2017-07-06 | Shenzhen Institutes Of Advanced Technology Chinese Academy Of Sciences | Three-dimensional point cloud model reconstruction method, computer readable storage medium and device |
CN107219483A (zh) * | 2017-04-22 | 2017-09-29 | 天津大学 | 一种基于扩散峰度成像的径向峰度各项异性定量方法 |
-
2017
- 2017-12-06 CN CN201711278418.9A patent/CN108022221A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102629367A (zh) * | 2012-01-17 | 2012-08-08 | 安徽建筑工业学院 | 基于KDTree的点云数据双边滤波去噪的方法 |
CN103530899A (zh) * | 2013-10-10 | 2014-01-22 | 浙江万里学院 | 基于几何特征性的点云简化方法 |
US20170193692A1 (en) * | 2015-12-30 | 2017-07-06 | Shenzhen Institutes Of Advanced Technology Chinese Academy Of Sciences | Three-dimensional point cloud model reconstruction method, computer readable storage medium and device |
CN106530238A (zh) * | 2016-09-21 | 2017-03-22 | 燕山大学 | 一种散乱点云的保特征滤波方法 |
CN107219483A (zh) * | 2017-04-22 | 2017-09-29 | 天津大学 | 一种基于扩散峰度成像的径向峰度各项异性定量方法 |
Non-Patent Citations (5)
Title |
---|
MIN KI PARK E.T.: "Multi-scale tensor voting for feature extraction from unstructured point clouds", 《GRAPHICAL MODELS》 * |
廖昌粟 等: "一种新的三维点云模型去噪光滑算法研究", 《计算机应用与软件》 * |
杨千里 等: "三维各向异性扩散滤波在地震数据处理中的应用", 《地球物理学进展》 * |
王丽坤: "三维各向异性扩散滤波方法研究与应用", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
章大勇,吴文启: "《激光雷达/惯性组合导航***的一致性与最优估计问题研究》", 28 February 2017, 北京:国防工业出版社 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109508489A (zh) * | 2018-11-07 | 2019-03-22 | 山东大学 | 一种各向异性多孔结构的建模方法及*** |
CN114428295A (zh) * | 2020-09-24 | 2022-05-03 | 中国石油化工股份有限公司 | 一种基于断层置信度参数控制的边缘保持扩散滤波方法 |
CN114428295B (zh) * | 2020-09-24 | 2024-03-29 | 中国石油化工股份有限公司 | 一种基于断层置信度参数控制的边缘保持扩散滤波方法 |
CN112381823A (zh) * | 2021-01-13 | 2021-02-19 | 北京智拓视界科技有限责任公司 | 用于图像的几何特征的提取方法和相关产品 |
CN112381823B (zh) * | 2021-01-13 | 2021-04-16 | 北京智拓视界科技有限责任公司 | 用于图像的几何特征的提取方法和相关产品 |
CN112837331A (zh) * | 2021-03-08 | 2021-05-25 | 电子科技大学 | 一种基于自适应形态重建模糊三维sar图像目标提取方法 |
CN113034638A (zh) * | 2021-03-22 | 2021-06-25 | 山东大学 | 一种基于定向平滑约束反演的场地污染物刻画方法及*** |
CN113902779A (zh) * | 2021-09-28 | 2022-01-07 | 燕山大学 | 一种基于张量投票方法的点云配准方法 |
CN113902779B (zh) * | 2021-09-28 | 2024-05-28 | 燕山大学 | 一种基于张量投票方法的点云配准方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108022221A (zh) | 一种各项异性扩散滤波的三维散乱点云平滑去噪方法 | |
Wang et al. | Design automation for customized apparel products | |
CN107067473A (zh) | 对3d建模对象进行重构 | |
CN110930334B (zh) | 一种基于神经网络的网格去噪方法 | |
US20180033181A1 (en) | Audio-based caricature exaggeration | |
JP2016157429A (ja) | 細分割曲面上における2d画像のエングレービング | |
Chen et al. | Non‐Oriented MLS Gradient Fields | |
Lee et al. | Adaptive tetrahedral element generation and refinement to improve the quality of bulk metal forming simulation | |
CN103530472A (zh) | 基于重要性采样的三维模型自动化简化方法 | |
CN106991204A (zh) | 对受约束的非对称细分网格的修改 | |
US20210141869A1 (en) | Automated Analysis of Mechanical Designs | |
JP2017504890A5 (zh) | ||
US20160140255A1 (en) | Method and apparatus for modeling deformable body by fusing surface particles and internal skeletal structure | |
Rohmer et al. | Folded paper geometry from 2D pattern and 3D contour | |
CN103077554B (zh) | 基于节单位的作物交互式设计方法和*** | |
Dabbaghchian et al. | Simplification of vocal tract shapes with different levels of detail | |
Wu et al. | Real‐Time Symmetry‐Preserving Deformation | |
CN103970952B (zh) | 一种面向三维制造技术的快速应力分析方法 | |
Li et al. | Four‐Dimensional Geometry Lens: A Novel Volumetric Magnification Approach | |
Du et al. | High fidelity virtual try-on network via semantic adaptation and distributed componentization | |
CN104637091B (zh) | 一种通过稀疏表示创建流形的方法 | |
CN107038741B (zh) | 三维渲染二维光影的方法 | |
Wu et al. | Diffusion equations over arbitrary triangulated surfaces for filtering and texture applications | |
Su et al. | Discrete Lie flow: A measure controllable parameterization method | |
EP2856433B1 (fr) | Procédé de traitement de données pour transformer une première surface 3d adaptée pour revêtir un premier objet en une deuxième surface 3d adaptée pour revêtir un deuxième objet, programme d'ordinateur et dispositif de traitement associés |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180511 |
|
RJ01 | Rejection of invention patent application after publication |