Do you want to publish a course? Click here

Fourier-Bessel rotational invariant eigenimages

95   0   0.0 ( 0 )
 Added by Zhizhen Zhao
 Publication date 2012
and research's language is English




Ask ChatGPT about the research

We present an efficient and accurate algorithm for principal component analysis (PCA) of a large set of two dimensional images, and, for each image, the set of its uniform rotations in the plane and its reflection. The algorithm starts by expanding each image, originally given on a Cartesian grid, in the Fourier-Bessel basis for the disk. Because the images are bandlimited in the Fourier domain, we use a sampling criterion to truncate the Fourier-Bessel expansion such that the maximum amount of information is preserved without the effect of aliasing. The constructed covariance matrix is invariant to rotation and reflection and has a special block diagonal structure. PCA is efficiently done for each block separately. This Fourier-Bessel based PCA detects more meaningful eigenimages and has improved denoising capability compared to traditional PCA for a finite number of noisy images.



rate research

Read More

The spherical Fourier-Bessel (SFB) decomposition is a natural choice for the radial/angular separation that allows optimal extraction of cosmological information from large volume galaxy surveys. In this paper we develop a SFB power spectrum estimator that allows the measurement of the largest angular and radial modes with the next generation of galaxy surveys. The code measures the pseudo-SFB power spectrum, and takes into account mask, selection function, pixel window, and shot noise. We show that the local average effect is significant only in the largest-scale mode, and we provide an analytical covariance matrix. By imposing boundary conditions at the minimum and maximum radius encompassing the survey volume, the estimator does not suffer from the numerical instabilities that have proven challenging in the past. The estimator is demonstrated on simplified Roman-like, SPHEREx-like, and Euclid-like mask and selection functions. For intuition and validation, we also explore the SFB power spectrum in the Limber approximation. We release the associated public code written in Julia.
We describe a novel scheme of detecting rotational anisotropy second harmonic generation (RA-SHG) signals using a lock-in amplifier referenced to a fast scanning RA-SHG apparatus. The method directly measures the $n^{th}$ harmonics of the scanning frequency corresponding to SHG signal components of $C_n$ symmetry that appear in a Fourier series expansion of a general RA-SHG signal. GaAs was used as a test sample allowing comparison of point-by-point averaging with the lock-in based method. When divided by the $C_infty$ signal component, the lock-in detected data allowed for both self-referenced determination of ratios of $C_n$ components of up to 1 part in 10$^4$ and significantly more sensitive measurement of the relative amount of different $C_n$ components when compared with conventional methods.
This paper makes two points. First, we show that the line-of-sight solution to cosmic microwave anisotropies in Fourier space, even though formally defined for arbitrarily large wavelengths, leads to position-space solutions which only depend on the sources of anisotropies inside the past light-cone of the observer. This happens order by order in a series expansion in powers of the visibility $gamma=e^{-mu}$, where $mu$ is the optical depth to Thompson scattering. We show that the CMB anisotropies are regulated by spacetime window functions which have support only inside the past light-cone of the point of observation. Second, we show that the Fourier-Bessel expansion of the physical fields (including the temperature and polarization momenta) is an alternative to the usual Fourier basis as a framework to compute the anisotropies. In that expansion, for each multipole $l$ there is a discrete tower of momenta $k_{i,l}$ (not a continuum) which can affect physical observables, with the smallest momenta being $k_{1,l} ~ l$. The Fourier-Bessel modes take into account precisely the information from the sources of anisotropies that propagates from the initial value surface to the point of observation - no more, no less. We also show that the physical observables (the temperature and polarization maps), and hence the angular power spectra, are unaffected by that choice of basis. This implies that the Fourier-Bessel expansion is the optimal scheme with which one can compute CMB anisotropies. (Abridged)
In a recent paper on exotic crossed products, we included a lemma concerning ideals of the Fourier-Stieltjes algebra. Buss, Echterhoff, and Willett have pointed out to us that our proof of this lemma contains an error. In fact, it remains an open question whether the lemma is true as stated. In this note we indicate how to contain the resulting damage. Our investigation of the above question leads us to define two properties emph{ordered} and emph{weakly ordered} for invariant ideals of Fourier-Stieltjes algebras, and we initiate a study of these properties.
High-precision cosmology requires the analysis of large-scale surveys in 3D spherical coordinates, i.e. spherical Fourier-Bessel decomposition. Current methods are insufficient for future data-sets from wide-field cosmology surveys. The aim of this paper is to present a public code for fast spherical Fourier-Bessel decomposition that can be applied to cosmological data or 3D data in spherical coordinates in other scientific fields. We present an equivalent formulation of the spherical Fourier-Bessel decomposition that separates radial and tangential calculations. We propose the use of the existing pixelisation scheme HEALPix for a rapid calculation of the tangential modes. 3DEX (3D EXpansions) is a public code for fast spherical Fourier-Bessel decomposition of 3D all-sky surveys that takes advantage of HEALPix for the calculation of tangential modes. We perform tests on very large simulations and we compare the precision and computation time of our method with an optimised implementation of the spherical Fourier-Bessel original formulation. For surveys with millions of galaxies, computation time is reduced by a factor 4-12 depending on the desired scales and accuracy. The formulation is also suitable for pre-calculations and external storage of the spherical harmonics, which allows for additional speed improvements. The 3DEX code can accommodate data with masked regions of missing data. 3DEX can also be used in other disciplines, where 3D data are to be analysed in spherical coordinates. The code and documentation can be downloaded at http://ixkael.com/blog/3dex.
comments
Fetching comments Fetching comments
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا