由斑点图重建天文图像的方法

首先介绍了斑点图形成机制以及获取方法,并借助计算机模拟了斑点图的生成。然后重点描述斑点图的图像重建方法,包括频率域和空间域重建方法。在频率域重建方法部分依次详细叙述了用于复原模信息的斑点干涉法和用于复原相位信息的由模复原相位法、Knox-Thompson法、斑点掩模法;空间域重建方法部分介绍了典型的LWH法和迭代位移叠加法;另外,盲消卷积法也为斑点图重建提供了新的思路,文章对此方法也做了介绍;特别地,还对国内在该领域的相关工作做了介绍;关于斑点图重建方法在天文学、天体物理学等领域的应用,文章做了有针对性的介绍;最后一章是斑点图图像重建方法的总结。     

Astronomical Image Co-adding Method based on Wide-field Image Deconvolution

Due to the atmospheric turbulence, the static aberration, tracking and pointing errors of telescopes, the point spread functions (PSFs) in different fields of view are different. Meanwhile, there are different PSFs in the images obtained by different telescopes. The quality of co-adding image is limited by the image with the poorest quality, and finally the resolution and sensitivity of the quad-channel telescope will also be affected. Dividing the image into some regions with the same type of PSF, and deconvolving these regions can improve the quality of the co-adding image. According to this theory, an image restoration algorithm based on the PSF clustering is proposed. Firstly, this paper makes the PSF clustering analysis by using Self-Organizing Maps, and makes the image segmentation based on the result of the PSF clustering analysis, then using the clustered PSFs to make deconvolutions on the sub-images. Then, the restored sub-images after deconvolution are joined together. Finally, by through the image registration and co-adding, the image with a high signal to noise ratio can be obtained. The result shows that the signal to noise ratio of the astronomical images are improved with our method, and the detection capability on faint stars is also improved.     

基于大视场反卷积的天文图像叠加方法

地基望远镜在成像过程中,由于受大气湍流、望远镜静态像差、跟踪误差、指向误差及视场变化的影响,不同视场区域的PSF (Point Spread Function)具有差异;同时,不同望远镜获取的图像PSF也存在差异.将多个望远镜获取的星象直接叠加至相同的区域后,图像质量受像质最差的望远镜限制,最终观测分辨率和灵敏度均会受到影响.通过图像复原,可以提高图像质量,进而提高叠加效果.根据该思路提出了1种基于PSF分区的迭代图像复原方法:该方法首先通过SOM (Self-organizing Maps)对PSF进行聚类分析,利用同类别PSF的平均PSF进行反卷积,再将反卷积结果按PSF聚类结果分割为不同大小的子图,最后将子图进行拼接.图像复原在提高图像质量的同时,降低了PSF不一致性对图像叠加带来的影响.将几个望远镜在同一时刻获取的图像经反卷积处理之后利用图像配准算法进行矫正并叠加,可获得高信噪比图像.对实际望远镜获取的数据处理后的结果表明:图像在进行复原和叠加过程中,星象目标信噪比不断提升,提高了成像系统对暗星的探测能力.     

Near Real-Time Image Reconstruction

In recent years, post-facto image-processing algorithms have been developed to achieve diffraction-limited observations of the solar surface. We present a combination of frame selection, speckle-masking imaging, and parallel computing which provides real-time, diffraction-limited, 256×256 pixel images at a 1-minute cadence. Our approach to achieve diffraction limited observations is complementary to adaptive optics (AO). At the moment, AO is limited by the fact that it corrects wavefront abberations only for a field of view comparable to the isoplanatic patch. This limitation does not apply to speckle-masking imaging. However, speckle-masking imaging relies on short-exposure images which limits its spectroscopic applications. The parallel processing of the data is performed on a Beowulf-class computer which utilizes off-the-shelf, mass-market technologies to provide high computational performance for scientific calculations and applications at low cost. Beowulf computers have a great potential, not only for image reconstruction, but for any kind of complex data reduction. Immediate access to high-level data products and direct visualization of dynamic processes on the Sun are two of the advantages to be gained.     

图象质量的测度

本文论述了图象质量的各种定量测度,包括连续的和离散的形式。着重讨论了基于视觉特性的图象质量测度,它们与主观评价具有很好的相关性;提出了用图象信息量作图象质量测度的方法。这些测度在图象压缩、图文传真中具有重要的应用。图象质量的测度是对图象质量的评价,它对复原算法优劣的评价是一个十分重要的依据。     

光纤星象切分器

本文介绍了为云台CCD—Coude光谱仪系统研制的星象切分器。计算表明,运用切分器可提高光力2个星等以上。文中讲述了原理并给出有关参数的选择。本星象切分器入射端面为直径1mm的圆面,出射端面为0.15的一线阵,共由38根光纤组成。     

图象复原的最佳准则

