CN106210448B - 一种视频图像抖动消除处理方法 - Google Patents
一种视频图像抖动消除处理方法 Download PDFInfo
- Publication number
- CN106210448B CN106210448B CN201610587481.XA CN201610587481A CN106210448B CN 106210448 B CN106210448 B CN 106210448B CN 201610587481 A CN201610587481 A CN 201610587481A CN 106210448 B CN106210448 B CN 106210448B
- Authority
- CN
- China
- Prior art keywords
- image
- block
- frame
- video
- data
- 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
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N5/00—Details of television systems
- H04N5/14—Picture signal circuitry for video frequency region
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N5/00—Details of television systems
- H04N5/14—Picture signal circuitry for video frequency region
- H04N5/21—Circuitry for suppressing or minimising disturbance, e.g. moiré or halo
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Picture Signal Circuits (AREA)
Abstract
本发明公开了一种视频图像抖动消除处理方法,其包括图像预处理、选取图像特征窗口、特征跟踪、轨迹平滑和图像置换等步骤实现图像的平稳输出,所述图像预处理通过时域滤波处理和空域滤波处理对图像进行降噪处理。本发明的视频图像抖动消除处理方法具有加图像输出平稳、图像质量高和实现性强的特点。
Description
技术领域
本发明涉及视频图像处理技术领域,具体是指一种视频图像抖动消除处理方法。
背景技术
网络视频图像是由多帧静止图像组合而来,并且静止图像间保持着某种相关性,这些相关的静止图像连接在一起,形成了具有实际意义的视频序列。消除数字视频抖动技术实质上是一种统计方法。在时间和空间上,视频序列通常包含统计冗余度。消除视频抖动技术所依赖的基本统计特性为像素之间的相关性,在各连续帧之间存在简单的相关性平移运动,消除视频抖动算法就是采用一些编码技术,以n*n像素的画面块为单位,有效的开发同一画面各附件像素之间的空间相关性。若附件帧中各像素间具有较大的相关性时,也就是说两个连续帧的内容很相似或相同时,就可以认为视频图像是相对稳定的,然而,若附近帧中各像素间具有较小的相关性时,也就是说两个连续帧的内容发生位置上的平移时,就可以认为发生了抖动,严重影响了视频的播放效果。
发明内容
本发明的目的是提供一种视频图像抖动消除处理方法,具有加图像输出平稳、图像质量高和实现性强的特点。
本发明可以通过以下技术方案来实现:
本发明公开了一种视频图像抖动消除处理方法,包括以下步骤:
选取图像特征窗口:在视频输入数据中的一帧图像选取具有变化结构模型的区间作为图的一系列像特征窗口。在视频选取特征窗口的过程,典型的特征点选择方法是基于灰度值的局部最大变化量,特征点出现在拐角处或不满足平滑约束条件的表面不连续处;拐点,角点也称为兴趣点,而特征点的跟踪就是在连续帧中兴趣点的匹配。为了对特征点进行跟踪,发现一帧图像中一系列特征窗口,在选择特征窗口过程中把那些比较重要的点形成图像的结构;由于在一帧图像中,一般只具有一个具有变化结构模型的区间,但是,同密度活线性密度的区域却很少,可以有效保证视频图像数据处理选取的准确性。
特征跟踪:在一帧图像中发现一系列特征窗口后,计算视频序列轨迹,从而在视频序列的2个基本相邻帧图像中寻找匹配点,具体算法如下:
I(x,y,t+)=I(x-Δx,y-Δy,t);
其中x,y表示图像窗口的位置。t表示某一时间捕捉到的第一帧图像;
轨迹平滑:在完成特征跟踪后,对视频序列轨迹分别在X方向和Y方向平滑原轨迹,具体平滑算法为:对X方向:x=m1t+c,对Y方向:y=m2t+c2;其中x=m1,c1,m2,c2是衡量,t表示时间。
利用特征跟踪来发现视频码流的轨迹,轨迹方向分别包括X方向和Y方向。因为任何复杂的轨迹因为X方向和Y方向的不同而分成X方向矢量和Y方向矢量。在第一帧中选择特征窗口,然后再下面N个连续图像中跟踪,在跟踪这N个连续图像后得到视频序列中不是很准确的原轨迹,因为一些干扰移动,这个轨迹有可能是混乱的,必须进行平滑化处理。
图像置换:按照如下算法对完成轨迹平滑后的轨迹进行置换处理即可实现图像的平稳输出:
其中Mk,Ni,S(k,l)表示子块中坐标为(k,l)处的像素亮度值;Sr(k+I,l+j)表示参考图像搜索窗中坐标为(k+I,l+j)处的像素亮度值。
为了产生平稳的视频,特征窗口的轨迹适应满足特征跟踪的两个公式,置换每一帧图像的结果是X方向和Y方向的干扰移动将会被抵消。但是如果被置换图像比置换图像大的多,图像的边界区域将会变成黑色。经过图像置换后,特征窗口可以比较平稳的方式运动,因此整个视频都是平稳的。
进一步地,在所述选取图像特征窗口之前,视频输入数据还包括图像预处理步骤,所述图像预处理步骤对图像实施降噪过程,排除图像干扰保留图像需要处理的部分,并过滤掉图像不需要的部分。视频图像降噪作为图像预处理技术,其表现往往决定后期的图像处理。降噪有一个固有的瓶颈,就是在降噪的同时会造成图像细节模糊与丢失,基于此,视频图像降噪包含两个方面的内容,图像噪声的滤除以及图像边缘特征的增强,有效保证了图像输出的质量。
进一步地,所述图像预处理步骤首先将待降噪的当前帧分为16×16大小的块,以块为单位进行后续处理;在划分块数据的同时将图中输入的参考帧分别置为前一帧和后一帧得到帧数据;然后将块数据和帧数据经过滤波处理然后进行平均,就能得到最终的降噪结果输出降噪的图像。
进一步地,所述滤波处理包括时域滤波和空域滤波。
进一步地,所述时域滤波处理的过程如下:首先将帧数据和块数据进行运动估计,结合MAD准则在前一帧中找到其匹配块;然后通过当前块和匹配块的共同信息检测当前运动强度系数k,最后根据运动强度k来确定当前块和匹配块时域平滑滤波强度,进行加权平均滤波。其算法数学描述如下:
Bout是当前块的降噪结果;Bt是当前块,Bt-1是运动估计在前一帧搜索到的匹配块;w是时域加权滤波系数;Th和Tl分别是根据经验设定的块运动强度高阈值和低阈值。
首先将帧数据和块数据进行运动估计,搜索匹配块;接着进行运动检测,将运动强度较小的数据先将时域加权滤波系数设定在某个范围内,然后在进行视频图像滤波时根据运动强度和噪声方差大小自适应的选取滤波系数。时域滤波算法是结合运动估计模块,根据运动强度确定当前块与匹配块时域滤波强度,从而沿着块运动轨迹进行自适应的加权平均滤波,滤波系数和阈值均随着噪声方差估计的大小自适应调节,以达到最佳的降噪效果。
进一步地,所述空域滤波处理的过程如下:将块数据进行噪声标准差估计,噪声标准差估计小的数据直接进行空域滤波处理算法处理;若噪声标准差偏大进行运动检测,选择运动强度较大的数据进行空域滤波算法处理:所述空域滤波算法如下:
NL(z)(i0,j0)=∑{i,j}∈Iω(i,j)z(i,j)
其中SAD(N(i,j),N(m,n))表示对(i,j)和(m,n)点所在的块进行SAD计,h表示系数调节,Z函数实现了归一化过程。
进一步地,所述噪声标准差估计是采用基于求块方差的加权平均方法进行的。因为噪声估计方法必须遵循的准则是尽量避免在噪声方差估计过程中受到视频图像中有用信息的影响,寻求视频图像的空白块或者较为平坦的地方。如果单纯的对选定的8×8大小宏块进行方差的计算并选出其中的最小值,可能会受到很对有用图像信息的影响,降低估计效果。因此,在对各个块进行方差计算之前先进行有用图像信息的滤除。具体实现方法为:先将8×8大小的宏块分割为四个4×4大小的子块,分别算出这四个4×4大小子块的方差,然后对着四个子块的方差进行加权平均,以加权平均后的方差值作为8*8宏块的方差估算值,最后从所有宏块中选出方差最小的一个值作为图像最终的方差估算值。
本发明一种视频图像抖动消除处理方法,具有如下的有益效果:
第一、图像输出平稳,在本发明中,消除视频图像中的抖动是利用帧与帧之间的相关性,分析出相邻帧,静止图像的抖动因子,从而消除视频图像中的抖动,其处理的过程实质是求解运动估计矢量和图像位移补偿的过程,运动估计矢量是通过估计视频图像在相邻帧之间的运动得到的,它是整个视频图像抖动消除算法的基础;抖动消除后的视频图像补偿技术由图像平移得到,并且视频图像边缘的补偿技术直接决定了抖动消除的效果;
第二、图像质量高,图像的降噪处理通过分析噪声主要来源,建立噪声模型,然后根据噪声的强弱分别进行时域滤波和空域滤波,根据环境自适应调节视频图像噪声大小,降噪后输出视频图像主观视觉质量较好,没有明显的块效应或图像质量下降;
第三、实现性强,图像抖动消除处理和图像降噪处理的步骤简单,算法简洁,处理流程化繁为简,具有较强的可实现性。
附图说明
图1为本发明一种视频图像抖动消除处理方法的流程图。
图2为本发明一种视频图像抖动消除处理方法的图像预处理的流程图。
具体实施方式
下面结合附图,对本发明的较优的实施例作进一步的详细说明。
如图1所示,本发明公开了一种视频图像抖动消除处理方法,包括以下步骤:
选取图像特征窗口:在视频输入数据中的一帧图像选取具有变化结构模型的区间作为图的一系列像特征窗口:
特征跟踪:在一帧图像中发现一系列特征窗口后,计算视频序列轨迹。具体算法如下:
I(x,y,t+)=I(x-Δx,y-Δy,t);
其中x,y表示图像窗口的位置。t表示某一时间捕捉到的第一帧图像;
轨迹平滑:在完成特征跟踪后,对视频序列轨迹分别在X方向和Y方向平滑原轨迹,具体平滑算法为:对X方向:x=m1t+c1,对Y方向:y=m2t+c2;
其中x=m1,c1,m2,c2是衡量,t表示时间;
图像置换:按照如下算法对完成轨迹平滑后的轨迹进行置换处理即可实现图像的平稳输出:
其中Mk,Ni,S(k,l)表示子块中坐标为(k,l)处的像素亮度值。Sr(k+I,l+j)表示参考图像搜索窗中坐标为(k+I,l+j)处的像素亮度值。
进一步地,在所述选取图像特征窗口之前,视频输入数据还包括图像预处理步骤,所述图像预处理步骤对图像实施降噪过程,排除图像干扰保留图像需要处理的部分,并过滤掉图像不需要的部分。
如图2所示,所述图像预处理步骤首先将待降噪的当前帧分为16×16大小的块,以块为单位进行后续处理;在划分块数据的同时将图中输入的参考帧分别置为前一帧和后一帧得到帧数据;然后将块数据和帧数据经过滤波处理然后进行平均,就能得到最终的降噪结果输出降噪的图像。所述滤波处理包括时域滤波和空域滤波。
如图2所示,所述时域滤波处理的过程如下:所述时域滤波处理的过程如下:首先将帧数据和块数据进行运动估计,结合MAD准则在前一帧中找到其匹配块;然后通过当前块和匹配块的共同信息检测当前运动强度系数k,最后根据运动强度k来确定当前块和匹配块时域平滑滤波强度,进行加权平均滤波。其算法数学描述如下:
Bout是当前块的降噪结果;Bt是当前块,Bt-1是运动估计在前一帧搜索到的匹配块;w是时域加权滤波系数。Th和Tl分别是根据经验设定的块运动强度高阈值和低阈值。
如图2所示,所述空域滤波处理的过程如下:将块数据进行噪声标准差估计,噪声标准差估计小的数据直接进行空域滤波处理算法处理;若噪声标准差偏大进行运动检测,选择运动强度较大的数据进行空域滤波算法处理:所述空域滤波算法如下:
NL(z)(i0,j0)=∑{i,j}∈Iω(i,j)z(i,j)
其中SAD(N(i,j),N(m,n))表示对(i,j)和(m,n)点所在的块进行SAD计,h表示系数调节,Z函数实现了归一化过程。所述噪声标准差估计是采用基于求块方差的加权平均方法进行的。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (4)
1.一种视频图像抖动消除处理方法,其特征在于包括以下步骤:
选取图像特征窗口:在视频输入数据中的一帧图像选取具有变化结构模型的区间作为图像的一系列特征窗口;
特征跟踪:在一帧图像中发现一系列特征窗口后,计算视频序列轨迹;具体算法如下:I(x,y,t+)=I(x-Δx,y-Δy,t);
其中x,y表示图像窗口的位置;t表示某一时间捕捉到的第一帧图像;
轨迹平滑:在完成特征跟踪后,对视频序列轨迹分别在X方向和Y方向平滑原轨迹,具体平滑算法为:对X方向:x=m1t+c1,对Y方向:y=m2t+c2;
其中m1,c1,m2,c2是恒量,t表示时间;
图像置换:按照如下算法对完成轨迹平滑后的轨迹进行置换处理即可实现图像的平稳输出:
其中S(k,l)表示子块中坐标为(k,l)处的像素亮度值;SR(k+i,l+j)表示参考图像搜索窗中坐标为(k+i,l+j)处的像素亮度值;
在所述选取图像特征窗口之前,视频输入数据还包括图像预处理步骤,所述图像预处理步骤对图像实施降噪过程,排除图像干扰保留图像需要处理的部分,并过滤掉图像不需要的部分;
所述图像预处理步骤首先将待降噪的当前帧分为16×16大小的块,以块为单位进行后续处理;在划分块数据的同时将图中输入的参考帧分别置为前一帧和后一帧得到帧数据;然后将块数据和帧数据经过滤波处理然后进行平均,就能得到最终的降噪结果输出降噪的图像;
所述滤波处理包括时域滤波和空域滤波。
4.根据权利要求3所述的视频图像抖动消除处理方法,其特征在于:所述噪声标准差估计是采用基于求块方差的加权平均方法进行的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610587481.XA CN106210448B (zh) | 2016-07-22 | 2016-07-22 | 一种视频图像抖动消除处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610587481.XA CN106210448B (zh) | 2016-07-22 | 2016-07-22 | 一种视频图像抖动消除处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106210448A CN106210448A (zh) | 2016-12-07 |
CN106210448B true CN106210448B (zh) | 2020-08-18 |
Family
ID=57491788
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610587481.XA Active CN106210448B (zh) | 2016-07-22 | 2016-07-22 | 一种视频图像抖动消除处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106210448B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3352454B1 (en) * | 2017-01-19 | 2019-01-02 | Axis AB | Method for fixed pattern noise reduction and use of such method |
CN108460731A (zh) * | 2017-03-16 | 2018-08-28 | 中国人民解放军海军航空工程学院青岛校区 | 一种消除联网飞行仿真中视景图像抖动的方法 |
CN108632501B (zh) * | 2017-03-23 | 2020-07-03 | 展讯通信(上海)有限公司 | 视频防抖方法及装置、移动终端 |
CN109859126B (zh) | 2019-01-17 | 2021-02-02 | 浙江大华技术股份有限公司 | 一种视频降噪方法、装置、电子设备及存储介质 |
CN110264431A (zh) * | 2019-06-29 | 2019-09-20 | 北京字节跳动网络技术有限公司 | 视频美化方法、装置及电子设备 |
CN110351453A (zh) * | 2019-08-16 | 2019-10-18 | 焦作大学 | 一种计算机视频数据处理方法 |
EP4113428A4 (en) * | 2020-02-25 | 2023-04-19 | Guangdong Oppo Mobile Telecommunications Corp., Ltd. | IMAGE PROCESSING METHOD AND APPARATUS AND STORAGE MEDIA |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101009021A (zh) * | 2007-01-25 | 2007-08-01 | 复旦大学 | 基于特征匹配与跟踪的视频稳定方法 |
CN101630488A (zh) * | 2009-08-05 | 2010-01-20 | 硅谷数模半导体(北京)有限公司 | 用于数字显示系统的自适应视频图像抖动方法和装置 |
CN101661611A (zh) * | 2009-09-25 | 2010-03-03 | 西安电子科技大学 | 基于贝叶斯非局部均值滤波器的实现方法 |
CN101771811A (zh) * | 2010-01-14 | 2010-07-07 | 北京大学 | 航空电子稳像器 |
CN101951464A (zh) * | 2010-05-19 | 2011-01-19 | 上海稳像信息技术有限公司 | 基于积分图特征块匹配的实时视频稳像方法 |
CN102202164A (zh) * | 2011-05-20 | 2011-09-28 | 长安大学 | 基于运动估计的道路视频稳定方法 |
CN102348046A (zh) * | 2010-07-30 | 2012-02-08 | 富士通株式会社 | 一种视频去抖动方法和视频去抖动装置 |
CN102385690A (zh) * | 2010-09-01 | 2012-03-21 | 汉王科技股份有限公司 | 基于视频图像的目标跟踪方法及系统 |
CN102999759A (zh) * | 2012-11-07 | 2013-03-27 | 东南大学 | 一种基于光流的车辆运动状态估计方法 |
CN103139568A (zh) * | 2013-02-05 | 2013-06-05 | 上海交通大学 | 基于稀疏度和保真度约束的视频稳像方法 |
CN103902972A (zh) * | 2014-03-21 | 2014-07-02 | 哈尔滨工程大学 | 一种水面移动平台视觉系统图像分析处理方法 |
CN104135597A (zh) * | 2014-07-04 | 2014-11-05 | 上海交通大学 | 一种视频抖动自动检测方法 |
-
2016
- 2016-07-22 CN CN201610587481.XA patent/CN106210448B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101009021A (zh) * | 2007-01-25 | 2007-08-01 | 复旦大学 | 基于特征匹配与跟踪的视频稳定方法 |
CN101630488A (zh) * | 2009-08-05 | 2010-01-20 | 硅谷数模半导体(北京)有限公司 | 用于数字显示系统的自适应视频图像抖动方法和装置 |
CN101661611A (zh) * | 2009-09-25 | 2010-03-03 | 西安电子科技大学 | 基于贝叶斯非局部均值滤波器的实现方法 |
CN101771811A (zh) * | 2010-01-14 | 2010-07-07 | 北京大学 | 航空电子稳像器 |
CN101951464A (zh) * | 2010-05-19 | 2011-01-19 | 上海稳像信息技术有限公司 | 基于积分图特征块匹配的实时视频稳像方法 |
CN102348046A (zh) * | 2010-07-30 | 2012-02-08 | 富士通株式会社 | 一种视频去抖动方法和视频去抖动装置 |
JP2012034361A (ja) * | 2010-07-30 | 2012-02-16 | Fujitsu Ltd | 手ぶれ補正方法及び手ぶれ補正装置 |
CN102385690A (zh) * | 2010-09-01 | 2012-03-21 | 汉王科技股份有限公司 | 基于视频图像的目标跟踪方法及系统 |
CN102202164A (zh) * | 2011-05-20 | 2011-09-28 | 长安大学 | 基于运动估计的道路视频稳定方法 |
CN102999759A (zh) * | 2012-11-07 | 2013-03-27 | 东南大学 | 一种基于光流的车辆运动状态估计方法 |
CN103139568A (zh) * | 2013-02-05 | 2013-06-05 | 上海交通大学 | 基于稀疏度和保真度约束的视频稳像方法 |
CN103902972A (zh) * | 2014-03-21 | 2014-07-02 | 哈尔滨工程大学 | 一种水面移动平台视觉系统图像分析处理方法 |
CN104135597A (zh) * | 2014-07-04 | 2014-11-05 | 上海交通大学 | 一种视频抖动自动检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106210448A (zh) | 2016-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106210448B (zh) | 一种视频图像抖动消除处理方法 | |
KR101830804B1 (ko) | 적응 필터링을 이용한 디지털 이미지 안정화 방법 | |
EP0512422A2 (en) | Method for temporally adaptive filtering of noisy image sequences | |
JP5107409B2 (ja) | 動き領域の非線形スムージングを用いた動き検出方法およびフィルタリング方法 | |
US20070070250A1 (en) | Methods for adaptive noise reduction based on global motion estimation | |
KR100985805B1 (ko) | 적응적인 칼만필터를 이용한 영상 안정화 장치 및 방법 | |
CN106331723B (zh) | 一种基于运动区域分割的视频帧率上变换方法及系统 | |
CN109963048B (zh) | 降噪方法、降噪装置及降噪电路系统 | |
Jin et al. | Quaternion-based impulse noise removal from color video sequences | |
JP2007181674A (ja) | ブロック整合及び動き補償補間を利用した映像形成方法 | |
CN111614965B (zh) | 基于图像网格光流滤波的无人机视频稳像方法及系统 | |
EP1859627B1 (en) | Global motion estimation | |
WO2014069103A1 (ja) | 画像処理装置 | |
CN108270945B (zh) | 一种运动补偿去噪方法及装置 | |
Philip et al. | A comparative study of block matching and optical flow motion estimation algorithms | |
Reeja et al. | Real time video denoising | |
JP4443767B2 (ja) | 動きベクトルの伝送コスト低減のための動き推定の方法 | |
Zhang et al. | Dehazing with improved heterogeneous atmosphere light estimation and a nonlinear color attenuation prior model | |
KR100772405B1 (ko) | 전역 움직임 추정에 근거한 적응성 잡음 감쇄 방법 및 영상처리 시스템 | |
Lu et al. | An artifact information based motion vector processing method for motion compensated frame interpolation | |
Xu et al. | Interlaced scan CCD image motion deblur for space-variant motion blurs | |
Chen et al. | Motion estimation for noise reduction in historical films: Mpeg encoding effects | |
Zlokolica et al. | Wavelet domain noise-robust motion estimation and noise estimation for video denoising | |
Talukdar et al. | A Novel Global Motion Estimation and Compensation Framework in Compressed Domain for Sign Language Videos | |
Yang et al. | An adaptive frame interpolation algorithm using statistic analysis of motions and residual energy |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20200721 Address after: Room 2b05, NANDA industry university research building, No. 6, Yuexing Third Road, high tech Zone, Yuehai street, Nanshan District, Shenzhen City, Guangdong Province Applicant after: Shenzhen smart city construction and Operation Management Co.,Ltd. Address before: 518000, Guangdong District, Longhua, Shenzhen province Longhua office, clear lake community, Silicon Valley power clear lake garden, A9 building, fourth floor Applicant before: NPE INTELLIGENT Co.,Ltd. |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant |