CN103076552B - 一种抑制偏置电压纹波和列输出放大器通道不均匀的方法 - Google Patents
一种抑制偏置电压纹波和列输出放大器通道不均匀的方法 Download PDFInfo
- Publication number
- CN103076552B CN103076552B CN201210080080.7A CN201210080080A CN103076552B CN 103076552 B CN103076552 B CN 103076552B CN 201210080080 A CN201210080080 A CN 201210080080A CN 103076552 B CN103076552 B CN 103076552B
- Authority
- CN
- China
- Prior art keywords
- signal
- row
- focal plane
- integration
- detector
- 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
Abstract
本发明公开了一种抑制偏置电压纹波和列输出放大器通道不均匀的方法,外界辐射入射到改进的探测器焦平面上由焦平面的光电转换效应将光能转换成电流的形式,即在标准的 象元的探测器焦平面阵列基础上做的列延拓,位于象元阵列周围,使原有的阵列大小变为,将蓝宝石作为列延拓背景象元的材料,整个象元阵列与读出电路连接;读出电路按照时序要求对每一帧的图像按行积分,对电压信号进行处理,使得最后的实际输出都保持在一个均匀的程度,从而有效的去除图像上的滚动条纹。
Description
技术领域
本发明属于微测辐射热计读出电路硬件定标技术,特别是一种抑制偏置电压纹波和列输出放大器通道不均匀的方法。
背景技术
红外辐射能够被探测是由于物质内部的相互作用,通过光子间交换的加热效应,或是通过直接电子吸收/内部光电子效应。前者,吸收辐射而加热,使温度升高,因此,测量出探测器随温度变化的参数,这样的器件称为平方律器件;后者,带隙半导体中的价带电子由于直接吸取光子而跃迁到导带上,因此能够检测出导带内电子数的增加,这类器件称为量子探测器,平方律和量子探测器之间最突出的差别就是后者必须制冷到低温才能工作。
对于不需要进行制冷的微测辐射热计,其工作状态受环境温度影响特别大,表现在整机结构中为各电路板上的IC器件的工作状态会随着外界环境温度的变化而发生些许的改变,特别是一些分立元件的电学特性,导致部分提供给探测器工作所需的偏置电压在一个平均值范围内发生上下微小的波动,形成一定的噪声,即如下式所示:
(1)
其中,为随时间变化的偏置电压实际值,为理论上应当输出的偏置电压标准值,为受环境温度影响的电压波动噪声。这种波动噪声会影响到探测器进行图像行积分时所采到的有效信号的大小,导致即使是观察均匀的背景,其图像上每行之间也存在明显的差异,这种差异就以条纹的形式表现出来,又因为是行积分,所以条纹会在图像上不停的滚动,大大降低了图像质量及信噪比(SNR)。
目前应用最广泛的去除这种由于偏置电压波动而引起的条纹噪声的办法是对于得到的每一行的数字信号进行多帧累加取均值的算法,即连续采集多帧内各行的数据并分别取其平均值,然后用各行数据减去该平均值的方法来抑制条纹的形成,这样做的结果虽然能够在一定程度上抑制条纹,但是对于外界照度很强时各行的均值差异较大,很多时候并不能完全去除该种条纹噪声。
发明内容
本发明的目的在于提供一种抑制偏置电压纹波和列输出放大器通道不均匀的方法,针对外界输出微测辐射热计焦平面的偏置电压上存在随时间变化波动的噪声而导致的图像行积分采样时列通道输出的不均匀带来的滚动条纹问题,能够有效的去除图像上的滚动条纹。
实现本发明目的的技术解决方案为:一种抑制偏置电压纹波和列输出放大器通道不均匀的方法,步骤如下:
第一步,外界辐射入射到改进的探测器焦平面上由焦平面的光电转换效应将光能转换成电流的形式,即在标准的象元的探测器焦平面阵列基础上做的列延拓,位于象元阵列周围,使原有的阵列大小变为,将蓝宝石作为列延拓背景象元的材料,整个象元阵列与读出电路连接,的取值是根据现有的焦平面来选取,其中为整数,N最大为1024,M最大为1024,n小于N;
第二步,读出电路按照时序要求对每一帧的图像按行积分:读出电路是将探测器阵列所采集的电信号,按时序输出,并对输出信号进行处理,输出N*M像素的图像;
第三步,对电压信号进行处理:
(1)读出电路计算延拓出来的n列背景象元上当前的入射辐射所引起的温度—电压信号的平均值;
(2)对一行上的所有N个像素积分完成后,将每个像素的积分值减去n列背景象元的电压平均值;
(3)通过探测器焦平面输出端口上的采样电路对减去“背景象元”平均值的信号进行采样并输出即得到了均匀的每一行像素的电压值,从而达到消除条纹的目标。
本发明与现有技术相比,其显著优点:通过采用对环境辐射温度不敏感的蓝宝石材料制作的背景象元来对读出电路做n元列延拓,通过对n列背景象元求和取均值的方式计算出标准值,并在输出前将各采样值与标准背景值取差,使得最后的实际输出都保持在一个均匀的程度,从而有效的去除图像上的滚动条纹。
下面结合附图对本发明作进一步详细描述。
附图说明
图1是探测器焦平面结构图。
图2是单个象元结构图。
图3是读出电路结构简图。
图4 是读出电路时序要求。
图5 是输出信号的采样波形。
图6 是元延拓方式。
图7是偏置电压噪声引起的含有滚动条纹的图像。
图8 是列输出不均匀性被抑制的无滚动条纹的图像。
具体实施方式
本发明抑制偏置电压纹波和列输出放大器通道不均匀的方法,步骤如下:
第一步,外界辐射入射到改进的探测器焦平面上由焦平面的光电转换效应将光能转换成电流的形式,即在标准的象元的探测器焦平面阵列基础上做的列延拓,位于象元阵列周围,使原有的阵列大小变为,将蓝宝石作为列延拓背景象元的材料,整个象元阵列与读出电路连接;的取值是根据现有的焦平面来选取,其中N、M、n为整数,N最大为1024,M最大为1024,n一般小于N,取其二分之一较为合适。下面内容,取N为320,M为240,n为160为例。
第二步,读出电路按照时序要求对每一帧的图像按行积分,
(1)读出电路:读出电路是将探测器阵列所采集的电信号,按如图4的时序输出,并对输出信号进行处理,读出电路包含与红外探测器阵列相应的单元传输电流阵列,每个单元传输电流阵列由偏置电路10、水平增益电路11、积分电路12、采样电路13、多路开关14依次连接组成,电流输入端9与每个单元传输电流阵列连接,电流信号从电流输入端9处输入,依次经过240个并联的偏置电路10进行调整,然后由水平增益电路11调整,再按行信号积分电路12进行积分,由采样电路13进行采样和多路开关14的选择,输出N*M像素的图像,如图3所示。
(2)时序要求,其中,MC为探测器主时钟,RESET为复位信号,INT为积分信号,SORTIE为视频信号,DATAVALID是数据有效信号,LIGNE1是第一行读出信号。
探测器主时钟MC的最高频率为5.5 MHz和50%的占空比,是整个电路的操作同步的基础。像素寻址是由主时钟控制的静态同步移位寄存器执行的。 像素信息是在每一个主时钟周期率的读出。
积分信号INT允许来自一行的焦平面的信号积分,当它是高水平的时候,于纹波操作模式,在INT阶段必须在每个发送行积分,每帧240次。
积分信号INT在主时钟的上升沿必须改变状态。它的周期必须大于340主时钟周期。积分时间等于INT为高水平的持续时间减去2个主时钟周期。最低积分时间是主时钟周期。最大积分时间为318个主时钟周期。
复位信号RESET,通过迫使焦平面第一行的信号积分,使读出电路的操作复位。在主时钟的上升沿过程中,复位信号必须改变其状态。每帧不得重复一次以上。
数据有效信号DATAVALID,为了获得读出信息,是在积分结束后自动设置到一个高电平。
第一行读出信号LIGNE1,自动设置在高电平,以获得来自第一行的FPA的信号输出。
第三步,对电压信号进行处理:
(1)读出电路计算延拓出来的n列背景象元上当前的入射辐射所引起的温度—电压信号的平均值;
(2)对一行上的所有N个像素积分完成后,将每个像素的积分值减去n列背景象元的电压平均值;
(3)通过探测器焦平面输出端口上的采样电路对减去“背景象元”平均值的信号进行采样并输出即得到了均匀的每一行像素的电压值,从而达到消除条纹的目标。
本发明从读出电路结构上来对探测器图像行积分时由于工作偏置电压的波动带来的条纹噪声进行抑制的方法,原理清晰,结构简单,且不占用后端的资源以及功耗。红外辐射的光能要转变为电能继而以电压的形式输出,并能被后端的处理电路接收并进行处理需要经过光电转换—电流电压转换两个过程,对于探测器来说,这两个过程分别发生在探测器内部的焦平面和读出电路中,焦平面接收外界辐射的光能,将其转变为电流的形式,并通过读出电路的积分过程将电流转变为电压形式输出,附图1、2为探测器焦平面及读出电路的结构图。
图1所示为整个焦平面的结构分布,红外辐射1辐照在焦平面2上,焦平面2通过铟柱3与下面的读出电路4连接,将吸收到的红外辐射转变为电流后通过读出电路4进行信号输出5。
图2是单一象元的结构图,对于个象元的焦平面阵列,就是由个图2中示意的结构组成了整个焦平面。红外辐射1透过衬底6依次辐照在焦平面N沟道7和焦平面P沟道8上,焦平面P沟道8通过铟柱3与下面的读出电路4连接,将吸收到的红外辐射转变为电流后通过读出电路4进行信号输出5。
本发明的方法是针对象元的微测辐射热计焦平面阵列做列延拓,用受辐射照射温度变化十分小的半导体材料作为“背景象元”,该背景象元搭载在读出电路上,位于标准的大小阵列周围,使原有的阵列大小变为,用图6表示。
本发明的工作原理是:首先分析了该条纹产生的原因,是由于每个像素工作时偏置电压上微小的波动导致采样值无法保持相等,而是在一个统计平均值的周围上下波动,从图5中得到:
(2)
其中为随时间变化的偏置电压实际值,为理论上应当输出的偏置电压标准值,为受环境温度影响的电压波动噪声。由2式可以得到对于一行像素,其总的输出响应Vout应该满足下式:
(3)
由于每行像素的统计平均值有所差别,导致显示在监视器上时出现滚动条纹,因此,从结构的角度分析,只要在每行像素积分完成后能够减去同一个背景值,那么就相当于将噪声的波形扰动直流分量去除掉,每个像素的实际值就能保持在理论范围的同一个均匀水平上,该统计平均值Vave可以用下式得到:
(4)
在图5中,令的标准像素值为,的背景象元的值为,则读出电路的总像素值V可用下式表示:
(5)
在完成了光电转换过程后,读出电路根据外部配置的读出时序的要求开始对外进行信号输出。读出电路按照行的方向逐行的对信号电流进行积分和采样,然后输出有效的图像信号,其读出结构简图如附图3所示。从图中可以看出偏置电压FID和EB直接与探测器的输出动态范围和直流电平有关,也就是说,这两个偏置电压对于条纹噪声的贡献最大,下面图4为单一象元在提供偏置电压的情况下的输出信号的采样情况图。从图3,4中可以对这种噪声波动引起的条纹做解释,当每行的电流信号被积分电容转变为电压信号后,由于偏置电压的噪声影响,其输出波形并不是各个时间相同的幅值,而是在一个平均值周围上下波动的波形,经由后面的采样输出放大器采样后就会得到一些同样在一个均值附近上下跳动的值,在一帧图像内,各行的均值都不相同的情况下就会产生这种条纹。
假设不考虑读出电路的KTC噪声以及1/f噪声,=,因此只要将式3,4两式相减,然后再将这个值输出就能够得到各行均匀的探测器信号,就能从根本上去除滚动条纹。
本发明的具体实现需要利用到两方面的控制,一是探测器内部读出电路上结构延拓;二是需要用后端的FPGA对其“背景象元”的求和取平均操作进行数字控制。要在探测器正常的行积分时间内完成对“背景象元”的求和取平均过程,并在一行像素积分完成的时候将其有效信号值与“背景象元”的平均值取差,然后再经过采样的过程对外输出。
“背景象元”需要采用对于外界辐射温度—电压变化不敏感的半导体材料,因为该延拓是直接连接在读出电路周围的,所以也会受到部分焦平面上外界辐射的影响,在选取材料时还需要特别注意,被选择的材料要与探测器象元的化学性质相配合,特别是与衬底的化学性质的配合,因此,本发明使用蓝宝石材料为背景象元的制作材料,蓝宝石材料目前已经被广泛的使用与制作或者经由掺杂后制作制冷型探测器焦平面的衬底材料,其受温度的影响非常小,可以保证外界辐射无论如何变化,其性质都保持稳定,即可以为后面的n列延拓的求和取平均提供十分稳定的输出值。
此外,读出电路的结构相比普通的微测辐射热计结构需要有所调整,在多出的n行,m列背景象元上需要有CMOS结构的时序控制设计,以便让后端的FPGA能够实现对前端的时序操作及控制。在图1中的信号输出接插件中需要有能够和外部的FPGA连接的背景象元时序控制端口,后端的FPGA除了提供探测器正常工作所需要的数字信号,包括时钟信号,积分信号等以外,还需要提供背景象元控制信号,该信号与行积分信号后端同步,即当行积分信号结束时,正常的N个像素的值已经得到的情况下,延拓后的n个背景象元的平均值也已经计算完毕,此时经过求差操作就可以得到均匀的信号值。
背景象元的求和方式等同于普通象元的积分累加方式,即,在提供一个工作电压和数字时序的情况下,每一行后面延拓的背景象元的值通过累加的方式合到一起,然后通过放大器的除法电路作用,将计算出来的电压值存在电容中,等到有效像元的积分完成后将电容中的电压和信号电压通过专门的减法电路取差,最后由连接外部输出端口的输出采样放大电路完成输出的工作,整个设计原理清晰,实现方法容易实现,由于是硬件实现方法,能够从根本上改善输出图像的质量,提高图像信噪比,相比于繁琐的软件算法,本发明不仅快速,而且能节省大量的系统资源,降低系统功耗,适合于长时间使用电池工作的方式,是一项十分具有工程意义的发明项目。
图6、图7为含有严重条纹的图像及使用该方法后的图像,通过对比可以发现,条纹已经几乎全部去除掉,证明了该方法对抑制由偏置电压噪声带来的列输出不均匀问题起到了很好的效果。
Claims (3)
1.一种抑制偏置电压纹波和列输出放大器通道不均匀的方法,其特征在于步骤如下:
第一步,外界辐射入射到改进的探测器焦平面上由焦平面的光电转换效应将光能转换成电流的形式,即在标准的N×M象元的探测器焦平面阵列基础上做n的列延拓,位于N×M象元阵列周围,使原有的阵列大小N×M变为(N+n)×M,将蓝宝石作为n列延拓背景象元的材料,整个(N+n)×M象元阵列与读出电路连接,N、M、n的取值是根据现有的焦平面来选取,其中N、M、n为整数,N最大为1024,M最大为1024,n小于N;
第二步,读出电路按照时序要求对每一帧的图像按行积分:读出电路是将探测器阵列所采集的电信号,按时序输出,并对输出信号进行处理,输出N×M像素的图像;
第三步,对电压信号进行处理:
(1)读出电路计算延拓出来的n列背景象元上当前的入射辐射所引起的温度—电压信号的平均值;
(2)对一行上的所有N个像素积分完成后,将每个像素的积分值减去n列背景象元的温度—电压信号的平均值;
(3)通过探测器焦平面输出端口上的采样电路对减去“背景象元”温度—电压信号的平均值的信号进行采样并输出即得到了均匀的每一行像素的电压值,从而达到消除条纹的目标。
2.根据权利要求1所述的抑制偏置电压纹波和列输出放大器通道不均匀的方法,其特征在于第二步中的读出电路包含与探测器阵列相应的单元传输电流阵列,电流输入端(9)与每个单元传输电流阵列连接,电流信号从电流输入端(9)处输入,依次经过240个并联的偏置电路(10)进行调整,然后由水平增益电路(11)调整,再按行信号积分电路(12)进行积分,由采样电路(13)进行采样和多路开关(14)的选择,输出N×M像素的图像。
3.根据权利要求1所述的抑制偏置电压纹波和列输出放大器通道不均匀的方法,其特征在于第二步的时序要求中,MC为探测器主时钟,RESET为复位信号,INT为积分信号,SORTIE为视频信号,DATAVALID是数据有效信号,LIGNE1是第一行读出信号:探测器主时钟MC的最高频率为5.5MHz和50%的占空比,是整个电路的操作同步的基础,像素寻址是由主时钟控制的静态同步移位寄存器执行的,像素信息是在每一个主时钟周期率的读出;
积分信号INT允许来自一行的焦平面的信号积分,当积分信号INT是高电平的时候,在INT阶段必须在逐个发送行积分,每帧240次;
积分信号INT在探测器主时钟MC的上升沿改变状态,积分信号INT的周期大于340个主时钟周期,积分时间等于INT为高电平的持续时间减去2个主时钟周期,最低积分时间是主时钟周期,最大积分时间为318个主时钟周期;
复位信号RESET,通过迫使焦平面第一行的信号积分,使读出电路的操作复位,在探测器主时钟MC的上升沿过程中,复位信号改变其状态,每帧不得重复一次以上;
数据有效信号DATAVALID在积分结束后自动设置到一个高电平;
第一行读出信号LIGNE1自动设置在高电平,以获得来自第一行的FPA的信号输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210080080.7A CN103076552B (zh) | 2012-03-23 | 2012-03-23 | 一种抑制偏置电压纹波和列输出放大器通道不均匀的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210080080.7A CN103076552B (zh) | 2012-03-23 | 2012-03-23 | 一种抑制偏置电压纹波和列输出放大器通道不均匀的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103076552A CN103076552A (zh) | 2013-05-01 |
CN103076552B true CN103076552B (zh) | 2015-05-13 |
Family
ID=48153145
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210080080.7A Expired - Fee Related CN103076552B (zh) | 2012-03-23 | 2012-03-23 | 一种抑制偏置电压纹波和列输出放大器通道不均匀的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103076552B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114205488B (zh) * | 2020-08-31 | 2023-10-27 | 京东方科技集团股份有限公司 | 图像传感电路、图像传感器以及图像传感电路的驱动方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5420419A (en) * | 1992-06-19 | 1995-05-30 | Honeywell Inc. | Camera for producing video output signal, infrared focal plane array package for such camera, and method and apparatus for generating video signals from passive focal plane array of elements on a semiconductor substrate |
CN1301080A (zh) * | 2000-12-28 | 2001-06-27 | 西安交通大学 | 一种用于抑制大功率直流电源纹波的方法及装置 |
US7105821B1 (en) * | 2003-12-15 | 2006-09-12 | Raytheon Company | Thermally stabilized radiation detector utilizing temperature controlled radiation filter |
CN101212172A (zh) * | 2006-12-25 | 2008-07-02 | 中兴通讯股份有限公司 | 开关型电源变换器限流状态输出电压低频纹波抑制电路 |
CN101951133A (zh) * | 2010-08-11 | 2011-01-19 | 京信通信系统(中国)有限公司 | 直放站供电系统上电瞬间的冲击电流抑制装置 |
-
2012
- 2012-03-23 CN CN201210080080.7A patent/CN103076552B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5420419A (en) * | 1992-06-19 | 1995-05-30 | Honeywell Inc. | Camera for producing video output signal, infrared focal plane array package for such camera, and method and apparatus for generating video signals from passive focal plane array of elements on a semiconductor substrate |
CN1301080A (zh) * | 2000-12-28 | 2001-06-27 | 西安交通大学 | 一种用于抑制大功率直流电源纹波的方法及装置 |
US7105821B1 (en) * | 2003-12-15 | 2006-09-12 | Raytheon Company | Thermally stabilized radiation detector utilizing temperature controlled radiation filter |
CN101212172A (zh) * | 2006-12-25 | 2008-07-02 | 中兴通讯股份有限公司 | 开关型电源变换器限流状态输出电压低频纹波抑制电路 |
CN101951133A (zh) * | 2010-08-11 | 2011-01-19 | 京信通信系统(中国)有限公司 | 直放站供电系统上电瞬间的冲击电流抑制装置 |
Non-Patent Citations (5)
Title |
---|
640×512制冷探测器积分时间与SNR分析;刘宁等;《红外与激光工程》;20120131;第41卷(第1期);第43-48页 * |
640×512红外焦平面探测器前端噪声分析及抑制技术;刘宁等;《红外技术》;20101031;第32卷(第10期);第572-575、582页 * |
Influencing factors of microscanning performance based on flat optical component;Xiubao Sui.etc;《CHINESE OPTICS LETTERS》;20110510;第052302-1-052302-4页 * |
基于BP神经网络的红外焦平面非均匀性校正技术;王娴雅等;《光子学报》;20090630;第38卷(第6期);第1504-1506页 * |
基于配准的红外焦平面阵列条纹非均匀性校正;任建乐等;《红外与毫米波学报》;20111231;第30卷(第6期);第499-502、506页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103076552A (zh) | 2013-05-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10616512B2 (en) | Systems, methods, and media for high dynamic range imaging using dead-time-limited single photon detectors | |
Kavusi et al. | Quantitative study of high-dynamic-range image sensor architectures | |
US10015429B2 (en) | Method and system for reducing noise in an image sensor using a parallel multi-ramps merged comparator analog-to-digital converter | |
EP2974280B1 (en) | Threshold-monitoring, conditional-reset image sensor | |
Fowler et al. | A 5.5 Mpixel 100 frames/sec wide dynamic range low noise CMOS image sensor for scientific applications | |
US8922685B2 (en) | Solid-state image sensing device and camera system with divided pixels | |
US9609243B2 (en) | Systems and methods for providing low-noise readout of an optical sensor | |
US20170195590A1 (en) | Method and system for reducing noise in an image sensor using a parallel multi-ramps merged comparator analog-to-digital converter | |
CN101272453A (zh) | 拖尾校正信号发生电路、方法、程序、拖尾校正电路和成像器件 | |
Bugalho et al. | Experimental results with TOFPET2 ASIC for time-of-flight applications | |
US20210377470A1 (en) | Digital pixel having high sensitivity and dynamic range | |
Daigle et al. | The darkest EMCCD ever | |
Daigle et al. | CCCP: a CCD controller for counting photons | |
CN103076552B (zh) | 一种抑制偏置电压纹波和列输出放大器通道不均匀的方法 | |
Liu et al. | A dual-exposure wide dynamic range CMOS image sensor with 12 bit column-parallel incremental sigma-delta ADC | |
TWI683577B (zh) | 具有電流鏡之緩衝直接注入讀出中的影像滯後緩解 | |
Ren et al. | Efficient single image stripe nonuniformity correction method for infrared focal plane arrays | |
Daigle et al. | Photon counting with an EMCCD | |
Weidenspointner et al. | Calibration of the non-linear system response of a prototype set-up of the DSSC detector for the European XFEL | |
Greffe et al. | Characterization of low light performance of a CMOS Sensor for ultraviolet astronomical applications | |
CN115278100B (zh) | 一种像素单元电路、信号采集装置和信号采集方法 | |
Guo et al. | Column readout circuit with improved offset mismatch and charge sharing for CMOS image sensor | |
Zhu et al. | Explosive growth of image sensors in smart government technology and economic Scale | |
Contreras et al. | A hardware Kalman-based offset estimator for nonuniformity correction on IRFPA | |
Kamberova et al. | The effect of radiometric correction on multicamera algorithms |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150513 Termination date: 20190323 |