CN1809837A - 图像处理装置、图像处理方法和程序及记录介质 - Google Patents
图像处理装置、图像处理方法和程序及记录介质 Download PDFInfo
- Publication number
- CN1809837A CN1809837A CNA2004800174722A CN200480017472A CN1809837A CN 1809837 A CN1809837 A CN 1809837A CN A2004800174722 A CNA2004800174722 A CN A2004800174722A CN 200480017472 A CN200480017472 A CN 200480017472A CN 1809837 A CN1809837 A CN 1809837A
- Authority
- CN
- China
- Prior art keywords
- image data
- value
- aforementioned
- conversion coefficient
- output
- 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
- 238000003672 processing method Methods 0.000 title claims description 11
- 238000006243 chemical reaction Methods 0.000 claims abstract description 59
- 230000002708 enhancing effect Effects 0.000 claims description 42
- 238000001914 filtration Methods 0.000 claims description 9
- 238000009499 grossing Methods 0.000 claims description 4
- 238000000034 method Methods 0.000 description 29
- 230000009466 transformation Effects 0.000 description 17
- 238000010586 diagram Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 8
- 230000000694 effects Effects 0.000 description 6
- 230000000873 masking effect Effects 0.000 description 6
- 238000009795 derivation Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000001965 increasing effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- UNPLRYRWJLTVAE-UHFFFAOYSA-N Cloperastine hydrochloride Chemical compound Cl.C1=CC(Cl)=CC=C1C(C=1C=CC=CC=1)OCCN1CCCCC1 UNPLRYRWJLTVAE-UHFFFAOYSA-N 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000003707 image sharpening Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004335 scaling law Methods 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 239000002023 wood Substances 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/73—Deblurring; Sharpening
- G06T5/75—Unsharp masking
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
- H04N1/40—Picture signal circuits
- H04N1/409—Edge or detail enhancement; Noise or error suppression
- H04N1/4092—Edge or detail enhancement
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
- Compression Or Coding Systems Of Tv Signals (AREA)
Abstract
提供可以同时实现噪声去除和轮廓增强的图像处理。高通滤波器(1)输出输入图像数据f(m,n)的高通分量h(m,n)。增强控制量导出部(2)对输入图像数据进行离散小波变换,从而求出增强控制量e(m,n)。乘法部(3)对每个像素求出并输出来自增强控制量导出部(2)的增强控制量e(m,n)与作为来自高通滤波器(1)的输出的高通分量h(m,n)的相乘值。放大部(5)使乘法部(3)的输出为常数倍(λ倍)。加法部(4)将来自乘法部(3)的相乘值与输入图像数据相加,输出输出图像数据。
Description
技术领域
本发明涉及一种图像处理装置、图像处理方法和程序及记录介质,特别是涉及一种去除图像中的噪声,并且增强轮廓部分的亮度对比度的图像处理装置、图像处理方法和程序及记录介质。
背景技术
以往,作为以图像增强为目的广泛利用的方法,通常有反锐化掩模(unsharp masking)。这是如下的方法:对于输入图像数据,利用滤波器求出高频分量,乘以决定增强控制量的任意常数,与输入图像数据相加,输出增强图像。但是,反锐化掩模法,在输入图像数据中叠加有噪声的情况下,不能避免噪声分量的放大。另外,在噪声叠加图像的增强中,原图像的噪声分量被放大,出现在增强结果中。
为了解决这一点,提出了对高频分量,局部控制增强控制量的方法。在非专利文献1中,提出了如下方法:对每个输入图像数据的像素,将从亮度值的水平方向的中央差分、垂直方向的中央差分的平方和所导出的亮度级与高频分量相乘,从而进行增强控制量控制。在非专利文献2中,对非专利文献1的方法进行改良,导入根据模糊规则进行的控制。另外,在专利文献1中,提出了为实现噪声去除而对于小振幅的高频分量使符号翻转与输入图像数据相加的噪声去除装置。
非专利文献1
G.Ramponi,“A cubic unsharp masking technique forcontrast enhancement,”Signal Processing,vol.67,pp.211-222,June 1998.
非专利文献2
木村誠聡、田口亮、村田裕,“フアジ一推論を用いた雑音が重畳された画像の強調に関する一手法”,信学論A、vol.J81-A、no.9、pp.1247-1256、1998年9月
专利文献1
日本特开2001-274995号公报
发明内容
然而,非专利文献1的方法是利用在图像的轮廓部分亮度变化大、只在轮廓部分使高频分量的振幅增大的方法,但是在轮廓部分的亮度变化小的情况下、以及在噪声振幅大的情况下效果较低。另外,不能够去除叠加在输入图像数据中的噪声。另外,非专利文献2是基于与非专利文献1相同的亮度级的方法,因此,具有与非专利文献1同样的问题。并且,专利文献1具有亮度变化小的轮廓被判断为噪声而被平滑化的问题。
在本发明中,鉴于以上点,目的在于提供可以同时实现噪声去除和轮廓增强的图像处理装置、图像处理方法和程序及记录介质。
根据本发明的第一解决方法,提供一种图像处理装置,去除输入图像的噪声且增强轮廓部分的对比度,其具备:
滤波器,使输入图像数据的高频分量通过;
导出部,对输入图像数据进行离散小波变换,从而求出在图像轮廓部分和噪声中成为不同的大小关系的第一及第二变换系数,根据第一变换系数的平方的值、第一及第二变换系数的积的值、以及预先决定的设定值,求出增强控制量;
乘法部,输出来自前述导出部的增强控制量与来自前述滤波器的输出的相乘值;
加法部,将来自前述乘法部的相乘值与输入图像数据相加而得到输出图像数据。
另外,前述导出部可以具备:
离散小波变换部,对输入图像数据进行离散小波变换,求出第一及第二变换系数;
第一电路,具有对第一变换系数进行平方的平方电路;
第二电路,具有将第一及第二变换系数相乘的乘法器;
设定部,计算并输出使前述第一电路的输出为预先决定的α倍的值、使前述第二电路的输出为预先决定的β倍的值、和预先决定的γ值的线性和。
根据本发明的第二解决方法,提供一种图像处理方法,去除输入图像的噪声且增强轮廓部分的对比度,其包含:
对输入图像数据进行离散小波变换,从而求出在图像轮廓部分和噪声中成为不同的大小关系的第一及第二变换系数,根据第一变换系数的平方的值、第一及第二变换系数的积的值、以及预先决定的设定值,求出增强控制量;
输出前述增强控制量与输入图像数据的高频分量的相乘值;
将前述相乘值和输入图像数据相加而得到输出图像数据的加法部。
根据本发明的第三解决方法,提供一种去除输入图像的噪声且增强轮廓部分的对比度的图像处理程序以及记录了该程序的计算机可读取的记录介质,用于使计算机执行以下步骤:
处理部从存储部或者输入部读取输入图像数据;
处理部对输入图像数据进行离散小波变换,从而求出在图像轮廓部分和噪声中成为不同的大小关系的第一及第二变换系数,根据第一变换系数的平方的值、第一及第二变换系数的积的值、以及预先决定的设定值,求出增强控制量;
处理部输出前述增强控制量与输入图像数据的高频分量的相乘值;
处理部将前述相乘值与输入图像数据相加而得到输出图像数据的加法部;
处理部将得到的输出图像数据存储到存储部以及/或者输出到输出部或显示部。
附图说明
图1是图像处理装置的结构图。
图2是高通滤波器的滤波系数的说明图。
图3是增强控制量导出部2的结构图。
图4是最大尺度为2时的小波变换的说明图。
图5是表示对一维信号进行小波变换的三个滤波器用的系数的例子的说明图。
图6是用于实现二维小波变换的滤波器组的结构图。
图7是表示离散二进小波变换的高通滤波器、低通滤波器的结构的图。
图8是一维信号的小波变换与小波变换的尺度间的积的说明图。
图9是表示输入图像数据例的图。
图10是表示根据作为现有方法的反锐化掩模法得到的处理结果的图。
图11是表示根据非专利文献1所示的方法得到的增强图像的处理结果的图。
图12是表示利用本发明同时进行了轮廓部分的增强和噪声去除的处理结果的图。
图13是与本实施方式相关的硬件的结构图。
图14是图像处理的流程图。
具体实施方式
1.图像处理装置
图1示出了图像处理装置的结构图。
本实施方式的图像处理装置具备高通滤波器1、增强控制量导出部2、乘法部3、加法部4、放大部5,去除输入图像的噪声且增强轮廓部分的对比度。
高通滤波器1,使输入图像数据f(m,n)的高频分量通过,输出高通分量h(m,n)。此外,作为输入图像数据f(m,n),高通滤波器1及增强控制量导出部2等处理中所需的、坐标(m,n)的对象像素的周围的像素也作为输入图像数据被适当输入。增强控制量导出部2,根据输入图像数据对每个像素输出增强控制量e(m,n)。增强控制量导出部2,对输入图像数据进行离散小波变换,从而求出在图像轮廓部分和噪声中成为不同的大小关系的第一及第二变换系数,根据第一变换系数的平方的值、第一及第二变换系数的积的值、以及预先决定的设定值,求出增强控制量e(m,n)。乘法部3,对每个像素求出并输出来自增强控制量导出部2的增强控制量e(m,n)和作为来自高通滤波器1的输出的高通分量h(m,n)的相乘值。另外,根据需要设置放大部5,从而可使乘法部3的输出为常数倍(λ倍),向加法部4输出λe(m,n)h(m,n)。该常数λ是决定增强程度的正的常数,可以在处理前决定。加法部4,将来自乘法部3的相乘值和输入图像数据相加,得到输出图像数据f(m,n)+λe(m,n)h(m,n),输出其结果。
高通滤波器1
图2是高通滤波器1的滤波系数的说明图。作为高通滤波器1的例子,有具有如图所示的滤波系数的拉普拉斯滤波器。
增强控制量导出部2
图3是增强控制量导出部2的结构图。在该增强控制量导出部2的例子中,从具有两个尺度的离散二进小波变换导出增强控制量e(m,n)。
增强控制量导出部2具备离散小波变换部21、第一及第二平方电路22及23、第一加法器24、第一及第二乘法器25及26、第二加法器27、设定部28、限幅器29。
离散小波变换部21,对输入图像数据进行离散小波变换,求出第一及第二尺度的各水平及垂直方向变换系数。第一平方电路22对第一尺度的水平方向的变换系数进行平方,第二平方电路23对第一尺度的垂直方向的变换系数进行平方。第一加法器24将第一及第二平方电路22及23的输出相加。另外,第一乘法器25将第一及第二尺度的水平方向的变换系数相乘,第二乘法器26将第一及第二尺度的垂直方向的变换系数相乘。第二加法器27将第一及第二乘法器25及26的输出相加。设定部28计算并输出使第一加法器24的输出为预定的α倍的值、使第二加法器27的输出为预定的β倍的值、和预定的γ值的线性和。限幅器29对计算的线性和的数值范围进行限制。
在本实施方式中,为了这样导出增强控制量e(m,n),求出从输入图像数据的离散二进小波变换得到的变换系数的平方、在不同的尺度间得到的变换系数的积、以及常数的线性和。在设定部28中,通过选择线性和的权重,在图像轮廓部分中增强控制量成为正值、在图像平坦部分中成为负值。将它输入到决定增强控制量的下限和上限的限幅器29中,使限幅器29的输出为常数倍之后,利用乘法部3,与从输入图像数据的滤波器输出得到的高频分量相乘。在限幅器29的输出为负值的情况下,即在图像平坦部分,从输入图像数据减去高频分量,装置作为平滑化滤波器工作;在限幅器29的输出为正值的情况下,即在图像轮廓部分,将高频分量与输入图像数据相加,因此作为图像增强滤波器动作。
离散小波变换部21
在此,对离散小波变换部21进行说明。
通常,离散二进小波变换由多个小波函数和图像的卷积运算进行定义。另外,小波变换由将小波函数作为滤波系数来保持的数字滤波器实现。
小波函数通过将基本小波函数向时间轴方向延伸为2j倍来定义。在此,j是大于等于1的整数,称为尺度。在将尺度j的最大值设为J的情况下,小波变换输出从尺度1到尺度J为止的J个尺度和与它们对应的变换系数。
图4示出了使最大尺度为2的情况下的一维小波变换的说明图。
离散二进小波变换,作为一个例子,可以由图中所示的高通滤波器和低通滤波器构成的滤波器组结构来实现。即,对于一维信号的小波变换,在设定为J=2的情况下,由该图所示的滤波器的级联连接来实现。在高通滤波器1中,使用由基本小波函数决定的滤波系数。在高通滤波器2中,使用在高通滤波器的滤波系数的取样期间插入零的滤波系数,由低通滤波器对它进行插值,从而导出2倍尺度的小波系数。
图5是表示对这些一维信号进行小波变换的3个滤波器用的系数的例子的说明图。图5(a)表示高通滤波器1的滤波系数的例子,图5(b)表示高通滤波器2的滤波系数的例子,图5(c)表示低通滤波器的滤波系数的例子。
在这些图中,将各图的中心作为输出的中心,示出了与周围像素相乘的系数。另外,在该例子中,对于对象像素f(n),f(N)(N=n-2,n-1,n,n+1,n+2)作为输入图像数据被输入并使用。
根据上图(a)、(b)、(c)所示的滤波系数,对于输入信号x(n),尺度1的小波系数Wf1(n)被计算为
Wf1(n)=x(n-1)-x(n+1),
低通滤波器输出Sf1(n)被计算为
由此,尺度2的小波系数Wf1(n)根据Sf1(n)被计算为
Wf2(n)=Sf1(n-2)-Sf1(n+2)。
图6示出了用于实现二维小波变换的滤波器组的结构图。
离散小波变换部21具备高通滤波器1_61、高通滤波器2_62、高通滤波器3_63、高通滤波器4_64、低通滤波器65。
当是作为二维信号的图像数据时,通过交替地重复一维滤波处理来实现。
高通滤波器1_61具有图5(a)的滤波系数,高通滤波器4_64具有图5(b)的滤波系数,对图像的水平方向的每一行进行一维滤波处理。高通滤波器2_62具有图5(a)的滤波系数,高通滤波器3_63具有图5(b)的滤波系数,对图像的垂直方向的每一行进行一维滤波处理。低通滤波器65通过如下实现:将根据图2(c)的滤波系数的滤波处理,对水平方向的每一行执行后,对垂直方向的每一行执行。
高通滤波器1_61输出水平方向的第一变换系数,高通滤波器2_62输出垂直方向的第一变换系数。高通滤波器3_63输出水平方向的第二变换系数,高通滤波器4_64输出垂直方向的第二变换系数。
图7是表示离散二进小波变换的高通滤波器、低通滤波器的其他结构的图。图(a)表示高通滤波器1_61,图(b)表示高通滤波器3_63,图(c)表示高通滤波器2_62,图(d)表示高通滤波器4_64,图(e)表示低通滤波器65。高通滤波器1_61和高通滤波器2_62,为了在轮廓部分使变换系数产生最大值,作为例子可以使用该图所示的滤波系数。另外,在低通滤波器65中,为了满足小波变换的相似定律,作为例子可以使用该图所示的滤波系数。
设定部28
下面,对增强控制量导出部2的设定部28进行说明。在设定部28中,设定常数α、β、γ,计算:尺度1的垂直方向系数的平方值与水平方向系数的平方值的和的α倍的值、将尺度1和尺度2的垂直方向系数相乘的值与将尺度1和尺度2的水平方向的系数相乘的值的和的β倍的值、以及常数γ的总和e。
通常,可知在轮廓部分中,尺度1的小波变换系数小于等于尺度2的小波变换系数。另外,相反对于高斯性噪声等噪声,可知尺度1的小波变换系数比尺度2的小波变换系数大。因此,作为例子设定为α=-1、β=1、γ=0,从而可以只在图像轮廓部分使增强控制量为正值、在平坦部分为负值。另外,为了对所有的像素实现均匀的对比度增强,设定为α=0、β=0、γ=1。相反,为了对所有的像素均匀地应用平滑化,设定为α=0、β=0、γ=-1。不进行噪声去除而实现抑制了噪声放大的图像增强时,设定为α=0、β=1、γ=0。可以根据增强控制部的三个参数设定,同时实现噪声去除和增强、以及抑制了噪声放大的增强等各种效果。
图8示出了一维信号的小波变换与小波变换的尺度间的积的说明图。
该图中,对于小波变换的尺度和信号的关系,将一维信号作为例子示出。
在上图中,作为输入信号,作为一个例子使用如下信号:在从15点到45点为止的亮度值为100、除此之外为0的波形中,叠加了噪声。用曲线图表示用于从该输入信号导出增强控制量的小波变换尺度1的平方与小波变换尺度1和尺度2的积。
仅由小波变换尺度1的平方值,难以区分噪声和图像轮廓,但是由两个尺度的积,由于仅对噪声的小波变换系数的振幅随着尺度的增加而减少的性质,轮廓部分作为极大值明确地被显示。对该两个特征量乘以系数α、β,加上常数γ,设定α、β、γ使增强控制量在平坦部分中为负值、在图像轮廓部分中为正值,由此成为在包含噪声的平坦部分中从原图像减去高通分量的动作,可以从输入图像数据去除包含在高通分量中的噪声分量。作为例子设定为(α,β,γ)=(-1,1,0)时,在平坦部分中产生噪声去除效果。
限幅器29
根据上述计算求出的总和e,在作为正值过大时引起过增强,另外,在作为负值过小时引起噪声产生,因此,利用限幅器29限制数值范围。预先决定U为正的常数、L为负的常数。在限幅器29的输入e超过U的情况下,限幅器29输出U,在不足L的情况下输出L,除此之外的情况下将输入值e直接输出。在使用上图所示的拉普拉斯滤波器的情况下,使L的值为-0.2/λ时,可以保证本装置在图像平坦部分作为平滑化滤波器进行动作。
输入图像数据和增强控制量的导出
以下说明输入图像数据的像素、和使用于导出每个像素的增强控制量的像素。
以f(m,n)表示输入图像数据的坐标(m,n)中的亮度值,对于坐标(m,n)中的小波变换结果以及增强控制量,也如图中以数列所示,决定于每个坐标(m,n)。
用于图像处理的输入图像数据,是对象像素f(m,n)的周围的像素,包含二维小波变换的计算所需的像素。这些输入图像数据,例如,预先存储到存储器中,根据需要由导出部读出并使用。例如,在本实施方式中,对于对象像素f(m,n),在高通及低通滤波器中,如由图7等所示,使用滤波处理所需的周围的像素数据计算离散二进小波变换。另外,也可以递归地利用输出图像数据。
另外,对于增强控制量的导出中使用的像素的关系,决定增强控制量的小波系数的导出如上,并且,对于增强控制量e(m,n)
e(m,n)=F[α〔(W1 H(m,n))2+(W1 V(m,n))2〕+β〔W1 H(m,n)W2 H(m,n)+W1 V(m,n)W2 V(m,n)〕+γ]
在此F[·]由表示限幅器29的输入输出关系的非线性函数和明示的式子示出。
在本发明中,增强控制量根据离散二进小波变换的不同尺度中的变换系数的大小关系,成为正值及负值,在正值的情况下表现出增强的效果,在负值的情况下表现出噪声去除的效果。离散二进小波变换系数的大小关系,不依赖于图像的亮度的振幅及对比度,在图像轮廓和噪声中成为不同的关系,因此,可以不受轮廓部分的对比度的大小、噪声振幅的影响,只增强图像轮廓部分。
2.图像处理的具体例
为了表示本实施方式的噪声去除及轮廓增强的效果,使用了图9所示的输入图像数据。图9是对某一图像应用了3×3像素的平均值滤波器后,叠加了方差为50的高斯噪声的图像。
图10是表示根据作为现有方法的反锐化掩模法得到的处理结果的图,图11是表示根据非专利文献1所示的方法得到的增强图像的处理结果的图。另一方面,图12是表示利用本发明同时进行了轮廓部分的增强和噪声去除的处理结果的图。为得到图12的结果,作为一个例子,设定常数λ为0.001、α=-1、β=1、γ=0、L=-0.2/λ。
比较三个增强图像时,图像轮廓部分的增强为相同程度,但在根据本发明的增强图像(图12)中,图像背景部分的噪声与其他方法相比减少,可以确认本发明的有效性。
3.图像处理程序
本发明的图像处理方法或者图像处理装置/系统,可通过用于使计算机执行该图像处理方法或者图像处理装置/系统的各步骤的图像处理程序、记录了图像处理程序的计算机可读取的记录介质、包含图像处理程序并可加载到计算机的内部存储器中的程序制品、包含该程序的服务器等的计算机等来提供。
图13是与本实施方式相关的硬件的结构图。
该硬件具有作为中央处理装置(CPU)的处理部101、输入部102、输出部103、显示部104及存储部105。另外,处理部101、输入部102、输出部103、显示部104及存储部105,以星状或者总线等适当的连接方法被连接。存储部105包含:输入图像文件151,存储有被进行图像处理的输入图像数据f(m,n);增强控制量文件152,存储有被进行计算的增强控制量e(m,n);输出图像文件153,存储有被进行了图像处理的输出图像数据。
图14表示图像处理的流程图。另外,各处理的细节与“1.图像处理装置”中说明的一样。
图像处理程序使计算机执行以下的处理,从而去除输入图像的噪声且增强轮廓部分的对比度。首先,处理部101从存储部105的输入图像文件151或者输入部102读取输入图像数据(步骤S 1)。处理部101通过对输入图像数据进行离散小波变换,求出在图像轮廓部分和噪声中成为不同的大小关系的第一及第二变换系数,并根据第一变换系数的平方的值、第一及第二变换系数的积的值、以及预先决定的设定值,求出增强控制量(步骤S2)。处理部101根据需要,将求出的增强控制量存储到增强控制量文件152中。处理部101,输出增强控制量与输入图像数据的高频分量的相乘值(步骤S3)。处理部101,将相乘值和输入图像数据相加,得到输出图像数据的加法部(步骤S4)。处理部101,将得到的输出图像数据存储到存储部105的输出图像文件153中以及/或者输出到输出部103或显示部104(步骤S5)。另外,处理部101也可以根据得到的输出图像数据,递归地计算以上的图像处理,进一步求出输出图像数据。
产业上的可利用性
根据本发明,具有噪声去除能力,但提出方法的运算量的增加,与非专利文献1的方法相比,只有用于导出尺度2的小波变换系数所需的低通滤波器和高通滤波器3、4的计算、以及用于导出增强控制量的计算(每个像素乘8次和加4次、阈值操作2次)。考虑到以与该运算量的增加相同程度的运算量来实现只保存图像轮廓部分的噪声去除一直以来很难,根据本发明,可以实现用于同时进行增强和噪声去除的运算量的减少。
根据本发明,通过变更增强控制量导出部2的常数(α,β,γ),不需要变更装置整体,可以实现增强、噪声去除或者同时实现两者,得到各种增强特性。
现有的技术(非专利文献1、2及专利文献1),利用依靠亮度振幅及差分的方法来进行噪声和图像轮廓的分离,因此,输入图像数据的对比度低的情况及噪声振幅大的情况下,对于噪声放大的抑制效果降低,根据本发明,为了根据离散二进小波变换系数的尺度间的大小关系来实现噪声和图像轮廓的分离,可以不依靠输入图像数据的对比度及噪声振幅,实现噪声去除和图像增强。
本发明特别适用于例如数字照相机、数字视频照相机、图像扫描仪中的输入图像的图像清晰化。
Claims (15)
1.一种图像处理装置,去除输入图像的噪声且增强轮廓部分的对比度,其具备:
滤波器,使输入图像数据的高频分量通过;
导出部,对输入图像数据进行离散小波变换,从而求出在图像轮廓部分和噪声中成为不同的大小关系的第一及第二变换系数,根据第一变换系数的平方的值、第一及第二变换系数的积的值、以及预先决定的设定值,求出增强控制量;
乘法部,输出来自前述导出部的增强控制量与来自前述滤波器的输出的相乘值;
加法部,将来自前述乘法部的相乘值与输入图像数据相加而得到输出图像数据。
2.根据权利要求1所述的图像处理装置,其特征在于,
前述导出部在图像平坦部分输出负值,利用前述加法部从输入图像数据减去高频分量,另一方面,前述导出部在图像轮廓部分输出正值,利用前述加法部将高频分量与输入图像数据相加。
3.根据权利要求1所述的图像处理装置,其特征在于,前述导出部,具备:
离散小波变换部,对输入图像数据进行离散小波变换,求出第一及第二变换系数;
第一电路,具有对第一变换系数进行平方的平方电路;
第二电路,具有将第一及第二变换系数相乘的乘法器;
设定部,计算并输出使前述第一电路的输出为预先决定的α倍的值、使前述第二电路的输出为预先决定的β倍的值、预先决定的γ值的线性和。
4.根据权利要求3所述的图像处理装置,其特征在于,
前述设定部,设定为α=-1、β=1、γ=0,从而使增强控制量在图像轮廓部分为正值、在平坦部分为负值。
5.根据权利要求3所述的图像处理装置,其特征在于,
前述设定部,设定为α=0、β=0、γ=1,从而对于所有的像素均匀地实现对比度增强。
6.根据权利要求3所述的图像处理装置,其特征在于,
前述设定部,设定为α=0、β=0、γ=-1,从而对于所有的像素均匀地应用平滑化。
7.根据权利要求3所述的图像处理装置,其特征在于,
前述设定部,设定为α=0、β=1、γ=0,从而实现不进行噪声去除而抑制了噪声放大的图像增强。
8.根据权利要求1所述的图像处理装置,其特征在于,
前述导出部具备离散小波变换部,该离散小波变换部具有:
第一高通滤波器,对输入图像数据的水平方向的每一行进行一维滤波处理,输出水平方向的第一变换系数;
第二高通滤波器,对输入图像数据的垂直方向的每一行进行一维滤波处理,输出垂直方向的第一变换系数;
低通滤波器,对输入图像数据水平方向的每一行执行且对垂直方向的每一行执行;
第三高通滤波器,对来自前述低通滤波器的输出的水平方向的每一行进行一维滤波处理,输出水平方向的第二变换系数;
第四高通滤波器,对来自前述低通滤波器的输出的垂直方向的每一行进行一维滤波处理,输出垂直方向的第二变换系数。
9.根据权利要求1所述的图像处理装置,其特征在于,
前述导出部还具备限幅器,该限幅器用于限制所计算的线性和的数值范围。
10.根据权利要求1所述的图像处理装置,其特征在于,
还具备放大部,该放大部使前述乘法部的输出为常数倍并输出到前述加法部,从而决定增强的程度。
11.一种图像处理方法,去除输入图像的噪声且增强轮廓部分的对比度,其包含:
对输入图像数据进行离散小波变换,从而求出在图像轮廓部分和噪声中成为不同的大小关系的第一及第二变换系数,根据第一变换系数的平方的值、第一及第二变换系数的积的值、以及预先决定的设定值,求出增强控制量;
输出前述增强控制量与输入图像数据的高频分量的相乘值;
将前述相乘值和输入图像数据相加而得到输出图像数据的加法部。
12.根据权利要求11所述的图像处理方法,其特征在于,还包含:
计算并输出使第一变换系数的平方的值为预先决定的α倍的值、使第一及第二变换系数的积的值为预先决定的β倍的值、以及预先决定的γ的值的线性和。
13.根据权利要求12所述的图像处理方法,其特征在于,
作为前述设定值,设定为α=-1、β=1、γ=0,从而使增强控制量在图像轮廓部分为正值、在平坦部分为负值。
14.一种图像处理程序,去除输入图像的噪声且增强轮廓部分的对比度,用于使计算机执行以下步骤:
处理部从存储部或者输入部读取输入图像数据;
处理部对输入图像数据进行离散小波变换,从而求出在图像轮廓部分和噪声中成为不同的大小关系的第一及第二变换系数,根据第一变换系数的平方的值、第一及第二变换系数的积的值、以及预先决定的设定值,求出增强控制量;
处理部输出前述增强控制量与输入图像数据的高频分量的相乘值;
处理部将前述相乘值与输入图像数据相加而得到输出图像数据的加法部;
处理部将得到的输出图像数据存储到存储部以及/或者输出到输出部或显示部。
15.一种记录介质,是记录了去除输入图像的噪声且增强轮廓部分的对比度的图像处理程序的、计算机可读取的记录介质,用于使计算机执行以下步骤:
处理部从存储部或者输入部读取输入图像数据;
处理部对输入图像数据进行离散小波变换,从而求出在图像轮廓部分和噪声中成为不同的大小关系的第一及第二变换系数,根据第一变换系数的平方的值、第一及第二变换系数的积的值、以及预先决定的设定值,求出增强控制量;
处理部输出前述增强控制量与输入图像数据的高频分量的相乘值;
处理部将前述相乘值与输入图像数据相加而得到输出图像数据的加法部;
处理部将得到的输出图像数据存储到存储部以及/或者输出到输出部或显示部。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP270889/2003 | 2003-07-04 | ||
JP2003270889A JP2005025696A (ja) | 2003-07-04 | 2003-07-04 | 画像処理装置、画像処理方法及びプログラム及び記録媒体 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN1809837A true CN1809837A (zh) | 2006-07-26 |
Family
ID=33562632
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2004800174722A Pending CN1809837A (zh) | 2003-07-04 | 2004-03-24 | 图像处理装置、图像处理方法和程序及记录介质 |
Country Status (4)
Country | Link |
---|---|
US (1) | US20070009171A1 (zh) |
JP (1) | JP2005025696A (zh) |
CN (1) | CN1809837A (zh) |
WO (1) | WO2005004057A1 (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102113307A (zh) * | 2008-07-30 | 2011-06-29 | 三菱电机株式会社 | 图像处理装置、方法以及图像显示装置 |
CN102792124A (zh) * | 2009-10-17 | 2012-11-21 | 荷马里光电子公司 | 增强的成像方法及装置 |
CN102844786A (zh) * | 2010-03-01 | 2012-12-26 | 夏普株式会社 | 图像放大装置、图像放大程序、存储图像放大程序的存储介质和显示装置 |
CN103076339A (zh) * | 2013-01-17 | 2013-05-01 | 杭州电子科技大学 | 一种基于限幅均值滤波的电子清纱器纱疵快速识别方法 |
CN104883980A (zh) * | 2012-11-01 | 2015-09-02 | 日立阿洛卡医疗株式会社 | 医用图像诊断装置以及医用图像生成方法 |
CN109118442A (zh) * | 2018-07-20 | 2019-01-01 | 湖北工业大学 | 一种基于Sobel算子滤波的图像增强算法 |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101278315A (zh) | 2005-09-30 | 2008-10-01 | 皇家飞利浦电子股份有限公司 | 图像细节增强 |
US7738713B2 (en) * | 2006-08-31 | 2010-06-15 | Ching-Wei Yeh | Method for processing digital image with discrete wavelet transform and apparatus for the same |
JPWO2009037817A1 (ja) * | 2007-09-19 | 2011-01-06 | パナソニック株式会社 | 輪郭補正装置、輪郭補正方法、映像表示装置 |
US8326075B2 (en) | 2008-09-11 | 2012-12-04 | Google Inc. | System and method for video encoding using adaptive loop filter |
US8780971B1 (en) * | 2011-04-07 | 2014-07-15 | Google, Inc. | System and method of encoding using selectable loop filters |
US8781004B1 (en) | 2011-04-07 | 2014-07-15 | Google Inc. | System and method for encoding video using variable loop filter |
US8780996B2 (en) | 2011-04-07 | 2014-07-15 | Google, Inc. | System and method for encoding and decoding video data |
US8849057B2 (en) * | 2011-05-19 | 2014-09-30 | Foveon, Inc. | Methods for digital image sharpening with noise amplification avoidance |
US8885706B2 (en) | 2011-09-16 | 2014-11-11 | Google Inc. | Apparatus and methodology for a video codec system with noise reduction capability |
US9131073B1 (en) | 2012-03-02 | 2015-09-08 | Google Inc. | Motion estimation aided noise reduction |
US9344729B1 (en) | 2012-07-11 | 2016-05-17 | Google Inc. | Selective prediction signal filtering |
KR20150090515A (ko) * | 2014-01-29 | 2015-08-06 | 삼성전자주식회사 | 고주파수 성분의 위상 변조를 통한 영상 질감 향상 방법 및 그 장치 |
US10102613B2 (en) | 2014-09-25 | 2018-10-16 | Google Llc | Frequency-domain denoising |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1166280A (ja) * | 1997-08-25 | 1999-03-09 | Shimadzu Corp | 医用画像処理装置 |
JPH11252368A (ja) * | 1998-03-04 | 1999-09-17 | Ricoh Co Ltd | 画像処理装置 |
EP1064782A1 (en) * | 1999-01-15 | 2001-01-03 | Koninklijke Philips Electronics N.V. | Sharpness enhancement |
JP4244094B2 (ja) * | 1999-02-17 | 2009-03-25 | 富士フイルム株式会社 | 画像処理方法および装置並びに記録媒体 |
US6665448B1 (en) * | 2000-09-29 | 2003-12-16 | Hewlett-Packard Development Company, L.P. | Selective smoothing and sharpening of images by generalized unsharp masking |
-
2003
- 2003-07-04 JP JP2003270889A patent/JP2005025696A/ja active Pending
-
2004
- 2004-03-24 WO PCT/JP2004/004085 patent/WO2005004057A1/ja active Application Filing
- 2004-03-24 US US10/563,284 patent/US20070009171A1/en not_active Abandoned
- 2004-03-24 CN CNA2004800174722A patent/CN1809837A/zh active Pending
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102113307A (zh) * | 2008-07-30 | 2011-06-29 | 三菱电机株式会社 | 图像处理装置、方法以及图像显示装置 |
CN102113307B (zh) * | 2008-07-30 | 2014-04-30 | 三菱电机株式会社 | 图像处理装置、方法以及图像显示装置 |
CN102792124A (zh) * | 2009-10-17 | 2012-11-21 | 荷马里光电子公司 | 增强的成像方法及装置 |
CN102792124B (zh) * | 2009-10-17 | 2016-12-28 | 荷马里光电子公司 | 增强的成像方法及装置 |
CN102844786A (zh) * | 2010-03-01 | 2012-12-26 | 夏普株式会社 | 图像放大装置、图像放大程序、存储图像放大程序的存储介质和显示装置 |
CN104883980A (zh) * | 2012-11-01 | 2015-09-02 | 日立阿洛卡医疗株式会社 | 医用图像诊断装置以及医用图像生成方法 |
CN103076339A (zh) * | 2013-01-17 | 2013-05-01 | 杭州电子科技大学 | 一种基于限幅均值滤波的电子清纱器纱疵快速识别方法 |
CN109118442A (zh) * | 2018-07-20 | 2019-01-01 | 湖北工业大学 | 一种基于Sobel算子滤波的图像增强算法 |
CN109118442B (zh) * | 2018-07-20 | 2021-08-03 | 湖北工业大学 | 一种基于Sobel算子滤波的图像增强方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2005004057A1 (ja) | 2005-01-13 |
US20070009171A1 (en) | 2007-01-11 |
JP2005025696A (ja) | 2005-01-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1809837A (zh) | 图像处理装置、图像处理方法和程序及记录介质 | |
CN107358586B (zh) | 一种图像增强方法、装置及设备 | |
CN101076079A (zh) | 一种视频信号图像增强方法及其设备 | |
US20100179974A1 (en) | Signal Processing Method for Hierarchical Empirical Mode Decomposition and Apparatus Therefor | |
CN1806257A (zh) | 图像处理设备、图像处理方法、用于图像处理方法的程序、以及记录有用于图像处理方法的程序的记录介质 | |
US8594456B2 (en) | Image denoising method | |
CN1130069C (zh) | 假轮廓校正装置和方法 | |
RU2011143859A (ru) | Устройство для совершенствования изображений, способ совершенствования изображений, программа для совершенствования изображений и устройство для обработки сигналов | |
CN1135501C (zh) | 图像处理方法和装置 | |
CN1774030A (zh) | 图像噪声去除方法 | |
US8406559B2 (en) | Method and system for enhancing image sharpness based on local features of image | |
CN1291354C (zh) | 图象处理装置、图象处理方法、存储介质和程序 | |
US8885965B2 (en) | Image denoising method | |
CN108961172B (zh) | 一种基于Gamma校正的图像对比度增强方法 | |
US20040189874A1 (en) | Image detail enhancement system | |
TWI392361B (zh) | 濾波系統及實施影像之有效率濾波之方法 | |
CN1723691A (zh) | 降噪方法和设备 | |
EP1569168B1 (en) | Finite impulse response filter method and apparatus | |
CN104077746A (zh) | 灰度图像处理方法及其装置 | |
CN107481203A (zh) | 一种图像导向滤波方法及计算设备 | |
JP2008500757A (ja) | 映像信号の鮮鋭度を強調するための方法およびシステム | |
US10395344B2 (en) | Image processing method | |
CN113269681B (zh) | 一种基于多算法融合的水下实时图像增强方法 | |
CN115330637A (zh) | 一种图像锐化方法、装置、计算设备及存储介质 | |
Negi et al. | Survey of various image enhancement techniques in spatial domain using MatLab |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |