CN113820360B - 一种基于正交相位编码线性调频的高分辨率光热脉冲压缩热成像检测方法 - Google Patents
一种基于正交相位编码线性调频的高分辨率光热脉冲压缩热成像检测方法 Download PDFInfo
- Publication number
- CN113820360B CN113820360B CN202110953038.0A CN202110953038A CN113820360B CN 113820360 B CN113820360 B CN 113820360B CN 202110953038 A CN202110953038 A CN 202110953038A CN 113820360 B CN113820360 B CN 113820360B
- Authority
- CN
- China
- Prior art keywords
- thermal
- signal
- pulse compression
- chirp
- resolution
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 26
- 230000006835 compression Effects 0.000 title claims abstract description 25
- 238000007906 compression Methods 0.000 title claims abstract description 25
- 238000001931 thermography Methods 0.000 title claims abstract description 25
- 230000005284 excitation Effects 0.000 claims abstract description 40
- 238000001914 filtration Methods 0.000 claims abstract description 33
- 238000012545 processing Methods 0.000 claims abstract description 7
- 238000000034 method Methods 0.000 claims abstract description 6
- 239000011158 industrial composite Substances 0.000 claims abstract description 5
- 239000000463 material Substances 0.000 claims abstract description 5
- 239000002131 composite material Substances 0.000 abstract description 11
- 230000007547 defect Effects 0.000 abstract description 3
- 238000012512 characterization method Methods 0.000 abstract description 2
- 238000012805 post-processing Methods 0.000 abstract description 2
- 238000000605 extraction Methods 0.000 abstract 1
- 229920000049 Carbon (fiber) Polymers 0.000 description 9
- 239000004917 carbon fiber Substances 0.000 description 9
- 230000000052 comparative effect Effects 0.000 description 8
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 8
- 239000004918 carbon fiber reinforced polymer Substances 0.000 description 7
- 210000001519 tissue Anatomy 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- 230000002401 inhibitory effect Effects 0.000 description 2
- 238000009659 non-destructive testing Methods 0.000 description 2
- 241001494479 Pecora Species 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000003365 glass fiber Substances 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 239000012720 thermal barrier coating Substances 0.000 description 1
- 230000003685 thermal hair damage Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/72—Investigating presence of flaws
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/20—Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J2005/0077—Imaging
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明提供一种基于正交相位编码线性调频的高分辨率光热脉冲压缩热成像检测方法,包括以下步骤:步骤10:将瞬时频率在单个子脉冲内呈线性变化,且在单个子脉冲间进行二元相位调制信号作为激励信号,发射所述激励信号对待测样品进行加热;步骤20:利用红外热像仪获取待测样品表面的热回波信号;步骤30:将热波回波信号进行匹配滤波处理,得到匹配滤波输出信号并利用先进后处理算法进行特征提取及缺陷定量化表征。本发明方法结合了单一的线性调频信号与单一的二进制相位编码信号的优点,具有较强的抗干扰性能和较好的脉冲压缩质量,可应用于复合材料、生物组织的高分辨率光热无损检测,因此对于工业复合材料、生物医学领域具有重要意义。
Description
技术领域
本发明涉及多物理场无损检测技术领域,尤其涉及一种基于正交相位编码线性调频的高分辨率光热脉冲压缩热成像检测方法。
背景技术
脉冲压缩热成像技术由于具有高信噪比、大的动态探测深度,已被用于碳纤维等工业复合材料、牙齿等生物组织的无损检测,该技术即使在仅使用低功率外部激励源情形下,也能显著提高信噪比和增大热成像的探测范围/深度分辨率,因此通常不会对待测样品表面产生热损伤。然而,距离/深度分辨率是制约脉冲压缩热成像技术发展的主要瓶颈,目前,可通过施加合适的激励波形和采用合适的后处理算法来进一步提高脉冲压缩热成像的分辨率,此外,所施加激励波形能够在实际应用中简易实施也是所要考虑的重要因素。因此,探索新的激励波形对于实现高分辨率脉冲压缩热成像具有重要的应用价值。
目前,脉冲压缩热成像技术中通常采用线性调频激励波形和相位调制巴克码激励波形。线性调频激励波形可实现动态探测,但其匹配滤波输出通常具有较高的旁瓣,不利于高分辨率脉冲压缩热成像;巴克码相位调制激励波形具有很好的抗噪性能,但其探测深度范围通常比较有限。
发明内容
本发明目的在于针对现有技术的缺陷,提供一种基于正交相位编码线性调频的高分辨率光热脉冲压缩热成像检测方法。
为解决上述技术问题,本发明提供技术方案如下:
一种基于正交相位编码线性调频的高分辨率光热脉冲压缩热成像检测方法,其特征在于:包括如下步骤:
步骤10:设计瞬时频率在单个子脉冲内呈线性变化,且在单个子脉冲间进行二元相位调制的正交相位编码线性调频信号作为激励光源的发射波形,并将其通过信号发生器向待测样品发射;
步骤20:获取待测样品表面的热回波信号;
步骤30:对所述步骤20中所获取的待测样品表面热回波信号进行匹配滤波处理,得到待测样品表面热回波信号的匹配滤波输出信号;
步骤40:根据所述步骤30中得到的匹配滤波输出信号,得到高分辨率的光热脉冲压缩热成像。
进一步的,在所述步骤30和步骤40之间还包括步骤35:利用窗函数抑制待测样品表面热回波信号的匹配滤波输出信号的旁瓣,以得到旁瓣更低的匹配滤波输出信号。
进一步的,所述步骤10中的正交相位编码线性调频激励波形信号s(t)为:
式中,fOPCLFM(t)=f1+(f2-f1)/T′·t表示瞬时频率,k=(f2-f1)/T′表示频率调制斜率,f1表示起始频率,f2表示终止频率,T′表示激励持续时间,C为二进制相位编码,m为单个子脉冲的数目,P表示码元系数;其中,m取13,此时C=[1,1,1,1,1,-1,-1,1,1,-1,1,-1,1],fc为中心频率。
进一步的,根据权利要求3所述一种基于正交相位编码线性调频的高分辨率光热脉冲压缩热成像检测方法,其特征在于:所述步骤30的匹配滤波输出信号R(τ)为:
式中,T为热回波信号,τ为时间延迟。
进一步的,所述步骤35中的窗函数为凯瑟窗函数。
进一步的,所述步骤10中的待测样品为工业复合材料或生物组织。
与现有技术相比,本发明的有益效果是:1、本发明的匹配滤波输出由于具有较低的旁瓣和较窄的主峰,结合了单一的线性调频信号与单一的二进制相位编码信号的优点,具有很好的抗噪性能,其探测深度范围较深,因此展现出很高的信噪比和深度分辨率。2、凯瑟(Kaiser)窗函数的引入可进一步抑制其匹配滤波输出的旁瓣,可实现高分辨率光热脉冲压缩热成像。3、本发明具有普适性,不仅可应用于碳纤维复合材料的光热无损检测,还可用于牙齿等生物组织的高分辨率光热无损表征。
附图说明
图1为本发明步骤原理图;
图2为碳纤维增强复合板的正交相位编码线性调频激励波形与理想状态下的热回波信号波形图;
图3为碳纤维增强复合板表面正交相位编码线性调频激励热回波信号的不加窗函数和加窗函数匹配滤波输出对比图;
图4为图3的局部放大图;
图5为碳纤维增强复合板的线性调频激励波形与理想状态下的热回波信号;
图6为碳纤维增强复合板表面线性调频激励热回波信号的不加窗函数和加窗函数匹配滤波输出对比图;
图7为图6的局部放大图;
图8为碳纤维增强复合板的巴克码相位调制激励波形与理想状态下的热回波信号;
图9为碳纤维增强复合板表面巴克码相位调制激励热回波信号的不加窗函数和加窗函数匹配滤波输出对比图;
图10为图9的局部放大图。
具体实施方式
为了加深本发明的理解,下面我们将结合附图对本发明作进一步详述,该实施例仅用于解释本发明,并不构成对本发明保护范围的限定。
如图1所示,一种基于正交相位编码线性调频的高分辨率光热脉冲压缩热成像检测方法,包括如下步骤:
步骤10:设计瞬时频率在单个子脉冲内呈线性变化,且在单个子脉冲间进行二元相位调制的正交相位编码线性调频信号作为激励光源的发射波形,并将其通过信号发生器向待测样品发射,该激励波形信号s(t)为:
式中,fOPCLFM(t)=f1+(f2-f1)/T′·t表示瞬时频率,k=(f2-f1)/T′表示频率调制斜率,f1表示起始频率,f2表示终止频率,T′表示激励持续时间,C为二进制相位编码,m为单个子脉冲的数目,P表示码元系数;其中,m取13,此时C=[1,1,1,1,1,-1,-1,1,1,-1,1,-1,1],fc为中心频率;
步骤20:获取待测样品表面的热回波信号;
步骤35:利用合适的窗函数进一步抑制待测样品表面热回波信号的匹配滤波输出信号的旁瓣,以得到旁瓣更低的匹配滤波输出信号。
步骤40:根据步骤35中得到的匹配滤波输出信号,得到高分辨率的光热脉冲压缩热成像。
优选地,步骤35中的窗函数为凯瑟窗函数;步骤10中的待测样品既可以为热障涂层类多层复合结构或玻璃纤维复合材料等工业复合材料,也可以是生物医学领域的牙齿、羊骨等生物组织。
下面结合具体的实施例进一步说明本发明的方法:
本实施例研究对象为单层厚度为0.2mm,总厚度为3.2mm的碳纤维增强聚合板,其热扩散率为α=6.16×10-7m2/s,热传导率为k=0.46W/mK。
步骤10:为了与13位巴克码激励波形相对应,中心频率fc=0.1Hz被定义,因此,线性调频波形与正交相位编码线性调频激励波形的起始频率f1=0.1×(1-0.33)Hz,终止频率f2=0.1×(1+0.33)Hz,激励持续时间T′=130s。利用式(1)设计瞬时频率在单个子脉冲内呈线性变化,且在单个子脉冲间进行二元相位调制的正交相位编码线性调频激励波形,如图2中实线所示,所获得的理想状态下的无噪声热回波信号如图2中虚线所示。
其中,fOPCLFM(t)=f1+(f2-f1)/T′·t表示瞬时频率,k=(f2-f1)/T′表示频率调制斜率,C为二进制相位编码,m为单个子脉冲的数目,P表示码元系数,其中m取13,此时,C=[1,11,1,1,-1,-1,1,1,-1,1,-1,1],
步骤20:获取碳纤维增强聚合板表面的热回波信号。
步骤30:利用式(2)对碳纤维增强聚合板表面热回波信号进行匹配滤波处理,得到碳纤维增强聚合板表面热回波信号的匹配滤波输出信号R(τ),如图3和图4中虚线所示。
步骤35:利用凯瑟(Kaiser)窗函数进一步抑制热回波信号的匹配滤波输出的旁瓣,以得到旁瓣更低的匹配滤波输出波形,如图3和图4中实线所示,其旁瓣水平由-47.52dB降为-49.56dB。
为了能够更好的说明上述实施例的优点,现设置以下两组对比例。
对比例1:
步骤1:将线性调频信号作为激励信号,向碳纤维增强聚合板发射激励信号加热;激励信号的波形如图5中实线所示;该线性调频激励波形的起始频率f1=0.1×(1-0.33)Hz,终止频率f2=0.1×(1+0.33)Hz,激励持续时间T′=130s;所获得的理想状态下的无噪声热回波信号如图5中虚线所示。
步骤2:利用式(2)将热回波信号进行匹配滤波处理,得到匹配滤波输出信号,如图6和图7中虚线所示。
步骤3:利用凯瑟(Kaiser)窗函数进一步抑制热回波信号的匹配滤波输出的旁瓣,以得到旁瓣更低的匹配滤波输出波形,如图6和图7中实线所示,其旁瓣水平由-30.56dB降为-31.17dB。
对比例2:
步骤1:将中心频率f=0.1Hz的13位相位调制巴克码波形作为激励信号,向碳纤维增强聚合板表面发射激励信号,对碳纤维增强聚合板进行加热;激励信号的波形如图8中实线所示,所获得的理想状态下的无噪声热回波信号如图8中虚线所示。
步骤2:利用式:2)将热回波信号进行匹配滤波处理,得到匹配滤波输出信号,如图9和图10中虚线所示。
步骤3:利用凯瑟窗函数进一步抑制热回波信号的匹配滤波输出的旁瓣,以得到旁瓣更低的匹配滤波输出波形,如图9和图10中实线所示,其旁瓣水平由-50.21dB降为-52.51dB。
对比实施例与对比例1:采用线性调频激励波形得到的热回波信号波形的旁瓣水平为-30.56dB,即使利用凯瑟窗函数进一步抑制热回波信号的匹配滤波输出的旁瓣后仍为-31.17dB,效果不是很理想;而采用本发明的方法得到的热回波信号波形的旁瓣水平为-47.52dB,进一步利用凯瑟窗函数后,其旁瓣水平再降低为-49.56dB,因此,本发明相比对比例1,其热回波信号的旁瓣明显更低,具有更强的抗噪性能。
对比实施例与对比例2:虽然采用巴克码相位调制激励波形的热回波信号相比本发明具有更好的旁瓣水平,但是巴克码相位调制激励波形对待测样品的探测深度范围有限,需要进行多次单频检测;而本发明波形可动态调制热扩散深度,因此可避免巴克码波形的多次单频检测,实现碳纤维增强复合板缺陷的一次性可靠检测,大大提高检测效率。
上述具体实施方式,仅为说明本发明的技术构思和结构特征,目的在于让熟悉此项技术的相关人士能够据以实施,但以上内容并不限制本发明的保护范围,凡是依据本发明的精神实质所作的任何等效变化或修饰,均应落入本发明的保护范围之内。
Claims (6)
1.一种基于正交相位编码线性调频的高分辨率光热脉冲压缩热成像检测方法,其特征在于:包括如下步骤:
步骤10:设计瞬时频率在单个子脉冲内呈线性变化,且在单个子脉冲间进行二元相位调制的正交相位编码线性调频信号作为激励光源的发射波形,并将其通过信号发生器向待测样品发射;
步骤20:获取待测样品表面的热回波信号;
步骤30:对所述步骤20中所获取的待测样品表面热回波信号进行匹配滤波处理,得到待测样品表面热回波信号的匹配滤波输出信号;
步骤40:根据所述步骤30中得到的匹配滤波输出信号,得到高分辨率的光热脉冲压缩热成像。
2.根据权利要求1所述一种基于正交相位编码线性调频的高分辨率光热脉冲压缩热成像检测方法,其特征在于:在所述步骤30和步骤40之间还包括步骤35,
步骤35:利用窗函数抑制待测样品表面热回波信号的匹配滤波输出信号的旁瓣,以得到旁瓣更低的匹配滤波输出信号。
5.根据权利要求4所述一种基于正交相位编码线性调频的高分辨率光热脉冲压缩热成像检测方法,其特征在于:所述步骤35中的窗函数为凯瑟窗函数。
6.根据权利要求5所述一种基于正交相位编码线性调频的高分辨率光热脉冲压缩热成像检测方法,其特征在于:所述步骤10中的待测样品为工业复合材料或生物组织。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110953038.0A CN113820360B (zh) | 2021-08-19 | 2021-08-19 | 一种基于正交相位编码线性调频的高分辨率光热脉冲压缩热成像检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110953038.0A CN113820360B (zh) | 2021-08-19 | 2021-08-19 | 一种基于正交相位编码线性调频的高分辨率光热脉冲压缩热成像检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113820360A CN113820360A (zh) | 2021-12-21 |
CN113820360B true CN113820360B (zh) | 2022-12-27 |
Family
ID=78913274
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110953038.0A Active CN113820360B (zh) | 2021-08-19 | 2021-08-19 | 一种基于正交相位编码线性调频的高分辨率光热脉冲压缩热成像检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113820360B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102011118697A1 (de) * | 2011-11-16 | 2013-05-16 | Carl Zeiss Optronics Gmbh | Bilderfassungssystem |
CN103901073A (zh) * | 2014-04-22 | 2014-07-02 | 哈尔滨工业大学 | 一种基于移相调频的光热成像检测方法 |
CN104535615A (zh) * | 2015-01-18 | 2015-04-22 | 何赟泽 | 一种互相关涡流热成像缺陷检测和层析成像方法及系统 |
CN104677944A (zh) * | 2015-03-25 | 2015-06-03 | 何赟泽 | 一种微波调频热波成像系统及方法 |
CN104713906A (zh) * | 2015-04-01 | 2015-06-17 | 何赟泽 | 一种微波锁相热成像系统及方法 |
CN106940332A (zh) * | 2017-03-31 | 2017-07-11 | 哈尔滨工业大学 | 一种基于相位编码调制方式的红外无损检测系统及方法 |
CN107490599A (zh) * | 2017-09-29 | 2017-12-19 | 电子科技大学 | 一种铅钢多层材料脱粘缺陷周期脉冲热成像检测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9638648B2 (en) * | 2012-03-29 | 2017-05-02 | General Electric Company | Flaw detection using transient thermography |
-
2021
- 2021-08-19 CN CN202110953038.0A patent/CN113820360B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102011118697A1 (de) * | 2011-11-16 | 2013-05-16 | Carl Zeiss Optronics Gmbh | Bilderfassungssystem |
CN103901073A (zh) * | 2014-04-22 | 2014-07-02 | 哈尔滨工业大学 | 一种基于移相调频的光热成像检测方法 |
CN104535615A (zh) * | 2015-01-18 | 2015-04-22 | 何赟泽 | 一种互相关涡流热成像缺陷检测和层析成像方法及系统 |
CN104677944A (zh) * | 2015-03-25 | 2015-06-03 | 何赟泽 | 一种微波调频热波成像系统及方法 |
CN104713906A (zh) * | 2015-04-01 | 2015-06-17 | 何赟泽 | 一种微波锁相热成像系统及方法 |
CN106940332A (zh) * | 2017-03-31 | 2017-07-11 | 哈尔滨工业大学 | 一种基于相位编码调制方式的红外无损检测系统及方法 |
CN107490599A (zh) * | 2017-09-29 | 2017-12-19 | 电子科技大学 | 一种铅钢多层材料脱粘缺陷周期脉冲热成像检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113820360A (zh) | 2021-12-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Misaridis et al. | Potential of coded excitation in medical ultrasound imaging | |
Smith et al. | Width-modulated square-wave pulses for ultrasound applications | |
CN1218174C (zh) | 用来对超声散射体进行成象的系统和方法 | |
Garcia-Rodriguez et al. | Application of Golay codes to improve the dynamic range in ultrasonic Lamb waves air-coupled systems | |
CN1210003C (zh) | 基于脉冲压缩的超声波收发装置 | |
CN101496728B (zh) | 超声频率复合成像方法与装置 | |
CN113820360B (zh) | 一种基于正交相位编码线性调频的高分辨率光热脉冲压缩热成像检测方法 | |
US20180156907A1 (en) | Echo measurement | |
CN105628795B (zh) | 利用频率相位复合编码激励信号进行焊缝超声检测的方法 | |
CN105866257B (zh) | 一种基于Golay互补卷积码的多元调制超声编码单次激励方法 | |
CN109639345B (zh) | 一种基于时域反射法tdr技术的线缆带宽测试方法 | |
Nowicki et al. | On the application of signal compression using Golay's codes sequences in ultrasound diagnostic | |
Nowicki et al. | Estimation of ultrasonic attenuation in a bone using coded excitation | |
Yichun et al. | Optimization design of NLFM signal and its pulse compression simulation | |
Svilainis et al. | Ultrasound transmission spectral compensation using arbitrary position and width pulse sets | |
White et al. | Pseudo random binary sequence techniques for ultrasonic NDE of high loss materials | |
Zanelli et al. | A robust hydrophone for HIFU metrology | |
Fan et al. | Coded excitation with Nonlinear Frequency Modulation Carrier in Ultrasound Imaging System | |
CN113358697B (zh) | 基于非线性调频的高分辨率光热脉冲压缩热成像检测方法 | |
Fu et al. | Barker coded excitation using LFM carrier for improving axial resolution in ultrasound imaging | |
Lahav et al. | FoCUS: Fourier-based coded ultrasound | |
Svilainis et al. | Application of arbitrary position and width pulse trains signals in ultrasonic imaging: correlation performance study | |
Litniewski et al. | Sound fields for coded excitations in water and tissue: experimental approach | |
CN110954907A (zh) | 一种应用于深海浅表的高精度沉积物声学探测装置 | |
Raman et al. | Pre-enhancement of chirp signal for inverse filtering in medical ultrasound |
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 |