A Regularized Image Restoration Algorithm for Lossy Compression in Astronomy

Astronomical images currently provide large amounts of data. Lossy compression algorithms have recently been developed for high compression ratios. These compression technique introduce distortion in the compressed images and for high compression ratios, a blocking effect appears. We propose a modified compression algorithm based on the hcompress scheme, and we introduce a new decompression method based on the regularization theory The image is restored scale by scale in a multiresolution scheme and the information lost during the compression is recovered by applying a Tikhonov regularization constraint. The experimental results show that the blocking effect is reduced and some measurements made on a simulated image show that the astrometric and the photometric properties of the restored images are improved.     

数字图象压缩方法在天文上应用的研究及实验

本文研究了数字图象压缩的方法、分类和原理以及该技术在天文领域应用的必要性。针对天文观测的特点和需要，经过分析比较，提出了天文图象压缩的可行性方案。通过应用计算机编程、压缩实验得到相应结果。文中对天文数字图象分别进行了以下压缩方法实验：无失真压缩中的哈夫曼编码和位平面编码；限失真压缩中的离散余弦变换编码及混合编码等。其压缩比分别可达到２—４．５（无失真压缩）和１０—３０（限失真压缩）。图像压缩所需时间随使用的计算机而异。文中详细列出各种压缩方法的实验结果。     

CLEAN算法在天文图像空域重建中的应用

将CLEAN算法用于天文图像空域重建方法 :迭代位移叠加法中的消卷积过程 ,由于它同样是在空间域中进行的迭代减法过程 ,避免了必须把图像变换到傅里叶空间频率域中 ,然后用除法消卷积的传统做法 ,因而满足了仅在空间域中重建目标像的需要 ,使迭代位移叠加法更加简捷。用该方法对天文目标进行的观测实验得到了很好的像复原结果。     

Fast compression and reconstruction of astronomical images based on compressed sensing

With the fast increase in the resolution of astronomical images, the question of how to process and transfer such large images has become a key issue in astronomy. We propose a new real-time compression and fast reconstruction algorithm for astronomical images based on compressive sensing techniques. We first reconstruct the original signal with fewer measurements, according to its compressibility. Then,based on the characteristics of astronomical images, we apply Daubechies orthogonal wavelets to obtain a sparse representation. A matrix representing a random Fourier ensemble is used to obtain a sparse representation in a lower dimensional space. For reconstructing the image, we propose a novel minimum total variation with block adaptive sensing to balance the accuracy and computation time. Our experimental results show that the proposed algorithm can efficiently reconstruct colorful astronomical images with high resolution and improve the applicability of compressed sensing.     

具有泊松噪声分布的图象复原方法

本文以ＣＣＤ图象为技术背景，指出ＣＣＤ图象噪声主要将以与信号有关的泊松噪声形式存在。根据图象形成的模型及图象统计特性，提出了最大似然（ＭＩ）法，最大后验概率（ＭＡＰ）法及最大熵（ＭＥ）复原方法。针对ＭＡＰ法提出了具体算法和参数估计方法。并提出应用图象分割法来提高计算机复原速度，节省存储空间。     

上海天文台1．56米望远镜斑点干涉成像实验进展

斑点干涉成像技术是克服大气湍流影响,提高地面大口径望远镜分辨本领的有效途径之一。该技术利用斑点相机拍摄一系列的短曝光像,使得大气湍流冻结,再经过图像处理获得高分辨率重建像。该技术设备简单,易于实现,很快在观测天文学中得到了广泛的应用,尤其是对双星的研究。首先回顾了天文高分辨率重建技术的发展,并介绍了相关研究成果。描述了几种典型的斑点干涉成像处理方法及其优缺点。对图像噪声类型及滤波方法进行了分析。在上海天文台1．56m望远镜上开展了双星斑点干涉观测实验,目标星等4～7mag,双星目标星等差小于2。分别采用斑点干涉术和迭代位移叠加法成功实现了双星目标的高分辨率成像,初步证明了在1．56m望远镜上进行斑点干涉成像实验,能够达到接近望远镜衍射极限的分辨率水平。     

光学系统参数确定带宽约束的天文图像盲反卷积

