Researcher profile

José Luis Romero

José Luis Romero contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
13topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

16 published item(s)

preprint2026arXiv

Model agnostic signal encoding by leaky integrate and fire, performance and uncertainty

Integrate-and-fire is a resource efficient time-encoding mechanism that summarizes into a signed spike train those time intervals where a signal's charge exceeds a certain threshold. We analyze the IF encoder in terms of a very general notion of approximate bandwidth, which is shared by most commonly-used signal models. This complements results on exact encoding that may be overly adapted to a particular signal model. We take into account, possibly for the first time, the effect of uncertainty in the exact location of the spikes (as may arise by decimation), uncertainty of integration leakage (as may arise in realistic manufacturing), and boundary effects inherent to finite periods of exposure to the measurement device. The analysis is done by means of a concrete bandwidth-based Ansatz that can also be useful to initialize more sophisticated model specific reconstruction algorithms, and uses the earth mover's (Wasserstein) distance to measure spike discrepancy.

preprint2022arXiv

Efficient computation of the zeros of the Bargmann transform under additive white noise

We study the computation of the zero set of the Bargmann transform of a signal contaminated with complex white noise, or, equivalently, the computation of the zeros of its short-time Fourier transform with Gaussian window. We introduce the adaptive minimal grid neighbors algorithm (AMN), a variant of a method that has recently appeared in the signal processing literature, and prove that with high probability it computes the desired zero set. More precisely, given samples of the Bargmann transform of a signal on a finite grid with spacing $δ$, AMN is shown to compute the desired zero set up to a factor of $δ$ in the Wasserstein error metric, with failure probability $O(δ^4 \log^2(1/δ))$. We also provide numerical tests and comparison with other algorithms.

preprint2022arXiv

Invertibility of frame operators on Besov-type decomposition spaces

We derive an extension of the Walnut-Daubechies criterion for the invertibility of frame operators. The criterion concerns general reproducing systems and Besov-type spaces. As an application, we conclude that $L^2$ frame expansions associated with smooth and fast-decaying reproducing systems on sufficiently fine lattices extend to Besov-type spaces. This simplifies and improves recent results on the existence of atomic decompositions, which only provide a particular dual reproducing system with suitable properties. In contrast, we conclude that the $L^2$ canonical frame expansions extend to many other function spaces, and, therefore, operations such as analyzing using the frame, thresholding the resulting coefficients, and then synthesizing using the canonical dual frame are bounded on these spaces.

preprint2022arXiv

Spectral-norm risk rates for multi-taper estimation of Gaussian processes

We consider the estimation of the covariance of a stationary Gaussian process on a multi-dimensional grid from observations taken on a general acquisition domain. We derive spectral-norm risk rates for multi-taper estimators. When applied to one dimensional acquisition intervals, these show that Thomson's classical multi-taper has optimal risk rates, as they match known benchmarks. We also extend existing lower risk bounds to multi-dimensional grids and conclude that multi-taper estimators associated with certain two-dimensional acquisition domains also have almost optimal risk rates.

preprint2022arXiv

The planar low temperature Coulomb gas: separation and equidistribution

We consider planar Coulomb systems consisting of a large number $n$ of repelling point charges in the low temperature regime, where the inverse temperature $β$ grows at least logarithmically in $n$ as $n \longrightarrow \infty$, i.e., $β\gtrsim \log n$. Under suitable conditions on an external potential we prove results to the effect that the gas is with high probability uniformly separated and equidistributed with respect to the corresponding equilibrium measure (in the given external field). Our results generalize earlier results about Fekete configurations, i.e., the case $β=\infty$. There are also several auxiliary results which could be of independent interest. For example, our method of proof of equidistribution (a variant of "Landau's method") works for general families of configurations which are uniformly separated and which satisfy certain sampling and interpolation inequalities.

preprint2022arXiv

Zeros of Gaussian Weyl-Heisenberg functions and hyperuniformity of charge

