CN104794692B - 一种图像去锯齿的系统 - Google Patents

一种图像去锯齿的系统 Download PDF

Info

Publication number
CN104794692B
CN104794692B CN201510179591.8A CN201510179591A CN104794692B CN 104794692 B CN104794692 B CN 104794692B CN 201510179591 A CN201510179591 A CN 201510179591A CN 104794692 B CN104794692 B CN 104794692B
Authority
CN
China
Prior art keywords
mrow
mtd
mtr
pixel
edge direction
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
CN201510179591.8A
Other languages
English (en)
Other versions
CN104794692A (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.)
Beijing Jilang Semiconductor Technology Co Ltd
Original Assignee
Institute of Automation of Chinese Academy of Science
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 Institute of Automation of Chinese Academy of Science filed Critical Institute of Automation of Chinese Academy of Science
Priority to CN201510179591.8A priority Critical patent/CN104794692B/zh
Publication of CN104794692A publication Critical patent/CN104794692A/zh
Application granted granted Critical
Publication of CN104794692B publication Critical patent/CN104794692B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Processing (AREA)
  • Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)

Abstract

本发明公开了一种图像去锯齿系统包括边缘检测单元、方向滤波器、小角度方向滤波器、方向滤波融合单元、结果融合单元,边缘检测单元依据输入图像输出像素的边缘方向和边缘置信度,方向滤波器依据输入图像和像素的边缘方向输出方向滤波结果,小角度方向滤波器依据输入图像和像素的边缘方向输出小角度方向滤波结果;方向滤波融合单元依据各像素的边缘方向将方向滤波器输出结果和小角度方向滤波器输出结果加权组合,输出方向滤波加权组合结果;结果融合单元依据各像素的边缘置信度、输入图像、方向滤波加权组合结果输出优化图像。本发明能很好地处理任意方向的方向滤波,并对低角度锯齿具有良好的效果。

Description

