CN108955571A - 双频外差与相移编码相结合的三维测量方法 - Google Patents

双频外差与相移编码相结合的三维测量方法 Download PDF

Info

Publication number
CN108955571A
CN108955571A CN201810470370.XA CN201810470370A CN108955571A CN 108955571 A CN108955571 A CN 108955571A CN 201810470370 A CN201810470370 A CN 201810470370A CN 108955571 A CN108955571 A CN 108955571A
Authority
CN
China
Prior art keywords
phase
striped
heterodyne
shift coding
shift
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.)
Granted
Application number
CN201810470370.XA
Other languages
English (en)
Other versions
CN108955571B (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.)
Nanchang Hangkong University
Original Assignee
Nanchang Hangkong University
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 Nanchang Hangkong University filed Critical Nanchang Hangkong University
Priority to CN201810470370.XA priority Critical patent/CN108955571B/zh
Publication of CN108955571A publication Critical patent/CN108955571A/zh
Application granted granted Critical
Publication of CN108955571B publication Critical patent/CN108955571B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/24Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/25Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures by projecting a pattern, e.g. one or more lines, moiré fringes on the object
    • G01B11/254Projection of a pattern, viewing through a pattern, e.g. moiré
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T9/00Image coding
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00Indexing scheme for image data processing or generation, in general
    • G06T2200/04Indexing scheme for image data processing or generation, in general involving 3D image data

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Multimedia (AREA)
  • Theoretical Computer Science (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了双频外差与相移编码相结合的三维测量方法,由外差原理、相移编码原理、及三维测量原理三大关键部分组成。本发明的优点是:(1)双频外差后的条纹的周期无需覆盖整个视场,也即,两条纹的节距相差可以相差较大,从而可以提高分辨被测细节的能力,同时更容易满足外差限制条件V0表示初始相位主值与外差后相位频率的比值;表示相位主值的误差;ΔΦ表示外差后的相位误差。即,对相位主值的精度有较高的要求)。(2)用N个周期的相位编码条纹就可以解出大于2N个周期的绝对相位,从而达到用频率较低的相位编码条纹,实现高精度的表面形貌测量。

Description

