首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 500 毫秒
1.
2.
A probabilistic technique for the joint estimation of background and sources with the aim of detecting faint and extended celestial objects is described. Bayesian probability theory is applied to gain insight into the co-existence of background and sources through a probabilistic two-component mixture model, which provides consistent uncertainties of background and sources. A multiresolution analysis is used for revealing faint and extended objects in the frame of the Bayesian mixture model. All the revealed sources are parametrized automatically providing source position, net counts, morphological parameters and their errors.
We demonstrate the capability of our method by applying it to three simulated data sets characterized by different background and source intensities. The results of employing two different prior knowledge on the source signal distribution are shown. The probabilistic method allows for the detection of bright and faint sources independently of their morphology and the kind of background. The results from our analysis of the three simulated data sets are compared with other source detection methods. Additionally, the technique is applied to ROSAT All-Sky Survey data.  相似文献   

3.
We present and test a new method for the reconstruction of cosmological initial conditions from a full-sky galaxy catalogue. This method, called ZTRACE, is based on a self-consistent solution of the growing mode of gravitational instabilities according to the Zel'dovich approximation and higher order in Lagrangian perturbation theory. Given the evolved redshift-space density field, smoothed on some scale, ZTRACE finds, via an iterative procedure, an approximation to the initial density field for any given set of cosmological parameters; real-space densities and peculiar velocities are also reconstructed. The method is tested by applying it to N -body simulations of an Einstein–de Sitter and an open cold dark matter universe. It is shown that errors in the estimate of the density contrast dominate the noise of the reconstruction. As a consequence, the reconstruction of real-space density and peculiar velocity fields using non-linear algorithms is little improved over those based on linear theory. The use of a mass-preserving adaptive smoothing, equivalent to a smoothing in Lagrangian space, allows an unbiased (although noisy) reconstruction of initial conditions, as long as the (linearly extrapolated) density contrast does not exceed unity. The probability distribution function of the initial conditions is recovered to high precision, even for Gaussian smoothing scales of ∼5  h −1 Mpc, except for the tail at δ ≥1. This result is insensitive to the assumptions of the background cosmology.  相似文献   

4.
An aliasing effect brought up by mass assignment onto Fast Fourier Transformation (FFT) grids may bias measurement of the power spectrum of large scale structures. In this paper, based on the Beylkin's unequally spaced FFT technique, we propose a new precise method to extract the true power spectrum of a large discrete data set. We compare the traditional mass assignment schemes with the new method using the Daub6 and the 3rd-order B-spline scaling functions. Our measurement of Poisson samples and samples of N-body simulations shows that the B-spline scaling function is an optimal choice for mass assignment in the sense that (1) it has a compact support in real space and thus yields an efficient algorithm (2) without any extra corrections. The Fourier space behavior of the 3rd-order B-spline scaling function enables it to be able to accurately recover the true power spectrum with errors less than 5% up to k < kN. It is expected that such a method can be applied to higher order statistics in Fourier space and will enable us to have a precision capture of the non-Gaussian features in the large scale structure of the universe.  相似文献   

5.
An algorithm is proposed for denoising the signal induced by cosmic strings in the cosmic microwave background. A Bayesian approach is taken, based on modelling the string signal in the wavelet domain with generalized Gaussian distributions. Good performance of the algorithm is demonstrated by simulated experiments at arcminute resolution under noise conditions including primary and secondary cosmic microwave background anisotropies, as well as instrumental noise.  相似文献   

6.
A strong correlation between the gamma-ray burster peak energy and the peak luminosity of the associated supernova was discovered by Li for four GRBs. Despite the fact that the formal significance level of the correlation is 0.3 per cent, the smallness of the data set requires careful further evaluation of the result. Subject to the assumption that the data are bivariate Gaussian, a 95 per cent confidence interval of  (−0.9972, 0.02)  for the correlation is derived. Using data from the literature, it is shown that the distribution of known peak GRB energies is not Gaussian if X-ray flashes are included in the sample. This leads to a proposed alternative to the bivariate Gaussian model, which entails describing the dependence between the two variables by a Gaussian copula. The copula is still characterized by a correlation coefficient. The Bayesian posterior distribution of the correlation coefficient is evaluated using a Markov chain Monte Carlo method. The mean values of the posterior distributions range from −0.33 to about zero, depending on the specifics of the supernova (SN) peak brightness distribution. The implication is that the existing data favour a modest correlation between the GRB peak energy and the SN peak brightness; confidence intervals are very wide and include zero.  相似文献   

7.
We discuss spherical needlets and their properties. Needlets are a form of spherical wavelets which do not rely on any kind of tangent plane approximation and enjoy good localization properties in both pixel and harmonic space; moreover needlet coefficients are asymptotically uncorrelated at any fixed angular distance, which makes their use in statistical procedures very promising. In view of these properties, we believe needlets may turn out to be especially useful in the analysis of cosmic microwave background (CMB) data on the incomplete sky, as well as of other cosmological observations. As a final advantage, we stress that the implementation of needlets is computationally very convenient and may rely completely on standard data analysis packages such as healp ix.  相似文献   

