RESUMEN

Addressing insufficient and irregular sampling is a difficult challenge in seismic processing and imaging. Recently, rank reduction methods have become popular in seismic processing algorithms for simultaneous denoising and interpolating. These methods are based on rank reduction of the trajectory matrices using truncated singular value decomposition (TSVD). Estimation of the ranks of these trajectory matrices depends on the number of plane waves in the processing window; however, for the more complicated data, the rank reduction method may fail or give poor results. In this paper, we propose an adaptive weighted rank reduction (AWRR) method that selects the optimum rank in each window automatically. The method finds the maximum ratio of the energy between two singular values. The AWRR method selects a large rank for the highly curved complex events, which leads to remaining residual errors. To overcome the residual errors, a weighting operator on the selected singular values minimizes the effect of noise projection on the signal projection. We tested the efficiency of the proposed method by applying it to both synthetic and real seismic data.

Asunto(s)

RESUMEN

To effectively find the direction of non-circular signals received by a uniform linear array (ULA) in the presence of non-negligible perturbations between array elements, i.e., mutual coupling, in colored noise, a direction of arrival (DOA) estimation approach in the context of high order statistics is proposed in this correspondence. Exploiting the non-circularity hidden behind a certain class of wireless communication signals to build up an augmented cumulant matrix, and carrying out a reformulation of the distorted steering vector to extract the angular information from the unknown mutual coupling, by exploiting the characteristic of mutual coupling, i.e., a limited operating range and an inverse relation of coupling effects to interspace, we develop a MUSIC-like estimator based on the rank-reduction (RARE) technique to directly determine directions of incident signals without mutual coupling compensation. Besides, we provide a solution to the problem of coherency between signals and mutual coupling between sensors co-existing, by selecting a middle sub-array to mitigate the undesirable effects and exploiting the rotation-invariant property to blindly separate the coherent signals into different groups to enhance the degrees of freedom. Compared with the existing robust DOA methods to the unknown mutual coupling under the framework of fourth-order cumulants (FOC), our work takes advantage of the larger virtual array and is able to resolve more signals due to greater degrees of freedom. Additionally, as the effective aperture is virtually extended, the developed estimator can achieve better performance under scenarios with high degree of mutual coupling between two sensors. Simulation results demonstrate the validity and efficiency of the proposed method.

RESUMEN

Model-based image reconstruction has improved contrast and spatial resolution in imaging applications such as magnetic resonance imaging and emission computed tomography. However, these methods have not succeeded in pulse-echo applications like ultrasound imaging due to the typical assumption of a finite grid of possible scatterer locations in a medium⁻an assumption that does not reflect the continuous nature of real world objects and creates a problem known as off-grid deviation. To cope with this problem, we present a method of dictionary expansion and constrained reconstruction that approximates the continuous manifold of all possible scatterer locations within a region of interest. The expanded dictionary is created using a highly coherent sampling of the region of interest, followed by a rank reduction procedure. We develop a greedy algorithm, based on the Orthogonal Matching Pursuit, that uses a correlation-based non-convex constraint set that allows for the division of the region of interest into cells of any size. To evaluate the performance of the method, we present results of two-dimensional ultrasound imaging with simulated data in a nondestructive testing application. Our method succeeds in the reconstructions of sparse images from noisy measurements, providing higher accuracy than previous approaches based on regular discrete models.

RESUMEN

Many modern statistical problems can be cast in the framework of multivariate regression, where the main task is to make statistical inference for a possibly sparse and low-rank coefficient matrix. The low-rank structure in the coefficient matrix is of intrinsic multivariate nature, which, when combined with sparsity, can further lift dimension reduction, conduct variable selection, and facilitate model interpretation. Using a Bayesian approach, we develop a unified sparse and low-rank multivariate regression method to both estimate the coefficient matrix and obtain its credible region for making inference. The newly developed sparse and low-rank prior for the coefficient matrix enables rank reduction, predictor selection and response selection simultaneously. We utilize the marginal likelihood to determine the regularization hyperparameter, so our method maximizes its posterior probability given the data. For theoretical aspect, the posterior consistency is established to discuss an asymptotic behavior of the proposed method. The efficacy of the proposed approach is demonstrated via simulation studies and a real application on yeast cell cycle data.

RESUMEN

Data from the Iowa mumps epidemic of 2006 were collected on a spatial lattice over a regular temporal interval. Without access to a person-to-person contact graph, it is sensible to analyze these data as homogenous within each areal unit and to use the spatial graph to derive a contact structure. The spatio-temporal partition is fine, and the counts of new infections at each location at each time are sparse. Therefore, we propose a spatial compartmental epidemic model with general latent time distributions (spatial PS SEIR) that is capable of smoothing the contact structure, while accounting for spatial heterogeneity in the mixing process between locations. Because the model is an extension of the PS SEIR model, it simultaneously handles non-exponentially distributed latent and infectious time distributions. The analysis within focuses on the progression of the disease over both space and time while assessing the impact of a large proportion of the infected people dispersing at the same time because of spring break and the impact of public awareness on the spread of the mumps epidemic. We found that the effect of spring break increased the mixing rate in the population and that the spatial transmission of the disease spreads across multiple conduits.

Asunto(s)

RESUMEN

Modern scientific research produces datasets of increasing size and complexity that require dedicated numerical methods to be processed. In many cases, the analysis of spectroscopic data involves the denoising of raw data before any further processing. Current efficient denoising algorithms require the singular value decomposition of a matrix with a size that scales up as the square of the data length, preventing their use on very large datasets. Taking advantage of recent progress on random projection and probabilistic algorithms, we developed a simple and efficient method for the denoising of very large datasets. Based on the QR decomposition of a matrix randomly sampled from the data, this approach allows a gain of nearly three orders of magnitude in processing time compared with classical singular value decomposition denoising. This procedure, called urQRd (uncoiled random QR denoising), strongly reduces the computer memory footprint and allows the denoising algorithm to be applied to virtually unlimited data size. The efficiency of these numerical tools is demonstrated on experimental data from high-resolution broadband Fourier transform ion cyclotron resonance mass spectrometry, which has applications in proteomics and metabolomics. We show that robust denoising is achieved in 2D spectra whose interpretation is severely impaired by scintillation noise. These denoising procedures can be adapted to many other data analysis domains where the size and/or the processing time are crucial.

Asunto(s)