基于EMCCD的成像系统的构建

为了探索对较暗目标进行成像的新方法，我们在TI公司生产的带有增益寄存器的EMCCD（Electron Multiply CCD）芯片基础上，构建了一套成像系统，主要用于对此芯片进行性能测试，为将来构建天文观测用相机做准备。本文主要介绍了EMCCD的结构特点、工作原理，并且详细介绍实验室采用TI公司的EMCCD搭建的成像系统，以及在此系统基础上对EMCCD进行地评测。     

Direct imaging with highly diluted apertures – II. Properties of the point spread function of a hypertelescope

In the future, optical stellar interferometers will provide true images thanks to larger number of telescopes and to advanced cophasing subsystems. These conditions are required to have sufficient resolution elements (resel) in the image and to provide direct images in the hypertelescope mode. It has already been shown that hypertelescopes provide snapshot images with a significant gain in sensitivity without inducing any loss of the useful field of view for direct imaging applications. This paper aims at studying the properties of the point spread functions of future large arrays using the hypertelescope mode. Numerical simulations have been performed and criteria have been defined to study the image properties. It is shown that the choice of the configuration of the array is a trade-off between the resolution, the halo level and the field of view. A regular pattern of the array of telescopes optimizes the image quality (low halo level and maximum encircled energy in the central peak), but decreases the useful field of view. Moreover, a non-redundant array is less sensitive to the space aliasing effect than a redundant array.     

Particle-mesh simulations of the Lyα forest

We have analysed short-exposure high-resolution images obtained on the William Herschel Telescope (WHT), Big Telescope Azimuthal (BTA) and CanadaFranceHawaii Telescope (CFHT) in order to investigate mechanical vibrations that are capable of reducing the high angular resolution allowed by speckle-interferometric methods. After filtering the photon noise we have computed power spectra of the image centres of gravity. In these spectra we have found vibration features that differ from one telescope to another and that vary in particular with the zenith angle. In the case of the WHT we discuss these results in terms of possible causes and in terms of degradation of the transfer function. We present some means to improve the data.     

基于USB30.0的EMCCD相机高速数据传输系统韵研究

自动成像协会（AutomatedImagingAssociation,AIA）在2013年初发布了用于高速图像数据传输的USB3Vision标准。首先简要介绍了基于该标准的EMCCD相机高速数据传输系统的设计方案,重点介绍了传输系统的构建和图像采集软件的设计。其中,传输系统的构建主要是在QUARTUSII的开发环境下,移植EMCCD相机数字控制器,使用VHDL语言编程设计一个控制器,产生USB3．0芯片USB3014的读写时序以及相关的逻辑信号,并且完成模拟图像产生和针对USB3Vision标准的数据传输格式转换的功能;图像采集软件的编写是在VS2010开发环境下利用CYPRESS公司提供的应用程序接151（API）,采用c＋＋语言实现。最后进行了模拟图像的采集实验,并且进行了误码率估算。     

Photon Spectroscopy with Imaging X-Ray Instruments

Individual X-ray photons in the keV energy range produce hundreds of photoelectrons in a single pixel of a CCD array detector. The number of photoelectrons produced is a linear function of the photon energy, allowing the measurement of spectral information with an imaging detector system. Most solar X-ray telescopes, such as Yohkoh/SXT and Hinode/XRT, use CCD detectors in an integrating mode and are designed to make temperature estimates from multiband filter photometry. We show how such instruments can be used in a new way to perform a limited type of this photon spectroscopy. By measuring the variance in intensity of a series of repeated images through a single filter of an X-ray source, the mean energy per detected photon can be determined. This energy is related to the underlying coronal spectrum, and hence it can be used to deduce the mean plasma temperature. We apply this technique to data from the Yohkoh Soft X-Ray Telescope and compare the temperatures obtained with this technique with the temperatures derived using the standard filter ratio method for a postflare loop system. Given the large dynamic range of the soft X-ray flux observed from the Sun, we describe the requirements for a future instrument that would be better suited to performing photon spectroscopy. B.J. Labonte deceased 24 October 2005.     

Digital stacking of photographic plates with SuperCOSMOS

Photographic Schmidt plates are among the most effective tools in wide-field astronomy. One of the principal difficulties of photographic plates when compared to modern detectors is lack of image depth. We present a technique for stacking plates digitized using the SuperCOSMOS microdensitometer, aimed at maximizing the signal-to-noise ratio in faint images. The efficacy of several image combination algorithms is tested by stacking plate frames in the presence of spurious images. We find that an 'average sigma clipping' type pixel rejection in conjunction with our weighting scheme is most effective in delivering a clean, high signal-to-noise ratio stack. The gain in limiting magnitude obtained from stacking is found to be consistent with that expected: Δ M ∼1.5 for a stack of 16 good-quality plates.     

