CN105547179B - 一种非球面方程的测量方法 - Google Patents

一种非球面方程的测量方法 Download PDF

Info

Publication number
CN105547179B
CN105547179B CN201510865940.1A CN201510865940A CN105547179B CN 105547179 B CN105547179 B CN 105547179B CN 201510865940 A CN201510865940 A CN 201510865940A CN 105547179 B CN105547179 B CN 105547179B
Authority
CN
China
Prior art keywords
mrow
msub
aspherical
rho
measurement
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
CN201510865940.1A
Other languages
English (en)
Other versions
CN105547179A (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.)
Beijing Guowang Optical Technology Co., Ltd.
Original Assignee
Changchun Institute of Optics Fine Mechanics and Physics 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 Changchun Institute of Optics Fine Mechanics and Physics of CAS filed Critical Changchun Institute of Optics Fine Mechanics and Physics of CAS
Priority to CN201510865940.1A priority Critical patent/CN105547179B/zh
Publication of CN105547179A publication Critical patent/CN105547179A/zh
Application granted granted Critical
Publication of CN105547179B publication Critical patent/CN105547179B/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/2441Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures using interferometry

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种非球面方程的测量方法,包括:1)利用移相干涉仪测量被测非球面的环形面形图,利用测距干涉仪测量被测非球面距离上一测量点的相对位移量;2)依据所述移相干涉仪测量的环形面形图,计算获得在测量点对应的“零条纹”半径;3)依据所述“零条纹”半径和所述测距干涉仪测量的相对位移量,计算获得被测非球面在测量点对应的矢高导数值;4)依据在测量点对应的矢高导数值,计算获得被测非球面方程的导数表达式;5)依据所述被测非球面方程的导数表达式,计算获得被测非球面方程。该方法具有测量方法简单,成本低,便于操作等优点。

Description

