Acquisition comprimée

L'acquisition comprimée (en anglais compressed sensing) est une technique permettant de trouver la solution la plus parcimonieuse d'un système linéaire sous-déterminé. Elle englobe non seulement les moyens pour trouver cette solution mais aussi les systèmes linéaires qui sont admissibles. En anglais, elle porte le nom de Compressive sensing, Compressed Sampling ou Sparse Sampling.

Histoire

modifier

Plusieurs disciplines scientifiques utilisent les techniques d'acquisition comprimée depuis très longtemps[réf. souhaitée], mais ce n'est que grâce aux articles d'Emmanuel Candès, Justin Romberg, Terence Tao et David Donoho que l'on a commencé à comprendre quels étaient les ensembles de mesures qui étaient admissibles. En particulier, on s'est aperçu que l'un des moyens de retrouver des signaux parcimonieux était la résolution de problèmes d'optimisation linéaires faisant intervenir la norme  [1]. En statistique, la méthode des moindres carrés a été complétée par la norme   introduite par Laplace. Après l'introduction de l'optimisation linéaire et de l'algorithme du simplexe de George Dantzig, la norme   a été utilisée en statistique numérique. En théorie statistique la norme   a été utilisée par George Brown et plus tard elle servit à définir des estimateurs non-biaisés de la médiane. Elle a été utilisée par Peter Huber et d'autres dans leurs travaux sur les statistiques robustes. Dans les années 1970, la norme   a été utilisée par les sismologues pour analyser les données de réflexions d'ondes entre les couches de l'écorce terrestre lorsqu'elles ne semblaient pas respecter les conditions du théorème d'échantillonnage de Nyquist-Shannon[2]. Toujours en traitement du signal, elle a été introduite en 1993 dans l'algorithme « Matching Pursuit » et dans l'estimateur Lasso par Robert Tibshirani en 1996[3] et la poursuite de base (basis pursuit) en 1998[4]. Des résultats théoriques permettaient de savoir quand ces algorithmes permettaient de retrouver des solutions éparses, mais le type et le nombre de mesures requis pour cela étaient sous-optimaux et l'algorithme d'acquisition comprimée a permis une amélioration sensible de ce domaine du traitement du signal. Vers 2004, Emmanuel Candès, Terence Tao et David Donoho découvrirent des résultats importants sur le nombre minimum de données nécessaires à la reconstruction d'une image qui apparut plus petit que le nombre minimum déterminé par le critère de Nyquist-Shannon[5],[6].

Principes

modifier

On considère un signal   à valeur réelle, à une dimension et de longueur finie, que l’on notera sous la forme d’un vecteur colonne de longueur  . Dans le cas d’images, ou de signaux de plus grande dimension, on peut agencer les données pour former un long vecteur.

On exprime ce signal dans une «base  », ce que l'on exprime par la relation linéaire  . On souhaite choisir la base telle que la majorité des coefficients de   soient nuls, c'est-à-dire telle que   soit parcimonieux dans la base  . Ces bases sont connues et habituellement utilisées pour comprimer le signal de départ. L’acquisition comprimée propose d'acquérir directement la version comprimée du signal ce qui évite de traiter les échantillons inutiles[7].

Le processus de mesure linéaire utilisé consiste à faire   produits scalaires, avec  , entre   et une collection de vecteurs  ,   allant de 1 à  . On obtient ainsi les échantillons de mesures  . Considérant la matrice de mesure   ayant pour colonnes les  , on peut alors écrire   comme :   (avec  ).

La conception de la matrice de mesure doit permettre de retrouver les signaux les plus parcimonieux. Pour cela, ces matrices de mesure doivent suivre certaines propriétés dont l'une est appelée RIP (Restricted Isometry Property). De plus   doit être incohérente par rapport à  . Il s’avère que ces propriétés sont satisfaites avec une grande probabilité simplement en choisissant aléatoirement  , par exemple en utilisant une distribution gaussienne. Bien que ces conditions RIP ou d’incohérence soient satisfaites pour certains ensembles de mesures, ce ne sont que des conditions nécessaires et non suffisantes. Il existe d'autres propriétés, aussi non suffisantes, qui permettent d'avoir besoin d'encore moins échantillons pour garantir une reconstruction de signal parcimonieux[8].