一种图像去锯齿的系统
技术领域
本发明涉及图像视频处理技术领域,尤其涉及对图像、视频、图形进行锯齿去除的领域。
背景技术
图像和视频处理中的缩放和去隔行等处理会给图像的边缘带来锯齿效应,即在源图像中看上去直的边缘在处理后的图像上,在边缘附近出现了锯齿效应。在图形渲染中,处理一条直线时,如果处理不当,渲染出的直线会出现锯齿。当图像出现锯齿时,需要对锯齿进行去除。
一种传统的去锯齿系统如图1所示,先检测图像的边缘,根据边缘方向在方向滤波器组中选取一个滤波器,进行沿边缘方向的滤波,最后利用边缘的置信度对源图像和方向滤波融合单元组合的结果进行组合,输出最终结果。这种去锯齿系统有几个缺点,第一个缺点是方向往往经过量化,量化成有限的几个的方向,每个方向对应一个沿方向的滤波器,这样对于落在量化方向中间的方向,只能选取就近的量化方向进行滤波,或者对邻近的几个量化方向滤波结果进行组合,准确性受到影响。第二个缺点是在处理低角度(即接近水平方向的角度)时,对于低于量化方向中最低角度的方向,无法有效处理该方向出现的锯齿。第三个缺点是在处理小角度(即接近水平方向的角度)时,参与滤波的像素相隔较远,容易造成错误。
发明内容
为了解决上述技术问题,本发明提出了一种图像去锯齿系统,该系统能很好地处理任意方向的方向滤波,并对低角度锯齿具有良好的效果。
本发明提出的一种图像去锯齿系统,包括边缘检测单元、方向滤波器、结果融合单元,边缘检测单元依据输入图像输出像素的边缘方向和边缘置信度,方向滤波器依据输入图像和像素的边缘方向输出方向滤波结果,结果融合单元依据像素的边缘置信度、输入图像和方向滤波结果输出优化图像,所述的方向滤波器根据待插值像素的边缘方向,计算出穿过该像素的沿边缘方向的一条线与邻域水平扫描线或邻域垂直扫描线的交点坐标,并利用一维插值的方法,插出位于交点的临时像素,再利用滤波器对沿边缘方向上的各个临时像素进行滤波,输出方向滤波结果。
为了实现更好的去锯齿效果,本发明的一种图像去锯齿系统还包括小角度方向滤波器、方向滤波融合单元;小角度方向滤波器依据输入图像和像素的边缘方向输出小角度方向滤波结果;方向滤波融合单元依据像素的边缘方向将方向滤波器输出结果和小角度方向滤波器输出结果加权组合,输出方向滤波加权组合结果;结果融合单元依据像素的边缘置信度、输入图像、方向滤波加权组合结果输出优化图像;
所述方向滤波器通过计算穿过待插值像素的沿边缘方向的一条线与邻域水平扫描线的交点坐标来得到位于交点的临时像素时,所述的小角度方向滤波器依据待插值像素的边缘方向,计算出穿过该像素的沿边缘方向的一条线与邻域垂直扫描线的交点坐标,并利用一维插值的方法,插出位于交点的临时像素;再利用一维滤波器对沿边缘方向上的各个临时像素进行滤波,输出小角度方向滤波结果。
所述方向滤波器通过计算穿过待插值像素的沿边缘方向的一条线与邻域垂直扫描线的交点坐标来得到位于交点的临时像素时,所述的小角度方向滤波器依据待插值像素的边缘方向,计算出穿过该像素的沿边缘方向的一条线与邻域水平扫描线的交点坐标,并利用一维插值的方法,插出位于交点的临时像素;再利用一维滤波器对沿边缘方向上的各个临时像素进行滤波,输出小角度方向滤波结果。
本发明具有如下有益效果
(1)由于不进行方向量化,可以处理任意方向的锯齿;
(2)可以很好地处理小角度的锯齿,可以处理的角度可以任意接近水平方向。
附图说明
图1传统的去锯齿系统框图;
图2本发明图像去锯齿系统框图;
图3本发明边缘检测单元框图;
图4归一化的边缘置信度计算曲线示意图;
图5本发明方向滤波器结构示意图;
图6本发明方向滤波示意图;
图7本发明在小角度进行方向滤波示意图;
图8本发明小角度方向滤波示意图;
图9本发明小角度方向滤波器结构示意图;
图10本发明方向滤波的融合权重计算曲线示意图;
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
如图2所示,本实施例的图像去锯齿系统包括边缘检测单元、方向滤波器、小角度方向滤波器、方向滤波融合单元、结果融合单元,边缘检测单元依据输入图像f输出像素的边缘方向D和边缘置信度Redge,方向滤波器依据输入图像f和像素的边缘方向D输出方向滤波结果fdh,小角度方向滤波器依据输入图像f和像素的边缘方向D输出小角度方向滤波结果fdl;方向滤波融合单元依据像素的边缘方向D将方向滤波器输出结果fdh和小角度方向滤波器输出结果fdl加权组合,输出方向滤波加权组合结果fd;结果融合单元依据像素的边缘置信度Redge、输入图像f、方向滤波加权组合结果fd输出优化图像f’。
如图3所示,边缘检测单元包括水平梯度计算模块、垂直梯度计算模块、局部协方差矩阵据算模块、特征值和特征向量计算模块、边缘方向计算模块、边缘置信度计算模块,水平梯度计算模块依据输入图像f计算图像像素的水平梯度gx,垂直梯度计算模块依据输入图像f计算图像像素的垂直梯度gy,局部协方差矩阵据算模块依据像素的水平梯度gx和垂直梯度gy计算像素的局部协方差Mc,特征值和特征向量计算模块依据局部协方差Mc计算像素局部协方差矩阵的特征值λ和特征向量e,边缘方向计算模块依据特征值λ计算像素的边缘方向d,边缘置信度计算模块用于依据特征向量e计算像素的边缘置信度Redge
其中水平梯度gx和垂直梯度gy可以用sobel梯度算子或其他梯度算子得到。局部协方差Mc的计算公式如式(1),假设用于计算局部协方差的窗口大小为(2K+1)*(2N+1),则位于(i,j)的像素其局部协方差为
协方差矩阵Mc特征值为λ1和λ2,λ1≥λ2,特征向量为eλ1和eλ2。则边缘方向D和较小特征值λ2对应的特征向量相同,即如式(2)所示
D=eλ2 (2)
边缘方向的置信度R计算如式(3)所示,其中Tr表示矩阵的迹,Det表示矩阵的行列式。
归一化的边缘置信度Redge计算如图4所示,其公式如图4所示
其中R_T1和R_T2为预设的两个阈值。
方向滤波器采用邻域水平扫描线和邻域垂直扫描线均可达到相应的滤波效果,本实施例采用邻域水平扫描线进行详细说明。所述的方向滤波器根据待处理像素的边缘方向D,计算出穿过该像素的沿边缘方向的一条线与邻域水平扫描线的交点坐标,并利用一维插值的方法,插出位于交点的临时像素,再利用滤波器对沿边缘方向上的各个临时像素进行滤波,输出方向滤波结果fdh
为方便叙述和理解,本实施例仅举穿过该像素的沿边缘方向的一条线与邻域水平扫描线的交点计算单元为5个的例子。
方向滤波示意图如图6所示,其中黑点为原始图像像素,P为待滤波像素,根据待滤波像素P的边缘方向,计算出穿过该像素的沿边缘方向的一条线Lp与邻域水平扫描线(如图6所示的H0,H1,H2,H3,H4)的交点坐标,并利用水平插值的方法,插出位于交点的临时像素(P0,P1,P2,P3,P4),如图6的五角星所示,再利用沿方向的一维滤波器对各个边缘方向上的临时像素进行滤波,得到最终的方向滤波结果。
如图5所示方向滤波器中的边缘方向与水平扫描线交点计算单元一共有m个(本实施例中m=5),其中边缘方向与水平扫描线交点计算单元n(n=0,1...,m-1)计算出边缘方向与水平扫描线Hn(n=0,1...,m-1)之间的交点Pn的坐标,即边缘方向与水平扫描线交点计算单元0计算出边缘方向与水平扫描线H0之间的交点P0的坐标,边缘方向与水平扫描线交点计算单元4计算出边缘方向与水平扫描线H4之间的交点P4的坐标。交点坐标计算过程如下,设待滤波的像素P的坐标为(i,j),设输入方向D=[v1,v2]T,则穿过该像素的沿边缘方向的一条线Lp与水平线Hn(n=0,1,2,3,4)的交点Pn(n=0,1,2,3,4)的坐标(Pn.y,Pn,x)如公式(5)所示
如图5所示方向滤波器中的一维插值单元接受图像输入,以及边缘方向与水平扫描线的交点计算单元计算出的交点坐标,用一维插值的方法插出边缘方向上的临时像素点Pn。一维插值利用水平邻近的几个原始像素的加权组合计算出待插值点的值。具体一维插值的方法可以选取三次插值,线性插值,及多相位滤波器插值方法。
如图5所示方向滤波器中的一维滤波单元接受5个位于边缘方向上的临时像素点Pn(n=0,1,2,3,4)的输入,对5个临时像素点进行一维滤波得到最终的滤波结果fdh,如式(6)所示
其中cn为预设的进行滤波的系数,可以取做cn=1/5.即均值滤波。
如果用图6所示的方向滤波方式对图像进行滤波,当像素的边缘方向为小角度时(即接近水平角度)时,边缘方向上的临时像素点P0,P1,P3,P4与待滤波像素P之间距离较远,如图7所示,相距较远的像素可能已经不属于待滤波像素所处的边缘,让非边缘的像素也参与平滑滤波,有可能出现错误。
为解决这个问题,本实施例采用如图8所示的小角度方向滤波方法,根据待滤波像素P的边缘方向D,此时计算出穿过该像素的沿边缘方向的一条线Lp与邻域垂直扫描线(如图8所示的V0,V1,V2,V3,V4)的交点坐标,并利用一维插值的方法,插出位于交点的临时像素P0、P1、P2、P3、P4(如图8的五角星所示),再利用一维滤波器对沿边缘方向上的各个临时像素进行滤波,输出小角度方向滤波结果fdl
小角度方向滤波器单元的结构图如图9所示,接受图像f和边缘检测单元计算的方向D的输入,输出方向滤波的结果fdl。为方便叙述和理解,本实施例仅举本实施例仅举穿过该像素的沿边缘方向的一条线与邻域水平扫描线的交点计算单元为5个的例子。
如图9所示小角度方向滤波器中的边缘方向与垂直扫描线交点计算单元一共有m个(本实施例中m=5),其中边缘方向与垂直扫描线交点计算单元n(n=0,1...,m-1)计算出边缘方向与垂直扫描线Vn(n=0,1...,m-1)之间的交点Pn的坐标,即边缘方向与垂直扫描线交点计算单元0计算出边缘方向与垂直扫描线V0之间的交点P0的坐标,边缘方向与垂直扫描线交点计算单元4计算出边缘方向与垂直扫描线V4之间的交点P4的坐标。交点坐标计算过程如下,设待滤波的像素P的坐标为(i,j),设输入方向D=[v1,v2]T,则穿过该像素的沿边缘方向的一条线Lp与垂直扫描线Vn(n=0,1,2,3,4)的交点Pn(n=0,1,2,3,4)的坐标(Pn.y,Pn,x)如公式(7)所示
如图9所示小角度方向滤波器中的一维插值单元接受图像输入,以及边缘方向与垂直扫描线的交点计算单元计算出的交点坐标,用一维插值的方法插出位于交点的临时像素点Pn。一维插值利用垂直方向上邻近的几个原始像素的加权组合计算出待插值点的值。具体一维插值的方法可以选取三次插值,线性插值,及多相位滤波器插值方法。
如图9所示小角度方向滤波器中的一维滤波单元接受5个位于边缘方向上的临时像素点Pn(n=0,1,2,3,4)的输入,对5个临时像素点进行一维滤波得到最终的滤波结果fdl,如式(8)所示
其中cn为预设的进行滤波的系数,可以取做|cn=1/5.即均值滤波。
如图2所示的图像去锯齿系统中的方向滤波融合单元接受方向滤波器的输入fdh和小角度方向滤波器的输入fdl,以及边缘方向D=[v1,v2]T输入,v1代表方向向量的水平分量,v2代表代表方向向量的垂直分量。根据边缘方向的角度进行加权组合,如式(9)所示,其中w为加权的权重
fd=(1-w)*fdh+w*fdl (9)
其中权重计算如图10所示,其计算公式如(10)所示
其中θ为边缘方向D对应的角度,θ_T1和θ_T2为预设的两个阈值。
边缘方向D对应的角度θ的计算公式如(11)所示
结果融合单元接受图像输入f以及方向滤波结果fd,以及边缘置信度Redge,根据边缘置信度对fd和f进行加权融合,如式(12)所示
fd=(1-Redge)*f+Redge*fd (12)

