CN101639938B - 一种基于双边滤波器和余量去卷积的图像复原方法 - Google Patents
一种基于双边滤波器和余量去卷积的图像复原方法 Download PDFInfo
- Publication number
- CN101639938B CN101639938B CN2009101023226A CN200910102322A CN101639938B CN 101639938 B CN101639938 B CN 101639938B CN 2009101023226 A CN2009101023226 A CN 2009101023226A CN 200910102322 A CN200910102322 A CN 200910102322A CN 101639938 B CN101639938 B CN 101639938B
- Authority
- CN
- China
- Prior art keywords
- image
- restored image
- ring
- fuzzy core
- algorithm
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000001914 filtration Methods 0.000 claims abstract description 17
- 238000011084 recovery Methods 0.000 claims description 19
- 230000008030 elimination Effects 0.000 claims description 4
- 238000003379 elimination reaction Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 description 20
- 238000010586 diagram Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 6
- 238000000605 extraction Methods 0.000 description 6
- 239000012467 final product Substances 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000002146 bilateral effect Effects 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000012804 iterative process Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 241000272525 Anas platyrhynchos Species 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000033001 locomotion Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000003716 rejuvenation Effects 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 235000014347 soups Nutrition 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明公开了一种基于双边滤波器和余量去卷积的图像复原算法,(1)根据待复原图像及其已知模糊核,采用标准RL算法做高次迭代,得到含有振铃的复原图像;(2)采用双边滤波器对含有振铃的复原图像实施滤波操作,得到消除振铃后的复原图像;(3)用已知模糊核对消除振铃后的复原图像进行卷积处理,得到重模糊图像;(4)利用步骤(3)所得的重模糊图像及待复原图像采用余量去卷积得到复原图像细节,并与消除振铃后的复原图像相加得到新的含有振铃的复原图像;(5)多次重复步骤(2)~步骤(4)得到最终复原图像。本发明方法能够有效抑制迭代算法复原图像中产生的中、高频振铃,保持丰富的图像细节,得到高质量的复原图像。
Description
技术领域
本发明涉及计算机图像处理技术,尤其涉及一种图像去模糊的方法。
背景技术
在拍照过程中,由于拍照者手抖动或拍摄物体高速运动等原因,造成图像模糊,为了避免产生模糊,通常的方法是缩短曝光时间或采用三脚架拍摄,然而前者在弱光照条件下会产生严重的噪声,后者则太笨重,于是如何利用软件复原模糊图像就显得尤为重要,是现代图像处理领域的热点。图像模糊的逆过程称为图像去卷积,其目的就是从模糊图像逆向得到清晰图像,它在天文拍摄、医学成像等领域具有重要意义。
图像去卷积起源于天文观测,较早的算法包括标准RL算法、维纳滤波等,后来又出现了总体变分正则化算法以及小波域的图像复原算法等,然而即使模糊核已知,由于吉布斯现象的存在将导致振铃效应,影响复原效果,如RL算法随着迭代次数的增加图像细节逐步得到恢复,振铃效应也会越来越明显,总体变分正则算法虽然在一定程度上压制了振铃,但损失了图像细节,因此如何在保持图像细节的同时压制振铃是图像复原的关键。
模糊核对图像复原至关重要,在某些情况下模糊核是可求的,如在天文拍摄时,可将星体视为点光源,其像即为模糊核,又如有些学者提出了基于自然图像梯度分布求模糊核、采用两幅图像求模糊核以及利用图像序列求模糊核的方法,这些方法为图像复原奠定了基础。
根据算法所需信息可将图像去卷积算法分为三类,即盲复原算法、已知模糊核的复原算法和基于多幅图像的复原算法。
盲复原是在模糊核未知的情况下仅从单幅图像复原出原图像的方法,由于已知信息不足,现有的盲复原算法通常假设待求图像和模糊核的分布服从某种规律,然后在贝叶斯框架下得到复原结果,所以如何对图像和模糊核的分布做出合理假设就成为图像复原的关键。近年来这方面的研究取得了一定进展,如先假设待求图像服从自然图像梯度分布规律,并在尺度空间上求得模糊核,然后用标准RL算法得到复原结果,该方法较传统算法的优势在于能够处理模糊核较复杂的情况;又如在最大后验框架下采用对图像噪声建模和引入局部约束压制振铃的方法可得到高质量的复原图像,优于传统的盲复原算法。
已知模糊核的复原算法又称为非盲复原算法,由于模糊核已知,求解相对简单,所以涌现出了众多算法,其中一些经典算法一直沿用至今,如RL、维纳滤波算法等,它们在运行效率,图像恢复质量上都有很好的表现,然而吉布斯现象导致的振铃效应一直成为困扰这些经典算法的难题,为了降低振铃效应的影响,科研工作者提出了许多新算法,如总体变分正则化算法,将稀疏自然图像约束与共轭梯度法相结合的算法,根据联合双边/双边滤波器的思想对标准RL算法进行改进,并在尺度空间上复原图像的算法等。此外,随着小波技术的发展,还出现了许多基于小波域的图像复原技术。
由于在实际工作中单幅图像包含的信息有限,所以要获得准确的模糊核并非易事,于是又出现了基于多幅图像的复原算法,如基于欠曝图像和模糊图像的复原算法;基于两幅模糊图像的复原算法以及基于序列图像的复原算法。
发明内容
本发明提出了一种基于双边滤波器和余量去卷积的图像复原算法,结合了标准RL算法、双边滤波器和余量去卷积算法的特点,通过多次滤波、重模糊和细节复原得到清晰图像,该算法能有效抑制迭代过程中产生的中、高频振铃,保持丰富的图像细节,步骤如下:
(1)根据待复原图像及其已知模糊核,采用标准RL算法做高次迭代,得到含有振铃的复原图像;
(2)采用双边滤波器对含有振铃的复原图像实施滤波操作,得到消除振铃后的复原图像;
(3)用步骤(1)所述的已知模糊核对消除振铃后的复原图像进行卷积处理,得到重模糊图像;
(4)利用步骤(3)所得的重模糊图像及待复原图像采用余量去卷积得到复原图像细节,将复原图像细节与步骤(2)的消除振铃后的复原图像相加得到新的含有振铃的复原图像,并将其作为步骤(2)中双边滤波器的输入图像;
(5)多次重复步骤(2)~步骤(4)得到最终复原图像。
图像模糊过程可用如下的数学模型表述:
在已知模糊核的图像复原算法中B和h是已知的,I是待求项,解决该问题的方法有很多种,例如共轭梯度法、标准RL算法和维纳滤波等,然而由于吉布斯现象的存在,所有这些算法都会在复原结果中引入振铃现象,振铃的频率不仅取决于模糊核和图像的相对大小,而且与模糊核的强度分布有重要关系,低频振铃会在图像背景中引入较宽的波纹,但不会对图像细节辨识造成太大影响,而中、高频的振铃将影响到图像细节,产生伪信息,如何在保持图像细节的情况下压制振铃是一个难点。
本发明是一种已知模糊核的图像复原算法,模糊核的获得有多种方法,如可根据镜头的成像特性推导出模糊核、根据序列图像求得模糊核等等,在此认为模糊核是已知的。模糊核已知后就可以利用本发明方法复原图像。
本发明方法由如下三部分核心步骤构成:
标准RL算法。
标准RL算法是一种经典的图像去卷积算法,它假设图像噪声服从泊淞分布规律,由此推导出图像概率表达式,然后再求出该表达式的最大值,即为复原结果,迭代过程如下:
其中,h*(i,j)=h(-j,-i),i,j表示以模糊核中心为坐标原点,模糊核中各像素的坐标值,t表示迭代次数,It表示图像复原过程中经过t次迭代后得到的复原图像。标准RL算法的特点在于它具有较强的细节复原能力,迭代次数越高,细节越丰富,同时振铃的影响也越来越严重,特别是当迭代次数高于某一阈值时,图像细节将基本保持不变,而振铃效应将愈加显著。
双边滤波器。
双边滤波器是一种非线性滤波器,具有很强的细节提取能力,它首先为含有振铃的复原图像中的每个像素规定一个邻域,滤波操作被限定在该邻域中,滤波结果不仅取决于邻域中心像素与相邻像素间的相对距离,还取决于它们的相对强度,其表达式如下:
其中,Ap f表示对含有振铃的复原图像实施滤波后所得图像中某点的像素值,利用上式得到图像中所有点的像素值后即得到了消除振铃后的复原图像; 为归一化系数;W表示含有振铃的复原图像中某个待处理像素的邻域,p表示该邻域中心坐标,p′表示邻域中其他像素的坐标,Ap表示该邻域中心像素值,Ap′表示邻域中其他像素的像素值,Gd和Gr都为高斯函数,其标准差分别为σd和σr,σd同时还限定了邻域的大小。从上式可以看出,双边滤波器首先为邻域W中的每个像素Ap′乘以一个权重Gd(p′-p)Gr(Ap-Ap′),然后求得加权和 再做归一化处理即乘以得到待处理像素滤波后的值,高斯函数Gd使得距离中心像素越远的像素对滤波结果的影响越小,同理Gr则使得邻域中与中心像素强度相差越大的点对结果的影响越小。双边滤波器具有很强的细节提取能力,只要调整σd和σr的值,便可以得到各种程度不同的滤波效果。
由于本发明方法中需要多次重复步骤(2)~步骤(4),随着重复次数的增加,复原图像中的振铃效应将越来越小,而双边滤波器的作用在于为下面所述的余量去卷积提供平滑图像,所以双边滤波器参数应随步骤(2)~步骤(4)重复次数的增加做适当减小,以减少细节提取量,只需保证每次滤波恰好除去输入图像中的振铃即可。
余量去卷积。
研究表明,由吉布斯现象产生的振铃效应的幅度正比于图像中的像素突变值,因此,在图像边缘振铃效应较明显,为了削弱振铃,首先在步骤(2)中用双边滤波器对步骤(1)或(4)中所得的含有振铃的复原图像进行滤波操作得到消除振铃后的复原图像NF,则剩余模糊图像 其中即为步骤(3)所述的重模糊图像。
由于ΔB中可能含有负值,应先为其加上偏置1,再使用标准RL算法进行迭代得到复原图像细节如下:
其中t表示迭代次数,ΔIt表示复原图像细节,由于该算法是在剩余模糊图上实施标准RL算法,所以突变点处的像素值已被削弱,因此它能够有效抑制振铃效应。
由于本发明方法中需要多次重复步骤(2)~步骤(4)才得到振铃较少、细节丰富的最终复原图像。在此过程中,步骤(4)中余量去卷积迭代次数将随步骤(2)~步骤(4)重复操作次数的增加做适当减少,但应高于某一阈值,这有利于抑制振铃效应。
本发明方法步骤(1)中利用了所述的标准RL算法的特点,通过高次迭代,获得了细节丰富的复原图像,同时结果中的振铃效应也较明显。
本发明方法步骤(2)~(4)中则利用了双边滤波器在提取细节方面的优势,通过参数调整,刚好去掉输入图像中的振铃,得到平滑图像NF,然后再将NF与模糊核h做卷积得到重模糊图像它与原模糊图的差值图像 即为剩余模糊图像,再采用余量去卷积进行高次迭代即可得到复原图像细节ΔIt,则通过步骤(2)~(4)的操作得到新的复原图像I=NF+ΔIt,此时余量去卷积算法已对振铃实施了明显压制,且图像细节得到了保留,为了进一步消除残余振铃,只需多次重复步骤(2)~(4)的操作即可。
本发明方法能够有效抑制迭代算法复原图像中产生的中、高频振铃,保持丰富的图像细节,得到高质量的复原图像。适用于卫星遥感、飞机航拍和医学图像处理等领域。
附图说明
图1(a)为处理前的模糊图像;
图1(b)为模糊核;
图2为图1(a)经过本发明步骤(1)处理后的效果示意图;
图3为图2经过本发明步骤(2)处理后的效果示意图;
图4为图3经过本发明步骤(3)处理后的效果示意图;
图5为图4经过本发明步骤(4)处理后的效果示意图;
图6为图5再经过4次本发明步骤(2)~步骤(4)的重复处理后的效果示意图;
图7为本发明方法流程图。
具体实施方式
如图1所示,本发明方法的已知信息为待复原模糊图像B和模糊核h,其中模糊图像B大小为512×512,模糊核h为高斯型,大小为10×10;首先根据模糊图像B和模糊核采用标准RL算法做高次迭代,得到含有振铃的复原图像。
迭代的过程如下:
其中:
h*(i,j)=h(-j,-i),i,j表示以模糊核中心为坐标原点,模糊核中各像素的坐标值;
t表示迭代次数;
It表示图像复原过程中经过t次迭代后得到的复原图像;
B表示待复原图像;
h表示待复原图像的模糊核;
鉴于标准RL算法的特点,该步骤中的迭代次数t应高于某一阈值(本实验为150次),这样做的目的是使复原结果中含有尽可能多的图像细节,当然该复原图像将受到振铃现象的影响,为了便于叙述,将该步骤称为“预操作”,图2为“预操作”结果。
然后采用双边滤波器对含有振铃的复原图像进行细节提取,双边滤波器有很强的细节提取能力,只要调节其参数就能得到不同程度的平滑图像NF,这里只需保证“预操作”结果中的振铃被刚好去掉即可,方法如下式:
其中:
Ap f表示对含有振铃的复原图像实施滤波后所得图像中某点的像素值;
W表示含有振铃的复原图像中某个待处理像素的邻域,p表示该邻域中心坐标,p′表示邻域中其他像素的坐标,Ap表示该邻域中心像素值,Ap′表示邻域中其他像素的像素值;
Gd和Gr都为高斯函数,其标准差分别为σd和σr,(σd=3,σr=0.3)。
结果如图3所示。
图4是对图3进行重模糊的结果,即用模糊核h与平滑图像NF进行卷积操作即可,然后求出剩余模糊图像 接着采用余量去卷积算法进行迭代得到图像细节ΔIt,
其中:
t表示迭代次数;
ΔIt表示复原图像细节;
B表示待复原图像;
h表示待复原图像的模糊核;
在余量去卷积过程中,同样应保证迭代次数t高于某一阈值(本实验t为50次)以便得到丰富的细节,则得到新的复原图像为I=NF+ΔIt,此时该复原图像中的振铃已得到明显抑制,如图5所示。
由于在标准RL算法中,由吉布斯现象产生的振铃效应的幅度正比于图像中的像素突变值,因此在图像的边缘点处将产生明显的振铃现象,而余量去卷积是在剩余图像上实施标准RL算法,原图中的像素突变值已被削弱,所以结果中的振铃被明显抑制,又因为迭代次数较高,所以图像细节得以保留。为了便于叙述,将上述双边滤波、平滑图像重模糊和余量去卷积合称为“操作”。
为了进一步削弱复原图像中的残余振铃,应对新的复原图像做多次“操作”,在多次“操作”过程中,应使双边滤波器的参数逐次减小,只要保证在每次操作中振铃被刚好消除即可,同时余量去卷积的次数也应逐次降低,但应高于某一阈值(本实验为50次),这样既有利于抑制振铃同时又能保持复原图像中的细节,图6是经过共五次“操作”所得的结果,明显优于图2。
本发明方法能够有效抑制迭代算法复原图像中产生的中、高频振铃,保持丰富的图像细节,得到高质量的复原图像。
Claims (3)
1.一种基于双边滤波器和余量去卷积的图像复原方法,其特征在于,包括如下步骤:
(1)根据待复原图像及其已知模糊核,采用标准RL算法做高次迭代,得到含有振铃的复原图像;
(2)采用双边滤波器对含有振铃的复原图像实施滤波操作,得到消除振铃后的复原图像;
(3)用步骤(1)所述的已知模糊核对消除振铃后的复原图像进行卷积处理,得到重模糊图像;
(4)利用步骤(3)所得的重模糊图像及待复原图像采用余量去卷积得到复原图像细节,将复原图像细节与步骤(2)的消除振铃后的复原图像相加得到新的含有振铃的复原图像,并将其作为步骤(2)中双边滤波器的输入图像;所述的利用步骤(3)所得的重模糊图像及待复原图像采用余量去卷积得到复原图像细节的方法如下式:
其中:
t表示迭代次数;
ΔIt表示复原图像细节;
B表示待复原图像;
h表示待复原图像的模糊核;
表示卷积运算;
h*(i,j)=h(-j,-i),i,j表示以模糊核中心为坐标原点,模糊核中各像素的坐标值;
(5)多次重复步骤(2)~步骤(4)得到最终复原图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101023226A CN101639938B (zh) | 2009-08-28 | 2009-08-28 | 一种基于双边滤波器和余量去卷积的图像复原方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101023226A CN101639938B (zh) | 2009-08-28 | 2009-08-28 | 一种基于双边滤波器和余量去卷积的图像复原方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101639938A CN101639938A (zh) | 2010-02-03 |
CN101639938B true CN101639938B (zh) | 2011-10-19 |
Family
ID=41614907
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009101023226A Expired - Fee Related CN101639938B (zh) | 2009-08-28 | 2009-08-28 | 一种基于双边滤波器和余量去卷积的图像复原方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101639938B (zh) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102542533A (zh) * | 2010-12-16 | 2012-07-04 | 方正国际软件(北京)有限公司 | 一种图像模糊方法及*** |
CN102131079B (zh) * | 2011-04-20 | 2012-08-22 | 浙江宇视科技有限公司 | 一种消除图像运动模糊的方法和装置 |
CN102142134A (zh) * | 2011-04-29 | 2011-08-03 | 中山大学 | 一种基于三维栅格平滑模型的图像细节增强方法 |
CN102147915B (zh) * | 2011-05-06 | 2012-06-27 | 重庆大学 | 一种权重的稀疏边缘正则化图像复原方法 |
CN102903077A (zh) * | 2011-07-26 | 2013-01-30 | 陈奋 | 一种快速图像去模糊算法 |
CN105096261B (zh) * | 2014-05-13 | 2018-04-17 | 北京大学 | 图像处理装置和图像处理方法 |
CN103971345B (zh) * | 2014-05-27 | 2018-03-02 | 上海电力学院 | 一种基于改进双边滤波的图像去噪方法 |
CN104008531B (zh) * | 2014-06-17 | 2017-02-08 | 中国电子科技集团公司第二十八研究所 | 一种基于混合型马尔科夫专家场的模糊图像盲复原方法 |
CN104933693A (zh) * | 2015-07-02 | 2015-09-23 | 浙江大学 | 一种用于有饱和像素场景的多幅图像复原方法 |
CN104966277B (zh) * | 2015-07-10 | 2017-11-07 | 浙江大学 | 一种基于长短曝光图像对的图像复原方法 |
CN106530381B (zh) * | 2016-10-19 | 2019-01-29 | 浙江大学 | 一种基于gpu加速的三维荧光显微图像的去卷积算法 |
CN106600543A (zh) * | 2016-11-10 | 2017-04-26 | 江苏信息职业技术学院 | 基于粒子群优化的图像复原方法 |
CN107123098A (zh) * | 2017-04-26 | 2017-09-01 | 航天科技控股集团股份有限公司 | 一种用于解决图像重影及折射不定状态的仪表图像处理方法及装置 |
CN107124531A (zh) * | 2017-05-26 | 2017-09-01 | 努比亚技术有限公司 | 一种图像处理方法及移动终端 |
CN107945127B (zh) * | 2017-11-27 | 2021-05-07 | 南昌大学 | 一种基于图像列灰度概率一致性的高速运动图像去模糊方法 |
CN108447041B (zh) * | 2018-01-30 | 2020-12-15 | 中国航天电子技术研究院 | 一种基于增强学习的多源图像融合方法 |
CN108510459A (zh) * | 2018-04-08 | 2018-09-07 | 哈尔滨理工大学 | 一种基于小波自适应阈值与双边滤波图像降噪算法 |
CN109029425B (zh) * | 2018-06-25 | 2020-07-31 | 中国科学院长春光学精密机械与物理研究所 | 一种采用区域滤波的模糊星图复原方法 |
CN110717870B (zh) * | 2019-06-29 | 2022-04-12 | 浙江大学 | Usb表面缺陷检测的图像预处理方法 |
CN112116544B (zh) * | 2020-10-08 | 2022-04-22 | 东南数字经济发展研究院 | 一种抵抗图像平滑滤波的预处理方法 |
CN112614080A (zh) * | 2021-01-06 | 2021-04-06 | 中国科学院微小卫星创新研究院 | 基于遥感图像序列的去模糊方法 |
CN113592728B (zh) * | 2021-07-01 | 2024-04-05 | 温州理工学院 | 一种图像复原方法、***、处理终端及计算机介质 |
-
2009
- 2009-08-28 CN CN2009101023226A patent/CN101639938B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101639938A (zh) | 2010-02-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101639938B (zh) | 一种基于双边滤波器和余量去卷积的图像复原方法 | |
Zheng et al. | Adaptive unfolding total variation network for low-light image enhancement | |
Dudhane et al. | Burst image restoration and enhancement | |
CN111028177B (zh) | 一种基于边缘的深度学习图像去运动模糊方法 | |
US10755390B2 (en) | Image deblurring method based on light streak information in an image | |
CN108550118B (zh) | 运动模糊图像的模糊处理方法、装置、设备及存储介质 | |
Lai et al. | Blur kernel estimation using normalized color-line prior | |
CN101930601B (zh) | 一种基于边缘信息的多尺度模糊图像盲复原方法 | |
CN103761710B (zh) | 基于边缘自适应的高效图像盲去模糊方法 | |
CN104599242B (zh) | 使用多尺度非局部正则的模糊核估计方法 | |
CN104091312B (zh) | 一种根据图像频谱信息提取模糊核先验的单透镜成像方法 | |
CN102254309A (zh) | 一种基于近红外图像的运动模糊图像去模糊方法和装置 | |
Schelten et al. | Interleaved regression tree field cascades for blind image deconvolution | |
CN107730466A (zh) | 基于结构自相似性和低秩矩阵表示的图像盲复原方法 | |
Zhang et al. | Deep motion blur removal using noisy/blurry image pairs | |
Wang et al. | An improved image blind deblurring based on dark channel prior | |
CN110852947B (zh) | 一种基于边缘锐化的红外图像超分辨方法 | |
CN111626943B (zh) | 基于一阶前向后向算法的全变分图像去噪方法 | |
Guo et al. | Image blind deblurring using an adaptive patch prior | |
Senshiki et al. | PSF estimation using total variation regularization and shock filter for blind deconvolution | |
US20240135496A1 (en) | System and method for burst image restoration and enhancement | |
CN103337055A (zh) | 一种基于梯度拟合的文本图像去模糊方法 | |
Esmaeilzehi et al. | DSegAN: A deep light-weight segmentation-based attention network for image restoration | |
Zhang et al. | Parametric model for image blur kernel estimation | |
Wang et al. | Efficient image deblurring via blockwise non-blind deconvolution algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20111019 Termination date: 20130828 |