CN107909555B - 一种保持锐度的网格噪声消除方法 - Google Patents

一种保持锐度的网格噪声消除方法 Download PDF

Info

Publication number
CN107909555B
CN107909555B CN201711201625.4A CN201711201625A CN107909555B CN 107909555 B CN107909555 B CN 107909555B CN 201711201625 A CN201711201625 A CN 201711201625A CN 107909555 B CN107909555 B CN 107909555B
Authority
CN
China
Prior art keywords
image
pixel
voting
sharpness
edge
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
Application number
CN201711201625.4A
Other languages
English (en)
Other versions
CN107909555A (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.)
China Daheng Group Inc Beijing Image Vision Technology Branch
Beijing Daheng Image Vision Co ltd
Original Assignee
China Daheng Group Inc Beijing Image Vision Technology Branch
Beijing Daheng Image Vision Co 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 China Daheng Group Inc Beijing Image Vision Technology Branch, Beijing Daheng Image Vision Co ltd filed Critical China Daheng Group Inc Beijing Image Vision Technology Branch
Priority to CN201711201625.4A priority Critical patent/CN107909555B/zh
Publication of CN107909555A publication Critical patent/CN107909555A/zh
Application granted granted Critical
Publication of CN107909555B publication Critical patent/CN107909555B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/73Deblurring; Sharpening
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20024Filtering details
    • G06T2207/20032Median filtering

Landscapes

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

Abstract

本发明公开了一种保持锐度的网格噪声消除方法,所述方法包括:基于对目标图像中各个像素是否属于平滑区域的判定,对目标图像中的平滑区域进行滤波处理。本发明采用投票法进行平滑区域筛选,从图像的一角为起点,逐个像素地对每个像素进行平滑区投票,并且将当前像素的投票值赋予其相关像素,以确定每个像素的投票值,因此,去除网格噪声的同时保持了图像的锐度,并且具有抗噪性。本发明参数很少,并且与图像亮度无关,对于不同网格性状的图像均有良好的处理效果。

Description

