Source author record

Bernhard G. Bodmann

Bernhard G. Bodmann appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

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

11 published item(s)

preprint2022arXiv

A scattering transform for graphs based on heat semigroups, with an application for the detection of anomalies in positive time series with underlying periodicities

This paper develops an adaptive version of Mallat's scattering transform for signals on graphs. The main results are norm bounds for the layers of the transform, obtained from a version of a Beurling-Deny inequality that permits to remove the nonlinear steps in the scattering transform. Under statistical assumptions on the input signal, the norm bounds can be refined. The concepts presented here are illustrated with an application to traffic counts which exhibit characteristic daily and weekly periodicities. Anomalous traffic patterns which deviate from these expected periodicities produce a response in the scattering transform.

preprint2022arXiv

Spikes, Roots, and Modulations: Phase Retrieval for Finitely-Supported Complex Measures

We study the recovery of a finitely supported distribution, a complex linear combination of Dirac measures, from intensity measurements. The distribution $μ=\sum_{j=1}^{s}c_{j}δ_{t_{j}}$ is given by a coefficient vector $c\in\mathbb{C}^s$ and its support $\{t_1, t_2, \dots, t_s\}$ is contained in $ [0,Λ]$ for some $Λ>0$. The intensity measurements evaluate (squared) magnitudes of a set of linear functionals applied to $μ$, obtained by sampling $\hat μ$, the Fourier transform of $μ$, or by evaluating differences between modulated samples. Following a strategy by Alexeev et al., the structure of the linear functionals, and hence of the non-linear magnitude measurement, is encoded with a graph, where the vertices represent point evaluations of $\hat μ$ at $\{v_1, v_2, \dots, v_n\} \subset [-Ω,Ω]$ and each edge represents a (modulated) difference between vertices incident with it. We show that a Ramanujan graph with degree $d \ge 3$ and $n>\frac{6(1 + 6 /\ln(s/ΛΩ)) s}{1-2\sqrt{d-1}/d}$ vertices provides $M=(d+1)n$ magnitudes that are sufficient for identifying the complex measure up to an overall unimodular multiplicative constant. At the cost of including an additional oversampling step and with an additional requirement that $n-1$ is prime, we construct an explicit recovery algorithm that is based on the Prony method.

preprint2016arXiv

Maximal Orthoplectic Fusion Frames from Mutually Unbiased Bases and Block Designs

The construction of optimal line packings in real or complex Euclidean spaces has shown to be a tantalizingly difficult task, because it includes the problem of finding maximal sets of equiangular lines. In the regime where equiangular lines are not possible, some optimal packings are known, for example, those achieving the orthoplex bound related to maximal sets of mutually unbiased bases. In this paper, we investigate the packing of subspaces instead of lines and determine the implications of maximality in this context. We leverage the existence of real or complex maximal mutually unbiased bases with a combinatorial design strategy in order to find optimal subspace packings that achieve the orthoplex bound. We also show that maximal sets of mutually unbiased bases convert between coordinate projections associated with certain balanced incomplete block designs and Grassmannian 2-designs. Examples of maximal orthoplectic fusion frames already appeared in the works by Shor, Sloane and by Zauner. They are realized in dimensions that are a power of four in the real case or a power of two in the complex case.

preprint2016arXiv

Region-of-Interest reconstruction from truncated cone-beam projections

Region-of-Interest (ROI) tomography aims at reconstructing a region of interest $C$ inside a body using only x-ray projections intersecting $C$ with the goal to reduce overall radiation exposure when only a small specific region of the body needs to be examined. We consider x-ray acquisition from sources located on a smooth curve $Γ$ in $\mathbb{R}^3$ verifying classical Tuy's condition. In this situation, the {\it non-trucated} cone-beam transform $D f$ of smooth densities $f$ admits an explicit inverse $Z$; however $Z$ cannot directly reconstruct $f$ from ROI-truncated projections. To deal with the ROI tomography problem, we introduce a novel reconstruction approach. For densities $f$ in $L^{\infty}(B)$ where $B$ is a bounded ball in $\mathbb{R}^3$, our method iterates an operator $U$ combining ROI-truncated projections, inversion by the operator $Z$ and appropriate regularization operators. Assuming only knowledge of projections corresponding to a spherical ROI $C \subset B$, given $ε>0$, we prove that if $C$ is sufficiently large our iterative reconstruction algorithm converges uniformly to an $ε$-accurate approximation of $f$, where the accuracy depends on the regularity of $f$ quantified in the Sobolev norm $W^5(B)$. This result shows the existence of a critical ROI radius ensuring the convergence of the ROI reconstruction algorithm to $ε$-accurate approximations of $f$. We numerically verified these theoretical results using simulated acquisition of ROI-truncated cone-beam projection data for multiple acquisition geometries. Numerical experiments indicate that the critical ROI radius is fairly small with respect to the support region~$B$.

preprint2015arXiv

Achieving the orthoplex bound and constructing weighted complex projective 2-designs with Singer sets

