Minimax optimal convex methods for Poisson inverse problems under q-ball sparsity

Abstract

In this paper, we study the minimax rates and provide an implementable convex algorithm for Poisson inverse problems under weak sparsity and physical constraints. In particular we assume the model yi Poisson(Taif*) for 1 ≤ i ≤ n where T ∈ R+ is the intensity, and we impose weak sparsity on f* ∈ Rp by assuming f* lies in an q-ball when rotated according to an orthonormal basis D ∈ Rp × p. In addition, since we are modeling real physical systems we also impose positivity and flux-preserving constraints on the matrix A = [a1, a2,...,an] and the function f*. We prove minimax lower bounds for this model which scale as Rq ( pT)1 - q/2 where it is noticeable that the rate depends on the intensity T and not the sample size n. We also show that a 1-based regularized least-squares estimator achieves this minimax lower bound, provided a suitable restricted eigenvalue condition is satisfied. Finally we prove that provided n ≥ K p where K = O(Rq ( pT)- q/2) represents an approximate sparsity level, our restricted eigenvalue condition and physical constraints are satisfied for random bounded ensembles. We also provide numerical experiments that validate our mean-squared error bounds. Our results address a number of open issues from prior work on Poisson inverse problems that focuses on strictly sparse models and does not provide guarantees for convex implementable algorithms.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…