本文论述了统计复原的各种最佳准则,这些准则曾广泛地应用于线性或非线性图象复原方法中;指出了最大信息法(MI)更适合复原天文图象。在图象复原领域内,引入了三种新的图象复原准则—最大图象信息法(MII)、最小视觉约方差法(MPMSE)和最小检测概率均方差法(MDP—MSE),前者基于模糊集原理,后两者基于人眼视觉特性。着重讨论了基于检测概率的复原准则(MDP—MSE)、最大信息复原准则(MI)。     

差分像运动视宁度优化监测法

大气视宁度是衡量台址大气光学品质的重要指标。差分像运动视宁度监测仪(Differential Image Motion Monitor,DIMM)被广泛应用于国内外天文选址的视宁度测量作业。介绍了一种优良的视宁度测量方法——差分像运动视宁度优化监测(Improved to Differential Image Motion Monitor,I-DIMM)法。首先对I-DIMM的结构设计和视宁度计算方法进行了详细描述;随后通过设置0. 36 m和0. 12 m两种口径望远镜进行视宁度的模拟测量,将I-DIMM测量结果与传统的DIMM测量的结果进行对比,均证明I-DIMM视宁度计算方法比DIMM更为精确;最后对模拟结果进行了分析,证明了I-DIMM相比于DIMM的优势。     

Multi-Scale Gaussian Normalization for Solar Image Processing

Extreme ultra-violet images of the corona contain information over a wide range of spatial scales, and different structures such as active regions, quiet Sun, and filament channels contain information at very different brightness regimes. Processing of these images is important to reveal information, often hidden within the data, without introducing artefacts or bias. It is also important that any process be computationally efficient, particularly given the fine spatial and temporal resolution of Atmospheric Imaging Assembly on the Solar Dynamics Observatory (AIA/SDO), and consideration of future higher resolution observations. A very efficient process is described here, which is based on localised normalising of the data at many different spatial scales. The method reveals information at the finest scales whilst maintaining enough of the larger-scale information to provide context. It also intrinsically flattens noisy regions and can reveal structure in off-limb regions out to the edge of the field of view. We also applied the method successfully to a white-light coronagraph observation.     

星像抖动的时域频谱

本文顾及湍流外尺度的影响导出了星像抖动频谱的普遍表达式，在均匀光路的情况下得到了分析表达式。理论结果与实验数据相当一致。同时讨论了测量星像抖动所需要的频带宽度以及有限带宽所产生的测量误差。     

CCD图像的一维定心方法

在天体测量学中，数字星象的高精度定位方法日益重要。本文利用边缘分布的方法将二维星像转换到一维，比较了高斯拟合法、修正矩法、中值法和寻导法四种一维定心算法。从精度上看，高斯拟合法最佳，修正矩次之，而寻导法最低。但从计算速度上看，高斯拟合法太慢了。因此作为一个满足高精度和高效率的定心方法，修正矩是最佳选择。     

消除CCD图像中宇宙射线的算法的比较

CCD天文图像在采集过程中会受到各种噪声的影响,其中宇宙射线噪声有时会严重影响到图像中的有用信息。研究如何有效识别和剔除宇宙射线噪声对于天文图像的信息提取是非常重要的。针对目前国际上较新的三种消除宇宙射线的算法:Laplacian边缘检测算法(A Laplacian Edge Detection Algorithm),基于直方图的快速算法(A Fast Algorithm BasedOn Histogram)以及万能噪声消除算法(A Universal Noise Removal Algorithm),采用由云南天文台1m望远镜拍摄的CCD图像进行了模拟实验和实际宇宙射线的处理。实验表明Laplacian边缘检测算法能相对准确地探测到恒星和星系图像上的宇宙射线。并对算法的效果和复杂度进行比较和分析。最后探讨了如何准确替代宇宙射线像素点的灰度值。     

Image recovering from the modulated CCD imagers

The CCD imagers have a spatial resolution comparable to the size of a constituent photosensitive cell. The modulation by means of a rotating Ronchi grating helps to exceed this limit. While the grating rotates in front of the CCD imager the counts of each cell are modulated thus producing a time-varying modulation pattern. The small systematic variations of the form of the pattern encode the fragment of the image pertaining to the cell. The images can be recovered from the modulation patterns by means of the least-squares as well as maximum entropy reconstruction technique. Both methods are able to decode the modulation patterns having high signal-to-noise ratio, however, maximum entropy reconstruction seems to be more robust. A gain in spatial resolution may be of the order of magnitude under favourable circumstances.     

LAMOST二维光纤光谱抽谱方法的研究

对LAMOST二维光纤光谱数据处理系统中的抽谱方法进行了阐述,分析了抽谱中关键参数采样点选取对结果的影响,并根据高斯分布函数的特性以及实验情况确定了采样点的选取范围,解决了抽谱结果出现负值的问题.进一步针对大噪声背景下抽谱存在误差较大的问题,提出了基于频域滤波思想的改进抽谱方法.首先利用快速傅立叶变换和低通滤波器将影响光纤实际轮廓的尖锐噪声滤除,然后再进行正常的抽谱.利用LAMOST二维光纤光谱数据处理系统提供的模拟数据进行了实验测试,结果表明改进抽谱方法的可行性与有效性.     

Review of Image Quality Measures for Solar Imaging

Observations of the solar photosphere from the ground encounter significant problems caused by Earth’s turbulent atmosphere. Before image reconstruction techniques can be applied, the frames obtained in the most favorable atmospheric conditions (the so-called lucky frames) have to be carefully selected. However, estimating the quality of images containing complex photospheric structures is not a trivial task, and the standard routines applied in nighttime lucky imaging observations are not applicable. In this paper we evaluate 36 methods dedicated to the assessment of image quality, which were presented in the literature over the past 40 years. We compare their effectiveness on simulated solar observations of both active regions and granulation patches, using reference data obtained by the Solar Optical Telescope on the Hinode satellite. To create images that are affected by a known degree of atmospheric degradation, we employed the random wave vector method, which faithfully models all the seeing characteristics. The results provide useful information about the method performances, depending on the average seeing conditions expressed by the ratio of the telescope’s aperture to the Fried parameter, \(D/r_{0}\). The comparison identifies three methods for consideration by observers: Helmli and Scherer’s mean, the median filter gradient similarity, and the discrete cosine transform energy ratio. While the first method requires less computational effort and can be used effectively in virtually any atmospheric conditions, the second method shows its superiority at good seeing (\(D/r_{0}<4\)). The third method should mainly be considered for the post-processing of strongly blurred images.     

JPEG2000 Image Compression on Solar EUV Images

For future solar missions as well as ground-based telescopes, efficient ways to return and process data have become increasingly important. Solar Orbiter, which is the next ESA/NASA mission to explore the Sun and the heliosphere, is a deep-space mission, which implies a limited telemetry rate that makes efficient onboard data compression a necessity to achieve the mission science goals. Missions like the Solar Dynamics Observatory (SDO) and future ground-based telescopes such as the Daniel K. Inouye Solar Telescope, on the other hand, face the challenge of making petabyte-sized solar data archives accessible to the solar community. New image compression standards address these challenges by implementing efficient and flexible compression algorithms that can be tailored to user requirements. We analyse solar images from the Atmospheric Imaging Assembly (AIA) instrument onboard SDO to study the effect of lossy JPEG2000 (from the Joint Photographic Experts Group 2000) image compression at different bitrates. To assess the quality of compressed images, we use the mean structural similarity (MSSIM) index as well as the widely used peak signal-to-noise ratio (PSNR) as metrics and compare the two in the context of solar EUV images. In addition, we perform tests to validate the scientific use of the lossily compressed images by analysing examples of an on-disc and off-limb coronal-loop oscillation time-series observed by AIA/SDO.     

一种强有力图象数据缓冲区的开发与使用

以往由于内存小,CCD的读出过程需小块读出,随时存盘;而且由于小块读出中间停顿,数据质量差;尽管现在超过1兆的微机已相当普遍,然而普通操作系统的内存管理仍局限于640k。目前CCD器件的尺寸越来越大,使这一矛盾变得越来越尖锐。针对这些问题,本文介绍一种用扩充内存作图象数据缓冲区的图象采集处理系统方案,并给出具体实现步骤。     

Image Quality in High-resolution and High-cadence Solar Imaging

Broad-band imaging and even imaging with a moderate bandpass (about 1 nm) provides a photon-rich environment, where frame selection (lucky imaging) becomes a helpful tool in image restoration, allowing us to perform a cost-benefit analysis on how to design observing sequences for imaging with high spatial resolution in combination with real-time correction provided by an adaptive optics (AO) system. This study presents high-cadence (160 Hz) G-band and blue continuum image sequences obtained with the High-resolution Fast Imager (HiFI) at the 1.5-meter GREGOR solar telescope, where the speckle-masking technique is used to restore images with nearly diffraction-limited resolution. The HiFI employs two synchronized large-format and high-cadence sCMOS detectors. The median filter gradient similarity (MFGS) image-quality metric is applied, among others, to AO-corrected image sequences of a pore and a small sunspot observed on 2017 June 4 and 5. A small region of interest, which was selected for fast-imaging performance, covered these contrast-rich features and their neighborhood, which were part of Active Region NOAA 12661. Modifications of the MFGS algorithm uncover the field- and structure-dependency of this image-quality metric. However, MFGS still remains a good choice for determining image quality without a priori knowledge, which is an important characteristic when classifying the huge number of high-resolution images contained in data archives. In addition, this investigation demonstrates that a fast cadence and millisecond exposure times are still insufficient to reach the coherence time of daytime seeing. Nonetheless, the analysis shows that data acquisition rates exceeding 50 Hz are required to capture a substantial fraction of the best seeing moments, significantly boosting the performance of post-facto image restoration.