大气湍流极大限制了地基大口径望远镜观测天文目标图像的空间分辨率.根据最大似然估计原理,提出了用实际光学带宽约束的可有效减小天文观测图像中大气湍流影响的盲反卷积方法,通过共轭梯度优化算法使卷积误差函数趋向最小.建立了望远镜光学系统参数和图像频域带宽的关系,采用变量正性约束、点扩散函数带宽有限约束,提高算法的收敛性.为避免图像处理中有效傅立叶变换频率超出截止频率,要求采集望远镜焦面图像时单个成像单元(如CCD像素单元)应小于四分之一衍射斑直径.算法中未用目标支持域约束,所提出的方法适用于全视场天文图像恢复.用计算机模拟和对实际天文目标双鱼座图像数据的恢复结果验证了所提出方法的有效性.     

Reconstruction of RHESSI Solar Flare Images with a Forward Fitting Method

We describe a forward-fitting method that has been developed to reconstruct hard X-ray images of solar flares from the Ramaty High-Energy Solar Spectroscopic Imager (RHESSI), a Fourier imager with rotation-modulated collimators that was launched on 5 February 2002. The forward-fitting method is based on geometric models that represent a spatial map by a superposition of multiple source structures, which are quantified by circular gaussians (4 parameters per source), elliptical gaussians (6 parameters), or curved ellipticals (7 parameters), designed to characterize real solar flare hard X-ray maps with a minimum number of geometric elements. We describe and demonstrate the use of the forward-fitting algorithm. We perform some 500 simulations of rotation-modulated time profiles of the 9 RHESSI detectors, based on single and multiple source structures, and perform their image reconstruction. We quantify the fidelity of the image reconstruction, as function of photon statistics, and the accuracy of retrieved source positions, widths, and fluxes. We outline applications for which the forward-fitting code is most suitable, such as measurements of the energy-dependent altitude of energy loss near the limb, or footpoint separation during flares. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1022469811115     

Synthetic Aperture Technic in Astronomy Using Slit Aperture Telescope

The interest in a Rotating Slit-Aperture Telescope (RSAT) among other synthetic aperture telescopes is its capability of being easily coupled with a spectrograph, in order to give reconstructed images of an astronomical object as a function of the light wavelength. Each colored image is comparable with the others for fruitful astrophysical applications. The principle of image reconstruction is well known: it consists of the inversion of the set of projections (Radon transform) given by the telescope during its rotation around its optical axis. A full coverage of the two dimensional Fourier plane can be obtained by rotating the SAT. This problem has led to intense developments for medical imaging (tomography). One of the main difficulties in the reconstruction process in space may come from the jitter of the rotation axis of the RSAT. A set of projections uncorrected for this jitter produces very fuzzy reconstructed images. An elegant solution to the necessary phasing between successive projections is proposed which makes use of a small auxilliary telescope, and some numerical simulations are presented. This revised version was published online in July 2006 with corrections to the Cover Date.     

Derivation of planetary topography using multi-image shape-from-shading

In many cases, the derivation of high-resolution digital terrain models (DTMs) from planetary surfaces using conventional digital image matching is a problem. The matching methods need at least one stereo pair of images with sufficient texture. However, many space missions provide only a few stereo images and planetary surfaces often possess insufficient texture.This paper describes a method for the generation of high-resolution DTMs from planetary surfaces, which has the potential to overcome the described problem. The suggested method, developed by our group, is based on shape-from-shading using an arbitrary number of digital optical images, and is termed “multi-image shape-from-shading” (MI-SFS). The paper contains an explanation of the theory of MI-SFS, followed by a presentation of current results, which were obtained using images from NASA's lunar mission Clementine, and constitute the first practical application with our method using extraterrestrial imagery. The lunar surface is reconstructed under the assumption of different kinds of reflectance models (e.g. Lommel-Seeliger and Lambert). The represented results show that the derivation of a high-resolution DTM of real digital planetary images by means of MI-SFS is feasible.     

图象质量的测度

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

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.     

Automated Processing of LASCO Coronal Images: Spurious Point-Source-Filtering and Missing-Blocks Correction

