CN110188321B - 一种基于神经网络算法的主次镜校准方法 - Google Patents

一种基于神经网络算法的主次镜校准方法 Download PDF

Info

Publication number
CN110188321B
CN110188321B CN201910430789.7A CN201910430789A CN110188321B CN 110188321 B CN110188321 B CN 110188321B CN 201910430789 A CN201910430789 A CN 201910430789A CN 110188321 B CN110188321 B CN 110188321B
Authority
CN
China
Prior art keywords
neural network
coefficients
primary
axis
detuning
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
CN201910430789.7A
Other languages
English (en)
Other versions
CN110188321A (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.)
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 CN201910430789.7A priority Critical patent/CN110188321B/zh
Publication of CN110188321A publication Critical patent/CN110188321A/zh
Application granted granted Critical
Publication of CN110188321B publication Critical patent/CN110188321B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/06Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons
    • G06N3/061Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons using biological neurons, e.g. biological neurons connected to an integrated circuit
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods

Abstract

本发明公开了一种基于神经网络算法的主次镜校准方法。首先通过仿真软件建立主次镜失调模型,再利用随机加入的失调量得到相应泽尼克多项式系数;其次将失调量与相应的泽尼克多项式系数组合成数据组,重复以上步骤,得到足够的数据组;然后将其作为训练样本,对神经网络进行训练;最后,再将代表待测系统的泽尼克多项式系数输入训练好的神经网络,进而得到主次镜镜头对准的失调量。本发明方法适用于多种类型的镜头对准失调量的计算,能有效提高望远镜系统光学镜头的装调效率,同时适用于主次镜的实时校准。

Description

一种基于神经网络算法的主次镜校准方法
技术领域
本发明属于离轴反射式系统主次镜装调领域,具体涉及一种基于神经网络算法的主次镜校准方法。
背景技术
离轴反射式光学系统具有无遮拦、视场大、结构紧凑及无色差等优点,在三维测绘、空间遥感、天文观测及多光谱热成像等领域得到了广泛的应用。目前离轴反射式系统主次镜对准误差的计算速度不快,因此,如何快速又准确地得出离轴反射式系统主次镜的失调量对于实时装调具有重要意义。
早期望远镜系统主次镜的装调主要依靠人为一次次地调整,精度不够而且需要反复调整,装调效果具有明显的随机性,而主次镜的失调量完全依据装调人员经验得到,并不能定量给出。随着光学及计算机等学科的发展与进步,逐渐发展出计算机辅助装调技术,该技术利用干涉仪器测量得到系统的干涉图,经过对其分析,得到对应泽尼克多项式的系数;再通过更改失调量,得到相应的一系列系数,进而利用这种方式拟合出失调量与泽尼克多项式系数的关系。该方法在一定失调区间范围内,能够得到波长级精确的失调量。在这种装调方法中,失调量的计算主要通过求解系统灵敏度矩阵来得到。该方法的优点是当失调量小时,能够得出精确的失调量。缺点是每次求解失调量都需要大量的计算,装调效率不高,不利于实时校准,同时当主次镜具有大失调量时,失调量计算误差大,不能应用于装调;之后发展出一种逆向优化法,为解决灵敏度矩阵的弊端,它以装调自由度为变量,以实测泽尼克系数与理想泽尼克系数差值为评价函数。它的优点是:当有大失调量时,它的求解精度高。缺点是不可用于多镜系统的装调;再后来又发展出一种微分波前采样法:该方法将代表系统波像差的失调量函数二阶展开,构成Hessian矩阵,通过求解以该矩阵为系数的方程组,得到失调量。它能够很好地解决非线性与耦合问题,但精度非常依赖测量精度;之前的这三种方法都是需要建立数值化模型,后来又出现了一种矢量差模型,它以旋转对称光学系统的像差理论为依据,即:在像面上整个光学系统的像差是各表面像差贡献之和,每个面的像差场中心以连接该表面的曲率中心和光瞳中心的直线为中心。它既可用来辅助偏心倾斜光学系统的设计,又可用来分析失调光学系统,但是它的普适性较低,并且校正模型的鲁棒性有待论证。
发明内容
本发明的目的在于提供一种基于神经网络算法的主次镜校准方法,减少因为解失调量而需要的大量计算对装调效率的影响,同时提高求解失调量的准确度,获得更加稳定的装调效果。
本发明采用的技术方案是:一种基于神经网络算法的主次镜校准方法,用于望远镜系统上的主次镜对准,其具体步骤是:
步骤一模型的建立:
首先在光学仿真软件中分别输入待测望远镜系统的结构参数,建立光学模型;
步骤二训练样本的获取:
先记录无失调量的前37项泽尼克系数Z0,如下所示:
Figure BDA0002068909860000021
其中θ1,θ2,θ3表示望远镜系统三个视场角;再加入失调误差;用矩阵Δ1中的四个元素表示四个变量,其式子表示如下:
Figure BDA0002068909860000022
再得到相应的前37项泽尼克多项式的系数矩阵Zi,其表示如下:
Figure BDA0002068909860000023
其中i=1,2…,表示第i次加入失调量误差;
重复以上步骤,得到N组系数组;
步骤三训练样本的获取:
得到系数组后,先将它减去Z0,得到一个新的矩阵ZΔi,表达式如下所示:
Figure BDA0002068909860000024
其中j=1,2…37;
将得到的系列矩阵组成一个集合A,其中A的表达式为:
A=(ZΔ1…ZΔi)
同样,用集合B表示一系列的失调量,其表达式为:
B=(Δ1…Δi)
分别将这两个集合结合起来,得到两个集合的映射F,即:F:A→B,具体表达式为:F:(ZΔ1…ZΔi)→(Δ1…Δi)
即组成一个完整的训练样本集;
步骤四神经网络的训练:
选取多层神经网络;
将得到的训练集输入到多层神经网络中对其进行训练,当损失函数的值达到极小值时,整个训练过程完成;
步骤五失调量的计算:
先将自准干涉仪对待测望远镜系统进行处理,得到泽尼克多项式的前37项系数;将得到的37项系数作为输入量输入到已经经过训练的神经网络,经过神经网络的处理后,便得到了系统的失调误差。
本发明的系统组成成分主要有:光学模型,神经网络,自准干涉仪以及待测光学系统。其中神经网络主要由多层输入层,多层隐含层以及一层输出层组成。
本发明的原理是:神经网络是一种模仿生物中枢神经结构和功能的数学模型,它用于对关系函数进行近似或估计。神经网络由神经元联结进行相关数学计算,同时,它能在外界信息的基础上改变自身内部结构,从而达到最佳近似效果。因此,它也是一种自适应的非线性统计性数据建模工具。
本发明与现有技术相比,具有以下优点:
(1)经过训练后的神经网络,能够快速准确地得出望远镜系统的失调量,提高系统的装调效率;
(2)神经网络的适用性范围广,针对不同的待测系统,只需要用新样本进行训练,训练完成后即可工作,可重复使用性高。
附图说明
图1为基于神经网络算法的主次镜失调量计算方法的流程图。
图2为主次镜对准失调误差定义图。
图3为神经网络层图。
具体实施方式
下面结合附图和具体实施步骤对本发明方案作进一步描述。
如图1中右边的流程图所示,基于神经网络算法的主次镜校准方法,包括如下步骤:
步骤一模型的建立:
首先在光学仿真软件中分别输入待测离轴光学系统的结构参数,建立光学模型;
步骤二训练样本的获取:
先记录没有添加任何对准误差时的前37项泽尼克系数Z0,如下所示:
Figure BDA0002068909860000041
其中θ1,θ2,θ3表示望远镜系统不同的视场角;
如图2所示,两块镜片的失调量主要包括X和Y轴上的离轴误差dr,Z轴方向上的离焦误差dr(因为离焦误差dz容易校准,因此通常无需特意求解),以及X和Y轴上的倾斜误差ψ;
然后分别对望远镜系统的次镜加入任意的四个自由度的对准失调误差;即沿X轴方向上的离轴量,沿Y轴方向上的离轴量,沿X轴方向上的倾斜量以及沿Y轴方向上的倾斜量;
用矩阵Δ1中的四个元素表示四个变量,其式子表示如下:
Figure BDA0002068909860000042
再从光学仿真软件中的分析模块中得到相应的三个不同视场角的前37项泽尼克多项式的系数矩阵Zi,其表示如下:
Figure BDA0002068909860000043
其中i=1,2…,表示第i次加入失调量误差;
重复以上步骤,得到N组泽尼克多项式的系数组;
步骤三训练样本的获取:
得到三个不同视场角度下的泽尼克多项式系数后,先将它减去无失调量的泽尼克多项式的系数,得到一个新的矩阵ZΔi,表达式如下所示:
Figure BDA0002068909860000044
其中j=1,2…37;
分别将得到的系列矩阵组成一个集合A,其中A的表达式如下所示:
A=(ZΔ1…ZΔi)
同样的,用集合B来表示一系列的失调量,其表达式如下所示:
B=(Δ1…Δi)
分别将其结合起来,得到两个集合的映射F,即:
F:A→B
具体表达式如下所示:
F:(ZΔ1…ZΔi)→(Δ1…Δi)
即组成一个完整的训练样本集;
步骤四神经网络的训练:
因为我们选取了三个不同的视场角,同时由于误差主要分成离轴误差和倾斜误差这两类,所以我们选取多层的神经网络;
其中多层神经网络主要包括两层隐含层以及一层输出层,其结构如图3所示;
将得到的训练集输入到多层神经网络中对其进行训练,当损失函数的值达到极小值时,整个训练过程才算完成,其流程图如图1中左边所示;
步骤五失调量的计算:
先用自准干涉仪对待测望远镜系统进行处理,得到泽尼克多项式的前37项系数;
将得到的37项系数作为输入量输入到已经经过训练的神经网络,经过神经网络的处理后,便得到了系统的失调误差。
经过训练后的神经网络,能够快速准确地得出望远镜系统的失调量,提高系统的装调效率,同时可适用于主次镜实时校准阶段;神经网络训练时间长短取决于电脑配置以及训练集的大小,但当神经网络训练好后,每次对失调量的计算不超过20ms。神经网络的普适性好,针对不同的待测系统,只需要用新样本进行训练,训练完成后即可工作,可重复使用性高。

