A half-quadratic block-coordinate descent method for spectral estimation
Signal Processing
| Hi-index | 35.68 |
Formulated as a linear inverse problem, spectral estimation is particularly underdetermined when only short data sets are available. Regularization by penalization is an appealing nonparametric approach to solve such ill-posed problems. Following Sacchi et al. (see ibid., vol.46, no.1, p.32-38, 1998), we first address line spectra recovering in this framework. Then, we extend the methodology to situations of increasing difficulty: the case of smooth spectra and the case of mixed spectra, i.e., peaks embedded in smooth spectral contributions. The practical stake of the latter case is very high since it encompasses many problems of target detection and localization from remote sensing. The stress is put on adequate choices of penalty functions: following Sacchi et al., separable functions are retained to retrieve peaks, whereas Gibbs-Markov potential functions are introduced to encode spectral smoothness. Finally, mixed spectra are obtained from the conjunction of contributions, each one bringing its own penalty function. Spectral estimates are defined as minimizers of strictly convex criteria. In the cases of smooth and mixed spectra, we obtain nondifferentable criteria. We adopt a graduated nondifferentiability approach to compute an estimate. The performance of the proposed techniques is tested on the well-known Kay and Marple (1982) example