Liver cancers give rise to a heavy burden of health care worldwide. Understanding the tumour microenvironment (TME) underpins the development of precision therapy. Single-cell RNA sequencing (scRNA-seq) technology has generated high-quality cell atlases of the TME, but its wider application faces enormous costs for various clinical circumstances. Fortunately, a variety of deconvolution algorithms can instead repurpose bulk RNA-seq data, alleviating the need for generating scRNA-seq datasets. In this study, we reviewed major public omics databases for relevance in this study and utilized 8 RNA-seq and 1 microarray datasets from clinical studies. To decipher the TME of liver cancer, we estimated the fractions of liver cell components by deconvoluting the samples with Cibersortx using three reference scRNA-seq atlases. We also confirmed that Cibersortx can accurately deconvolute cell types/subtypes of interest. Compared with non-tumorous liver, liver cancers showed multiple decreased cell types forming normal liver microarchitecture, as well as elevated cell types involved in fibrogenesis, abnormal angiogenesis and disturbed immune responses. Survival analysis shows that the fractions of five cell types/subtypes significantly correlated with patient outcomes, indicating potential therapeutic targets. Therefore, deconvolution of bulk RNA-seq data with scRNA-seq atlas references can be a useful tool to help understand the TME.

Decision making is a complex process involving various parts of the brain which are active during different times. It is challenging to measure externally the exact instance when any given region becomes active during the decision-making process. Here we try to extract and visualize the dynamic functional brain activation information from the observed fMRI data. We propose the use of a regularized deconvolution model to simultaneously map various activation regions within the brain and track how different activation regions changes with time. Thus, providing both spatial and temporal brain activation information. The activation information can then be further analyzed as per need and requirements. The proposed technique was validated using simulated data and then applied to a simple decision-making task for identification of various brain regions involved in different stages of decision making. The visualization aspect of the algorithm allows us to actually see the flow of activation (and deactivation) in form of a motion picture. The dynamic estimate may aid in understanding the causality of activation between various brain regions in a better way.

In order to solve the problems of waveform distortion and signal delay by many physical and electrical systems with linear low-pass transfer characteristics with multiple complex poles, a general digital-signal-processing (DSP)-based method of real-time recovery of the original source waveform from the distorted output waveform is proposed. From the convolution kernel representation of a multiple-pole low-pass transfer function with an arbitrary denominator polynomial with real valued coefficients, it is shown that the source waveform can be accurately recovered in real time using a particular moving average algorithm with real-valued DSP computations only, even though some or all of the poles are complex. The proposed digital signal recovery method is DC-accurate and unaffected by initial conditions, transient signals, and resonant amplitude enhancement. This method can be applied to most sensors and amplifiers operating close to their frequency response limits or around their resonance frequencies to accurately deconvolute the multiple-pole characteristics and to improve the overall performances of data acquisition systems and digital feedback control systems.

To overcome the temperature effect of NaI(Tl) detectors for energy spectrometry without additional hardware, a new correction method was put forward based on pulse deconvolution, trapezoidal shaping and amplitude correction, named DTSAC. To verify this method, actual pulses from a NaI(Tl)-PMT detector were measured at different temperatures from -20℃ to 50℃. Pulse processing and spectrum synthesis showed that the position drift of the 662 keV peak was less than 3 keV, and the corresponding resolution at 662 keV of the sum spectra ranges from 6.9% to 10.60% with trapezoidal width set as 1000 ns to 100 ns. The DTSAC method corrects temperature effect by pulse processing, and need no reference peak, reference spectrum and additional circuits. The method solves the problem of correction of pulse shape and pulse amplitude at the same time, and can be used even at high counting rate.

A wealth of Full Waveform (FW) LiDAR data are available to the public from different sources, which is poised to boost the extensive application of FW LiDAR data. However, we lack a handy and open source tool that can be used by potential users for processing and analyzing FW LiDAR data. To this end, we introduce waveformlidar, an R package dedicated to FW LiDAR processing, analysis and visualization as a solution to the constraint. Specifically, this package provides several commonly used waveform processing methods such as Gaussian, adaptive Gaussian and Weibull decompositions, and deconvolution approaches (Gold and Richard-Lucy (RL)) with users’ customized settings. In addition, we also developed functions to derive commonly used waveform metrics for characterizing vegetation structure. Moreover, a new way to directly visualize FW LiDAR data is developed through converting waveforms into points to form the Hyper Point cloud (HPC), which can be easily adopted and subsequently analyzed with existing discrete-return LiDAR processing tools such as LAStools and FUSION. Basic explorations of the HPC such as 3D voxelization of the HPC and conversion from original waveforms to composite waveforms are also available in this package. All of these functions are developed based on small-footprint FW LiDAR data, but they can be easily transplanted to the large footprint FW LiDAR data such as Geoscience Laser Altimeter System (GLAS) and Global Ecosystem Dynamics Investigation (GEDI) data analysis. It is anticipated that these functions will facilitate the widespread use of FW LiDAR and be beneficial for better estimating biomass and characterizing vegetation structure at various scales. The package and code examples can be found at https://github.com/tankwin08/waveformlidar.