8.
We present a Bayesian approach to modelling galaxy clusters using multi-frequency pointed observations from telescopes that exploit the Sunyaev–Zel'dovich effect. We use the recently developed multinest technique to explore the high-dimensional parameter spaces and also to calculate the Bayesian evidence. This permits robust parameter estimation as well as model comparison. Tests on simulated Arcminute Microkelvin Imager observations of a cluster, in the presence of primary CMB signal, radio point sources (detected as well as an unresolved background) and receiver noise, show that our algorithm is able to analyse jointly the data from six frequency channels, sample the posterior space of the model and calculate the Bayesian evidence very efficiently on a single processor. We also illustrate the robustness of our detection process by applying it to a field with radio sources and primordial CMB but no cluster, and show that indeed no cluster is identified. The extension of our methodology to the detection and modelling of multiple clusters in multi-frequency SZ survey data will be described in a future work.  相似文献   

9.
Although tracking with a passive optical telescope is a powerful technique for space debris observation, it is limited by its sensitivity to dynamic background noise. Traditionally, in the field of astronomy, static background subtraction based on a median image technique has been used to extract moving space objects prior to the tracking operation, as this is computationally efficient. The main disadvantage of this technique is that it is not robust to variable illumination conditions. In this article,we propose an approach for tracking small and dim space debris in the context of a dynamic background via one of the optical telescopes that is part of the space surveillance network project, named the AsiaPacific ground-based Optical Space Observation System or APOSOS. The approach combines a fuzzy running Gaussian average for robust moving-object extraction with dim-target tracking using a particlefilter-based track-before-detect method. The performance of the proposed algorithm is experimentally evaluated, and the results show that the scheme achieves a satisfactory level of accuracy for space debris tracking.  相似文献   

10.
We adapt the smooth tests of goodness-of-fit developed by Rayner and Best to the study of the non-Gaussianity of interferometric observations of the cosmic microwave background (CMB). The interferometric measurements (visibilities) are transformed into signal-to-noise ratio eigenmodes, and then the method is applied directly in Fourier space. This transformation allows us to perform the analysis in different subsets of eigenmodes according to their signal-to-noise ratio level. The method can also deal with non-uniform or incomplete coverage of the UV plane. We explore here two possibilities: we analyse either the real and imaginary parts of the complex visibilities (Gaussianly distributed under the Gaussianity hypothesis) or their phases (uniformly distributed under the Gaussianity hypothesis). The power of the method in discriminating between Gaussian and non-Gaussian distributions is studied by using several kinds of non-Gaussian simulations. On the one hand, we introduce a certain degree of non-Gaussianity directly into the Fourier space using the Edgeworth expansion, and afterwards the desired correlation is introduced. On the other hand, we consider interferometric observations of a map with topological defects (cosmic strings). To these previous non-Gaussian simulations we add different noise levels and quantify the required signal-to-noise ratio necessary to achieve a detection of these non-Gaussian features. Finally, we have also studied the ability of the method to constrain the so-called non-linear coupling constant f NL using χ2 simulations. The whole method is illustrated here by application to simulated data from the Very Small Array interferometer.  相似文献   

11.
天体光谱分类是天文学研究的重要内容之一,其关键是从光谱数据中选择和提取对分类识别最有效的特征构建特征空间.提出一种新的基于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维傅里叶谱图像,谱图像构成了新的光谱数据特征和特征空间,新的特征对于光谱数据分类是有效的.此方法是对光谱分类的一种全新尝试,对海量天体光谱的分类和挖掘处理有一定的开创意义.  相似文献   

12.
We present a simple and intuitive approximation for solving the perturbation theory (PT) of small cosmic fluctuations. We consider only the spherically symmetric or monopole contribution to the PT integrals, which yields the exact result for tree-graphs (i.e. at leading order). We find that the non-linear evolution in Lagrangian space is then given by a simple local transformation over the initial conditions, although it is not local in Euler space. This transformation is found to be described by the spherical collapse (SC) dynamics, as it is the exact solution in the shearless (and therefore local) approximation in Lagrangian space. Taking advantage of this property, it is straightforward to derive the one-point cumulants, ξJ, for both the unsmoothed and smoothed density fields to arbitrary order in the perturbative regime. To leading-order this reproduces, and provides us with a simple explanation for, the exact results obtained by Bernardeau. We then show that the SC model leads to accurate estimates for the next corrective terms when compared with the results derived in the exact perturbation theory making use of the loop calculations. The agreement is within a few per cent for the hierarchical ratios S J  = ξ J J −12. We compare our analytic results with N -body simulations, which turn out to be in very good agreement up to scales where σ ≈ 1. A similar treatment is presented to estimate higher order corrections in the Zel'dovich approximation. These results represent a powerful and readily usable tool to produce analytical predictions that describe the gravitational clustering of large-scale structure in the weakly non-linear regime.  相似文献   