Claims (6)

1.一种基于神经网络算法的主次镜校准方法,用于望远镜系统主次镜装调过程中对光学镜头对准,其特征在于,包括步骤如下:
步骤一模型的建立:
首先在光学仿真软件中分别输入待测离轴光学系统的结构参数,建立光学模型;
步骤二训练样本的获取:
1)先记录没有添加任何对准误差时的前37项泽尼克系数
Figure FDA0003619346240000011
其中θ1,θ2,θ3分别表示望远镜系统三个不同的视场角;
2)分别对望远镜系统的次镜加入任意的四个自由度的对准失调误差,即沿X轴方向上的离轴量,沿Y轴方向上的离轴量,沿X轴方向上的倾斜量以及沿Y轴方向上的倾斜量,它们分别用以下四个变量表示:
Figure FDA0003619346240000012
再从光学仿真软件中的分析模块中得到相应的三个视场角的前37项泽尼克多项式的系数:
Figure FDA0003619346240000013
其中i=1,2…,表示第i次加入失调量误差;
重复步骤2),得到N组泽尼克多项式的系数组;
步骤三训练样本的获取:
得到了三个不同视场角度下的泽尼克多项式系数后,先将它减去无失调量的泽尼克多项式的系数:
Figure FDA0003619346240000014
其中i=1,2…;
分别将得到的系列矩阵组成一个集合A,其中A的表达式为:A=(ZΔ1…ZΔi),同样,用集合B表示一系列的失调量,其表达式为:B=(Δ1…Δi);
分别将这两个集合结合起来,得到两个集合的映射F,即:F:A→B,具体表达式为:F:(ZΔ1…ZΔi)→(Δ1…Δi);
从而得到一个完整的训练样本集;
步骤四神经网络的训练:
因为选取了三个不同视场角的泽尼克多项式系数,且误差主要由离轴误差和倾斜误差组成,所以选用多层的神经网络;
将得到的训练集输入到多层神经网络中对其进行训练,当损失函数的值达到极小值时,整个训练过程才完成;
步骤五失调量的计算:
先将自准干涉仪对待测望远镜系统进行处理,得到泽尼克多项式的前37项系数;
将得到的37项系数作为输入量输入到已经经过训练的神经网络,经过神经网络的处理后,便得到了系统的失调误差,用于主次镜的实时校准。
2.根据权利要求1所述的一种基于神经网络算法的主次镜校准方法,其特征在于,步骤二中的X,Y轴上的离轴量指的是镜片顶点偏离光轴在X,Y轴上的距离,其中坐标系是以主镜的顶点为坐标原点,水平光线传播方向为Z轴,坐标系方向符合右手定则。
3.根据权利要求1所述的一种基于神经网络算法的主次镜校准方法,其特征在于,步骤二中的泽尼克多项式的系数指泽尼克标准系数。
4.根据权利要求1所述的一种基于神经网络算法的主次镜校准方法,其特征在于,步骤三中先将它减去无失调量的泽尼克多项式的系数指有失调误差的泽尼克多项式系数减去无失调误差泽尼克多项式系数。
5.根据权利要求1所述的一种基于神经网络算法的主次镜校准方法,其特征在于,步骤四中的损失函数指的是目标函数,此处指交叉熵损失函数,它的形式是:J(W,b,a,y)=-[ylna+(1-y)ln(1-a)]。
6.根据权利要求1所述的一种基于神经网络算法的主次镜校准方法,其特征在于,步骤五中的自准干涉仪对待测望远镜系统进行处理指的是利用干涉仪对系统发出标准球面波,在经过光学系统后成为平面波;在系统前放置一块平面镜,则由系统传播过来的平面波经过反射后,又通过光学系统,进入干涉仪与标准球面波发生干涉现象,由干涉仪可以直接得到多项泽尼克多项式系数。
CN201910430789.7A 2019-05-22 2019-05-22 一种基于神经网络算法的主次镜校准方法 Active CN110188321B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910430789.7A CN110188321B (zh) 2019-05-22 2019-05-22 一种基于神经网络算法的主次镜校准方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910430789.7A CN110188321B (zh) 2019-05-22 2019-05-22 一种基于神经网络算法的主次镜校准方法