Visual Evoked Potentials (VEPs) are used in clinical applications in ophthalmology, neurology and extensively in brain computer interface (BCI) research. BCI literature covers steady state VEP (SSVEP) and code modulated VEP (c-VEP) BCIs along with sophisticated methods to improve information transfer rates (ITR). There is a gap of knowledge regarding the VEP adaptation dynamics, physiological generation mechanisms and relation with BCI performance. A simple dual display VEP switch was developed to test signatures elicited by non-isochronic, non-singular, low jitter stimuli at the rates of 10, 32, 50 and 70 reversals per second (rps). Non-isochronic, low-jitter stimulation elicits Quasi-Steady-State VEPs (QSS-VEPs) that are utilized for simultaneous generation of transient VEP and QSS-VEP. QSS-VEP is a special case of c-VEPs and it is assumed that it shares the similar generators of the SSVEPs. Eight subjects were recorded and the performance of the overall system was analyzed by means of Receiver Operating Characteristic (ROC) curves, accuracy plots and ITRs. In summary QSS-VEPs performed better than transient VEPs. It was found that in general 32rps stimulation had the highest ROC area, accuracy and ITRs in general. To investigate the reasons behind this, adaptation dynamics of transient VEPs and QSS-VEPs at all four rates were analyzed and speculated. Moreover, QSS-VEPs were found to lead to higher accuracy by the template matching compared to SSVEPs at 10rps and 32rps.

In the last two decades many detailed full transcriptomic studies on complex biological samples have been published and included in large gene expression repositories. These studies primarily provide bulk expression signal for each sample, including multiple cell-types mixed within the global signal. The cellular heterogeneity in these mixtures does not allow the activity of specific genes in specific cell types to be identified. Therefore, inferring relative cellular composition is a very powerful tool to achieve a more accurate molecular profiling of complex biological samples. In recent decades, computational techniques have been developed to solve this problem by applying cell deconvolution methods, designed to decompose cell mixtures into their cellular components and calculate the relative proportions of these elements. Some of them, only calculate the cell proportions (supervised methods), while other deconvolution algorithms can also identify the gene signatures specific for each cell type (unsupervised methods). In these work, five deconvolution methods (CIBERSORT, FARDEEP, DECONICA, LINSEED and ABIS) were implemented and used to analyze blood and immune cells, and also cancer cells, in complex mixture samples (using three bulk expression datasets). Our study provides three analytical tools (corrplots, cell-signature plots and bar-mixture plots) that allow a thorough comparative analysis of the cell mixture data. The work indicates that CIBERSORT is a robust method optimized for the identification of immune cell-types, but not as efficient in the identification of cancer cells. We also observed that LINSEED is a very powerful unsupervised method that provides precise and specific gene signatures for each of the five main immune cell-types analyzed: neutrophils and monocytes (in the myeloid lineage), B-cells, NK cells and T-cells (in the lymphoid lineage).

The retrieval of tree height percentiles with satellite lidar waveforms over the mountainous ar-eas is greatly challenged due to slanted and rugged topographic relief of underlying ground. A novel algorithm is proposed to extract the height percentiles based on the target response waveform (TRW) that is resolved by using a Richardson–Lucy deconvolution algorithm with an adaptive iterative time. The height metrics of the energy percentiles of 25%, 50%, 75% and 95% for the cumulative TRW distribution are determined by their vertical distances relative to the ground elevation in this study. To validate the proposed algorithm, we select the received waveforms of the Global Ecosystem Dynamics Investigation (GEDI) lidar over the Pahvant Mountains of central Utah, USA. The results reveal that the resolved TRWs resemble closely to the actual target response waveforms from the coincident airborne lidar data, with the mean values of the coefficient of correlation, total bias and root mean square error (RMSE) taking val-ue of 0.92, 0.0813 and 0.0016. In addition, the accuracies of the derived height percentiles from the proposed algorithm have been greatly improved compared with the conventional Gaussian de-composition method and the slope-adaptive waveform metrics method. The RMSE values de-crease by the mean values of 1.96 m and 2.72 m, respectively. It demonstrates that the proposed algorithm has good potential in the extraction of forest structure parameter over rugged moun-tainous areas.