双频外差与相移编码相结合的三维测量方法
技术领域
本发明涉及一种三维测量方法,尤其涉及多频外差和相移编码三维测量方法。
背景技术
条纹投影轮廓术因其高速、高精度、非接触、全场测量、快速信息获取等优点,在三维测量中有重要意义,已广泛应用于工业制造、文物保护、医疗等诸多领域。本发明的三维测量系统如图1所示,包括DLP投影仪1、CCD2、工作站 3、测量支架4、参考平面5和待测物体6;DLP投影仪1和CCD2放在测量支架4上;DLP投影仪1、CCD2分别通过数据线连接工作站3;待测物体6放在参考平面5上;工作站3内包含图像采集卡、投影软件、测量软件。DLP投影仪1将带有特征信息的条纹聚焦投射到被测物体6表面,由CCD2采集条纹信息,经过工作站3处理后提取出特征信息,并按照特定算法进行三维重建。DLP投影仪1光轴和CCD2光轴相交于O点。DLP投影仪1和CCD2为同一高度,它们之间的距离为d,它们到参考平面的距离为l0。被测物体6的高度计算公式为:
其中fy为参考平面上的正弦条纹在竖直方向上的分频率,为物体表面图像和参考平面图像对应点的连续相位差。
随着科技的发展,人们对物体的三维数据要求不断提高。通过对国内外研究现状及发展动向分析研究,国内外研究者开始对条纹的类型、反射率、强度、对比度等进行研究,以期扩大测量对象的范围,以及解决测量结果不理想的问题。由于多频外差方法有较高的测量精度,但是容易出现跳变误差,且随着外差的次数增加越容易出现跳变误差,而且对选用的频率有一定的限制,为了减少外差次数和放宽对频率的限制,同时,相移编码条纹可以求出比自身周期大的条纹级次,因此想到本发明的思路。
发明内容
本发明的目的在于提出一种双频外差与相移编码相结合的三维测量方法,此方法能够提高双频外差的测量精度,而且能以低周期相位编码获得高周期的绝对相位。
双频外差与相移编码相结合的三维测量方法,(1)双频外差原理:外差原理是运用光调制技术以及电子相位计测量经过调制的光信号的相位,将频率为f1和f2的相位函数φ1(x,y)和φ2(x,y)经过叠加得到一种合适的频率为f 其相位函数φ12(x,y)。
多频外差原理基本思想是利用合成的条纹图使其周期能覆盖整个视场,然后以合成的相位图的相位值来解得原始条纹的相位。此方法利用合成的低频的测量结果来确定原始的高频测量时相位展开值,从而就能够得到真实相位,为了在全场范围内准确地进行相位展开,必须选择合适的f1和f2值,使得f=1。而本发明无需f等于1,使得f1和f2值的选取可以根据实验需要进行取值。
(2)相移编码条纹原理:将相移编码条纹的周期设置为双频外差后条纹周期的二分之一,即相移编码周期N=1/2f。利用相移编码条纹,求出双频外差后的相位的条纹级次k,由条纹级次k将展开成连续相位Φ12(x)。
具体包括以下步骤:
1)生成两种频率的条纹:根据实验需求选择条纹的频率f1和f2,设置具有相移功能的正弦条纹I1(x,y)和I2(x,y),其强度表示为:
I1(x,y)=A(x,y)+B(x,y)cos[φ1(x,y)+π(n-1)/2] (1)
I2(x,y)=A(x,y)+B(x,y)cos[φ2(x,y)+π(n-1)/2] (2)
其中,A(x,y)是背景光强,B(x,y)是调制光强,φ1(x,y)和φ2(x,y)是待求相位,n根据相移步数N而改变,取大于0小于N的整数,本发明以N=4进行四步相移为例,则(n=1,2,3,4)。
2)编写相移编码条纹:编码条纹的相位平均分配于2π的长度段中,像素点(x, y)的阶梯编码相位可表示为:
其中[x/p]=k(x,y)为通过取整运算得到的理想条纹级次,N为相移编码条纹总的周期数目N=1/2λ,p为条纹间距(即每个周期的像素个数)。将阶梯型编码相位加到相移条纹中得到相移编码条纹:
3)用投影仪依次投射上述生成的正弦条纹和相移编码条纹于参考面和物体上,用CCD采集这三组图片,送入计算机中。
4)由两组四步正弦条纹可以得到包裹相位将包裹相位经外差处理得到包裹相位
5)通过相移编码条纹可以得到的条纹级次k(x,y)。
6)利用条纹级次k(x,y)将包裹相位进行相位展开得到连续相位φ12(x,y),展开公式为:
7)用进行相位展开得到φ1(x,y)和φ2(x,y),展开公式为:
其中α、α1、α2的相位角。
8)根据上述方法,进行相位进展开,分别得到两个频率参考面的连续相位α1和α2及带物体的两个频率的连续相位β1和β2
9)通过公式:得到连续相位差,从而利用相位-高度公式:
最后得到两个物体表面每一点的高度信息,可以互相校正得到更精确的高度信息。
本发明的优点是:1.双频外差后的条纹的周期无需覆盖整个视场,也即,两条纹的节距相差可以相差较大,从而可以提高分辨被测细节的能力,同时更容易满足外差限制条件(V0表示初始相位主值与外差后相位频率的比值;表示相位主值的误差;ΔΦ表示外差后的相位误差。即,对相位主值的精度有较高的要求)。2.用N个周期的相位编码条纹就可以解出大于2N个周期的绝对相位,从而达到用频率较低的相位编码条纹,实现高精度的表面形貌测量。
附图说明
图1为本发明的三维测量系统示意图。
图2为本发明的双频条纹。
图3为本发明的相移编码条纹。
图4为本发明的双频外差原理图。
图5为本发明的文本流程图。
图6为包裹相位φ12(x,y)和条纹级次k(x,y)。
图7为相位展开过程。
图8为展开相位φ1(x,y)。
附图标记说明:1.DLP投影仪、2.CCD、3.工作站、4.测量支架、5.参考平面、 6.待测物体。
具体实施方式
下面结合附图对本发明的优选实施例进行详细阐述,以使本发明的优点和特征能更易于被本领域技术人员理解,从而对本发明的保护范围做出更为清楚明确的界定。
具体包括以下步骤:
1)生成两种频率的条纹:根据实验需求选择条纹的频率f1和f2,设置具有相移功能的正弦条纹I1(x,y)和I2(x,y),如图2所示,其强度表示为:
I1(x,y)=A(x,y)+B(x,y)cos[φ1(x,y)+π(n-1)/2] (1)
I2(x,y)=A(x,y)+B(x,y)cos[φ2(x,y)+π(n-1)/2] (2)
其中,A(x,y)是背景光强,B(x,y)是调制光强,φ1(x,y)和φ2(x,y)是待求相位,n根据相移步数N而改变,取大于0小于N的整数,本发明以N=4进行四步相移为例,则(n=1,2,3,4)。
2)编写相移编码条纹:编码条纹的相位平均分配于2π的长度段中,像素点(x, y)的阶梯编码相位可表示为:
其中[x/p]=k(x,y)为通过取整运算得到的理想条纹级次,N为相移编码条纹总的周期数目N=1/2λ,p为条纹间距(即每个周期的像素个数)。将阶梯型编码相位加到相移条纹中得到相移编码条纹:
3)用投影仪依次投射上述生成的正弦条纹和相移编码条纹于参考面和物体上,用CCD采集这三组图片,送入计算机中。
4)由两组四步正弦条纹可以得到包裹相位将包裹相位经外差处理得到包裹相位
5)通过相移编码条纹可以得到的条纹级次k(x,y)。
6)利用条纹级次k(x,y)将包裹相位进行相位展开得到连续相位φ12(x,y),展开公式为:
7)用φ12(x,y)对进行相位展开得到φ1(x,y)和φ2(x,y),展开公式为:
其中α、α1、α2的相位角。
8)根据上述方法,进行相位进展开,分别得到两个频率参考面的连续相位α1和α2及带物体的两个频率的连续相位β1和β2
9)通过公式:得到连续相位差,从而利用相位-高度公式:
最后得到两个物体表面每一点的高度信息,可以互相校正得到更精确的高度信息。
实施例1
1、投影计算机生成的两种频率的正弦条纹(如图2)和相位编码条纹(如图3)到参考平面和待测物体上,由相机捕获调制后的条纹图。以下以平面仿真为例进行相位展开说明。
2、用计算机处理相机捕获四步相移正弦条纹图和相移条纹图,得到包裹相位φ12(x,y)和条纹级次k(x,y)(如图6)。
3、利用和条纹级次k(x,y)通过公式(5)用计算机对进行相位展开得到展开相位φ12(x,y)。
4、利用公式(6)用计算机进行处理得到的展开相位φ1(x,y)(如图8)。
5、重复2-4对带物体的图像进行数据处理,得到带物体的展开相位,与参考面的展开相位作差,得到相位差通过公式(8)得到物体高度信息。
不局限于此,任何不经过创造性劳动想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书所限定的保护范围为准。