Publications (2)

Publication Number Publication Date
CN110188321A CN110188321A (zh) 2019-08-30
CN110188321B true CN110188321B (zh) 2022-07-19

Family

ID=67717283

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910430789.7A Active CN110188321B (zh) 2019-05-22 2019-05-22 一种基于神经网络算法的主次镜校准方法

Country Status (1)

Country Link
CN (1) CN110188321B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110705040A (zh) * 2019-09-11 2020-01-17 中国科学院光电技术研究所 一种基于泽尼克多项式系数与最小二乘法求取主次镜失调误差量的方法
CN111209689B (zh) * 2020-02-14 2022-07-26 北京理工大学 非零干涉非球面测量回程误差去除方法及装置
CN111399210B (zh) * 2020-04-22 2021-05-04 中国科学院长春光学精密机械与物理研究所 大口径大视场望远镜粗对准方法、装置、设备及存储介质
CN111811785B (zh) * 2020-07-21 2021-07-06 中国科学院长春光学精密机械与物理研究所 大动态范围望远镜像差的检测方法、装置、设备及介质
CN111812836A (zh) * 2020-07-21 2020-10-23 中国科学院长春光学精密机械与物理研究所 大口径大视场相机的标校装调方法、装置、设备及介质
CN113283162B (zh) * 2021-04-30 2022-10-28 晋城鸿智纳米光机电研究院有限公司 塑料光学镜片面形误差预测方法、电子装置及存储介质
CN113326653B (zh) * 2021-05-17 2022-03-29 中国工程物理研究院激光聚变研究中心 一种基于神经网络的离轴三反光学系统的设计方法
CN114859552B (zh) * 2022-06-07 2023-09-19 中国科学院光电技术研究所 一种基于深度学习的离轴望远镜低阶像差校正方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103134660A (zh) * 2013-01-30 2013-06-05 中国科学院光电技术研究所 基于像散分解获得望远镜主次镜对准误差的方法
CN103744173A (zh) * 2014-01-24 2014-04-23 中国科学院光电技术研究所 一种基于光斑清晰度函数的望远镜次镜位置校正方法
CN106646867A (zh) * 2016-12-16 2017-05-10 中国科学院光电研究院 一种深紫外光学系统共焦对准装置与方法
CN106873152A (zh) * 2017-01-11 2017-06-20 浙江大学 一种基于机器学习的高速像差校正方法
CN109596227A (zh) * 2018-12-06 2019-04-09 浙江大学 一种卷积神经网络先验增强的光学元件中频误差的相位恢复检测系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8616712B2 (en) * 2011-03-24 2013-12-31 University Of Rochester Nonsymmetric optical system and design method for nonsymmetric optical system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103134660A (zh) * 2013-01-30 2013-06-05 中国科学院光电技术研究所 基于像散分解获得望远镜主次镜对准误差的方法
CN103744173A (zh) * 2014-01-24 2014-04-23 中国科学院光电技术研究所 一种基于光斑清晰度函数的望远镜次镜位置校正方法
CN106646867A (zh) * 2016-12-16 2017-05-10 中国科学院光电研究院 一种深紫外光学系统共焦对准装置与方法
CN106873152A (zh) * 2017-01-11 2017-06-20 浙江大学 一种基于机器学习的高速像差校正方法
CN109596227A (zh) * 2018-12-06 2019-04-09 浙江大学 一种卷积神经网络先验增强的光学元件中频误差的相位恢复检测系统

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Active optical alignment of off-axis telescopes based on nodal aberration theory;Xiaobin Zhang;《Optics Express》;20150917;第24卷(第23期);26392-26413 *
基于人工神经网络方法的自由曲面光学系统装调;王钰 等;《光 学 学 报》;20131231;第33卷(第12期);81-86 *
改进的灵敏度矩阵法在离轴望远镜装调中的应用;顾志远 等;《 光学精密工程》;20150915;2595-2604 *
离轴望远镜失调误差检测技术研究;刘柱;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20210815;C030-12 *

