Functional Spectral Imaging by Ultrasound (FSIU): A Spectral-Theoretic Basis for Functional Ultrasound

TL;DR

This work addresses what tissue properties can be inferred from spectral data and how to recover them without ionizing radiation. It develops Functional Spectral Imaging (FSI) based on eigenvalue perturbations of a self-adjoint elliptic operator L = -D(x)∇^2 + γ(x), with first-order Hadamard formulas δλ_i = ∫Ω δγ(x)(ψ_i^0)^2 dx − ∫Ω (∇ψ_i^0)^T δD(x) ∇ψ_i^0 dx and a Fréchet-derivative framework D F_(0,0)[δD, δγ], enabling variational inversion with explicit gradients. The study demonstrates in silico that retaining ~10–15 low-order modes preserves ≈85% of anomaly contrast while suppressing noise, achieving submillimetric localization (~0.1–0.3 mm) and mg-scale detectability under ideal noise, and introduces a spectral-entropy index to differentiate compact from diffuse inclusions. Concluding, the operator-spectral FSI provides a rigorous, non-ionizing avenue for localized functional imaging with clear pathways to phantom validation and eventual in vivo studies, and suggests integration with ultrasound hardware for portable point-of-care applications.

Abstract

Functional Spectral Imaging (FSI) models image formation as the recovery of tissue surrogates such as density and stiffness from spectral perturbations of a self-adjoint elliptic operator. Rather than relying on reflectivity or relaxation kinetics, FSI tracks shifts of a truncated set of eigenmodes under controlled excitation, providing a non-ionizing and operator-theoretic route to contrast. Tissue heterogeneity is modeled as a small perturbation of L = -div(D grad) + gamma, with first-order Hadamard formulas linking local contrasts to eigenvalue shifts. Frechet derivatives and their adjoints yield gradients for variational inversion, stabilized by Tikhonov or total-variation regularization and modal truncation. Finite-element simulations show submillimetric localization (about 0.1-0.3 mm) and milligram-scale detectability (thresholds near 1 mg) under ideal noise. Retaining 10-15 modes preserves about 85 percent of anomaly contrast while suppressing noise. A spectral-entropy index separates compact from diffuse inclusions and acts as a morphology surrogate. FSI thus provides a mathematically controlled, non-ionizing framework for localized functional imaging, motivating validation in physical phantoms and in vivo studies.

Figures (7)

• Figure 1: Conceptual scheme of the domain $\Omega$ with Dirichlet boundary conditions and operator $L$. The eigenfunctions $\psi_i$ form an orthogonal basis of $L^2(\Omega)$.
• Figure 2: Eigenvalue spectrum schematic. Well-separated low-order modes are robust and informative, whereas high-order modes become dense and noise-sensitive (Weyl asymptotics).
• Figure 3: Eigenvalue spectrum trade-off and modal robustness. Low-order modes are robust yet less sensitive; high-order modes become noise-sensitive as spectral gaps shrink. Insets: a compact, low-sensitivity mode (left) and a highly oscillatory, noise-sensitive mode (right).
• Figure 4: Modal contribution spectrum for a synthetic anomaly. Each bar displays the weighted perturbation $\alpha_i\delta\lambda_i$ for mode $i$. Dominant contributions arise from $i=1$–12, consistent with the SNR-based truncation rule.
• Figure 5: Resolution benchmark from controlled numerical simulations (*in silico*). By tracking eigenmode perturbations, FSI resolves inclusions smaller than $0.3$ mm and localizes them within $0.15$ mm. The mechanism rests on mode-shift sensitivity to local changes in density or stiffness, rather than on a propagating wavelength scale.