CN103257033A - 基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法 - Google Patents

基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法 Download PDF

Info

Publication number
CN103257033A
CN103257033A CN2013101784952A CN201310178495A CN103257033A CN 103257033 A CN103257033 A CN 103257033A CN 2013101784952 A CN2013101784952 A CN 2013101784952A CN 201310178495 A CN201310178495 A CN 201310178495A CN 103257033 A CN103257033 A CN 103257033A
Authority
CN
China
Prior art keywords
sub
aperture
window glass
interferometer
scanning motion
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
Application number
CN2013101784952A
Other languages
English (en)
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.)
Institute of Optics and Electronics of CAS
Original Assignee
Institute of Optics and Electronics of CAS
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 Institute of Optics and Electronics of CAS filed Critical Institute of Optics and Electronics of CAS
Priority to CN2013101784952A priority Critical patent/CN103257033A/zh
Publication of CN103257033A publication Critical patent/CN103257033A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明涉及基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法,根据干涉仪的口径和窗口玻璃的口径计算检测窗口玻璃所需子孔径个数,调整干涉仪与平面镜到自准直使干涉仪上条纹为零条纹;将四维扫描运动平台放置在干涉仪与平面镜之间,四维扫描运动平台带动窗口玻璃扫描运动,调整窗口玻璃与平面镜使干涉仪上的条纹为达到最少;将窗口玻璃调整到第一个子孔径位置,干涉仪检测窗口玻璃并记录此时子孔径透射波前数据;将窗口玻璃调整到第二个子孔径位置,干涉仪检测窗口玻璃并记录此时子孔径透射波前数据,直到窗口玻璃所有子孔径波前数据被检测并记录下来;将所有子孔径波前数据通过算法拼接,得到窗口玻璃全口径透射波前。

Description

