CN105890538A - 三表面干涉式高精度曲面轮廓测量系统及方法 - Google Patents
三表面干涉式高精度曲面轮廓测量系统及方法 Download PDFInfo
- Publication number
- CN105890538A CN105890538A CN201410856974.XA CN201410856974A CN105890538A CN 105890538 A CN105890538 A CN 105890538A CN 201410856974 A CN201410856974 A CN 201410856974A CN 105890538 A CN105890538 A CN 105890538A
- Authority
- CN
- China
- Prior art keywords
- interference
- wedge
- wave number
- centerdot
- curved surface
- 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.)
- Pending
Links
Landscapes
- Length Measuring Devices By Optical Means (AREA)
- Instruments For Measurement Of Length By Optical Means (AREA)
Abstract
本发明公开了一种三表面干涉式高精度曲面轮廓测量系统及方法。这是一种利用迈克尔逊干涉原理,获取被测零件表面三维轮廓信息的测量系统。首先用计算机控制半导体激光器进行波数扫描,与此同时CCD相机连续拍摄不同波数下的干涉图像。通过光路中光楔前后表面的干涉图像进行波数在线监测。然后将干涉图像每个像素沿时间轴傅立叶变换,在被测零件曲面轮廓和光楔前表面干涉信号峰值处,提取卷绕相位信息。解卷绕后,得到被测零件表面三维轮廓。本发明的曲面三维轮廓测量精度为±10nm,稳定可靠,不需经常实时效验,不需参考曲面,同时保持很高的测量精度。
Description
技术领域
本发明提出一种三表面干涉式高精度曲面轮廓测量系统及方法,适用于精密测试、数控设备等领域,属于物体表面三维轮廓测量领域。
背景技术
高精度曲面轮廓测量技术在光学零件制造、模具加工等领域有着广泛的应用。目前,这种技术主要基于对比测量法原理,测量精度可达到亚微米级,但是需要一个高精度对比测量参考面,而这对于用户来讲,不易满足,且测量成本较高。
波数扫描干涉测量技术是以激光波数扫描干涉方法为基础进行测量的一门技术。它利用CCD相机和计算机,将多幅干涉图像进行快速处理,获取被测零件的信息。其非接触式的测量方法和基于傅里叶变换的频率域数据解调方法确保了系统的稳定性和测量的超高精度,可用于长度、角度、轮廓、位移、介质折射率变化和振动等方面的几何量测量。
本专利设计了一种三表面干涉式高精度曲面轮廓测量系统。此系统基于迈克尔逊干涉原理,首先利用可调谐激光进行波数扫描,与此同时CCD相机拍摄多幅三表面干涉图像。接着通过傅立叶变换、相位解卷绕技术处理获得被测零件的表面轮廓信息。该系统具有非接触、超高精度、可靠性高和成本低等特点,特别适用于机械加工后的质量检验应用中。
发明内容
本发明提供了一种三表面干涉式高精度曲面轮廓测量系统及方法。
本发明通过如下技术方案实现:
三表面干涉式高精度曲面轮廓测量系统如图1所示,包括激光控制器(1)、半导体激光器(2)、准直透镜(3)、分光镜(4)、偏振片(5)、被测零件(6)、光楔(7)、CCD相机(8)、计算机(9)。
测量过程如下:先将被测零件(6)固定,然后打开激光器(2),利用激光控制器(1)进行波数扫描,其扫描时间、扫描步距等由计算机(9)发出指令控制。激光器(2)发出的激光被分光镜(4)分为两束,其中一束光经过偏振片(5),在被测零件(6)表面产生反射光;另一束光在光楔(7)前后表面也产生反射光。这三个表面的反射光经分光镜(4)的透射和反射作用后,相互叠加形成干涉图像,由CCD相机(8)实时拍摄,最后在计算机(9)上将存储的干涉图像数据进行分析处理。
三表面干涉式高精度曲面轮廓测量系统数据处理方法的具体步骤如下:
1)激光波数扫描:激光控制器(1)通过调节半导体激光器(2)壳体的工作温度,对半导体激光器(2)的输出光进行波数扫描,使输出光的波数k跟随时间t变化,其函数关系可表示为:
其中Δk为波数扫描最大变化量;T为波数扫描周期;k0为起始波数。
2)三表面干涉光强叠加的相位差:CCD相机(8)采集到一组干涉条纹图像,其光强为:
φpq(x,y)=2·Λpq(x,y)·k0, (4)
Λpq(x,y)=npq·zpq(x,y), (5)
φpq,zpq和npq分别为面p与面q(p,q=1,2,3)反射光之间的相位差,距离和折射率;fpq为面p与面q反射光之间的干涉频率;Λpq为面p与面q之间的光程差;(x,y)为被测零件表面轮廓的空间坐标,如图1所示。
3)半导体激光器(2)输出光的波数在线监测:如图1所示,三个干涉面分别为光楔前后表面S1、S2和被测曲面S3。其中光楔前表面S1和后表面S2彼此线性倾斜,用于监测半导体激光器输出光的波数。单幅三表面干涉图像经二维傅立叶变换后,得到幅频和相频图。对于每幅干涉图像均作上述处理,然后在非直流分量的峰值频率处提取相角,得到对应每个波数下光楔前后表面干涉的卷绕相位。通过解卷绕运算,得到每个波数下光楔干涉信号的解卷绕相位φuw 12(t),离散化后为φuw 12(n),得到波数扫描干涉的波数序列k(n):
其中Λ120为光楔前后表面在(x=0mm,y=0mm)点的光程,N为CCD相机在波数扫描中拍摄干涉图像的总数目。
4)利用随机采样离散傅立叶变换获取被测零件表面的三维轮廓:CCD相机上拍摄的三表面干涉图像光强序列的随机采样离散傅里叶变换为
其中w()为波数扫描过程中采样的窗口函数。
将图像各个像素点的干涉光强序列进行随机采样离散傅立叶变换后,幅频特性上有3个峰值,分别对应光楔前后表面S1和S2干涉、S1和被测曲面S3干涉、S2和被测曲面S3干涉。由于光楔前表面S1是平整的,其平整度<λ/20,因此各个像素点的光楔前表面S1和被测曲面S3干涉峰值对应的卷绕相位φ13(x,y)包含被测曲面的高度信息。将光楔前表面S1和被测曲面S3干涉卷绕相位解卷绕后,就可得到被测曲面S3的三维轮廓信息:
其中φuw 13(x,y)为光楔前表面S1和被测曲面S3干涉的解卷绕相位图;z13(x,y)为S1和S3之间的高度差。由于光楔前表面S1平整且严格垂直z方向,Δz13(x,y)=Δz3(x,y),后者就是被测零件的曲面三维轮廓。
附图说明
图1三表面干涉式高精度曲面轮廓测量系统;包括激光控制器(1)、半导体激光器(2)、准直透镜(3)、分光镜(4)、偏振片(5)、被测零件(6)、光楔(7)、CCD相机(8)、计算机(9)。
图2三表面干涉图像经二维傅立叶变换后的幅频特性图。
图3三表面干涉图像经二维傅立叶变换后的相频特性图。
图4波数扫描对应的光楔卷绕相位序列。
图5波数扫描对应的光楔解卷绕相位序列。
图6干涉图像在x=0,y=0处的干涉光强序列。
图7光楔前表面S1和被测曲面S3干涉的卷绕相位图。
图8被测曲面S3的三维轮廓。
具体实施方式
下面结合实例和附图对本发明作进一步说明,但不应限制本发明的保护范围。
激光控制器(1)使用ILX Lightwave LDC-3724,半导体激光器(2)使用Toptica公司的DFB半导体激光器LD-0855-0150-DFB-1,中心波数7.39×106m-1,在波数扫描过程中,波数无模跳,扫描范围为4.13×103m-1。光楔(7)中心厚度6mm,倾角6′,CCD相机(8)使用PCO公司的PCO1600。计算机(9)使用Windows XP 32位系统,2G内存。三表面干涉式高精度曲面轮廓测量系统搭建在大恒光电公司的防振平台DH-OTB-1200-800-100上。
被测零件(6)为苹果iPad Air的LOGO,将被测零件安装固定以后,计算机控制半导体激光器进行波数扫描,同时CCD相机拍摄N张干涉图像,如图1所示,三个干涉面分别为光楔前后面S1、S2和被测曲面S3。
单幅三表面干涉图像经二维傅立叶变换后,幅频和相频特性图如图2,3所示。对于每幅干涉图像,在非直流分量的峰值频率处(图2中点①)提取相角,得到对应每个激光波数下光楔前后表面干涉的卷绕相位,如图4所示。通过解卷绕运算,得到每个激光波数下光楔干涉信号的解卷绕相位,离散化后为φuw 12(n),如图5所示。
图6为x=0,y=0像素点的光强序列,灰度值0~255表示光强。将干涉图像各个像素点的光强序列进行随机采样离散傅立叶变换后,幅频特性上有3个峰值,分别对应光楔前后表面S1和S2干涉、S1和被测曲面S3干涉、S2和S3干涉。提取光楔前表面S1和被测曲面S3干涉峰值处的卷绕相位,如图7所示。将上述卷绕相位解卷绕后,根据公式(8)求出被测曲面S3的三维轮廓信息,如图8所示。
本发明的优点在于:
(1)曲面三维轮廓测量精度为±10nm。
(2)系统稳定可靠,鲁棒性好。
(3)不需经常实时效验,不需参考曲面,对不同曲面轮廓仍保持很高的测量精度。
综上所述,在三维曲面轮廓测量领域,本发明所提出的三表面干涉式高精度曲面轮廓测量系统及方法具有广阔的应用前景。
Claims (6)
1.一种三表面干涉式高精度曲面轮廓测量系统,包括激光控制器(1)、半导体激光器(2)、准直透镜(3)、分光镜(4)、偏振片(5)、被测零件(6)、光楔(7)、CCD相机(8)、计算机(9)。
2.根据权利要求1中所述的三表面干涉式高精度曲面轮廓测量系统,其特征在于:激光控制器(1)控制半导体激光器(2)壳体的温度线性变化,使半导体激光器(2)输出光的波数单调变化。
3.根据权利要求1中所述的三表面干涉式高精度曲面轮廓测量系统,其特征在于:光路为迈克尔逊干涉式,一个臂安装被测零件(6),另一臂安装光楔(7)。
4.根据权利要求1、2、3中所述的三表面干涉式高精度曲面轮廓测量系统,其特征在于:光楔(7)的前表面为迈克尔逊干涉的参考面。
5.根据权利要求1、2、3中所述的三表面干涉式高精度曲面轮廓测量系统,其特征在于:光楔(7)的前后表面用于实时监测半导体激光器(2)输出光的波数。
6.一种三表面干涉式高精度曲面轮廓测量方法,其具体步骤如下:
1)激光波数扫描:激光控制器(1)通过调节半导体激光器(2)壳体的工作温度,对半导体激光器(2)输出光进行波数扫描,使输出光的波数k跟随时间t变化,其函数关系可表示为:
其中Δk为波数扫描最大变化量;T为波数扫描周期;k0为起始波数。
2)三表面干涉光强叠加的相位差:CCD相机(8)采集到一组干涉条纹图像,其光强可表示为:
φpq(x,y)=2·Λpq(x,y)·k0
Λpq(x,y)=npq·zpq(x,y)
φpq,zpq和npq分别为面p与面q(p,q=1,2,3)反射光之间的相位差,距离和折射率;fpq为面p与面q反射光之间的干涉频率;Λpq为面p与面q之间的光程差;(x,y)为被测零件表面轮廓的空间坐标。
3)半导体激光器(2)输出光的波数在线监测:三个干涉面分别为光楔(7)的前后表面S1、S2和被测曲面S3。其中光楔(7)的前表面S1和后表面S2彼此线性倾斜,用于监测半导体激光器(2)输出光的波数。单幅三表面干涉图像经二维傅立叶变换后,得到幅频和相频图;对于每幅干涉图像均作上述处理,在非直流分量的峰值频率处提取相角,得到对应每个波数下光楔(7)前后表面干涉的卷绕相位;通过解卷绕运算,得到每个波数下光楔(7)干涉信号的解卷绕相位φuw 12(t),离散化后为φuw 12(n),得到波数扫描干涉的波数序列k(n):
其中Λ120为光楔(7)前后表面在x=0mm,y=0mm处的光程,N为CCD相机(8)在波数扫描中拍摄图像的总数目。
4)利用随机采样离散傅立叶变换获取被测零件(6)表面的三维轮廓:CCD相机(8)上拍摄的三表面干涉图像光强序列的随机采样离散傅里叶变换为
其中w()为波数扫描过程中采样的窗口函数。
将图像各个像素点的干涉光强序列进行随机采样离散傅立叶变换后,幅频特性上有3个峰值,分别对应为光楔(7)前后表面S1和S2干涉、S1和被测曲面S3干涉、S2和S3干涉。由于光楔(7)的前表面S1是平整的,其平整度<λ/20,因此图像各个像素点的光楔(7)前表面S1和被测曲面S3干涉峰值对应的卷绕相位φ13(x,y)包含被测曲面S3高度信息。将光楔(7)前表面S1和被测曲面S3干涉的卷绕相位解卷绕后,就可得到被测曲面S3的三维轮廓信息:
其中φuw 13(x,y)为光楔(7)前表面S1和被测曲面S3干涉的解卷绕相位图;z13(x,y)为光楔(7)前表面S1和被测曲面S3之间的高度差。由于光楔(7)的前表面S1平整且严格垂直z方向,Δz13(x,y)=Δz3(x,y),后者就是被测零件(6)的曲面三维轮廓。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410856974.XA CN105890538A (zh) | 2014-12-30 | 2014-12-30 | 三表面干涉式高精度曲面轮廓测量系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410856974.XA CN105890538A (zh) | 2014-12-30 | 2014-12-30 | 三表面干涉式高精度曲面轮廓测量系统及方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105890538A true CN105890538A (zh) | 2016-08-24 |
Family
ID=56700945
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410856974.XA Pending CN105890538A (zh) | 2014-12-30 | 2014-12-30 | 三表面干涉式高精度曲面轮廓测量系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105890538A (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108917641A (zh) * | 2018-05-15 | 2018-11-30 | 广东工业大学 | 基于激光器波数合成的样件内部轮廓检测方法及系统 |
CN109990729A (zh) * | 2019-01-15 | 2019-07-09 | 广东工业大学 | 一种基于电流波数扫描的透镜三维轮廓测量装置及方法 |
CN110011178A (zh) * | 2019-03-28 | 2019-07-12 | 广东工业大学 | 电流-温控半导体激光器的光学干涉光源装置及其构成的测量系统 |
CN110260812A (zh) * | 2019-05-15 | 2019-09-20 | 广东工业大学 | 一种基于欠定盲源分离双通道光学三维干涉方法及系统 |
CN110658781A (zh) * | 2019-09-25 | 2020-01-07 | 华南理工大学 | 一种基于xyz轴滑台的曲面力控制轮廓跟随装置及方法 |
CN112066909A (zh) * | 2020-08-24 | 2020-12-11 | 南京理工大学 | 一种基于倾斜平面高精度提取的抗振动干涉测量方法 |
CN112762817A (zh) * | 2020-12-23 | 2021-05-07 | 广东工业大学 | 倾斜斐索波数扫描干涉仪 |
CN113237432A (zh) * | 2021-04-19 | 2021-08-10 | 武汉光谷航天三江激光产业技术研究院有限公司 | 激光焊接熔深提取方法及设备 |
CN113251945A (zh) * | 2021-05-17 | 2021-08-13 | 东北大学秦皇岛分校 | 一种线轮廓成像装置的解调方法及成像装置 |
CN113654482A (zh) * | 2021-08-30 | 2021-11-16 | 东北大学秦皇岛分校 | 一种基于色差和谱域干涉的光学3d成像装置及方法 |
CN118067017A (zh) * | 2023-12-15 | 2024-05-24 | 惠然科技有限公司 | 用于检测半导体高度的方法、高度检测系统及相关产品 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1724550B1 (en) * | 2005-05-18 | 2009-01-07 | Mitutoyo Corporation | Interferometer and shape measuring method |
KR20090025971A (ko) * | 2007-09-07 | 2009-03-11 | 한국표준과학연구원 | 형상 측정장치 및 그 방법 |
CN101788275A (zh) * | 2010-02-01 | 2010-07-28 | 天津大学 | 利用波长作相移获取表面三维形貌的方法 |
CN103267494A (zh) * | 2013-05-20 | 2013-08-28 | 湖北工业大学 | 一种表面形貌干涉测量的方法及装置 |
CN103528524A (zh) * | 2013-09-16 | 2014-01-22 | 广东工业大学 | 透视测量树脂基复合材料内部离面位移场分布的装置及方法 |
CN103727901A (zh) * | 2014-01-14 | 2014-04-16 | 中国科学院长春光学精密机械与物理研究所 | 基于波长移相法检测平面间平行度的方法 |
-
2014
- 2014-12-30 CN CN201410856974.XA patent/CN105890538A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1724550B1 (en) * | 2005-05-18 | 2009-01-07 | Mitutoyo Corporation | Interferometer and shape measuring method |
KR20090025971A (ko) * | 2007-09-07 | 2009-03-11 | 한국표준과학연구원 | 형상 측정장치 및 그 방법 |
CN101788275A (zh) * | 2010-02-01 | 2010-07-28 | 天津大学 | 利用波长作相移获取表面三维形貌的方法 |
CN103267494A (zh) * | 2013-05-20 | 2013-08-28 | 湖北工业大学 | 一种表面形貌干涉测量的方法及装置 |
CN103528524A (zh) * | 2013-09-16 | 2014-01-22 | 广东工业大学 | 透视测量树脂基复合材料内部离面位移场分布的装置及方法 |
CN103727901A (zh) * | 2014-01-14 | 2014-04-16 | 中国科学院长春光学精密机械与物理研究所 | 基于波长移相法检测平面间平行度的方法 |
Non-Patent Citations (2)
Title |
---|
任寰等: "多表面干涉下的光学元件面形检测", 《光学 精密工程》 * |
马力等: "基于波长移相调谐的多表面干涉研究", 《光学与光电技术》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108917641A (zh) * | 2018-05-15 | 2018-11-30 | 广东工业大学 | 基于激光器波数合成的样件内部轮廓检测方法及系统 |
CN109990729A (zh) * | 2019-01-15 | 2019-07-09 | 广东工业大学 | 一种基于电流波数扫描的透镜三维轮廓测量装置及方法 |
US11092511B2 (en) | 2019-01-15 | 2021-08-17 | Guangdong University Of Technology | Device and method for measuring lens contour based on laser wave number scanning |
CN110011178A (zh) * | 2019-03-28 | 2019-07-12 | 广东工业大学 | 电流-温控半导体激光器的光学干涉光源装置及其构成的测量系统 |
US11108212B2 (en) | 2019-03-28 | 2021-08-31 | Guangdong University Of Technology | Optical interference light source device of current-temperature controlled semiconductor laser and measurement system including the same |
US11060849B2 (en) | 2019-05-15 | 2021-07-13 | Guandong University Of Technology | Dual-channel optical three-dimensional interference method and system based on underdetermined blind source separation |
CN110260812B (zh) * | 2019-05-15 | 2021-02-09 | 广东工业大学 | 一种基于欠定盲源分离双通道光学三维干涉方法及系统 |
CN110260812A (zh) * | 2019-05-15 | 2019-09-20 | 广东工业大学 | 一种基于欠定盲源分离双通道光学三维干涉方法及系统 |
CN110658781A (zh) * | 2019-09-25 | 2020-01-07 | 华南理工大学 | 一种基于xyz轴滑台的曲面力控制轮廓跟随装置及方法 |
CN112066909A (zh) * | 2020-08-24 | 2020-12-11 | 南京理工大学 | 一种基于倾斜平面高精度提取的抗振动干涉测量方法 |
CN112066909B (zh) * | 2020-08-24 | 2022-04-08 | 南京理工大学 | 一种基于倾斜平面高精度提取的抗振动干涉测量方法 |
CN112762817A (zh) * | 2020-12-23 | 2021-05-07 | 广东工业大学 | 倾斜斐索波数扫描干涉仪 |
CN113237432A (zh) * | 2021-04-19 | 2021-08-10 | 武汉光谷航天三江激光产业技术研究院有限公司 | 激光焊接熔深提取方法及设备 |
CN113251945A (zh) * | 2021-05-17 | 2021-08-13 | 东北大学秦皇岛分校 | 一种线轮廓成像装置的解调方法及成像装置 |
CN113251945B (zh) * | 2021-05-17 | 2022-07-12 | 东北大学秦皇岛分校 | 一种线轮廓成像装置的解调方法及成像装置 |
CN113654482A (zh) * | 2021-08-30 | 2021-11-16 | 东北大学秦皇岛分校 | 一种基于色差和谱域干涉的光学3d成像装置及方法 |
CN118067017A (zh) * | 2023-12-15 | 2024-05-24 | 惠然科技有限公司 | 用于检测半导体高度的方法、高度检测系统及相关产品 |
CN118067017B (zh) * | 2023-12-15 | 2024-09-06 | 惠然科技有限公司 | 用于检测半导体高度的方法、高度检测系统及相关产品 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105890538A (zh) | 三表面干涉式高精度曲面轮廓测量系统及方法 | |
Huang et al. | Review of phase measuring deflectometry | |
De Groot | A review of selected topics in interferometric optical metrology | |
CN109163672B (zh) | 一种基于白光干涉零光程差位置拾取算法的微观形貌测量方法 | |
Zhang et al. | Phase measuring deflectometry for obtaining 3D shape of specular surface: a review of the state-of-the-art | |
Chen et al. | Overview of 3-D shape measurement using optical methods | |
CN102788558B (zh) | 散斑相关和散斑干涉相结合的三维变形测量系统及方法 | |
CN103528524B (zh) | 透视测量树脂基复合材料内部离面位移场分布的方法 | |
CN104729419B (zh) | 透视测量高聚物材料内部热变形场分布的装置及方法 | |
CN104296698A (zh) | 一种超高精度的光学表面平整度测量方法 | |
Huang et al. | Fast full-field out-of-plane deformation measurement using fringe reflectometry | |
CN103630086A (zh) | 一种基于单色ccd的双波长同时相移干涉测量方法 | |
Dhanasekar et al. | Digital speckle interferometry for assessment of surface roughness | |
CN105091781B (zh) | 一种单桢干涉条纹图测量光学表面的方法与装置 | |
JP2020517911A (ja) | スペクトル制御干渉法による曲率半径測定 | |
CN109059787B (zh) | 一种基于横向剪切干涉的玻璃板厚度分布测量方法及系统 | |
CN109000781A (zh) | 一种结构微振动线域测量装置及方法 | |
Huang et al. | Fast template matching method in white-light scanning interferometry for 3D micro-profile measurement | |
CN114440789A (zh) | 旋转体速度、距离和三维形貌同步干涉测量方法及系统 | |
CN106123805A (zh) | 基于白光扫描干涉的镀膜器件三维形貌测量方法 | |
Bergström et al. | Shape verification using dual-wavelength holographic interferometry | |
Carvalho et al. | Deflectometry in cylindrical coordinates using a conical mirror: principles and proof of concept | |
US20050177339A1 (en) | Precision surface measurement | |
Frade et al. | In situ 3D profilometry of rough objects with a lateral shearing interferometry range finder | |
CN103376072B (zh) | 数字全息干涉和变频率投影条纹复合测量系统及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160824 |
|
WD01 | Invention patent application deemed withdrawn after publication |