Wide-field tracking with zenith-pointing telescopes

Equipped with a suitable optical relay system, telescopes employing low-cost fixed primary mirrors could point and track while delivering high-quality images to a fixed location. Such an optical tracking system would enable liquid-mirror telescopes to access a large area of sky and employ infrared detectors and adaptive optics. Such telescopes could also form the elements of an array in which light is combined either incoherently or interferometrically. Tracking of an extended field requires correction of all aberrations including distortion, field curvature and tilt. A specific design is developed that allows a 10-m liquid-mirror telescope to track objects for as long as 30 min and to point as far as 4° from the zenith, delivering a distortion-free diffraction-limited image to a stationary detector, spectrograph or interferometric beam combiner.     

Detection of Faint Asteroids Based on Image Shifting and Stacking Method

In order to improve the ability to find the faint and small celestial bodies in the solar system, a method of shifting and stacking images which improves the detection efficiency of faint moving objects is applied to process the sequential optical images. This method determines the existence of moving objects by using the method of false position to pre-estimate the apparent velocities of moving objects, then determines iteratively the accurate positions of moving objects based on the SNR (Signal-to-Noise Ratio) and elongation of stellar image. Using the sequential images of the China Near Earth Object Survey Telescope (CNEOST), we carry out an experiment and succeed in detecting asteroids fainter than 21 magnitude which are invisible on a single image. Thus, the feasibility of this method is verified.     

Photon counting strategies with low-light-level CCDs

Low light level charge-coupled devices (L3CCDs) have recently been developed, incorporating on-chip gain. They may be operated to give an effective readout noise of much less than one electron by implementing an on-chip gain process allowing the detection of individual photons. However, the gain mechanism is stochastic and so introduces significant extra noise into the system. In this paper we examine how best to process the output signal from an L3CCD so as to minimize the contribution of stochastic noise, while still maintaining photometric accuracy.
We achieve this by optimizing a transfer function that translates the digitized output signal levels from the L3CCD into a value approximating the photon input as closely as possible by applying thresholding techniques. We identify several thresholding strategies and quantify their impact on the photon counting accuracy and the effective signal-to-noise ratio.
We find that it is possible to eliminate the noise introduced by the gain process at the lowest light levels. Reduced improvements are achieved as the light level increases up to about 20 photon pixel⁻¹ and above this there is negligible improvement. Operating L3CCDs at very high speeds will keep the photon flux low, giving the best improvements in signal-to-noise ratio.     

Instrument and data analysis challenges for imaging spectropolarimetry