基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法
技术领域
本发明属于光学检测领域,涉及一种基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法,可用于检测大口径窗口玻璃透射波前。
背景技术
窗口玻璃连接并隔离光学系统与外界环境,广泛的应用于空间观测光学系统、红外光学系统等各种光学仪器和设备中,窗口玻璃隔离光学系统内外压差和温度变化,在野外还可以使光学系统不受雨点、沙粒等侵蚀,窗口玻璃不仅起着保护光学系统的作用,同时还参与成像,是光学系统不可或缺的重要组成部分。
由于窗口玻璃在光学系统中参与成像,窗口玻璃透射波前的好坏直接影响光学系统成像质量,因此需要对窗口玻璃透射波前进行标定和检测。对小口径窗口玻璃透射波前检测方法是采用干涉仪与窗口玻璃、平面镜组成子准直光路,通过干涉仪直接检测。随着制造和科学技术的发展,光学系统口径越来越大,窗口玻璃的口径也随之越来越大。现有干涉仪的口径有限(≤600mm),同时加工大口径、高精度平面镜也存在着价格贵、周期长的问题。当窗口玻璃的口径大于干涉仪口径时,就不能对窗口玻璃透射波前进行检测。
发明内容
针对现有仪器和方法不能检测大口径窗口玻璃透射波前问题,本发明提供一种新的检测方法,用于对大口径窗口玻璃透射波前的检测。该方法采用小口径干涉仪和小口径平面镜检测大口径窗口玻璃透射波前,降低了检测成本,拓展了干涉仪动态范围,提高了检测的分辨力。
本发明解决上述技术问题,采取基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法的技术方案是:利用由干涉仪、窗口玻璃、四维扫描运动平台、平面镜组成检测系统,并通过以下步骤完成对窗口玻璃透射波前的测量步骤如下:
步骤S1:根据干涉仪的口径和窗口玻璃的口径计算检测窗口玻璃所需子孔径个数,以其中任意一个子孔径中心坐标为原点,计算其余子孔径中心坐标,并对子孔径进行编号;
步骤S2:调整干涉仪与平面镜到自准直使干涉仪上条纹为零条纹;
步骤S3:将四维扫描运动平台放置在干涉仪与平面镜之间,将窗口玻璃安装在四维扫描运动平台上;
步骤S4:四维扫描运动平台带动窗口玻璃扫描运动,调整窗口玻璃与平面镜使干涉仪上的条纹为达到最少;
步骤S5:通过四维扫描运动平台的运动,将窗口玻璃调整到第一个子孔径位置,干涉仪检测窗口玻璃并记录此时子孔径透射波前数据;
步骤S6:通过四维扫描运动平台的运动,将窗口玻璃调整到第二个子孔径位置,干涉仪检测窗口玻璃并记录此时子孔径透射波前数据;
步骤S7:重复步骤S5、步骤S6直到窗口玻璃所有子孔径波前数据被检测并记录下来;
步骤S8:将所有子孔径波前数据通过算法拼接,得到窗口玻璃全口径透射波前,完成对窗口玻璃的透射波前检测。
本发明的优点是:
(1)、本发明基于子孔径拼接干涉检测技术,精度高、可以突破干涉仪口径的限制检测大口径窗口玻璃透射波前。
(2)、本发明采用小口径干涉仪和平面镜检测大口径窗口玻璃透射波前,降低了检测成本。
附图说明
图1为基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法流程图。
图2为基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法原理示意图。
图3为基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法的子孔径规划图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法,该方法利用干涉仪与窗口玻璃、平面镜组成自准直光路对窗口玻璃透射波前进行测量,干涉仪发出的平行光经过窗口玻璃后又经平面镜反射回干涉仪形成干涉,通过干涉对窗口玻璃的该部分子孔径波前进行检测,窗口玻璃安装在四维扫描运动平台上,由于干涉仪口径比窗口玻璃口径小,每次只能检测窗口玻璃部分透射子孔径波前,每次子孔径波前检测完成后,由四维调整运动平台将窗口玻璃移动到下一个子孔径位置,当所有子孔径波前检测完成后,通过算法将所有子孔径波前拼接,完成对整个窗口玻璃透射波前的检测。本发明突破现有干涉仪口径的限制,可以检测口径比干涉仪大的窗口玻璃透射波前。
结合附图1和图2对本发明作进一步说明,图2中,干涉仪11的口径为600mm,窗口玻璃12口径为1000mm,四维扫描运动平台13可在水平(x方向)、垂直(y方向)两方向运动,同时可以在倾斜、俯仰两方向进行角度调整,平面镜14的口径为600mm。如图1示出基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法的流程图,该方法可通过以下步骤完成对窗口玻璃2透射波前的测量:
步骤S1:根据干涉仪11口径和窗口玻璃12口径计算检测窗口玻璃12所需子孔径个数,以其中任意一个子孔径坐标为原点,计算其余子孔径坐标,并对子孔径进行编号;
其中,步骤S1中计算子孔径个数的步骤如下:设定子孔径坐标坐标是以任意一个包含窗口玻璃12有效透射波前数据的子孔径为初始子孔径,编号为第一子孔径1,以(x1,y1)作为第一子孔径1中心坐标位置,其中x1,y1的坐标点分别取0;对第一子孔径1进行位置标定。以第一子孔径1圆心(0,0)为坐标原点建立坐标系,在四维扫描运动平台左右运动为x方向,在四维扫描运动平台x方向向左运动为负,向右运动为正,四维扫描运动平台上下运动为y方向,在四维扫描运动平台y方向向上运动为正,向下运动为负;所有子孔径数据区域大小相同,覆盖的面积均为s;任意方向移动第一子孔径1的中心坐标得到第二子孔径2,第二子孔径2中心坐标为(x2,y2)与第一子孔径1相邻并且有一重叠区域,重叠区域面积大小为s2,获取s2与s之比等于18%;任意方向移动第二子孔2径得到第三子孔径3,第三子孔径3中心坐标为(x3,y3)与第二子孔径2相邻且有一重叠区域,重叠区域面积大小为s3,获取s3与s之比等于18%;依次将子孔径移到其相邻位置得到下一个子孔径,记录下一个子孔径位置坐标,并保证子该孔径与其相邻子孔径即下一个子孔径有一重叠区域,且重叠区域面积与s之比等于18%,重复上述步骤直到所有子孔径能完全覆盖窗口玻璃12口径。
图3为基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法的子孔径规划图,以任意一个子孔径中心坐标为原点(0,0),并将该子孔径标记为第一子孔径1,第二子孔径2与第一子孔径1相邻,重叠区域波前数据面积为子孔径波前数据面积的18%,可计算得到第二子孔径2中心坐标为(424,0)。第三子孔径3与第二子孔径2相邻,重叠区域波前数据面积为子孔径波前数据面的18%,可计算第三子孔径3的中心坐标(212,367)。依次可以计算其余子孔径中心坐标,第四子孔径4中心坐标(-212,367),第五子孔径5中心坐标(-424,0),第六子孔径6中心坐标(-212,-367),第七子孔径7中心坐标(212,-367);
步骤S2:调整干涉仪11与平面镜14到自准直使干涉仪11上的条纹在5根以内;其中,所述调整干涉仪11与平面镜14到自准直是使干涉仪11的条纹数小于5根,再在干涉仪11与平面镜14中间加入四维扫描运动平台13,再接着调整窗口玻璃12和平面镜14,使干涉仪11上条纹小于5根。
步骤S3:将四维扫描运动平台13放置在干涉仪11与平面镜14之间,将窗口玻璃12安装在四维扫描运动平台13上;
步骤S4:将四维扫描运动平台13带动窗口玻璃12扫描运动,调整窗口玻璃12与平面镜14使干涉仪11上的条纹为达到最少;
步骤S5:通过四维扫描运动平台13的运动,将窗口玻璃12调整到第一子孔径1位置(0,0),干涉仪11检测窗口玻璃12并记录此时第一子孔径透射波前数据w1;
步骤S6:通过四维扫描运动平台13的运动,将窗口玻璃12调整到第二子孔径2位置(424,0),干涉仪11检测窗口玻璃12并记录此时第二子孔径透射波前数据w2;
步骤7S:重复步骤S5、步骤S6直到窗口玻璃12所有子孔径波前数据被检测并记录下来,第三子孔径3到第七子孔径7波前数据分别表示为w3、w4、w5、w6、w7;
步骤S8:将所有子孔径波前数据通过算法拼接,得到窗口玻璃12透射波前,完成对窗口玻璃12的透射波前检测。
其中,所有子孔径波前数据通过子孔径拼接算法拼接(具体算法见:子孔径拼接干涉检测大口径平面波前,汪利华,吴时彬等,光电工程,2009,36(6):126-130),具体的拼接方法如下:
w1~w7为干涉仪11直接检测的子孔径波前数据,检测过程中由于四维扫描运动平台13的运动,会给真实的子孔径波前数据引入倾斜和平移调整误差,而每个子孔径倾斜和平移调整误差也都不相同,干涉仪11直接检测的子孔径波前数据wi可表示为:
wi=Pi+Txi·x+Tyi·y+w0i       (1)
式(1)中i=1,2…7,x、y分别是指子孔径波前数据x方向和y方向的坐标值,Pi为第i个子孔径平移调整系数,Txi、Tyi分别为x方向和y方向上倾斜误差调整系数。wi为直接检测的子孔径波前数据,w0i为第i个子孔径真实波前数据,w0i在相邻子孔径重叠区域数据应相等,因此将相邻子孔径重叠区域数据相减,即可得到子孔径间相对调整误差系数:
wi-wj=Pij+Txij·x+Tyij·y        (2)
式(2)中j=1,2…7(j≠i),wj为直接检测的第j个子孔径波前数据,子孔径i和子孔径j相对平移Pij和倾斜误差调整系数Txij、Tyij
通过式(2)可求出每个子孔径相对倾斜和平移调整误差,得到真实的子孔径波前数据,再将真实的子孔径波前数据插值到窗口玻璃12对应的区域,即可得到窗口玻璃12全口径透射波前,从而对窗口玻璃12透射波前的检测。
本发明未详细阐述的部分属于本领域的公知技术。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解想到的变换或替换,都应涵盖在本发明的包含范围之内。