Equiangular tight frames are examples of Grassmannian line packings for a Hilbert space. More specifically, according to a bound by Welch, they are minimizers for the maximal magnitude occurring among the inner products of all pairs of vectors in a unit-norm frame. This paper is dedicated to packings in the regime in which the number of frame vectors precludes the existence of equiangular frames. The orthoplex bound then serves as an alternative to infer a geometric structure of optimal designs. We construct frames of unit-norm vectors in $K$-dimensional complex Hilbert spaces that achieve the orthoplex bound. When $K-1$ is a prime power, we obtain a tight frame with $K^2+1$ vectors and when $K$ is a prime power, with $K^2+K-1$ vectors. In addition, we show that these frames form weighted complex projective 2-designs that are useful additions to maximal equiangular tight frames and maximal sets of mutually unbiased bases in quantum state tomography. Our construction is based on Singer's family of difference sets and the related concept of relative difference sets.

preprint2014arXiv

Algorithms and error bounds for noisy phase retrieval with low-redundancy frames

The main objective of this paper is to find algorithms accompanied by explicit error bounds for phase retrieval from noisy magnitudes of frame coefficients when the underlying frame has a low redundancy. We achieve these goals with frames consisting of $N=6d-3$ vectors spanning a $d$-dimensional complex Hilbert space. The two algorithms we use, phase propagation or the kernel method, are polynomial time in the dimension $d$. To ensure a successful approximate recovery, we assume that the noise is sufficiently small compared to the squared norm of the vector to be recovered. In this regime, the error bound is inverse proportional to the signal-to-noise ratio. Upper and lower bounds on the sample values of trigonometric polynomials are a central technique in our error estimates.

preprint2014arXiv

Frame potentials and the geometry of frames

This paper concerns the geometric structure of optimizers for frame potentials. We consider finite, real or complex frames and rotation or unitarily invariant potentials, and mostly specialize to Parseval frames, meaning the frame potential to be optimized is a function on the manifold of Gram matrices belonging to finite Parseval frames. Next to the known classes of equal-norm and equiangular Parseval frames, we introduce equidistributed Parseval frames, which are more general than the equiangular type but have more structure than equal-norm ones. We also provide examples where this class coincides with that of Grassmannian frames, the minimizers for the maximal magnitude among inner products between frame vectors. These different types of frames are characterized in relation to the optimization of frame potentials. Based on results by Lojasiewicz, we show that the gradient descent for a real analytic frame potential on the manifold of Gram matrices belonging to Parseval frames always converges to a critical point. We then derive geometric structures associated with the critical points of different choices of frame potentials. The optimal frames for families of such potentials are thus shown to be equal-norm, or additionally equipartitioned, or even equidistributed.

preprint2013arXiv

Gabor Shearlets

In this paper, we introduce Gabor shearlets, a variant of shearlet systems, which are based on a different group representation than previous shearlet constructions: they combine elements from Gabor and wavelet frames in their construction. As a consequence, they can be implemented with standard filters from wavelet theory in combination with standard Gabor windows. Unlike the usual shearlets, the new construction can achieve a redundancy as close to one as desired. Our construction follows the general strategy for shearlets. First we define group-based Gabor shearlets and then modify them to a cone-adapted version. In combination with Meyer filters, the cone-adapted Gabor shearlets constitute a tight frame and provide low-redundancy sparse approximations of the common model class of anisotropic features which are cartoon-like functions.

preprint2013arXiv

Random fusion frames are nearly equiangular and tight

This paper demonstrates that random, independently chosen equi-dimensional subspaces with a unitarily invariant distribution in a real Hilbert space provide nearly tight, nearly equiangular fusion frames. The angle between a pair of subspaces is measured in terms of the Hilbert-Schmidt inner product of the corresponding orthogonal projections. If the subspaces are selected at random, then a measure concentration argument shows that these inner products concentrate near an average value. Overwhelming success probability for near tightness and equiangularity is guaranteed if the dimension of the subspaces is sufficiently small compared to that of the Hilbert space and if the dimension of the Hilbert space is small compared to the sum of all subspace dimensions.

preprint2013arXiv

Stable phase retrieval with low-redundancy frames

We investigate the recovery of vectors from magnitudes of frame coefficients when the frames have a low redundancy, meaning a small number of frame vectors compared to the dimension of the Hilbert space. We first show that for vectors in d dimensions, 4d-4 suitably chosen frame vectors are sufficient to uniquely determine each signal, up to an overall unimodular constant, from the magnitudes of its frame coefficients. Then we discuss the effect of noise and show that 8d-4 frame vectors provide a stable recovery if part of the frame coefficients is bounded away from zero. In this regime, perturbing the magnitudes of the frame coefficients by noise that is sufficiently small results in a recovery error that is at most proportional to the noise level.

preprint2010arXiv

Spanning and independence properties of frame partitions

We answer a number of open problems in frame theory concerning the decomposition of frames into linearly independent and/or spanning sets. We prove that in finite dimensional Hilbert spaces, Parseval frames with norms bounded away from 1 can be decomposed into a number of sets whose complements are spanning, where the number of these sets only depends on the norm bound. We also prove, assuming the Kadison-Singer conjecture is true, that this holds for infinite dimensional Hilbert spaces. Further, we prove a stronger result for Parseval frames whose norms are uniformly small, which shows that in addition to the spanning property, the sets can be chosen to be independent, and the complement of each set to contain a number of disjoint, spanning sets.