BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Hau-Tieng Wu

Hau-Tieng Wu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAmath.STStatistics TheoryApplicationsMachine Learningeess.SPMethodologyphysics.data-anNumerical Analysismath.SPphysics.med-phQuantitative Methodsmath.CAmath.OCmath.PRphysics.opticsSoundComputational Engineering, Finance, and Scienceeess.IVHuman-Computer Interactionmath.CVmath.FAMultimediaNeurons and Cognition

Catalog footprint

What is connected

46works
24topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

46 published item(s)

preprint2026arXiv

On spectral interference of the short-time Fourier transform and its nonlinear variations

Spectral interference, the frequency counterpart of the beating phenomenon in the time domain, can severely distort time-frequency representations (TFRs) in physical applications. We study this phenomenon for the short-time Fourier transform (STFT) with a Gaussian window and for nonlinear refinements based on the reassignment method, with an emphasis on the synchrosqueezing transform (SST). Working with a two-component harmonic model, we quantify when STFT can (and cannot) resolve two nearby frequencies: a sharp transition occurs at a critical gap that scales inversely to kernel bandwidth and depends explicitly on the amplitude ratio. Below this threshold, the spectrogram ridges undergo bifurcation and form repeating time-frequency bubbles, which we describe asymptotically and, in the balanced-amplitude case, approximate closely by ellipses. We then analyze the STFT phase, showing a canonical winding behavior, and relate the complex-valued SST reassignment map to a holomorphic structure via the Bargmann transform. In the two-component setting the reassignment rule admits an explicit Mobius-geometry description, sending frequency lines to circular arcs in the complex plane. Finally, viewing SST and reassignment through a measure mapping perspective, we derive small-kernel asymptotics that explain when reassignment sharpens energy and when it produces distorted or misleading TFRs; we also introduce a generalized synchrosqueezing framework that isolates the role of STFT weighting and clarifies how alternative choices can mitigate interference in certain regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Probabilistic Analysis of Scalogram Ridges in Signal Processing

While ridges in the scalogram, determined by the squared modulus of analytic wavelet transform (AWT), is a widely accepted concept and utilized in nonstationary time series analysis, their behavior in noisy environments remains underexplored. Our object is to provide a theoretical foundation for scalogram ridges by defining ridges as a potentially set-valued random process connecting local maxima of the scalogram along the scale axis and analyzing their properties when the signal fulfills the adaptive harmonic model and is contaminated by stationary Gaussian noise. In addition to establishing several key properties of the AWT for random processes, we investigate the probabilistic characteristics of the resulting random ridge points in the scalogram. Specifically, we establish the uniqueness property of the ridge point at individual time instances and prove the upper hemicontinuity of the ridge random process. Furthermore, we derive bounds on the probability that the deviation between the ridges of noisy and clean signals exceeds a specified threshold, and these bounds depend on the signal-to-noise ratio. To achieve these ridge deviation results, we derive maximal inequalities for the complex modulus of nonstationary Gaussian processes, leveraging classical tools such as the Borell-TIS inequality and Dudley's theorem, which might be of independent interest.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Efficient Artifacts Removal for Adaptive Deep Brain Stimulation and a Temporal Event Localization Analysis

Adaptive deep brain stimulation (aDBS) leverages symptom-related biomarkers to deliver personalized neuromodulation therapy, with the potential to improve treatment efficacy and reduce power consumption compared to conventional DBS. However, stimulation-induced signal contamination remains a major technical barrier to advancing its clinical application. Existing artifact removal strategies, both front-end and back-end, face trade-offs between artifact suppression and algorithmic flexibility. Among back-end algorithms, Shrinkage and Manifold-based Artifact Removal using Template Adaptation (SMARTA) has shown promising performance in mitigating stimulus artifacts with minimal distortion to local field potentials (LFPs), but its high computational demand and inability to handle transient direct current (DC) artifacts limit its use in real-time applications. To address this, we developed SMARTA+, a computationally efficient extension of SMARTA capable of suppressing both stimulus and transient DC artifacts while supporting flexible algorithmic design. We evaluated SMARTA+ using semi-real aDBS data and real data from Parkinson's disease patients. Compared to SMARTA and other established methods, SMARTA+ achieved comparable or superior artifact removal while significantly reducing computation time. It preserved spectral and temporal structures, ranging from beta band to high-frequency oscillations, and demonstrated robustness across diverse stimulation protocols. Temporal event localization analysis further showed improved accuracy in detecting beta bursts. These findings support SMARTA+ as a promising tool for advancing real-time, closed-loop aDBS systems.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Design a Metric Robust to Complicated High Dimensional Noise for Efficient Manifold Denoising

In this manuscript, we propose an efficient manifold denoiser based on landmark diffusion and optimal shrinkage under the complicated high dimensional noise and compact manifold setup. It is flexible to handle several setups, including the high ambient space dimension with a manifold embedding that occupies a subspace of high or low dimensions, and the noise could be colored and dependent. A systematic comparison with other existing algorithms on both simulated and real datasets is provided. This manuscript is mainly algorithmic and we report several existing tools and numerical results. Theoretical guarantees and more comparisons will be reported in the official paper of this manuscript.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Disentangling modes with crossover instantaneous frequencies by synchrosqueezed chirplet transforms, from theory to application