Spatial biology is a rapidly growing research field which focuses on the transcriptomic or proteomic profiling of single cells within tissues with preserved spatial information. Imaging-based spatial transcriptomics uses epifluorescence microscopy, which has shown remarkable results for identification of multiple targets in situ. Nonetheless, the number of genes that can be reliably visualized is limited by the diffraction of light. Here, we investigate the effect of structured illumination (SIM), a super-resolution microscopy approach, over the performance of single gene transcript detection in spatial transcriptomics experiments. We performed direct mRNA-targeted hybridization in situ sequencing for multiple genes in mouse coronal brain tissue sections. We evaluated spot detection performance in widefield and confocal images versus those with SIM in combination with 20X, 25X and 60X objectives. In general, SIM increases the detection efficiency of gene transcripts spots compared to widefield and confocal modes. For each case, the specific fold increase in localizations is dependent on gene transcripts density and the numerical aperture of the objective used, which showed to play an important role especially for densely clustered spots. Taken together, our results suggest that SIM has the capacity to improve spot detection and overall data quality in spatial transcriptomics.

Blind deconvolution (BD) is one of the effective methods that help pre-process vibration signals and assist in bearing fault diagnosis. Currently, most BD methods design an optimization criterion and use frequency or time domain information independently to optimize a deconvolution filter. It recovers weak periodic impulses related to incipient faults. However, the random noise interference may cause the optimizer to overfit. The time-domain-based BD methods tend to extract fault-unrelated single peak impulse, and the frequency-domain-based BD methods tend to retain the maximum energy frequency component, which will lose the fault-related harmonics frequency components. To solve the above issue, we propose a hybrid criterion that combines the kurtosis for time domain optimization and the $G-l_1/l_2$ norm for the frequency domain. These two criteria are monotonically increasing and decreasing, so they mutually constrain to avoid overfitting. After that, we design a multi-task one-dimensional convolutional neural network with time and frequency branches to achieve an optimal solution for this hybrid criterion. The multi-task neural network realizes the simultaneous optimization of two domains. Experimental results show that our proposed method outperforms other state-of-the-art methods.

Photoacoustic spectroscopy allows the identification of specific molecules in gases. We evaluate the spectral resolution and detection limits for a PAS system in the broadband infrared wavelength region 3270 nm ≲ λ ≲ 3530 nm driven by a continuous wave optical parametric oscillator with P ≈ 1.26 W by measuring the absorption of diluted propane, ethane and methane test gases at low concentrations c ~ 100 ppm for ~1350 discrete wavelengths λ_i. The resulting spectra I_PAS(λ_i) were compared to the high resolution cross section data σ_FTIR obtained by Fourier Transform Infrared Spectroscopy from the HITRAN database. Deviations as little as 7.1(6)% for propane, 8.7(11)% for ethane and 15.0(14)% for methane with regard to the average uncertainty between I_PAS(λ_i) and the expected reference values based on σ_FTIR were recorded. The wavelengths λ_res of the characteristic absorption lines can be pinpointed with a high relative accuracy <5 × 10⁻⁵ corresponding to a resolution of λ_res ~ 0.16 nm. Detection limits range between 7.1 ppb (ethane) to 13.6 ppb (methane) coinciding with high experimental signal-to-noise ratios. Moreover, using EUREQA, an artificial intelligence program, simulated mixed gas samples at low limits of detection could be deconvoluted. These results justify a further development of PAS technology to support, e.g., biomedical research.

Cancer immunotherapy reaches a remarkable achievement in various cancer types and brings new possibilities to improve cancer patients’ long-term survival. However, outcomes vary from case to case, and the present protocol benefits a small fraction of patients. One notable factor is the tumor microenvironment, especially the immune cell components, that may reflect the immune response's status quo on site. Thus, understanding the content of infiltrating immune cells in tumors is not only for research interesting but also a crucial subject toward precision medicine. We implement an algorithm for resolving relative proportions of twenty-one immune cell subclasses from a human tissue profiled transcriptome by microarray technology to reach the goal above. By selecting gene features and then adopting ?-Support Vector Regression, we can construct a deconvolution model and resolve the immune cell context. The excellent consistency between the estimated values and the correct immune-cell composition further demonstrates this approach provides a more natural alternative to revealing samples' immune cell content and reliable results like recent single-cell technologies. Based on this algorithm, the web-based deconvolution tool implemented named mySORT provides a user-friendly interface for estimating the immune cell content by uploading gene expression profiling. We also present comprehensive visualization 2D/3D plots in mySORT so that users can easily make a comparison between different samples. Finally, we synthesized pseudo-bulk expression data from single-cell transcriptomic datasets of 17 melanoma and 16 head and neck cancer patients. The deconvolution results of microarray-based data in the previous study and synthetic pseudo-bulk data all proved the excellent performance of mySORT. We believe that mySORT can help researchers in all fields easily understand complex immune microenvironment. The website of mySORT is freely accessible on https://symbiosis.iis.sinica.edu.tw/mySORT/.

Detecting gear rim fatigue cracks using vibration signal analysis is often a challenging task, which typically requires a series of signal processing steps to detect and enhance fault features. This task becomes even harder in helicopter planetary gearboxes due to the complex interactions between different gear sets, and the presence of vibration from sources other than the planetary gear set. In this paper, we propose an effectual processing algorithm to isolate and enhance rim crack features and to trend crack growth in planet gears. The algorithm is based on using cepstrum editing (or liftering) of the hunting tooth synchronous averaged signals (angular domain) to extract harmonics and sidebands of the planet gears, and low-pass filtering and minimum entropy deconvolution (MED) to enhance extracted fault features. The algorithm has been successfully applied to a vibration dataset collected from a planet gear rim crack propagation test undertaken in the Helicopter Transmission Test Facility (HTTF) at DSTG Melbourne. In this test, a seeded notch generated by an electric discharge machine (EDM) was used to initiate a fatigue crack that propagated through the gear rim body over 94 load cycles. The proposed algorithm demonstrated a successful isolation of incipient fault features and provided a reliable trending capability to monitor crack progression. Results of a comparative analysis showed that the proposed algorithm outperformed the traditional signal processing approach.

The performance of bearings plays a pivotal role in determining the dependability and security of rotating machinery. In intricate systems demanding exceptional reliability and safety, the ability to accurately forecast fault occurrences during operation holds profound significance. Such predictions serve as invaluable guides for crafting well-considered reliability strategies and executing maintenance practices aimed at enhancing reliability. In order to ensure the reliability of bearing operation, this article investigates the application of three advanced techniques—Maximum Correlation Kurtosis Deconvolution (MCKD), Multi-Scale Permutation Entropy (MPE), and Long Short-Term Memory (LSTM) recurrent neural networks—for the prediction of the remaining useful life (RUL) of rolling bearings. Each technique's principles, methodologies, and applications are comprehensively reviewed, offering insights into their respective strengths and limitations. Case studies and experimental evaluations are presented to assess their performance in RUL prediction. Findings reveal that MCKD enhances fault signatures, MPE captures complexity, and LSTM excels in modeling temporal patterns. The root mean square error of the prediction results is 0.007. The fusion of these techniques offers a comprehensive approach to RUL prediction, leveraging their unique attributes for more accurate and reliable predictions.

Traditional discrimination techniques for Euler deconvolution use only the color spectrums of structural indexes, without considering the spatial distribution characteristics and inherent relationships among the Euler solutions to separate adjacent causative sources. In the present study, a new approach was developed for discriminating uncorrelated Euler solutions from coherent solutions based on the focusing levels indicated by the probability density distributions generated using multivariate kernel density estimations (KDE). A novel multiple coverage technique was proposed by using a series of different sized moving windows over gridded gravity data, which formed tight clusters of Euler solutions for different sized causative sources. The results of the probability density distributions were obtained using a 3-D KDE method for the Euler solution subsets {x, y, z} of synthetical models, and real data from a survey conducted in British Columbia (Canada) which had successfully established more credible and meaningful geological models when compared with three other subsets.

Properties of time domain parameters of the vibration signal have been extensively studied for the fault diagnosis of rolling element bearings (REB). Parameters like kurtosis and Envelope Harmonic-to-Noise Ratio are most widely applied in this field and some important progress has been made. However, since only one-sided information is contained in these parameters respectively, problems still exist in practice when the signals collected are of complicated structure and/or contaminated by strong background noises. A new parameter, named Shock pulse index (SPI), is proposed in this paper. It integrates the mutual advantage of both parameters above and can help effectively identify fault related impulse components under the interference of strong background noises, unrelated harmonic components and random impulses. The SPI optimizes the parameters of Maximum Correlated Kurtosis Deconvolution (MCKD), which is used to filter the signals under consideration. Finally, the interested transient information contained in the filtered signal can be highlighted through demodulation with Teager Energy Operator (TEO). Fault related impulse components can therefore be extracted accurately. Simulations and experiment analyses verify the effectiveness and correctness of the SPI.