Claims (5)

1.一种基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法,其特征是:利用由干涉仪、窗口玻璃、四维扫描运动平台、平面镜组成检测系统,并通过以下步骤完成对窗口玻璃透射波前的测量步骤如下:
步骤S1:根据干涉仪的口径和窗口玻璃的口径计算检测窗口玻璃所需子孔径个数,以其中任意一个子孔径中心坐标为原点,计算其余子孔径中心坐标,并对子孔径进行编号;
步骤S2:调整干涉仪与平面镜到自准直使干涉仪上条纹为零条纹;
步骤S3:将四维扫描运动平台放置在干涉仪与平面镜之间,将窗口玻璃安装在四维扫描运动平台上;
步骤S4:四维扫描运动平台带动窗口玻璃扫描运动,调整窗口玻璃与平面镜使干涉仪上的条纹为达到最少;
步骤S5:通过四维扫描运动平台的运动,将窗口玻璃调整到第一个子孔径位置,干涉仪检测窗口玻璃并记录此时子孔径透射波前数据;
步骤S6:通过四维扫描运动平台的运动,将窗口玻璃调整到第二个子孔径位置,干涉仪检测窗口玻璃并记录此时子孔径透射波前数据;
步骤S7:重复步骤S5、步骤S6直到窗口玻璃所有子孔径波前数据被检测并记录下来;
步骤S8:将所有子孔径波前数据通过算法拼接,得到窗口玻璃全口径透射波前,完成对窗口玻璃的透射波前检测。
2.根据权利要求1所述的窗口玻璃透射波前检测方法,其特征在于:步骤S1中计算子孔径个数的步骤如下:设定子孔径坐标是以任意一个包含窗口玻璃有效透射波前数据的子孔径为初始子孔径,编号为第一子孔径,以(x1,y1)作为第一子孔径中心坐标位置,其中x1,y1的坐标点分别取0对第一子孔径进行位置标定,以第一子孔径圆心(0,0)为坐标原点建立坐标系,在四维扫描运动平台沿水平方向左右运动为x方向,在四维扫描运动平台x方向向左运动为负,向右运动为正;四维扫描运动平台上下运动为y方向,在四维扫描运动平台y方向向上运动为正,向下运动为负;所有子孔径数据区域大小相同,覆盖的面积均为s;在任意方向移动第一子孔径的中心坐标得到第二子孔径,第二子孔径中心坐标为(x2,y2)与第一子孔径相邻并且有一重叠区域,重叠区域面积大小为s2,获取s2与s之比等于18%;任意方向移动第二子孔径得到第三子孔径,第三子孔径中心坐标为(x3,y3)与第二子孔径相邻且有一重叠区域,重叠区域面积大小为s3,获取s3与s之比等于18%;依次将子孔径移到其相邻位置得到下一个子孔径,记录下一个子孔径位置坐标,并保证子该孔径与其相邻子孔径即下一个子孔径有一重叠区域,且重叠区域面积与s之比等于18%,重复上述步骤直到所有子孔径能完全覆盖窗口玻璃口径。
3.根据权利要求1所述的窗口玻璃透射波前检测方法,其特征在于:所述调整干涉仪与平面镜到自准直是使干涉仪的条纹数小于5根,再在干涉仪与平面镜中间加入四维扫描运动平台,再接着调整窗口玻璃和平面镜,使干涉仪上条纹小于5根。
4.根据权利要求1所述的窗口玻璃透射波前检测方法,其特征在于:所有子孔径波前数据通过子孔径拼接算法拼接。
5.根据权利要求1所述的窗口玻璃透射波前检测方法,其特征在于:干涉仪口径小于窗口玻璃口径。
CN2013101784952A 2013-05-15 2013-05-15 基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法 Pending CN103257033A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2013101784952A CN103257033A (zh) 2013-05-15 2013-05-15 基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2013101784952A CN103257033A (zh) 2013-05-15 2013-05-15 基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法