一种保持锐度的网格噪声消除方法
技术领域
本发明涉及一种图像质量提升技术,针对于目前含有网格模板噪声的光学传感器,该方法能够去除图像中网格噪声的同时保留图像边缘的锐度。
背景技术
网格噪声的产生是因为图像在水平和垂直方向存在周期性的亮度差异,这些差异具体表现有横竖条纹。同时具有横竖条纹的图像即为网格。由于图像缩放算法的原因,当图像缩小时网格会更加明显。
针对于含有网格噪声的图像,目前有三种方式进行处理。
1、提升相机的出图质量,对光学传感器进行改进。此方法成本较高,不易于实现。
2、对图像进行平滑处理,采用均值滤波的方式消除像素之间的固定差异,此方法会降低图像的锐度,尤其对于边缘部分,会变的模糊。
3、对图像首先进行平滑区域判断,再对平滑区域进行加权平均处理。如专利200810105709.2中所示,其首先计算周围八邻域相互之间的亮度差异,再与设定好的阈值进行比较来进行平滑区域判定,最后对平滑区域处的像素点与邻域内的全部或部分像素进行加权平均,结果作为更新后像素点值。此方法需要很多阈值和权值,普适性差,对于噪声仅依靠阈值来进行处理。阈值较小则会导致噪点被判定为非平滑区域。阈值较大则会导致弱边缘被判定为平滑区域。如果不同位置的网格亮度差异较大,便很难有一个良好的阈值筛选出所有的平滑区域进行处理。
现有技术中还有一些复杂的算法来去除网格噪声或在网格噪声去除过程中进行边缘检测,但是现有的去除网格噪声或边缘检测算法都会消耗大量的运算资源,导致图像处理的速度无法满足高速运行状况的应用。
发明内容
针对含有网格噪声的光学传感器及其产生的具有网格噪声的图像,本发明设计了一种具有普适性、抗噪性的方法,在保持原有图像锐度的同时消除网格噪声。
这里所提到的普适性是指针对于网格性状不同的光学传感器同样具有良好的处理效果。(网格性状不同指的是网格在周期、幅值和位置三个方面存在差异)
这里所提到的抗噪性是指针对于噪声情况不同的光学传感器同样具有良好的处理效果。(噪声会影响平滑区域的判定)
这里所提到的保持图像的锐度是指本发明对图像锐度的影响最小化,远小于其他的网格噪声去除方法。本发明的方法在网格噪声去除过程中对图像边缘不进行处理,因为边缘的变化趋势会掩盖网格噪声。
具体而言,本发明提供一种保持锐度的网格噪声消除方法,其特征在于,所述方法包括:基于对目标图像中各个像素是否属于平滑区域的判定,对目标图像中的平滑区域进行滤波处理。
优选地,所述方法包括下述步骤:
(1)、从图像的一角为起点,逐个像素地对每个像素进行平滑区投票,并且将当前像素的投票值赋予其相关像素,以确定每个像素的投票值;
(2)、读取或设定目标图像的投票值阈值,并根据投票值阈值对所述目标图像进行平滑区判定,以确定目标图像中每个像素是否属于边缘区;
(3)、对图像中属于边缘区的像素不进行滤波,对图像中不属于边缘区的像素进行滤波。
优选地,所述步骤(1)包括:
对于目标点a及其横向相关点b、c以及纵向相关点为d、f,如果(a-b)*(b-c)>0,则判定横向存在边缘,如果(a-d)*(d-f)>0,则判定纵向存在边缘。
优选地,在所述步骤(3)中对像素进行中值滤波。
优选地,所述方法还包括基于目标图像的复杂度进行图像或图像区域的分类。
优选地,所述方法还包括基于图像噪声水平调整投票阈值。
优选地,所述方法还包括,从图像的一角为起点,向图像的对角方向,逐个像素地对每个像素进行平滑区投票。该投票包括分别沿横向和纵向两个方向确定目标像素的相关像素,并且将相关像素的投票值赋予目标像素。
本发明的方法将目标点在横向和/或纵向的相关点的边缘检测结果赋予目标点,既简单易行,又准确可靠。
本发明在去除网格的同时保持了图像的锐度,并且具有抗噪性。对于不同性状的网格噪声具有普适性。并且,本发明的网格噪声消除方法对运算资源的要求低,能够适应工业生产中对图像处理的高速检测需求,在同样的检测效率情况下,能够有效降低系统资源成本。更为重要的是,本发明的网格噪声消除方法不依赖于图像亮度,不需要计算邻域的亮度差异。
附图说明
图1为本发明方法中边缘检测部分示意图;
图2本发明方法中边缘判定部分第一种投票方案示意图;
图3本发明方法中边缘判定部分第二种投票方案示意图;
图4本发明方法中所采用的中值滤波领域选取示意图;
图5为含有微弱网格的原图片;
图6利用均值滤波方法对微弱网格图像处理后的图片;
图7利用本发明方法对微弱网格图像处理后的图片。
图8-10分别为原图的垂直清晰度、现有方法处理后的垂直清晰度以及采用本发明方法处理后的垂直清晰度的视图。
图11中示出了细节较多的区域和细节较少的区域。
具体实施方式
以下结合附图及其实施例对本发明进行详细说明,但并不因此将本发明的保护范围限制在实施例描述的范围之中。
本发明采用2×2的中值滤波对图像进行非线性处理,为保持场景中边缘的锐度,在滤波之前进行平滑区域判定,对于非平滑区域将不进行中值滤波处理。
1非平滑区域投票(包括边缘检测过程和投票过程)。
1.1边缘检测
采用当前点右下方3×3范围内横纵两个方向进行边缘检测。
如图1所示,当前点为a。
横向相关点为b、c,如果(a-b)*(b-c)>0,
则判定横向存在边缘,否则为非边缘,即平滑区域。
纵向相关点为d,f,如果(a-d)*(d-f)>0,
则判定纵向存在边缘,否则为非边缘,即平滑区域。
这个处理的顺序是从左上角的开始,因为Sensor的数据是逐行读出的,在具体的实现中首先缓存两行,当第三行来第一个数据时,对第一行的第一个数据进行边缘判定,并且记录下两个方向上的检测情况。这样会节省资源。之后就可以根据上面和左面的四个相关点的结算结果进行投票了。简单说就是当前点的相关点,上面两个,左面两个边缘检测完成后会记录下来,当前点直接根据相关点和本身的边缘检测结果进行投票。
下述步骤1.2中对当前点的投票判定是在前两行和左面两个点检测完成后进行的。
1.2非平滑区的投票
如图2所示,当前点为C。判定相关点为A、B、E、D。这四个点在进行边缘检测时(见1.1)与C点相关,因此这四个点的边缘检测结果会关系到C的判定结果。ABCDE五个点共同为C点是否处在平滑区域进行投票。
针对于不同的场景,共有两种投票方式,本发明将场景分为两类场景,第一类场景为简单场景、第二类场景为复杂场景。对于场景的分类可以人工进行,也可以通过软件实现,比如,对图像进行一个特征提取,特征数量超过一定阈值的则定义为复杂场景。此外,这里的场景分类不仅限于对不同图像进行分类,还可以指的是对图像的不同区域也进行分类。
1)对于第一类图像或图像区域,每次与C相关的边缘检测计1票。
针对于细节较少的场景,即简单场景。图像的边缘结构比较明显,结构比较清晰,主要是低频和中频信息。
如图2所示,AB的纵向边缘检测、ED的横向边缘检测和C点的横纵双方向边缘检测会与C点相关(共计6次检测)。检测为边缘则给C点投1票。则C点最多获得6票,最少为0票。即,对于此种场景,进行目标相关点逐次计票。
2)对于第二类图像或图像区域,每个与C相关点的边缘检测计1票。
第二种投票方式主要是针对于细节较多的第二类场景,即复杂场景,采用该投票方式可以保护与边缘相邻的点,使之不进行中值滤波。
如下图11所示,对于细节较多的场景,边缘之间互相影响,第二种投票方式,不是相关边缘检测的次数结果计1票,而是相关点的判定情况计1票,这样如果相关点判定为边缘,当前点也会提高判定为边缘的概率。因为当前点的检测结果的权重降低了,只占了一票。边缘的相邻点更大概率被判定为边缘。举个例子,当前的上方只是横向判定为边缘,但是会为当前点是边缘投一票(虽然纵向没有判定为边缘)。
如图3所示,ABCDE 5个点在边缘检测时与C点相关。每个点在边缘检测时,横纵任一方向检测为边缘则认为此点为边缘,此时给C点投1票。则C点最多获得5票,最少获得0票。
2平滑区域判定和中值滤波
2.1平滑区域判定
对于每一个像素点,当票数大于判定阈值时,此点则被判定为非平滑区域,当票数小于等于判定阈值时,此点则被判定为平滑区域。
判定阈值的选取:阈值的范围为0到最大票数。(第一种投票方式最大为6票,第二种投票方式最大为5票。)推荐此阈值设定为2,因为处于边缘处的像素点至少会被投票2次。对于噪声比较大的图像,可以适当提高阈值。
2.2中值滤波
对于判定为平滑区域的点,对其进行2×2的中值滤波。将当前点的亮度值更新为该点右下方2×2邻域窗口内的所有像素点灰度值的中值。如图4所示,A点为当前点,四个点的亮度值由大到小的顺序下为A>B>C>D,则A点被替换为
Figure GDA0001544973960000071
本发明采用投票法进行平滑区域筛选,去除网格噪声的同时保持了图像的锐度,并且具有抗噪性。
下面给出了对图5中所示图像采用本发明的滤波方法和采用均值滤波方法进行滤波处理后的处理结果对比
一、边缘像素数量进行对比
表2-1图像边缘处灰度值由10%上升到90%范围内包含的像素数量
Figure GDA0001544973960000072
从表格中的对比可以看出,采用本实施例中的算法,处理后像素数量相差不大,或者减小,说明锐度保持不变或者轻微升高。相比之下,采用现有的均值滤波会使边缘处像素数量增加,锐度下降。
二、清晰度对比
表3-1为图像在MTF50处能够达到的分辨率。(这个是表示图像锐度的一个值。在表示相机图像解析力时,通常采用MTF50或者MTF50P。MTF50表示的是MTF为最大值的50%(即MTF=0.5)时,对应的空间频率)
Figure GDA0001544973960000081
其中LW/PH为线宽每图像高,代表着图像的清晰程度,值越高图像清晰度越好。
从表格中可以看出,采用本发明的噪声消除方法处理后清晰度基本不变,均值滤波后下降明显。
本发明参数很少,并且与图像亮度无关,对于不同网格性状的图像均有良好的处理效果。
如图5至图7中所示,图5为含有微弱网格的图片,图6为采用现有的均值滤波方法处理后的图像,而图7则是利用本发明方法对图5中的微弱网格图像处理后的图片。原图中缩小至51%后观察,网格明显,而采用本发明的方法去网格后网格消失。
图8-10分别为原图的垂直清晰度、现有方法处理后的垂直清晰度以及采用本发明方法处理后的垂直清晰度的视图。
图8、图9和图10示出了imatest的计算结果,表3-1和2-1的数据可以在图中看到。另外图中的MTF曲线的形状能表示我们的效果比中值滤波好,我们基本和原图一致,中值滤波的形状明显陡峭。
虽然,现有技术中也存在着一些边缘检测算法,但是现有的去除网格噪声或边缘检测算法都会消耗大量的运算资源,导致图像处理的速度无法满足高速运行状况的应用,而本发明的基于投票的边缘检测方法CPU占用资源非常低,运算量远低于现有的边缘检测算法。边缘检测是去网格算法的一部分,Sobel算子每个像素点需要11次运算,本文提到的算法需要8次运算,可以明显降低运算量。
以上所述仅为本发明的较佳实施例,并非对本发明做任何形式上的限制,凡在本发明的精神和原则之内,依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化与修饰,均仍属于本发明的保护范围之内。
虽然上面结合本发明的优选实施例对本发明的原理进行了详细的描述,本领域技术人员应该理解,上述实施例仅仅是对本发明的示意性实现方式的解释,并非对本发明包含范围的限定。实施例中的细节并不构成对本发明范围的限制,在不背离本发明的精神和范围的情况下,任何基于本发明技术方案的等效变换、简单替换等显而易见的改变,均落在本发明保护范围之内。

