CN102494638B - 三个编码周期的三基色三角形相移三维信息获取方法 - Google Patents
三个编码周期的三基色三角形相移三维信息获取方法 Download PDFInfo
- Publication number
- CN102494638B CN102494638B CN201110431404.2A CN201110431404A CN102494638B CN 102494638 B CN102494638 B CN 102494638B CN 201110431404 A CN201110431404 A CN 201110431404A CN 102494638 B CN102494638 B CN 102494638B
- Authority
- CN
- China
- Prior art keywords
- coding
- point
- code period
- sampled point
- primary
- 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
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提出了一种三个编码周期的三基色三角形相移三维信息获取方法，本发明方法通过选择三个基色空间中不同的编码周期，在保留相位编码方法对投影角连续划分的优点同时，对相位编码法的测量范围拓展了第三个编码周期倍。本发明是一种提高相移编码结构光测量范围，减小测量误差，提高测量精度，克服被测物表面反射率不一致和环境光对测量影响的基于三个编码周期三基色三角形相移采样点三维信息获取技术。
Description
技术领域
本发明涉及三维信息获取方法，具体涉及一种三基色编码周期互质且基色等级对称分布的相移结构光采样点三维信息获取技术。
背景技术
数量众多的光学非接触三维信息获取技术被广泛的应用于现代工农业生产过程中，其中结构光技术由于其主动测量原理，采样密度可调，测量精度高等优点被广泛应用于逆向工程、考古、医学、服饰、人体测量等需要恢复物体的三维形状的领域中。结构光技术是利用普通或者激光作为光源，将已知的图案投射到被测物的表面，根据通过成像装置从另一角度获取的图像，利用三角原理获得采样点的三维信息。
结构光技术根据投影图案的不同可以分为单点结构光、单线结构光、多线结构光和编码结构光等。其中单点和单线结构光由于其需要高精度的扫面机械及高精度的定位装置，已经不能引起广大研究者的兴趣。多线结构光由于使用时必须确定条纹对应的问题，也不适用于实际测量中。但编码结构光由于其条纹次序明确，不需考虑其他几何约束的优点已经成为目前该领域的研究热点。编码结构光技术向被测物投射编码图案，测量速度由于不需扫描装置而被大大提高了。在通过预先标定明确系统结构参数的条件下，确定被测物采样点并将其与条纹图像、编码图案对应起来，是编码结构光目前主要的研究问题。
编码结构光根据编码方法可分为时间编码、空间编码和直接编码、三者各具优点。相对于空间编码和直接编码，时间编码具有容易识别、准确度高等优点而具有广阔的应用前景。时间编码又可分为数字编码和连续编码两种技术，连续编码方法具有测量空间编码连续的优点，抗干扰能力强的优点，已经被越来越多研究者所重视。
目前连续编码方法主要分为正弦相移编码法和强度比编码法。正弦相移编码法具有采样密度高，在单独一个周期内，相位值连续且唯一，但是在多周编码中存在着周期数判定不准确，容易产生多义性，且该技术涉及大量的正切及反正切计算，计算量大，这些缺点限制了该技术的应用范围。强度比法恰恰弥补了正弦相移法计算量大的缺点。该技术利用线性变化光照下每一像素值与持续光照下像素值的比率，即强度比，替代正弦相位法中的相位值，具有速度快，计算量小的优点，但该技术不可避免的存在着解包裹时，采样点不唯一的问题。
大衍求一术被引入到连续正弦编码结构光解码中，目前编码方法采用两个编码频率对被测物进行采样，解码时除了求解同余方程外，还需要查询事先标定好的误差表进行修正，这类技术的采样精度没有明显的提高。
三周期正弦编码法也被用来进行正弦相移的解包裹算法，但是该方法不可避免正弦相移方法中计算量大，且投射编码图案众多的问题。
发明内容
本发明的目的是提供一种提高相移编码结构光测量范围，减小测量误差，提高测量精度，克服被测物表面反射率不一致和环境光对测量影响的三个编码周期的三基色三角形相移三维信息获取方法。
本发明的目的通过以下的步骤实现：
(1) 选择编码周期，制作编码图案。根据投射源及测量范围在三个基色空间上分别确定彼此互质的编码周期。编码图案是由投射源投射到被测物体表面的主动信息，是由在编码图案上，按X和Y方向顺序排列的编码点的基色等级构成。根据投射源与条纹图像获取装的几何位置确定编码方向，在另一个方向上是编码方向上编码点基色等级的重复。编码周期是指编码点的基色等级从某一个基色值历过基色值最高、基色值最低，且将此基色等级重复三次时，沿编码方向上编码点的总数。在一个周期内，基色等级的编码是线性的，并且关于1/2周期点对称。为了提高编码图像的解码时抗干扰能力和测量精度，对于一个基色空间采用相同的编码周期，制作三幅沿着编码方向依次仅移动1/3周期的编码图案，总计3幅基于三基色编码的彩色相移编码图案。
依次将3幅彩色编码相移图案经由投射源投影到被测物体表面，平面的编码图案受到被测物体表面三维信息的调制，产生变形，利用条纹图像获取装置，获取包含有被测物体及变形编码图案的3幅彩色条纹图像。
对同一个采样点的三幅同一基色内的基色数据进行处理，计算得到采样点在此基色空间内对应于编码图案的相对编码点，得到采样点在该及色空间内的余数方程。针对同一采样点在其余2个不同基色空间的其余幅条纹图像进行相同处理，总计得到3个余数方程，联立得到采样点对应绝对编码点的同余方程，采用大衍求一术确定同余方程的乘率，进而等到该同余方程的解，即采样点对应的绝对编码点，进而获得采样点对应的投影角度。
在确定采样点的对应的投影角度后，还需根据投射源与条纹图像获取装置的几何关系和三角原理确定采样点的三维信息。
本发明具有以下有益效果：
1.本发明方法通过选择三个基色空间中不同的编码周期，在保留相位编码方法对投影角连续划分的优点同时，对相位编码法的测量范围拓展了第三个编码周期倍。
本发明方法在相同测量范围内，可以选择三种周期较小的基色等级线性编码，这样减少了在同一幅值范围内基色等级分布的点数，在同等情况下，与大周期的彩色编码相比，基色等级分布的减少更加明显，出现误差的可能性大为降低。
本发明方法需要投射3幅编码图案，有效地减少了投射图案的数量，提高了三维信息获取的速度。
本发明方法对同一个基色空间，投射三幅移动1/3周期的编码图案，可以提高抗干扰能力，同时减少了环境光照和被测物表面反射系数不一致对三维信息获取结果的影响。
本发明方法采用编码点作为确定投影角的中间变量，可以减少计算量，进一步提高了三维信息获取的速度。
本发明方法在解码时通过三基色空间基色等级进行编码，扩大了绝对编码点的表示范围，将投影角进一步细分了。
本发明方法在解码时，不依靠条纹图像中采样点对应像点周围的像素点的数据，仅仅依靠本身的灰度信息，从原理上截断了误差传递的可能，避免了累计大误差的出现。
本发明提出的三基色编码周期互质且基色等级对称分布的编码解码技术即可解决投影角连续的划分问题，又可以提高测量准确度和采样密度。
附图说明
图1为三个编码周期的三基色三角形相移三维信息获取技术的实施流程框图；
图2为三角形基色等级编码排列序号；
图3为三角形基色等级编码的三个编码周期基色等级示意图；
图4为三基色三角形基色等级编码投射图案；
图5为解码后对应三个周期的相对编码点；
图6为求解同于方程后确定的绝对编码点。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步的说明：
本发明方法在实施中需要投射源，条纹图像获取装置和用于控制及计算的计算机。用来测量被测物表面被采样点的三维信息。其中被采样点的三维信息包括，在以投射源投射光学系统光学中心为原点，以原点与条纹图像获取装置照相光学系统光学中心连线方向为X轴方向，以投射源投射光学系统的光轴方向为Z轴方向的左手空间坐标系(世界坐标系)中，被采样点距离三个坐标轴的距离x，y，z，或者根据(x，y，z)及必要的空间解析几何计算方法，可以获取的诸如被采样点处得法向量、曲率或结合其他被采样点可以计算的被测物表面的诸如弯曲、倾斜、扭曲、翘曲之类的信息。
图1是三个编码周期的三基色三角形相移三维信息获取技术的实施流程框图。在使用前需要对测量系统包括的投射器、条纹图像获取装置和两者之间的距离进行标定，本领域的技术人员不难理解这一步骤是利用三角法进行三维信息获取的必要步骤。
按3个编码周期生成基色等级的三角形编码，每个编码周期各生成沿编码方向移动1/3周期的编码图案，共计3幅彩色编码图案，利用投射器，向被测物投射3幅三基色彩色编码编码图案，虽然此处编码图案的生产和编码图案的投射是顺序描述的，但本领域的技术人员不难理解可以通过同步技术实现在前1幅编码图案投射的时间内生成下1幅编码图案，以供下一个投射时间使用。
通过条纹图像获取装置，采集包含由被测物表面所产生变形的条纹及环境背景的条纹图像，由分析、控制和处理装置对条纹图像进行预处理后，确定条纹图像中采样点对应于编码图案中的相对编码点，建立该基色空间下相同编码周期的余数方程，3个基色空间内的3个余数方程联立得到同余方程组，根据大衍求一术确定乘率，得到采样点在编码图案中的绝对编码点，确定采样点对应的投影角度，根据三角原理获得采样点的三维信息。
图2是由计算机生成的编码图案族中某一个基色空间基色等级编码排列的序号。计算机控制投射源，按照设计的投射顺序和投影时间，依次将图2所示的三基色三角形基色等级编码的编码图案投射到被测物表面。编码公式如式所示
式中I
max
和I
min
是灰度编码时赋予编码点的最大基色等级级和最小基色等级级，mod(n,T
i
)表示编码点的在编码图案中同一行内的序号对周期T
i
取余得到余数。
图3是本发明方法采用的三基色三角形基色等级编码的基色等级示意图。其中每一幅编码图案根据投射源最大投射分辨力，按照行的编码方向，以周期T
i
(i=1，2，3)重复编码，赋予行方向编码点以具体的基色等级，列方向的编码点的基色等级分别重复同一列中上一行编码点的基色等级。本领域的技术人员将会理解，编码图案是列方向平行，行方向基色等级变化的彩色条纹编码图案。
图4是本发明方法实际应用时通过投射源向被测物表面投射的编码图案。原理上在一个编码周期T
i
(即同一个基色空间)内采用一幅三角形基色等级编码图案和一幅最大基色等级编码图案/或一幅最小基色等级编码图案，也可用于计算被测物表面采样点的三维信息，但是为了克服环境光及被测物表面采样点表面反射系数不同对测量的影响，采用投射3幅灰度级变化的编码图案，可以有效的提高采样点灰度计算时的鲁棒性及抗干扰的能力。
计算机生成的编码图案是在三个基色空间采用互质的编码周期，且不同编码图案中相同基色空间的编码依次沿编码方向上移动1/3周期组成，总计3幅彩色编码图案；每一个基色空间内的编码图案是三角形编码；编码周期是指以任意一点为起点，沿着编码方向基色等级级由0到255变化一次且该点基色等级出现三次时编码点的前一个编码点与起点之间的编码点的总数。三个编码周期对应的三角形基色等级编码周期是任意两个周期都为互质的，此处互质是指两个编码周期的最大公约数为1。
条纹图像获取装置所获取的条纹图像包括：被测物在多幅编码图案照射时，物体表面在多个照明条件下的反射光；通过转换电路，转换为能被识别并可处理的3幅条纹图像；其获取条纹图像的时间及方式由计算机根据编码图案的周期及编码图案的投射时间计算确定。
对同一基色空间的3幅条纹图像进行基色等级计算，确定采样点在该基色空间内编码图案中对应相对编码点r(m，n)，求解公式如式(2)所示。
式中I
1
、I
2
和I
3
是同一编码周期获取的三幅条纹图像中，同一采样点的不同基色等级，T为对应的编码周期。
根据编码周期T
i
、相对编码点r
i
(m，n)，构造余数a
i
，建立相对编码点r
i
(m，n)，在编码周期T
i
内的余数方程r
i
(m，n)≡a
i
(mod T
i
)；其中i表示三个编码周期的编号，i=1,2,3；
图5是以编码周期T
1
=3，T
2
=4，T
3
=5为例，条纹图像采样点，在不同编码周期的编码图案照射条件下，求解对应编码周期内余数方程可能的相对编码点值。
如果仅一个编码周期进行编解码的话，那么超出一个周期以外的采样点对应的编码点不唯一。
可以看出利用两个编码周期可以确定两个周期乘积个绝对编码点，但是此时为了获得到更高的采样精度，必须采用以下几种方法：
1. 增大每一个编码周期，以满足在整个测量空间内采样点与编码点唯一对应，但是编码周期的增大，不可避免的使采样点灰度信息获取容易受到其他采样点反射灰度信息的影响，同时容易产生误差；
2. 在不改变同一周期内灰度等级数量的前提下，只能增加编码图案中，同一灰度编码点的个数，这样在原理上降低了采样精度；
3. 额外增加其他约束条件，为了获取到与本方法相同的采样精度，必须与投射同一灰度编码点个数相同的线移编码图案，大大超出本方法投射的编码图案数量，降低了整体测量速度。
两个编码周期利用了两个基色空间，组成的颜色不完整，且获取的数据存在1/3的冗余。
建立关于采样点在3个基色空间编码图案照射条件下的投射源对应的绝对编码点x和3个编码周期的同余方程组：x(m，n)≡a
i
(mod T
i
)；
计算得到采样点对应于编码图案上的绝对编码点x(m，n)：根据大衍求一术确定同余方程组的乘率K
i
，利用公式(3)进而等到该同余方程组的解；
根据同余方程组的解，可以确定采样点在编码图案中对应的绝对编码点。
图6是是以编码周期T
1
=3，T
2
=4，T
3
=5为例，条纹图像采样点，在不同编码周期的编码图案照射条件下，求解对应编码周期内余数方程可能的相对编码点值。
确定采样点在编码图案中对应的绝对编码点x(m，n)后，根据预先标定的投射源投影光学系统的投射范围角度、条纹图像获取装置的照相光学系统的投影光学系统与照相光学系统之间的实际距离，计算得到采样点的三维信息。
本发明方法应用系统可以由一个投射编码图案的DLP投射器、一个或多个采集编码图像的数字相机、一个产生编码图案、进行系统控制和数据处理的笔记本电脑所组成，其特点是不含运动件、图案编码方便，而且所有部件均可在市场上购得。系统标定通过控制带有特征点的标准平面在视场内移动，采用照相测量方法完成。现有工艺水平完全可以满足原理样机和标定装置的加工要求。
Claims (1)
1.一种三个编码周期的三基色三角形相移三维信息获取方法，其特征在于它通过以下的步骤实现：
(1) 选择编码周期，制作编码图案：根据投射源及测量范围在三个基色空间上分别确定彼此互质的编码周期；编码图案是由投射源投射到被测物体表面的主动信息，是由在编码图案上，按X和Y方向顺序排列的编码点的基色等级构成；根据投射源与条纹图像获取装的几何位置确定编码方向，在另一个方向上是编码方向上编码点基色等级的重复；编码周期是指编码点的基色等级从某一个基色值历过基色值最高、基色值最低，且将此基色等级重复三次时，沿编码方向上编码点的总数；在一个周期内，基色等级的编码是线性的，并且关于1/2周期点对称；一个基色空间采用相同的编码周期；选择三个基色空间中不同的编码周期，在保留相位编码方法对投影角连续划分的优点同时，对相位编码法的测量范围拓展了第三个编码周期倍；三个编码周期对应的梯形基色等级编码周期是任意两个周期都为互质的，此处互质是指两个编码周期的最大公约数为1；选择三种周期中周期小的基色等级线性编码，这样减少了在同一幅值范围内基色等级分布的点数，在同等情况下，与大周期的彩色编码相比，基色等级分布的减少更加明显，出现误差的可能性大为降低；为了提高编码图像的解码时抗干扰能力和测量精度，对于一个基色空间采用相同的编码周期，制作三幅沿着编码方向依次仅移动1/3周期的编码图案，总计3幅基于三基色编码的彩色相移编码图案；
(2) 依次将3幅彩色编码相移图案经由投射源投影到被测物体表面，平面的编码图案受到被测物体表面三维信息的调制，产生变形，利用条纹图像获取装置，获取包含有被测物体及变形编码图案的3幅彩色条纹图像；
(3) 对同一个采样点的三幅同一基色内的基色数据进行处理，计算得到采样点在此基色空间内对应于编码图案的相对编码点，得到采样点在该基色空间内的余数方程；确定采样点在该基色空间内编码图案中对应相对编码点r(m，n)，求解公式如式(1)所示； 
公式（1）；式中I1、I2和I3是同一编码周期获取的三幅条纹图像中，同一采样点的不同基色等级，T为对应的编码周期；根据编码周期Ti、相对编码点ri(m，n)，构造余数ai，建立相对编码点ri(m，n)，在编码周期Ti内的余数方程ri(m，n)≡ai(mod Ti)；条纹图像采样点，在不同编码周期的编码图案照射条件下，求解对应编码周期内余数方程的相对编码点值；针对同一采样点在其余2个不同基色空间的其余幅条纹图像进行相同处理，总计得到3个余数方程，联立得到采样点对应绝对编码点的同余方程，采用大衍求一术确定同余方程的乘率，进而得到该同余方程的解，即采样点对应的绝对编码点，建立关于采样点在3个基色空间编码图案照射条件下的投射源对应的绝对编码点x和3个编码周期的同余方程组：x(m，n)≡ai(mod Ti)；根据大衍求一术确定同余方程组的乘率Ki，利用公式(2)进而得到该同余方程组的解，计算得到采样点对应于编码图案上的绝对编码点x(m，n)；确定采样点在编码图案中对应的绝对编码点x(m，n)后，进而获得采样点对应的投影角度；公式（2）；
(4) 在确定采样点的对应的投影角度后，还需根据投射源与条纹图像获取装置的几何关系，预先标定的投射源投影光学系统的投射范围角度、条纹图像获取装置的照相光学系统的投影光学系统与照相光学系统之间的实际距离和三角原理；确定采样点的三维信息。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110431404.2A CN102494638B (zh) | 2011-12-21 | 2011-12-21 | 三个编码周期的三基色三角形相移三维信息获取方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110431404.2A CN102494638B (zh) | 2011-12-21 | 2011-12-21 | 三个编码周期的三基色三角形相移三维信息获取方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102494638A CN102494638A (zh) | 2012-06-13 |
| CN102494638B true CN102494638B (zh) | 2014-06-11 |
Family
ID=46186482
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201110431404.2A Expired - Fee Related CN102494638B (zh) | 2011-12-21 | 2011-12-21 | 三个编码周期的三基色三角形相移三维信息获取方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102494638B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103383249B (zh) * | 2013-07-12 | 2016-08-10 | 西安交通大学 | 灰度条纹投影光强非线性校正方法及基于该方法的相位校正方法 |
| CN108955574B (zh) * | 2018-07-09 | 2020-04-28 | 广东工业大学 | 一种三维测量方法及系统 |
| CN114659466B (zh) * | 2022-03-30 | 2023-04-28 | 西安交通大学 | 一种基于事件相机的边移结构光三维测量方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101738172A (zh) * | 2008-11-05 | 2010-06-16 | 哈尔滨理工大学 | 基于绿条纹分割的高采样密度彩色结构光三维测量方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4032603B2 (ja) * | 2000-03-31 | 2008-01-16 | コニカミノルタセンシング株式会社 | ３次元計測装置 |
-
2011
- 2011-12-21 CN CN201110431404.2A patent/CN102494638B/zh not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101738172A (zh) * | 2008-11-05 | 2010-06-16 | 哈尔滨理工大学 | 基于绿条纹分割的高采样密度彩色结构光三维测量方法 |
Non-Patent Citations (7)
| Title |
|---|
| JP特开2001-280916A 2001.10.10 |
| 三原色格雷码结合梯形相移三维测量方法;饶勇 等;《哈尔滨理工大学学报》;20110630;第16卷(第3期);77-81 * |
| 单鹂娜 等.彩色梯度相移与格雷码结合的三维测量方法.《光学精密工程》.2010,第18卷(第11期),2497-2504. |
| 单鹂娜.组合式时间编码结构光三维测量方法与技术.《全国博士学位论文数据库 信息科技辑》.2011,(第6期),第6-10页. |
| 彩色梯度相移与格雷码结合的三维测量方法;单鹂娜 等;《光学精密工程》;20101130;第18卷(第11期);2497-2504 * |
| 组合式时间编码结构光三维测量方法与技术;单鹂娜;《全国博士学位论文数据库 信息科技辑》;20110615(第6期);第6-10页 * |
| 饶勇 等.三原色格雷码结合梯形相移三维测量方法.《哈尔滨理工大学学报》.2011,第16卷(第3期),77-81. |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102494638A (zh) | 2012-06-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102494637B (zh) | 三个编码周期的三基色梯形相移三维信息获取方法 | |
| CN102519390B (zh) | 三编码周期灰度梯形相移结构光三维信息获取方法 | |
| Wang | Review of real-time three-dimensional shape measurement techniques | |
| CN101881605B (zh) | 基于相位编码技术的光学三维测量方法 | |
| CN107036556B (zh) | 基于分段量化相位编码的结构光三维测量方法 | |
| CN101441066B (zh) | 彩色条纹编码的相位去包裹方法 | |
| CN110174079B (zh) | 一种基于四步相移编码型面结构光的三维重建方法 | |
| CN105844633B (zh) | 基于De序列和相位编码的单帧结构光深度获取方法 | |
| CN102261896A (zh) | 一种基于相位测量的物体三维形貌测量方法及系统 | |
| CN104197861A (zh) | 基于结构光灰度向量的三维数字成像方法 | |
| CN110849290B (zh) | 基于形态学操作的分段量化编码强度的三维测量方法 | |
| CN106546193B (zh) | 一种高反射物体表面三维测量方法和系统 | |
| CN109945802A (zh) | 一种结构光三维测量方法 | |
| CN105303572B (zh) | 基于主被动结合的深度信息获取方法 | |
| Yu et al. | Unequal-period combination approach of gray code and phase-shifting for 3-D visual measurement | |
| CN102519396B (zh) | 三个灰度对称线性编码周期的采样点三维信息获取方法 | |
| CN111207692A (zh) | 一种改进的分段阶梯相位编码三维测量方法 | |
| Liu et al. | Background and amplitude encoded fringe patterns for 3D surface-shape measurement | |
| CN108596008A (zh) | 针对三维人脸测量的面部抖动补偿方法 | |
| CN102494638B (zh) | 三个编码周期的三基色三角形相移三维信息获取方法 | |
| CN112945089A (zh) | 一种基于条纹宽度调制的结构光编码方法 | |
| CN113763540A (zh) | 一种基于散斑条纹混合调制的三维重建方法及设备 | |
| CN107421467B (zh) | 结构光三维测量方法及系统 | |
| CN113345039B (zh) | 三维重建量化结构光相位图像编码方法 | |
| US10801834B2 (en) | Fringe projection for determining topography of a body |
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: 20140611 Termination date: 20161221 |