CN107784638A - 一种优化的东巴古籍图像增强方法 - Google Patents
一种优化的东巴古籍图像增强方法 Download PDFInfo
- Publication number
- CN107784638A CN107784638A CN201711031437.1A CN201711031437A CN107784638A CN 107784638 A CN107784638 A CN 107784638A CN 201711031437 A CN201711031437 A CN 201711031437A CN 107784638 A CN107784638 A CN 107784638A
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- msub
- image
- 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
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000005457 optimization Methods 0.000 title claims abstract description 15
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 18
- 230000002708 enhancing effect Effects 0.000 claims abstract description 10
- 238000006243 chemical reaction Methods 0.000 claims abstract description 8
- 230000008569 process Effects 0.000 claims abstract description 4
- 230000009466 transformation Effects 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000011426 transformation method Methods 0.000 claims description 4
- 206010003694 Atrophy Diseases 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 230000037444 atrophy Effects 0.000 claims description 3
- 230000001419 dependent effect Effects 0.000 claims description 3
- 244000062793 Sorghum vulgare Species 0.000 claims 1
- 235000019713 millet Nutrition 0.000 claims 1
- 235000015170 shellfish Nutrition 0.000 claims 1
- 230000008901 benefit Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000005728 strengthening Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000004807 localization Effects 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/10004—Still image; Photographic image
-
- 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/20—Special algorithmic details
- G06T2207/20024—Filtering details
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明涉及一种优化的东巴古籍图像增强方法,其步骤:对原始东巴古籍图像作图像类型的转换,将RGB图像转换为灰度图像;对灰度图像进行小波分解,选取小波分解后的db4小波基将图像进行二层分解,得到对应的小波低频系数和高频系数;对小波二层分解后得到的高频系数采用改进的阈值滤波算法处理,实现对小波分解后的高频系数的去噪处理;对小波二级分解后得到的低频系数进行模糊集增强处理,得到增强的低频子带系数;对处理后的低频系数和高频系数利用现有技术中的小波反变换进行重构,完成对东巴古籍图像的增强。本发明能在增强图像细节的同时抑制噪声,并能有效提高图像对比度。
Description
技术领域
本发明涉及一种图像处理领域,特别是关于一种优化的东巴古籍图像增强方法。
背景技术
图像增强是东巴经典古籍数字化处理过程中极其重要的环节之一。其基本思想是:通过采用一系列技术去改善图像的质量和视觉效果,突出图像中感兴趣的特征,获取图像中有价值的信息,使图像更方便人和机器进行分析处理。目前,频域法和空域法是两种主要的图像增强方法。频域法图像增强包括Curvelet变换、Contourlet变换、小波变换和傅里叶变换等方法。空域法图像增强包括多尺度Retinex算法、直方图均衡化、中值滤波和灰度变换等方法。虽然国内外研究图像增强算法的学者,无论是从数学理论方面、图像视觉效果方面,还是从图像特征方面,都提出了很多种算法,但面对不同的对象,其都存在一定的局限性。
图像增强技术的关键在于如何有效的改善图像质量和视觉效果的同时更好的保留图像的边缘和细节信息,现常用的两种图像增强方法为:1)直方图均衡:该方法通常用来增加许多图像的全局对比度,尤其是当图像的有用数据的对比度相当接近的时候。通过这种方法,亮度可以更好地在直方图上分布,从而可以用于增强局部的对比度而不影响整体的对比度。但是,该方法会增加背景噪声的对比度并且降低有用信号的对比度,也可能导致过增强与细节丢失问题(如图4所示)。2)小波阈值滤波:小波滤波与传统滤波相比,因小波变换具有时频局部化特性和多分辨的特性,使其在去噪方面有独特的优势,但其处理后的图像清晰度和对比度都很低,存在细节丢失的情况(如图5所示)。
发明内容
针对上述问题,本发明的目的是提供一种优化的东巴古籍图像增强方法,其能在增强图像细节的同时抑制噪声,并能有效提高图像对比度。
为实现上述目的,本发明采取以下技术方案:一种优化的东巴古籍图像增强方法,其特征在于包括以下步骤:(1)对原始东巴古籍图像作图像类型的转换,将RGB图像转换为灰度图像;(2)对灰度图像进行小波分解,选取小波分解后的db4小波基将图像进行二层分解,得到对应的小波低频系数和高频系数;(3)对小波二层分解后得到的高频系数采用改进的阈值滤波算法处理,得到改进的贝叶斯萎缩阈值,实现对小波分解后的高频系数的去噪处理;(4)对小波二级分解后得到的低频系数进行模糊集增强处理,得到增强的低频子带系数;(5)对处理后的低频系数和高频系数利用现有技术中的小波反变换进行重构,完成对东巴古籍图像的增强。
进一步,所述步骤(3)中,改进的阈值滤波算法具体步骤如下:(3.1)采用贝叶斯萎缩法估计噪声与信号的阈值;(3.2)计算第k尺度、第s个高频子带的噪声标准差(3.3)计算无噪声标准方差,即子带的信号标准差(3.4)在步骤(3.1)中的噪声与信号的阈值表达式中加上一个自适应变量ω,得到改进的贝叶斯萎缩阈值(3.5)将步骤(3.4)中得到的阈值代入阈值函数W(x),实现对小波分解后的高频系数进行去噪处理。
进一步,所述步骤(3.1)中,假设小波变换得到的高频系数总体上服从广义高斯分布,那么根据贝叶斯准则,得到贝叶斯萎缩阈值为:
式中,表示子带的信号标准差;为第k尺度、第s个高频子带的噪声标准差。
进一步,所述步骤(3.4)中,改进的贝叶斯萎缩阈值
其中,LJ表示第J层子带的长度。
进一步,所述步骤(3.5)中,将改进的贝叶斯萎缩阈值代入阈值函数W(x)有:
式中,sgn(x)表示返回函数,即x>0,sgn(x)=1,x=0,sgn(x)=0,x<0,sgn(x)=-1;x表示原系数;a表示引入的参数。
进一步,所述步骤(4)中,具体步骤如下:(4.1)定义大小为M×N的一幅图像的模糊矩阵A为:
其中,μij/xij表示图像中点(i,j)的系数值xij具有某种特征的程度为μij;(4.2)将低频子带系数通过隶属度函数映射到模糊空间:
其中;T(xij)表示低频子带(i,j)处系数值xij在模糊空间内构成的矩阵;xij表示低频子带(i,j)处系数值;xmin表示低频子带系数中的最小值、xmax表示低频子带系数中的最大值;(4.3)对矩阵T(xij)通过非线性变换实现模糊增强处理,得到新的模糊特征平面;(4.4)通过对新的模糊特征平面进行逆变换,得到增强的低频子带系数x′ij。
进一步,所述步骤(4.3)中,通过一个对比度增强算子来实现非线性变换方法增强图像,其表达式如下:
其中,μ′ij表示增强后的点(i,j)处具有某种特征的程度;T(μij)表示表示在新的模糊平面点(i,j)处具有某种特征的程度;T为低频子带中的阈值;当μij>T时,增强隶属度μij;反之,μij<T时则压缩隶属度μij。
进一步,所述步骤(4.4)中,采用的逆变换式如下:
x′ij=T-1(μij),
x′ij=(xmax-xmin)×μ′ij+xmin。
本发明由于采取以上技术方案,其具有以下优点:1、本发明结合了模糊集增强和小波阈值滤波的优点,不仅增强了图像细节和提高了对比度,而且取得了很好的降噪效果。2、本发明在软阈值函数和硬阈值函数之间引入参数a,针对不同的图像,通过调整参数a的大小,来获得合适的阈值函数。适当缩小了含噪系数与阈值的差值,优化了传统算法。3、本发明在模糊空间,通过一个对比度增强算子以此提高图像对比度不但运算速度快,而且可以保存图像中的低灰度信息。
附图说明
图1是本发明的整体流程示意图;
图2是本发明实施例中一幅原始输入的东巴古籍图像;
图3是对图2中的东巴古籍图像采用本发明方法进行图像增强显示的结果;
图4是对图2中的东巴古籍图像采用采用现有直方图均衡进行图像增强显示的结果;
图5是对图2中的东巴古籍图像采用现有小波阈值滤波进行图像增强显示的结果。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明提供一种优化的东巴古籍图像增强方法,其包括以下步骤:
(1)如图2所示,对原始东巴古籍图像作图像类型的转换,将RGB图像转换为灰度图像;
(2)对灰度图像进行小波分解,选取小波分解后的db4小波基将图像进行二层分解,得到对应的小波低频系数和高频系数;
(3)对小波二层分解后得到的高频系数采用改进的阈值滤波算法处理,得到改进的Bayes(贝叶斯)萎缩阈值,即噪声与信号的阈值,进而抑制噪声。
改进的阈值滤波算法具体步骤如下:
(3.1)采用贝叶斯萎缩法估计噪声与信号的阈值。假设小波变换得到的高频系数总体上服从广义高斯分布,那么根据贝叶斯准则,得到贝叶斯萎缩阈值为:
式中,表示子带的信号标准差;为第k尺度、第s个高频子带的噪声标准差。
(3.2)计算第k尺度、第s个高频子带的噪声标准差为:
其中,为第k尺度、第s个高频子带(i,j)位置上的高频系数;median()表示求解该子带系数的中值,其中,下标n表示高频系数的数量。
(3.3)计算无噪声标准方差,即子带的信号标准差为:
其中,为相应子带系数的方差。
(3.4)采用局部阈值法,依据分解层数的不同及其变化规律,设定不同的局部阈值。在步骤(3.1)中的Bayes(贝叶斯)萎缩阈值表达式中加上一个自适应变量ω,得到改进的贝叶斯萎缩阈值
其中,LJ表示第J层子带的长度。
(3.5)将步骤(3.4)中得到的阈值代入阈值函数W(x),实现对小波分解后的高频系数进行去噪处理:
式中,W(x)表示阈值函数;sgn(x)表示返回函数,即x>0,sgn(x)=1,x=0,sgn(x)=0,x<0,sgn(x)=-1;x表示原系数;a表示引入的参数,针对不同的图像,通过调整a的大小,来获得合适的阈值函数。
当原始图像较平滑,背景部分较多而细节部分较少时,a在[0.5,1]区域内取值;当原图纹理复杂,小波系数高频部分较重要时,a取低于0.5的值。
(4)对小波二级分解后得到的低频系数进行模糊集增强处理,得到增强的低频子带系数;
具体步骤如下:
(4.1)定义大小为M×N的一幅图像的模糊矩阵A为:
其中,μij/xij表示图像中点(i,j)的系数值xij具有某种特征的程度为μij。
(4.2)将低频子带系数通过隶属度函数映射到模糊空间:
其中;T(xij)表示低频子带(i,j)处系数值xij在模糊空间内构成的矩阵;xij表示低频子带(i,j)处系数值;xmin表示低频子带系数中的最小值、xmax表示低频子带系数中的最大值。
(4.3)在模糊空间,为了提高图像对比度,通过非线性变换方法增强图像。即对矩阵T(xij)通过非线性变换实现模糊增强处理,得到新的模糊特征平面;
通过一个对比度增强算子来实现非线性变换方法增强图像,其表达式如下:
其中,μ′ij表示增强后的点(i,j)处具有某种特征的程度;T(μij)表示表示在新的模糊平面点(i,j)处具有某种特征的程度;T为低频子带中的阈值,由于在进行步骤(2)小波分解后,低频子带包含了大量的图像基本信息,且基本滤除了噪声,故此处阈值T选用原始贝叶斯萎缩阈值对结果影响较小。当μij>T时,增强隶属度μij;反之,μij<T时则压缩隶属度μij。
(4.4)通过对新的模糊特征平面进行逆变换,得到增强的低频子带系数x′ij。
采用的逆变换式如下:
x′ij=T-1(μij),
x′ij=(xmax-xmin)×μ′ij+xmin。
(5)对处理后的低频系数和高频系数利用现有技术中的小波反变换进行重构,完成对东巴古籍图像的增强。
综上所述,本发明在使用时,相比于现有技术中采用直方图均衡对东巴古籍图像进行图像增强效果(如图4所示),或采用小波阈值滤波对东巴古籍图像进行图像增强的效果(如图5所示),本发明的东巴古籍图像增强效果较好(如图3所示),本发明不仅增强了图像细节和提高了对比度,而且取得了很好的降噪效果。
上述各实施例仅用于说明本发明,各部件的结构、尺寸、设置位置及形状都是可以有所变化的,在本发明技术方案的基础上,凡根据本发明原理对个别部件进行的改进和等同变换,均不应排除在本发明的保护范围之外。
Claims (8)
1.一种优化的东巴古籍图像增强方法,其特征在于包括以下步骤:
(1)对原始东巴古籍图像作图像类型的转换,将RGB图像转换为灰度图像;
(2)对灰度图像进行小波分解,选取小波分解后的db4小波基将图像进行二层分解,得到对应的小波低频系数和高频系数;
(3)对小波二层分解后得到的高频系数采用改进的阈值滤波算法处理,得到改进的贝叶斯萎缩阈值,实现对小波分解后的高频系数的去噪处理;
(4)对小波二级分解后得到的低频系数进行模糊集增强处理,得到增强的低频子带系数;
(5)对处理后的低频系数和高频系数利用现有技术中的小波反变换进行重构,完成对东巴古籍图像的增强。
2.如权利要求1所述的一种优化的东巴古籍图像增强方法,其特征在于:所述步骤(3)中,改进的阈值滤波算法具体步骤如下:
(3.1)采用贝叶斯萎缩法估计噪声与信号的阈值;
(3.2)计算第k尺度、第s个高频子带的噪声标准差
(3.3)计算无噪声标准方差,即子带的信号标准差
(3.4)在步骤(3.1)中的噪声与信号的阈值表达式中加上一个自适应变量ω,得到改进的贝叶斯萎缩阈值
(3.5)将步骤(3.4)中得到的阈值代入阈值函数W(x),实现对小波分解后的高频系数进行去噪处理。
3.如权利要求2所述的一种优化的东巴古籍图像增强方法,其特征在于:所述步骤(3.1)中,假设小波变换得到的高频系数总体上服从广义高斯分布,那么根据贝叶斯准则,得到贝叶斯萎缩阈值为:
<mrow>
<msubsup>
<mi>T</mi>
<mi>k</mi>
<mi>s</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mi>n</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>,</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>,</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
式中,表示子带的信号标准差;为第k尺度、第s个高频子带的噪声标准差。
4.如权利要求2所述的一种优化的东巴古籍图像增强方法,其特征在于:所述步骤(3.4)中,改进的贝叶斯萎缩阈值
<mrow>
<msubsup>
<mi>T</mi>
<mi>k</mi>
<msup>
<mi>s</mi>
<mo>&prime;</mo>
</msup>
</msubsup>
<mo>=</mo>
<mi>&omega;</mi>
<mfrac>
<mrow>
<msubsup>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mi>n</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>,</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>,</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中,LJ表示第J层子带的长度。
5.如权利要求2所述的一种优化的东巴古籍图像增强方法,其特征在于:所述步骤(3.5)中,将改进的贝叶斯萎缩阈值代入阈值函数W(x)有:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>W</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>sgn</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>*</mo>
<mi>max</mi>
<mo>{</mo>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mi>x</mi>
<mo>|</mo>
<mo>-</mo>
<msubsup>
<mi>aT</mi>
<mi>k</mi>
<msup>
<mi>s</mi>
<mo>&prime;</mo>
</msup>
</msubsup>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>0</mn>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>a</mi>
<mo>&Element;</mo>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
式中,sgn(x)表示返回函数,即x>0,sgn(x)=1,x=0,sgn(x)=0,x<0,sgn(x)=-1;x表示原系数;a表示引入的参数。
6.如权利要求1所述的一种优化的东巴古籍图像增强方法,其特征在于:所述步骤(4)中,具体步骤如下:
(4.1)定义大小为M×N的一幅图像的模糊矩阵A为:
<mrow>
<mi>A</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&mu;</mi>
<mn>11</mn>
</msub>
<mo>/</mo>
<msub>
<mi>x</mi>
<mn>11</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&mu;</mi>
<mrow>
<mn>1</mn>
<mi>N</mi>
</mrow>
</msub>
<mo>/</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mi>N</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>M</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>/</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>M</mi>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>M</mi>
<mi>N</mi>
</mrow>
</msub>
<mo>/</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>M</mi>
<mi>N</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,μij/xij表示图像中点(i,j)的系数值xij具有某种特征的程度为μij;
(4.2)将低频子带系数通过隶属度函数映射到模糊空间:
<mrow>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>T</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
其中;T(xij)表示低频子带(i,j)处系数值xij在模糊空间内构成的矩阵;xij表示低频子带(i,j)处系数值;xmin表示低频子带系数中的最小值、xmax表示低频子带系数中的最大值;
(4.3)对矩阵T(xij)通过非线性变换实现模糊增强处理,得到新的模糊特征平面;
(4.4)通过对新的模糊特征平面进行逆变换,得到增强的低频子带系数x′ij。
7.如权利要求6所述的一种优化的东巴古籍图像增强方法,其特征在于:所述步骤(4.3)中,通过一个对比度增强算子来实现非线性变换方法增强图像,其表达式如下:
<mrow>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<mi>T</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mfrac>
<msubsup>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mi>T</mi>
</mfrac>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo><</mo>
<mi>T</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>T</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>></mo>
<mi>T</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,μ′ij表示增强后的点(i,j)处具有某种特征的程度;T(μij)表示表示在新的模糊平面点(i,j)处具有某种特征的程度;T为低频子带中的阈值;当μij>T时,增强隶属度μij;反之,μij<T时则压缩隶属度μij。
8.如权利要求6所述的一种优化的东巴古籍图像增强方法,其特征在于:所述步骤(4.4)中,采用的逆变换式如下:
x′ij=T-1(μij),
x′ij=(xmax-xmin)×μ′ij+xmin。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711031437.1A CN107784638A (zh) | 2017-10-27 | 2017-10-27 | 一种优化的东巴古籍图像增强方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711031437.1A CN107784638A (zh) | 2017-10-27 | 2017-10-27 | 一种优化的东巴古籍图像增强方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107784638A true CN107784638A (zh) | 2018-03-09 |
Family
ID=61432331
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711031437.1A Pending CN107784638A (zh) | 2017-10-27 | 2017-10-27 | 一种优化的东巴古籍图像增强方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107784638A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108961179A (zh) * | 2018-06-19 | 2018-12-07 | 上海中和软件有限公司 | 一种医学图像后处理系统及其使用方法 |
CN109615590A (zh) * | 2018-11-14 | 2019-04-12 | 江苏科技大学 | 基于模糊算法与分数阶微分算法的声呐图像增强方法 |
CN110009583A (zh) * | 2019-03-28 | 2019-07-12 | 上海交通大学 | 基于灰度均衡和nsct分解的红外图像预处理方法 |
CN110097518A (zh) * | 2019-04-28 | 2019-08-06 | 东软医疗系统股份有限公司 | 图像去噪方法、装置及终端设备 |
CN110390652A (zh) * | 2019-07-24 | 2019-10-29 | 东北大学 | 一种模糊集与分数阶微分相结合的医疗图像增强方法 |
CN112862716A (zh) * | 2021-02-07 | 2021-05-28 | 江苏科技大学 | 基于双树复数小波和模糊理论的声呐图像目标增强方法 |
CN112991196A (zh) * | 2021-02-02 | 2021-06-18 | 武汉科技大学 | 一种回转窑火焰图像频域去噪方法 |
CN113763267A (zh) * | 2021-08-26 | 2021-12-07 | 长春理工大学 | 一种基于nsct图像融合的强散射环境下图像复原方法 |
CN114202491A (zh) * | 2021-12-08 | 2022-03-18 | 深圳市研润科技有限公司 | 一种增强光学图像的方法及系统 |
CN115187465A (zh) * | 2022-05-30 | 2022-10-14 | 中国矿业大学 | 基于nsct和图像匹配的浮选泡沫图像增强方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101799916A (zh) * | 2010-03-16 | 2010-08-11 | 刘国传 | 基于贝叶斯估计的生物芯片图像小波去噪方法 |
CN102314675A (zh) * | 2011-09-27 | 2012-01-11 | 西安电子科技大学 | 基于小波高频的贝叶斯去噪方法 |
CN106408539A (zh) * | 2016-09-26 | 2017-02-15 | 微美光速资本投资管理(北京)有限公司 | 一种基于模糊逻辑的增强现实图像再现方法 |
-
2017
- 2017-10-27 CN CN201711031437.1A patent/CN107784638A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101799916A (zh) * | 2010-03-16 | 2010-08-11 | 刘国传 | 基于贝叶斯估计的生物芯片图像小波去噪方法 |
CN102314675A (zh) * | 2011-09-27 | 2012-01-11 | 西安电子科技大学 | 基于小波高频的贝叶斯去噪方法 |
CN106408539A (zh) * | 2016-09-26 | 2017-02-15 | 微美光速资本投资管理(北京)有限公司 | 一种基于模糊逻辑的增强现实图像再现方法 |
Non-Patent Citations (5)
Title |
---|
GENG QIAOMAN ET AL: "Application of adaptive median filter and wavelet transform to dongba manuscript images denoising", 《2017 IEEE 13TH INTERNATIONAL CONFERENCE ON ELECTRONIC MEASUREMENT & INSTRUMENTS》 * |
全永奇 等: "模糊集与非线性增益相结合的自适应图像增强算法", 《计算机应用研究》 * |
刘兴淼 等: "基于小波变换与模糊理论的图像增强算法研究", 《弹箭与制导学报》 * |
李金 等: "一种改进的小波域贝叶斯图像去噪方法", 《模式识别与仿真》 * |
郭海涛 等: "贝叶斯改进阈值超声图像去噪方法", 《应用声学》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108961179A (zh) * | 2018-06-19 | 2018-12-07 | 上海中和软件有限公司 | 一种医学图像后处理系统及其使用方法 |
CN108961179B (zh) * | 2018-06-19 | 2022-10-18 | 上海中和软件有限公司 | 一种医学图像后处理系统及其使用方法 |
CN109615590A (zh) * | 2018-11-14 | 2019-04-12 | 江苏科技大学 | 基于模糊算法与分数阶微分算法的声呐图像增强方法 |
CN110009583A (zh) * | 2019-03-28 | 2019-07-12 | 上海交通大学 | 基于灰度均衡和nsct分解的红外图像预处理方法 |
CN110097518A (zh) * | 2019-04-28 | 2019-08-06 | 东软医疗系统股份有限公司 | 图像去噪方法、装置及终端设备 |
CN110097518B (zh) * | 2019-04-28 | 2022-12-27 | 东软医疗系统股份有限公司 | 图像去噪方法、装置及终端设备 |
CN110390652A (zh) * | 2019-07-24 | 2019-10-29 | 东北大学 | 一种模糊集与分数阶微分相结合的医疗图像增强方法 |
CN112991196A (zh) * | 2021-02-02 | 2021-06-18 | 武汉科技大学 | 一种回转窑火焰图像频域去噪方法 |
CN112862716A (zh) * | 2021-02-07 | 2021-05-28 | 江苏科技大学 | 基于双树复数小波和模糊理论的声呐图像目标增强方法 |
CN113763267A (zh) * | 2021-08-26 | 2021-12-07 | 长春理工大学 | 一种基于nsct图像融合的强散射环境下图像复原方法 |
CN114202491A (zh) * | 2021-12-08 | 2022-03-18 | 深圳市研润科技有限公司 | 一种增强光学图像的方法及系统 |
CN115187465A (zh) * | 2022-05-30 | 2022-10-14 | 中国矿业大学 | 基于nsct和图像匹配的浮选泡沫图像增强方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107784638A (zh) | 一种优化的东巴古籍图像增强方法 | |
CN101950414B (zh) | 自然图像非局部均值去噪方法 | |
CN101944230B (zh) | 基于多尺度的自然图像非局部均值去噪方法 | |
CN112233026A (zh) | 一种基于多尺度残差注意力网络的sar图像去噪方法 | |
CN103761719B (zh) | 一种基于邻域相关性的自适应小波阈值去噪方法 | |
CN103700072A (zh) | 一种基于自适应小波阈值和双边滤波器的图像去噪方法 | |
CN107464226A (zh) | 一种基于改进二维经验模态分解算法的图像去噪方法 | |
CN107909558A (zh) | 一种基于无监督学习的非局部均值图像去噪方法 | |
CN107610074B (zh) | 一种提高遥感图像质量的方法 | |
CN106570843A (zh) | 一种自适应小波阈值函数图像噪声抑制方法 | |
CN104881847A (zh) | 一种基于小波分析和伪彩色处理的比赛视频图像增强方法 | |
CN107784639A (zh) | 一种无人机遥感图像改进多边滤波去噪方法 | |
Gupta et al. | A novel approach for image denoising and performance analysis using SGO and APSO | |
CN104732498B (zh) | 一种基于非下采样Contourlet变换的阈值化图像去噪方法 | |
CN111242854A (zh) | 一种图像去噪方法 | |
CN104616259B (zh) | 一种噪声强度自适应的非局部均值图像去噪方法 | |
CN102222327A (zh) | 基于Treelet变换和最小均方误差估计的图像去噪方法 | |
CN101957984B (zh) | 基于非局部萎缩因子参数估计的图像去噪方法 | |
CN106067164A (zh) | 基于自适应小波域处理的彩色图像对比度增强算法 | |
CN104240208A (zh) | 非制冷红外焦平面探测器图像细节增强方法 | |
CN109003247A (zh) | 彩色图像混合噪声的去除方法 | |
CN113204051A (zh) | 一种基于变分模态分解的低秩张量地震数据去噪方法 | |
CN105023257B (zh) | 基于N‑Smoothlets的图像去噪方法 | |
CN106023103A (zh) | 一种基于精确局部方差先验建模的自适应正交小波图像去噪方法 | |
CN107230191B (zh) | 一种基于结构相似性的非局部均值去噪优化方法 |
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: 20180309 |
|
RJ01 | Rejection of invention patent application after publication |