Claims (4)

1.一种保持锐度的网格噪声消除方法,其特征在于,所述方法包括:基于对目标图像中各个像素是否属于平滑区域的判定,对目标图像中的平滑区域进行滤波处理 ,其中,所述方法包括下述步骤:
(1)、从图像的左上角为起点,逐个像素地对每个像素进行平滑区投票,并且将当前像素的投票值赋予其相关像素,以确定每个像素的投票值,首先缓存两行像素,当第三行来第一个数据时,对第一行的第一个数据进行边缘判定,并且记录下两个方向上的检测情况,之后就根据上面和左面的四个相关点的结算结果进行投票;
(2)、读取或设定目标图像的投票值阈值,并根据投票值阈值对所述目标图像进行平滑区判定,以确定目标图像中每个像素是否属于边缘区;
(3)、对图像中属于边缘区的像素不进行滤波,对图像中不属于边缘区的像素进行滤波,
所述步骤(1)包括:
对于目标点a及其横向相关点b、c以及纵向相关点为d、f,如果(a-b)*(b-c)>0,则判定横向存在边缘,如果(a-d)*(d-f)>0,则判定纵向存在边缘。
2.根据权利要求1所述的保持锐度的网格噪声消除方法,其特征在于,在所述步骤(3)中对像素进行中值滤波。
3.根据权利要求1所述的保持锐度的网格噪声消除方法,其特征在于,所述方法还包括基于目标图像的复杂度进行图像或图像区域的分类。
4.根据权利要求1所述的保持锐度的网格噪声消除方法,其特征在于,所述方法还包括基于图像噪声水平调整投票值阈值。
CN201711201625.4A 2017-11-27 2017-11-27 一种保持锐度的网格噪声消除方法 Active CN107909555B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711201625.4A CN107909555B (zh) 2017-11-27 2017-11-27 一种保持锐度的网格噪声消除方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711201625.4A CN107909555B (zh) 2017-11-27 2017-11-27 一种保持锐度的网格噪声消除方法

