技术领域
[0001] 本
发明涉及光学与天文学领域，特别涉及一种基于二阶压缩感知的超灵敏光谱成像天文望远镜及方法。
背景技术
[0002] 天文望远镜是观测天体的重要工具，可以毫不夸大地说，没有望远镜的诞生和发展，就没有现代天文学。随着望远镜在各方面性能的改进和提高，天文学也正经历着巨大的飞跃，迅速推进着人类对
宇宙的认识。
[0003] 按工作波段的不同，天文望远镜可分为
光学望远镜和射电望远镜。其中光学望远镜主要以可见光为工作波段，根据使用地点的不同，可以分为地面天文望远镜和空间天文望远镜。由于光学系统的不同，又可以分为反射望远镜，折射望远镜，折反射望远镜等类型；射电望远镜主要以
无线电波为工作波段。目前地面观测绝大多数处于凝聚态的天体(
恒星等)仍以光学波段观测为主要手段，这是由于：大多数恒星等天体
温度范围从数千度到数万度，
辐射集中在光学波段；携带大量天体物理信息的谱线，主要集中于可见区；大气在可见区有良好的透射。
[0004] 在天文观测中，光谱信息的获取具有重要的意义，这是因为天文学中大量的信息能以光谱的形式表现出来。第一，对宇宙和
星系的研究。宇宙的诞生、星系的形成等前沿问题都建立在对星系物理的研究
基础之上。研究宇宙大尺度结构依赖于星系红移巡天的工作。获取星系的光谱就能得到星系的红移，进而知道它的距离，由此获得星系的三维分布，这样就可以了解整个宇宙空间的结构，同时可以研究包括星系的形成、演化在内的宇宙大尺度结构和星系物理。获取星系的光谱是进行这一工作最基础的需要。第二，对恒星和
银河系的结构特征的研究。由于不同元素具有不同的特征谱线，通过一颗恒星的光谱，可以分析出其元素构成和含量等化学组成，可以分析出其
密度、温度等物理条件，还可以测量出其运动速度和运行轨迹等。研究了不同种类的恒星的分布，可以研究出
银河系的结构和银河系的形成。第三，对宇宙生命的研究。通过恒星或行星的光谱，可以研究其表面
水分和
氧气的含量，以确定是否存在
生物的可能。因此，在天文学中对光谱的研究具有重要而不可替代的作用。
[0005] 然而，天文望远镜要想同时获得天文图像和天文光谱信息是十分困难的，其中最主要的困难是维度的问题。二维的天文图像与一维的光谱信息共有三维信息，按照传统的信息获取方式，则需要具有三个维度的探测器，而目前这显然是无法实现的，因此现有的大量的天文望远镜只能分别获得天文图像信息或天文光谱信息，而无法同时获得两方面的信息。一种解决方法是在普通的天文望远镜上通过二维探测器获取图像信息，再通过滤光片等方式滤出某一感兴趣的波段的光
信号进行成像，这样可以获得单一波段的光谱成像，而要获得多波段或全波段的光谱成像只能通过改变滤光系统进行重复测量，而获得不同波段的图像。这种光谱成像的方式需要通过扫描光谱实现，要获得高
分辨率的光谱，必然会带来巨大的时间成本，而且本质上仍然无法实现天文图像信息和天文光谱信息的同时获取。
[0006] 灵敏度是天文望远镜非常重要的指标，因为天文望远镜灵敏度提高，就能够看到更暗更远的天体，这等同于能够看到更早期的宇宙，这对于研究宇宙的起源等人类关心的最基本问题具有重要意义。在天文光谱探测中，由于只获取单一波段的信息，与全波段成像时相比
光信号的强度大大减弱，因此对灵敏度的要求更高。天文光谱望远镜的灵敏度提高，就可以将光谱测量时的
波长分得更细，获得更高的光谱分辨率。因此天文成像和天文光谱成像的发展需要更高灵敏度的天文望远镜。
[0007] 目前天文望远镜灵敏度的提高主要通过口径的增大来实现，望远镜的口径越大，集光能
力越强，灵敏度也会越高，因此现代天文望远镜的口径制作地越来越大。但是，随着望远镜口径的增大，一系列的技术问题接踵而来。例如，口径5米的海尔望远镜曾经是世界上最大的天文望远镜，它的镜头自重达14.5吨，可动部分的重量为530吨，而后来建成的6米口径天文望远镜更是重达800吨。一方面，望远镜的自重过大会使镜头
变形相当明显，另一方面，镜体温度不均也令镜面产生畸变，进而影响成像
质量。从制造方面看，传统方法制造望远镜的
费用几乎与口径的平方或立方成正比，所以制造更大口径的望远镜在性能和费用上都受到极大限制。
[0008] 影响天文望远镜灵敏度的另一个重要因素在于光学探测器的性能，高灵敏度的探测器必然可以有效提高天文望远镜的灵敏度。基于盖革工作模式的
雪崩光电
二极管(APD)可以探测到单个
光子的
能量，是理论上灵敏度最高的探测器，也称为单光子探测器。其他高灵敏度探测器还包括
光电倍增管(PMT)，其灵敏度可达到几个或几十个光子。现阶段我国只有单点的APD探测器，由于制作工艺的限制还没有生产阵列APD的能力；国际上可用的阵列APD最大
像素也只有128×128，远远达不到获得高分辨率天文图像的需求。另外，PMT由于工作机理的原因也没有阵列探测器。对高灵敏度探测器像素数不足问题的解决，一种办法是使用点探测器进行扫描实现成像，这样带来的问题是扫描探测器会耗费大量的时间，大大降低图像获取速度，同时图像不同
位置的信息探测时间产生差异，扫描期间的图像偏移会造成成像分辨率的下降。另一种办法是将大量点探测器拼成阵列进行探测，但是要获得足够的分辨率，需要数量极其庞大的单点探测器，如要获得1024×768像素的图像则需要大约80万的点探测器，造成极高的成本，并且点探测器拼接会存在严重的占空比问题，造成集光效果的下降，进而影响望远镜的灵敏度。因此，利用
现有技术这些高灵敏度的探测器都无法解决光谱成像天文望远镜存在的探测维度的问题，无法同时获取天文图像和天文光谱信息。
[0009] 综上所述，现有的光谱成像天文望远镜存在图像和光谱探测信息维度过大、探测器维度不足的问题，而且无法实现高灵敏度的探测。由于工作原理的局限，传统的光谱成像天文望远镜在实现多维探测和提高探测灵敏度的途径上存在制约，天体物理的发展亟需灵敏度更高的光谱成像天文望远镜。
发明内容
[0010] 本发明的目的在于以二阶压缩感知为基础，克服现有技术中的光谱成像天文望远镜在多维探测和灵敏度上的不足，从而提供一种基于二阶压缩感知的超灵敏光谱成像天文望远镜。
[0011] 为了实现上述目的，本发明提供了一种基于二阶压缩感知的超灵敏光谱成像天文望远镜，包括光学单元I和电学单元II；其中，光学单元I包括天文望远镜镜头1、第一空间光
调制器2、
准直部件3、光谱分光部件4、光谱会聚部件5、第二空间光调制器6、收集部件7；电学单元II包括单光子点探测器8、计数器9、随机数发生器10、控
制模块11、数据包
存储器12以及压缩感知模块13；其中，所述准直部件3包括收集透镜3_1、光阑3_2、
准直透镜3_3；所述随机数发生器10包括第一随机数发生器10_1、第二随机数发生器10_2；
[0012] 从天体传播而来的单光子水平的光信号由所述天文望远镜镜头1收集，并成像到所述第一空间光调制器2上；所述第一空间光调制器2对成像在其表面的天文图像进行随机调制，以随机概率将图像上不同位置的光反射到所述准直部件3方向；所述第一空间光调制器2随机反射的光首先由所述收集透镜3_1会聚到所述光阑3_2，限制光斑尺寸，形成近似点
光源，然后经过所述准直透镜3_3准直形成平行光，照射在所述光谱分光部件4上；所述光谱分光部件4将不同波长的光向不同方向出射；经过所述光谱会聚部件5后不同波长的光会聚到所述光谱会聚部件5焦平面上不同位置，形成光谱带；所述第二空间光调制器6放置在所述光谱会聚部件5的焦平面位置，对其表面的光谱带进行随机调制，以随机概率将光谱带上不同位置的光反射到所述收集部件7方向；所述收集部件7收集所述第二空间光调制器6反射而来的光信号，由所述电学单元II中单光子点探测器8探测；
[0013] 所述第一随机数发生器10_1产生随机数用于控制所述第一空间光调制器2、所述第二随机数发生器10_2产生随机数用于控制所述第二空间光调制器6；所述第一空间光调制器2和所述第二空间光调制器6根据该随机数实现对光信号的随机调制；所述单光子点探测器8探测待测极弱光中的单光子，将采集到的单光子信号转换成脉冲形式的
电信号后输出；所述计数器9记录所述单光子点探测器8发出的代表单光子数目的电脉冲个数；所述
控制模块11对整个超灵敏光谱成像天文望远镜进行控制协调，包括对各部件的工作控制和同步脉冲触发信号发射，确保所述计数器9、所述第一空间光调制器2和所述第二空间光调制器6同步工作；所述计数器9所记录的单光子数目和所述第一随机数发生器10_1、所述第二随机数发生器10_2生成的随机矩阵全部存入所述数据包存储器12中；所述压缩感知模块13利用所述数据包存储器12中的单光子数目以及对应的随机矩阵，并选取稀疏基对天文光谱图像进行重建，得到极弱光水平的天文光谱图像。
[0014] 上述技术方案中，所述第一随机数发生器10_1、第二随机数发生器10_2用于生成二值伯努利分布的散斑或二值非均匀分布的散斑，二值由0和1组成；其中，当第二随机数发生器10_2生成二值伯努利分布的散斑时，需使第一
帧的散斑全1，而伯努利分布由Walsh或Hadamard或noiselet变换获得；当第二随机数发生器10_2生成二值非均匀分布的散斑时，每帧散斑中1的个数需远小于0的个数，且1在每帧散斑的空间分布上是随机的；第一随机数发生器10_1生成的散斑不做限定。
[0015] 上述技术方案中，所述天文望远镜镜头1采用下列任意一种天文望远镜类型的镜头：反射式天文望远镜，包括
牛顿式、卡塞格林式、格里式；折射式天文望远镜，包括伽利略望远镜、开普勒望远镜；折反射式天文望远镜，包括施密特-卡塞格林式、
马克苏托夫-卡塞格林式；多镜面望远镜；双筒望远镜；也包括应用于卫星、空间站之上的空间天文望远镜。
[0016] 上述技术方案中，所述第一空间光调制器2或第二空间光调制器6采用数字微镜器件实现。
[0017] 上述技术方案中，所述准直部件3中的收集透镜3_1，准直透镜3_3通过透镜或凹面镜实现；光阑3_2通过狭缝或小孔实现。
[0018] 上述技术方案中，所述光谱分光部件4包括用于将不同波长的光分开的色散分光部件，所述色散分光部件采用包括光栅、棱镜在内的具有分光能力的器件实现。
[0019] 上述技术方案中，所述光谱分光部件4还包括用于滤除不进行光谱成像的波长的光的预滤光部件，所述预滤光部件由滤光片实现。
[0020] 上述技术方案中，所述光谱会聚部件5由透镜或凹面镜实现；所述光谱会聚部件5将不同波长的光按波长从小到大依次透射到所述第二空间光调制器6的不同像素上。
[0021] 上述技术方案中，所述第一空间光调制器2稳定于每一次随机调制时间内时，所述第二空间光调制器6进行若干次随机调制。
[0022] 上述技术方案中，所述收集部件7由透镜或凹面镜实现。
[0023] 上述技术方案中，所述单光子点探测器8采用盖革模式
雪崩二极管或光电倍增管实现。
[0024] 上述技术方案中，所述控制模块11确保所述计数器9和所述第一空间光调制器2、所述第二空间光调制器6之间同步工作包括：所述第一空间光调制器2进行若干次随机调制；当所述第一空间光调制器2稳定于每一次随机调制时间内，所述第二空间光调制器6进行若干次随机调制；所述第二空间光调制器6每进行一次随机调制，所述计数器9分别累积所述单光子点探测器8发出的代表单光子数目的电脉冲个数，直到所述第二空间光调制器6进行下一次随机调制，所述计数器9将所述第二空间光调制器6稳定于一次随机调制时间内的光子计数传输至数据包存储器12，并将计数清零，开始下一次计数。
[0025] 上述技术方案中，所述压缩感知模块13采用下列
算法中的任意一种实现压缩感知：匹配
跟踪算法MP、
正交匹配跟踪算法OMP、基跟踪算法BP、贪心重建算法、LASSO、LARS、GPSR、贝叶斯估计算法、magic、IST、TV、StOMP、CoSaMP、LBI、SP、l1_ls、smp算法、SpaRSA算法、TwIST算法、l0重建算法、l1重建算法、l2重建算法；稀疏基采用离散余弦变换基、小波基、傅里叶变换基、梯度基、gabor变换基中的任意一种；当所测天文图像本身具有很好的稀疏性时，不通过稀疏基的变化，直接对原始信号进行重建。
[0026] 本发明还提供了基于二阶压缩感知的超灵敏光谱成像天文望远镜所实现的天文光谱图像获取方法，包括：
[0027] 步骤1)光信号获取的步骤：
[0028] 从天体传播而来的单光子水平的光信号由所述天文望远镜镜头1收集，并成像到所述第一空间光调制器2上；所述第一空间光调制器2对成像在其表面的天文图像进行随机调制，以随机概率将图像上不同位置的光反射到所述准直部件3方向；所述第一空间光调制器2随机反射的光首先由所述收集透镜3_1会聚到所述光阑3_2，限制光斑尺寸，形成近似点光源，然后经过所述准直透镜3_3准直形成平行光，照射在所述光谱分光部件4上；所述光谱分光部件4将不同波长的光向不同方向出射；经过所述光谱会聚部件5后不同波长的光会聚到所述光谱会聚部件5焦平面上不同位置，形成光谱带；所述第二空间光调制器6放置在所述光谱会聚部件5的焦平面位置，对其表面的光谱带进行随机调制，以随机概率将光谱带上不同位置的光反射到所述收集部件7方向；所述收集部件7收集所述第二空间光调制器6反射而来的光信号，由所述电学单元II中单光子点探测器8探测；
[0029] 步骤2)光学调制与单光子探测、计数同步工作的步骤；
[0030] 所述第一随机数发生器10_1产生随机数用于控制所述第一空间光调制器2、所述第二随机数发生器10_2产生随机数用于控制所述第二空间光调制器6；所述第一空间光调制器2根据随机数在图像维度上实现对光信号的随机调制，所述第二空间光调制器6根据随机数在光谱维度上实现对光信号的随机调制；所述单光子点探测器8探测待测极弱光中的单光子，将采集到的单光子信号转换成脉冲形式的电信号后输出；所述计数器9记录所述单光子点探测器8发出的代表单光子数目的电脉冲个数；所述控制模块11对整个超灵敏光谱成像天文望远镜进行控制协调，包括对各部件的工作控制和同步脉冲触发信号发射，确保所述计数器9、所述第一空间光调制器2和所述第二空间光调制器6同步工作；
[0031] 步骤3)单光子数目与随机矩阵预处理的步骤；
[0032] 对第一随机数发生器10_1所对应的随机矩阵不做预处理；
[0033] 当第二随机数发生器10_2生成二值伯努利分布的散斑时，若第1帧为全1，将计数器9上对应某个波长的一个像素上的所有的单光子数目组成一个列向量y，维度为m×1，m为总的测量数，将第二空间光调制器6上的随机矩阵记作A，维度为m×n，n为总的信号长度，令第一帧所对应的单光子数目为y1，则使2y-y1作为新的单光子数目，2A-1作为新的随机矩阵；
[0034] 当第二随机数发生器10_2生成二值非均匀分布的散斑时，跳过此步骤3)；
[0035] 步骤4)压缩感知光谱图像恢复的步骤；
[0036] 所述计数器9所记录的单光子数目和所述第一随机数发生器10_1、所述第二随机数发生器10_2生成的随机矩阵全部存入所述数据包存储器12中；所述压缩感知模块13首先利用所述数据包存储器12中的单光子数目以及对应的所述第二随机数发生器10_2生成的随机矩阵进行压缩感知算法重建，获得所述第一空间光调制器2各次随机调制下的光谱曲线；然后利用各条光谱曲线中代表某一波长的光子数值与对应的所述第一随机数发生器10_1生成的随机矩阵进行压缩感知算法重建，获得某一波长的天文图像；分别对不同波长的天文图像进行重建，得到极弱光水平的天文光谱图像。
[0037] 上述技术方案中，在步骤1)之前还包括对所述第二空间光调制器6上各个像素对应波长进行标定的步骤；该步骤包括：
[0038] 选取几个特定波长的
激光器发射特定波长的光或者用滤光片从宽谱光源滤出某些特定波长的光，然后将这些特别波长的光从天文望远镜镜头照射进光学系统，通过对所述第二空间光调制器6的调制状态进行控制，利用所述单光子点探测器8测量所述第二空间光调制器6上各个像素的光子数，光子数分布最大值的像素位置对应这些特定波长；其他像素位置对应的波长可按照线性分布计算得出。
[0039] 上述技术方案中，在步骤1)之前还包括减少仪器噪声的步骤；该步骤包括：对仪器进行密闭封装，或提高光学部件的透过率，或提高仪器内部的清洁度，或提高光谱分光系统4的效率，或提高单光子点探测器8的包括探测效率、暗计数在内的参数，或提高仪器
稳定性。
[0040] 上述技术方案中，在步骤1)之前还包括采用主动光学或自适应光学提高图像
信噪比的步骤；其中，所述主动光学通过促动器主动改变主镜镜面的形状，修正由于重力、温度和
风力造成的镜面本身的形变对成像带来的影响，减少由此产生的光学畸变；所述自适应光学需要首先检测波前扭曲情况，然后通过安装在望远镜焦面后方的携带有促动器的小型可变形镜面对波前实时进行矫正，从而修复大气
湍流等因素对光波波前的扭曲。
[0041] 本发明的优点在于：
[0042] 1、本发明利用了最新的数学研究成果——压缩感知理论，只需要单光子点探测器即可获取一维光谱、二维图像共三维信息，实现高分辨率的天文光谱图像观测，解决了现阶段高灵敏度光谱成像中信息维度过高，探测器维度不足的问题；
[0043] 2、天文光谱图像获取过程中单光子点探测器位置可保持位置固定，既不需对图像进行扫描，也不需对光谱进行扫描，减少了机械移动产生的误差；每一次测量过程单光子点探测器均会获取天文图像整体区域以及光谱整体范围的信息，保证了图像不同位置和不同波长的信息获取的均匀性；
[0044] 3、每次测量中空间光调制器(图像)2和空间光调制器(光谱)6上的光信号各有一半收集到单光子点探测器，因此每次测量单光子点探测器都可获取整幅图像所有波长光信号的1/4，远高于对图像进行扫描或对光谱进行扫描时单次测量的光信号强度，同样远高于使用单光子线阵或阵列探测器进行图像和光谱探测时单个像素上的光信号强度，是一种高通量、高信噪比的测量方式，由此允许光谱成像时获取更小尺度的图像信息和光谱信息，可实现高图像分辨率、高光谱分辨率的高灵敏度光谱图像探测；
[0045] 4、压缩感知理论允许亚
采样的采样数，本发明中在图像探测和光谱探测两阶段均利用了压缩感知理论，在两阶段均实现了亚采样，测量次数远小于单光子点探测器扫描模式实现图像探测和光谱探测的测量次数可以利用更短的时间获取天文光谱图像；
[0046] 5、本发明利用单光子点探测器实现了远高于现有天文望远镜的灵敏度，从根本上解决以往从提高望远镜口径的途径提高天文望远镜灵敏度的方式，不需要超大口径的望远镜镜头即可实现高灵敏度的天文图像探测，适当尺寸的望远镜口径可以提高镜头的均匀性和光学、力学性能，提高成像
精度和分辨率；
[0047] 6、本发明中的基于二阶压缩感知的超灵敏光谱成像天文望远镜可以广泛应用在地面、空间等工作条件下的天文望远镜，对于天文学、
宇宙学、天体物理等领域的发展有重要作用。
附图说明
[0048] 图1是本发明的超灵敏天文望远镜的结构示意图；
[0049] 图面说明
[0050] I 光学单元
[0051] 1天文望远镜镜头 2 空间光调制器(图像)
[0052] 3 准直部件3_1收集透镜
[0053] 3_2光阑 3_3准直透镜
[0054] 4光谱分光部件 5光谱会聚部件
[0055] 6 空间光调制器(光谱) 7收集部件
[0056] II电学单元
[0057] 8单光子点探测器 9计数器
[0058] 10随机数发生器10_1随机数发生器(图像)
[0059] 10_2随机数发生器(光谱) 11控制模块
[0060] 12数据包存储器13压缩感知模块
[0061] 图2是数字微镜器件中单个微镜的反射机制描述图。
[0062] 图3是利用单光子线阵探测器上不同像素的光子数实现光谱成像的示意图；
[0063] 图面说明
[0064] 6单光子线阵探测器 7计数器
[0065] 8随机数发生器
具体实施方式
[0066] 现结合附图对本发明作进一步的描述。
[0067] 本发明的基于二阶压缩感知的超灵敏光谱成像天文望远镜利用了压缩感知(Compressive Sensing，简称CS)原理，所述的压缩感知原理是由Donoho、Tao和Candès等人提出的一个全新数学理论。按照压缩感知，通过对信号进行随机采样的方式，可以利用远低于奈奎斯特/香农采样定理要求的采样数实现对信号信息的采样，并通过数学算法完美地恢复原始信号，且具有很高的鲁棒性。压缩感知主要分为三步骤：压缩采样、
稀疏变换与算法重建；其中，压缩采样，是指以少于信号数量的测量数对信号进行采样的过程y＝Ax，其中x为待测信号，A为测量矩阵，y为测量值。同时对信号的线性随机采样可以对探测维度进行压缩，只需要低于原始信号维度探测器即可获取信号的线性
叠加信息。所述的稀疏变换是选取适当的稀疏基Ψ，使得x经Ψ作用所得值x’是稀疏的，即x在Ψ
框架下可稀疏表达；所述的算法重建是在已知测量值y、测量矩阵A和稀疏基Ψ的条件下求解y＝AΨx'+e的过程，最后再由 反演出x。
[0068] 参考图1，本发明的基于压缩感知原理的超灵敏光谱成像天文望远镜包括光学单元I和电学单元II；其中，光学单元I包括天文望远镜镜头1、第一空间光调制器2、准直部件3、光谱分光部件4、光谱会聚部件5、第二空间光调制器6、收集部件7；电学单元II包括单光子点探测器8、计数器9、随机数发生器10、控制模块11、数据包存储器12以及压缩感知模块13；
[0069] 在光学单元I中，从天体传播而来的单光子水平的光信号由天文望远镜镜头1收集，并成像到第一空间光调制器2上，天文望远镜镜头的成像面大小应与第一空间光调制器2的有效面积相当，使得第一空间光调制器2的有效面积上完全
覆盖图像信息，同时天文望远镜镜头所成图像不会超出第一空间光调制器2的有效面积之外；第一空间光调制器2对成像在其表面的天文图像进行随机调制，以随机概率将图像上不同位置的光反射到准直部件3方向；准直部件3包括收集透镜3_1，光阑3_2，准直透镜3_3；第一空间光调制器2随机反射的光首先由收集透镜3_1会聚到光阑3_2，限制光斑尺寸，形成近似点光源，然后经过准直透镜3_3准直形成平行光，照射在光谱分光部件4上；光谱分光部件4将不同波长的光向不同方向出射；经过光谱会聚部件5后不同波长的光会聚到光谱会聚部件5焦平面上的不同位置，形成光谱带；第二空间光调制器6放置在光谱会聚部件5的焦平面位置，对其表面的光谱带进行随机调制，以随机概率将光谱带上不同位置的光反射到收集部件7方向；收集部件7收集第二空间光调制器6反射而来的光信号，由电学单元II中单光子点探测器8探测；
[0070] 在电学单元II中，所述随机数发生器10包括第一随机数发生器10_1，第二随机数发生器10_2，分别产生随机数用于控制所述第一空间光调制器2和所述第二空间光调制器6；所述第一空间光调制器2根据随机数在图像维度上实现对光信号的随机调制，所述第二空间光调制器6根据该随机数在光谱维度上实现对光信号的随机调制；所述单光子点探测器8探测待测极弱光中的单光子，将采集到的单光子信号转换成脉冲形式的电信号后输出；所述计数器9记录所述单光子点探测器8发出的代表单光子数目的电脉冲个数；所述控制模块11对整个超灵敏光谱成像天文望远镜进行控制协调，包括对各部件的工作控制和同步脉冲触发信号发射，确保所述计数器9、所述第一空间光调制器2和所述第二空间光调制器6同步工作；所述计数器9所记录的单光子数目和所述第一随机数发生器10_1、所述第二随机数发生器10_2生成的随机矩阵全部存入所述数据包存储器12中；所述压缩感知模块13利用所述数据包存储器12中的单光子数目以及对应的随机矩阵，并选取适当的稀疏基对天文光谱图像进行重建，得到极弱光水平的天文光谱图像。
[0071] 以上是对本发明的基于二阶压缩感知的超灵敏光谱成像天文望远镜的总体结构的描述，下面对基于二阶压缩感知的超灵敏光谱成像天文望远镜中各个部件的具体实现做进一步的描述。
[0072] 所述天文望远镜镜头1用于收集从天体发射并传播到望远镜位置的光子信号，并对天体进行成像。天文望远镜的成像分辨率和像差、色差等图像质量主要由天文望远镜镜头决定。天文望远镜镜头的结构可以采用下列任意一种天文望远镜类型的镜头：反射式天文望远镜，包括牛顿式、卡塞格林式、格里式等；折射式天文望远镜，包括伽利略望远镜、开普勒望远镜等；折反射式天文望远镜，包括施密特-卡塞格林式、马克苏托夫-卡塞格林式等；多镜面望远镜；双筒望远镜；也包括应用于卫星、空间站之上的空间天文望远镜。
[0073] 所述第一空间光调制器2和所述第二空间光调制器6均属于空间光调制器(SLM)，它能将信息加载于一维或两维的光场上。这类器件可在随时间变化的电驱动信号或其他信号的控制下，改变空间上光分布的振幅或强度、
相位、偏振态以及波长，或者把非相干光转
化成相干光，是实时光学信息处理、光学计算、光学神经网络和自适应光学等现代光学领域的关键器件，其种类有很多种，主要有数字微镜器件(Digital Micro-mirror Device，简称DMD)、
液晶光
阀、毛玻璃等。在本
实施例中，所述SLM为数字微镜器件，包括微镜阵列和集成
电路部分。在其他实施例中，也可以是其它类型的SLM。第一空间光调制器2和所述第二空间光调制器6可以是相同类型的器件，如都采用DMD，也可以是不同类型的器件，如第一空间光调制器2采用DMD，第二空间光调制器6采用液晶光阀。
[0074] 本实施例中所采用的DMD是包含有大量安装在
铰链上的微镜的阵列(主流的DMD由1024×768的阵列构成)，每一镜片的尺寸为14μm×14μm，并可对每个像素上的光实现独立控制。通过对每一个镜片下的存储单元以二进制信号进行
电子化寻址，便可让每个镜片在静电作用下向两侧翻转10～12°左右(本实施例中为+12°和-12°)，把这两种状态记为1和0，分别对应“开”和“关”，当镜片不工作时，它们处于0°的“停泊”状态。
[0075] 在图2中，对DMD中的单个微镜的反射机制做了描述。图中矩形表示DMD微镜，0°位置为微镜初始位置。图中标出了微镜处于初始位置时的法线方向，及光线入射、出射方向。当微镜处于+12°翻转状态时，微镜顺
时针旋转+12°，法线方向也随之旋转+12°。
根据反射定律，反射光将顺时针旋转24°；同理，当微镜处于-12°翻转状态时，反射光将逆时针旋转24°。因此，两个方向的反射光成48°夹
角。当准直部件3或收集部件7处于+12°或-12°反射方向时，不会收集到向另一个方向反射的光子，即可实现随机将DMD上不同位置的光收集进光路。
[0076] 所述准直部件3用于准直第一空间光调制器2随机调制的光，使其成为平行光并提供给所述光谱分光部件4，光的平行度越高，光谱分光的分辨率越高。所述收集透镜3_1将光会聚到所述光阑3_2，限制光斑尺寸，形成近似点光源，然后经过所述准直透镜3_3准直形成平行光。所述收集透镜3_1，所述准直透镜3_3通过透镜或凹面镜实现；所述光阑3_2通过狭缝或小孔实现。
[0077] 所述光谱分光部件4包括色散分光部件和预滤光部件。色散分光部件用于将不同波长的光分开。平行光照在色散分光部件后，不同波长的光会以不同角度透射或反射。色散分光部件采用光栅、棱镜等具有分光能力的器件实现，本实施例中，色散分光部件采用闪耀光栅实现。预滤光部件用于在光照射在光谱分光器件前首先滤出需要探测的波长的光，滤除其他不进行光谱成像的波长的光，可以减小光路中的噪声。预滤光部件由滤光片实现。作为一种可选的实现方式，所述光谱分光部件4仅包括色散分光部件，不包含预滤光部件。
[0078] 所述光谱会聚部件5用于会聚所述光谱分光部件4色散后的光。以相同方向入射到所述光谱会聚部件5的光会聚到其焦平面上相同的点，不同方向入射的光会聚到所述光谱会聚部件5焦平面上不同的点，因此所述光谱会聚部件5将不同波长的光按波长从小到大依次排列在焦平面上。所述光谱会聚部件5由透镜或凹面镜实现。本实施例中，所述光谱会聚部件5由透镜实现。
[0079] 所述第二空间光调制器6放置在所述光谱会聚部件5的焦平面上，因此不同波长的光按波长从小到大依次排列在所述第二空间光调制器6上不同像素上。波长的排列方向可以与所述第二空间光调制器6上像素排列方向平行，也可沿对角线的方向，及其他与所述第二空间光调制器6上像素排列方向成任意角度的方向。本实施例中，波长的排列方向与所述第二空间光调制器6上像素排列方向平行。
[0080] 所述第一空间光调制器2稳定于每一次随机调制时间内时，所述第二空间光调制器6进行若干次随机调制。所述第一空间光调制器2处于不同随机调制状态时，所述第二空间光调制器6的随机调制序列可以是相同的或者不同的。本实施例中，所述第一空间光调制器2处于每一次随机调制状态时，所述第二空间光调制器6的随机调制序列都是相同的。例如，空间光调制器2的随机调制是1010，当空间光调制器2处于1状态时，空间光调制器6进行多次调制，如110011；空间光调制器2处于0状态时，空间光调制器6仍然进行调制110011。
[0081] 所述收集部件7由透镜或凹面镜实现。本实施例中，所述收集部件7采用透镜实现。
[0082] 所述单光子点探测器8采用盖革模式雪崩二极管或光电倍增管实现。本实施例中，所述单光子点探测器8采用盖革模式雪崩二极管实现。
[0083] 所述控制模块11实现对各部件的工作控制和触发脉冲控制，确保所述计数器9和所述第一空间光调制器2、所述第二空间光调制器6之间同步工作包括：所述第一空间光调制器2进行若干次随机调制；当所述第一空间光调制器2稳定于每一次随机调制时间内时，所述第二空间光调制器6进行若干次随机调制；所述第二空间光调制器6每进行一次随机调制，所述计数器9分别累积所述单光子点探测器8发出的代表单光子数目的电脉冲个数，直到所述第二空间光调制器6进行下一次随机调制，所述计数器9再将所述第二空间光调制器6稳定于一次随机调制时间内的光子计数传输至数据包存储器12，并将计数清零，开始下一次计数。
[0084] 所述压缩感知模块13利用所述数据包存储器12中的单光子数目以及对应的随机矩阵，并选取适当的稀疏基对不同波长的天文图像进行重建，得到极弱光水平的天文光谱图像。该模块仅需各波长下天文图像和光谱带的少量线性随机投影便可重建出天文光谱图像，并可利用矩阵填充理论弥补天文光谱图像中的信息缺失。其中，所述的稀疏变换是选取合适的Ψ，使得天文图像x可在Ψ框架下稀疏表达。压缩感知模块13采用下列算法中的任意一种实现压缩感知：匹配跟踪算法MP、正交匹配跟踪算法OMP、基跟踪算法BP、贪心重建算法、LASSO、LARS、GPSR、贝叶斯估计算法、magic、IST、TV、StOMP、CoSaMP、LBI、SP、l1_ls、smp算法、SpaRSA算法、TwIST算法、l0重建算法、l1重建算法、l2重建算法。稀疏基采用离散余弦变换基、小波基、傅里叶变换基、梯度基、gabor变换基中的任意一种。当所测天文图像本身具有很好的稀疏性时，可以不通过稀疏基的变化，直接对原始信号进行重建。
[0085] 图3描述了利用单光子点探测器8探测得到的光子数和第一空间光调制器2、第二空间光调制器6上的随机矩阵实现光谱成像的过程。图中的2表示第一空间光调制器，6表示第二空间光调制器，8表示单光子点探测器，9表示计数器。
[0086] 第一空间光调制器2共进行N次随机调制，对应N个随机矩阵n1，n2…nN，在第一空间光调制器2的一次随机调制时间内，第二空间光调制器6进行M次随机调制，对应M个随机矩阵m1，m2…mM。第二空间光调制器6固定于每次随机调制时间内，单光子点探测器8探测到的光子数由计数器9记录。第一空间光调制器2每一次随机调制时间内，第二空间光调制器6的M次随机调制对应M个计数值，得到一条光强曲线I。第一空间光调制器2的N次随机调制对应N条光强曲线I1，I2…IN。
[0087] 光谱图像重建时，首先分别利用1条光强曲线的M个光子计数值与第二空间光调制器6的M个随机矩阵m1，m2…mM进行压缩感知计算，得到第一空间光调制器2固定于一次随机调制时，第二空间光调制器6上的光强分布，由于第二空间光调制器6上每个像素点对应波长从小到大的光信号，即得到了一条光谱曲线，其光谱段个数与第二空间光调制器6上一行的像素数K相等。由N条光强曲线即可得到N条光谱曲线，与第一空间光调制器2上的N次随机调制对应。
[0088] 取N条光谱曲线上某一相同位置的光强，得到表示某一波长λk的光信号的强度值N个，分别对应第一空间光调制器2上的N次随机调制。利用这N个光子计数值和N个随机矩阵n1，n2…nN进行压缩感知计算，即可得到该波长λk的图像。利用N条光谱曲线上各个位置的光强，与第一空间光调制器2上的N个随机矩阵n1，n2…nN进行压缩感知计算，可以得到各个波长的图像，实现对天文目标的光谱成像。
[0089] 以上是对本发明的基于二阶压缩感知的超灵敏光谱成像天文望远镜的结构说明。下面对该光谱成像天文望远镜的工作过程进行描述。
[0090] 本发明的基于二阶压缩感知的超灵敏光谱成像天文望远镜在工作时包括以下步骤：
[0091] 步骤1)光信号获取的步骤：
[0092] 从天体传播而来的单光子水平的光信号由所述天文望远镜镜头1收集，并成像到所述第一空间光调制器2上；所述第一空间光调制器2对成像在其表面的天文图像进行随机调制，以随机概率将图像上不同位置的光反射到所述准直部件3方向；所述准直部件3包括收集透镜3_1，光阑3_2，准直透镜3_3；所述第一空间光调制器2随机反射的光首先由所述收集透镜3_1会聚到所述光阑3_2，限制光斑尺寸，形成近似点光源，然后经过所述准直透镜3_3准直形成平行光，照射在所述光谱分光部件4上；所述光谱分光部件4将不同波长的光向不同方向出射；经过所述光谱会聚部件5后不同波长的光会聚到所述光谱会聚部件5焦平面上不同位置，形成光谱带；所述第二空间光调制器6放置在所述光谱会聚部件5的焦平面位置，对其表面的光谱带进行随机调制，以随机概率将光谱带上不同位置的光反射到所述收集部件7方向；所述收集部件7收集所述第二空间光调制器6反射而来的光信号，由所述电学单元II中单光子点探测器8探测；
[0093] 步骤2)光学调制与单光子探测、计数同步工作的步骤；
[0094] 所述随机数发生器10包括第一随机数发生器10_1，第二随机数发生器10_2，分别产生随机数用于控制所述第一空间光调制器2和所述第二空间光调制器6；所述第一空间光调制器2根据随机数在图像维度上实现对光信号的随机调制，所述第二空间光调制器6根据该随机数在光谱维度上实现对光信号的随机调制；所述单光子点探测器8探测待测极弱光中的单光子，将采集到的单光子信号转换成脉冲形式的电信号后输出；所述计数器9记录所述单光子点探测器8发出的代表单光子数目的电脉冲个数；所述控制模块11对整个超灵敏光谱成像天文望远镜进行控制协调，包括对各部件的工作控制和同步脉冲触发信号发射，确保所述计数器9、所述第一空间光调制器2和所述第二空间光调制器6同步工作；
[0095] 步骤3)单光子数目与随机矩阵预处理的步骤；
[0096] 对第一随机数发生器(10_1)所对应的随机矩阵不做预处理；
[0097] 当第二随机数发生器(10_2)生成二值伯努利分布的散斑时，若第1帧为全1，将计数器(9)上对应某个波长的一个像素上的所有的单光子数目组成一个列向量y，维度为m×1，m为总的测量数，将第二空间光调制器(6)上的随机矩阵记作A，维度为m×n，n为总的信号长度，令第一帧所对应的单光子数目为y1，则使2y-y1作为新的单光子数目，2A-1作为新的随机矩阵；
[0098] 当第二随机数发生器(10_2)生成二值非均匀分布的散斑时，跳过此步骤3)；
[0099] 步骤4)压缩感知光谱图像恢复的步骤；
[0100] 所述计数器9所记录的单光子数目和所述第一随机数发生器10_1、所述第二随机数发生器10_2生成的随机矩阵全部存入所述数据包存储器12中；所述压缩感知模块13首先利用所述数据包存储器12中的单光子数目以及对应的所述第二随机数发生器10_2生成的随机矩阵进行压缩感知算法重建，获得所述第一空间光调制器2各次随机调制下的光谱曲线；然后利用各条光谱曲线中代表某一波长的光子数值与对应的所述第一随机数发生器10_1生成的随机矩阵进行压缩感知算法重建，获得某一波长的天文图像；分别对不同波长的天文图像进行重建，得到极弱光水平的天文光谱图像。
[0101] 作为一种优选实现方式，在另一个实施例中，在步骤1)之前还包括对所述第二空间光调制器6上各个像素对应波长进行标定的操作。在标定时，首先选取几个特定波长的激光器，或者用滤光片从宽谱光源滤出某些特定波长的光，然后分别将特别波长的光从天文望远镜镜头照射进光学系统，通过对所述第二空间光调制器6的调制状态控制，利用所述单光子点探测器8测量所述第二空间光调制器6上各个像素的光子数，光子数分布最大值的像素位置即对应这些特定波长。其他像素位置对应的波长可按照线性分布计算得出。
[0102] 作为一种优选实现方式，在又一个实施例中，在步骤1)之前还包括有减少仪器噪声的操作。仪器噪声来源包含环境噪声、光学噪声、电学噪声等。压缩感知采样中，信息存在于探测值的涨落中，若仪器噪声的涨落淹没了信号的涨落，则压缩感知算法失效；若仪器噪声的涨落小于或远小于信号的涨落，则能较好甚至完美重建图像。因此，减少仪器噪声有助于提高成像质量。减少仪器噪声的方式有多种，如对仪器进行密闭封装，以遮挡外部环境光信号进入光学系统和探测系统；提高光学部件的透过率，提高仪器内部的清洁度，减少光信号的衰减和散射；提高光谱分光部件4的效率；提高单光子点探测器8的探测效率、暗计数等参数；提高仪器稳定性，减少仪器震动对成像分辨率的影响。
[0103] 作为一种优选实现方式，在又一个实施例中，在步骤1)之前还包括有利用主动光学、自适应光学提高光谱成像信噪比的操作。主动光学是一种为消除望远镜的光学系统及
支架受重力、温度、
风力等影响引起的变形而采用的波面校正技术。通过促动器主动改变主镜镜面的形状，可以修正由于重力、温度和风力造成的镜面本身的形变对成像带来的影响，减少由此产生的光学畸变。自适应光学是一种补偿由大气湍流或其他因素造成的成像过程中波前畸变的技术。自适应光学需要首先检测波前扭曲情况，然后通过安装在望远镜焦面后方的携带有促动器的小型可变形镜面对波前实时进行矫正，从而修复大气湍流等因素对光波波前的扭曲。对天文望远镜镜头1根据主动光学或自适应光学的要求进行设计，可以有效提高天文望远镜镜头1的成像质量，进而提高超灵敏光谱成像天文望远镜获得的天文图像质量。
[0104] 最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行
修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的
权利要求范围当中。
