2D and 3D prestack seismic data regularization using an accelerated sparse time-invariant Radon transform

We present a fast and robust sparse time-invariant Radon transform (RT) in the mixed frequency-time domain. In the proposed method, the 2D alternating split Bregman algorithm is introduced to solve the optimization problem. Two proximity operators are applied to promote the robustness and sparseness of the RT in the time domain, while the forward and inverse RTs are computed in the frequency domain. With high efficiency and easy implementation, the proposed method is resistant to the outliers. The synthetic and real data examples demonstrate the better performance of the proposed method when applied on the noise-contaminated data, compared with the traditional non-robust RTs.

In the early stage of acoustic full-waveform inversion (AFWI), it is important to exploit the long-wavelength features of the gradient and suppress its short-wavelength features to update the background velocity. However, due to strong near-offset PP reflections and multiples within the pressure data, conventional AFWI sometimes primarily reconstructs the short-wavelength features of a given model and fails to converge on a reliable subsurface P-wave velocity model. Therefore, this study, we propose the usage of pseudo-horizontal particle acceleration (pseudo-ax) data for AFWI, rather than the original pressure data. Because the amplitudes of PP reflections are implicitly weighted according to the reflection angle in pseudo-ax data, the near-offset PP reflections and multiples of the pseudo-ax data are much weaker than those of the pressure data. As a result, AFWI using pseudo-ax data focuses on reconstructing the longer-wavelength features of a given model in the early iterations, and then gradually reconstructs short-wavelength features as the inversion process continues. Using synthetic and real data examples for the Volve oilfield in the North Sea, we examined the effects and feasibility of this strategy. In the synthetic data example, the proposed strategy rapidly reduced far-offset data misfits related to long-wavelength feature updates and achieved smaller model misfits than the conventional strategy. The real data example showed that the velocity model reconstructed using the proposed strategy flattened some of the reflectors at the angle-domain common-image gathers (ADCIGs) better than the velocity model obtained using the conventional strategy. These findings demonstrate that our strategy converges closer to the real P-wave velocity model than the conventional strategy by building a hierarchical velocity model from long- to short-wavelength structures.

In transmission X-ray microscopy (TXM) systems, the rotation of a scanned sample might be restricted to a limited angular range to avoid collision with other system parts or high attenuation at certain tilting angles. Image reconstruction from such limited angle data suffers from artifacts because of missing data. In this work, deep learning is applied to limited angle reconstruction in TXMs for the first time. With the challenge to obtain sufficient real data for training, training a deep neural network from synthetic data is investigated. In particular, U-Net, the state-of-the-art neural network in biomedical imaging, is trained from synthetic ellipsoid data and multi-category data to reduce artifacts in filtered back-projection (FBP) reconstruction images. The proposed method is evaluated on synthetic data and real scanned chlorella data in 100° limited angle tomography. For synthetic test data, U-Net significantly reduces the root-mean-square error (RMSE) from 2.55 × 10-3 µm-1 in the FBP reconstruction to 1.21 × 10-3 µm-1 in the U-Net reconstruction and also improves the structural similarity (SSIM) index from 0.625 to 0.920. With penalized weighted least-square denoising of measured projections, the RMSE and SSIM are further improved to 1.16 × 10-3 µm-1 and 0.932, respectively. For real test data, the proposed method remarkably improves the 3D visualization of the subcellular structures in the chlorella cell, which indicates its important value for nanoscale imaging in biology, nanoscience and materials science.

Seismic inverse problems aim to infer the properties of subsurface geology, such as elastic and petrophysical properties. Existing seismic inversion methods for the joint estimation of these properties are mainly based either on Gassmann theory for prestack seismic data processed with stochastic optimisation techniques or on the Wyllie formula for poststack seismic data processed by deterministic optimisation techniques. The purpose of this study is to develop a strategy for the joint estimation of elastic and petrophysical properties from prestack seismic data based on Gassmann equations with deterministic optimisation techniques. Given poor-quality prestack seismic data, two regularisation parameters are introduced to control the trade-off between fidelity to the data and the smoothness of the solution. An appropriate linearised system of equations for the joint model update is derived from Newton’s method, which fits seismic data, the description of the rock physics medium and prior information, simultaneously. We show the preliminary results obtained with the proposed framework for synthetic and real data examples.