Publications (2)

Publication Number Publication Date
CN107909555A CN107909555A (zh) 2018-04-13
CN107909555B true CN107909555B (zh) 2020-06-02

Family

ID=61848725

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711201625.4A Active CN107909555B (zh) 2017-11-27 2017-11-27 一种保持锐度的网格噪声消除方法

Country Status (1)

Country Link
CN (1) CN107909555B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109493357B (zh) * 2018-11-16 2021-06-18 浙江海洋大学 一种数码影像近景图像的边界提取方法
CN110327033B (zh) * 2019-04-04 2022-05-03 浙江工业大学 一种基于深度神经网络的心肌梗死心电图的筛查方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101150658A (zh) * 2007-08-22 2008-03-26 北京中星微电子有限公司 一种坏点自适应的网格噪声消除装置及方法
CN101272503A (zh) * 2008-04-30 2008-09-24 北京中星微电子有限公司 一种保持图像清晰度的网格噪声消除方法及其装置
CN101793501A (zh) * 2010-04-14 2010-08-04 华中科技大学 一种基于图像的输电线覆冰状态检测方法
CN101794446A (zh) * 2010-02-11 2010-08-04 东南大学 一种图像角点的线搜索式检测方法
CN103530878A (zh) * 2013-10-12 2014-01-22 北京工业大学 一种基于融合策略的边缘提取方法
AU2015283079A1 (en) * 2014-06-30 2016-11-24 Ventana Medical Systems, Inc. Detecting edges of a nucleus using image analysis

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8320670B1 (en) * 2008-08-29 2012-11-27 Adobe Systems Incorporated Hough transform method for linear ribbon and circular ring detection in the gradient domain

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101150658A (zh) * 2007-08-22 2008-03-26 北京中星微电子有限公司 一种坏点自适应的网格噪声消除装置及方法
CN101272503A (zh) * 2008-04-30 2008-09-24 北京中星微电子有限公司 一种保持图像清晰度的网格噪声消除方法及其装置
CN101794446A (zh) * 2010-02-11 2010-08-04 东南大学 一种图像角点的线搜索式检测方法
CN101793501A (zh) * 2010-04-14 2010-08-04 华中科技大学 一种基于图像的输电线覆冰状态检测方法
CN103530878A (zh) * 2013-10-12 2014-01-22 北京工业大学 一种基于融合策略的边缘提取方法
AU2015283079A1 (en) * 2014-06-30 2016-11-24 Ventana Medical Systems, Inc. Detecting edges of a nucleus using image analysis

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Road Network Extraction using Edge Detection and Spatial Voting;Beril Sırmacek 等;《Pattern Recognition》;20100826;第3113-3116页 *
基于排序统计的图像边缘增强滤波;李向吉 等;《天津大学学报》;19991130;第32卷(第6期);第687页第1节第1段 *
基于边缘连接的快速直线段检测算法;董晶 等;《光学学报》;20130331;第33卷(第3期);第0315003-1至0315003-8页 *