Analysis of signals with oscillatory modes with crossover instantaneous frequencies is a challenging problem in time series analysis. One way to handle this problem is lifting the 2-dimensional time-frequency representation to a 3-dimensional representation, called time-frequency-chirp rate (TFC) representation, by adding one extra chirp rate parameter so that crossover frequencies are disentangled in higher dimension. The chirplet transform is an algorithm for this lifting idea, which leads to a TFC representation. However, in practice, we found that it has a strong ``blurring'' effect in the chirp rate axis, which limits its application in real-world data. Moreover, to our knowledge, we have limited mathematical understanding of the chirplet transform in the literature. Motivated by the need for the real-world data analysis, in this paper, we propose the synchrosqueezed chirplet transform (SCT) that enhances the TFC representation given by the chirplet transform. The resulting concentrated TFC representation has high contrast so that one can better distinguish different modes with crossover instantaneous frequencies. The basic idea is to use the phase information in the chirplet transform to determine a reassignment rule that sharpens the TFC representation determined by the chirplet transform. We also analyze the chirplet transform and provide theoretical guarantees of SCT.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Eigenvector Phase Retrieval: Recovering eigenvectors from the absolute value of their entries

We consider the eigenvalue problem $Ax = λx$ where $A \in \mathbb{R}^{n \times n}$ and the eigenvalue is also real $λ\in \mathbb{R}$. If we are given $A$, $λ$ and, additionally, the absolute value of the entries of $x$ (the vector $(|x_i|)_{i=1}^n$), is there a fast way to recover $x$? In particular, can this be done quicker than computing $x$ from scratch? This may be understood as a special case of the phase retrieval problem. We present a randomized algorithm which provably converges in expectation whenever $λ$ is a simple eigenvalue. The problem should become easier when $|λ|$ is large and we discuss another algorithm for that case as well.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Predicting Trust Using Automated Assessment of Multivariate Interactional Synchrony

Diverse disciplines are interested in how the coordination of interacting agents' movements, emotions, and physiology over time impacts social behavior. Here, we describe a new multivariate procedure for automating the investigation of this kind of behaviorally-relevant "interactional synchrony", and introduce a novel interactional synchrony measure based on features of dynamic time warping (DTW) paths. We demonstrate that our DTW path-based measure of interactional synchrony between facial action units of two people interacting freely in a natural social interaction can be used to predict how much trust they will display in a subsequent Trust Game. We also show that our approach outperforms univariate head movement models, models that consider participants' facial action units independently, and models that use previously proposed synchrony or similarity measures. The insights of this work can be applied to any research question that aims to quantify the temporal coordination of multiple signals over time, but has immediate applications in psychology, medicine, and robotics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Scalability and robustness of spectral embedding: landmark diffusion is all you need

While spectral embedding is a widely applied dimension reduction technique in various fields, so far it is still challenging to make it scalable to handle ``big data''. On the other hand, the robustness property is less explored and there exists only limited theoretical results. Motivated by the need of handling such data, recently we proposed a novel spectral embedding algorithm, which we coined Robust and Scalable Embedding via Landmark Diffusion (ROSELAND). In short, we measure the affinity between two points via a set of landmarks, which is composed of a small number of points, and ``diffuse'' on the dataset via the landmark set to achieve a spectral embedding. Roseland can be viewed as a generalization of the commonly applied spectral embedding algorithm, the diffusion map (DM), in the sense that it shares various properties of DM. In this paper, we show that Roseland is not only numerically scalable, but also preserves the geometric properties via its diffusion nature under the manifold setup; that is, we theoretically explore the asymptotic behavior of Roseland under the manifold setup, including handling the U-statistics-like quantities, and provide a $L^\infty$ spectral convergence with a rate. Moreover, we offer a high dimensional noise analysis and show that Roseland is robust to noise. We also compare Roseland with other existing algorithms with numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Spatiotemporal Analysis Using Riemannian Composition of Diffusion Operators

Multivariate time-series have become abundant in recent years, as many data-acquisition systems record information through multiple sensors simultaneously. In this paper, we assume the variables pertain to some geometry and present an operator-based approach for spatiotemporal analysis. Our approach combines three components that are often considered separately: (i) manifold learning for building operators representing the geometry of the variables, (ii) Riemannian geometry of symmetric positive-definite matrices for multiscale composition of operators corresponding to different time samples, and (iii) spectral analysis of the composite operators for extracting different dynamic modes. We propose a method that is analogous to the classical wavelet analysis, which we term Riemannian multi-resolution analysis (RMRA). We provide some theoretical results on the spectral analysis of the composite operators, and we demonstrate the proposed method on simulations and on real data.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An Efficient Forecasting Approach to Reduce Boundary Effects in Real-Time Time-Frequency Analysis

Time-frequency (TF) representations of time series are intrinsically subject to the boundary effects. As a result, the structures of signals that are highlighted by the representations are garbled when approaching the boundaries of the TF domain. In this paper, for the purpose of real-time TF information acquisition of nonstationary oscillatory time series, we propose a numerically efficient approach for the reduction of such boundary effects. The solution relies on an extension of the analyzed signal obtained by a forecasting technique. In the case of the study of a class of locally oscillating signals, we provide a theoretical guarantee of the performance of our approach. Following a numerical verification of the algorithmic performance of our approach, we validate it by implementing it on biomedical signals.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Convergence of Graph Laplacian with kNN Self-tuned Kernels

