CN113160088B - 一种基于Sobel算子及图像熵的散斑干涉相位图滤波评价方法 - Google Patents
一种基于Sobel算子及图像熵的散斑干涉相位图滤波评价方法 Download PDFInfo
- Publication number
- CN113160088B CN113160088B CN202110492892.1A CN202110492892A CN113160088B CN 113160088 B CN113160088 B CN 113160088B CN 202110492892 A CN202110492892 A CN 202110492892A CN 113160088 B CN113160088 B CN 113160088B
- Authority
- CN
- China
- Prior art keywords
- filtering
- image
- interference phase
- speckle interference
- gradient
- 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
- 238000001914 filtration Methods 0.000 title claims abstract description 121
- 238000011156 evaluation Methods 0.000 title claims abstract description 16
- 238000010587 phase diagram Methods 0.000 claims abstract description 35
- 238000010586 diagram Methods 0.000 claims abstract description 27
- 230000009466 transformation Effects 0.000 claims abstract description 20
- 238000009499 grossing Methods 0.000 claims abstract description 9
- 230000009467 reduction Effects 0.000 claims abstract description 7
- 230000008859 change Effects 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 23
- 230000008569 process Effects 0.000 abstract description 7
- 238000005259 measurement Methods 0.000 abstract description 3
- 238000011158 quantitative evaluation Methods 0.000 abstract 1
- 239000002131 composite material Substances 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
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
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Processing (AREA)
Abstract
本发明涉及一种基于Sobel算子及图像熵的散斑干涉相位图滤波评价方法,所述方法包括对散斑干涉相位图进行正、余弦变化,并对其进行均值滤波平滑处理,而后利用得到的一次滤波后的所述正弦图和余弦图进行反正切变换,获取一次滤波后的散斑干涉相位图,之后与Sobel算子进行卷积得到散斑干涉相位图的梯度图,通过对所述梯度图求取一维图像熵,最后重复滤波过程,得到多个所述一维图像熵,并进行统计绘制成折线图,当折线图变化达到平台期,滤波停止,获得滤波后的降噪散斑干涉相位图。本发明在单一散斑相位图持续滤波过程中避免了欠滤波或过滤波的问题,对滤波进行定量评价,控制滤波过程,提高了数字散斑测量智能化程度。
Description
技术领域
本发明涉及激光无损检测和光学图像处理领域,具体涉及一种基于Sobel算子及图像熵的散斑干涉相位图滤波评价方法。
背景技术
近年来,随着航空航天的发展,各种复合材料获得广泛应用。然而,这些复合材料在加工制造和服役过程中由于变形、位移等缺陷导致构件性能显著下降,因此服役时间减少。为此,需要对这些复合材料进行检测及评估,为后续抢修提供方案,保障航空航天设备的健康运行。
研究表明当复合材料被与相干的光照射时,物体表面发生形变或者位移时,物面的形变就转化成像面上散斑的相位变化,这种现象称为散斑干涉。数字散斑干涉(DigitalSpeckle Pattern Interferometry,DSPI)是一种非接触式全场测量技术,它可以对复合材料的位移、形变、表面缺陷等进行测量。采集的散斑相位图由于大量噪声干扰,导致散斑相位图信噪比低、相位测量灵敏度低,无法获得测量物体的精确信息。因此,有必要对散斑相位图进行处理,提高图像的信噪比,去除噪声干扰。
散斑干涉相位图常用滤波方法有均值滤波法、中值滤波法、小波变换法。均值滤波法是一种线性滤波法,利用模板进行计算获得新像素点的灰度值,新的灰度值是模板内以目标像素点为中心周围的若干个像素灰度值的平均值。均值滤波法在滤除噪声的同时破坏了图像的细节部分,导致处理后的图像变得模糊,使相位突变的信息丢失;
中值滤波法是一种非线性滤波法,将目标像素点的邻域内各像素点灰度值的中值用来代替该像素点的灰度值。与线性滤波方法相比,中值滤波方法在降噪的同时可以保留边缘信息,且方法简单,但对高斯噪声滤波效果不明显,平滑效果较差;
小波变换法是对图像进行变换,将图像从空域变换到频域,在频域对图像进行降噪处理,小波变换法是将图像分解成不同的频带,利用预先设置的阈值获取有用频带,然后进行重建获取降噪后的图像,该方法在使用时较为繁琐。
发明内容
本发明为有效解决对数字散斑相位图中大量噪声干扰时易出现的欠滤波和过滤波的问题,提供了一种基于Sobel算子及图像熵的散斑干涉相位图滤波评价方法,采用Sobel算子对散斑相位图进行处理,无需参数设置,可以对图片的滤波次数进行判断,获得精确的测量信息。
为了实现上述目的,本发明提出了一种基于Sobel算子及图像熵的散斑干涉相位图滤波评价方法,所述方法包括以下步骤:
步骤1:对散斑干涉相位图进行正、余弦变化,得到散斑干涉相位图的正弦图和余弦图;
对所述正弦图和余弦图分别进行均值滤波平滑处理,得到一次滤波后的正弦图和余弦图;
步骤2:利用一次滤波后的所述正弦图和余弦图进行反正切变换,得到一次滤波后的散斑干涉相位图;
步骤3:基于一次滤波后的所述散斑干涉相位图利用Sobel算子进行卷积,得到一次滤波后的散斑干涉相位图的梯度图;
步骤4:对所述梯度图求取一维图像熵,得到所述一次滤波后的散斑干涉相位图的梯度图的一维图像熵;
步骤5:重复步骤1~4,得到多个所述一维图像熵,并进行统计绘制成折线图,当折线图变化达到平台期,滤波停止,获得滤波后的降噪散斑干涉相位图;
步骤3具体包括:
将一次滤波后的所述散斑干涉相位图分别与水平方向上的Sobel算子模板和垂直方向上的Sobel算子模板进行卷积,获得其水平方向上的梯度图和垂直方向上的梯度图;
分别对水平方向上的梯度图和垂直方向上的梯度图求模并相加,获得一次滤波后的散斑干涉相位图的梯度图;
通过公式(6)~(9)得到一次滤波后的散斑干涉相位图的梯度图:
其中:gx和gy分别为水平和垂直方向的Sobel算子,为使计算方便,近似取其梯度图的模为水平分量模与垂直分量模之和;
在x方向上的梯度的模为:
在y方向上的梯度的模为:
将x方向和y方向上的梯度相加求和,获得滤波后图像整体的梯度灰度值矩阵为:
步骤4具体包括:
统计一次滤波后的所述散斑干涉相位图像梯度图中各灰度值出现的概率,并对其取对数;
针对各灰度值,令其出现概率与出现概率的对数相乘,并求和取相反数;
如公式(10)所示:
其中,pi表示灰度值为i的点在S(f(n+1)(x,y))中出现的概率。
进一步地,步骤1具体包括:
步骤1.1:针对图像每个像素点的灰度值,分别进行正弦变换和余弦变换,计算出原图像的正弦变换图像和余弦变换图像,如公式(1)和(2)所示:
其中,-π<f(n)(x,y)<+π,f(n)为散斑干涉相位图;
步骤1.2:通过均值滤波平滑处理对每一个像素点根据给定的窗口大小,以该像素点为中心,计算窗口内所有像素点的灰度值均值,以均值替换原本像素点灰度值,如公式(3)和(4)所示:
进一步地,步骤2具体包括:
将一次滤波后的正弦图与余弦图逐像素点进行相除,对相除后的图像进行反正切变换,获得一次滤波后的散斑干涉相位图,如公式(5)所示:
进一步地,步骤5具体包括:
重复对图像进行滤波,统计历次滤波后梯度图像的一维图像熵,并画出折线图;
当第N次滤波得到的一维图像熵小于第N-1次滤波得到的一维图像熵,且第N次滤波得到的一维图像熵小于第N+1滤波得到的一维图像熵时,滤波停止;
如公式(11)所示:
H(S(f(n-1)(x,y)))>H(S(f(n)(X,y)));
H(S(f(n)(x,y)))<H(S(f(n+1)(x,y))) (11);
获取最后一次滤波的图像,即滤波后的降噪散斑干涉相位图。
通过上述技术方案,本发明的有益效果为:
本发明首先对散斑干涉相位图进行正、余弦变化,得到散斑干涉相位图的正弦图和余弦图,接着对所述正弦图和余弦图分别进行均值滤波平滑处理,得到一次滤波后的正弦图和余弦图,利用一次滤波后的所述正弦图和余弦图进行反正切变换,得到一次滤波后的散斑干涉相位图,并基于一次滤波后的所述散斑干涉相位图利用Sobel算子进行卷积,得到一次滤波后的散斑干涉相位图的梯度图,进而对所述梯度图求取一维图像熵,得到所述一次滤波后的散斑干涉相位图的梯度图的一维图像熵,最后重复步骤1~4,得到多个所述一维图像熵,并进行统计绘制成折线图,当折线图变化达到平台期,滤波停止,获得滤波后的降噪散斑干涉相位图。
本发明针对Sobel算子与图像卷积的特点,提出的求取散斑相位图梯度图的方法,对于滤波效果的评价较为直观;
对图像熵代表信息量的特点,本发明提出的计算滤波过程中散斑相位图一维图像熵的的方法,获得持续滤波过程中一维图像熵变化的折线图,能自动判断滤波是否完成;
本发明能够自适应的对散斑干涉相位图片进行降噪,避免传统滤波法繁琐的滤波次数设置。
附图说明
图1为本发明一种基于Sobel算子及图像熵的散斑干涉相位图滤波评价方法的流程图;
图2为本发明一种基于Sobel算子及图像熵的散斑干涉相位图滤波评价方法的原理结构图;
图3为本发明一种基于Sobel算子及图像熵的散斑干涉相位图滤波评价方法中包含噪声的散斑干涉相位图;
图4为利用本发明一种基于Sobel算子及图像熵的散斑干涉相位图滤波评价方法后得到的降噪散斑干涉相位图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
如图1所示,一种基于Sobel算子及图像熵的散斑干涉相位图滤波评价方法,所述方法包括:
步骤1:对散斑干涉相位图进行正、余弦变化,得到散斑干涉相位图的正弦图和余弦图;
对所述正弦图和余弦图分别进行均值滤波平滑处理,得到一次滤波后的正弦图和余弦图;
步骤2:利用一次滤波后的所述正弦图和余弦图进行反正切变换,得到一次滤波后的散斑干涉相位图;
步骤3:基于一次滤波后的所述散斑干涉相位图利用Sobel算子进行卷积,得到一次滤波后的散斑干涉相位图的梯度图;
步骤4:对所述梯度图求取一维图像熵,得到所述一次滤波后的散斑干涉相位图的梯度图的一维图像熵;
步骤5:重复步骤1~4,得到多个所述一维图像熵,并进行统计绘制成折线图,当折线图变化达到平台期,滤波停止,获得滤波后的降噪散斑干涉相位图;
步骤3具体包括:
将一次滤波后的所述散斑干涉相位图分别与水平方向上的Sobel算子模板和垂直方向上的Sobel算子模板进行卷积,获得其水平方向上的梯度图和垂直方向上的梯度图;
分别对水平方向上的梯度图和垂直方向上的梯度图求模并相加,获得一次滤波后的散斑干涉相位图的梯度图;
通过公式(6)~(9)得到一次滤波后的散斑干涉相位图的梯度图:
其中:gx和gy分别为水平和垂直方向的Sobel算子,为使计算方便,近似取其梯度图的模为水平分量模与垂直分量模之和;
在x方向上的梯度的模为:
在y方向上的梯度的模为:
将x方向和y方向上的梯度相加求和,获得滤波后图像整体的梯度灰度值矩阵为:
步骤4具体包括:
统计一次滤波后的所述散斑干涉相位图像梯度图中各灰度值出现的概率,并对其取对数;
针对各灰度值,令其出现概率与出现概率的对数相乘,并求和取相反数;
如公式(10)所示:
其中,pi表示灰度值为i的点在S(f(n+1)(x,y))中出现的概率。
本实施例实现了在包含噪声的相位图前提下可以对散斑干涉相位图进行均值滤波,根据Sobel算子求得的梯度图,计算梯度图的一维图像熵,统计梯度图一维图像熵的变化情况,获得滤波后的散斑相位图。
实施例2
基于上述实施例对步骤1进行优化具体的:
步骤1.1:针对图像每个像素点的灰度值,分别进行正弦变换和余弦变换,计算出原图像的正弦变换图像和余弦变换图像,如公式(1)和(2)所示:
其中,-π<f(n)(x,y)<+π,f(n)为散斑干涉相位图;
步骤1.2:通过均值滤波平滑处理对每一个像素点根据给定的窗口大小,以该像素点为中心,计算窗口内所有像素点的灰度值均值,以均值替换原本像素点灰度值,如公式(3)和(4)所示:
从而得到一次滤波后的正弦图和余弦图。
实施例3
基于上述实施例1对步骤2进行优化具体的:
步骤2具体包括:
将一次滤波后的正弦图与余弦图逐像素点进行相除,对相除后的图像进行反正切变换,获得一次滤波后的散斑干涉相位图,如公式(5)所示:
实施例4
基于上述多个实施例,为了实现自适应的对散斑干涉相位图片进行降噪,对步骤5进行优化,具体包括:
重复对图像进行滤波,统计历次滤波后梯度图像的一维图像熵,并画出折线图;
当第N次滤波得到的一维图像熵小于第N-1次滤波得到的一维图像熵,且第N次滤波得到的一维图像熵小于第N+1滤波得到的一维图像熵时,滤波停止;
如公式(11)所示:
H(S(f(n-1)(x,y)))>H(S(f(n)(x,y)));
H(S(f(n)(x,y)))<H(S(f(n+1)(x,y))) (11);
获取最后一次滤波的图像,即滤波后的降噪散斑干涉相位图。
如图2~4所示通过本发明方法,实现了对包含大量噪声的散斑干涉相位图进行滤波,无需提前设置滤波次数,通过Sobel算子与散斑相位图像卷积获得梯度图,根据梯度图一维图像熵判断滤波完成情况,降低散斑相位图像中的噪声。
以上所述之实施例,只是本发明的较佳实施例而己,并非限制本发明的实施范围,故凡依本发明专利范围所述的构造、特征及原理所做的等效变化或修饰,均应包括于本发明申请专利范围内。
Claims (4)
1.一种基于Sobel算子及图像熵的散斑干涉相位图滤波评价方法,其特征在于,所述方法包括:
步骤1:对散斑干涉相位图进行正、余弦变化,得到散斑干涉相位图的正弦图和余弦图;
对所述正弦图和余弦图分别进行均值滤波平滑处理,得到一次滤波后的正弦图和余弦图;
步骤2:利用一次滤波后的所述正弦图和余弦图进行反正切变换,得到一次滤波后的散斑干涉相位图;
步骤3:基于一次滤波后的所述散斑干涉相位图利用Sobel算子进行卷积,得到一次滤波后的散斑干涉相位图的梯度图;
步骤4:对所述梯度图求取一维图像熵,得到所述一次滤波后的散斑干涉相位图的梯度图的一维图像熵;
步骤5:重复步骤1~4,得到多个所述一维图像熵,并进行统计绘制成折线图,当折线图变化达到平台期,滤波停止,获得滤波后的降噪散斑干涉相位图;
步骤3具体包括:
将一次滤波后的所述散斑干涉相位图分别与水平方向上的Sobel算子模板和垂直方向上的Sobel算子模板进行卷积,获得其水平方向上的梯度图和垂直方向上的梯度图;
分别对水平方向上的梯度图和垂直方向上的梯度图求模并相加,获得一次滤波后的散斑干涉相位图的梯度图;
通过公式(6)~(9)得到一次滤波后的散斑干涉相位图的梯度图:
其中:gx和gy分别为水平和垂直方向的Sobel算子,为使计算方便,近似取其梯度图的模为水平分量模与垂直分量模之和;
在x方向上的梯度的模为:
在y方向上的梯度的模为:
将x方向和y方向上的梯度相加求和,获得滤波后图像整体的梯度灰度值矩阵为:
步骤4具体包括:
统计一次滤波后的所述散斑干涉相位图像梯度图中各灰度值出现的概率,并对其取对数;
针对各灰度值,令其出现概率与出现概率的对数相乘,并求和取相反数;
如公式(10)所示:
其中,pi表示灰度值为i的点在S(f(n+1)(x,y))中出现的概率。
4.根据权利要求1所述的一种基于Sobel算子及图像熵的散斑干涉相位图滤波评价方法,其特征在于,步骤5具体包括:
重复对图像进行滤波,统计历次滤波后梯度图像的一维图像熵,并画出折线图;
当第N次滤波得到的一维图像熵小于第N-1次滤波得到的一维图像熵,且第N次滤波得到的一维图像熵小于第N+1滤波得到的一维图像熵时,滤波停止;
如公式(11)所示:
H(S(f(n-1)(x,y)))>H(S(f(n)(X,y)));
H(S(f(n)(x,y)))<H(S(f(n+1)(x,y))) (11);
获取最后一次滤波的图像,即滤波后的降噪散斑干涉相位图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110492892.1A CN113160088B (zh) | 2021-04-30 | 2021-04-30 | 一种基于Sobel算子及图像熵的散斑干涉相位图滤波评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110492892.1A CN113160088B (zh) | 2021-04-30 | 2021-04-30 | 一种基于Sobel算子及图像熵的散斑干涉相位图滤波评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113160088A CN113160088A (zh) | 2021-07-23 |
CN113160088B true CN113160088B (zh) | 2022-08-12 |
Family
ID=76873553
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110492892.1A Active CN113160088B (zh) | 2021-04-30 | 2021-04-30 | 一种基于Sobel算子及图像熵的散斑干涉相位图滤波评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113160088B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114018150B (zh) * | 2021-11-08 | 2023-08-15 | 河北荣众信息科技有限公司 | 一种数字全息显微测量装置及其控制方法 |
CN114299402B (zh) * | 2022-03-07 | 2022-05-20 | 成都理工大学 | 隐患点自动识别方法、电子设备及计算机可读存储介质 |
CN116429002A (zh) * | 2023-03-01 | 2023-07-14 | 北京航空航天大学 | 一种强粘接、耐高温散斑制备方法及制备系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104902267A (zh) * | 2015-06-08 | 2015-09-09 | 浙江科技学院 | 一种基于梯度信息的无参考图像质量评价方法 |
CN107292866A (zh) * | 2017-05-17 | 2017-10-24 | 浙江科技学院 | 一种基于相对梯度的图像质量客观评价方法 |
CN107563969A (zh) * | 2017-08-03 | 2018-01-09 | 天津大学 | 基于变分模态分解的dspi相位滤波方法 |
CN108648154A (zh) * | 2018-04-27 | 2018-10-12 | 合肥工业大学 | 相位图的滤波评价方法 |
CN110766627A (zh) * | 2019-10-16 | 2020-02-07 | 北京信息科技大学 | 一种散斑干涉图像的降噪方法和装置 |
CN112381809A (zh) * | 2020-11-19 | 2021-02-19 | 天津中科虹星科技有限公司 | 一种基于相对熵的虹膜图像清晰度评价方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3955899B2 (ja) * | 2002-09-20 | 2007-08-08 | 国立大学法人埼玉大学 | 電子的スペックル干渉法を用いた変形計測方法および装置 |
-
2021
- 2021-04-30 CN CN202110492892.1A patent/CN113160088B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104902267A (zh) * | 2015-06-08 | 2015-09-09 | 浙江科技学院 | 一种基于梯度信息的无参考图像质量评价方法 |
CN107292866A (zh) * | 2017-05-17 | 2017-10-24 | 浙江科技学院 | 一种基于相对梯度的图像质量客观评价方法 |
CN107563969A (zh) * | 2017-08-03 | 2018-01-09 | 天津大学 | 基于变分模态分解的dspi相位滤波方法 |
CN108648154A (zh) * | 2018-04-27 | 2018-10-12 | 合肥工业大学 | 相位图的滤波评价方法 |
CN110766627A (zh) * | 2019-10-16 | 2020-02-07 | 北京信息科技大学 | 一种散斑干涉图像的降噪方法和装置 |
CN112381809A (zh) * | 2020-11-19 | 2021-02-19 | 天津中科虹星科技有限公司 | 一种基于相对熵的虹膜图像清晰度评价方法 |
Non-Patent Citations (9)
Title |
---|
"An Image Quality Assessment Algorithm Based on Feature Selection";Ting Lu等;《IScIDE 2012》;20131231;第289-297页 * |
"Nonreference Image Quality Evaluation Algorithm Based on Wavelet Convolutional Neural Network and Information Entropy";Jinhua Liu等;《Entropy》;20191031;第1-18页 * |
"The simulation of adaptive optical image even and pulse noise and research of image quality evaluation";Wen Changli等;《International Symposium on Photoelectronic Detection and Imaging 2013》;20131231;第1-7页 * |
"基于信息熵和细节方差均值与背景方差均值比的无参考图像锐化结果评价";张玉娟等;《哈尔滨师范大学自然科学学报》;20190131;第36-40页 * |
"基于平滑样条拟合的相位图像滤波评价方法";赵琪涵等;《光学学报》;20180831;第1-7页 * |
"基于梯度误差电阻抗断层成像技术的图像评价";吕婧华等;《中国组织工程研究与临床康复》;20081223;第10277-10280页 * |
"基于能量的DSPI相位图正余弦降噪";肖启阳等;《光学精密工程》;20180430;第936-943页 * |
"数字图像相关法中散斑图的质量评价";冯益春等;《计算机与数字工程》;20200731;第1743-1747页 * |
"数字散斑干涉相位信息处理关键技术研究";肖启阳;《万方数据企业知识服务平台》;20210204;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113160088A (zh) | 2021-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113160088B (zh) | 一种基于Sobel算子及图像熵的散斑干涉相位图滤波评价方法 | |
CN101661611B (zh) | 基于贝叶斯非局部均值滤波器的实现方法 | |
CN109191387B (zh) | 一种基于巴特沃斯滤波器的红外图像去噪方法 | |
CN101901476A (zh) | 基于NSCT域边缘检测和Bishrink模型的SAR图像去噪方法 | |
Xiong et al. | Application of histogram equalization for image enhancement in corrosion areas | |
CN106169174B (zh) | 一种图像放大方法 | |
CN102968770A (zh) | 噪声消除方法及装置 | |
CN109118434A (zh) | 一种图像预处理方法 | |
CN107194889A (zh) | 分块双边全变分正则化图像噪声消除方法 | |
CN111986098A (zh) | 一种含固定背景的被动式太赫兹图像增强方法 | |
Ma et al. | An automatic detection method of Mura defects for liquid crystal display | |
CN113899349B (zh) | 海浪参数检测方法、设备及存储介质 | |
CN104616259A (zh) | 一种噪声强度自适应的非局部均值图像去噪方法 | |
CN103455986B (zh) | 基于分数阶微分梯度的随机噪声点检测方法 | |
CN116012265B (zh) | 基于时空域自适应滤波的红外视频降噪方法和装置 | |
CN116883987A (zh) | 一种用于变电站无人巡检的指针仪表读数识别方法 | |
CN111899196A (zh) | 一种基于经典复原算法的叶片缺陷运动模糊图像复原方法 | |
CN107907542B (zh) | 一种ivmd及能量估计相结合的dspi相位滤波方法 | |
CN103325090A (zh) | 一种相干斑抑制方法和装置 | |
CN115829967A (zh) | 一种工业金属表面缺陷图像去噪和增强方法 | |
CN115937302A (zh) | 结合边缘保持的高光谱图像亚像元定位方法 | |
Paul et al. | MR image enhancement using an extended neighborhood filter | |
Dasgupta | Comparative analysis of non-blind deblurring methods for noisy blurred images | |
Daniel et al. | Retinal image enhancement using wavelet domain edge filtering and scaling | |
Mar et al. | Guided filter smoothing for third order edge mask |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230821 Address after: No. 217-225, Zone A3, Yilian Building Materials Home Furnishing Hardware City, Yandu Street, Ye County, Pingdingshan City, Henan Province, 467200 Patentee after: Henan Haolun cloud Technology Co.,Ltd. Address before: Shunhe District in Henan province 475001 City Minglun Street No. 85 Patentee before: Henan University |