Also Published As

Publication number Publication date
CN107909555A (zh) 2018-04-13

Similar Documents

Publication Publication Date Title
CN102479383B (zh) 椒盐噪声去除方法及装置
CN109658424B (zh) 一种改进的鲁棒二维otsu阈值图像分割方法
CN106599783B (zh) 一种视频遮挡检测方法及装置
JP6505237B2 (ja) 画像処理装置
CN115063421B (zh) 极片区域检测方法及系统及装置及介质及缺陷检测方法
CN107292828B (zh) 图像边缘的处理方法和装置
CN111145105B (zh) 一种图像快速去雾方法、装置、终端及存储介质
CN109727202B (zh) 气体红外视频增强方法、装置及存储介质
CN109583435A (zh) 一种自适应投票车道线检测方法
CN114127784A (zh) 用于生成相机流的掩模的方法、计算机程序产品和计算机可读介质
CN107909555B (zh) 一种保持锐度的网格噪声消除方法
CN100367770C (zh) 一种去除视频孤立噪声点的方法
CN103996177A (zh) 一种无参考检测消除雪花噪声算法
CN114004850A (zh) 一种天空分割方法、图像去雾方法、电子设备和存储介质
CN115829956A (zh) 一种基于fpga的低照度视频增强硬件实现方法
CN104010114A (zh) 视频去噪方法和装置
CN115294149A (zh) 一种天文图像背景提取方法及系统
CN103514588A (zh) 图像增强方法和系统
CN106846262B (zh) 去除蚊式噪声的方法及系统
CN106874882B (zh) 路面病害边缘检测方法和装置
CN111161211A (zh) 一种图像检测方法及装置
CN115035311A (zh) 基于可见光与热红外融合的托辊检测方法
CN113781328A (zh) 一种sigma图像滤波方法及系统
KR101179500B1 (ko) 적응적 잡음 제거 방법 및 장치
CN115393782A (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
GR01 Patent grant
GR01 Patent grant