Kernelized Gram matrix $W$ constructed from data points $\{x_i\}_{i=1}^N$ as $W_{ij}= k_0( \frac{ \| x_i - x_j \|^2} {σ^2} )$ is widely used in graph-based geometric data analysis and unsupervised learning. An important question is how to choose the kernel bandwidth $σ$, and a common practice called self-tuned kernel adaptively sets a $σ_i$ at each point $x_i$ by the $k$-nearest neighbor (kNN) distance. When $x_i$'s are sampled from a $d$-dimensional manifold embedded in a possibly high-dimensional space, unlike with fixed-bandwidth kernels, theoretical results of graph Laplacian convergence with self-tuned kernels have been incomplete. This paper proves the convergence of graph Laplacian operator $L_N$ to manifold (weighted-)Laplacian for a new family of kNN self-tuned kernels $W^{(α)}_{ij} = k_0( \frac{ \| x_i - x_j \|^2}{ ε\hatρ(x_i) \hatρ(x_j)})/\hatρ(x_i)^α\hatρ(x_j)^α$, where $\hatρ$ is the estimated bandwidth function {by kNN}, and the limiting operator is also parametrized by $α$. When $α= 1$, the limiting operator is the weighted manifold Laplacian $Δ_p$. Specifically, we prove the point-wise convergence of $L_N f $ and convergence of the graph Dirichlet form with rates. Our analysis is based on first establishing a $C^0$ consistency for $\hatρ$ which bounds the relative estimation error $|\hatρ - \barρ|/\barρ$ uniformly with high probability, where $\barρ = p^{-1/d}$, and $p$ is the data density function. Our theoretical results reveal the advantage of self-tuned kernel over fixed-bandwidth kernel via smaller variance error in low-density regions. In the algorithm, no prior knowledge of $d$ or data density is needed. The theoretical results are supported by numerical experiments on simulated data and hand-written digit image data.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Prenatal stress perturbs fetal iron homeostasis in a sex-specific manner

What is the influence of chronic maternal prenatal stress (PS) on fetal iron homeostasis? In a prospective case-control study in 164 pregnant women, we show that cord blood transferrin saturation is lower in male stressed neonates. The total effect of PS exposure on fetal ferritin revealed a decrease of 15.4% compared with controls. Electrocardiogram-based Fetal Stress Index (FSI) identified affected fetuses non-invasively during the third trimester of gestation. FSI-based timely detection of fetuses affected by PS can support early individualized iron supplementation and neurodevelopmental follow-up to prevent long-term sequelae due to PS-exacerbated impairment of the iron homeostasis.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A persistent homology approach to heart rate variability analysis with an application to sleep-wake classification