Accelerated parabolic Radon domain 2D adaptive multiple subtraction with fast iterative shrinkage thresholding algorithm and its application in parabolic Radon domain hybrid demultiple method

Infrared (IR) images are essential to improve the visibility of dark or camouflaged objects. Object recognition and segmentation based on a neural network using IR images provide more accuracy and insight than color visible images. But the bottleneck is the amount of relevant IR images for training. It is difficult to collect real-world IR images for special purposes, including space exploration, military and fire-fighting applications. To solve this problem, we created color visible and IR images using a Unity-based 3D game editor. These synthetically generated color visible and IR images were used to train cycle consistent adversarial networks (CycleGAN) to convert visible images to IR images. CycleGAN has the advantage that it does not require precisely matching visible and IR pairs for transformation training. In this study, we discovered that additional synthetic data can help improve CycleGAN performance. Neural network training using real data (N = 20) performed more accurate transformations than training using real (N = 10) and synthetic (N = 10) data combinations. The result indicates that the synthetic data cannot exceed the quality of the real data. Neural network training using real (N = 10) and synthetic (N = 100) data combinations showed almost the same performance as training using real data (N = 20). At least 10 times more synthetic data than real data is required to achieve the same performance. In summary, CycleGAN is used with synthetic data to improve the IR image conversion performance of visible images.

We address the directional derivative trace transform (DDTT) for seismic wavefield separation. The goal is to separate the different seismic waves, which is a practical requirement for seismic exploration and geophysics. The DDTT is based on the trace transform, which is a generalization of the Radon transform (RT). It calculates the integral of a directional derivative along the trace line. The directional derivative is the pointwise rate of change of a function in a certain direction. Therefore, the DDTT better reflects the property of linear wavefields and other wavefields than the F-K filtering and RT methods. The DDTT domain consists of two parts: One part mostly represents linear wavefields such as surface waves, and the other part mostly represents the reflection signals. Based on this property, a parametric method is proposed for separating the ground roll better. In order to return data to the time-offset domain, we derive the inverse DDTT using the properties of the Fourier transform and the Hilbert transform. We show the potential of our method for the removal of ground roll on synthetic and real data examples.

Ultra-low dose abdominal MDCT: Using a knowledge-based Iterative Model Reconstruction technique for substantial dose reduction in a prospective clinical study

Deep-learning-based methods have been successfully applied to seismic data random noise attenuation. Among them, the supervised deep-learning-based methods dominate the unsupervised ones. The supervised methods need accurate noise-free data as training labels. However, the field seismic data cannot meet this requirement. To circumvent it, some researchers utilized realistic-looking synthetic data or denoised results via conventional methods as labels. The former ones encounter the problem of weak generalization ability because it requires the same distribution of test and training data. The latter ones encounter the issue of insufficient denoising ability because its denoising ability is difficult to significantly exceed the conventional methods which were used to generate labels. To avoid preparing noise-free labels, we propose a novel deep learning framework for attenuating random noise of prestack seismic data in an unsupervised manner. The prestack seismic data, such as common-reflection-point (CRP) gathers and common-midpoint (CMP) gathers after normal moveout (NMO) correction, have high self-similarity. It is because their events are coherent in the time–space domain and approximately horizontal from shallow to deep layers. The generator convolutional neural network (GCN) first learns self-similar features before any learning. The useful signals are more self-similar than random noise, which is incoherent and randomly distributed. Therefore, the GCN extracts features of useful signals before random noise. We select the specific training iteration and adopt the early stopping strategy to suppress random noise. Both synthetic and field prestack seismic data examples demonstrate the validity of our methods.

Coronary artery calcium (CAC) scoring with computed tomography (CT) is an established tool for quantifying calcified atherosclerotic plaque burden. Despite the widespread use of novel image reconstruction techniques in CT, the effect of iterative model reconstruction on CAC score remains unclear. We sought to assess the impact of iterative model based reconstruction (IMR) on coronary artery calcium quantification as compared to the standard filtered back projection (FBP) algorithm and hybrid iterative reconstruction (HIR). In addition, we aimed to simulate the impact of iterative reconstruction techniques on calcium scoring based risk stratification of a larger asymptomatic population. We studied 63 individuals who underwent CAC scoring. Images were reconstructed with FBP, HIR and IMR and CAC scores were measured. We estimated the cardiovascular risk reclassification rate of IMR versus HIR and FBP in a larger asymptomatic population (n = 504). The median CAC scores were 147.7 (IQR 9.6-582.9), 107.0 (IQR 5.9-526.6) and 115.1 (IQR 9.3-508.3) for FBP, HIR and IMR, respectively. The HIR and IMR resulted in lower CAC scores as compared to FBP (both p < 0.001), however there was no difference between HIR and IMR (p = 0.855). The CAC score decreased by 7.2 % in HIR and 7.3 % in IMR as compared to FBP, resulting in a risk reclassification rate of 2.4 % for both HIR and IMR. The utilization of IMR for CAC scoring reduces the measured calcium quantity. However, the CAC score based risk stratification demonstrated modest reclassification in IMR and HIR versus FBP.

After multiple prediction, adaptive multiple subtraction is essential for the success of multiple removal. The 3D blind separation of convolved mixtures (3D BSCM) method, which is effective in conducting adaptive multiple subtraction, needs to solve an optimization problem containing L1-norm minimization constraints on primaries by the iterative reweighted least-squares (IRLS) algorithm. The 3D BSCM method can better separate primaries and multiples than the 1D/2D BSCM method and the method with energy minimization constraints on primaries. However, the 3D BSCM method has high computational cost because the IRLS algorithm achieves nonquadratic optimization with an LS optimization problem solved in each iteration. In general, it is good to have a faster 3D BSCM method. To improve the adaptability of field data processing, the fast iterative shrinkage thresholding algorithm (FISTA) is introduced into the 3D BSCM method. The proximity operator of FISTA can solve the L1-norm minimization problem efficiently. We demonstrate that our FISTA-based 3D BSCM method achieves similar accuracy of estimating primaries as that of the reference IRLS-based 3D BSCM method. Furthermore, our FISTA-based 3D BSCM method reduces computation time by approximately 60% compared with the reference IRLS-based 3D BSCM method in the synthetic and field data examples.

Seismic derivable elastic attributes, e.g., elastic impedance, lambda-rho, mu-rho, and Poisson impedance (PI), are routinely being used for reservoir characterization practice. These attributes could be derived from inverted VP, VS, and density, and usually indicate high sensitivity to reservoir lithology and fluid. Due to the high sensitivity of such elastic attributes, errors or measurement noise associated with the acquisition, processing, and inversion of prestack seismic data will propagate through the inversion products, and will lead to even larger errors in the computed attributes. To solve this problem, we have developed a two-step cascade workflow that combines linear inversion and nonlinear optimization techniques for the improved estimation of elastic attributes and better prediction and delineation of reservoir lithology and fluids. The linear inversion in the first step is an inversion scheme with a sparseness assumption, based on L1-norm regularization. This step is used to select the major reflective layer locations, followed in the second step by a nonlinear optimization process with the predefined layer structure. The combination of these two procedures produces a reasonable blocky earth model with consistent elastic properties, including the ones that are sensitive to reservoir lithology and fluid change, and thus provides an accurate approach for seismic reservoir characterization. Using PI, as one of the target elastic attributes, as an example, this workflow has been successfully applied to synthetic and field data examples. The results indicate that our workflow improves the estimation of elastic attributes from the noisy prestack seismic data and may be used for the identification of the reservoir lithology and fluid.