Claims (4)

1.双频外差与相移编码相结合的三维测量方法,其特征在于,包括以下步骤:
1)生成两种频率的条纹:根据实验需求选择条纹的频率f1和f2,设置具有相移功能的正弦条纹I1(x,y)和I2(x,y),其强度表示为:
I1(x,y)=A(x,y)+B(x,y)cos[φ1(x,y)+π(n-1)/2] (1)
I2(x,y)=A(x,y)+B(x,y)cos[φ2(x,y)+π(n-1)/2] (2)
2)编写相移编码条纹:编码条纹的相位平均分配于2π的长度段中,像素点(x,y)的阶梯编码相位可表示为:
将(4)中阶梯型编码相位加到相移条纹中得到相移编码条纹:
3)用投影仪依次投射上述生成的正弦条纹和相移编码条纹于参考面和物体上,用CCD采集调制的相位编码条纹图片和相位编码条纹的强度分布图,送入计算机中;
4)由两组四步正弦条纹可以得到包裹相位将包裹相位经外差处理得到包裹相位
5)通过相移编码条纹可以得到的条纹级次k(x,y)。
6)利用条纹级次k(x,y)将包裹相位进行相位展开得到连续相位φ12(x,y),展开公式为:
7)用φ12(x,y)对进行相位展开得到φ1(x,y)和φ2(x,y),展开公式为:
α、α1、α2的相位角;
8)根据上述方法,进行相位进展开,分别得到两个频率参考面的连续相位α1和α2及带物体的两个频率的连续相位β1和β2
9)通过公式:得到连续相位差,从而利用相位-高度公式:
最后得到两个物体表面每一点的高度信息,可以互相校正得到更精确的高度信息。
2.根据权利要求1所述的双频外差与相移编码相结合的三维测量方法,其特征在于,步骤1)中,A(x,y)是背景光强,B(x,y)是调制光强,φ1(x,y)和φ2(x,y)是待求相位,n根据相移步数N而改变,取大于0小于N的整数。
3.根据权利要求1所述的双频外差与相移编码相结合的三维测量方法,其特征在于,步骤2)中,[x/p]=k(x,y)为通过取整运算得到的理想条纹级次,N为相移编码条纹总的周期数目N=1/2f,p为条纹间距,即每个周期的像素个数。
4.根据权利要求1所述的双频外差与相移编码相结合的三维测量方法,其特征在于,步骤7)中
CN201810470370.XA 2018-05-16 2018-05-16 双频外差与相移编码相结合的三维测量方法 Active CN108955571B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810470370.XA CN108955571B (zh) 2018-05-16 2018-05-16 双频外差与相移编码相结合的三维测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810470370.XA CN108955571B (zh) 2018-05-16 2018-05-16 双频外差与相移编码相结合的三维测量方法

