Source-linked AI summary

Data-Driven Sparse Sensor Placement for Reconstruction

Krithika Manohar, Bingni W. Brunton, J. Nathan Kutz, Steven L. Brunton

arXiv:1701.07569v2math.OCeess.SY

TL;DR

The paper addresses how to place few sensors for reconstructing high-dimensional signals when useful structure is available in the data. It develops data-driven sensing with tailored features and compares it with universal compressed sensing, finding more accurate reconstruction for optimized tailored sensing while recognizing compressed sensing’s wider signal coverage.

  • Problem

    The paper asks how to design a small number of optimized measurements for reconstructing high-dimensional signals when sensor cost and estimation latency matter.

  • Method

    The paper extracts tailored low-rank features from training data and uses sparse sensing algorithms, including QR pivoting, to select point sensors.

  • Results

    Optimized sensing in a tailored basis typically provides more accurate signal reconstruction than compressed sensing.

  • Takeaways & Limitations

    When signal type information is available, principled sensor selection can reduce the required number of sensors and downstream computation.

  • Takeaways & Limitations

    Compressed sensing provides more flexibility when the signal is unknown, whereas tailored sensing relies on available information about the signal type.

Abstract

from arXiv · show

Optimal sensor placement is a central challenge in the design, prediction, estimation, and control of high-dimensional systems. High-dimensional states can often leverage a latent low-dimensional representation, and this inherent compressibility enables sparse sensing. This article explores optimized sensor placement for signal reconstruction based on a tailored library of features extracted from training data. Sparse point sensors are discovered using the singular value decomposition and QR pivoting, which are two ubiquitous matrix computations that underpin modern linear dimensionality reduction. Sparse sensing in a tailored basis is contrasted with compressed sensing, a universal signal recovery method in which an unknown signal is reconstructed via a sparse representation in a universal basis. Although compressed sensing can recover a wider class of signals, we demonstrate the benefits of exploiting known patterns in data with optimized sensing. In particular, drastic reductions in the required number of sensors and improved reconstruction are observed in examples ranging from facial images to fluid vorticity fields. Principled sensor placement may be critically enabling when sensors are costly and provides faster state estimation for low-latency, high-bandwidth control. MATLAB code is provided for all examples.

Compressed sensing: Random measurements in a universal basis

Compressed sensing reconstructs signals from low-dimensional measurements by exploiting sparsity in a universal basis. Random measurements and sparse recovery can reduce sampling requirements, but reconstruction depends on signal sparsity and measurement conditions.

  • Signal sparsity: Signals such as images and audio are highly compressible because only a few Fourier or wavelet modes may carry most information.This sparsity underlies compression methods and allows accurate reconstruction from a reduced representation.
  • Measurement and recovery: Compressed sensing collects low-dimensional measurements and infers sparse coefficients instead of first acquiring and discarding a high-dimensional signal.The measurements are used to solve for a sparse coefficient vector in a universal basis.
  • Measurement and recovery: Recovery seeks the sparsest solution to an underdetermined linear system, typically using ℓ1 minimization under conditions such as incoherence and the restricted isometry property.These conditions help preserve the geometry of sparse vectors and support recovery with high probability.
  • Contrast with tailored sensing: Data-driven sensing can use as few optimized samples as intrinsic rank r, often much smaller than sparsity K, but it applies to signals represented by known tailored patterns.Tailored sensor placement uses training-data features and can reduce sensors and downstream computation for structured physical signals.
  • Measurement conditions: Compressed sensing generally requires measurements scaling with K log(n/K), where K is sparsity, while less coherent measurements reduce the required count.Random projections are incoherent, but their incoherence is not optimal and random sensing can be difficult for physical systems.

Proper orthogonal decomposition