一种非球面方程的测量方法
技术领域
本发明涉及光学检测领域,特别提供了一种非球面方程的测量方法。
背景技术
在光学设计中,单个球面可以供优化的自由度只有曲率半径;而非球面除了顶点曲率半径之外,还有二次曲面常数和高阶项系数。由于非球面比球面拥有更多设计的自由度,所以在很多光学***中,都普遍采用非球面元件来减小***的复杂度,并提高***的成像质量。
当前,对于非球面的检测,大多都集中在非球面的面形上,即采用拼接法(环形子孔径拼接或者圆形子孔径拼接)或零位补偿法(补偿镜法或计算全息法)来对非球面检测检测,很少关注对非球面方程的检测。
因此,研发一种非球面方程的检测方法,成为人们亟待解决的问题。
发明内容
鉴于此,本发明的目的在于提供一种非球面方程的测量方法,以实现在完全未知非球面方程的情况下测量获得非球面方程。
本发明提供的技术方案具体为,一种非球面方程的测量方法,其特征在于,包括:
1)利用移相干涉仪测量被测非球面的环形面形图,利用测距干涉仪测量被测非球面距离上一环形面形图测量点的相对位移量;
2)依据所述移相干涉仪测量的环形面形图,计算获得在测量点对应的“零条纹”半径;
3)依据所述“零条纹”半径、所述测距干涉仪测量的相对位移量以及非球面方程的导数表达式,采用迭代优化法计算出被测非球面的方程。
优选,所述依据所述移相干涉仪测量的环形面形图,计算获得在测量点对应的“零条纹”半径,包括:
如果环形干涉条纹与图像中心点不连通,在环形面形图上采用数值微分法计算出满足阈值要求的坐标,而后采用最小二乘法计算获得在测量点对应的“零条纹”半径;
如果环形干涉条纹与图像中心点连通,采用Zernike的离焦和球差拟合法,计算获得在测量点对应的“零条纹”半径。
进一步优选,所述采用Zernike的离焦和球差拟合法,计算获得在测量点对应的“零条纹”半径公式为:
其中,a1表示离焦系数,a2表示球差系数,且所述a1和a2可通过对环形面形图进行Zernike拟合获得。
进一步优选,所述非球面方程的导数表达式为:
其中,c表示最佳拟合球曲率,ρ表示“零条纹”半径,ρmax表示非球面径向坐标最大值,M表示整数,bn表示多项式{Pn(x)}的系数,且多项式{Pn(x)}满足:
进一步优选,M=6。
进一步优选,所述迭代优化法具体为:
1)利用测距干涉仪测量“猫眼”位置与所述环形面形图中间环带位置的距离,并将测量得到的距离作为所述被测非球面方程中最佳拟合球曲率c的初值,并依据最佳拟合球曲率c的初值,估算出对应的矢高值z(ρi);
2)依据最佳拟合球曲率c和矢高值z(ρi),利用公式(a)和最小二乘法,计算获得多项式{Pn(x)}的系数bn
其中,ρi表示在测量点对应的“零条纹”半径,si表示测距干涉仪测量的相对位移量,ρ表示“零条纹”半径;
3)依据系数bn获得被测非球面方程的导数具体表达式,并依据所述具体表达式计算获得非球面方程;
4)依据所述非球面方程,重新计算获得最佳拟合球曲率c以及矢高值z(ρi),将计算获得的最佳拟合球曲率c与步骤2)代入的最佳拟合球曲率c相减,获得差值;
5)将所述差值与阈值比较,当大于阈值时,重复步骤2)~步骤4),直至所述差值小于等于阈值。
本发明提供的非球面方程的测量方法,可实现在完全未知非球面方程的情况下测量获得非球面方程,丰富了光学元件的可测量参数,具有测量方法简单,成本低,便于操作等优点。
附图说明
图1为非球面方程的测量***示意图;
图2为仿真生成的环形干涉图;
图3为计算获得“零条纹”位置示意图;
图4为当干涉条纹与图像中心不连通时,满足设定阈值的“零条纹”位置分布图;
图5为利用图2中“零条纹”位置计算出的最佳拟合圆;
图6为当干涉条纹与图像中心连通时,仿真生成的干涉图;
图7为图6对应的面形图;
图8为某一测量位置实际测量的环形干涉图;
图9为图8对应的面形图;
图10为测长干涉仪实际测量的相对位移与计算出的“零条纹”的位置关系;
图11为最终获得的非球面方程计算出的矢高;
图12为轮廓测量的矢高与发明中论述方法测量的矢高的偏差。
具体实施方式
下面以具体的实施方案对本发明进行进一步解释,但是并不用于限制本发明的保护范围。
为了实现在完全未知非球面方程的情况下测量获得非球面方程,本实施方案提供了一种非球面方程的测量方法,其中,所用的测量装置参见图1,包括:稳频激光器S1、准直扩束***S2、分光棱镜S3、汇聚镜组S4、移相参考面S5、测距干涉仪S6、非球面S7、夹持机构S8、多维调整台S9、成像镜S10、探测器S11、图像采集单元S12,其中,S1、S2、S3、S4、S5、S10、S11和S12构成了移相干涉仪。
具体的测量过程如下:
准备工作:选择恰当的球面参考镜,将被测非球面与移相干涉仪放置同心,将测距干涉仪置零;
1)沿移相干涉仪的轴向移动被测非球面,利用移相干涉仪测量被测非球面的环形面形图,利用测距干涉仪测量被测非球面轴向移动距离。
2)依据移相干涉仪测量的环形面形图,计算获得在测量点对应的“零条纹”半径;
3)依据“零条纹”半径、测距干涉仪测量的轴向移动距离以及非球面方程的导数表达式,采用迭代优化的计算方法获得被测非球面方程。
在本实施方中,步骤2)的计算方法如下:
当移相干涉仪发出的球面波与被测非球面径向任一位置法线垂直时,会出现所示的环形干涉图如图2所示,计算示意图如图3所示。P为非球面面上的任意一点,非球面上过P点的法线与光轴相交于点F;若P点处于“零条纹”位置,则P点与F点的坐标可以表示为(ρi,z(ρi))和(0,si),且满足:
其中ρi表示所采集的第i环带所对应的非球面径向位置,z(ρi)表示此位置处非球面的矢高,z'(ρi)表示此位置处矢高的导数,si表示测距干涉仪测量的相对位移。
为了便于计算,公式(1)可以写成:
并且满足:
在面形图上,“零条纹”位置满足:
其中F(x,y)表示测量的面形。
具体计算ρi(“零条纹”半径)可以分为两种情况:
情况一:当环形面形图中干涉条纹与图像中心不连通时,为了便于计算,公式(4)可以写成:
对面形图采用数值微分法,设置合适的阈值,可以得到“零条纹”位置分布图,如图4所示,采用最小二乘法进行拟合,即可得到“零条纹”的中心与半径,如图5所示。
情况二:当环形面形图中干涉条纹与图像中心连通时(此时移相干涉仪焦点接近非球面顶点球球心),所测量的面形主要由离焦和球差项组成,如图6、图7所示,即:
F(x,y)≈a1(2ρ2-1)+a2(6ρ4-6ρ+1) (6)
其中a1和a2分别表示离焦和球差系数,其可以通过对环形面形图进行Zernike拟合获得。
将公式(6)带入公式(4),可得:
在本实施方案中,步骤3)的计算优化方法为:
为了提高测量精度,被测非球面方程的导数表达式描述为:
其中,c为最佳拟合球曲率,ρmax为非球面径坐标最大值,M可以根据实际的需要进行整数取值,本实施方案中M取6,bn表示多项式{Pn(x)}的系数,且多项式{Pn(x)}满足:
根据公式(8),非球面方程可以表示为
其中Tn(x)=∫Pn(x)dx,n≥0。
令公式(3)与公式(10)相等,为了便于迭代优化,可以写为:
具体的迭代优化过程如下:
(Ⅰ)利用测距干涉仪测量“猫眼”位置与所述环形面形图中间环带位置的距离,并将测量得到的距离作为所述被测非球面方程中最佳拟合球曲率c的初值,并依据最佳拟合球曲率c的初值,估算出被测非球面的矢高z(ρi);
(Ⅱ)将测距干涉仪测量的相对位移、步骤2)中获得的“零条纹”半径以及最佳拟合球曲率c和矢高z(ρi),带入公式(11),采用最小二乘法计算出拟合系数{bn};
(Ⅲ)利用公式(10)计算出非球面的方程;
(Ⅳ)根据非球面方程重新计算最佳拟合球的曲率c和矢高z(ρi),将计算获得的最佳拟合球曲率与步骤(Ⅱ)代入的最佳拟合球曲率相减,获得差值;
(Ⅴ)将所述差值与阈值比较,当大于阈值时,重复步骤(Ⅱ)~步骤(Ⅳ),直至所述差值小于等于阈值,此时,获得的非球面方程即为被测非球面方。
通过上述方法针对一高次偶次非球面进行测量,此非球面参数如表1所示:
表1非球面参数
利用步骤1)~步骤3)完成了对非球面方程的测量,结果如图8~11所示,图8为某一测量位置的环形干涉图,图9为与干涉图相对应的相位图,图10为测长干涉仪测量的相对位移与“零条纹”的位置,图11利用测量获得的非球面方程计算出的非球面矢高。
为了验证此方法的可靠性,此非球面也采用轮廓仪进行了矢高的测量,二者测量的偏差如图12所示,从图中可以看出,二者测量的PV值偏差约为1.6μm。