Publications (2)

Publication Number Publication Date
CN108955571A true CN108955571A (zh) 2018-12-07
CN108955571B CN108955571B (zh) 2019-07-09

Family

ID=64499399

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810470370.XA Active CN108955571B (zh) 2018-05-16 2018-05-16 双频外差与相移编码相结合的三维测量方法

Country Status (1)

Country Link
CN (1) CN108955571B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109883353A (zh) * 2019-02-19 2019-06-14 西安理工大学 一种刚性运动物体的动态三维测量方法
CN109974626A (zh) * 2019-04-08 2019-07-05 四川大学 一种基于相移量编码条纹级次的结构光三维测量方法
CN110081817A (zh) * 2019-04-24 2019-08-02 杭州光粒科技有限公司 消除背景光的方法、装置、计算机设备和存储介质
CN110132430A (zh) * 2019-03-29 2019-08-16 黑龙江科技大学 相移法两级编码高精度绝对相位获取方法
CN110207621A (zh) * 2019-06-04 2019-09-06 重庆大学 一种基于变相移的四步相移的面结构光解相方法
CN110230997A (zh) * 2019-06-04 2019-09-13 江南大学 一种基于改进单调法的阴影区相位噪声校正方法
CN112330771A (zh) * 2020-10-29 2021-02-05 南京理工大学智能计算成像研究院有限公司 一种用于快速三维测量的pcps条纹图案设计方法
CN113465545A (zh) * 2021-06-30 2021-10-01 北京航空航天大学 一种基于高速led阵列的三维测量系统及其测量方法
US20220364853A1 (en) * 2019-10-24 2022-11-17 Shining 3D Tech Co., Ltd. Three-Dimensional Scanner and Three-Dimensional Scanning Method

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101655360A (zh) * 2009-07-01 2010-02-24 四川大学 一种可绝对位相展开的十六步双频光栅相移轮廓术
CN102519390A (zh) * 2011-12-21 2012-06-27 哈尔滨理工大学 三编码周期灰度梯形相移结构光三维信息获取方法
JP5032943B2 (ja) * 2007-11-06 2012-09-26 パナソニック株式会社 3次元形状計測装置及び3次元形状計測方法
CN102853783A (zh) * 2012-09-18 2013-01-02 天津工业大学 一种高精度多波长三维测量方法
CN103759650A (zh) * 2014-02-07 2014-04-30 成都天拓众成科技有限公司 一种基于双频编码的快速三维测量方法
DE10362329B4 (de) * 2003-01-31 2015-07-30 Sirona Dental Systems Gmbh Aufnahmeverfahren für ein Bild eines Aufnahmeobjekts und Aufnahmevorrichtung
CN106840038A (zh) * 2017-01-25 2017-06-13 哈尔滨理工大学 一种双频相移与格雷码组合方法
US9690092B2 (en) * 2013-06-28 2017-06-27 Intel Corporation MEMS scanning mirror light pattern generation
CN107063128A (zh) * 2016-04-29 2017-08-18 华南师范大学 一种双频相移三维测量方法及系统

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10362329B4 (de) * 2003-01-31 2015-07-30 Sirona Dental Systems Gmbh Aufnahmeverfahren für ein Bild eines Aufnahmeobjekts und Aufnahmevorrichtung
JP5032943B2 (ja) * 2007-11-06 2012-09-26 パナソニック株式会社 3次元形状計測装置及び3次元形状計測方法
CN101655360A (zh) * 2009-07-01 2010-02-24 四川大学 一种可绝对位相展开的十六步双频光栅相移轮廓术
CN102519390A (zh) * 2011-12-21 2012-06-27 哈尔滨理工大学 三编码周期灰度梯形相移结构光三维信息获取方法
CN102853783A (zh) * 2012-09-18 2013-01-02 天津工业大学 一种高精度多波长三维测量方法
US9690092B2 (en) * 2013-06-28 2017-06-27 Intel Corporation MEMS scanning mirror light pattern generation
CN103759650A (zh) * 2014-02-07 2014-04-30 成都天拓众成科技有限公司 一种基于双频编码的快速三维测量方法
CN107063128A (zh) * 2016-04-29 2017-08-18 华南师范大学 一种双频相移三维测量方法及系统
CN106840038A (zh) * 2017-01-25 2017-06-13 哈尔滨理工大学 一种双频相移与格雷码组合方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
孙伟文等: "用于实时三维成像的双频结构光编解码方法", 《强激光与粒子束》 *
郭进等: "基于多频外差相移结构光的三维测量系统", 《工业仪表与自动化装置》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109883353A (zh) * 2019-02-19 2019-06-14 西安理工大学 一种刚性运动物体的动态三维测量方法
CN109883353B (zh) * 2019-02-19 2021-02-12 西安理工大学 一种刚性运动物体的动态三维测量方法
CN110132430A (zh) * 2019-03-29 2019-08-16 黑龙江科技大学 相移法两级编码高精度绝对相位获取方法
CN109974626B (zh) * 2019-04-08 2021-02-09 四川大学 一种基于相移量编码条纹级次的结构光三维测量方法
CN109974626A (zh) * 2019-04-08 2019-07-05 四川大学 一种基于相移量编码条纹级次的结构光三维测量方法
CN110081817A (zh) * 2019-04-24 2019-08-02 杭州光粒科技有限公司 消除背景光的方法、装置、计算机设备和存储介质
CN110207621A (zh) * 2019-06-04 2019-09-06 重庆大学 一种基于变相移的四步相移的面结构光解相方法
CN110230997A (zh) * 2019-06-04 2019-09-13 江南大学 一种基于改进单调法的阴影区相位噪声校正方法
CN110230997B (zh) * 2019-06-04 2020-04-21 江南大学 一种基于改进单调法的阴影区相位噪声校正方法
US20220364853A1 (en) * 2019-10-24 2022-11-17 Shining 3D Tech Co., Ltd. Three-Dimensional Scanner and Three-Dimensional Scanning Method
US12007224B2 (en) * 2019-10-24 2024-06-11 Shining 3D Tech Co., Ltd. Three-dimensional scanner and three-dimensional scanning method
CN112330771A (zh) * 2020-10-29 2021-02-05 南京理工大学智能计算成像研究院有限公司 一种用于快速三维测量的pcps条纹图案设计方法
CN112330771B (zh) * 2020-10-29 2023-12-15 南京理工大学智能计算成像研究院有限公司 一种用于快速三维测量的pcps条纹图案设计方法
CN113465545A (zh) * 2021-06-30 2021-10-01 北京航空航天大学 一种基于高速led阵列的三维测量系统及其测量方法