Proper orthogonal decomposition (POD) extracts a low-dimensional, data-tailored representation of high-dimensional states, enabling reconstruction and computationally efficient estimation from sparse point measurements.

  • Uses and assumptions: The POD basis must be tailored using training data, so its effectiveness depends on the available data and chosen rank.Selecting the intrinsic rank without magnifying noise is difficult, and traveling-wave data require particular care because the coefficients vary over time.
  • POD construction: The resulting coefficients can be lifted back to the full state through a linear combination of POD modes.This provides a reduced representation of x while retaining a route to full-state reconstruction.
  • POD construction: POD uses singular value decomposition to obtain orthonormal eigenmodes and low-dimensional coefficients from state snapshots.For a data matrix X, the eigenmodes are the left singular vectors, and the coefficients represent states in the POD subspace.
  • Uses and assumptions: POD reduces high-dimensional systems for numerical time-stepping, parameter estimation, and control.The reduced representation is used to streamline tasks that would otherwise be expensive in the full state space.
  • Sparse reconstruction: Sensor placement selects rows of a tailored basis so Θ = CΨr is well-conditioned for recovering low-rank coefficients.QR pivoting provides an approximate greedy selection of point sensors, including oversampled cases with p > r.
  • Sparse reconstruction: Selected point observations can reconstruct the full state by estimating POD coefficients with the pseudoinverse of the measurement-basis operator.The reconstruction is expressed as a = Θ†y, where y contains the selected observations.

Comparison of methods

The paper compares data-driven QR sensor placement in tailored bases with compressed sensing, interpolation alternatives, and random sampling across fluid, image, and environmental examples. Tailored QR sensors exploit low-rank structure to achieve accurate reconstruction with fewer sensors, while performance depends on sufficient training data and an appropriate basis rank.

  • Method comparison: Across examples, tailored QR sensors outperform random sensors and can reduce the number of sensors required for a given reconstruction performance.The reported benefits include lower sensor cost and reduced computation latency when training data sufficiently represent the reconstruction task.
  • Other tailored bases: For polynomial interpolation, QR-selected pivots overcome the Gibbs phenomenon associated with ill-conditioned equispaced points.The tailored basis is the Vandermonde matrix of polynomial features, with QR pivots supplying the interpolation points.
  • Flow past a cylinder: Using modes beyond r > 42 causes overfitting because QR pivoting begins selecting sensors associated with uninformative modes.The reconstruction accuracy therefore stops improving once the informative POD modes have been included.
  • Comparison with DEIM: Oversampled QR with p = 2r achieves a 4x reduction in reconstruction error relative to DEIM, while QR and DEIM are comparable at p = r = 40.The figure reports that QR-based sampling remains within a constant factor of the optimal POD approximation using full states.
  • Method comparison: QR pivoting selects optimized sensors from tailored bases, whereas compressed sensing uses random measurements in a universal basis.The paper frames compressed sensing as more flexible for unknown signal structure, while data-driven sensing requires samples tied to the data’s intrinsic rank.

Discussion

The paper argues that sparse sensing has no all-purpose strategy: data-driven tailored sensing favors scoped problems with training data, while random universal sensing is more general. Across applications, tailored sensing typically improves reconstruction accuracy, reduces sensor counts, and speeds computation.

  • The study demonstrates practical implementations of several sparse sensing algorithms on relevant real-world examples.
  • No all-purpose strategy exists for sparse sensing of high-dimensional systems.
  • Data-driven tailored sensing and random universal sensing form two fundamental categories.The choice depends on training data, task scope, sensor cost, and required latency.
  • Tailored sensing produces optimized sensors for tightly scoped reconstruction when sufficient training data is available.
  • Random sensing requires more sensors but makes fewer assumptions about the underlying signal, making it more general.
  • Tailored sensing typically provides more accurate reconstruction than random measurements while reducing the number of sensors by about a factor of two.
  • Tailored-basis sensor selection is computationally efficient and simple to implement, whereas compressed sensing generally requires a costly iterative algorithm.
  • Reducing measurements to key pixels can speed computation and reduce power consumption in high-dimensional sensing tasks.

Reproducible Research

The authors provide MATLAB materials to reproduce the manuscript’s results, including datasets or links to public data and scripts for recreating result figures.

  • MATLAB code and datasets or public-data links are provided to support reproduction of the results.
  • MATLAB scripts are included to recreate the figures of results.
Loading 1701.07569v2…