Claims (4)

1.一种非球面方程的测量方法,其特征在于,包括:
1)利用移相干涉仪测量被测非球面的环形面形图,利用测距干涉仪测量被测非球面距离上一环形面形图测量点的相对位移量;
2)依据所述移相干涉仪测量的环形面形图,计算获得在测量点对应的“零条纹”半径;
3)依据所述“零条纹”半径、所述测距干涉仪测量的相对位移量以及非球面方程的导数表达式,采用迭代优化法计算出被测非球面的方程;
所述依据所述移相干涉仪测量的环形面形图,计算获得在测量点对应的“零条纹”半径,包括:
如果环形干涉条纹与图像中心点不连通,在环形面形图上采用数值微分法计算出满足阈值要求的坐标,而后采用最小二乘法计算获得在测量点对应的“零条纹”半径;
如果环形干涉条纹与图像中心点连通,采用Zernike的离焦和球差拟合法,计算获得在测量点对应的“零条纹”半径;
所述采用Zernike的离焦和球差拟合法,计算获得在测量点对应的“零条纹”半径公式为:
<mrow> <msub> <mi>&amp;rho;</mi> <mi>i</mi> </msub> <mo>=</mo> <msqrt> <mfrac> <mrow> <mn>3</mn> <msub> <mi>a</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>a</mi> <mn>1</mn> </msub> </mrow> <mrow> <mn>6</mn> <msub> <mi>a</mi> <mn>2</mn> </msub> </mrow> </mfrac> </msqrt> </mrow>
其中,a1表示离焦系数,a2表示球差系数,且所述a1和a2可通过对环形面形图进行Zernike拟合获得。
2.按照权利要求1所述非球面方程的测量方法,其特征在于,所述非球面方程的导数表达式为:
<mrow> <msup> <mi>z</mi> <mo>&amp;prime;</mo> </msup> <mrow> <mo>(</mo> <mi>&amp;rho;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>c</mi> <mi>&amp;rho;</mi> </mrow> <msqrt> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mi>c</mi> <mn>2</mn> </msup> <msup> <mi>&amp;rho;</mi> <mn>2</mn> </msup> </mrow> </msqrt> </mfrac> <mo>+</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>M</mi> </munderover> <msub> <mi>b</mi> <mi>n</mi> </msub> <msub> <mi>P</mi> <mi>n</mi> </msub> <mrow> <mo>(</mo> <mi>&amp;rho;</mi> <mo>/</mo> <msub> <mi>&amp;rho;</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow>
其中,c表示最佳拟合球曲率,ρ表示“零条纹”半径,ρmax表示非球面径向坐标最大值,M表示整数,bn表示多项式{Pn(x)}的系数,且多项式{Pn(x)}满足:
3.按照权利要求2所述非球面方程的测量方法,其特征在于:M=6。
4.按照权利要求2所述非球面方程的测量方法,其特征在于,所述迭代优化法具体为:
1)利用测距干涉仪测量“猫眼”位置与所述环形面形图中间环带位置的距离,并将测量得到的距离作为所述被测非球面方程中最佳拟合球曲率c的初值,并依据最佳拟合球曲率c的初值,估算出对应的矢高值z(ρi);
2)依据最佳拟合球曲率c和矢高值z(ρi),利用公式(a)和最小二乘法,计算获得多项式{Pn(x)}的系数bn
<mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mi>M</mi> </munderover> <msub> <mi>b</mi> <mi>n</mi> </msub> <msub> <mi>P</mi> <mi>n</mi> </msub> <mrow> <mo>(</mo> <mi>&amp;rho;</mi> <mo>/</mo> <msub> <mi>&amp;rho;</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <msub> <mi>&amp;rho;</mi> <mi>i</mi> </msub> <mrow> <msub> <mi>s</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>z</mi> <mrow> <mo>(</mo> <msub> <mi>&amp;rho;</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>-</mo> <mfrac> <mrow> <mi>c</mi> <mi>&amp;rho;</mi> </mrow> <msqrt> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mi>c</mi> <mn>2</mn> </msup> <msup> <mi>&amp;rho;</mi> <mn>2</mn> </msup> </mrow> </msqrt> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mi>a</mi> <mo>)</mo> </mrow> </mrow>
其中,ρi表示在测量点对应的“零条纹”半径,si表示测距干涉仪测量的相对位移量,ρ表示“零条纹”半径;
3)依据系数bn获得被测非球面方程的导数具体表达式,并依据所述具体表达式计算获得非球面方程;
4)依据所述非球面方程,重新计算获得最佳拟合球曲率c以及矢高值z(ρi),将计算获得的最佳拟合球曲率c与步骤2)代入的最佳拟合球曲率c相减,获得差值;
5)将所述差值与阈值比较,当大于阈值时,重复步骤2)~步骤4),直至所述差值小于等于阈值。
CN201510865940.1A 2015-12-01 2015-12-01 一种非球面方程的测量方法 Active CN105547179B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510865940.1A CN105547179B (zh) 2015-12-01 2015-12-01 一种非球面方程的测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510865940.1A CN105547179B (zh) 2015-12-01 2015-12-01 一种非球面方程的测量方法

