Non-stationary predictive filtering for seismic random noise suppression - A tutorial

Predictive filtering in the frequency domain is one of the most widely used denoising algorithms in the seismic data processing workflow. Predictive filtering is based on the assumption of linear/planar events in the time-space domain. In traditional predictive filtering method, the predictive filter is fixed across the spatial dimension, which cannot deal with the spatial variation of seismic data well. To handle the curving events, the predictive filter is either applied in local windows or extended to a non-stationary version. The regularized non-stationary autoregression (RNAR) method can be treated as a non-stationary extension of the traditional predictive filtering, where the predictive filter coefficients are variable in different space locations. The highly under-determined inverse problem is solved by shaping regularization with a smoothness constraint in space. We further extend the RNAR method to a more general case, where we can apply more constraints to the filter coefficients according to the features of seismic data. First, apart from the smoothness in space, we also apply a smoothing constraint in frequency, considering the coherency of the coefficients in the frequency dimension. Secondly, we apply a frequency dependent smoothing radius along the space dimension to better take advantage of the non-stationarity of seismic data in the frequency axis, and to better deal with noise. The proposed method is validated via several synthetic and field data examples.

Spatially and Spectrally Concatenated Neural Networks for Efficient Lossless Compression of Hyperspectral Imagery

To achieve efficient lossless compression of hyperspectral images, we design a concatenated neural network, which is capable of extracting both spatial and spectral correlations for accurate pixel value prediction. Unlike conventional neural network based methods in the literature, the proposed neural network functions as an adaptive filter, thereby eliminating the need for pre-training using decompressed data. To meet the demand for low-complexity onboard processing, we use a shallow network with only two hidden layers for efficient feature extraction and predictive filtering. Extensive simulations on commonly used hyperspectral datasets and the standard CCSDS test datasets show that the proposed approach attains significant improvements over several other state-of-the-art methods, including standard compressors such as ESA, CCSDS-122, and CCSDS-123.

Automatically steering cardiac catheters in vivo with respiratory motion compensation

A robotic system for automatically navigating ultrasound (US) imaging catheters can provide real-time intra-cardiac imaging for diagnosis and treatment while reducing the need for clinicians to perform manual catheter steering. Clinical deployment of such a system requires accurate navigation despite the presence of disturbances including cyclical physiological motions (e.g., respiration). In this work, we report results from in vivo trials of automatic target tracking using our system, which is the first to navigate cardiac catheters with respiratory motion compensation. The effects of respiratory disturbances on the US catheter are modeled and then applied to four-degree-of-freedom steering kinematics with predictive filtering. This enables the system to accurately steer the US catheter and aim the US imager at a target despite respiratory motion disturbance. In vivo animal respiratory motion compensation results demonstrate automatic US catheter steering to image a target ablation catheter with 1.05 mm and 1.33° mean absolute error. Robotic US catheter steering with motion compensation can improve cardiac catheterization techniques while reducing clinician effort and X-ray exposure.

Three-parameter prestack seismic inversion based on L1-2 minimization

Prestack inversion has become a common approach in reservoir prediction. At present, the critical issue in the application of seismic inversion is the estimation of elastic parameters in the thin layers and weak reflectors. To improve the resolution and the accuracy of the inversion results, we introduced the difference of [Formula: see text] and [Formula: see text] norms as a nearly unbiased approximation of the sparsity of a vector, denoted as the [Formula: see text] norm, to the prestack inversion. The nonconvex penalty function of the [Formula: see text] norm can be decomposed into two convex subproblems via the difference of convex algorithm, and each subproblem can be solved efficiently by the alternating direction method of multipliers. Compared with the [Formula: see text] norm regularization, the [Formula: see text] minimization can reconstruct reflectivities more accurately. In addition, the [Formula: see text]-[Formula: see text] predictive filtering was introduced to guarantee the lateral continuity of the location and the amplitude of the reflectivity series. The generalized linear inversion and [Formula: see text]-[Formula: see text] predictive filtering are combined for stable elastic impedance inversion results, and three parameters of P-wave velocity, S-wave velocity, and density can be inverted with the Bayesian linearized amplitude variation with offset inversion. The inversion results of synthetic and real seismic data demonstrate that the proposed method can effectively improve the resolution and accuracy of the inversion results.