We study Gaussian random functions on the complex plane whose stochastics are invariant under the Weyl-Heisenberg group (twisted stationarity). The theory is modeled on translation invariant Gaussian entire functions, but allows for non-analytic examples, in which case winding numbers can be either positive or negative. We calculate the first intensity of zero sets of such functions, both when considered as points on the plane, or as charges according to their phase winding. In the latter case, charges are shown to be in a certain average equilibrium independently of the particular covariance structure (universal screening). We investigate the corresponding fluctuations, and show that in many cases they are suppressed at large scales (hyperuniformity). This means that universal screening is empirically observable at large scales. We also derive an asymptotic expression for the charge variance. As a main application, we obtain statistics for the zero sets of the short-time Fourier transform of complex white noise with general windows, and also prove the following uncertainty principle: the expected number of zeros per unit area is minimized, among all window functions, exactly by generalized Gaussians. Further applications include poly-entire functions such as covariant derivatives of Gaussian entire functions.

preprint2021arXiv

On dual molecules and convolution-dominated operators

We show that sampling or interpolation formulas in reproducing kernel Hilbert spaces can be obtained by reproducing kernels whose dual systems form molecules, ensuring that the size profile of a function is fully reflected by the size profile of its sampled values. The main tool is a local holomorphic calculus for convolution-dominated operators, valid for groups with possibly non-polynomial growth. Applied to the matrix coefficients of a group representation, our methods improve on classical results on atomic decompositions and bridge a gap between abstract and concrete methods.

preprint2021arXiv

Sampling the flow of a bandlimited function

We analyze the problem of reconstruction of a bandlimited function $f$ from the space-time samples of its states $f_t=ϕ_t\ast f$ resulting from the convolution with a kernel $ϕ_t$. It is well-known that, in natural phenomena, uniform space-time samples of $f$ are not sufficient to reconstruct $f$ in a stable way. To enable stable reconstruction, a space-time sampling with periodic nonuniformly spaced samples must be used as was shown by Lu and Vetterli. We show that the stability of reconstruction, as measured by a condition number, controls the maximal gap between the spacial samples. We provide a quantitative statement of this result. In addition, instead of irregular space-time samples, we show that uniform dynamical samples at sub-Nyquist spatial rate allow one to stably reconstruct the function $\widehat f$ away from certain, explicitly described blind spots. We also consider several classes of finite dimensional subsets of bandlimited functions in which the stable reconstruction is possible, even inside the blind spots. We obtain quantitative estimates for it using Remez-Turán type inequalities. En route, we obtain a Remez-Turán inequality for prolate spheroidal wave functions. To illustrate our results, we present some numerics and explicit estimates for the heat flow problem.

preprint2020arXiv

Balian-Low type theorems on homogeneous groups

We prove strict necessary density conditions for coherent frames and Riesz sequences on homogeneous groups. Let $N$ be a connected, simply connected nilpotent Lie group with a dilation structure (a homogeneous group) and let $(π, \mathcal{H}_π)$ be an irreducible, square-integrable representation modulo the center $Z(N)$ of $N$ on a Hilbert space $\mathcal{H}_π$ of formal dimension $d_π$. If $g \in \mathcal{H}_π$ is an integrable vector and the set $\{ π(λ)g : λ\in Λ\}$ for a discrete subset $Λ\subseteq N / Z(N)$ forms a frame for $\mathcal{H}_π$, then its density satisfies the strict inequality $D^-(Λ)> d_π$, where $D^-(Λ)$ is the lower Beurling density. An analogous density condition $D^+(Λ) < d_π$ holds for a Riesz sequence in $\mathcal{H}_π$ contained in the orbit of $(π, \mathcal{H}_π)$. The proof is based on a deformation theorem for coherent systems, a universality result for $p$-frames and $p$-Riesz sequences, some results from Banach space theory, and tools from the analysis on homogeneous groups.

preprint2020arXiv

Filtering with Wavelet Zeros and Gaussian Analytic Functions

We present the continuous wavelet transform (WT) of white Gaussian noise and establish a connection to the theory of Gaussian analytic functions. Based on this connection, we propose a methodology that detects components of a signal in white noise based on the distribution of the zeros of its continuous WT. To illustrate that the continuous theory can be employed in a discrete setting, we establish a uniform convergence result for the discretized continuous WT and apply the proposed method to a variety of acoustic signals.