Also Published As

Publication number Publication date
CN108955571B (zh) 2019-07-09

Similar Documents

Publication Publication Date Title
CN108955571A (zh) 双频外差与相移编码相结合的三维测量方法
CN109489585B (zh) 基于改进多频条纹结构光的三维测量方法
Song et al. A new phase unwrapping algorithm based on three wavelength phase shift profilometry method
CN109186476A (zh) 一种彩色结构光三维测量方法、装置、设备及存储介质
CN103697815B (zh) 基于相位编码的混频结构光三维信息获取方法
CN107036556B (zh) 基于分段量化相位编码的结构光三维测量方法
CN109945802B (zh) 一种结构光三维测量方法
CN110006365B (zh) 基于二维查找表的相位展开方法、装置及电子设备
CN110174079A (zh) 一种基于四步相移编码型面结构光的三维重建方法
CN104197861A (zh) 基于结构光灰度向量的三维数字成像方法
CN110692084B (zh) 用于导出场景的拓扑信息的装置和机器可读存储介质
CN103791856A (zh) 基于四幅光栅条纹图像的相位求解与去包裹方法
CN109631798A (zh) 一种基于π相移方法的三维面形垂直测量方法
CN104154879B (zh) 一种非均匀条纹分段生成方法
CN112504165A (zh) 一种基于双边滤波优化的复合立体相位展开方法
CN105303572A (zh) 基于主被动结合的深度信息获取方法
Lv et al. An improved phase-coding method for absolute phase retrieval based on the path-following algorithm
CN109974624A (zh) 一种基于多频相移的减少投影图像数量的方法
CN104697468B (zh) 基于倾斜相位编码条纹的光学三维测量方法
Shi et al. A phase unwrapping method suitable for high frequency fringe based on edge feature
CN107421467B (zh) 结构光三维测量方法及系统
CN106840039B (zh) 一种基于复合光栅投影的快速三维测量方法
Song et al. Phase unwrapping method combining dual-frequency heterodyne with dual complementary phase encoding
Duan et al. Quasi-periodic phase coding for long-depth-range 3-D measurements of moving objects
CN111947600A (zh) 基于相位级次代价滤波的鲁棒立体相位展开方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant