CN106017305B - 一种基于差分进化算法的相位解包裹方法 - Google Patents

一种基于差分进化算法的相位解包裹方法 Download PDF

Info

Publication number
CN106017305B
CN106017305B CN201610298044.6A CN201610298044A CN106017305B CN 106017305 B CN106017305 B CN 106017305B CN 201610298044 A CN201610298044 A CN 201610298044A CN 106017305 B CN106017305 B CN 106017305B
Authority
CN
China
Prior art keywords
vector
differential evolution
phase
follows
evolution 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.)
Active
Application number
CN201610298044.6A
Other languages
English (en)
Other versions
CN106017305A (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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong 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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201610298044.6A priority Critical patent/CN106017305B/zh
Publication of CN106017305A publication Critical patent/CN106017305A/zh
Application granted granted Critical
Publication of CN106017305B publication Critical patent/CN106017305B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00—Measuring instruments characterised by the use of optical techniques
    • G01B9/02—Interferometers
    • G01B9/02083—Interferometers characterised by particular signal processing and presentation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B9/00—Measuring instruments characterised by the use of optical techniques
    • G01B9/02—Interferometers
    • G01B9/02055—Reduction or prevention of errors; Testing; Calibration

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Complex Calculations (AREA)
  • Instruments For Measurement Of Length By Optical Means (AREA)

Abstract

本发明公开了一种基于差分进化算法的相位解包裹方法，其核心是利用有限项正交多项式(Zernike)多项式全局逼近真实相位，并利用差分进化算法求解待定系数。为了排除被大噪声污染的数据点影响最终的解包结果，采用加权Zernike多项式拟合的方法。本发明将Zernike多项式作为基底函数对测量得到的离散数据进行拟合，保证了拟合结果的精确性，利用差分进化算法(Differential Evolution，DE)动态随机搜索，以求解待定系数，具有较强的全局收敛性和鲁棒性。

Description

一种基于差分进化算法的相位解包裹方法
【技术领域】
本发明属于光干涉计量技术领域，涉及到差分进化算法(DifferentialEvolution，DE)，具体是一种基于差分进化算法的相位解包裹方法。
【背景技术】
干涉图的处理是光干涉计量中的关键技术。而相位解包裹是快速傅里叶变换(Fast Fourier Transformation，FFT)、移相、空间投影等技术引入反正切函数计算的干涉处理技术的关键。绝大多数的相位都是通过反正切函数求出，而反正切函数的值域为(-π,π]，这样求解出的相位被包裹在(-π,π]之间，相位解包就是把包裹在(-π,π]之间的相位还原成真实相位的过程。
为了很好地实现解包裹运算，国内外学者已经提出了很多相位解包裹算法，归纳起来主要是两种，一种是寻求积分结果与路径无关条件的路径跟踪算法，另一种是寻求满足最小范数解的路径无关相位解包裹算法。除此之外，还有全局算法以及分区域算法，各种方法都有其优缺点和适用性。
但是上述算法在大噪声下的鲁棒性都较差，因此，实有必要提供一种可以解决鲁棒性的技术方案。
【发明内容】
本发明的目的是改进现有技术存在的问题，采用加权Zernike多项式拟合的方法，这样被大噪声污染的数据点就不会影响最终的解包结果，提高了该方法对噪声的鲁棒性。具体地说，本发明提供一种基于差分进化算法的相位解包裹方法，该方法相比于基于微分Zernike多项式拟合的相位解包裹方法，不需要计算相位梯度和多项式的一阶导数。本发明首次将差分进化算法应用于光干涉计量领域的相位解包裹中。
为实现上述目的，本发明采用以下技术方案：
一种基于差分进化算法的相位解包裹方法，包括以下步骤：
(1)采用加权Zernike多项式拟合的方法建立实际干涉图的多项式，以剔除噪声污染的数据点，建立的多项式为：
其中，(x,y)为第i个像素的坐标，M为像素的数量，ci为多项式的系数，C＝【c1、c2、c3……cM】，W(x,y)为坐标(x,y)下的权重，Zi(x,y)为Zernike多项式；
(2)建立拟合系数的目标函数
对步骤(1)的多项式取余弦，建立多项式拟合系数C的最小化约束函数：
其中，L为干涉图有效孔径，λ为约束因子，γ为估计干涉图与实际干涉图的相关系数；
(3)利用DE算法对多项式拟合系数C的最小化约束函数进行全局优化，求解拟合系数，从而得到真实的解包相位。
进一步，假设相邻两个像素之间的相位之差的绝对值不超过π，则步骤(1)中的权重W(x,y)根据以下公式计算：
进一步，步骤(2)中的估计干涉图与实际干涉图的相关系数γ定义如下：
进一步，所述步骤(3)的具体步骤包括以下：
假设待优化的变量维数为D，进化代数为NP，那么进化的第G代其解向量Ci,G为：
初始种群在初始化时均匀随机覆盖整个搜索空间，则搜索空间的上下边界分别表示为：
初始种群第i个个体的第j个参数按照如下公式产生：
其中，rand(0,1)表示均匀分布在[0,1]之间的随机数；
得到初始种群之后，对其进行变异、交叉和选择，直到搜索到约束函数的全局极小值C。
进一步，采用DE/rand/1变异策略将当前种群的目标向量Ci,G进行变异得到突变向量Vi,G，所述的DE/rand/1变异策略表示为：
其中，分别为第一次变异时第一个个体、第二个个体、第三个个体的目标向量，正整数索引随机产生，且满足Cbest,G是第G代最优的个体，常数缩放因子F∈[0,2]控制差分向量对变异个体的影响大小；
突变向量Vi,G表示为:
进一步，对初始种群进行交叉的方法为：将每一对目标向量Ci,G和其相应的变异向量Vi,G进行交叉操作产生试验向量Ui,G，该试验向量Ui,G表示为:
Ui,G定义如下：
式中，表示第j对目标向量Ci,G和其相应的变异向量Vi,G进行交叉操作产生试验向量，交叉常数CR∈[0,1]为自定义常数，其大小控制从变异向量Vi,G拷贝比例的大小；如果新产生的试验向量Ui,G中的某些值超过设定值的上下边界，那么将其均匀随机地在其限定的边界内进行初始化，以符合要求。
进一步，计算所有试验向量Ui,G对应的目标函数f的值，并进行选择操作的具体方法为：将试验向量Ui,G所对应的目标函数值f(Ui,G)与当前种群中相应的目标向量Ci,G所对应的目标函数值f(Ci,G)进行比较，如果试验向量函数值f(Ui,G)小于或等于目标向量函数值f(Ci,G)，将试验向量体替代目标向量进入下一代种群，否则，目标向量仍然留在种群中进入下一代。
与现有技术相比，本发明至少具有以下有益效果：本发明在相位拟合时，采用加权Zernike多项式拟合的方法，这样被大噪声污染的数据点就不会影响最终的解包结果，消除了大噪声情况下拟合失真的问题，提高了该方法对噪声的鲁棒性。相比于基于微分Zernike多项式拟合的相位解包裹方法，本发明不需要计算相位梯度和多项式的一阶导数。本发明首次将差分进化算法应用于光干涉计量领域的相位解包裹中。
【附图说明】
图1是本发明的实施路线图；
图2是干涉图在不同噪声水平下的解包结果；
图3是点衍射干涉测量实验中采用5步相移法获取的5幅干涉条纹图；
图4是由干涉条纹图获取的包裹相位；
图5是对获取的包裹相位采用本发明的解包裹方法得到最终面型；
图6是ZYGO干涉仪的测量结果；
图7是本发明的解包结果与ZYGO干涉仪的解包结果之间的点对点残差。
【具体实施方式】
光学表面检测大多数情况中，被测光学表面或光学系统的出射波面总是趋于光滑且连续的，这样的波面函数一定可以表示成一个完备的基底函数的线性组合。因此常用Zernike多项式作为基底函数对测量得到的离散数据进行拟合，把实际波面或面形表示为Zernike多项式各项的线性组合，然后全局逼近真实相位，并利用优化算法求解待定系数。
而差分进化算法(Differential Evolution，DE)最早由Storn和Price于1997年提出，是一种简单有效基于群体智能的随机并行优化算法，通过模仿生物群体内个体间的合作与竞争产生的启发式群体智能来指导优化搜索，对于定义在连续空间上的全局优化问题非常有效。DE算法特有的记忆能力使其可以动态跟踪当前的搜索情况，以调整其搜索策略，实现自适应寻优，因此具有较强的全局收敛性和鲁棒性。
本发明采用加权Zernike多项式拟合的方法，这样被大噪声污染的数据点就不会影响最终的解包结果，提高了该方法对噪声的鲁棒性。
本发明所采用的技术方案包括以下步骤：
步骤1：剔除噪声污染的数据点
参见图2并结合图1，表示的是不同噪声水平下的解包结果，从上到下噪声量增多，由图看得出，在大噪声情况下拟合失真的问题，被大噪声污染的数据点会影响到最终的解包结果，对此，本发明采用加权Zernike多项式拟合的方法来消除噪声的影响；也就是将原有的多项式(式(1))改写为式(2)：
改写为:
其中，(x,y)为第i个像素的坐标，M为像素的数量，ci为多项式的系数，C＝【c1、c2、c3……cM】，W(x,y)为坐标(x,y)下的权重，Zi(x,y)为Zernike多项式。代表i个像素的真实相位。
首先假设相邻两个像素之间的相位之差的绝对值不超过π，则W(x,y)如下式所示：
代表坐标点(x,y)处的相位，代表坐标点(x+1,y)处的相位，代表坐标点(x,y+1)处的相位，abs表示取绝对值。
这样被大噪声污染的数据点就不会影响最终的解包结果，从而提高了该方法对噪声的鲁棒性。
步骤2：建立拟合系数的目标函数
对相位拟合等式两边同时取余弦，估计多项式拟合系数C，得到最小化约束函数f(C)；
对相位拟合等式两边同时取余弦，可得：
为了估计多项式拟合系数C，可以得到最小化约束函数：
其中L为干涉图有效孔径，λ为约束因子，γ为估计干涉图与实际干涉图的相关系数，定义如下：
步骤3：求解拟合系数
利用DE算法完成多模函数f(C)的全局优化；具体为：
假设待优化的变量维数为D，进化代数为NP，那么进化的第G代其解向量Ci,G可以表示为：
初始种群在初始化时应该尽可能地均匀随机覆盖整个搜索空间，则搜索空间的上下边界分别表示为：
初始种群(G＝0)第i个个体的第j个参数可以按照如下公式产生：
其中rand(0,1)表示均匀分布在[0,1]之间的随机数。
在得到初始种群之后，对其进行变异(Mutation)、交叉(Crossover)和选择(Selection)三步操作：
(a)第一步是变异操作，DE算法根据变异策略将当前种群的目标向量Ci,G进行变异得到突变向量Vi,G；常见的变异策略很多，此处采用经典的DE/rand/1变异策略，即：`
式中，正整数索引可以随机产生，且满足Cbest,G是第G代最优的个体，常数缩放因子F∈[0,2]，控制差分向量对变异个体的影响大小。突变向量Vi,G可以表示为:
(b)第二步是将每一对目标向量Ci,G和其相应的变异向量Vi,G进行交叉操作产生试验向量(trialvector)Ui,G；可以表示为:
DE算法中，采用二项交叉产生Ui,G,其定义如下：
式中，交叉常数CR∈[0,1]为用户自定义常数，其大小控制从变异向量Vi,G拷贝比例的大小。如果新产生的试验向量Ui,G中的某些值超过设定的上下边界，那么可以将其均匀随机地在其限定的边界内进行初始化，以符合要求。
(c)第三步需要计算所有试验向量Ui,G对应的目标函数f的值，并进行选择操作。其具体规则是将试验向量Ui,G所对应的目标函数值f(Ui,G)与当前种群中相应的目标向量Ci,G所对应的目标函数值f(Ci,G)进行比较，如果试验向量函数值f(Ui,G)小于或等于目标向量函数值f(Ci,G)，将试验向量体替代目标向量进入下一代种群。否则，目标向量仍然留在种群中进入下一代。该选择操作可表示为：
步骤4：反复进行变异、交叉和选择这三个操作，设定终止条件，搜索得到约束函数的全局极小值C。将C代入加权Zernike多项式,得到真实解包相位。
本发明采用点衍射干涉仪对一球面进行了测量，并把测量结果与ZYGO公司的GPI系列的干涉仪测量结果做了对比。在点衍射干涉测量中，图3是实验中采用5步相移法获取的5幅干涉条纹图，图4是获取的包裹相位。图5是对获取的包裹相位采用本发明的解包裹方法得到最终面型，图6为ZYGO干涉仪的测量结果，图7为本发明的解包结果与ZYGO干涉仪的解包结果之间的点对点残差。从图中可以看出，基于差分进化算法的解包结果与ZYGO干涉仪的结果残差很小。
与现有技术比较，本发明具有以下有益效果：
本发明在相位拟合时，采用加权Zernike多项式拟合的方法，这样被大噪声污染的数据点就不会影响最终的解包结果，消除了大噪声情况下拟合失真的问题，提高了该方法对噪声的鲁棒性。本发明基于差分进化算法的相位解包裹方法，该方法相比于基于微分Zernike多项式拟合的相位解包裹方法，不需要计算相位梯度和多项式的一阶导数。本发明首次将差分进化算法应用于光干涉计量领域的相位解包裹中。

Claims (7)

1.一种基于差分进化算法的相位解包裹方法，用于光干涉计量领域的相位解包裹，其特征在于：包括以下步骤：
(1)采用加权Zernike多项式拟合的方法建立实际干涉图的多项式，以剔除噪声污染的数据点，建立的多项式为：
其中，(x,y)为第i个像素的坐标，M为像素的数量，ci为多项式的系数，C＝【c1、c2、c3……cM】，W(x,y)为坐标(x,y)下的权重，Zi(x,y)为Zernike多项式；
(2)建立拟合系数的目标函数
对步骤(1)的多项式取余弦，建立多项式拟合系数C的最小化约束函数：
其中，L为干涉图有效孔径，λ为约束因子，γ为估计干涉图与实际干涉图的相关系数；
(3)利用DE算法对多项式拟合系数C的最小化约束函数进行全局优化，求解拟合系数，从而得到真实的解包相位。
2.根据权利要求1所述的一种基于差分进化算法的相位解包裹方法，其特征在于：假设相邻两个像素之间的相位之差的绝对值不超过π，则步骤(1)中的权重W(x,y)根据以下公式计算：
3.根据权利要求1所述的一种基于差分进化算法的相位解包裹方法，其特征在于：步骤(2)中的估计干涉图与实际干涉图的相关系数γ定义如下：
4.根据权利要求1所述的一种基于差分进化算法的相位解包裹方法，其特征在于：所述步骤(3)的具体步骤包括以下：
假设待优化的变量维数为D，进化代数为NP，那么进化的第G代其目标向量Ci,G为：
初始种群在初始化时均匀随机覆盖整个搜索空间，则搜索空间的上下边界分别表示为：
初始种群第i个个体的第j个参数按照如下公式产生：
其中，rand(0,1)表示均匀分布在[0,1]之间的随机数；
得到初始种群之后，对其进行变异、交叉和选择，直到搜索到约束函数的全局极小值C。
5.根据权利要求4所述的一种基于差分进化算法的相位解包裹方法，其特征在于：采用DE/rand/1变异策略将当前种群的目标向量Ci,G进行变异得到突变向量Vi,G，所述的DE/rand/1变异策略表示为：
其中，分别为第一次变异时第一个个体、第二个个体、第三个个体的目标向量，正整数索引随机产生，且满足常数缩放因子F∈[0,2]控制差分向量对变异个体的影响大小；
突变向量Vi,G表示为:
6.根据权利要求5所述的一种基于差分进化算法的相位解包裹方法，其特征在于：对初始种群进行交叉的方法为：将每一对目标向量Ci,G和其相应的突变向量Vi,G进行交叉操作产生试验向量Ui,G，该试验向量Ui,G表示为:
Ui,G定义如下：
式中，表示第j对目标向量Ci,G和其相应的突变向量Vi,G进行交叉操作产生试验向量，交叉常数CR∈[0,1]为自定义常数，其大小控制从突变向量Vi,G拷贝比例的大小；如果新产生的试验向量Ui,G中的某些值超过设定值的上下边界，那么将其均匀随机地在其限定的边界内进行初始化，以符合要求。
7.根据权利要求6所述的一种基于差分进化算法的相位解包裹方法，其特征在于：计算所有试验向量Ui,G对应的目标函数f的值，并进行选择操作的具体方法为：将试验向量Ui,G所对应的目标函数值f(Ui,G)与当前种群中相应的目标向量Ci,G所对应的目标函数值f(Ci,G)进行比较，如果试验向量函数值f(Ui,G)小于或等于目标向量函数值f(Ci,G)，将试验向量体替代目标向量进入下一代种群，否则，目标向量仍然留在种群中进入下一代。
CN201610298044.6A 2016-05-06 2016-05-06 一种基于差分进化算法的相位解包裹方法 Active CN106017305B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610298044.6A CN106017305B (zh) 2016-05-06 2016-05-06 一种基于差分进化算法的相位解包裹方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610298044.6A CN106017305B (zh) 2016-05-06 2016-05-06 一种基于差分进化算法的相位解包裹方法

Publications (2)

Publication Number Publication Date
CN106017305A CN106017305A (zh) 2016-10-12
CN106017305B true CN106017305B (zh) 2019-05-21

Family

ID=57081931

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610298044.6A Active CN106017305B (zh) 2016-05-06 2016-05-06 一种基于差分进化算法的相位解包裹方法

Country Status (1)

Country Link
CN (1) CN106017305B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111397645A (zh) * 2020-04-06 2020-07-10 华中科技大学 一种相位差分解调方法及系统

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107590538B (zh) * 2017-08-28 2021-04-27 南京航空航天大学 一种基于在线序列学习机的危险源识别方法
CN111461224B (zh) * 2020-04-01 2022-08-16 西安交通大学 一种基于残差自编码神经网络的相位数据解包裹方法
CN112179269B (zh) * 2020-08-18 2021-08-13 西安交通大学 一种基于Zernike多项式与WKF的相位解包裹方法
CN113029042B (zh) * 2021-05-25 2021-08-03 四川大学 一种高温熔融态金属表面形貌动态测量装置及方法
CN113899305B (zh) * 2021-09-30 2023-08-22 广东技术师范大学 一种改进的相移相位测量方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1912541A (zh) * 2006-08-30 2007-02-14 天津大学 微纳结构3d轮廓测量中基于模板的相位解包裹方法
CN102012668A (zh) * 2010-09-16 2011-04-13 昆明理工大学 一种基于剪切原理的相位解包裹方法
CN102779207A (zh) * 2012-06-19 2012-11-14 北京航空航天大学 基于OpenCL的并行差分进化算法的翼型优化设计方法
KR101459534B1 (ko) * 2013-01-11 2014-11-10 연세대학교 산학협력단 저밀도 패리티 검사 부호의 차수 분포 설계 방법 및 장치
CN105516044A (zh) * 2015-12-01 2016-04-20 西南石油大学 一种基于差分进化算法的ofdm系统峰均比抑制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1912541A (zh) * 2006-08-30 2007-02-14 天津大学 微纳结构3d轮廓测量中基于模板的相位解包裹方法
CN102012668A (zh) * 2010-09-16 2011-04-13 昆明理工大学 一种基于剪切原理的相位解包裹方法
CN102779207A (zh) * 2012-06-19 2012-11-14 北京航空航天大学 基于OpenCL的并行差分进化算法的翼型优化设计方法
KR101459534B1 (ko) * 2013-01-11 2014-11-10 연세대학교 산학협력단 저밀도 패리티 검사 부호의 차수 분포 설계 방법 및 장치
CN105516044A (zh) * 2015-12-01 2016-04-20 西南石油大学 一种基于差分进化算法的ofdm系统峰均比抑制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《基于圆载频干涉图相位解调技术的点衍射干涉术》;李金鹏 等;《光学学报》;20150131;第1-9页

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111397645A (zh) * 2020-04-06 2020-07-10 华中科技大学 一种相位差分解调方法及系统
CN111397645B (zh) * 2020-04-06 2020-12-18 华中科技大学 一种相位差分解调方法及系统

Also Published As

Publication number Publication date
CN106017305A (zh) 2016-10-12

Similar Documents

Publication Publication Date Title
CN106017305B (zh) 一种基于差分进化算法的相位解包裹方法
CN111323830B (zh) 一种基于大地电磁和直流电阻率数据的联合反演方法
Gu et al. Local stereo matching with adaptive support-weight, rank transform and disparity calibration
Zhao et al. Initial guess by improved population-based intelligent algorithms for large inter-frame deformation measurement using digital image correlation
CN110245802A (zh) 基于改进梯度提升决策树的卷烟空头率预测方法及系统
CN106295880A (zh) 一种电力系统多目标无功优化的方法及系统
CN103714577B (zh) 一种适用于带纹理模型的三维模型简化方法
CN107907856A (zh) 一种基于虚拟参考标签的rfid定位方法及系统
CN109409614A (zh) 一种基于贝叶斯正则化神经网络的电力负荷预测方法
CN107832259A (zh) 一种基于时间序列和卡尔曼滤波的负荷预测方法
CN108460462A (zh) 一种基于区间参数优化的区间神经网络学习方法
CN110233484A (zh) 考虑补偿效益的多目标无功优化配置方法
CN103645632A (zh) 一种基于广义模糊 cmac 的模型参考自适应控制方法
CN108629809A (zh) 一种精确高效的立体匹配方法
CN105405100B (zh) 一种稀疏驱动sar图像重建正则化参数自动选择方法
CN109981195A (zh) 无线信号强度的处理方法及装置
Mole et al. Growing self-organizing surface map: Learning a surface topology from a point cloud
CN109698505A (zh) 大电网静态电压稳定在线防控的调控量化映射计算方法
CN113408622A (zh) 计及特征量信息表达差异的非侵入式负荷辨识方法及系统
CN110942504B (zh) 一种众核平台上面向规则网格问题的结构化着色方法
CN106054665B (zh) 一种大规模光伏逆变器系统的分群等值建模方法
CN109308360A (zh) 信道参数寻找方法、装置、存储介质和计算机设备
CN109323677A (zh) 改进布谷鸟搜索算法的圆度误差评定算法
Choi et al. Style-guided and disentangled representation for robust image-to-image translation
CN109004916A (zh) 量子状态滤波器及相关方法

Legal Events

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