We report on automated procedures for correcting the images of the LASCO coronagraph for i) spurious quasi-point-sources such as the impacts of cosmic rays, stars, and planets, and ii) the absence of signal due to transmission errors or dropouts, which results in blocks of missing information in the images. Correcting for these undesirable artifacts is mandatory for all quantitative works on the solar corona that require data inversion and/or long series of images, for instance. The nonlinear filtering of spike noise or point-like objects is based on mathematical morphology and implements the procedure opening by morphological reconstruction. However, a simple opening filter is applied whenever the fractional area of corrupted pixels exceeds 50 % of the original image. We describe different strategies for reconstructing the missing information blocks. In general, it is possible to implement the method of averaged neighbors using the two images obtained immediately before and after the corrupted image. For the other cases, and in particular when missing blocks overlapped in three images, we developed an original procedure of weighted interpolation along radial profiles from the center of the Sun that intercept the missing block(s). This procedure is also adequate for the saturated images of bright planets (such as Venus) that bleed along the neighboring pixels. Missing blocks in polarized images may generally be reconstructed using the associated unpolarized image of the same format. But in the case of overlapping missing blocks, we implemented our procedure of weighted interpolation. All tests performed on numerous LASCO-C2 images at various periods of solar activity (i.e. varying complexity of the structure of the corona) demonstrate the excellent performance of these new procedures, with results vastly superior to the methods implemented so far in the pipeline-processing of the LASCO images.     

基于CCD星图的空间目标精确光学定位方法研究

通过光学定位方法来确定空间目标的轨道是实现空间目标监测的重要手段之一。介绍一套小型的空间目标光学定位系统,针对该系统获取的CCD星图,提出Top-hat变换和中值滤波相结合的算法去除星图噪声,并提出基于梯度Hough变换的算法完成对CCD星图星点和轨迹的同时精确定位,最后利用星角距的旋转不变性和球面几何原理,实现对空间目标的精确定位。对分辨率为752×580的CCD实拍星图的应用研究结果表明,星图去噪算法和星点、轨迹同时精确定位算法是有效的,且空间目标的定位精度也可达到角秒量级。     

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

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

一种新的基于2维傅里叶谱图像的恒星光谱特征提取方法和深度网络分类应用

天体光谱分类是天文学研究的重要内容之一,其关键是从光谱数据中选择和提取对分类识别最有效的特征构建特征空间.提出一种新的基于2维傅里叶谱图像的特征提取方法,并应用于LAMOST (the Large Sky Area Multi-Object Fiber Spectroscopic Telescope)恒星光谱数据的分类研究中.光谱数据来源于LAMOST Data Release 5(DR5),选取30000条F、 G和K型星光谱数据,利用短时傅里叶变换(Short-Time Fourier Transform, STFT)将1维光谱数据变换成2维傅里叶谱图像,对得到的2维傅里叶谱图像采用深度卷积网络模型进行分类,得到的分类准确率是92.90%.实验结果表明通过对LAMOST恒星光谱数据进行STFT可得到光谱的2维傅里叶谱图像,谱图像构成了新的光谱数据特征和特征空间,新的特征对于光谱数据分类是有效的.此方法是对光谱分类的一种全新尝试,对海量天体光谱的分类和挖掘处理有一定的开创意义.     

Fast computation of 2D transfer functions from adaptive optics data

The use of atmospheric transfer functions is common in image reconstruction techniques such as speckle interferometry to calibrate the Fourier amplitudes of the reconstructed images. Thus, an accurate model is needed to ensure proper photometry in the reconstruction. The situation complicates when adaptive optics (AO) are used during data acquisition. I propose a novel technique to derive two‐dimensional transfer functions from data collected using AO simultaneously with the observations. The technique is capable to compute the relevant transfer functions within a short time for the prevailing atmospheric conditions and AO performance during data acquisition (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

图象复原的最佳准则

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

斑点全息术原理及其计算机模拟实验

介绍一种天文高分辨率像复原的方法——斑点全息术,并给出了人造目标和真实大气的模拟计算结果,最后指出此方法的应用局限及克服局限性的可能途径。     

用CCD照相多幅叠加处理检测弱星的试验

本文利用CCD动态范围大且便于数值处理的优点,观测月晕中(强背景上)的弱源。在背景尚未使CCD片子饱和时,即停止露光,存入计算机,多次露光,然后在计算机内多幅叠加,减去背景,显出弱源。通过观测试验和定量估算,结果表明,曝光重叠法,用10幅图叠加时能提高信噪比2.65倍。实测与理想情况符合较好。