The next generation of solar telescopes will enable us to resolve the fundamental scales of the solar atmosphere, i.e., the pressure scale height and the photon mean free path. High‐resolution observations of small‐scale structures with sizes down to 50 km require complex post‐focus instruments, which employ adaptive optics (AO) and benefit from advanced image restoration techniques. The GREGOR Fabry‐Pérot Interferometer (GFPI) will serve as an example of such an instrument to illustrate the challenges that are to be expected in instrumentation and data analysis with the next generation of solar telescopes (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

明安图射电频谱日像仪图像的强度定标方法

明安图射电频谱日像仪(Mingantu Spectral Radioheliograph, MUSER)能够在0.4--15GHz超宽频带内实现高时间、高空间、高频率分辨率的太阳射电成像. 而射电亮温度是描述太阳物理过程的一个重要的参数, 在研究不同射电辐射机制、太阳磁场以及太阳爆发过程中非热粒子加速等问题上有着非常重要的作用, 因此必须对MUSER观测的图像进行亮温度定标. 将介绍一种适用于射电日像仪图像强度定标的方法. 太阳射电图像中包含着太阳圆盘的结构信息, 利用射电日像仪短基线的可视度函数拟合第一类贝塞尔函数, 可以得到图像中宁静太阳圆盘的射电半径和强度, 再利用瑞利-金斯定律和每天的太阳射电流量可以计算得到每天图像的定标因子$G_c$, 从而实现对MUSER图像强度的定标. 将该方法应用到MUSER的实际观测数据中, 包括宁静太阳和有太阳射电爆发等不同的情况, $G_c$的误差基本不超过10%, 得到的宁静太阳亮温度与其他宁静太阳的结果具有较高的相关性, 表明了此方法的可行性和有效性.     

High‐order adaptive optical system for Big Bear Solar Observatory

A high‐order Adaptive Optical (AO) system for the 65 cm vacuum telescope of the Big Bear Solar Observatory (BBSO) is presented. The Coudé‐exit of the telescope has been modified to accommodate the AO system and two imaging magnetograph systems for visible‐light and near infrared (NIR) observations. A small elliptical tip/tilt mirror directs the light into an optical laboratory on the observatory's 2^nd floor just below the observing floor. A deformable mirror (DM) with 77 mm diameter is located on an optical table where it serves two wave‐front sensors (WFS), a correlation tracker (CT) and Shack‐Hartman (SH) sensor for the high‐order AO system, and the scientific channels with the imaging magnetographs. The two‐axis tip/tilt platform has a resonance frequency around 3.3 kHz and tilt range of about 2 mrad, which corresponds to about 25″ in the sky. Based on 32 × 32 pixel images, the CT detects image displacements between a reference frame and real‐time frames at a rate of 2 kHz. High‐order wave‐front aberrations are detected in the SH WFS channel from slope measurements derived from 76 sub‐apertures, which are recorded with 1,280 × 1,024 pixel Complex Metal Oxide Semiconductor (CMOS) camera manufactured by Photobit camera. In the 4 × 4 pixel binning mode, the data acquisition rate of the CMOS device is more than 2 kHz. Both visible‐light and NIR imaging magnetographs use Fabry‐Pérot etalons in telecentric configurations for two‐dimensional spectro‐polarimetry. The optical design of the AO system allows using small aperture prefilters, such as interference or Lyot filters, and 70 mm diameter Fabry‐Pérot etalons covering a field‐of‐view (FOV) of about 180″ × 180″.     

Transmission profile of the Dutch Open Telescope Hα Lyot filter

Accurate knowledge of the spectral transmission profile of a Lyot filter is important, in particular in comparing observations with simulated data. The paper summarizes available facts about the transmission profile of the Dutch Open Telescope (DOT) Hα Lyot filter pointing to a discrepancy between sidelobe‐free Gaussian‐like profile measured spectroscopically and signatures of possible leakage of parasitic continuum light in DOT Hα images. We compute wing‐to‐center intensity ratios resulting from convolutions of Gaussian and square of the sinc function with the Hα atlas profile and compare them with the ratios derived from observations of the quiet Sun chromosphere at disk center. We interpret discrepancies between the anticipated and observed ratios and the sharp limb visible in the DOT Hα image as an indication of possible leakage of parasitic continuum light. A method suggested here can be applied also to indirect testing of transmission profiles of other Lyot filters. We suggest two theoretical transmission profiles of the DOT Hα Lyot filter which should be considered as the best available approximations. Conclusive answer can only be given by spectroscopic re‐measurement of the filter. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

From planes to spheres: about gravitational lens magnifications

We discuss the classic theorem according to which a gravitational lens always produces at least one image with a magnification greater than unity. This theorem seems to contradict the conservation of total flux from a lensed source. The standard solution to this paradox is based on the exact definition of the reference 'unlensed' situation, in which the lens mass can be either removed or smoothly redistributed.
We calculate magnifications and amplifications (in photon number and energy flux density) for general lensing scenarios not limited to regions close to the optical axis. In this way the formalism is naturally extended from tangential planes for the source and lensed images to complete spheres. We derive the lensing potential theory on the sphere and find that the Poisson equation is modified by an additional source term that is related to the mean density and to the Newtonian potential at the positions of observer and source. This new term generally reduces the magnification, to below unity far from the optical axis, and ensures conservation of the total photon number received on a sphere around the source.
This discussion does not affect the validity of the focusing theorem , in which the unlensed situation is defined to have an unchanged affine distance between source and observer. The focusing theorem does not contradict flux conservation, because the mean total magnification (or amplification) directly corresponds to different areas of the source (or observer) sphere in the lensed and unlensed situation. We argue that a constant affine distance does not define an astronomically meaningful reference.
By exchanging source and observer, we confirm that magnification and amplification differ according to Etherington's reciprocity law, so that surface brightness is no longer strictly conserved. At this level we also have to distinguish between different surface brightness definitions that are based on photon number, photon flux and energy flux.     

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)     

Estimating statistics of sky brightness using radio interferometric observations