Publications (2)

Publication Number Publication Date
CN105547179A CN105547179A (zh) 2016-05-04
CN105547179B true CN105547179B (zh) 2018-01-26

Family

ID=55826539

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510865940.1A Active CN105547179B (zh) 2015-12-01 2015-12-01 一种非球面方程的测量方法

Country Status (1)

Country Link
CN (1) CN105547179B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105910552B (zh) * 2016-05-16 2018-07-03 中国科学院西安光学精密机械研究所 一种光学***立式装检装置及装检方法
CN106404354A (zh) * 2016-10-11 2017-02-15 中国科学院长春光学精密机械与物理研究所 非球面补偿器透射波前方程的测量装置和方法
WO2018068199A1 (zh) * 2016-10-11 2018-04-19 中国科学院长春光学精密机械与物理研究所 非球面补偿器透射波前方程的测量装置和方法
CN108955532B (zh) * 2018-08-23 2020-05-05 中国科学院上海光学精密机械研究所 用于绝对检测的大口径光学平面镜的旋转装置
CN109737892B (zh) * 2019-01-08 2020-06-23 北京卫星环境工程研究所 基于区域定位拟合算法的数字莫尔移相干涉面形测量方法
CN113011439A (zh) * 2021-03-19 2021-06-22 中国科学院长春光学精密机械与物理研究所 一种干涉条纹位移信息提取方法
CN112923871B (zh) * 2021-03-31 2021-12-28 中国科学院长春光学精密机械与物理研究所 一种自由曲面反射镜曲率半径检测装置及方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101949691A (zh) * 2010-09-07 2011-01-19 中国科学院长春光学精密机械与物理研究所 非零位补偿浅度光学非球面面形检测方法
CN103335613A (zh) * 2013-07-18 2013-10-02 中国科学院光电技术研究所 一种大口径非球面主镜检测装置与方法
CN104034262A (zh) * 2014-06-25 2014-09-10 中国科学院光电技术研究所 一种非球面镜二次常数测量方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101949691A (zh) * 2010-09-07 2011-01-19 中国科学院长春光学精密机械与物理研究所 非零位补偿浅度光学非球面面形检测方法
CN103335613A (zh) * 2013-07-18 2013-10-02 中国科学院光电技术研究所 一种大口径非球面主镜检测装置与方法
CN104034262A (zh) * 2014-06-25 2014-09-10 中国科学院光电技术研究所 一种非球面镜二次常数测量方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于曲面拟合的离轴非球面镜顶点半径计算方法;郭玲玲等;《红外与激光工程》;20140831;第43卷(第8期);第2694-2697页 *
斐索式红外干涉仪研制和测试技术研究;袁群等;《中国激光》;20110831;第38卷(第8期);第0808001-1-0808001-6页 *
用计算全息图对非球面的高精度检测与误差评估;高松涛等;《光学学报》;20130630;第33卷(第6期);第0612003-1-0612003-7页 *