With palm-top devices increasingly popular, security for both the data contained in the devices and e-commerce transactions through them becomes an issue. Iris texture is a nature-given password that has advantages over alphanumeric password and other biometric features. Existing iris recognition systems however are all designed for either mainframe or desktop applications. This paper describes an iris recognition algorithm that is compact and simple enough to be implemented on palm-top devices. The algorithm exploits a simple transformation – Radon Transform – to encode iris images. Experimental results on real image data are shown to illustrate how the proposed system performs in comparison with an existing Gabor wavelet coefficient based algorithm.

Accurately modelled computer-generated data can be used in place of real-world signals for the design, test and validation of signal processing techniques in situations where real data is difficult to obtain. Bio-signal processing researchers interested in working with fNIRS data are restricted due to the lack of freely available fNIRS data and by the prohibitively expensive cost of fNIRS systems. We present a simplified mathematical description and associated MATLAB implementation of model-based synthetic fNIRS data which could be used by researchers to develop fNIRS signal processing techniques. The software, which is freely available, allows users to generate fNIRS data with control over a wide range of parameters and allows for fine-tuning of the synthetic data. We demonstrate how the model can be used to generate raw fNIRS data similar to recorded fNIRS signals. Signal processing steps were then applied to both the real and synthetic data. Visual comparisons between the temporal and spectral properties of the real and synthetic data show similarity. This paper demonstrates that our model for generating synthetic fNIRS data can replicate real fNIRS recordings.

In comparison with conventional filtered backprojection (FBP) algorithms for x-ray computed tomography (CT) image reconstruction, statistical algorithms directly incorporate the random nature of the data and do not assume CT data are linear, noiseless functions of the attenuation line integral. Thus, it has been hypothesized that statistical image reconstruction may support a more favorable tradeoff than FBP between image noise and spatial resolution in dose-limited applications. The purpose of this study is to evaluate the noise-resolution tradeoff for the alternating minimization (AM) algorithm regularized using a nonquadratic penalty function. Idealized monoenergetic CT projection data with Poisson noise were simulated for two phantoms with inserts of varying contrast (7%-238%) and distance from the field-of-view (FOV) center (2-6.5 cm). Images were reconstructed for the simulated projection data by the FBP algorithm and two penalty function parameter values of the penalized AM algorithm. Each algorithm was run with a range of smoothing strengths to allow quantification of the noise-resolution tradeoff curve. Image noise is quantified as the standard deviation in the water background around each contrast insert. Modulation transfer functions (MTFs) were calculated from six-parameter model fits to oversampled edge-spread functions defined by the circular contrast-insert edges as a metric of local resolution. The integral of the MTF up to 0.5 1p/mm was adopted as a single-parameter measure of local spatial resolution. The penalized AM algorithm noise-resolution tradeoff curve was always more favorable than that of the FBP algorithm. While resolution and noise are found to vary as a function of distance from the FOV center differently for the two algorithms, the ratio of noises when matching the resolution metric is relatively uniform over the image. The ratio of AM-to-FBP image variances, a predictor of dose-reduction potential, was strongly dependent on the shape of the AM's nonquadratic penalty function and was also strongly influenced by the contrast of the insert for which resolution is quantified. Dose-reduction potential, reported here as the fraction (%) of FBP dose necessary for AM to reconstruct an image with comparable noise and resolution, for one penalty parameter value of the AM algorithm was found to vary from 70% to 50% for low-contrast and high-contrast structures, respectively, and from 70% to 10% for the second AM penalty parameter value. However, the second penalty, AM-700, was found to suffer from poor low-contrast resolution when matching the high-contrast resolution metric with FBP. The results of this simulation study imply that penalized AM has the potential to reconstruct images with similar noise and resolution using a fraction (10%-70%) of the FBP dose. However, this dose-reduction potential depends strongly on the AM penalty parameter and the contrast magnitude of the structures of interest. In addition, the authors' results imply that the advantage of AM can be maximized by optimizing the nonquadratic penalty function to the specific imaging task of interest. Future work will extend the methods used here to quantify noise and resolution in images reconstructed from real CT data.