CN111207910B - 基于色散条纹斜率分析的拼接镜共位相误差校正方法 - Google Patents
基于色散条纹斜率分析的拼接镜共位相误差校正方法 Download PDFInfo
- Publication number
- CN111207910B CN111207910B CN202010047451.6A CN202010047451A CN111207910B CN 111207910 B CN111207910 B CN 111207910B CN 202010047451 A CN202010047451 A CN 202010047451A CN 111207910 B CN111207910 B CN 111207910B
- Authority
- CN
- China
- Prior art keywords
- dispersion
- slope
- fringe
- mirrors
- sub
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
- G01M11/0242—Testing optical properties by measuring geometrical properties or aberrations
- G01M11/0271—Testing optical properties by measuring geometrical properties or aberrations by using interferometric methods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
- G01M11/0242—Testing optical properties by measuring geometrical properties or aberrations
- G01M11/0257—Testing optical properties by measuring geometrical properties or aberrations by analyzing the image formed by the object to be tested
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
基于色散条纹斜率分析的拼接镜共位相误差校正方法,属于光学检测领域,包括:构建检测平台;采集色散条纹图像,采用条纹斜率粗提取方法计算条纹斜率;若条纹斜率绝对值大于设定的阈值T1,则进行粗调整并计算两个拼接子镜之间的平移误差,重复步骤二;若条纹斜率绝对值小于设定的阈值T1,则进行精调整并重新计算条纹斜率;若条纹斜率绝对值小于设定的阈值T2,则迭代调整结束,否则重复步骤四。本发明利用子镜平移误差与色散条纹斜率之间存在定量化解析关系,通过分析色散条纹斜率确定子镜间平移误差的方向与大小,并进行迭代调整,直至满足迭代终止条件时结束。通过本发明的多次迭代调整之后可以使子镜间残留平移误差足够小。
Description
技术领域
本发明属于光学检测技术领域,具体涉及一种基于色散条纹斜率分析的拼接镜共位相误差校正方法。
背景技术
空间大口径望远镜对于我国科学进步与社会发展具有重要意义,其不仅可作为研究宇宙起源、天体形成与演化、暗物质与暗能量本质、地外生命***等天文学热点问题的关键设备,还能够用于高分辨率对地遥感,在灾害监测、资源勘探、地图测绘以及军事等方面发挥重要作用。
传统采用单块镜面的空间大口径望远镜在镜面加工、检测以及整个镜体的运输、发射等方面存在很大困难,采用拼接式主镜结构已成为新一代空间大口径望远镜发展的重要趋势。然而,拼接式主镜各子镜之间的共位相误差需要严格校正,若要实现衍射极限成像分辨率,拼接式主镜各子镜之间的共位相误差需要小于1/40波长。
色散条纹传感技术为共位相误差检测的常用方法之一。其利用宽光谱光源与色散元件实现宽谱段干涉,提高了对于子镜间平移误差的检测范围。与哈特曼相位传感器、四棱锥传感器以及泽尼克相衬传感器等相比,其不仅检测范围不受2π不确定性影响,而且不需要严格的光瞳对准或焦点对准,有较好的应用前景。色散条纹技术需要特定的色散条纹处理方法,从色散条纹中提取共位相误差信息。然而,当前的色散条纹处理方法普遍依赖于标定,包括波长标定、光源光谱标定、探测器位姿标定等。对于拼接式空间望远镜,这些复杂的标定过程将大大提高***造价与***复杂度,降低可靠性。
发明内容
本发明的目的是提供一种基于色散条纹斜率分析的拼接镜共位相误差校正方法,通过充分利用二维色散条纹斜率信息并从色散条纹自身寻找基准,省去了复杂的标定环节,更加适合在轨应用。
本发明为解决技术问题所采用的技术方案如下:
本发明的基于色散条纹斜率分析的拼接镜共位相误差校正方法,包括以下步骤:
步骤一、构建检测平台;
步骤二、采集色散条纹图像,采用条纹斜率粗提取方法计算条纹斜率;
步骤三、若条纹斜率绝对值大于设定的阈值T1,则进行粗调整并计算两个拼接子镜之间的平移误差,重复步骤二;
步骤四、若条纹斜率绝对值小于设定的阈值T1,则进行精调整并重新计算条纹斜率;
步骤五、若条纹斜率绝对值小于设定的阈值T2,则迭代调整结束,否则重复步骤四。
进一步的,步骤三中,采用图像二值化、连通域分析和方位角计算方法进行粗调整。
进一步的,步骤四中,采用低通滤波和数据拟合方法进行精调整。
进一步的,步骤一中,所述检测平台包括:宽光谱光源、准直镜、分光棱镜、透镜、两个拼接子镜、光阑、波前传感模块、两个压电陶瓷促动器、两个PI六足调整机构;两个拼接子镜分别连接两个压电陶瓷促动器,两个压电陶瓷促动器分别安装在两个PI六足调整机构上;
宽光谱光源发出的光经准直镜准直后变为平行光,平行光依次经分光棱镜、透镜作用后变为会聚光,会聚光经两个拼接子镜反射后变为反射光,反射光经透镜作用后再次变为平行光,平行光经分光棱镜反射后形成相位不同的两部分波前,这两部分波前分别经光阑的两个孔后再经波前传感模块进行成像,形成色散条纹图像。
进一步的,所述波前传感模块由色散棱栅、成像透镜和CCD探测器组成,两部分波前分别经光阑的两个孔后再经色散棱栅色散后形成色散条纹,经成像透镜成像在CCD探测器像面,形成色散条纹图像。
进一步的,步骤二具体包括以下步骤:
采集的色散条纹图像上不同位置(x,y)的光强分布公式如下:
其中,a、b和d分别表示光阑的宽度、长度以及两孔中心间距,x表示色散条纹图像上不同位置的横坐标,y表示色散条纹图像上不同位置的纵坐标,λ表示与所采集色散条纹图像相对应的波长,(λmin,λmax)表示与所采集色散条纹图像相对应的波长范围,f表示成像透镜等效焦距,φ(λ)表示色散棱栅的色散角,δ表示两个拼接子镜反射波面的光程差,即两个拼接子镜之间的平移误差;
二维色散条纹图像取得暗条纹的条件为:
λ与x之间的关系为:
λ(x)=λ0+C0x (3)
其中,λ0为宽光谱中心波长,C0为色散棱栅的线色散倒数;
进一步的,步骤三具体包括以下步骤:
若条纹斜率绝对值大于设定的阈值T1,则进行粗调整并计算两个拼接子镜之间的平移误差:在色散条纹中心线附近即y≈0处,条纹斜率的大小与两个拼接子镜反射波面的光程差δ正相关,且条纹斜率符号与光程差符号一一对应,根据条纹斜率的方向与大小,计算两个拼接子镜反射波面的光程差δ,然后重复步骤二。
进一步的,步骤四具体包括以下步骤:
若条纹斜率绝对值小于设定的阈值T1,则截取一系列不同波长位置的色散条纹数据,对每一列色散条纹数据进行低通滤波,滤掉毛刺,得到每一列色散条纹数据的峰值位置,同时取能量最强的三个峰值进行抛物线拟合,记录所拟合的抛物线顶点的位置,最终得到一系列峰值点与抛物线顶点;将所有的峰值点拟合出一条直线,将所有的抛物线顶点拟合出一条直线,峰值线相对于抛物线顶点拟合线的斜率即为色散条纹斜率的精确计算值。
本发明的有益效果是:
本发明在迭代调整过程中,采用粗、精两级迭代调整方式。在粗调整阶段,依靠图像处理手段对条纹斜率进行快速近似估算,主要利用图像二值化、连通域分析、方位角计算等方法;在精调整阶段,主要通过数据拟合的手段准确计算色散条纹相对于色散方向的倾斜,即该斜率计算过程需要从条纹自身寻找基准(色散方向),然后准确计算条纹相对于色散方向的倾斜,在此过程中需要利用低通滤波、数据拟合等方法。
本发明依据的基本原理是子镜平移误差与色散条纹斜率之间存在定量化解析关系,通过分析色散条纹斜率确定子镜间平移误差的方向与大小,并进行迭代调整,直至满足迭代终止条件(子镜间平移误差足够小)时结束。通过本发明的多次迭代调整之后可以使子镜间残留平移误差足够小,本发明校正精度高,适用于拼接式空间大口径天文望远镜共位相误差的检测。
附图说明
图1为检测平台结构组成示意图。
图2为原始色散条纹、对比度强化之后的色散条纹以及二值化色散条纹图像示意图。
图3为不同平移误差所对应的色散条纹图像及相应的二值化图像。
图4为经过粗调整后的色散条纹图像。
图5为得到的一系列不同波长位置的色散条纹数据。
图6为低通滤波结果。
图7为数据拟合结果。
图中,1、宽光谱光源,2、准直镜,3、分光棱镜,4、光阑,5、色散棱栅,6、成像透镜,7、CCD探测器,8、拼接子镜,9、透镜。
具体实施方式
本发明的基于色散条纹斜率分析的拼接镜共位相误差校正方法,主要包括以下步骤:
步骤一、根据如图1所示的色散条纹实现原理构建检测平台。检测平台主要包括:宽光谱光源1、准直镜2、分光棱镜3、透镜9、两个拼接子镜8、光阑4、波前传感模块、两个压电陶瓷促动器、两个PI六足调整机构;两个拼接子镜分别连接两个压电陶瓷促动器,两个压电陶瓷促动器分别安装在两个PI六足调整机构上。
宽光谱光源1发出的光经准直镜2准直后变为平行光,平行光依次经分光棱镜3、透镜9作用后变为会聚光,会聚光经两个拼接子镜8(两个拼接子镜8之间存在平移误差)反射后变为反射光,反射光经透镜9作用后再次变为平行光,平行光经分光棱镜3反射后形成相位不同的两部分波前,这两部分波前分别经光阑4的两个孔后再经波前传感模块进行成像,形成色散条纹图像。色散条纹本质上就是相同光程差条件下不同波长的干涉图像按照色散方向排布的结果。
其中,波前传感模块由色散棱栅5、成像透镜6和CCD探测器7组成,两部分波前分别经光阑4的两个孔后再经色散棱栅5色散后形成色散条纹,经成像透镜6成像在CCD探测器7像面,形成色散条纹图像。
步骤二、采集色散条纹图像,采用条纹斜率粗提取方法计算条纹斜率,以双矩形光阑所形成的色散条纹为例,则条纹斜率的计算过程具体包括以下步骤:
采集的色散条纹图像上不同位置(x,y)的光强分布公式如下:
其中,a、b和d分别表示光阑的宽度、长度以及两孔中心间距,x表示色散条纹图像上不同位置的横坐标,y表示色散条纹图像上不同位置的纵坐标,λ表示与所采集色散条纹图像相对应的波长,(λmin,λmax)表示与所采集色散条纹图像相对应的波长范围,f表示成像透镜等效焦距,φ(λ)表示色散棱栅的色散角,δ表示两个拼接子镜反射波面的光程差,即两个拼接子镜之间的平移误差;
二维色散条纹图像取得暗条纹的条件为:
λ与x之间的关系为:
λ(x)=λ0+C0x (3)
其中,λ0为宽光谱中心波长,C0为色散棱栅的线色散倒数;
步骤三、若条纹斜率绝对值大于设定的阈值T1,则采用图像二值化、连通域分析和方位角计算方法,同时利用PI六足调整机构进行粗调整并计算两个拼接子镜之间的平移误差:在色散条纹中心线附近即y≈0处,条纹斜率的大小与两个拼接子镜反射波面的光程差δ正相关,且条纹斜率符号与光程差符号一一对应,根据条纹斜率的方向与大小,计算两个拼接子镜反射波面的光程差δ,然后重复步骤二。
当两个拼接子镜间平移误差量级较大时,本发明通过亮条纹斜率近似计算暗条纹斜率,采用图像处理手段,首先对色散条纹进行对比度强化,然后进行图像二值化,然后以二值化色散条纹图像中面积最大的连通域的斜率作为条纹斜率的估计值。原始色散条纹、对比度强化之后的色散条纹以及二值化色散条纹图像如图2所示。与不同平移误差相对应的色散条纹图像二值化中间结果如图3所示,图3a为平移误差≈50μm,图3b为平移误差≈30μm,图3c为平移误差≈20μm,图3d为平移误差≈10μm。
利用以上粗调整方法可有效降低两个拼接子镜之间的平移误差,色散条纹也相应被调至近似水平。当计算的色散条纹斜率小于一定阈值时,将进入精调整步骤,在此步骤中需要进一步对色散条纹斜率进行精确计算。
步骤四、若条纹斜率绝对值小于设定的阈值T1,则采用低通滤波和数据拟合方法,同时利用压电陶瓷促动器进行精调整并重新计算条纹斜率:经过粗调整,色散条纹已经接近水平,如图4所示,此时,则截取一系列不同波长位置的色散条纹数据(为一维数据,对应单波长双孔干涉强度分布),如图5所示,对每一列色散条纹数据进行低通滤波,滤掉毛刺。在此基础上,一方面,得到每一列色散条纹数据的最大值(峰值)位置;另一方面,同时取能量最强的三个峰值进行抛物线拟合,记录所拟合的抛物线顶点的位置。最终,每一列色散条纹数据都能得到上述两个位置,于是便可以得到一系列峰值点E与抛物线顶点F,如图6所示;然后将所有的峰值点拟合出一条直线,将所有的抛物线顶点拟合出一条直线,如图7所示,峰值线相对于抛物线顶点拟合线的斜率即为色散条纹斜率的精确计算值。
在精调整过程中,可以进行步进调整,例如每次根据计算的条纹斜率符号使压电陶瓷促动器向相应方向步进50nm。
步骤五、在精调整迭代调整过程中,若条纹斜率绝对值小于设定的阈值T2,则迭代调整结束,否则重复步骤四。
本发明利用双频激光干涉仪搭建了干涉检验光路,其发出的球面波经过两个拼接子镜反射之后照射到球面镜上,光线被反射回干涉仪,携带共位相误差的实际波前与参考波前发生干涉,检测两个拼接子镜之间的共位相误差。检验结果能够证明,在不需要任何标定的情况下,利用本发明的方法能够将子镜间平移误差校正至100nm以内。本发明能够省去复杂标定过程,特别适合在轨应用。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (5)
1.基于色散条纹斜率分析的拼接镜共位相误差校正方法,其特征在于,包括以下步骤:
步骤一、构建检测平台;
所述检测平台包括:宽光谱光源、准直镜、分光棱镜、透镜、两个拼接子镜、光阑、波前传感模块、两个压电陶瓷促动器、两个PI六足调整机构;两个拼接子镜分别连接两个压电陶瓷促动器,两个压电陶瓷促动器分别安装在两个PI六足调整机构上;
宽光谱光源发出的光经准直镜准直后变为平行光,平行光依次经分光棱镜、透镜作用后变为会聚光,会聚光经两个拼接子镜反射后变为反射光,反射光经透镜作用后再次变为平行光,平行光经分光棱镜反射后形成相位不同的两部分波前,这两部分波前分别经光阑的两个孔后再经波前传感模块进行成像,形成色散条纹图像;
所述波前传感模块由色散棱栅、成像透镜和CCD探测器组成,两部分波前分别经光阑的两个孔后再经色散棱栅色散后形成色散条纹,经成像透镜成像在CCD探测器像面,形成色散条纹图像;
步骤二、采集色散条纹图像,采用条纹斜率粗提取方法计算条纹斜率,具体包括以下步骤:
采集的色散条纹图像上不同位置(x,y)的光强分布公式如下:
其中,a、b和d分别表示光阑的宽度、长度以及两孔中心间距,x表示色散条纹图像上不同位置的横坐标,y表示色散条纹图像上不同位置的纵坐标,λ表示与所采集色散条纹图像相对应的波长,λmin,λmax表示与所采集色散条纹图像相对应的波长范围,f表示成像透镜等效焦距,φ(λ)表示色散棱栅的色散角,δ表示两个拼接子镜反射波面的光程差,即两个拼接子镜之间的平移误差;
二维色散条纹图像取得暗条纹的条件为:
λ与x之间的关系为:
λ(x)=λ0+C0x (3)
其中,λ0为宽光谱中心波长,C0为色散棱栅的线色散倒数;
步骤三、若条纹斜率绝对值大于设定的阈值T1,则进行粗调整并计算两个拼接子镜之间的平移误差,重复步骤二;
步骤四、若条纹斜率绝对值小于设定的阈值T1,则进行精调整并重新计算条纹斜率;
步骤五、若条纹斜率绝对值小于设定的阈值T2,则迭代调整结束,否则重复步骤四。
2.根据权利要求1所述的基于色散条纹斜率分析的拼接镜共位相误差校正方法,其特征在于,步骤三中,采用图像二值化、连通域分析和方位角计算方法进行粗调整。
3.根据权利要求1所述的基于色散条纹斜率分析的拼接镜共位相误差校正方法,其特征在于,步骤四中,采用低通滤波和数据拟合方法进行精调整。
5.根据权利要求4所述的基于色散条纹斜率分析的拼接镜共位相误差校正方法,其特征在于,步骤四具体包括以下步骤:
若条纹斜率绝对值小于设定的阈值T1,则截取一系列不同波长位置的色散条纹数据,对每一列色散条纹数据进行低通滤波,滤掉毛刺,得到每一列色散条纹数据的峰值位置,同时取能量最强的三个峰值进行抛物线拟合,记录所拟合的抛物线顶点的位置,最终得到一系列峰值点与抛物线顶点;将所有的峰值点拟合出一条直线,将所有的抛物线顶点拟合出一条直线,峰值线相对于抛物线顶点拟合线的斜率即为色散条纹斜率的精确计算值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010047451.6A CN111207910B (zh) | 2020-01-16 | 2020-01-16 | 基于色散条纹斜率分析的拼接镜共位相误差校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010047451.6A CN111207910B (zh) | 2020-01-16 | 2020-01-16 | 基于色散条纹斜率分析的拼接镜共位相误差校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111207910A CN111207910A (zh) | 2020-05-29 |
CN111207910B true CN111207910B (zh) | 2021-05-25 |
Family
ID=70784758
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010047451.6A Active CN111207910B (zh) | 2020-01-16 | 2020-01-16 | 基于色散条纹斜率分析的拼接镜共位相误差校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111207910B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111551351B (zh) * | 2020-06-09 | 2021-08-03 | 中国科学院长春光学精密机械与物理研究所 | 一种相邻拼接镜间piston误差检测*** |
CN113324540B (zh) * | 2021-06-03 | 2023-07-07 | 合肥工业大学 | 一种基于扫描振镜的拼接镜位姿测试***及方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5317389A (en) * | 1989-06-12 | 1994-05-31 | California Institute Of Technology | Method and apparatus for white-light dispersed-fringe interferometric measurement of corneal topography |
CN107450176B (zh) * | 2017-09-21 | 2019-10-29 | 中国科学院上海天文台 | 一种空间稀疏孔径望远镜共相控制装置及控制方法 |
CN108845415B (zh) * | 2018-04-03 | 2021-02-26 | 华东交通大学 | 一种应用于拼接镜的粗共相调节的方法与检测*** |
CN108827596B (zh) * | 2018-04-03 | 2021-06-29 | 华东交通大学 | 一种应用于分块式拼接望远镜共相检测方法与装置 |
CN108519671B (zh) * | 2018-04-19 | 2020-09-18 | 中国科学院上海天文台 | 一种拼接望远镜***相位平移误差的闭环校正控制方法 |
CN109556513B (zh) * | 2018-11-23 | 2020-05-15 | 中国科学院长春光学精密机械与物理研究所 | 一种相邻拼接镜检测方法及检测*** |
-
2020
- 2020-01-16 CN CN202010047451.6A patent/CN111207910B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111207910A (zh) | 2020-05-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Rousset | Wave-front sensors | |
Mazoyer et al. | Estimation and correction of wavefront aberrations using the self-coherent camera: laboratory results | |
CN111207910B (zh) | 基于色散条纹斜率分析的拼接镜共位相误差校正方法 | |
CN111220361B (zh) | 一种测量微透镜阵列焦距的方法 | |
CN110146258B (zh) | 一种Poisson noise模型下对扩展目标成像时的相位恢复方法 | |
US20140285877A1 (en) | Calibratable beam shaping system and method | |
CN111694016B (zh) | 一种非干涉合成孔径超分辨成像重构方法 | |
US20060082882A1 (en) | Achromatic imaging lens with extended depth of focus | |
Potier et al. | Increasing the raw contrast of VLT/SPHERE with the dark hole technique-I. Simulations and validation on the internal source | |
Janin-Potiron et al. | The self-coherent camera as a focal plane fine phasing sensor | |
Li et al. | Self-measurements of point-spread function for remote sensing optical imaging instruments | |
Zepp et al. | Simulation-based design optimization of the holographic wavefront sensor in closed-loop adaptive optics | |
CN114076670A (zh) | 一种拼接主镜共相误差检测方法、***及储存介质 | |
Righini et al. | Imaging performance of multi-etalon bidimensional spectrometers | |
Volin | Portable snapshot infrared imaging spectrometer | |
CN115901192A (zh) | 具有实时对准功能的光学***波前拼接检测方法及装置 | |
Xie et al. | Restoration of sparse aperture images using spatial modulation diversity technology based on a binocular telescope testbed | |
CN113534436B (zh) | 基于七孔径的稀疏合成孔径成像***及其相位校正方法 | |
Sidick et al. | Studies of the effects of optical system errors on the HCIT contrast performance | |
Smith et al. | Coarse phasing of segmented mirrors using a dispersed fringe sensor | |
Guillard et al. | Optical performance of the JWST/MIRI flight model: characterization of the point spread function at high resolution | |
CN114323280B (zh) | 一种基于曲率传感的光谱分析装置及其方法 | |
Mugnier et al. | Inversion in optical imaging through atmospheric turbulence | |
Denolle et al. | Real-time focal-plane wavefront sensing for compact imaging phased-array telescopes: numerical and experimental demonstration | |
CN113203485B (zh) | 一种通过单次曝光实现轴向相位差波前重建的装置及方法 |
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 |