CN104217416B - 灰度图像处理方法及其装置 - Google Patents

灰度图像处理方法及其装置 Download PDF

Info

Publication number
CN104217416B
CN104217416B CN201310211284.4A CN201310211284A CN104217416B CN 104217416 B CN104217416 B CN 104217416B CN 201310211284 A CN201310211284 A CN 201310211284A CN 104217416 B CN104217416 B CN 104217416B
Authority
CN
China
Prior art keywords
pixel
value
edge pixel
edge
parameter
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.)
Expired - Fee Related
Application number
CN201310211284.4A
Other languages
English (en)
Other versions
CN104217416A (zh
Inventor
伍健荣
付轩
东明浩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to CN201310211284.4A priority Critical patent/CN104217416B/zh
Priority to JP2014106793A priority patent/JP6314654B2/ja
Priority to US14/290,169 priority patent/US9299163B2/en
Priority to EP14170582.2A priority patent/EP2819091B1/en
Publication of CN104217416A publication Critical patent/CN104217416A/zh
Application granted granted Critical
Publication of CN104217416B publication Critical patent/CN104217416B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/90—Determination of colour characteristics
    • 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
    • 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/20021—Dividing image into blocks, subimages or windows
    • 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/20048—Transform domain processing
    • G06T2207/20052—Discrete cosine transform [DCT]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Processing (AREA)
  • Facsimile Image Signal Circuits (AREA)
  • Image Analysis (AREA)

Abstract

本发明实施例提供一种灰度图像处理方法及其装置。该方法包括：计算原始灰度图像中边缘像素的均值和标准差，根据均值和标准差确定边缘像素的噪声方差；对边缘像素的第一邻近区域的图像进行二维DCT变换；根据该噪声方差，对二维DCT变换后的图像数据进行二维维纳滤波；接着进行二维DCT逆变换，获得滤波后的灰度图像；从滤波后的灰度图像中提取边缘像素的滤波后的像素值，将边缘像素的滤波后的像素值与该边缘像素的原始像素值进行加权求和，获得边缘像素的处理后的像素值。能够使得图像的边缘更加自然，并且使得图像的线边缘更加平直。

Description

灰度图像处理方法及其装置
技术领域
本发明涉及图像处理领域，特别涉及一种灰度图像处理方法及其装置。
背景技术
目前，随着数码相机、扫描仪等电子设备的流行与发展，很容易获得数码图像。然而，任何实体传感器（例如热传感器、电传感器或其他种类的传感器）都会受到一定程度噪声的影响，这些噪声将会影响信号测量的真实性，这样，获得的任何图像数据都包含了信号和噪声。各种各样与图像相关的应用，例如医学图像分析、图像分割、以及物体检测等，通常都需要有效的噪声抑制方法以进一步获得可靠的结果。因此，图像过滤已经成为图像处理以及计算机视觉中最为重要和广泛的研究课题之一。灰度图像是重要且广泛使用的图像类型之一，它的噪声抑制方法是非常重要的。但是，现有的滤波方法会使图像的边缘看起来不自然，并且，平直的边缘也可能看起来凹凸不平。
应该注意，上面对技术背景的介绍只是为了方便对本发明的技术方案进行清楚、完整的说明，并方便本领域技术人员的理解而阐述的。不能仅仅因为这些方案在本发明的背景技术部分进行了阐述而认为上述技术方案为本领域技术人员所公知。
发明内容
本发明实施例的目的在于提供一种灰度图像处理方法及其装置，能够使得图像的边缘更加自然，并且使得图像的线边缘更加平直。
根据本发明实施例的一个方面，提供了一种灰度图像处理方法，其中，所述方法包括：计算原始灰度图像中需要处理的每一个边缘像素的均值和标准差，根据所述均值和标准差确定每一个边缘像素的预先设定的噪声方差；对每一个边缘像素的第一邻近区域的原始灰度图像进行二维离散余弦变换；根据所述噪声方差，对二维离散余弦变换后的图像数据进行二维维纳滤波；对二维维纳滤波后的图像数据进行二维离散余弦逆变换，获得所述第一邻近区域的滤波后的灰度图像；从所述第一邻近区域的滤波后的灰度图像中提取每一个边缘像素的滤波后的像素值，将每一个边缘像素的滤波后的像素值与所述边缘像素的原始像素值进行加权求和，获得每一个边缘像素的处理后的像素值。
根据本发明实施例的另一个方面，提供了一种灰度图像处理装置，其中，所述装置包括：计算单元，所述计算单元用于计算原始灰度图像中需要处理的每一个边缘像素的均值和标准差；噪声方差确定单元，所述噪声方差确定单元用于根据所述均值和标准差确定每一个边缘像素的预先设定的噪声方差；二维离散余弦变换单元，所述二维离散余弦变换单元用于对每一个边缘像素的第一邻近区域的原始灰度图像进行二维离散余弦变换；滤波单元，所述滤波单元用于根据所述噪声方差，对二维离散余弦变换后的图像数据进行二维维纳滤波；二维离散余弦逆变换单元，所述二维离散余弦逆变换单元用于对二维维纳滤波后的图像数据进行二维离散余弦逆变换，获得所述第一邻近区域的滤波后的灰度图像；加权求和单元，所述加权求和单元用于从所述第一邻近区域的滤波后的灰度图像中提取每一个边缘像素的滤波后的像素值，将每一个边缘像素的滤波后的像素值与所述边缘像素的原始像素值进行加权求和，获得每一个边缘像素的处理后的像素值。
本发明实施例的有益效果在于：图像的边缘更加自然，并且，图像的线边缘更加平直。
参照后文的说明和附图，详细公开了本发明的特定实施方式，指明了本发明的原理可以被采用的方式。应该理解，本发明的实施方式在范围上并不因而受到限制。在所附权利要求的精神和条款的范围内，本发明的实施方式包括许多改变、修改和等同。
针对一种实施方式描述和/或示出的特征可以以相同或类似的方式在一个或更多个其它实施方式中使用，与其它实施方式中的特征相组合，或替代其它实施方式中的特征。
应该强调，术语“包括/包含”在本文使用时指特征、整件、步骤或组件的存在，但并不排除一个或更多个其它特征、整件、步骤或组件的存在或附加。
附图说明
参照以下的附图可以更好地理解本发明的很多方面。附图中的部件不是成比例绘制的，而只是为了示出本发明的原理。为了便于示出和描述本发明的一些部分，附图中对应部分可能被放大或缩小。在本发明的一个附图或一种实施方式中描述的元素和特征可以与一个或更多个其它附图或实施方式中示出的元素和特征相结合。此外，在附图中，类似的标号表示几个附图中对应的部件，并可用于指示多于一种实施方式中使用的对应部件。在附图中：
图1是本发明实施例1的灰度图像处理方法的流程图；
图2是本发明实施例1的根据边缘像素的均值和标准差确定噪声方差的方法的流程图；
图3是本发明实施例2的灰度图像处理装置的结构示意图；
图4是本发明实施例2的噪声方差确定单元的结构示意图；
图5是本发明实施例3的灰度图像处理方法的流程图；
具体实施方式
参照附图，通过下面的说明书，本发明的前述以及其它特征将变得明显。在说明书和附图中，具体公开了本发明的实施方式，其表明了其中可以采用本发明的原则的部分实施方式，应了解的是，本发明不限于所描述的实施方式，相反，本发明包括落入所附权利要求的范围内的全部修改、变型以及等同物。
目前，现有的滤波方法会使图像的边缘看起来不自然，并且，平直的边缘也可能看起来凹凸不平。本申请的发明人在实现本发明的过程中发现，可以通过检测边缘像素并进行平滑处理的方法使得图像边缘看起来更自然。
对于边缘像素的检测，现有的方法例如有索贝尔法（Sobel method）、普瑞维特法（Prewitt method）、以及梯度法（Gradient method）等。其中，
索贝尔法基于水平和垂直方向上邻近像素的权值差异来检测边缘像素，其中，边缘像素的权值差异比其他像素的权值差异大得多。然而，该检测结果并不精确，并且，检测结果的范围太小，例如，仅为一个像素宽度。
普瑞维特法基于水平和垂直方向上邻近像素的强度差异来检测边缘像素，其中，边缘像素的强度差异比其他像素的权值差异大得多。但是，该方法存在与上述索贝尔法相同的问题。
梯度法基于梯度的强度值检测边缘像素，其中，边缘像素的强度值比其他像素的强度值大的多。但是，该方法对于图像噪声过于敏感。
对于平滑处理，现有的方法例如有高斯滤波法（Gaussian filter method）、非局部平均法（non-local means method）、以及中值滤波法（median filter method）等。其中，
高斯滤波法利用高斯滤波来平滑图像，计算量较小且效率较高。但是，该方法使得图像模式且不自然。
非局部平均法利用非局部平均滤波来平滑图像，滤波结果比高斯滤波法的滤波结果更加自然，但是，图像边缘的凹凸不平感并没有被有效抑制。
中值滤波法与高斯滤波法相似，不同之处在于该方法将当前像素的强度替换为邻近区域的中值。但是，该方法存在的问题与高斯滤波法相同。
因此，现有的方法并不能使得图像的边缘自然，且使得图像的线边缘平直。
以下结合附图对本发明的非线性补偿装置及其方法进行详细说明。
实施例1
图1是本发明实施例1的灰度图像处理方法的流程图。如图1所示，该方法包括：
步骤101：计算原始灰度图像中需要处理的每一个边缘像素的均值和标准差，根据该均值和标准差确定每一个边缘像素的预先设定的噪声方差；
步骤102：对每一个边缘像素的第一邻近区域的原始灰度图像进行二维离散余弦变换（2D-Discrete Cosine Transform，2D-DCT）；
步骤103：根据该噪声方差，对二维离散余弦变换后的图像数据进行二维维纳滤波；
步骤104：对二维维纳滤波后的图像数据进行二维离散余弦逆变换，获得该第一邻近区域的滤波后的灰度图像；
步骤105：从该第一邻近区域的滤波后的灰度图像中提取每一个边缘像素的滤波后的像素值，将每一个边缘像素的滤波后的像素值与该边缘像素的原始像素值进行加权求和，获得每一个边缘像素的处理后的像素值。
由上述实施例可知，通过根据边缘像素的均值和标准差确定噪声方差，并且在二维DCT域中基于该噪声方差进行二维维纳滤波，能够使得图像的边缘更加自然，并且使得图像的线边缘更加平直。
在步骤101中，对于如何确定需要处理的边缘像素，可使用现有技术中的任一种方法。在本发明中，例如可将原始灰度图像中的像素的差值和标准差分别与预先设定的第一阈值和第二阈值进行比较，当该像素的差值大于该第一阈值且标准差大于该第二阈值时，则将该像素作为需要处理的边缘像素。
其中，对于计算像素的差值的方法，可使用现有技术中的任一种方法。在本发明中，可利用下述公式（1）计算像素(i,j)的差值：
diff(i,j)=max(|x(i,j)-x(k,m)|)，k∈[i-Rd,i+Rd]，m∈[j-Rd,j+Rd] （1）
其中，diff(i,j)表示像素(i,j)的差值，x(i,j)表示像素(i,j)的像素值，x(k,m)表示像素(k,m)的像素值，像素(k,m)是像素(i,j)的第二邻近区域内的像素，Rd表示该第二邻近区域的半径，Rd为大于等于1的整数，i和k为1～A的整数，j和m为1～B的整数，A和B分别是原始灰度图像的像素的行数和列数。
例如，对于像素(2,2)，可取Rd=1。那么，该第二邻近区域可表示为：
此时，像素(2,2)的差值可利用下述公式（2）计算：
diffex=max(|x(2,2)-x(k,m)|) k,m=1,2,3 （2）
其中，diffex表示像素(2,2)的差值，x(2,2)表示像素(2,2)的像素值，x(k,m)表示像素(k,m)的像素值，像素(k,m)是像素(2,2)的该第二邻近区域内的像素。
在本实施例中，对于计算像素的均值和标准差的方法，可使用现有技术中的任一种方法。在本发明中，可利用下述公式（3）和（4）计算像素(i,j)的均值和标准差：
其中，mean(i,j)表示像素(i,j)的均值，std(i,j)表示像素(i,j)的标准差，x(k,m)表示像素(k,m)的像素值，像素(k,m)是像素(i,j)的第三邻近区域内的像素，Rsm表示该第三邻近区域的半径，Rsm为大于等于1的整数，i和k为1～A的整数，j和m为1～B的整数，A和B分别是原始灰度图像的像素的行数和列数。
例如，对于像素(5,5)，可取Rsm=4。那么，该第三邻近区域可表示为：
此时，像素(5,5)的均值和标准差可利用下述公式（5）和（6）计算：
其中，meanex表示像素(5,5)的均值，stdex表示像素(5,5)的标准差，x(k,m)表示像素(k,m)的像素值，像素(k,m)是像素(5,5)的第三邻近区域内的像素。
这样，通过以上方法计算能够获得像素(i,j)的差值、均值和标准差。
在本实施例中，可根据实际需要预先设定第一阈值和第二阈值，当像素(i,j)的差值大于该第一阈值，并且像素(i,j)的标准差大于该第二阈值时，确定该像素(i,j)为边缘像素，需要对其进行处理。其中，该第一阈值例如可根据实际需要选择小于25的某个数值，该第二阈值例如可根据实际需要选择范围在10～30的某个数值。
在本实施例中，可基于实际需要，根据边缘像素的均值和标准差确定噪声方差，图2是本实施例的根据边缘像素的均值和标准差确定噪声方差的方法的流程图，但本发明不限于这种方法。如图2所示，该方法包括：
步骤201：当该标准差大于预先设定的第二阈值且该均值大于预先设定的第三阈值时，将该噪声方差设为第一参数；
步骤202：当该标准差大于该第二阈值且小于等于预先设定的第四阈值时，将该噪声方差设为第二参数；
步骤203：当该标准差大于该第四阈值时，将该噪声方差设为第三参数；
其中，第一参数、第二参数和第三参数满足以下关系：
第一参数＜第三参数＜第二参数。
其中，第二阈值的范围同上所述，第三阈值例如可根据实际需要选择范围在100以上的某个数值，第四阈值大于第二阈值，该第四阈值例如可根据实际需要选择范围在50～100的某个数值。
其中，该第一参数、第二参数和第三参数的选择，只要满足第一参数＜第三参数＜第二参数的关系即可，并不严格的限制。例如，第一参数可根据实际需要选择范围在10～50的某个数值，第二参数可根据实际需要选择范围在50～200的某个数值，第三参数可根据实际需要选择范围在30～100的某个数值。
在步骤102中，对每一个边缘像素的第一邻近区域的原始灰度图像进行二维离散余弦变换（2D-Discrete Cosine Transform，2D-DCT）。其中，对图像进行二维离散余弦变换属于现有技术，例如，本发明中对于边缘像素的第一邻近区域的原始灰度图像进行二维离散余弦变换可使用下述公式（7）和（8）进行变换：
其中，D(u,v)表示原始灰度图像中边缘像素的第一邻近区域内的像素(x,y)经过二维离散余弦变换后的图像数据，I(x,y)表示像素(x,y)的像素值，a(0),a(u),a(v)为变换系数，x,u=0,1,...M-1,y,v=0,1,...N-1，M和N表示该第一邻近区域的长度方向和宽度方向的像素个数，M和N均为大于1的整数，优选为大于等于9的整数。
通过上述公式（7）和（8），获得了边缘像素(x,y)的第一邻近区域的经过二维离散余弦变换后的图像数据。
在步骤103中，根据步骤101中确定的噪声方差，对二维离散余弦变换后的图像数据进行二维维纳滤波。
在本实施例中，可使用现有的二维维纳滤波方法对二维离散余弦变换后的图像数据进行二维维纳滤波。
例如，可使用下述公式（9）至（11）实现二维维滤波：
其中，mean(i,j)表示该第一邻近区域内的像素(i,j)的均值，x(k,m)表示像素(k,m)的像素值，像素(k,m)是像素(i,j)的第四邻近区域内的像素，var(i,j)表示像素(i,j)的方差，Wn(i,j)表示像素(i,j)经过二维维纳滤波后的图像数据，nv表示噪声方差，M'和N'表示该第四邻近区域的长度方向和宽度方向的像素个数，M'和N'均为大于1的整数。
对于二维离散余弦变换后的图像数据进行二维维纳滤波，即在DCT域中进行二维维纳滤波可在上述公式（11）的基础上，利用下述公式（12）进行：
其中，Wn(u,v)表示在DCT域中对应于像素(i,j)的图像数据经过二维维纳滤波后的图像数据，d(u,v)表示DCT域中对应于像素(i,j)的图像数据。
其中，噪声方差nv已经在上述步骤101中确定，根据该确定的nv，在DCT域中对于边缘像素(x,y)的第一邻近区域中的图像数据进行二维维纳变换，获得二维维纳滤波后的图像数据。
在步骤104中，对二维维纳滤波后的图像数据进行二维离散余弦逆变换，获得该边缘像素的第一邻近区域的滤波后的灰度图像。
其中，二维离散余弦逆变换可使用现有技术，例如，与上述公式（7）和（8）相对应，可使用下述公式（13）和（14）进行二维离散余弦逆变换：
其中，ID(x,y)表示边缘像素的第一邻近区域内像素(x,y)经过滤波后的像素值，D(u,v)表示表示原始灰度图像中边缘像素的第一邻近区域内的像素(x,y)经过二维离散余弦变换后的图像数据，a(0),a(u),a(v)为变换系数，x,u=0,1,...M-1,y,v=0,1,...N-1，M和N表示该第一邻近区域的长度方向和宽度方向的像素个数，M和N均为大于1的整数，优选为大于等于9的整数。
这样，就获得了获得一个边缘像素的该第一邻近区域的所有像素的滤波后的像素值，从而获得该边缘像素的第一邻近区域的滤波后的灰度图像。
在步骤105中，从一个边缘像素的该第一邻近区域的滤波后的灰度图像中提取该边缘像素的滤波后的像素值，将该边缘像素的滤波后的像素值与该边缘像素的原始像素值进行加权求和，获得该边缘像素的处理后的像素值。对于每一个边缘像素，重复以上步骤，则获得每一个边缘像素的处理后的像素值，从而获得处理后的灰度图像。
其中，可利用下述公式（15）进行像素值的加权求和：
ProImg=a*DCT-WN+(1-a)*OriImg （15）
其中，ProImg表示该边缘像素处理后的像素值，DCT-WN表示该边缘像素的滤波后的像素值，OriImg表示该边缘像素的原始像素值，a为预先设定的加权系数，0＜a≤1，其数值可根据实际需要而定。
由上述实施例可知，通过根据边缘像素的均值和标准差确定噪声方差，并且在二维DCT域中基于该噪声方差进行二维维纳滤波，能够使得图像的边缘更加自然，并且使得图像的线边缘更加平直。
实施例2
图3是本发明实施例2的灰度图像处理装置300的结构示意图，对应于实施例1的灰度图像处理方法。如图3所示，该装置300包括：计算单元301、噪声方差确定单元302、二维离散余弦变换单元303、滤波单元304、二维离散余弦逆变换单元305以及加权求和单元306。其中，
计算单元301用于计算原始灰度图像中需要处理的每一个边缘像素的均值和标准差；
噪声方差确定单元302用于根据该均值和标准差确定每一个边缘像素的预先设定的噪声方差；
二维离散余弦变换单元303用于对每一个边缘像素的第一邻近区域的原始灰度图像进行二维离散余弦变换；
滤波单元304用于根据该噪声方差，对二维离散余弦变换后的图像数据进行二维维纳滤波；
二维离散余弦逆变换单元305用于对二维维纳滤波后的图像数据进行二维离散余弦逆变换，获得该第一邻近区域的滤波后的灰度图像；
加权求和单元306用于从该第一邻近区域的滤波后的灰度图像中提取每一个边缘像素的滤波后的像素值，将每一个边缘像素的滤波后的像素值与该边缘像素的原始像素值进行加权求和，获得每一个边缘像素的处理后的像素值。
在本实施例中，该灰度图像处理装置300还可包括：
边缘像素确定单元307，该边缘像素确定单元307用于确定需要处理的边缘像素。
在本实施例中，该边缘像素确定单元307为可选部件，在图3中用虚线表示。
在本实施例中，噪声方差确定单元302用于根据该均值和标准差确定每一个边缘像素的预先设定的噪声方差。图4是本实施例的噪声方差确定单元302的结构示意图，但本发明不限于这种结构。如图4所示，噪声方差确定单元302包括判断单元401和确定单元402，其中，
判断单元401用于判断是否该标准差大于预先设定的第二阈值且该均值大于预先设定的第三阈值、或者是否该标准差大于该第二阈值且小于等于预先设定的第四阈值、或者是否该标准差大于该第四阈值；
确定单元402用于在该标准差大于预先设定的第二阈值且该均值大于预先设定的第三阈值时，将该噪声方差设为第一参数；在所述标准差大于该第二阈值且小于等于预先设定的第四阈值时，将该噪声方差设为第二参数；在该标准差大于所述第四阈值时，将该噪声方差设为第三参数；
其中，第一参数、第二参数和第三参数满足以下关系：
第一参数＜第三参数＜第二参数。
在本实施例中，确定需要处理的边缘像素的方法、计算边缘像素的均值和标准差的方法、根据该均值和标准差确定噪声方差的方法、进行二维离散余弦变换及逆变换、在二维DCT域进行二维维纳滤波的方法以及对像素值进行加权求和的方法均与实施例1中的记载相同，此处不再重复。
由上述实施例可知，通过根据边缘像素的均值和标准差确定噪声方差，并且在二维DCT域中基于该噪声方差进行二维维纳滤波，能够使得图像的边缘更加自然，并且使得图像的线边缘更加平直。
实施例3
图5是本发明实施例3的灰度图像处理方法的流程图。如图5所示，该方法包括：
步骤501：计算原始灰度图像中的一个像素的差值diff、标准差std和均值mean；
其中，计算差值、标准差和均值的方法与实施例1中的记载相同，此处不再重复。
步骤502：将该像素的差值diff和标准差std分别与预先设定的第一阈值Thre_d和第二阈值Thre_sl进行比较，当该差值diff大于该第一阈值Thre_d且标准差std大于第二阈值Thre_sl时，将该像素确定为需要处理的边缘像素，进入步骤503；否则，当前像素处理结束；
其中，第一阈值和第二阈值的设定方法与实施例1中的记载相同，此处不再重复。
步骤503：将该边缘像素的标准差和均值分别与预先设定的该第二阈值Thre_sl以及第三阈值Thre_m进行比较，当该标准差大于该第二阈值，并且，该均值大于该第三阈值时，进入步骤504；否则，进入步骤505；
其中，第三阈值的设定方法与实施例1中的记载相同，此处不再重复。
步骤504：将噪声方差设为第一参数，进入步骤509；
步骤505：将该边缘像素的标准差分别与预先设定的该第二阈值以及第四阈值进行比较，当该标准差大于该第二阈值Thre_sl且小于等于该第四阈值Thre_sh时，进入步骤506，否则，进入步骤507；
其中，第四阈值的设定方法与实施例1中的记载相同，此处不再重复。
步骤506：将噪声方差设为第二参数，进入步骤509；
步骤507：当该边缘像素的标准差大于预先设定的该第四阈值时，进入步骤508；否则，当前像素处理结束；
步骤508：将噪声方差设为第三参数，进入步骤509；
其中，步骤504、506、508中的第一参数、第二参数和第三参数满足以下关系：
第一参数＜第三参数＜第二参数。
其中，第一参数、第二参数和第三参数的设定方法与实施例1中的记载相同，此处不再重复。
步骤509：对该边缘像素的第一邻近区域的原始灰度图像进行二维离散余弦变换；
其中，进行二维离散余弦变换的方法与实施例1中的记载相同，此处不再重复。
步骤510：根据该噪声方差，对二维离散余弦变换后的图像数据进行二维维纳滤波；
其中，在DCT域中进行二维维纳滤波的方法与实施例1中的记载相同，此处不再重复。
步骤511：对二维维纳滤波后的图像数据进行二维离散余弦逆变换，获得该边缘的第一邻近区域的滤波后的灰度图像；
其中，进行二维离散余弦逆变换的方法与实施例1中的记载相同，此处不再重复。
步骤512：从该边缘的第一邻近区域的滤波后的灰度图像中提取该边缘像素的滤波后的像素值，将该边缘像素的滤波后的像素值与该边缘像素的原始像素值进行加权求和，获得该边缘像素的处理后的像素值。
其中，像素值加权求和的方法与实施例1中的记载相同，此处不再重复。
重复上述步骤501至512，即能够获得所有需要处理的边缘像素的处理后的像素值，从而获得处理后的灰度图像。
由上述实施例可知，通过根据边缘像素的均值和标准差确定噪声方差，并且在二维DCT域中基于该噪声方差进行二维维纳滤波，能够使得图像的边缘更加自然，并且使得图像的线边缘更加平直。
本发明以上的装置和方法可以由硬件实现，也可以由硬件结合软件实现。本发明涉及这样的计算机可读程序，当该程序被逻辑部件所执行时，能够使该逻辑部件实现上文该的装置或构成部件，或使该逻辑部件实现上文该的各种方法或步骤。
本发明还涉及用于存储以上程序的存储介质，如硬盘、磁盘、光盘、DVD、flash存储器等。
以上结合具体的实施方式对本发明进行了描述，但本领域技术人员应该清楚，这些描述都是示例性的，并不是对本发明保护范围的限制。本领域技术人员可以根据本发明的精神和原理对本发明做出各种变型和修改，这些变型和修改也在本发明的范围内。
附记
附记1、一种灰度图像的处理方法，其中，所述方法包括：
计算原始灰度图像中需要处理的每一个边缘像素的均值和标准差，根据所述均值和标准差确定每一个边缘像素的预先设定的噪声方差；
对每一个边缘像素的第一邻近区域的原始灰度图像进行二维离散余弦变换；
根据所述噪声方差，对二维离散余弦变换后的图像数据进行二维维纳滤波；
对二维维纳滤波后的图像数据进行二维离散余弦逆变换，获得所述第一邻近区域的滤波后的灰度图像；
从所述第一邻近区域的滤波后的灰度图像中提取每一个边缘像素的滤波后的像素值，将每一个边缘像素的滤波后的像素值与所述边缘像素的原始像素值进行加权求和，获得每一个边缘像素的处理后的像素值。
附记2、根据附记1所述的方法，其中，所述方法还包括：确定需要处理的边缘像素；
其中，所述确定需要处理的边缘像素包括：将原始灰度图像中的像素的差值大于预先设定的第一阈值且标准差大于预先设定的第二阈值的像素作为需要处理的边缘像素。
附记3、根据附记2所述的方法，其中，
利用以下的公式（1）计算所述像素的差值：
diff(i,j)=max(|x(i,j)-x(k,m)|)，k∈[i-Rd,i+Rd]，m∈[j-Rd,j+Rd] (1)
其中，diff(i,j)表示像素(i,j)的差值，x(i,j)表示像素(i,j)的像素值，x(k,m)表示像素(k,m)的像素值，像素(k,m)是像素(i,j)的第二邻近区域内的像素，Rd表示该第二邻近区域的半径，i和k为1～A的整数，j和m为1～B的整数，A和B分别是原始灰度图像的像素的行数和列数。
附记4、根据附记1所述的方法，其中，所述根据均值和标准差选择预先设定的噪声方差包括：
当所述标准差大于预先设定的第二阈值且所述均值大于预先设定的第三阈值时，将所述噪声方差设为第一参数；
当所述标准差大于所述第二阈值且小于等于预先设定的第四阈值时，将所述噪声方差设为第二参数；
当所述标准差大于所述第四阈值时，将所述噪声方差设为第三参数；
其中，所述第一参数、第二参数和第三参数满足以下关系：
所述第一参数＜所述第三参数＜所述第二参数。
附记5、根据附记1所述的方法，其中，所述将每一个边缘像素的滤波后的像素值与所述边缘像素的原始像素值进行加权求和，获得每一个边缘像素的处理后的像素值包括：
利用以下的公式（2）进行所述加权求和：
ProImg=a*DCT-WN+(1-a)*OriImg （2）
其中，ProImg表示所述边缘像素处理后的像素值，DCT-WN表示所述边缘像素的滤波后的像素值，OriImg表示所述边缘像素的原始像素值，a为预先设定的加权系数，0＜a≤1。
附记6、根据附记1-5中的任一项附记所述的方法，其中，所述第一邻近区域是每一个边缘像素周围的尺寸为M×N的区域，M表示所述第一邻近区域高度方向的像素个数，N表示所述第一邻近区域宽度方向的像素个数，M和N均为大于等于9的整数。
附记7、一种灰度图像的处理装置，其中，所述装置包括：
计算单元，所述计算单元用于计算原始灰度图像中需要处理的每一个边缘像素的均值和标准差；
噪声方差确定单元，所述噪声方差确定单元用于根据所述均值和标准差确定每一个边缘像素的预先设定的噪声方差；
二维离散余弦变换单元，所述二维离散余弦变换单元用于对每一个边缘像素的第一邻近区域的原始灰度图像进行二维离散余弦变换；
滤波单元，所述滤波单元用于根据所述噪声方差，对二维离散余弦变换后的图像数据进行二维维纳滤波；
二维离散余弦逆变换单元，所述二维离散余弦逆变换单元用于对二维维纳滤波后的图像数据进行二维离散余弦逆变换，获得所述第一邻近区域的滤波后的灰度图像；
加权求和单元，所述加权求和单元用于从所述第一邻近区域的滤波后的灰度图像中提取每一个边缘像素的滤波后的像素值，将每一个边缘像素的滤波后的像素值与所述边缘像素的原始像素值进行加权求和，获得每一个边缘像素的处理后的像素值。
附记8、根据附记7所述的装置，其中，所述装置还包括：
边缘像素确定单元，所述边缘像素确定单元用于确定需要处理的边缘像素；
其中，所述边缘像素确定单元将原始灰度图像中的像素的差值和标准差分别大于预先设定的第一阈值和第二阈值的像素作为需要处理的边缘像素。
附记9、根据附记8所述的装置，其中，
所述计算单元利用以下的公式（1）计算所述像素的差值：
diff(i,j)=max(|x(i,j)-x(k,m)|)，k∈[i-Rd,i+Rd]，m∈[j-Rd,j+Rd] (1)
其中，diff(i,j)表示像素(i,j)的差值，x(i,j)表示像素(i,j)的像素值，x(k,m)表示像素(k,m)的像素值，像素(k,m)是像素(i,j)的第二邻近区域内的像素，Rd表示该第二邻近区域的半径，i和k为1～A的整数，j和m为1～B的整数，A和B分别是原始灰度图像的像素的行数和列数。
附记10、根据附记7所述的装置，其中，所述噪声方差确定单元包括判断单元和确定单元，其中，
所述判断单元用于判断是否所述标准差大于预先设定的第二阈值且所述均值大于预先设定的第三阈值、或者是否所述标准差大于所述第二阈值且小于等于预先设定的第四阈值、或者是否所述标准差大于所述第四阈值；
所述确定单元用于在所述标准差大于预先设定的第二阈值且所述均值大于预先设定的第三阈值时，将所述噪声方差设为第一参数；在所述标准差大于所述第二阈值且小于等于预先设定的第四阈值时，将所述噪声方差设为第二参数；在所述标准差大于所述第四阈值时，将所述噪声方差设为第三参数；
其中，所述第一参数、第二参数和第三参数满足以下关系：
所述第一参数＜所述第三参数＜所述第二参数。
附记11、根据附记7所述的装置，其中，所述加权求和单元利用以下的公式（2）进行所述加权求和：
ProImg=a*DCT-WN+(1-a)*OriImg （2）
其中，ProImg表示所述边缘像素处理后的像素值，DCT-WN表示所述边缘像素的滤波后的像素值，OriImg表示所述边缘像素的原始像素值，a为预先设定的加权系数，0＜a≤1。
附记12、根据附记7-11中的任一项附记所述的装置，其中，所述第一邻近区域是每一个边缘像素周围的尺寸为M×N的区域，M表示所述第一邻近区域高度方向的像素个数，N表示所述第一邻近区域宽度方向的像素个数，M和N均为大于等于9的整数。

Claims (8)

1.一种灰度图像的处理方法，其中，所述方法包括：
计算原始灰度图像中需要处理的每一个边缘像素的均值和标准差，根据所述均值和标准差确定每一个边缘像素的预先设定的噪声方差；
对每一个边缘像素的第一邻近区域的原始灰度图像进行二维离散余弦变换；
根据所述噪声方差，对二维离散余弦变换后的图像数据进行二维维纳滤波；
对二维维纳滤波后的图像数据进行二维离散余弦逆变换，获得所述第一邻近区域的滤波后的灰度图像；
从所述第一邻近区域的滤波后的灰度图像中提取每一个边缘像素的滤波后的像素值，将每一个边缘像素的滤波后的像素值与所述边缘像素的原始像素值进行加权求和，获得每一个边缘像素的处理后的像素值；
其中，根据均值和标准差确定所述噪声方差包括：
当所述标准差大于预先设定的第二阈值且所述均值大于预先设定的第三阈值时，将所述噪声方差设为第一参数；
当所述标准差大于所述第二阈值且小于等于预先设定的第四阈值时，将所述噪声方差设为第二参数；
当所述标准差大于所述第四阈值时，将所述噪声方差设为第三参数；
其中，所述第一参数、第二参数和第三参数满足以下关系：
所述第一参数＜所述第三参数＜所述第二参数。
2.根据权利要求1所述的方法，其中，所述方法还包括：确定需要处理的边缘像素；
其中，所述确定需要处理的边缘像素包括：将原始灰度图像中的像素的差值大于预先设定的第一阈值且标准差大于预先设定的第二阈值的像素作为需要处理的边缘像素。
3.根据权利要求2所述的方法，其中，
利用以下的公式(1)计算所述像素的差值：
diff(i,j)＝max(|x(i,j)-x(k,m)|)，k∈[i-Rd,i+Rd]，m∈[j-Rd,j+Rd] (1)
其中，diff(i,j)表示像素(i,j)的差值，x(i,j)表示像素(i,j)的像素值，x(k,m)表示像素(k,m)的像素值，像素(k,m)是像素(i,j)的第二邻近区域内的像素，Rd表示该第二邻近区域的半径，i和k为1～A的整数，j和m为1～B的整数，A和B分别是原始灰度图像的像素的行数和列数。
4.根据权利要求1所述的方法，其中，所述将每一个边缘像素的滤波后的像素值与所述边缘像素的原始像素值进行加权求和，获得每一个边缘像素的处理后的像素值包括：
利用以下的公式(2)进行所述加权求和：
ProImg＝a*DCT-WN+(1-a)*OriImg (2)
其中，ProImg表示所述边缘像素处理后的像素值，DCT-WN表示所述边缘像素的滤波后的像素值，OriImg表示所述边缘像素的原始像素值，a为预先设定的加权系数，0＜a≤1。
5.根据权利要求1-4中的任一项所述的方法，其中，所述第一邻近区域是每一个边缘像素周围的尺寸为M×N的区域，M表示所述第一邻近区域高度方向的像素个数，N表示所述第一邻近区域宽度方向的像素个数，M和N均为大于等于9的整数。
6.一种灰度图像的处理装置，其中，所述装置包括：
计算单元，所述计算单元用于计算原始灰度图像中需要处理的每一个边缘像素的均值和标准差；
噪声方差确定单元，所述噪声方差确定单元用于根据所述均值和标准差确定每一个边缘像素的预先设定的噪声方差；
二维离散余弦变换单元，所述二维离散余弦变换单元用于对每一个边缘像素的第一邻近区域的原始灰度图像进行二维离散余弦变换；
滤波单元，所述滤波单元用于根据所述噪声方差，对二维离散余弦变换后的图像数据进行二维维纳滤波；
二维离散余弦逆变换单元，所述二维离散余弦逆变换单元用于对二维维纳滤波后的图像数据进行二维离散余弦逆变换，获得所述第一邻近区域的滤波后的灰度图像；
加权求和单元，所述加权求和单元用于从所述第一邻近区域的滤波后的灰度图像中提取每一个边缘像素的滤波后的像素值，将每一个边缘像素的滤波后的像素值与所述边缘像素的原始像素值进行加权求和，获得每一个边缘像素的处理后的像素值；
其中，所述噪声方差确定单元包括判断单元和确定单元，所述判断单元包括第一判断单元、第二判断单元和第三判断单元，所述确定单元包括第一确定单元、第二确定单元和第三确定单元，其中，
所述第一判断单元用于判断是否所述标准差大于预先设定的第二阈值且所述均值大于预先设定的第三阈值；
所述第一确定单元用于在所述标准差大于预先设定的第二阈值且所述均值大于预先设定的第三阈值时，将所述噪声方差设为第一参数；
所述第二判断单元用于判断是否所述标准差大于所述第二阈值且小于等于预先设定的第四阈值；
所述第二确定单元用于在所述标准差大于所述第二阈值且小于等于预先设定的第四阈值时，将所述噪声方差设为第二参数；
所述第三判断单元用于判断是否所述标准差大于所述第四阈值；
所述第三确定单元用于在所述标准差大于所述第四阈值时，将所述噪声方差设为第三参数；
其中，所述第一参数、第二参数和第三参数满足以下关系：
所述第一参数＜所述第三参数＜所述第二参数。
7.根据权利要求6所述的装置，其中，所述装置还包括：
边缘像素确定单元，所述边缘像素确定单元用于确定需要处理的边缘像素；
其中，所述边缘像素确定单元将原始灰度图像中的像素的差值和标准差分别大于预先设定的第一阈值和第二阈值的像素作为需要处理的边缘像素。
8.根据权利要求6所述的装置，其中，所述加权求和单元利用以下的公式(2)进行所述加权求和：
ProImg＝a*DCT-WN+(1-a)*OriImg (2)
其中，ProImg表示所述边缘像素处理后的像素值，DCT-WN表示所述边缘像素的滤波后的像素值，OriImg表示所述边缘像素的原始像素值，a为预先设定的加权系数，0＜a≤1。
CN201310211284.4A 2013-05-31 2013-05-31 灰度图像处理方法及其装置 Expired - Fee Related CN104217416B (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201310211284.4A CN104217416B (zh) 2013-05-31 2013-05-31 灰度图像处理方法及其装置
JP2014106793A JP6314654B2 (ja) 2013-05-31 2014-05-23 グレースケール画像処理方法及び装置
US14/290,169 US9299163B2 (en) 2013-05-31 2014-05-29 Method and apparatus for processing edges of a gray image
EP14170582.2A EP2819091B1 (en) 2013-05-31 2014-05-30 Method and apparatus for processing a gray image

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310211284.4A CN104217416B (zh) 2013-05-31 2013-05-31 灰度图像处理方法及其装置

Publications (2)

Publication Number Publication Date
CN104217416A CN104217416A (zh) 2014-12-17
CN104217416B true CN104217416B (zh) 2017-09-15

Family

ID=50841642

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310211284.4A Expired - Fee Related CN104217416B (zh) 2013-05-31 2013-05-31 灰度图像处理方法及其装置

Country Status (4)

Country Link
US (1) US9299163B2 (zh)
EP (1) EP2819091B1 (zh)
JP (1) JP6314654B2 (zh)
CN (1) CN104217416B (zh)

Families Citing this family (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104792263B (zh) 2015-04-20 2018-01-05 合肥京东方光电科技有限公司 确定显示母板的待检测区域的方法和装置
JP6729227B2 (ja) * 2016-09-14 2020-07-22 コニカミノルタ株式会社 画像処理装置、画像形成装置及びプログラム
CN106856001A (zh) * 2016-12-08 2017-06-16 彭志勇 灰度图像的保存方法和灰度图像原始像素的获取方法
CN107271046B (zh) * 2017-06-12 2019-08-23 中国科学院长春光学精密机械与物理研究所 一种基于局部统计相关性的非均匀性校正方法及装置
CN109523483A (zh) * 2018-11-14 2019-03-26 北京奇艺世纪科技有限公司 一种图像去雾方法及装置
CN109741267B (zh) * 2018-12-05 2023-04-25 西安电子科技大学 基于三边滤波和神经网络的红外图像非均匀性校正方法
CN109934262B (zh) * 2019-01-31 2023-08-22 平安科技（深圳）有限公司 图片差异性判断方法、装置、计算机设备和存储介质
CN110335198B (zh) * 2019-07-08 2021-06-08 威创集团股份有限公司 一种图像处理方法及系统
CN110390649B (zh) * 2019-07-16 2023-03-24 西安石油大学 一种用于油气管道焊缝图像降噪的方法
CN111179295B (zh) * 2019-12-23 2023-04-04 青海大学 改进的二维Otsu阈值图像分割方法和系统
CN111833341B (zh) * 2020-07-22 2024-08-13 浙江大华技术股份有限公司 确定图像中条纹噪声的方法及装置
CN112017109B (zh) * 2020-08-24 2023-03-21 内蒙古农业大学 一种在线铁谱视频图像气泡的消除方法
CN112487929B (zh) * 2020-11-25 2024-10-15 深圳市云希谷科技有限公司 儿童绘本的图像识别方法、装置、设备及存储介质
CN114065798A (zh) * 2021-02-23 2022-02-18 杭州博工科技有限公司 基于机器识别的视觉识别方法及装置
CN113129235B (zh) * 2021-04-22 2024-10-08 深圳市深图医学影像设备有限公司 一种医学图像噪声抑制算法
CN113902627A (zh) * 2021-08-24 2022-01-07 北京临近空间飞行器系统工程研究所 一种面向rd时频数据的数据增强系统
CN114399688A (zh) * 2021-12-15 2022-04-26 中国资源卫星应用中心 一种遥感图像刃边区域的选取方法及装置
CN114663451B (zh) * 2022-02-28 2024-10-18 惠州学院 一种可处理噪声和纹理的边缘检测方法和系统
CN114742727B (zh) * 2022-03-31 2023-05-05 南通电博士自动化设备有限公司 一种基于图像平滑的噪声处理方法及系统
WO2024031579A1 (zh) * 2022-08-11 2024-02-15 上海玄戒技术有限公司 图像降噪方法、装置及芯片
CN117710220A (zh) * 2022-08-31 2024-03-15 北京字跳网络技术有限公司 图像处理方法、装置、设备及介质
CN115870815B (zh) * 2022-12-28 2023-06-16 中山超精科技有限公司 镜片抛光检测方法及装置
CN115690430B (zh) * 2022-12-29 2023-03-28 尚特杰电力科技有限公司 光伏组件边缘特征图补齐方法、系统、设备及存储介质
CN116309579B (zh) * 2023-05-19 2023-08-01 惠州市宝惠电子科技有限公司 一种用图像处理的变压器焊缝质量检测方法
CN118014882B (zh) * 2024-04-08 2024-06-04 汉中精测电器有限责任公司 一种基于图像处理的扭矩传感器的钻孔定位方法
CN118096754B (zh) * 2024-04-26 2024-07-30 广州拓威天海国际物流有限公司 一种国际货运服务标准价格监测方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5729631A (en) * 1993-11-30 1998-03-17 Polaroid Corporation Image noise reduction system using a wiener variant filter in a pyramid image representation
US7590302B1 (en) * 2004-11-01 2009-09-15 Pixelworks, Inc. Image edge enhancement system and method
CN102006484A (zh) * 2009-08-31 2011-04-06 索尼公司 图像处理装置、图像处理方法以及程序

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10136258A (ja) * 1996-10-25 1998-05-22 Sony Corp 輪郭抽出装置及び輪郭抽出方法
KR100497398B1 (ko) * 2003-07-16 2005-06-23 삼성전자주식회사 영상 잡음 제거방법 및 장치
ITVA20060060A1 (it) * 2006-10-06 2008-04-07 St Microelectronics R&D Ltd Metodo e relativo dispositivo per stimare il livello di rumore gaussiano bianco che corrompe un'immagine digitale
JP2012256202A (ja) * 2011-06-09 2012-12-27 Sony Corp 画像処理装置、および画像処理方法、並びにプログラム
JP5766077B2 (ja) * 2011-09-14 2015-08-19 キヤノン株式会社 ノイズ低減のための画像処理装置及び画像処理方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5729631A (en) * 1993-11-30 1998-03-17 Polaroid Corporation Image noise reduction system using a wiener variant filter in a pyramid image representation
US7590302B1 (en) * 2004-11-01 2009-09-15 Pixelworks, Inc. Image edge enhancement system and method
CN102006484A (zh) * 2009-08-31 2011-04-06 索尼公司 图像处理装置、图像处理方法以及程序

Also Published As

Publication number Publication date
CN104217416A (zh) 2014-12-17
EP2819091A2 (en) 2014-12-31
EP2819091A3 (en) 2015-04-15
US20140355877A1 (en) 2014-12-04
US9299163B2 (en) 2016-03-29
EP2819091B1 (en) 2019-05-01
JP2014235736A (ja) 2014-12-15
JP6314654B2 (ja) 2018-04-25

Similar Documents

Publication Publication Date Title
CN104217416B (zh) 灰度图像处理方法及其装置
Maini et al. Study and comparison of various image edge detection techniques
CN110264426B (zh) 图像畸变校正方法和装置
KR102144994B1 (ko) 영상의 노이즈를 저감하는 방법 및 이를 이용한 영상 처리 장치
JP2015225665A (ja) 画像ノイズ除去方法及び画像ノイズ除去装置
CN104794685B (zh) 一种实现图像去噪的方法及装置
WO2016197571A1 (en) Image interpolation device and method thereof
CN101499164B (zh) 一种基于单张低分辨率图像的图像插值重构方法
Lo et al. Joint trilateral filtering for depth map super-resolution
EP2816526B1 (en) Filtering method and apparatus for recovering an anti-aliasing edge
EP1958158A2 (en) Method for detecting streaks in digital images
EP2724316B1 (en) Method and device for processing of an image
Cho et al. Dictionary-based anisotropic diffusion for noise reduction
Kang et al. Fusing soft-decision-adaptive and bicubic methods for image interpolation
KR101262164B1 (ko) 저 해상도 깊이 맵으로부터 고 해상도 깊이 맵을 생성하는 방법 및 이를 기록한 기록매체
CN109325909B (zh) 一种图像放大方法和图像放大装置
WO2017128646A1 (zh) 一种图像处理的方法及装置
CN110363723B (zh) 改进图像边界效果的图像处理方法及装置
CN112150532A (zh) 图像处理的方法、装置、电子设备和计算机可读介质
Wang et al. Remote sensing image magnification study based on the adaptive mixture diffusion model
Storozhilova et al. 2.5 D extension of neighborhood filters for noise reduction in 3D medical CT images
WO2009090641A1 (en) Multi-scale representation of an out of focus image
Yu et al. Removing blur kernel noise via a hybrid ℓ p norm
JP4930638B2 (ja) 画像補正装置および画像補正方法
Jayamanmadharao et al. IMPLEMENTATION OF OBJECT ORIENTED APPROCH FOR NOISE REDUCTION USING FUZZY LOGIC

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170915

CF01 Termination of patent right due to non-payment of annual fee