Publications (1)

Publication Number Publication Date
CN103257033A true CN103257033A (zh) 2013-08-21

Family

ID=48961090

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2013101784952A Pending CN103257033A (zh) 2013-05-15 2013-05-15 基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法

Country Status (1)

Country Link
CN (1) CN103257033A (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103575233A (zh) * 2013-11-20 2014-02-12 西安工业大学 大口径大相对孔径抛物面反射镜面形误差的检测方法
CN103714260A (zh) * 2014-01-13 2014-04-09 成都颐泰合同能源管理有限公司 窗玻璃太阳辐射能流计算方法
CN107796597A (zh) * 2017-09-18 2018-03-13 北京空间机电研究所 一种光学系统波前的子孔径反演与实施方法
CN107990839A (zh) * 2017-11-14 2018-05-04 中国科学院上海光学精密机械研究所 长条状光学元件面形的快速拼接检测装置和测量方法
CN109855743A (zh) * 2019-01-04 2019-06-07 北方工业大学 双频激光外差干涉相位测量大尺寸光学平面的装置及方法
CN117249912A (zh) * 2023-11-20 2023-12-19 苏州致将智能光电有限公司 大口径光学元件的检测方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19800844A1 (de) * 1998-01-13 1999-07-15 Johannes Prof Dr Schwider Shack-Hartmann-Sensor mit orts-varianter Linsenanordnung
JP2011133232A (ja) * 2009-12-22 2011-07-07 Nikon Corp 光学装置
CN102507155A (zh) * 2011-11-03 2012-06-20 中国科学院光电技术研究所 一种检测大口径光学系统波前的装置
CN102620683A (zh) * 2012-03-31 2012-08-01 中国科学院长春光学精密机械与物理研究所 子孔径拼接检测非球面调整误差补偿方法
CN102980532A (zh) * 2012-12-25 2013-03-20 中国科学院长春光学精密机械与物理研究所 采用三坐标测量仪拼接测量大口径非球面面形的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19800844A1 (de) * 1998-01-13 1999-07-15 Johannes Prof Dr Schwider Shack-Hartmann-Sensor mit orts-varianter Linsenanordnung
JP2011133232A (ja) * 2009-12-22 2011-07-07 Nikon Corp 光学装置
CN102507155A (zh) * 2011-11-03 2012-06-20 中国科学院光电技术研究所 一种检测大口径光学系统波前的装置
CN102620683A (zh) * 2012-03-31 2012-08-01 中国科学院长春光学精密机械与物理研究所 子孔径拼接检测非球面调整误差补偿方法
CN102980532A (zh) * 2012-12-25 2013-03-20 中国科学院长春光学精密机械与物理研究所 采用三坐标测量仪拼接测量大口径非球面面形的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
吴晶等: "基于Shack-Hartmann的子孔径拼接波前检验技术", 《光电工程》, vol. 38, no. 1, 31 January 2011 (2011-01-31), pages 39 - 43 *
张蓉竹等: "子孔径拼接干涉检测实验研究", 《光学技术》, vol. 30, no. 2, 31 March 2004 (2004-03-31), pages 173 - 175 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103575233A (zh) * 2013-11-20 2014-02-12 西安工业大学 大口径大相对孔径抛物面反射镜面形误差的检测方法
CN103714260A (zh) * 2014-01-13 2014-04-09 成都颐泰合同能源管理有限公司 窗玻璃太阳辐射能流计算方法
CN103714260B (zh) * 2014-01-13 2017-01-04 成都颐泰合同能源管理有限公司 窗玻璃太阳辐射能流计算方法
CN107796597A (zh) * 2017-09-18 2018-03-13 北京空间机电研究所 一种光学系统波前的子孔径反演与实施方法
CN107990839A (zh) * 2017-11-14 2018-05-04 中国科学院上海光学精密机械研究所 长条状光学元件面形的快速拼接检测装置和测量方法
CN109855743A (zh) * 2019-01-04 2019-06-07 北方工业大学 双频激光外差干涉相位测量大尺寸光学平面的装置及方法
CN117249912A (zh) * 2023-11-20 2023-12-19 苏州致将智能光电有限公司 大口径光学元件的检测方法及系统
CN117249912B (zh) * 2023-11-20 2024-02-13 苏州致将智能光电有限公司 大口径光学元件的检测方法及系统

Similar Documents

Publication Publication Date Title
CN103257033A (zh) 基于子孔径拼接干涉检测的窗口玻璃透射波前检测方法
CN109242908B (zh) 用于水下双目视觉测量系统的标定方法
CN105157606B (zh) 非接触式复杂光学面形高精度三维测量方法及测量装置
Luo et al. A simple calibration procedure for structured light system
CN106248352B (zh) 一种大口径平面镜拼接检测对准方法
CN105091792A (zh) 一种标定多光轴光学系统光轴平行度的装置及其标定方法
CN204313798U (zh) 一种激光束现场标定装置
CN102654391A (zh) 基于光束平差原理的条纹投影三维测量系统及其标定方法
CN104391366A (zh) 一种太赫兹波段离轴三反射镜系统及其装调方法
CN102243068B (zh) 一种子孔径拼接中系统误差的修正方法
CN106840023A (zh) 大口径复杂曲面光学参数精确测试和标定装置及方法
CN102937421A (zh) 旋转轴对称光学非球面的实时检测方法
CN108061514B (zh) 一种利用轴向扫描光干涉法检测非球面的动态建模方法
CN110428471B (zh) 一种针对光学自由曲面子孔径偏折测量的精确自定位方法
US20230194240A1 (en) Coordinate measuring system
CN103134443B (zh) 一种大口径大径厚比反射镜面形自准直检测装置及方法
CN108286950A (zh) 一种反射镜面形的在线检测方法
CN106705888A (zh) 干涉检测中的ccd坐标系与镜面坐标系非线性关系标定方法
CN110260817A (zh) 基于虚拟标志点的复杂曲面偏折测量自定位方法
CN104089583B (zh) 一种光学系统波前的子孔径反演方法
CN102589472A (zh) 在球面面形干涉检测中高精度消除调整误差的方法
KR20110065365A (ko) 비구면체 측정 방법 및 장치
CN109827607A (zh) 线结构光焊缝跟踪传感器的标定方法及装置
CN104581150B (zh) 一种定位与误差补偿方法
JP6223368B2 (ja) 光学要素の幾何学的構造を測定する方法及びツール

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20130821