一种基于哈特曼波前传感器的模式波前复原方法
【专利摘要】本发明公开了一种基于哈特曼波前传感器的模式波前复原方法，将哈特曼波前传感器子孔径内波前的斜率和二次曲率数据同时作为有效信息，根据各阶像差模式分布构成模式斜率曲率关系矩阵，在获得实际待测波前在各子孔径中的斜率和二次曲率数据之后，即可利用最小二乘算式计算待测波前中各阶像差模式系数数据，最终利用各阶像差模式分布以及对应系数数据重构波前；本发明利用哈特曼波前传感器子孔径内本被忽略的二次曲率信息，充分挖掘有效数据，缓解波前复原能力受限于子孔径采样频率的限制，在相同子孔径分割条件下，较传统模式波前复原方法，有望实现更多的复原阶数和更高的复原精度。
【专利说明】一种基于哈特曼波前传感器的模式波前复原方法
【技术领域】
[0001] 本发明属于光学信息测量【技术领域】，涉及一种通过哈特曼波前传感器复原入射光 束波前的方法，尤其涉及一种新型的基于子孔径内二次曲率信息的哈特曼波前传感器模式 波前复原方法。
【背景技术】
[0002] 哈特曼波前传感器作为一种简洁高效的光波波前像差测量装置，已经被广泛的应 用于光学检测、自适应光学、眼科医学、激光波前诊断、激光通讯等诸多领域。目前常见的哈 特曼波前传感器主要采用R. V. Shack等人在1971年提出的改进型结构，参见"Lenticular Hartmann Screen" [Platt B.C and R.V. Shack. [J]. Opt. Sci.Newsl. 5, 15-16, 1971]，主要 采用微透镜阵列对光波波前分割取样，形成许多子孔径，并将入射的光波分别汇聚到阵列 型光电探测器（常用CCD或CMOS相机）上，阵列型光电探测器靶面上则形成与子孔径一一 对应的光斑阵列，通过记录各个子光斑质心位置的偏移变化，就可通过计算处理获得所需 的波前相位测量数据。
[0003] 哈特曼波前传感器的原理结构决定了其对波前像差的重构属于一种间接测量， 其直接测量的数据是各子光斑的质心偏移，或者说是每个子孔径内波前的倾斜斜率值。 只有结合特定的波前复原方法，才能重构整个口径的入射光波波前。目前主流的哈特曼 波前传感器波前复原方法有区域法（参见"Wave-front estimation from wave-front slope measurements，'，W.H. Southwell. [J].J0SA 70 (8) ,998-1006, 1980)和模式法（参 见"Wave-front reconstruction using a Shack-Hartmann sensor"，R. G. Lane and M. Tallon. [J].Appl. Opt. 31 (32) ,6902-6908, 1992)两种。然而，不论采用何种波前复原方 法复原波前，均将各个子孔径内的波前看作是只含倾斜像差的简单平面，复原利用的原始 信息都是各个子孔径内的两维斜率数据。这意味着哈特曼波前传感器对入射波前的采样是 与微透镜阵列密度严格一致，即每个子孔径一个波前斜率采样点（在一个坐标方向上），因 此复原方法的复原能力直接受限于子孔径密度，尤其在子孔径数目少、传感器对波前采样 稀疏的情况下，只能复原出入射波前的低阶像差或者低频信息，难以复原更多的有效信息。 但在探测弱信号目标的波前相位信息时，入射光能非常有限，为提高光电探测器上的信噪 t匕，实现哈特曼波前传感器的有效探测，采用稀疏的子孔径分割是很有必要的。因而，若能 改造目前的复原方法，突破仅对子孔径内波前斜率信息进行采样的限制条件，使之在同等 子孔径密度条件下复原更多的波前细节，或者能保持同等复原精度条件下，却以更为稀疏 的子孔径分割实现，将对进一步提升哈特曼波前传感器的探测能力，扩宽其应用领域的广 度和深度均具有重要的意义。
【发明内容】
[0004] 为了克服现有技术的不足，充分利用哈特曼波前传感器子孔径内波前的信息量， 为此，本发明的目的是提供一种基于哈特曼波前传感器的模式波前复原方法，将子孔径内 波前的二次曲率和斜率同时作为波前复原方法的有效信息，改造传统的模式波前复原方 法，充分提高哈特曼波前传感器的采样率，从而提升其波前复原能力。
[0005] 为实现所述目的，本发明提供一种基于哈特曼波前传感器的模式波前复原方法， 其特征在于通过以下步骤实现波前复原：
[0006] 步骤S1 :根据哈特曼波前传感器子孔径排布以及波前描述所采用的像差模式，计 算各阶像差模式被子孔径分割后，落在各子孔径中的平均斜率和二次曲率大小；
[0007] 步骤S2:构建哈特曼波前传感器子孔径斜率和二次曲率矩阵D;
[0008] 步骤S3:计算哈特曼波前传感器子孔径斜率和二次曲率矩阵D的逆矩阵R，即为基 于哈特曼波前传感器子孔径斜率和二次曲率信息的模式复原矩阵；
[0009] 步骤S4:在哈特曼波前传感器获得待测波前分布于各子孔径的斜率和二次曲率 的原始数据向量G之后，即可计算待测波前中的各阶像差模式系数向量
[0010] A=R?G;
[0011] 步骤S5 :根据各阶像差模式分布，利用获得的像差模式系数信息，最终复原待测 波前分布。
[0012] 本发明与现有技术相比有如下优点：突破传统哈特曼波前传感器模式复原方法对 初始数据的需求限制，将各子孔径内波前的二次曲率引入哈特曼波前传感器模式波前复原 方法，充分增加波前复原计算的有效信息；本发明能够有效缓解哈特曼波前传感器的测量 精度严重受限于子孔径排布密度的技术现状，与现有技术相比，能够在同等子孔径密度条 件下复原更多的波前细节，提升波前复原精度，或者能在保持同等复原精度条件下，以更为 稀疏的子孔径排布密度实现。本发明充分提高了对哈特曼波前传感器对子孔径内波前有效 信息的挖掘和利用，突破了传统哈特曼波前传感器波前复原方法的复原能力严重受限于子 孔径采样频率的限制，有助于提升波前复原精度或者降低对子孔径排布密度的要求，有望 应用于弱光波前探测、眼科医学等相关领域。
【专利附图】
【附图说明】
[0013] 图1为本发明的哈特曼波前传感器基本结构图；
[0014] 图2为本发明所述的基于哈特曼波前传感器的模式波前复原方法流程图；
[0015] 图3为本发明所述方法与传统方法波前复原重构结果对比图；
[0016] 图4为本发明所述方法与传统方法模式系数复原结果对比图。
【具体实施方式】
[0017] 为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照 附图，对本发明进一步详细说明。
[0018]如图1是典型的哈特曼波前传感器结构，主要由缩束系统1、微透镜阵列2和光电 探测探测器阵列3构成，其中缩束系统1由两个或者两组焦距不同的透镜同轴共焦组成，缩 束系统1是光学匹配结构，主要将入射光束的尺寸缩小到与微透镜阵列2 口径匹配，即缩束 后光束的尺寸不大于微透镜阵列2的口径，微透镜阵列2由6X6个微透镜排列构成，置于 缩束系统1之后，每个微透镜将光波分割取样并分别聚焦到位于微透镜阵列2焦面的光电 探测器阵列3上，并在其靶面上形成光斑阵列图像。传统的哈特曼波前传感器模式复原方 法要求探测每个光斑的质心偏移，得到每个子孔径内波前斜率信息，然后利用事先算好的 斜率关于模式系数的复原矩阵计算得到入射光束中各阶像差模式系数，最终用各阶像差模 式重构波前。本发明所述模式复原方法不仅需要每个子孔径内波前斜率信息，而且同时采 用每个子孔径内二次曲率信息，事先构建复原矩阵的时候同样引入子孔径二次曲率信息， 也就是说将哈特曼波前传感器子孔径内二次曲率作为有效数据，使之参与波前复原的整个 运算过程，如图2流程图所示，结合实施例的具体步骤包括：
[0019] 步骤S1 :将哈特曼波前传感器子孔径被波前分布看作具有如下形式：
[0020] Wfsnb =axx+avy+x2 +axvxy+ayly2, (1)
[0021] 其中x、y为直角坐标系下坐标，ax、ay为沿x、y坐标轴方向的平均斜率， 0^为二次曲率系数。根据上述哈特曼波前传感器子孔径排布以及前65阶Zernike像差模 式，计算各阶Zernike像差模式分布被子孔径分割后，落在各子孔径中的平均斜率和二次 曲率大小。
[0022] 步骤S2 :构建哈特曼波前传感器子孔径斜率和二次曲率矩阵D :
[0023]
【权利要求】
1. 一种基于哈特曼波前传感器的模式波前复原方法，其特征在于：通过以下步骤实现 波前复原： 步骤Sl:根据哈特曼波前传感器子孔径排布以及波前描述所采用的像差模式，计算各 阶像差模式被子孔径分割后，落在各子孔径中的平均斜率和二次曲率大小； 步骤S2 :构建哈特曼波前传感器子孔径斜率和二次曲率矩阵D; 步骤S3 :计算哈特曼波前传感器子孔径斜率和二次曲率矩阵D的逆矩阵R，即为基于哈 特曼波前传感器子孔径斜率和二次曲率信息的模式复原矩阵； 步骤S4 :在哈特曼波前传感器获得待测波前分布于各子孔径的斜率和二次曲率的原 始数据向量G之后，即可计算待测波前中的各阶像差模式系数向量： A=R·G; (1) 步骤S5 :根据各阶像差模式分布，利用获得的像差模式系数信息，最终复原待测波前 分布。
2. 根据权利要求1所述的一种基于哈特曼波前传感器的模式波前复原方法，其特征在 于：所述步骤Sl中的各子孔径中的平均斜率和二次曲率，是将子孔径内波前分布看作具有 如下形式： wLllb = axx+ayv + αχ2 χ2 +αχ1 .τν+ a ^y1, (2) 其中Wfsub表示子孔径内波前分布，x、y为直角坐标系下坐标，ax、ay为沿x、y坐标轴方 向的平均斜率，、\v、为二次曲率系数。
3. 根据权利要求1所述的一种基于哈特曼波前传感器模式波前复原方法，其特征在 于：所述步骤S2中的矩阵D具有如下形式：
式中、<1,分别表示第1阶像差模式在第k个子孔径内波 前分布中沿X方向平均斜率、沿y方向平均斜率、X2项曲率系数、xy项曲率系数和y2项曲 率系数，1取值为1至n,k取值为1至m,η为波前复原像差模式阶数,m为哈特曼波前传感 器有效子孔径数目。
4. 根据权利要求1所述的一种基于哈特曼波前传感器的模式波前复原方法，其特征在 于：所述步骤S4中的原始数据矩阵G为列向量，具有如下形式：
式中Tf、if、分别表示待测波前在第k个子孔径内波前分布沿叉方 向平均斜率、沿y方向平均斜率、X2项曲率系数、Xy项曲率系数和y2项曲率系数，k取值为 1至m，m为哈特曼波前传感器有效子孔径数目。
5. 根据权利要求1或3所述的一种基于哈特曼波前传感器的模式波前复原方法，其特 征在于：所述像差模式可以是Zernike多项式、Legendre多项式以及其他正交完备函数序 列。
【文档编号】G06F19/00GK104239740SQ201410506076
【公开日】2014年12月24日 申请日期:2014年9月26日 优先权日:2014年9月26日
【发明者】王帅, 杨平, 许冰, 何星, 刘文劲, 董理治 申请人:中国科学院光电技术研究所