Magnetic Resonance Imaging (MRI) is widely used in clinics and research due to its accurate disease identification and non-invasive nature. MRI is based on Nuclear Magnetic Resonance (NMR), which is also extensively employed in various fields. NMR and its modern versions involve intricate pre-processing steps before data analysis. These steps are initially processed in the time domain and subsequently in the frequency domain. While our previous review focused on time domain pre-processing (https://www.preprints.org/manuscript/202310.2032/v1), this review delves into the mathematical and statistical aspects of frequency domain pre-processing. We discuss essential pre-processing steps like phase error correction, baseline correction, solvent filtering, calibration and alignment, reference deconvolution, binning/bucketing and peak picking, peak fitting/deconvolution and compound identification, integration and quantification, normalization and transformation. Furthermore, we offer practical recommendations for each step.

A refractive lens is one of the simplest, cost-effective and easily available imaging elements. With a spatially incoherent illumination, a refractive lens can faithfully map every object point to an image point in the sensor plane, when the object and image distances satisfy the imaging conditions. However, static imaging is limited to the depth of focus, beyond which the point-to-point mapping can be only obtained by changing either the location of the lens or the imaging sensor. In this study, the depth of focus of a refractive lens in static mode has been expanded using a recently developed computational reconstruction method, Lucy-Richardson-Rosen algorithm (LRRA). The technique consists of three steps. In this first step, the point spread functions (PSFs) were recorded along different depths and stored in the computer as PSF library. In the next step, the object intensity distribution was recorded. The LRRA was then applied to deconvolve the object information from the recorded intensity distributions in the final step. The results of LRRA were compared against two well-known reconstruction methods namely Lucy-Richardson algorithm and non-linear reconstruction.

Deep Learning is a recent and important addition to the computational toolbox available for image reconstruction in fluorescence microscopy. We review state-of-the-art applications such as image restoration, super-resolution, and light-field imaging, and discuss how the latest Deep Learning research can be applied to other image reconstruction tasks such as structured illumination, spectral deconvolution, and sample stabilisation. Despite its successes, Deep Learning also poses significant challenges, has often misunderstood capabilities, and overlooked limits. We will address key questions, such as: What are the challenges in obtaining training data? Can we discover structures not present in the training data? And, what is the danger of inferring unsubstantiated image details?

This paper realizes early bearing fault warning through bearing fault time series prediction, and proposes a bearing fault time series prediction model based on optimized maximum correlation kurtosis deconvolution (MCKD) and long short-term memory (LSTM) recurrent neural network to ensure bearings operation reliability. The model is based on lifecycle vibration signal of the bearing, to begin, the cuckoo search (CS) is utilized to optimize the parameter filter length L and deconvolution period T of MCKD, taking into account the influence and periodicity of the bearing time series, the fault impact component of the optimized MCKD deconvolution time series is improved. Then select the LSTM learning rate α depending on deconvolution time series. Finally, the dataset obtained through various preprocessing approaches are used to train and predict the LSTM model. The average prediction accuracy of the optimized MCKD-LSTM model is 26 percent higher than that of the original time series, proving the efficiency of this method, and the prediction results track the real fault data well, according to the XI'AN JIAOTONG University XJTU-SY bearing dataset.

Synthetic Aperture RADAR (SAR) is a radar imaging technique in which the relative motion of the sensor is used to synthesize a very long antenna and obtain high spatial resolution. Standard SAR raw data processing techniques assume uniform motion of the satellite (or aerial vehicle) and a fixed antenna beam pointing sideway orthogonally to the motion path, assumed rectilinear. Despite SAR data processing is a well established imaging technology that has become fundamental in several fields and applications, in this paper a novel approach has been used to exploit coherent illumination, demonstrating the possibility of extracting a large part of the ancillary data information from the raw data itself, to be used in the focusing procedure. In this work an effort has been carried out to try to focus the raw SAR complex data matrix without the knowledge of anyof the parameters needed in standard focusing procedures as Range Doppler (RD) algorithm, ω - K algorithm and Chirp Scaling (CS) algorithm. All the literature references regarding the algorithms needed to obtain a precise image from raw data use such parameters that refer both to the SAR system acquisition geometry and its radiometric specific parameters. In [12], authors introduced a preliminary work dealing with this problem and able to obtain, in the presence of a strong point scatterer in the observed scene, good quality images, if compared to the standard processing techniques. In this work the proposed technique is described and performances parameters are extracted to compare the proposed approach to RD.