13.
A model of the gravitationally evolved dark matter distribution, in the Eulerian space, is developed. It is a simple extension of the excursion set model that is commonly used to estimate the mass function of collapsed dark matter haloes. In addition to describing the evolution of the Eulerian space distribution of the haloes, the model allows one to describe the evolution of the dark matter itself. It can also be used to describe density profiles, on scales larger than the virial radius of these haloes, and to quantify the way in which matter flows in and out of Eulerian cells. When the initial Lagrangian space distribution is white noise Gaussian, the model suggests that the Inverse Gaussian distribution should provide a reasonably good approximation to the evolved Eulerian density field, in agreement with numerical simulations. Application of this model to clustering from more general Gaussian initial conditions is discussed at the end.  相似文献   

14.
With increasingly large data sets, weak lensing measurements are able to measure cosmological parameters with ever-greater precision. However, this increased accuracy also places greater demands on the statistical tools used to extract the available information. To date, the majority of lensing analyses use the two-point statistics of the cosmic shear field. These can be either studied directly using the two-point correlation function or in Fourier space, using the power spectrum. But analysing weak lensing data inevitably involves the masking out of regions, for example to remove bright stars from the field. Masking out the stars is common practice but the gaps in the data need proper handling. In this paper, we show how an inpainting technique allows us to properly fill in these gaps with only   N log  N   operations, leading to a new image from which we can compute straightforwardly and with a very good accuracy both the power spectrum and the bispectrum. We then propose a new method to compute the bispectrum with a polar fft algorithm, which has the main advantage of avoiding any interpolation in the Fourier domain. Finally, we propose a new method for dark matter mass map reconstruction from shear observations, which integrates this new inpainting concept. A range of examples based on 3D N -body simulations illustrates the results.  相似文献   

15.
16.
The maximum-entropy method is applied to the problem of reconstructing the projected mass density of a galaxy cluster using its gravitational lensing effects on background galaxies. We demonstrate the method by reconstructing the mass distribution in a model cluster using simulated shear and magnification data to which Gaussian noise is added. The mass distribution is reconstructed directly and the inversion is regularized using the entropic prior for this positive additive distribution. For realistic noise levels, we find that the method faithfully reproduces the main features of the cluster mass distribution not only within the observed field but also slightly beyond it. We estimate the uncertainties in the reconstruction by calculating an analytic approximation to the covariance matrix of the reconstruction values of each pixel. This result is compared with error estimates derived from Monte Carlo simulations for different noise realizations and found to be in good agreement.  相似文献   

17.
Finding outlier light curves in catalogues of periodic variable stars   总被引:2,自引:0,他引:2  
We present a methodology to discover outliers in catalogues of periodic light curves. We use a cross-correlation as the measure of 'similarity' between two individual light curves, and then classify light curves with lowest average 'similarity' as outliers. We performed the analysis on catalogues of periodic variable stars of known type from the MACHO and OGLE projects. This analysis was carried out in Fourier space and we established that our method correctly identifies light curves that do not belong to those catalogues as outliers. We show how an approximation to this method, carried out in real space, can scale to large data sets that will be available in the near future such as those anticipated from the Panoramic Survey Telescope & Rapid Response System (Pan-STARRS) and Large Synoptic Survey Telescope (LSST).  相似文献   

18.
The Fourier transform of cosmological density perturbations can be represented in terms of amplitudes and phases for each Fourier mode. We investigate the phase evolution of these modes using a mixture of analytical and numerical techniques. Using a toy model of one-dimensional perturbations evolving under the Zel'dovich approximation as an initial motivation, we develop a statistic that quantifies the information content of the distribution of phases. Using numerical simulations beginning with more realistic Gaussian random-phase initial conditions, we show that the information content of the phases grows from zero in the initial conditions, first slowly and then rapidly when structures become non-linear. This growth of phase information can be expressed in terms of an effective entropy. Gaussian initial conditions are a maximum entropy realization of the initial power spectrum; gravitational evolution decreases the phase entropy. We show that our definition of phase entropy results in a statistic that explicitly quantifies the information stored in the phases of density perturbations (rather than their amplitudes), and that this statistic displays interesting scaling behaviour for self-similar initial conditions.  相似文献   

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

20.
The Voigt–Hjerting function is fundamental in order to correctly model the profiles of absorption lines imprinted on the spectra of bright background sources by intervening absorbing systems. In this work, we present a simple analytic approximation to this function in the context of absorption-line profiles of intergalactic H  i absorbers. Using basic calculus tools, we derive an analytic expression for the Voigt–Hjerting function that contains only fourth-order polynomial and Gaussian functions. In connection with the absorption coefficient of intergalactic neutral hydrogen, this approximation is suitable for modelling Voigt profiles with an accuracy of 10−4 or better for an arbitrary wavelength baseline, for column densities up to   N H I= 10 22 cm−2  , and for damping parameters   a ≲ 10−4  , that is, the entire range of parameters characteristic to all Lyman transitions arising in a variety of H  i absorbing systems such as Lyman α (Lyα) forest clouds, Lyman limit systems and damped Lyα systems. We hence present an approximation to the Voigt–Hjerting function that is both accurate and flexible to implement in various types of programming languages and machines, and with which Voigt profiles can be calculated in a reliable and very simple manner.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号