preprint2020arXiv

Multitaper estimation on arbitrary domains

Multitaper estimators have enjoyed significant success in estimating spectral densities from finite samples using as tapers Slepian functions defined on the acquisition domain. Unfortunately, the numerical calculation of these Slepian tapers is only tractable for certain symmetric domains, such as rectangles or disks. In addition, no performance bounds are currently available for the mean squared error of the spectral density estimate. This situation is inadequate for applications such as cryo-electron microscopy, where noise models must be estimated from irregular domains with small sample sizes. We show that the multitaper estimator only depends on the linear space spanned by the tapers. As a result, Slepian tapers may be replaced by proxy tapers spanning the same subspace (validating the common practice of using partially converged solutions to the Slepian eigenproblem as tapers). These proxies may consequently be calculated using standard numerical algorithms for block diagonalization. We also prove a set of performance bounds for multitaper estimators on arbitrary domains. The method is demonstrated on synthetic and experimental datasets from cryo-electron microscopy, where it reduces mean squared error by a factor of two or more compared to traditional methods.

preprint2020arXiv

The Nyquist sampling rate for spiraling curves

We consider the problem of reconstructing a compactly supported function from samples of its Fourier transform taken along a spiral. We determine the Nyquist sampling rate in terms of the density of the spiral and show that below this rate spirals suffer from an approximate form of aliasing. This sets a limit to the amount of undersampling that compressible signals admit when sampled along spirals. More precisely, we derive a lower bound on the condition number for the reconstruction of functions of bounded variation, and for functions that are sparse in the Haar wavelet basis.

preprint2019arXiv

Strict density inequalities for sampling and interpolation in weighted spaces of holomorphic functions

Answering a question of Lindholm, we prove strict density inequalities for sampling and interpolation in Fock spaces of entire functions in several complex variables defined by a plurisubharmonic weight. In particular, these spaces do not admit a set that is simultaneously sampling and interpolating. To prove optimality of the density conditions, we construct sampling sets with a density arbitrarily close to the critical density. The techniques combine methods from several complex variables (estimates for $\bar \partial$) and the theory of localized frames in general reproducing kernel Hilbert spaces (with no analyticity assumed). The abstract results on Fekete points and deformation of frames may be of independent interest.

preprint2018arXiv

Sharp results on sampling with derivatives in shift-invariant spaces and multi-window Gabor frames

We study the problem of sampling with derivatives in shift-invariant spaces generated by totally-positive functions of Gaussian type or by the hyperbolic secant. We provide sharp conditions in terms of weighted Beurling densities. As a by-product we derive new results about multi-window Gabor frames with respect to vectors of Hermite functions or totally positive functions.

preprint2017arXiv

Computing reconstructions from nonuniform Fourier samples: Universality of stability barriers and stable sampling rates

We study the problem of recovering an unknown compactly-supported multivariate function from samples of its Fourier transform that are acquired nonuniformly, i.e. not necessarily on a uniform Cartesian grid. Reconstruction problems of this kind arise in various imaging applications, where Fourier samples are taken along radial lines or spirals for example. Specifically, we consider finite-dimensional reconstructions, where a limited number of samples is available, and investigate the rate of convergence of such approximate solutions and their numerical stability. We show that the proportion of Fourier samples that allow for stable approximations of a given numerical accuracy is independent of the specific sampling geometry and is therefore universal for different sampling scenarios. This allows us to relate both sufficient and necessary conditions for different sampling setups and to exploit several results that were previously available only for very specific sampling geometries. The results are obtained by developing: (i) a transference argument for different measures of the concentration of the Fourier transform and Fourier samples; (ii) frame bounds valid up to the critical sampling density, which depend explicitly on the sampling set and the spectrum. As an application, we identify sufficient and necessary conditions for stable and accurate reconstruction of algebraic polynomials or wavelet coefficients from nonuniform Fourier data.