1,989 research outputs found
Measuring the galaxy power spectrum and scale-scale correlations with multiresolution-decomposed covariance -- I. method
We present a method of measuring galaxy power spectrum based on the
multiresolution analysis of the discrete wavelet transformation (DWT). Since
the DWT representation has strong capability of suppressing the off-diagonal
components of the covariance for selfsimilar clustering, the DWT covariance for
popular models of the cold dark matter cosmogony generally is diagonal, or
(scale)-diagonal in the scale range, in which the second scale-scale
correlations are weak. In this range, the DWT covariance gives a lossless
estimation of the power spectrum, which is equal to the corresponding Fourier
power spectrum banded with a logarithmical scaling. In the scale range, in
which the scale-scale correlation is significant, the accuracy of a power
spectrum detection depends on the scale-scale or band-band correlations. This
is, for a precision measurements of the power spectrum, a measurement of the
scale-scale or band-band correlations is needed. We show that the DWT
covariance can be employed to measuring both the band-power spectrum and second
order scale-scale correlation. We also present the DWT algorithm of the binning
and Poisson sampling with real observational data. We show that the alias
effect appeared in usual binning schemes can exactly be eliminated by the DWT
binning. Since Poisson process possesses diagonal covariance in the DWT
representation, the Poisson sampling and selection effects on the power
spectrum and second order scale-scale correlation detection are suppressed into
minimum. Moreover, the effect of the non-Gaussian features of the Poisson
sampling can be calculated in this frame.Comment: AAS Latex file, 44 pages, accepted for publication in Ap
A Multichannel Spatial Compressed Sensing Approach for Direction of Arrival Estimation
The final publication is available at http://link.springer.com/chapter/10.1007%2F978-3-642-15995-4_57ESPRC Leadership Fellowship EP/G007144/1EPSRC Platform Grant EP/045235/1EU FET-Open Project FP7-ICT-225913\"SMALL
Quasi-local evolution of cosmic gravitational clustering in the weakly non-linear regime
We investigate the weakly non-linear evolution of cosmic gravitational
clustering in phase space by looking at the Zel'dovich solution in the discrete
wavelet transform (DWT) representation. We show that if the initial
perturbations are Gaussian, the relation between the evolved DWT mode and the
initial perturbations in the weakly non-linear regime is quasi-local. That is,
the evolved density perturbations are mainly determined by the initial
perturbations localized in the same spatial range. Furthermore, we show that
the evolved mode is monotonically related to the initial perturbed mode. Thus
large (small) perturbed modes statistically correspond to the large (small)
initial perturbed modes. We test this prediction by using QSO Ly
absorption samples. The results show that the weakly non-linear features for
both the transmitted flux and identified forest lines are quasi-localized. The
locality and monotonic properties provide a solid basis for a DWT
scale-by-scale Gaussianization reconstruction algorithm proposed by Feng & Fang
(Feng & Fang, 2000) for data in the weakly non-linear regime. With the
Zel'dovich solution, we find also that the major non-Gaussianity caused by the
weakly non-linear evolution is local scale-scale correlations. Therefore, to
have a precise recovery of the initial Gaussian mass field, it is essential to
remove the scale-scale correlations.Comment: 22 pages, 13 figures. Accepted for publication in the Astrophysical
Journa
A survey of parallel algorithms for fractal image compression
This paper presents a short survey of the key research work that has been undertaken in the application of parallel algorithms for Fractal image compression. The interest in fractal image compression techniques stems from their ability to achieve high compression ratios whilst maintaining a very high quality in the reconstructed image. The main drawback of this compression method is the very high computational cost that is associated with the encoding phase. Consequently, there has been significant interest in exploiting parallel computing architectures in order to speed up this phase, whilst still maintaining the advantageous features of the approach. This paper presents a brief introduction to fractal image compression, including the iterated function system theory upon
which it is based, and then reviews the different techniques that have been, and can be, applied in order to parallelize the compression algorithm
On the efficient Monte Carlo implementation of path integrals
We demonstrate that the Levy-Ciesielski implementation of Lie-Trotter
products enjoys several properties that make it extremely suitable for
path-integral Monte Carlo simulations: fast computation of paths, fast Monte
Carlo sampling, and the ability to use different numbers of time slices for the
different degrees of freedom, commensurate with the quantum effects. It is
demonstrated that a Monte Carlo simulation for which particles or small groups
of variables are updated in a sequential fashion has a statistical efficiency
that is always comparable to or better than that of an all-particle or
all-variable update sampler. The sequential sampler results in significant
computational savings if updating a variable costs only a fraction of the cost
for updating all variables simultaneously or if the variables are independent.
In the Levy-Ciesielski representation, the path variables are grouped in a
small number of layers, with the variables from the same layer being
statistically independent. The superior performance of the fast sampling
algorithm is shown to be a consequence of these observations. Both mathematical
arguments and numerical simulations are employed in order to quantify the
computational advantages of the sequential sampler, the Levy-Ciesielski
implementation of path integrals, and the fast sampling algorithm.Comment: 14 pages, 3 figures; submitted to Phys. Rev.
Scalar and vector modulation instabilities induced by vacuum fluctuations in fibers: numerical study
We study scalar and vector modulation instabilities induced by the vacuum
fluctuations in birefringent optical fibers. To this end, stochastic coupled
nonlinear Schrodinger equations are derived. The stochastic model is equivalent
to the quantum field operators equations and allow for dispersion,
nonlinearity, and arbitrary level of birefringence. Numerical integration of
the stochastic equations is compared to analytical formulas in the case of
scalar modulation instability and non depleted pump approximation. The effect
of classical noise and its competition with vacuum fluctuations for inducing
modulation instability is also addressed.Comment: 33 pages, 5 figure
A Multiresolution Census Algorithm for Calculating Vortex Statistics in Turbulent Flows
The fundamental equations that model turbulent flow do not provide much
insight into the size and shape of observed turbulent structures. We
investigate the efficient and accurate representation of structures in
two-dimensional turbulence by applying statistical models directly to the
simulated vorticity field. Rather than extract the coherent portion of the
image from the background variation, as in the classical signal-plus-noise
model, we present a model for individual vortices using the non-decimated
discrete wavelet transform. A template image, supplied by the user, provides
the features to be extracted from the vorticity field. By transforming the
vortex template into the wavelet domain, specific characteristics present in
the template, such as size and symmetry, are broken down into components
associated with spatial frequencies. Multivariate multiple linear regression is
used to fit the vortex template to the vorticity field in the wavelet domain.
Since all levels of the template decomposition may be used to model each level
in the field decomposition, the resulting model need not be identical to the
template. Application to a vortex census algorithm that records quantities of
interest (such as size, peak amplitude, circulation, etc.) as the vorticity
field evolves is given. The multiresolution census algorithm extracts coherent
structures of all shapes and sizes in simulated vorticity fields and is able to
reproduce known physical scaling laws when processing a set of voriticity
fields that evolve over time
Time scales in nuclear giant resonances
We propose a general approach to characterise fluctuations of measured cross
sections of nuclear giant resonances. Simulated cross sections are obtained
from a particular, yet representative self-energy which contains all
information about fragmentations. Using a wavelet analysis, we demonstrate the
extraction of time scales of cascading decays into configurations of different
complexity of the resonance. We argue that the spreading widths of collective
excitations in nuclei are determined by the number of fragmentations as seen in
the power spectrum. An analytic treatment of the wavelet analysis using a
Fourier expansion of the cross section confirms this principle. A simple rule
for the relative life times of states associated with hierarchies of different
complexity is given.Comment: 5 pages, 4 figure
PET Imaging of Post-infarct Myocardial Inflammation.
Funder: Department of HealthPurpose of reviewTo examine the use of positron emission tomography (PET) for imaging post-infarct myocardial inflammation and repair.Recent findingsDysregulated immune responses after myocardial infarction are associated with adverse cardiac remodelling and an increased likelihood of ischaemic heart failure. PET imaging utilising novel tracers can be applied to visualise different components of the post-infarction inflammatory and repair processes. This approach could offer unique pathophysiological insights that could prove useful for the identification and risk-stratification of individuals who would ultimately benefit most from emerging immune-modulating therapies. PET imaging could also bridge the clinical translational gap as a surrogate measure of drug efficacy in early-stage clinical trials in patients with myocardial infarction. The use of hybrid PET/MR imaging, in particular, offers the additional advantage of simultaneous in vivo molecular imaging and detailed assessment of myocardial function, viability and tissue characterisation. Further research is needed to realise the true clinical translational value of PET imaging after myocardial infarction
- …