La dernière étape du processus d’acquisition comprimée est la reconstruction du signal de départ. On connaît pour cela les   valeurs de  , la matrice de mesure utilisée ainsi que la base  . L’algorithme de reconstruction cherche à retrouver les coefficients de  , il est ensuite simple grâce à la connaissance de   de retrouver  , le signal de départ.

Il existe une infinité de solutions   à l’équation  , cependant on cherche la solution minimisant une certaine norme. La norme   mesure l’énergie du signal, c’est pourquoi en faisant une  -minimisation, on ne retrouvera quasiment jamais un résultat qui est K-parcimonieux. La norme   mesure la parcimonie du signal (on compte le nombre d’éléments non nuls), l’optimisation   tel que   donne un bon résultat. L’inconvénient est que ce problème n’est pas calculable numériquement. Il a été démontré que l’optimisation basée sur la norme   permet de retrouver exactement un signal K-parcimonieux, de plus le problème de la  -minimisation est convexe[9].

Sa résolution peut se réduire en un problème d'optimisation linéaire, plus connu sous le nom de poursuite de base (basis pursuit en anglais).

Applications

modifier

Une première application importante est en tomographie, en particulier pour l'IRM et la tomodensitométrie[10],[11],[12].

Il est également envisagé d'appliquer l'acquisition comprimée pour le radar[13] et l'imagerie radar[14]

En 2013, les laboratoires Bell annoncent avoir conçu une caméra sans lentilles, fonctionnant selon le principe de l'acquisition comprimée, avec un capteur ponctuel unique placé derrière un réseau d'orifices programmables[15].

Notes et références

modifier
(en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Compressed sensing » (voir la liste des auteurs).

Références

modifier
  1. List of L1 regularization ideas from Vivek Goyal, Alyson Fletcher, Sundeep Rangan, The Optimistic Bayesian: Replica Method Analysis of Compressed Sensing
  2. (en) Brian Hayes, « The Best Bits », American Scientist, (version du sur Internet Archive)
  3. The Lasso page, at Robert Tibshirani's homepage. "Regression shrinkage and selection via the lasso". J. Royal. Statist. Soc B., Vol. 58, No. 1, pages 267-288
  4. "Atomic decomposition by basis pursuit", by Scott Shaobing Chen, David L. Donoho, Michael A. Saunders. SIAM Journal on Scientific Computing
  5. E. J. Candès, J. Romberg and T. Tao. Stable signal recovery from incomplete and inaccurate measurements. Comm. Pure Appl. Math., 59 1207-1223 (www-stat.stanford.edu) [PDF]
  6. Donoho, D. L., Compressed Sensing, IEEE Transactions on Information Theory, V. 52(4), 1289–1306, 2006 (ieeexplore.ieee.org)
  7. Holtz, Olga V. An Introduction to Compressive Sensing. Stanford : s.n., 2009
  8. Measurement matrices and admissibility conjditioons in compressive sensing - a comprehensive list - - en anglais-
  9. Baraniuk, R.G., "Compressive Sensing [Lecture Notes]", Signal Processing Magazine, IEEE, vol.24, no.4, p. 118,121, July 2007
  10. Compressive sampling makes medical imaging safer
  11. (en) « Radiation dose reduction in computed tomography: techniques and future perspective », Imaging Med,‎ (lire en ligne)
  12. Pouryazdanpanah Kermani, Ali. "Sparse MRI and CT Reconstruction." (2017).
  13. (en) M.A. Herman et T. Strohmer, « High-Resolution Radar via Compressed Sensing », IEEE Transactions on Signal Processing, vol. 57,‎ , pp. 2275-2284
  14. (en) Potter, L.C., et al., et al., « Sparsity and Compressed Sensing in Radar Imaging. », Proceedings of the IEEE, vol. 98,‎ , pp. 1006-1020
  15. (en) Bell labs invents lensless camera

Voir aussi

modifier

Articles connexes

modifier

Lien externe

modifier