2. 中国科学院大学, 北京 100049
2. University of Chinese Academy of Sciences, Beijing 100049, China
1 m新真空太阳望远镜(New Vacuum Solar Telescope, NVST)是口径为1 m的地平式真空太阳望远镜,主要对太阳进行高分辨成像和光谱观测,自2010年投入使用以来,已取得了大量的太阳高分辨观测数据[1-2]。由于1 m太阳望远镜观测目标是强热源,虽然其光机结构经过特殊设计和处理[3],但整个系统的稳定性还是受太阳辐射变化的影响,目前主要表现为在连续观测时间内焦点的实时变化,当变化超出系统的焦深时,就会产生离焦像差,从而影响观测数据的质量[4-5],准确实时地探测1 m太阳望远镜的焦点变化并进行实时调整,将有助于提高望远镜的观测效率和数据质量。
如何实现光学系统的自动调焦已有不少学者进行了研究,文[6]提出了几种测距方法和图像清晰度的评价函数,而大型天文望远镜离焦主要是由于自身光学系统发生改变,这些调焦方法并不适用。文[7]基于图像处理采用搜索焦点的方法在南京大学浦口天文台的65 cm望远镜上以恒星为观测目标进行了自动调焦实验。文[8]提出了绝对方差的对焦评价函数,并在怀柔太阳观测基地用搜索焦点的方法以全日面像为目标进行了自动对焦实验。但这两个调焦系统不适用于1 m太阳望远镜,原因主要有两方面:(1)1 m太阳望远镜观测目标是局部日面像,其图像清晰度评价函数不适用;(2) 这两个调焦系统都是采用移动CCD的方法调节焦点,而1 m太阳望远镜终端成像系统包括多通道高分辨成像系统和多波段光谱仪,因此,1 m太阳望远镜调节焦点需要移动主光路中的M3成像镜。文[9]采用焦点扫描法,连续采集焦点前后不同扫描位置的局部太阳像,然后计算特定环上的功率谱,并进行多次平均以消除大气影响,最后通过高斯拟合功率谱与扫描位置的关系曲线从而获得焦点位置,文[9]对该方法进行了大量实验,充分论证了方法的稳定性和可靠性。本文的研究工作正是基于此方法开展的,结合1 m太阳望远镜的系统结构,为望远镜研制了焦点实时探测系统。
本文第1节根据1 m太阳望远镜的光学系统从理论和实测两方面分析了望远镜焦点位置随温度的变化;第2节根据文[9]中的焦点扫描探测法和1 m太阳望远镜的系统结构,给出了1 m太阳望远镜焦点实时探测系统总体结构,并分析了焦点探测算法对系统参数的具体要求;第3节详细设计了1 m太阳望远镜焦点实时探测系统的硬件和软件系统;第4节测试系统的重复精度和系统性能;第5节对本文研究工作进行了总结和展望。
1 1 m太阳望远镜焦点变化分析1 m太阳望远镜整体(包括建筑)和主光学系统如图 1,望远镜观测时,圆顶完全移开,除望远镜整体直接受太阳照射外,周围环境也随太阳照射而温度升高。为减小对温度的敏感性,M1与M2采用热膨胀系数很小的铟钢连接以控制M1与M2之间距离的变化,在焦点F1的焦面视场光阑处设计了热控系统。
根据文[5]的研究结果,焦点F3的位置变化与镜筒温度变化关系近似描述为(1) 式:
$ \Delta {F_3}/\Delta T = \left[{\left( {{L_{12}}{\beta _2} + {L_{24}}} \right){\alpha _1} + {L_{43}}{\alpha _2}} \right]{\beta _3}, $ | (1) |
其分析结果为:ΔF3/ΔT=2.57 mm/℃。
M3镜的位置变化量与焦点F3变化量近似是线性关系:
$ \Delta {M_3} = \Delta {F_3}/{\alpha _{\rm{t}}}, $ | (2) |
其中,αt=8.67是终端光学系统的轴向放大率,因此,1 m太阳望远镜焦点调节是通过调整M3的位置实现的,M3的调节量也同样反应望远镜的焦点变化。由(1) 式和(2) 式可得调节量与镜筒内温度变化的关系为ΔM3/ΔT=0.30 mm/℃。由此计算M3调节量随温度的理论变化值,并实测望远镜打开圆顶开始观测后温度上升较快的2小时内镜筒温度与焦点位置,理论值与实测值如图 2。
从理论计算与实测结果看,理论值能大致反应焦点随镜筒内温度的变化趋势,但二者差别较大,实测数据表明,焦点变化速度较理论分析的快,主要是理论分析只考虑了温度因素,而且有一定的简化,实测时,由于无法对真空镜筒内系统各点温度实测,只能以镜筒臂上的温度近似代替。总之,1 m太阳望远镜在开始观测的一段时间内,焦点变化较快,需要对焦点做频繁的调整。
打开圆顶开始观测以后,望远镜受太阳照射和环境温度的影响,主镜筒内热量累积使望远镜产生离焦。离焦图像与聚焦图像对比(λ=705.8 nm)如图 3。
根据文[5]的研究结果,离焦像差小于0.2 λ时,可以忽略对高分辨重建的影响,离焦像差与焦点F3的位置变化量的关系为
$ \omega = \frac{{\Delta {F_3}}}{{8\bar \lambda {{\left( {f/d} \right)}^2}}}. $ | (3) |
从(3) 式可以看出,离焦像差变化率与波长成反比,1 m太阳望远镜多通道高分辨成像系统有5个观测通道,最短观测波长为λ=393.3 nm,根据(3) 式,该波长对焦点变化最为敏感。根据(1) 和(3) 式计算λ=393.3 nm时,离焦像差随温度的变化率为Δω/ΔT=0.403 λ/℃。
根据文[5]的测量结果,在温度变化较快的时段,1小时内温度变化4.5℃,即产生近1.8 λ的离焦像差变化,因此,在焦点变化快的时候,1小时内焦点需要调节10次左右,这就要求焦点探测系统能够在6分钟以内完成一次焦点探测和焦点调节。
以上结果表明,虽然1 m太阳望远镜光机系统经过优化设计,但F3焦点仍随温度变化。该焦点探测系统研制之前,望远镜都是通过人为判断焦点变化进行调焦,这种方法效率和准确度都不高。因此,1 m太阳望远镜焦点实时探测系统是不可或缺的。
2 总体结构设计文[9]中的扫描焦点探测法要求在焦点前后位置连续扫描记录多帧图像和每帧图像对应的位置,同时为消除湍流大气的影响,需要连续扫描多组图像进行平均。1 m太阳望远镜多通道高分辨成像系统已投入观测的有Hα(λ=656.3 nm)和TiO(λ=705.8 nm)两个波段,在TiO成像通道进行分光,设计一个同步的成像系统进行图像采集,经后续算法处理后,得到焦点的位置,该位置信息转换为M3的调节量后,可以提醒观测人员进行调焦,也可以输送给调焦系统,形成闭环的自动调焦系统。1 m太阳望远镜焦点探测系统的总体结构设计如图 4。针对该算法,1 m太阳望远镜高分辨率成像系统和光学结构的具体要求和特点,对望远镜焦点探测系统各参数的详细要求如表 1。
参数名称 | 物理意义 | 数值 |
工作波长/nm | 图像采集系统工作中心波长 | 705.8 |
视场/′ | 图像采集系统的工作视场 | 1 × 1 |
焦比 | 焦点探测系统的光学系统焦比 | 25.9 |
每帧图像曝光/ms | 每帧图像的曝光时间 | <30 |
图像采集频率/Hz | 图像采集的频率 | 5 |
所需图像组数/组 | 一次焦点探测要求的图像组数 | 10 |
每组图像帧数/帧 | 每组图像(完成一次焦点扫描)所需的图像帧数 | 50 |
位置分辨率/μm | 扫描系统位置分辨率 | 1 |
位置精度/μm | 扫描系统的位置精度 | 5 |
扫描步长/μm | 扫描时连续两帧图之间的距离 | 758 |
步长运动时间/s | 扫描时连续两帧图之间的平台运动时间 | 0.2 |
扫描步数/步 | 每组图像直接扫描点,与每组图像帧数严格一致 | 50 |
扫描范围/mm | 扫描时,相机在焦点前后运动的距离 | 15.0 |
一次焦点时间/min | 完成一次焦点探测所需要的时间 | 6 |
根据(3) 式,离焦像差要求等于0.2 λ,λ=705.8 nm时,焦点探测系统的焦比f/d为25.9,ΔF3=758 μm,对应采集图像时的步长要小于758 μm,焦点扫描探测算法最后一步是高斯拟合得出焦点位置,为取得较好的拟合效果,扫描范围需要大于±2 λ,由(3) 式计算ΔF3=15.15 mm,因此,扫描范围需要大于15 mm。同时,因为要扫描多组进行平均,所以多组之间对应位置点的重复精度也是非常关键的指标。
3 软硬系统设计 3.1 硬件系统设计为了实现图像采集系统的扫描,在调焦成像系统的焦平面安放一个电动位移台,位移台上安装高精度的位置检测系统。平台上安放采集相机,图像采集和平台控制通过一台计算机完成,如图 5。
为满足表 1的指标要求,相机与位移台的选型和主要参数见表 2。系统安装时,位移台的运动方向要与光轴方向保持一致。根据表 1和表 2,位移台WN230 TA100 M的最大运动范围是10.0 cm,其位置检测系统的最小分辨率是1 μm,这相当于位移台的运动方向与光轴方向相差15.4′时带来的误差,这就要求位移台的运动方向与光轴方向相差小于15.4′。实际系统装调时,在10.0 cm范围移动位移台,用CCD记录穿过光轴中心的激光光斑的质心最大平移,然后折算成位移台的运动方向与光轴方向的夹角,将这一误差调节在5′以内,完全满足15.4′的要求。相机和位移台在焦点前后的位置与其它高分辨成像通道采用齐焦进行定标,保证焦点总能处在系统扫描范围内。
相机型号与主要参数 | 位移台型号与主要参数 | |||
厂家信息 | JAI集团 | 厂家信息 | 乐创自动化公司 | |
相机型号 | BM-141 GE | 位移台型号 | WN230 TA100 M | |
芯片类型 | CCD | 运动范围 | 10 cm | |
靶面大小 | 1392×1040 pixels | 编码器分辨率 | 1 μm | |
像元尺寸 | 6.45 μm | 控制电机类型 | 步进电机 | |
AD位数 | 8/12位 | 电机控制器型号 | MPC08 | |
读出速度 | 30帧/秒 | 系统控制模式 | 位置式/速度式 | |
曝光时间范围 | 63μs~33 ms | 零点类型 | 光耦 | |
读出接口 | 以太网 | 限位开关类型 | 机械式接近开关 |
软件系统主要完成相机控制与数据采集、相机与位移台的同步控制和数据处理并计算焦点位置等功能,软件系统的界面和工作流程分别如图 6和图 7。软件在WIN7平台上采用VC++开发,系统采用多线程实现[10],主要包括3个线程:(1) 系统主线程,包括界面和消息响应的实现;(2) 图像采集、平移台控制和数据处理线程;(3) 图像显示线程。系统完成一次检测焦点需要6 min。
4 位移台和焦点探测系统性能测试 4.1 位移台性能测试位移台性能测试主要测试位移台的重复移动精度是否满足要求。平台移动有两种工作模式:位置模式和速度模式,位置工作模式是给位移台一个指定的扫描起点和步进距,等位移台完成一步运动后保存平台的位置并采集一帧图像;速度模式是给位移台一个恒定的运动速度,移动过程中连续采集图像并记录平台的位置。为平滑大气的影响,需要连续采集多组图像进行平均,因此,组与组之间的位置精度需要满足算法要求,即平台的重复移动精度。测试结果如图 8。
测试时需满足连续两帧之间的位置差小于ΔF(757 μm),因此,对位置模式进行测试时,步长设定为300 μm,扫描一次采集50个点的位置数据,来回扫描10次,对应每个位置都采集了10个数据,然后计算10个位置数据均方根值,得到一个点的位置重复精度,逐个计算50个点位置均方根值,如图 8(a)。50个点的位置均方根值中,最大为3.03 μm,最小为0.78 μm,其重复移动精度满足系统精度要求。对速度模式进行测试时,位移台移动速度为2500 μm/s,每组同样采集50个点的位置数据,来回扫描10次,然后计算50个点的位置均方根值,如图 8(b)。50个点的位置均方根值中,最大为8.24 μm,最小为5.72 μm。通过比较可以看出,速度模式下位置数据的均方根值较位置模式下的均方根值大,因此,位移台工作在位置模式的重复移动精度较工作在速度模式的重复移动精度高,系统更稳定。
4.2 焦点探测系统测试焦点探测系统实验主要测试整个系统对焦点的探测精度,包括位移台、CCD采集系统、图像处理、焦点探测算法的一个综合测试。实验过程如下:
(1) 将1 m太阳望远镜对准日面上的任意观测区域,采集调焦系统中CCD的平场和暗场;
(2) 通过1 m太阳望远镜的M3调焦系统,人为随机调节M3位置,即改变望远镜的焦点,同时记录M3调节量,并根据(2) 式转换为焦点探测系统中的焦点变化量;
(3) 打开焦点探测系统,测量并记录一次焦点变化量;
(4) 反复进行(2) 和(3) 操作,多次测量;
(5) 计算多次测量中焦点实测变化量与给定变化量的偏差,并统计均方根值。
图 9是进行47次实验得到的焦点探测偏差分布图,均方根值是490 μm。
从焦深角度看,当离焦距离处在焦深内时系统成像是清晰的,而焦点探测的均方根值小于半焦深时则认为是可信的。1 m太阳望远镜的焦深计算公式如下:
$ \delta = \pm 2\bar \lambda {\left( {f/d} \right)^2}. $ | (4) |
把系统的焦比f/d=25.9和波长λ=705.8 nm代入(4) 式得到半焦深为947 μm,焦点探测系统偏差的均方根值为490 μm,因此从焦深角度,该调焦精度能够满足要求。
从天文高分辨统计重建的角度看,探测精度应小于0.2 λ。把f/d=25.9和λ=705.8 nm,ΔF3=490 μm代入(3) 式得到焦点探测偏差带来的离焦像差为0.128 λ,因此从天文高分辨统计重建的角度来看也是可行的。
根据对系统的探测精度和时间性能的测试,系统能有效探测1 m太阳望远镜的焦点变化,系统于2016年8月份投入1 m太阳望远镜的常规观测中使用,使用调焦系统的观测数据可参考1 m太阳望远镜的数据网站。系统投入使用后,还可以有效监测1 m太阳望远镜的焦点变化情况,图 10是1 m太阳望远镜在4天中的焦点调节(M3位置变化)情况,从图 10调焦量随时间的变化可以看出,在望远镜对准太阳后,焦点变化较快,需要不断调焦,开始观测后两个小时左右,焦点变化趋于稳定,在中午11:00到13:00,焦点变化在一个方向达到最大,随后往反方向变化,但变化相对于刚开机观测时缓慢,也需要不断调焦,这是由于太阳高度引起的辐射变化也是造成望远镜光机系统温度变化的原因之一。
5 总结本文针对1 m太阳望远镜在观测过程中焦点实时变化的问题,根据文[9]提出的焦点探测方法并结合1 m太阳望远镜的光机结构,设计了焦点探测系统。系统经过测试并已经投入使用,3个月运行结果表明:系统能准确地探测1 m太阳望远镜的焦点变化,并提供给观测人员焦点调焦量,提高了望远镜的观测效率和数据质量。下一步工作需要改造1 m太阳望远镜的焦点调节部分,与焦点探测系统形成闭环控制,进一步提高望远镜的工作效率。
[1] | Liu Zhong, Xu Jun, Gu Bozhong, et al. New vacuum solar telescope and observations with high resolution[J]. Research in Astronomy and Astrophysics, 2014, 14(6): 705–718. DOI: 10.1088/1674-4527/14/6/009 |
[2] | Xu Zhi, Jin Zhenyu, Xu Fangyu, et al. Primary observations of solar filaments using the multi-channel imaging system of the New Vacuum Solar Telescope[C]//Proceedings of the International Astronomical Union. 2013:117-120. http://journals.cambridge.org/abstract_S1743921313010831 |
[3] |
顾伯忠, 左恒. 云南天文台1m红外太阳望远镜的主镜热力学分析[J]. 天文研究与技术——国家天文台台刊, 2008, 5(1): 83–90 Gu Bozhong, Zuo Heng. The thermal analysis of the primary mirror of 1m infrared solar telescope at Yunnan Observatory[J]. Astronomical Research & Technology——Publications of National Astronomical Observatories of China, 2008, 5(1): 83–90. |
[4] |
向永源, 刘忠, 金振宇, 等. 高分辨率太阳图像重建方法[J]. 天文学进展, 2016, 34(1): 94–110 Xiang Yongyuan, Liu Zhong, Jin Zhenyu, et al. High resolution solar image reconstruction[J]. Progress in Astronomy, 2016, 34(1): 94–110. |
[5] |
方玉亮, 金振宇, 刘忠, 等. 一米新真空太阳望远镜离焦对高分辨太阳观测图像重建的影响[J]. 天文研究与技术, 2015, 12(2): 183–188 Fang Yuliang, Jin Zhenyu, Liu Zhong, et al. A study of influences of defocus aberrations on high-resolution image reconstruction for data from the new vacuum solar telescope of the YNAO[J]. Astronomical Research & Technology, 2015, 12(2): 183–188. |
[6] |
梁翠萍, 李清安, 乔彦峰, 等. 简析光学系统的自动调焦方法[J]. 电光与控制, 2006, 13(6): 93–96 Liang Cuiping, Li Qingan, Qiao Yanfeng, et al. On auto-focusing technology of optical system[J]. Electronics Optical & Control, 2006, 13(6): 93–96. |
[7] |
李晓燕, 朱庆生. 一种基于图像清晰度评价的天文望远镜自动调焦系统[J]. 天文研究与技术——国家天文台台刊, 2008, 5(3): 294–298 Li Xiaoyan, Zhu Qingsheng. An automatic focusing system of astronomical telescope based on image definition evaluation[J]. Astronomical Research & Technology——Publications of National Astronomical Observatories of China, 2008, 5(3): 294–298. |
[8] | 郭晶晶. 太阳望远镜的高精度自动导行和自动调焦方法的研究[D]. 昆明: 昆明理工大学, 2015. http://cdmd.cnki.com.cn/Article/CDMD-10674-1015641185.htm |
[9] | 方玉亮. NVST焦点探测[D]. 昆明: 中国科学院云南天文台, 2014. http://ir.ynao.ac.cn/handle/114a53/4911 |
[10] |
张宏, 黄小诚. VC++语言多线程编程及其实现[J]. 科技信息, 2008(7): 64–65 Zhang Hong, Huang Xiaocheng. Multithread programming in VC++[J]. Science & Technology Information, 2008(7): 64–65. |