Claims (10)

1.一种图像去锯齿系统,包括边缘检测单元、方向滤波器、结果融合单元,所述边缘检测单元依据输入图像输出像素的边缘方向和边缘置信度,所述方向滤波器依据所述输入图像和像素的边缘方向输出方向滤波结果,所述结果融合单元依据所述像素的边缘置信度、输入图像和方向滤波结果输出优化图像,其特征在于,所述的方向滤波器根据待处理像素的边缘方向,计算出穿过该像素的沿边缘方向的一条线与邻域水平扫描线或邻域垂直扫描线的交点坐标,并利用一维插值的方法,插出位于交点的临时像素,再利用滤波器对沿边缘方向上的各个临时像素进行滤波,输出方向滤波结果;
所述系统还包括小角度方向滤波器、方向滤波融合单元;所述小角度方向滤波器依据所述输入图像和像素的边缘方向输出小角度方向滤波结果;所述方向滤波融合单元依据所述像素的边缘方向将所述方向滤波器输出结果和所述小角度方向滤波器输出结果加权组合,输出方向滤波加权组合结果;所述结果融合单元依据所述像素的边缘置信度、输入图像、方向滤波加权组合结果输出优化图像;
所述方向滤波器通过计算穿过所述待处理像素的沿边缘方向的一条线与邻域水平扫描线的交点坐标来得到位于交点的临时像素时,所述的小角度方向滤波器依据所述待处理像素的边缘方向,计算出穿过该像素的沿边缘方向的一条线与邻域垂直扫描线的交点坐标,并利用一维插值的方法,插出位于交点的临时像素;再利用一维滤波器对沿边缘方向上的各个临时像素进行滤波,输出小角度方向滤波结果;
所述方向滤波器通过计算穿过所述待处理像素的沿边缘方向的一条线与邻域垂直扫描线的交点坐标来得到位于交点的临时像素时,所述的小角度方向滤波器依据所述待处理像素的边缘方向,计算出穿过该像素的沿边缘方向的一条线与邻域水平扫描线的交点坐标,并利用一维插值的方法,插出位于交点的临时像素;再利用一维滤波器对沿边缘方向上的各个临时像素进行滤波,输出小角度方向滤波结果。
2.如权利要求1所述的一种图像去锯齿系统,其特征在于,所述边缘检测单元包括水平梯度计算模块、垂直梯度计算模块、局部协方差矩阵据计算模块、特征值和特征向量计算模块、边缘方向计算模块、边缘置信度计算模块,所述水平梯度计算模块依据所述输入图像计算图像像素的水平梯度,所述垂直梯度计算模块依据所述输入图像计算图像像素的垂直梯度,所述局部协方差矩阵据计算模块依据所述像素的水平梯度和垂直梯度计算所述像素的局部协方差,所述特征值和特征向量计算模块依据所述局部协方差计算像素局部协方差矩阵的特征值和特征向量,所述边缘方向计算模块依据所述特征值计算像素的边缘方向,所述边缘置信度计算模块用于依据所述特征向量计算像素的边缘置信度。
3.如权利要求2所述的一种图像去锯齿系统,其特征在于,所述局部协方差的计算公式为
<mrow> <msub> <mi>M</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mrow> <msubsup> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mo>-</mo> <mi>N</mi> </mrow> <mi>N</mi> </msubsup> <msubsup> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mo>-</mo> <mi>K</mi> </mrow> <mi>K</mi> </msubsup> <msup> <mrow> <mo>(</mo> <msub> <mi>g</mi> <mi>x</mi> </msub> <mo>(</mo> <mrow> <mi>i</mi> <mo>+</mo> <mi>k</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mi>n</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mtd> <mtd> <mrow> <msubsup> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mo>-</mo> <mi>N</mi> </mrow> <mi>N</mi> </msubsup> <msubsup> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mo>-</mo> <mi>K</mi> </mrow> <mi>K</mi> </msubsup> <msub> <mi>g</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mi>k</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>*</mo> <msub> <mi>g</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mi>k</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mo>-</mo> <mi>N</mi> </mrow> <mi>N</mi> </msubsup> <msubsup> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mo>-</mo> <mi>K</mi> </mrow> <mi>K</mi> </msubsup> <msub> <mi>g</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mi>k</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>*</mo> <msub> <mi>g</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mi>k</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <msubsup> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mo>-</mo> <mi>N</mi> </mrow> <mi>N</mi> </msubsup> <msubsup> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mo>-</mo> <mi>K</mi> </mrow> <mi>K</mi> </msubsup> <msup> <mrow> <mo>(</mo> <msub> <mi>g</mi> <mi>y</mi> </msub> <mo>(</mo> <mrow> <mi>i</mi> <mo>+</mo> <mi>k</mi> <mo>,</mo> <mi>j</mi> <mo>+</mo> <mi>n</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中,所述局部协方差的窗口大小为(2K+1)*(2N+1),Mc(i,j)为坐标为(i,j)的像素的局部协方差,gx为像素的水平梯度,gy为像素的垂直梯度;
所述边缘方向的计算公式为D=eλ2,其中,所述D为边缘方向,eλ2为协方差矩阵的较小特征值对应的特征向量;
所述边缘置信度的计算公式为
<mrow> <msub> <mi>R</mi> <mrow> <mi>e</mi> <mi>d</mi> <mi>g</mi> <mi>e</mi> </mrow> </msub> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mtable> <mtr> <mtd> <mrow> <mi>i</mi> <mi>f</mi> </mrow> </mtd> <mtd> <mrow> <mi>R</mi> <mo>&gt;</mo> <mi>R</mi> <mo>_</mo> <mi>T</mi> <mn>2</mn> </mrow> </mtd> </mtr> </mtable> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mtable> <mtr> <mtd> <mrow> <mi>i</mi> <mi>f</mi> </mrow> </mtd> <mtd> <mrow> <mi>R</mi> <mo>&lt;</mo> <mi>R</mi> <mo>_</mo> <mi>T</mi> <mn>1</mn> </mrow> </mtd> </mtr> </mtable> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mi>R</mi> <mo>-</mo> <mi>R</mi> <mo>_</mo> <mi>T</mi> <mn>1</mn> </mrow> <mrow> <mi>R</mi> <mo>_</mo> <mi>T</mi> <mn>2</mn> <mo>-</mo> <mi>R</mi> <mo>_</mo> <mi>T</mi> <mn>1</mn> </mrow> </mfrac> </mtd> <mtd> <mrow> <mi>e</mi> <mi>l</mi> <mi>s</mi> <mi>e</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中,所述Redge为归一化的边缘方向的置信度,R_T1和R_T2为预设的两个阈值,R为边缘方向置信度,其中λ1和λ2分别为所述协方差矩阵Mc的两个特征值为,且λ1≥λ2
4.如权利要求3所述的一种图像去锯齿系统,其特征在于,所述的一维插值的方法可以采用三次插值方法、线性插值方法或多相位滤波器插值方法。
5.如权利要求4所述的一种图像去锯齿系统,其特征在于,待处理像素的沿边缘方向的一条线与邻域水平扫描线交点坐标的计算公式为:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>-</mo> <mi>m</mi> <mo>.</mo> <mi>y</mi> <mo>=</mo> <mi>i</mi> <mo>-</mo> <mi>m</mi> <mo>;</mo> </mrow> </mtd> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>-</mo> <mi>m</mi> <mo>.</mo> <mi>x</mi> <mo>=</mo> <mi>j</mi> <mo>+</mo> <mi>m</mi> <mo>*</mo> <mfrac> <mrow> <mi>v</mi> <mn>1</mn> </mrow> <mrow> <mi>v</mi> <mn>2</mn> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>-</mo> <mn>2.</mn> <mi>y</mi> <mo>=</mo> <mi>i</mi> <mo>-</mo> <mn>2</mn> <mo>;</mo> </mrow> </mtd> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>-</mo> <mn>2.</mn> <mi>x</mi> <mo>=</mo> <mi>j</mi> <mo>+</mo> <mn>2</mn> <mo>*</mo> <mfrac> <mrow> <mi>v</mi> <mn>1</mn> </mrow> <mrow> <mi>v</mi> <mn>2</mn> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>-</mo> <mn>1.</mn> <mi>y</mi> <mo>=</mo> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>;</mo> </mrow> </mtd> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>-</mo> <mn>1.</mn> <mi>x</mi> <mo>=</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> <mo>*</mo> <mfrac> <mrow> <mi>v</mi> <mn>1</mn> </mrow> <mrow> <mi>v</mi> <mn>2</mn> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>.</mo> <mi>y</mi> <mo>=</mo> <mi>i</mi> <mo>;</mo> </mrow> </mtd> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>.</mo> <mi>x</mi> <mo>=</mo> <mi>j</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>+</mo> <mn>1.</mn> <mi>y</mi> <mo>=</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>;</mo> </mrow> </mtd> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>+</mo> <mn>1.</mn> <mi>x</mi> <mo>=</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>*</mo> <mfrac> <mrow> <mi>v</mi> <mn>1</mn> </mrow> <mrow> <mi>v</mi> <mn>2</mn> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>+</mo> <mn>2.</mn> <mi>y</mi> <mo>=</mo> <mi>i</mi> <mo>+</mo> <mn>2</mn> <mo>;</mo> </mrow> </mtd> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>+</mo> <mn>2.</mn> <mi>x</mi> <mo>=</mo> <mi>j</mi> <mo>-</mo> <mn>2</mn> <mo>*</mo> <mfrac> <mrow> <mi>v</mi> <mn>1</mn> </mrow> <mrow> <mi>v</mi> <mn>2</mn> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>+</mo> <mi>m</mi> <mo>.</mo> <mi>y</mi> <mo>=</mo> <mi>i</mi> <mo>+</mo> <mi>m</mi> <mo>;</mo> </mrow> </mtd> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>+</mo> <mi>m</mi> <mo>.</mo> <mi>x</mi> <mo>=</mo> <mi>j</mi> <mo>-</mo> <mi>m</mi> <mo>*</mo> <mfrac> <mrow> <mi>v</mi> <mn>1</mn> </mrow> <mrow> <mi>v</mi> <mn>2</mn> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced>
其中(Pn.y,Pn.x)为坐标为(i,j)的待处理像素P的边缘方向与待处理像素所在的水平扫描线交点,m为大于等于2的自然数,v1代表边缘方向向量的水平分量,v2代表代表边缘方向向量的垂直分量。
6.如权利要求4所述的一种图像去锯齿系统,其特征在于,待处理像素的沿边缘方向的一条线与邻域垂直扫描线交点坐标的计算公式为:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>-</mo> <mi>m</mi> <mo>.</mo> <mi>y</mi> <mo>=</mo> <mi>i</mi> <mo>-</mo> <mi>m</mi> <mo>*</mo> <mfrac> <mrow> <mi>v</mi> <mn>2</mn> </mrow> <mrow> <mi>v</mi> <mn>1</mn> </mrow> </mfrac> <mo>;</mo> </mrow> </mtd> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>-</mo> <mi>m</mi> <mo>.</mo> <mi>x</mi> <mo>=</mo> <mi>j</mi> <mo>+</mo> <mi>m</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>-</mo> <mn>2.</mn> <mi>y</mi> <mo>=</mo> <mi>i</mi> <mo>-</mo> <mn>2</mn> <mo>*</mo> <mfrac> <mrow> <mi>v</mi> <mn>2</mn> </mrow> <mrow> <mi>v</mi> <mn>1</mn> </mrow> </mfrac> <mo>;</mo> </mrow> </mtd> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>-</mo> <mn>2.</mn> <mi>x</mi> <mo>=</mo> <mi>j</mi> <mo>+</mo> <mn>2</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>-</mo> <mn>1.</mn> <mi>y</mi> <mo>=</mo> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>*</mo> <mfrac> <mrow> <mi>v</mi> <mn>2</mn> </mrow> <mrow> <mi>v</mi> <mn>1</mn> </mrow> </mfrac> <mo>;</mo> </mrow> </mtd> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>-</mo> <mn>1.</mn> <mi>x</mi> <mo>=</mo> <mi>j</mi> <mo>+</mo> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>.</mo> <mi>y</mi> <mo>=</mo> <mi>i</mi> <mo>;</mo> </mrow> </mtd> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>.</mo> <mi>x</mi> <mo>=</mo> <mi>j</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>+</mo> <mn>1.</mn> <mi>y</mi> <mo>=</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>*</mo> <mfrac> <mrow> <mi>v</mi> <mn>2</mn> </mrow> <mrow> <mi>v</mi> <mn>1</mn> </mrow> </mfrac> <mo>;</mo> </mrow> </mtd> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>+</mo> <mn>1.</mn> <mi>x</mi> <mo>=</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>+</mo> <mn>2.</mn> <mi>y</mi> <mo>=</mo> <mi>i</mi> <mo>+</mo> <mn>2</mn> <mo>*</mo> <mfrac> <mrow> <mi>v</mi> <mn>2</mn> </mrow> <mrow> <mi>v</mi> <mn>1</mn> </mrow> </mfrac> <mo>;</mo> </mrow> </mtd> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>+</mo> <mn>2.</mn> <mi>x</mi> <mo>=</mo> <mi>j</mi> <mo>-</mo> <mn>2</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>+</mo> <mi>m</mi> <mo>.</mo> <mi>y</mi> <mo>=</mo> <mi>i</mi> <mo>+</mo> <mi>m</mi> <mo>*</mo> <mfrac> <mrow> <mi>v</mi> <mn>2</mn> </mrow> <mrow> <mi>v</mi> <mn>1</mn> </mrow> </mfrac> <mo>;</mo> </mrow> </mtd> <mtd> <mrow> <mi>P</mi> <mi>n</mi> <mo>+</mo> <mi>m</mi> <mo>.</mo> <mi>x</mi> <mo>=</mo> <mi>j</mi> <mo>-</mo> <mi>m</mi> </mrow> </mtd> </mtr> </mtable> </mfenced>
其中(Pn.y,Pn.x)为坐标为(i,j)的待处理像素P的边缘方向与待处理像素所在垂直扫描线交点,m为大于等于2的自然数,v1代表边缘方向向量的水平分量,v2代表代表边缘方向向量的垂直分量。
7.如权利要求5或6所述的一种图像去锯齿系统,其特征在于,对各个临时像素进行滤波的滤波器为一维滤波器。
8.如权利要求7所述的一种图像去锯齿系统,其特征在于,方向滤波加权组合结果为fd=(1-w)*fdh+w*fdl,其中fd为方向滤波加权组合结果,w为加权组合权重,fdh方向滤波结果,fdl为小角度方向滤波结果。
9.如权利要求8所述的一种图像去锯齿系统,其特征在于,方向滤波融合单元中的加权组合权重w的计算公式为
<mrow> <mi>w</mi> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mtable> <mtr> <mtd> <mrow> <mi>i</mi> <mi>f</mi> </mrow> </mtd> <mtd> <mrow> <mi>&amp;theta;</mi> <mo>&lt;</mo> <mi>&amp;theta;</mi> <mo>_</mo> <mi>T</mi> <mn>1</mn> </mrow> </mtd> </mtr> </mtable> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mtable> <mtr> <mtd> <mrow> <mi>i</mi> <mi>f</mi> </mrow> </mtd> <mtd> <mrow> <mi>&amp;theta;</mi> <mo>&gt;</mo> <mi>&amp;theta;</mi> <mo>_</mo> <mi>T</mi> <mn>2</mn> </mrow> </mtd> </mtr> </mtable> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mi>&amp;theta;</mi> <mo>_</mo> <mi>T</mi> <mn>2</mn> <mo>-</mo> <mi>&amp;theta;</mi> </mrow> <mrow> <mi>&amp;theta;</mi> <mo>_</mo> <mi>T</mi> <mn>2</mn> <mo>-</mo> <mi>&amp;theta;</mi> <mo>_</mo> <mi>T</mi> <mn>1</mn> </mrow> </mfrac> </mtd> <mtd> <mrow> <mi>e</mi> <mi>l</mi> <mi>s</mi> <mi>e</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中θ为边缘方向D对应的角度,θ_T1和θ_T2为预设的两个阈值。
10.如权利要求9所述的一种图像去锯齿系统,其特征在于,优化图像的加权融合公式为
f′=(1-Redge)*f+Redge*fd
其中f′为优化图像,Redge为归一化的边缘方向的置信度,f为输入图像,fd为方向滤波加权组合结果。
CN201510179591.8A 2015-04-16 2015-04-16 一种图像去锯齿的系统 Active CN104794692B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510179591.8A CN104794692B (zh) 2015-04-16 2015-04-16 一种图像去锯齿的系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510179591.8A CN104794692B (zh) 2015-04-16 2015-04-16 一种图像去锯齿的系统

Publications (2)

Publication Number Publication Date
CN104794692A CN104794692A (zh) 2015-07-22
CN104794692B true CN104794692B (zh) 2017-08-25

Family

ID=53559472

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510179591.8A Active CN104794692B (zh) 2015-04-16 2015-04-16 一种图像去锯齿的系统

Country Status (1)

Country Link
CN (1) CN104794692B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10410326B2 (en) 2015-04-16 2019-09-10 Institute Of Automation Chinese Academy Of Sciences Image anti-aliasing system
CN106886981B (zh) * 2016-12-30 2020-02-14 中国科学院自动化研究所 基于边缘检测的图像边缘增强方法和系统
CN107705256B (zh) * 2017-09-13 2021-03-16 西南交通大学 一种基于接触网语义的前向车载视频图像增强方法
CN109544459A (zh) * 2017-09-21 2019-03-29 中国电信股份有限公司 图像锯齿处理方法和装置以及计算机可读存储介质
CN107993213B (zh) * 2017-12-14 2021-08-27 天津鹿鹿通健康科技有限公司 一种去锯齿图像缩放器及其去锯齿方法
CN109993693B (zh) * 2017-12-29 2023-04-25 澜至电子科技(成都)有限公司 用于对图像进行插值的方法和装置
US10825154B2 (en) 2019-03-19 2020-11-03 Apple Inc. Directional bilateral filtering with improved noise reduction along edges
CN113837955B (zh) * 2021-08-17 2024-08-02 每平每屋(上海)科技有限公司 图像抗锯齿处理方法及电子设备
CN115984091A (zh) * 2021-10-14 2023-04-18 北京字跳网络技术有限公司 图像处理方法、装置、电子设备及可读存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101777176A (zh) * 2009-01-13 2010-07-14 北京大学 一种去除挂网图像锯齿的方法及装置
CN101790069A (zh) * 2010-03-09 2010-07-28 周艇 一种基于图像边缘方向的尺度变换方法
CN102025960A (zh) * 2010-12-07 2011-04-20 浙江大学 一种基于自适应插值的运动补偿去隔行方法
WO2012004830A1 (ja) * 2010-07-06 2012-01-12 三菱電機株式会社 画像処理装置
CN102750668A (zh) * 2012-07-04 2012-10-24 西南交通大学 结合局部方向特征的数字图像三倍插值放大方法
CN103093418A (zh) * 2013-02-21 2013-05-08 深圳市晶日盛科技有限公司 一种改进的数字图像缩放方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101777176A (zh) * 2009-01-13 2010-07-14 北京大学 一种去除挂网图像锯齿的方法及装置
CN101790069A (zh) * 2010-03-09 2010-07-28 周艇 一种基于图像边缘方向的尺度变换方法
WO2012004830A1 (ja) * 2010-07-06 2012-01-12 三菱電機株式会社 画像処理装置
CN102025960A (zh) * 2010-12-07 2011-04-20 浙江大学 一种基于自适应插值的运动补偿去隔行方法
CN102750668A (zh) * 2012-07-04 2012-10-24 西南交通大学 结合局部方向特征的数字图像三倍插值放大方法
CN103093418A (zh) * 2013-02-21 2013-05-08 深圳市晶日盛科技有限公司 一种改进的数字图像缩放方法

Also Published As

Publication number Publication date
CN104794692A (zh) 2015-07-22

Similar Documents

Publication Publication Date Title
CN104794692B (zh) 一种图像去锯齿的系统
CN102132323B (zh) 用于自动图像矫直的系统和方法
EP2962084B1 (en) Method and computer-readable medium for particle size distribution determination
CN102156972B (zh) 图像倾斜纠正方法及系统
Bertin et al. Field application of close‐range digital photogrammetry (CRDP) for grain‐scale fluvial morphology studies
CN113167569B (zh) 体积测量装置、系统、方法和程序
CN116165216B (zh) 一种液晶显示屏微划伤瑕疵3d检测方法、系统及计算设备
WO2016165097A1 (zh) 一种图像去锯齿的系统
DE60309412T2 (de) Verfahren zur Ermittlung des Verkippwinkels und der Lage eines Dokuments in einem abgescannten Bild
CN108445009B (zh) 一种太阳能电池板裂纹检测方法
CN108171674B (zh) 针对任意视角投影仪图像的视觉矫正方法
CN105678719A (zh) 一种全景拼接接缝处平滑方法及装置
CN106327576B (zh) 一种城市场景重建方法及系统
CN107729814A (zh) 一种检测车道线的方法及装置
CN105678720A (zh) 一种全景拼接判断图像匹配方法及装置
CN106204565A (zh) 一种通过斑块特征配准的图像清晰度检测装置
CN110197489B (zh) 一种无人机降落区域的选择方法和装置
CN108460348A (zh) 基于三维模型的道路目标检测方法
CN101132477A (zh) 数字图像处理方法及其装置
CN105469401B (zh) 一种基于计算机视觉的船槽定位方法
CN102737230A (zh) 一种基于方向场估计的非局部均值滤波方法
CN102073997B (zh) 小型文档图像透视复原方法
CN102831575A (zh) 图像处理装置和方法以及程序
CN106447660B (zh) 图片检测方法和装置
Gao et al. Gradient based interpolation for division of focal plane polarization imaging sensors

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate 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: 20171211

Address after: 102412 Beijing City, Fangshan District Yan Village Yan Fu Road No. 1 No. 11 building 4 layer 402

Patentee after: Beijing Si Lang science and Technology Co.,Ltd.

Address before: 100080 Zhongguancun East Road, Beijing, No. 95, No.

Patentee before: Institute of Automation, Chinese Academy of Sciences

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220119

Address after: 519031 room 532, building 18, No. 1889, Huandao East Road, Hengqin District, Zhuhai City, Guangdong Province

Patentee after: Zhuhai Jilang Semiconductor Technology Co.,Ltd.

Address before: 102412 room 402, 4th floor, building 11, No. 1, Yanfu Road, Yancun Town, Fangshan District, Beijing

Patentee before: Beijing Si Lang science and Technology Co.,Ltd.

CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: Room 701, 7th Floor, Building 56, No. 2, Jingyuan North Street, Beijing Economic and Technological Development Zone, Daxing District, Beijing 100176 (Beijing Pilot Free Trade Zone High-end Industry Zone Yizhuang Group)

Patentee after: Beijing Jilang Semiconductor Technology Co., Ltd.

Address before: 519031 room 532, building 18, No. 1889, Huandao East Road, Hengqin District, Zhuhai City, Guangdong Province

Patentee before: Zhuhai Jilang Semiconductor Technology Co.,Ltd.