Persistent homology (PH) is a recently developed theory in the field of algebraic topology to study shapes of datasets. It is an effective data analysis tool that is robust to noise and has been widely applied. We demonstrate a general pipeline to apply PH to study time series; particularly the instantaneous heart rate time series for the heart rate variability (HRV) analysis. The first step is capturing the shapes of time series from two different aspects -- {the PH's and hence persistence diagrams of its} sub-level set and Taken's lag map. Second, we propose a systematic {and computationally efficient} approach to summarize persistence diagrams, which we coined {\em persistence statistics}. To demonstrate our proposed method, we apply these tools to the HRV analysis and the sleep-wake, REM-NREM (rapid eyeball movement and non rapid eyeball movement) and sleep-REM-NREM classification problems. The proposed algorithm is evaluated on three different datasets via the cross-database validation scheme. The performance of our approach is better than the state-of-the-art algorithms, and the result is consistent throughout different datasets.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Airflow recovery from thoracic and abdominal movements using Synchrosqueezing Transform and Locally Stationary Gaussian Process Regression

Airflow signal encodes rich information about respiratory system. While the gold standard for measuring airflow is to use a spirometer with an occlusive seal, this is not practical for ambulatory monitoring of patients. Advances in sensor technology have made measurement of motion of the thorax and abdomen feasible with small inexpensive devices, but estimation of airflow from these time series is challenging. We propose to use the nonlinear-type time-frequency analysis tool, synchrosqueezing transform, to properly represent the thoracic and abdominal movement signals as the features, which are used to recover the airflow by the locally stationary Gaussian process. We show that, using a dataset that contains respiratory signals under normal sleep conditions, an accurate prediction can be achieved by fitting the proposed model in the feature space both in the intra- and inter-subject setups. We also apply our method to a more challenging case, where subjects under general anesthesia underwent transitions from pressure support to unassisted ventilation to further demonstrate the utility of the proposed method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An Adaptive QRS Detection Algorithm for Ultra-Long-Term ECG Recordings

Background: Accurate detection of QRS complexes during mobile, ultra-long-term ECG monitoring is challenged by instances of high heart rate, dramatic and persistent changes in signal amplitude, and intermittent deformations in signal quality that arise due to subject motion, background noise, and misplacement of the ECG electrodes. Purpose: We propose a revised QRS detection algorithm which addresses the above-mentioned challenges. Methods and Results: Our proposed algorithm is based on a state-of-the-art algorithm after applying two key modifications. The first modification is implementing local estimates for the amplitude of the signal. The second modification is a mechanism by which the algorithm becomes adaptive to changes in heart rate. We validated our proposed algorithm against the state-of-the-art algorithm using short-term ECG recordings from eleven annotated databases available at Physionet, as well as four ultra-long-term (14-day) ECG recordings which were visually annotated at a central ECG core laboratory. On the database of ultra-long-term ECG recordings, our proposed algorithm showed a sensitivity of 99.90% and a positive predictive value of 99.73%. Meanwhile, the state-of-the-art QRS detection algorithm achieved a sensitivity of 99.30% and a positive predictive value of 99.68% on the same database. The numerical efficiency of our new algorithm was evident, as a 14-day recording sampled at 200 Hz was analyzed in approximately 157 seconds. Conclusions: We developed a new QRS detection algorithm. The efficiency and accuracy of our algorithm makes it a good fit for mobile health applications, ultra-long-term and pathological ECG recordings, and the batch processing of large ECG databases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Convergence analysis of Adaptive Locally Iterative Filtering and SIFT method

Adaptive Local Iterative Filtering (ALIF) is a currently proposed novel time-frequency analysis tool. It has been empirically shown that ALIF is able to separate components and overcome the mode-mixing problem. However, so far its convergence is still an open problem, particularly for highly nonstationary signals, due to the fact that the kernel associated with ALIF is non-translational invariant, non-convolutional and non-symmetric. Our first contribution in this work is providing a convergence analysis of ALIF. From the practical perspective, ALIF depends on a robust frequencies estimator, based on which the decomposition can be achieved. Our second contribution is proposing a robust and adaptive decomposition method for noisy and nonstationary signals, which we coined the Synchrosqueezing Iterative Filtering Technique (SIFT). In SIFT, we apply the synchrosqueezing transform to estimate the instantaneous frequency, and then apply the ALIF to decompose a signal. We show numerically the ability of this new approach in handling highly nonstationary signals.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical computation of triangular complex spherical designs with small mesh ratio

This paper provides triangular spherical designs for the complex unit sphere $Ω^d$ by exploiting the natural correspondence between the complex unit sphere in $d$ dimensions and the real unit sphere in $2d-1$. The existence of triangular and square complex spherical $t$-designs with the optimal order number of points is established. A variational characterization of triangular complex designs is provided, with particular emphasis on numerical computation of efficient triangular complex designs with good geometric properties as measured by their mesh ratio. We give numerical examples of triangular spherical $t$-designs on complex unit spheres of dimension $d=2$ to $6$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the behavior of $1$-Laplacian Ratio Cuts on nearly rectangular domains

Given a connected set $Ω_0 \subset \mathbb{R}^2$, define a sequence of sets $(Ω_n)_{n=0}^{\infty}$ where $Ω_{n+1}$ is the subset of $Ω_n$ where the first eigenfunction of the (properly normalized) Neumann $p-$Laplacian $ -Δ^{(p)} ϕ= λ_1 |ϕ|^{p-2} ϕ$ is positive (or negative). For $p=1$, this is also referred to as the Ratio Cut of the domain. We conjecture that, unless $Ω_0$ is an isosceles right triangle, these sets converge to the set of rectangles with eccentricity bounded by 2 in the Gromov-Hausdorff distance as long as they have a certain distance to the boundary $\partial Ω_0$. We establish some aspects of this conjecture for $p=1$ where we prove that (1) the 1-Laplacian spectral cut of domains sufficiently close to rectangles of a given aspect ratio is a circular arc that is closer to flat than the original domain (leading eventually to quadrilaterals) and (2) quadrilaterals close to a rectangle of aspect ratio $2$ stay close to quadrilaterals and move closer to rectangles in a suitable metric. We also discuss some numerical aspects and pose many open questions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the spectral property of kernel-based sensor fusion algorithms of high dimensional data

We apply local laws of random matrices and free probability theory to study the spectral properties of two kernel-based sensor fusion algorithms, nonparametric canonical correlation analysis (NCCA) and alternating diffusion (AD), for two simultaneously recorded high dimensional datasets under the null hypothesis. The matrix of interest is the product of the kernel matrices associated with the databsets, which may not be diagonalizable in general. We prove that in the regime where dimensions of both random vectors are comparable to the sample size, if NCCA and AD are conducted using a smooth kernel function, then the first few nontrivial eigenvalues will converge to real deterministic values provided the datasets are independent Gaussian random vectors. Toward the claimed result, we also provide a convergence rate of eigenvalues of a kernel affinity matrix.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Can a composite heart rate variability biomarker shed new insights about autism spectrum disorder in school-aged children?

High-frequency heart rate variability (HRV) has identified parasympathetic nervous system alterations in autism spectrum disorder (ASD). In a cohort of school-aged children with and without ASD, we test a set of alternative linear and nonlinear HRV measures, including phase rectified signal averaging, applied to a segment of resting ECG, for associations with ASD vs. other psychiatric conditions. Using machine learning, we identify HRV measures derived from time, frequency, and geometric signal-analytical domains that (1) identify children with ASD relative to peers with receiver operating curve area of .89, and (2) differentiate such children from those with conduct problems or depression. Despite the small cohort and lack of prospective external validation, these preliminary results warrant larger prospective validation studies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Carrier frequencies, holomorphy and unwinding

We prove that functions of intrinsic-mode type (a classical models for signals) behave essentially like holomorphic functions: adding a pure carrier frequency $e^{int}$ ensures that the anti-holomorphic part is much smaller than the holomorphic part $ \| P_{-}(f)\|_{L^2} \ll \|P_{+}(f)\|_{L^2}.$ This enables us to use techniques from complex analysis, in particular the \textit{unwinding series}. We study its stability and convergence properties and show that the unwinding series can stabilize and show that the unwinding series can provide a high resolution time-frequency representation, which is robust to noise.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Convex Optimization approach to signals with fast varying instantaneous frequency

Motivated by the limitation of analyzing oscillatory signals composed of multiple components with fast-varying instantaneous frequency, we approach the time-frequency analysis problem by optimization. Based on the proposed adaptive harmonic model, the time-frequency representation of a signal is obtained by directly minimizing a functional, which involves few properties an "ideal time-frequency representation" should satisfy, for example, the signal reconstruction and concentrative time frequency representation. FISTA (Fast Iterative Shrinkage-Thresholding Algorithm) is applied to achieve an efficient numerical approximation of the functional. We coin the algorithm as {\it Time-frequency bY COnvex OptimizatioN} (Tycoon). The numerical results confirm the potential of the Tycoon algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Extract fetal ECG from single-lead abdominal ECG by de-shape short time Fourier transform and nonlocal median

The multiple fundamental frequency detection problem and the source separation problem from a single-channel signal containing multiple oscillatory components and a nonstationary noise are both challenging tasks. To extract the fetal electrocardiogram (ECG) from a single-lead maternal abdominal ECG, we face both challenges. In this paper, we propose a novel method to extract the fetal ECG signal from the single channel maternal abdominal ECG signal, without any additional measurement. The algorithm is composed of three main ingredients. First, the maternal and fetal heart rates are estimated by the de-shape short time Fourier transform, which is a recently proposed nonlinear time-frequency analysis technique; second, the beat tracking technique is applied to accurately obtain the maternal and fetal R peaks; third, the maternal and fetal ECG waveforms are established by the nonlocal median. The algorithm is evaluated on a simulated fetal ECG signal database ({\em fecgsyn} database), and tested on two real databases with the annotation provided by experts ({\em adfecgdb} database and {\em CinC2013} database). In general, the algorithm could be applied to solve other detection and source separation problems, and reconstruct the time-varying wave-shape function of each oscillatory component.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Heart beat classification from single-lead ECG using the Synchrosqueezing Transform

The processing of ECG signal provides a wealth of information on cardiac function and overall cardiovascular health. While multi-lead ECG recordings are often necessary for a proper assessment of cardiac rhythms, they are not always available or practical, for example in fetal ECG applications. Moreover, a wide range of small non-obtrusive single-lead ECG ambulatory monitoring devices are now available, from which heart rate variability (HRV) and other health-related metrics are derived. Proper beat detection and classification of abnormal rhythms is important for reliable HRV assessment and can be challenging in single-lead ECG monitoring devices. In this manuscript, we modeled the heart rate signal as an adaptive non-harmonic model and used the newly developed synchrosqueezing transform (SST) to characterize ECG patterns. We show how the proposed model can be used to enhance heart beat detection and classification between normal and abnormal rhythms. In particular, using the Massachusetts Institute of Technology-Beth Israel Hospital (MIT-BIH) arrhythmia database and the Association for the Advancement of Medical Instrumentation (AAMI) beat classes, we trained and validated a support vector machine (SVM) classifier on a portion of the annotated beat database using the SST-derived instantaneous phase, the R-peak amplitudes and R-peak to R-peak interval durations, based on a single ECG lead. We obtained sensitivities and positive predictive values comparable to other published algorithms using multiple leads and many more features.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Minimum-latency Time-frequency Analysis Using Asymmetric Window Functions

We study the real-time dynamics retrieval from a time series via the time-frequency (TF) analysis with the minimal latency guarantee. While different from the well-known intrinsic latency definition in the filter design, a rigorous definition of intrinsic latency for different time-frequency representations (TFR) is provided, including the short time Fourier transform (STFT), synchrosqeezing transform (SST) and reassignment method (RM). To achieve the minimal latency, a systematic method is proposed to construct an asymmetric window from a well-designed symmetric one based on the concept of minimum-phase, if the window satisfies some weak conditions. We theoretically show that the TFR determined by SST with the constructed asymmetric window does have a smaller intrinsic latency. Finally, the music onset detection problem is studied to show the strength of the proposed algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Modeling the pulse signal by wave-shape function and analyzing by synchrosqueezing transform

We apply the recently developed adaptive non-harmonic model based on the wave-shape function, as well as the time-frequency analysis tool called synchrosqueezing transform (SST) to model and analyze oscillatory physiological signals. To demonstrate how the model and algorithm work, we apply them to study the pulse wave signal. By extracting features called the spectral pulse signature, {and} based on functional regression, we characterize the hemodynamics from the radial pulse wave signals recorded by the sphygmomanometer. Analysis results suggest the potential of the proposed signal processing approach to extract health-related hemodynamics features.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Rank-1 accelerated illumination recovery in scanning diffractive imaging by transparency estimation

We consider the problem of blind ptychography, that is the joint estimation of an unknown object and an illumination function from diffraction intensity measurements. In ptychography, diffraction measurements from neighboring regions of the same object are related to each other by a pairwise relationship between overlapping frames. When the illumination is well known, the relationship among frames is given by a linear projection operator. We propose a power iteration-projection algorithm that minimizes the global pairwise discrepancy among frames. We accelerate the convergence of power method by subtracting the estimated localized average transparency of the unknown object. The method is effective for weakly scattering and low contrast objects or piecewise smooth specimens.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Sleep Apnea Detection Based on Thoracic and Abdominal Movement Signals of Wearable Piezo-Electric Bands

Physiologically, the thoracic (THO) and abdominal (ABD) movement signals, captured using wearable piezo-electric bands, provide information about various types of apnea, including central sleep apnea (CSA) and obstructive sleep apnea (OSA). However, the use of piezo-electric wearables in detecting sleep apnea events has been seldom explored in the literature. This study explored the possibility of identifying sleep apnea events, including OSA and CSA, by solely analyzing {one or both the THO and ABD signals. An adaptive non-harmonic model was introduced to model the THO and ABD signals, which allows us to design features for sleep apnea events. To confirm the suitability of the extracted features, a support vector machine was applied to classify three categories -- normal and hypopnea, OSA, and CSA. According to a database of} 34 subjects, the overall classification accuracies were on average $75.9\%\pm 11.7\%$ and $73.8\%\pm 4.4\%$, respectively, based on the cross validation. When the features determined from the THO and ABD signals were combined, the overall classification accuracy became $81.8\%\pm 9.4\%$. These features were applied for designing a state machine for online apnea event detection. Two event-by-event accuracy indices, S and I, were proposed for evaluating the performance {of the state machine. For the same database, the} S index was $84.01\%\pm 9.06\%$, and the I index was $77.21\%\pm 19.01\%$. The results indicate the considerable potential of applying the proposed algorithm to clinical examinations for both screening and homecare purposes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Vector Nonlocal Euclidean Median: Principal Bundle Captures The Nature of Patch Space

We extensively study the rotational group structure inside the patch space by introducing the fiber bundle structure. The rotational group structure leads to a new image denoising algorithm called the \textit{vector non-local Euclidean median} (VNLEM). The theoretical aspect of VNLEM is studied, which explains why the VNLEM and traditional non-local mean/non-local Euclidean median (NLEM) algorithm work. The numerical issue of the VNLEM is improved by taking the orientation feature in the commonly applied scale-invariant feature transform (SIFT), and a theoretical analysis of the robustness of the orientation feature in the SIFT is provided. The VNLEM is applied to an image database of 1,361 images and compared with the NLEM. Different image quality assessments based on the error-sensitivity or the human visual system are applied to evaluate the performance. The results confirmed the potential of the VNLEM algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Wave-shape function analysis -- when cepstrum meets time-frequency analysis

We propose to combine cepstrum and nonlinear time-frequency (TF) analysis to study mutiple component oscillatory signals with time-varying frequency and amplitude and with time-varying non-sinusoidal oscillatory pattern. The concept of cepstrum is applied to eliminate the wave-shape function influence on the TF analysis, and we propose a new algorithm, named de-shape synchrosqueezing transform (de-shape SST). The mathematical model, adaptive non-harmonic model, is introduced and the de-shape SST algorithm is theoretically analyzed. In addition to simulated signals, several different physiological, musical and biological signals are analyzed to illustrate the proposed algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spectral Convergence of the connection Laplacian from random samples

Spectral methods that are based on eigenvectors and eigenvalues of discrete graph Laplacians, such as Diffusion Maps and Laplacian Eigenmaps are often used for manifold learning and non-linear dimensionality reduction. It was previously shown by Belkin and Niyogi \cite{belkin_niyogi:2007} that the eigenvectors and eigenvalues of the graph Laplacian converge to the eigenfunctions and eigenvalues of the Laplace-Beltrami operator of the manifold in the limit of infinitely many data points sampled independently from the uniform distribution over the manifold. Recently, we introduced Vector Diffusion Maps and showed that the connection Laplacian of the tangent bundle of the manifold can be approximated from random samples. In this paper, we present a unified framework for approximating other connection Laplacians over the manifold by considering its principle bundle structure. We prove that the eigenvectors and eigenvalues of these Laplacians converge in the limit of infinitely many independent random samples. We generalize the spectral convergence results to the case where the data points are sampled from a non-uniform distribution, and for manifolds with and without boundary.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

When interpolation-induced reflection artifact meets time-frequency analysis

While extracting the temporal dynamical features based on the time-frequency analyses, like the reassignment and synchrosqueezing transform, attracts more and more interest in bio-medical data analysis, we should be careful about artifacts generated by interpolation schemes, in particular when the sampling rate is not significantly higher than the frequency of the oscillatory component we are interested in. In this study, we formulate the problem called the reflection effect and provide a theoretical justification of the statement. We also show examples in the anesthetic depth analysis with clear but undesirable artifacts. The results show that the artifact associated with the reflection effect exists not only theoretically but practically. Its influence is pronounced when we apply the time-frequency analyses to extract the time-varying dynamics hidden inside the signal. In conclusion, we have to carefully deal with the artifact associated with the reflection effect by choosing a proper interpolation scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Assess Sleep Stage by Modern Signal Processing Techniques

In this paper, two modern adaptive signal processing techniques, Empirical Intrinsic Geometry and Synchrosqueezing transform, are applied to quantify different dynamical features of the respiratory and electroencephalographic signals. We show that the proposed features are theoretically rigorously supported, as well as capture the sleep information hidden inside the signals. The features are used as input to multiclass support vector machines with the radial basis function to automatically classify sleep stages. The effectiveness of the classification based on the proposed features is shown to be comparable to human expert classification -- the proposed classification of awake, REM, N1, N2 and N3 sleeping stages based on the respiratory signal (resp. respiratory and EEG signals) has the overall accuracy $81.7\%$ (resp. $89.3\%$) in the relatively normal subject group. In addition, by examining the combination of the respiratory signal with the electroencephalographic signal, we conclude that the respiratory signal consists of ample sleep information, which supplements to the information stored in the electroencephalographic signal.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Connection graph Laplacian methods can be made robust to noise

Recently, several data analytic techniques based on connection graph laplacian (CGL) ideas have appeared in the literature. At this point, the properties of these methods are starting to be understood in the setting where the data is observed without noise. We study the impact of additive noise on these methods, and show that they are remarkably robust. As a by-product of our analysis, we propose modifications of the standard algorithms that increase their robustness to noise. We illustrate our results in numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Graph connection Laplacian and random matrices with random blocks

Graph connection Laplacian (GCL) is a modern data analysis technique that is starting to be applied for the analysis of high dimensional and massive datasets. Motivated by this technique, we study matrices that are akin to the ones appearing in the null case of GCL, i.e the case where there is no structure in the dataset under investigation. Developing this understanding is important in making sense of the output of the algorithms based on GCL. We hence develop a theory explaining the behavior of the spectral distribution of a large class of random matrices, in particular random matrices with random block entries of fixed size. Part of the theory covers the case where there is significant dependence between the blocks. Numerical work shows that the agreement between our theoretical predictions and numerical simulations is generally very good.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Real-time dynamics acquisition from irregular samples -- with application to anesthesia evaluation

The first objective of this paper is to introduce a unified approach to the D/A conversion, a real-time algorithm referred to as {\it blending operator}, based on spline functions of arbitrarily desired order, to interpolate the irregular data samples, while preserving all polynomials of the same spline order, with assured maximum order of approximation. This helps remove the two main obstacles for adapting the recently proposed time-frequency analysis technique {\it Synchrosqueezing transform} (SST) to irregular data samples in order to allow online computation. Secondly, for real-time dynamic information extraction from an oscillatory signal via SST, a family of vanishing-moment and minimum-supported spline-wavelets (to be called VM wavelets) are introduced for on-line computation of the CWT and its derivative. The second objective of this paper is to apply the proposed real-time algorithm and VM wavelets to clinical applications, particularly to the study of the "anesthetic depth" of a patient during surgery, with emphasis on analyzing two dynamic quantities: the "instantaneous frequencies" and the "non-rhythmic to rhythmic ratios" of the patient's respiration, based on a one-lead electrocardiogram (ECG) signal.It is envisioned that the proposed algorithm and VM wavelets should enable real-time monitoring of "anesthetic depth", during surgery, from the respiration signal via ECG measurement.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A New Test for One-Way ANOVA with Functional Data and Application to Ischemic Heart Screening

We propose and study a new global test, namely the $F_{\max}$-test, for the one-way ANOVA problem in functional data analysis. The test statistic is taken as the maximum value of the usual pointwise $F$-test statistics over the interval the functional responses are observed. A nonparametric bootstrap method is employed to approximate the null distribution of the test statistic and to obtain an estimated critical value for the test. The asymptotic random expression of the test statistic is derived and the asymptotic power is studied. In particular, under mild conditions, the $F_{\max}$-test asymptotically has the correct level and is root-$n$ consistent in detecting local alternatives. Via some simulation studies, it is found that in terms of both level accuracy and power, the $F_{\max}$-test outperforms the Globalized Pointwise F (GPF) test of \cite{Zhang_Liang:2013} when the functional data are highly or moderately correlated, and its performance is comparable with the latter otherwise. An application to an ischemic heart real dataset suggests that, after proper manipulation, resting electrocardiogram (ECG) signals can be used as an effective tool in clinical ischemic heart screening, without the need of further stress tests as in the current standard procedure.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Augmented projections for ptychographic imaging

Ptychography is a popular technique to achieve diffraction limited resolution images of a two or three dimensional sample using high frame rate detectors. We introduce a relaxation of common projection algorithms to account for instabilities given by intensity and background fluctuations, position errors, or poor calibration using multiplexing illumination. This relaxation introduces an additional phasing optimization at every step that enhances the convergence rate of common projection algorithms. Numerical tests exhibit the exact recovery of the object and the noise when there is high redundancy in the data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonparametric and adaptive modeling of dynamic seasonality and trend with heteroscedastic and dependent errors

Seasonality (or periodicity) and trend are features describing an observed sequence, and extracting these features is an important issue in many scientific fields. However, it is not an easy task for existing methods to analyze simultaneously the trend and {\it dynamics} of the seasonality such as time-varying frequency and amplitude, and the {\it adaptivity} of the analysis to such dynamics and robustness to heteroscedastic, dependent errors is not guaranteed. These tasks become even more challenging when there exist multiple seasonal components. We propose a nonparametric model to describe the dynamics of multi-component seasonality, and investigate the recently developed Synchrosqueezing transform (SST) in extracting these features in the presence of a trend and heteroscedastic, dependent errors. The identifiability problem of the nonparametric seasonality model is studied, and the adaptivity and robustness properties of the SST are theoretically justified in both discrete- and continuous-time settings. Consequently we have a new technique for de-coupling the trend, seasonality and heteroscedastic, dependent error process in a general nonparametric setup. Results of a series of simulations are provided, and the incidence time series of varicella and herpes zoster in Taiwan and respiratory signals observed from a sleep study are analyzed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Instantaneous frequency and wave shape functions (I)

Although one can formulate an intuitive notion of instantaneous frequency, generalizing "frequency" as we understand it in e.g. the Fourier transform, a rigorous mathematical definition is lacking. In this paper, we consider a class of functions composed of waveforms that repeat nearly periodically, and for which the instantaneous frequency can be given a rigorous meaning. We show that Synchrosqueezing can be used to determine the instantaneous frequency of functions in this class, even if the waveform is not harmonic, thus generalizing earlier results for cosine wave functions. We also provide real-life examples and discuss the advantages, for these examples, of considering such non-harmonic waveforms.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Local Linear Regression on Manifolds and its Geometric Interpretation

High-dimensional data analysis has been an active area, and the main focuses have been variable selection and dimension reduction. In practice, it occurs often that the variables are located on an unknown, lower-dimensional nonlinear manifold. Under this manifold assumption, one purpose of this paper is regression and gradient estimation on the manifold, and another is developing a new tool for manifold learning. To the first aim, we suggest directly reducing the dimensionality to the intrinsic dimension $d$ of the manifold, and performing the popular local linear regression (LLR) on a tangent plane estimate. An immediate consequence is a dramatic reduction in the computation time when the ambient space dimension $p\gg d$. We provide rigorous theoretical justification of the convergence of the proposed regression and gradient estimators by carefully analyzing the curvature, boundary, and non-uniform sampling effects. A bandwidth selector that can handle heteroscedastic errors is proposed. To the second aim, we analyze carefully the behavior of our regression estimator both in the interior and near the boundary of the manifold, and make explicit its relationship with manifold learning, in particular estimating the Laplace-Beltrami operator of the manifold. In this context, we also make clear that it is important to use a smaller bandwidth in the tangent plane estimation than in the LLR. Simulation studies and the Isomap face data example are used to illustrate the computational speed and estimation accuracy of our methods.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Synchrosqueezing-based Recovery of Instantaneous Frequency from Nonuniform Samples

We propose a new approach for studying the notion of the instantaneous frequency of a signal. We build on ideas from the Synchrosqueezing theory of Daubechies, Lu and Wu and consider a variant of Synchrosqueezing, based on the short-time Fourier transform, to precisely define the instantaneous frequencies of a multi-component AM-FM signal. We describe an algorithm to recover these instantaneous frequencies from the uniform or nonuniform samples of the signal and show that our method is robust to noise. We also consider an alternative approach based on the conventional, Hilbert transform-based notion of instantaneous frequency to compare to our new method. We use these methods on several test cases and apply our results to a signal analysis problem in electrocardiography.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Synchrosqueezing algorithm for time-varying spectral analysis: robustness properties and new paleoclimate applications

We analyze the stability properties of the Synchrosqueezing transform, a time-frequency signal analysis method that can identify and extract oscillatory components with time-varying frequency and amplitude. We show that Synchrosqueezing is robust to bounded perturbations of the signal and to Gaussian white noise. These results justify its applicability to noisy or nonuniformly sampled data that is ubiquitous in engineering and the natural sciences. We also describe a practical implementation of Synchrosqueezing and provide guidance on tuning its main parameters. As a case study in the geosciences, we examine characteristics of a key paleoclimate change in the last 2.5 million years, where Synchrosqueezing provides significantly improved insights.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Vector Diffusion Maps and the Connection Laplacian

We introduce {\em vector diffusion maps} (VDM), a new mathematical framework for organizing and analyzing massive high dimensional data sets, images and shapes. VDM is a mathematical and algorithmic generalization of diffusion maps and other non-linear dimensionality reduction methods, such as LLE, ISOMAP and Laplacian eigenmaps. While existing methods are either directly or indirectly related to the heat kernel for functions over the data, VDM is based on the heat kernel for vector fields. VDM provides tools for organizing complex data sets, embedding them in a low dimensional space, and interpolating and regressing vector fields over the data. In particular, it equips the data with a metric, which we refer to as the {\em vector diffusion distance}. In the manifold learning setup, where the data set is distributed on (or near) a low dimensional manifold $\MM^d$ embedded in $\RR^{p}$, we prove the relation between VDM and the connection-Laplacian operator for vector fields over the manifold.

0 citations0 reviewsNo code linkedOpen access