Also Published As

Publication number Publication date
CN105547179A (zh) 2016-05-04

Similar Documents

Publication Publication Date Title
CN105547179B (zh) 一种非球面方程的测量方法
CN103591888B (zh) 大口径离轴非球面光学元件几何参数的测算方法
CN102620683B (zh) 子孔径拼接检测非球面调整误差补偿方法
CN103926058B (zh) 非球面检测中采用自准平面镜测量光轴的方法
CN104142129B (zh) 离轴三反非球面***凸非球面次镜面形拼接检测方法
CN103134660B (zh) 基于像散分解获得望远镜主次镜对准误差的方法
CN106404354A (zh) 非球面补偿器透射波前方程的测量装置和方法
CN103471522B (zh) 检测范围广的凹非球面的实时检测方法
CN103234480A (zh) 一种环形凸非球面的快速面形检测方法
CN106468544B (zh) 基于光电自准直仪的卫星高精度测角方法
CN103969787A (zh) 一种离轴四反镜头的初装定位方法
CN102937421B (zh) 旋转轴对称光学非球面的实时检测方法
CN103557791B (zh) 一种大口径非球面主镜二次常数测量装置与方法
Li et al. Improved zonal integration method for high accurate surface reconstruction in quantitative deflectometry
CN106595529B (zh) 基于虚拟牛顿环的大曲率半径非零位干涉测量方法及装置
CN103292738A (zh) 一种球面面形误差绝对检测方法
CN105318847A (zh) 基于***建模的非球面非零位环形子孔径拼接方法
JP2013186024A (ja) 非球面形状計測方法、非球面形状計測装置、光学素子加工装置および光学素子
CN105444693A (zh) 一种浅度非球面的面形误差测量方法
CN101819017B (zh) 大口径非球面反射镜顶点曲率半径的检测装置和方法
CN105157572A (zh) 一种用于非球面环形子孔径拼接的中心偏移误差补偿方法
CN105352451B (zh) 一种基于可变形镜的准万能补偿镜及设计方法
Lyu et al. Robust and accurate measurement of optical freeform surfaces with wavefront deformation correction
CN104034262B (zh) 一种非球面镜二次常数测量方法
CN103196391A (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
TR01 Transfer of patent right

Effective date of registration: 20190905

Address after: Room 601-10, 6th floor, No. 2, Jingyuan Beijie, Beijing Economic and Technological Development Zone, Daxing District, Beijing, 100176

Patentee after: Beijing Guowang Optical Technology Co., Ltd.

Address before: No. 3888, Southeast Lake Road, Changchun Economic and Technological Development Zone, 130000 Jilin Province

Patentee before: Changchun Inst. of Optics and Fine Mechanics and Physics, Chinese Academy of Sci

TR01 Transfer of patent right