Also Published As

Publication number Publication date
CN110188321A (zh) 2019-08-30

Similar Documents

Publication Publication Date Title
CN110188321B (zh) 一种基于神经网络算法的主次镜校准方法
CN101290218B (zh) 用于校正非球面非零位检测时的原理误差的方法
EP2369319B1 (en) Aspheric object measuring method and apparatus
CN105423948A (zh) 采用变形镜的拼接干涉检测非球面面形的装置及方法
US10323938B2 (en) Method for calibrating a measuring device
CN103134660B (zh) 基于像散分解获得望远镜主次镜对准误差的方法
CN112926192B (zh) 一种离轴自由曲面光学系统偏振像差分析方法
CN111009014A (zh) 通用成像模型正交分光成像位姿传感器的标定方法
CN115290006A (zh) 一种用于反射光路光轴对准和面形曲率检测的系统及方法
CN112394508A (zh) 一种基于二阶灵敏度矩阵法的装调方法
CN110705040A (zh) 一种基于泽尼克多项式系数与最小二乘法求取主次镜失调误差量的方法
JP4340625B2 (ja) 光学検査方法および装置
JP2007303900A (ja) 眼鏡レンズの光学性能評価方法
CN112525071B (zh) 一种抑制大口径干涉仪中光学材料非均匀性影响的方法
CN109002567B (zh) 基于实际光学系统装调检测光路的计算机辅助装调方法
KR20110065365A (ko) 비구면체 측정 방법 및 장치
CN115423881A (zh) 一种基于神经网络和光斑形态的主次镜快速校准方法
JP2000097663A (ja) 干渉計
CN111399210B (zh) 大口径大视场望远镜粗对准方法、装置、设备及存储介质
CN115127473A (zh) 一种大口径光学元件面形检测装置及检测图像畸变校正方法
CN110793465B (zh) 一种微透射元件多面形大动态范围同步测量方法
CN109579739A (zh) 一种离轴折反式部分补偿器系统及设计方法
CN110307820B (zh) 全反长焦敏捷相机二次曲面实测参数获取方法
CN114137736B (zh) 大口径光学元件装调装置及其装调方法
Wang et al. Method for testing freeform surfaces based on a Shack-Hartmann sensor with plane wavefront scanning and stitching

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