Radio interferometers are used to construct high resolution images of the sky at radio frequencies and are the key instruments for accessing the statistical properties of the evolution of neutral hydrogen over cosmic time. Here we use simulated observations of the model sky to assess the efficacy of different estimators of the large-scale structure and power spectrum of the sky brightness distribution. We find that while the large-scale distribution can be reasonably estimated using the reconstructed image from interferometric data, estimates of the power spectrum of the intensity fluctuations calculated from the image are generally biased. This bias is found to be more pronounced for diffuse emission. The visibility based power spectrum estimator, however, gives an unbiased estimate of the true power spectrum. This work demonstrates that for an observation with diffuse emission the reconstructed image can be used to estimate the large-scale distribution of the intensity, while to estimate the power spectrum, visibility based methods should be preferred.With the upcoming experiments aimed at measuring the evolution of the power spectrum of the neutral hydrogen distribution, this is a very important result.     

Adaptive binning of X-ray galaxy cluster images

We present a simple method for adaptively binning the pixels in an image. The algorithm groups pixels into bins of size such that the fractional error on the photon count in a bin is less than or equal to a threshold value, and the size of the bin is as small as possible. The process is particularly useful for generating surface brightness and colour maps, with clearly defined error maps, from images with a large dynamic range of counts, for example X-ray images of galaxy clusters. We demonstrate the method in application to data from Chandra ACIS-S and ACIS-I observations of the Perseus cluster of galaxies. We use the algorithm to create intensity maps, and colour images that show the relative X-ray intensities in different bands. The colour maps can later be converted, through spectral models, into maps of physical parameters, such as temperature, column density, etc. The adaptive binning algorithm is applicable to a wide range of data, from observations or numerical simulations, and is not limited to two-dimensional data.     

Comparing modal noise and FRD of circular and non‐circular cross‐section fibres

Modal noise is a common source of noise introduced to the measurements by optical fibres and is particularly important for fibre‐fed spectroscopic instruments, especially for high‐resolution measurements. This noise source can limit the signal‐to‐noise ratio and jeopardize photon‐noise limited data. The subject of the present work is to compare measurements of modal noise and focal‐ratio degradation (FRD) for several commonly used fibres. We study the influence of a simple mechanical scrambling method (excenter) on both FRD and modal noise. Measurements are performed with circular and octagonal fibres from Polymicro Technology (FBP‐Series) with diameters of 100, 200, and 300μm and for square and rectangular fibres from CeramOptec, among others. FRD measurements for the same sample of fibres are performed as a function of wavelength. Furthermore, we replaced the circular fibre of the STELLA‐échelle‐spectrograph (SES) in Tenerife with an octagonal and found a SNR increase by a factor of 1.6 at 678 nm. It is shown in the laboratory that an excenter with a large amplitude and low frequency will not influence the FRD but will reduce modal noise rather effectively by up to 180%. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Exact optics – VI. Schmidt cameras and prime correctors

The methods of exact optics are here extended to cover the design of lens and mirror systems. Power series are given for the shapes of both Schmidt corrector plates and their complementary almost spherical mirrors. Their off-axis images are analysed.
A meniscus glass corrector extends the prime focus field of an f /3.333 reflector with subarcsecond images by a factor greater than 10.     

TRIFFID observations of the cores of the three globular clusters M15, M92 and NGC 6712

We present the results of TRIFFID simultaneous V - and B -band observations of the cores of the globular clusters M15, M92 and NGC 6712. A variability search of their dense centres was made feasible through performing post-exposure image sharpening on the images, increasing the image resolution by a factor of ∼2. The isis implementation of the image subtraction technique developed by Alard & Lupton was then used to detect flux variations in our image sets. We have obtained periods for all observable variables (in our field of view) in NGC 6712 and we have found two new RR Lyrae variables (an RRab and an RRc). We have confirmed three variables in our field of view of the M92. For M15, we detect 48 variables in our field of view, 23 of which are new discoveries. We obtain periods and amplitudes for all variables and classify new ones based on the light-curve shape, the most significant period and the mean magnitude in the V band. Among the detected RR Lyrae we find 19 RRc, 12 RRab and two RRd types. In the subsequent analysis we find a marked increase in RRc over RRab variables in the core. In a refined procedure to search for fainter objects we find no dwarf novae in our field of view of M15. Simulations performed on the data set to quantify our sensitivity to such objects indicate that an upper limit of 10 dwarf novae (at 92 per cent probability) exist in our field of view. The implications this result has on globular clusters are discussed.