Brevet FR3077908A1 — publié le 20190816

PROCEDE D'ESTIMATION PARAMETRIQUE D'UN OBJET A ESTIMER DANS UNE IMAGE NUMERIQUE, ET PROCEDE DE SUPPRESSION DE L'OBJET DANS L'IMAGE NUMERIQUE L'invention porte sur un procédé d'estimation paramétrique d'un objet à estimer dans une image numérique, et sur un procédé de suppression de l'objet dans l'image numérique. Elle s'applique notamment aux images numériques obtenues par imagerie à rayons X, mais peut être étendue à tout type d'imageur numérique, par exemple, à rayonnement infrarouge, ou encore dans le spectre visible. Dans la présente demande, trois composantes peuvent être présentes dans l'image numérique : 1. Le contenu ou fond clinique, correspondant à l'objet réel 15 imagé, que l'on cherche à visualiser. Il peut s'agir par exemple d'un patient, ou d'une zone d'un patient, en imagerie à rayons X. 2. Le bruit d'acquisition, correspondant aux perturbations aléatoires de la formation des images. Il peut s'agir du bruit électronique, auquel il faut ajouter, en imagerie à rayons X, le bruit photonique, dû au faible 20 nombre de photons présents dans l'image et qui suit une statistique de Poisson. 3. Les objets que l'on cherche à estimer, et qui peuvent s'avérer gênants pour la visualisation du contenu. 25 L'image numérique peut alors s'écrire de deux manières différentes selon que le bruit est traité indépendamment des deux autres composantes : Image = Objet + Contenu + bruit Ou considéré comme appartenant au contenu Image = Objet + Contenu 2 La figure 1 représente un exemple d'image comprenant un objet à estimer. En l'occurrence, cette figure représente une image radiologique dont le contenu est une zone d'un patient. Afin d'améliorer le contraste de l'image, et plus généralement la qualité de l'image radiologique, il 5 est courant pour le praticien d'utiliser une grille antidiffusante. Située entre le patient et le détecteur, elle permet de diminuer la contribution du rayonnement diffusé et causé par le patient. Sur la figure 1, les variations horizontales sont ainsi dues à la présence de la grille antidiffusante ; ces variations peuvent masquer une information importante, et conduire éventuellement à un 10 mauvais diagnostic. Un traitement d'image doit donc être appliqué, afin de retirer ces variations horizontales de l'image numérique, sans toutefois modifier le contenu. Un tel traitement comporte toutefois deux principales difficultés. Il faut tout d'abord pouvoir enlever uniquement l'objet sachant que l'image contient également du bruit, et que le contenu peut être un signal très fort et imprévisible. En effet, le contenu peut fortement varier selon la zone du patient à imager, et pour une même zone, le contenu peut varier d'un patient à l'autre. L'autre difficulté provient du fait que l'objet à estimer n'est pas parfaitement défini, à savoir que plusieurs paramètres de l'objet peuvent ne pas être connus. La figure 2 représente un autre exemple d'objet à estimer. Il s'agit d'un carré CR superposé à l'image IM, qui gêne la visualisation correcte du contenu. Ni l'amplitude de l'objet (la valeur du pixel, par exemple le niveau de gris pour une image de radiologie), ni sa taille ne sont connues. Actuellement, deux types de solutions permettent de traiter un signal afin d'en retirer des composantes non désirées. Le premier regroupe les filtres classiques, notamment le filtre adapté, le filtre de Wiener ou encore le filtre de Kalman. Ces solutions se fondent sur une estimation et sur une modélisation statistique du contenu et du signal gênant l'estimation du contenu (du bruit généralement aléatoire, décorrélé du contenu). Dans notre cas, le contenu (le patient) n'est pas modélisable et il représente pourtant la majorité du signal présent dans 3 l'image. De telles solutions peuvent donc être implémentées dans des zones dites « plates » de l'image où le patient n'est pas présent, mais demeurent sous-optimales lorsque le patient est présent. Ces solutions ont de plus l'inconvénient de modifier plus ou moins fortement le contenu de l'image lors 5 de la suppression de l'objet, par exemple si l'image comprend une composante qui ressemble à l'objet à supprimer. En référence à la figure 1, s'il y a, dans le contenu acquis par le détecteur, un élément ressemblant à la grille antidiffusante, cet élément sera potentiellement supprimé en cas de filtrage adapté. 10 Le deuxième type de solutions regroupe les méthodes d'estimation spectrale, notamment la méthode IAA (Iterative Adaptive Approach), décrite notamment dans l'article « Source localization and sensing: A nonparametric Iterative Adaptive Approach based on weighted least squares » (Yardibi T. et al., Aerospace and Electronic Systems, IEEE 15 Transactions on 46.1 2010), et la méthode SLIM (Sparse Learning via Iterative Minimization), décrite notamment dans l'article « Sparse Leaming via Iterative Minimization with application to MIMO Radar Imaging » (Tan X. et al., IEEE Transactions on Signal Processing, vol. 59, no. 3, March 2011). Ces méthodes sont adaptées au domaine des radars et des télécommunications, 20 pour lesquels le signal est unidimensionnel, mais ne peuvent en l'état pas être transposées au cas d'un signal bidimensionnel, et être ainsi utilisées pour du traitement d'image. Par ailleurs, dans la méthode SLIM, le bruit est supposé être blanc. Ce postulat ne peut être utilisé pour une image acquise par un détecteur, où le bruit est propre à chaque technique d'imagerie. En particulier, 25 en imagerie à rayons X, le bruit photonique n'est pas blanc. Enfin, les méthodes précitées ne considèrent que deux signaux, à savoir le contenu et le bruit. Elles ne permettent donc pas d'effectuer une estimation paramétrique d'un objet gênant la visualisation du contenu, et distinct du bruit d'acquisition. 30 L'invention vise donc à obtenir un procédé d'estimation paramétrique d'un objet à estimer dans une image numérique, distinct du bruit 4 d'acquisition, ne nécessitant de connaissance exacte préalable ni du contenu, ni de l'objet à estimer, ni du bruit d'acquisition. L'invention vise également un procédé de suppression de l'objet dans l'image numérique, qui ne modifie pas, dans l'image, le contenu correspondant à l'objet réel imagé. 5 Un objet de l'invention est donc un procédé d'estimation paramétrique d'un objet à estimer dans une image numérique représentant un contenu réel imagé, comprenant au moins : a) une étape initiale comprenant la création d'un dictionnaire de 10 composantes de contenu et la création d'un dictionnaire de composantes d'objet, les composantes de contenu et les composantes d'objet ayant les mêmes dimensions que l'image numérique ; b) une étape de détermination conjointe de l'amplitude de chacune des composantes de contenu du dictionnaire de composantes de contenu 15 et des composantes d'objet du dictionnaire de composantes d'objet présentes dans l'image numérique ; c) une étape de détermination de la valeur d'au moins un paramètre caractérisant l'objet à estimer à partir de l'amplitude de chacune des composantes d'objet. 20 Avantageusement, l'étape c) de détermination de la valeur d'au moins un paramètre caractérisant l'objet à estimer comprend les sous-étapes de : - détermination d'un sous-ensemble de composantes d'objet de 25 plus forte amplitude ; - détermination de la valeur dudit paramètre en fonction d'une combinaison linéaire des composantes d'objet dudit sous-ensemble. 30 Avantageusement, la création du dictionnaire de contenu comprend la génération d'un ensemble d'images de signaux sinusoïdaux bidimensionnels. 5 En variante, la création du dictionnaire de contenu comprend la génération d'un ensemble d'images représentatives de la fonction de transfert de modulation du capteur formant l'image numérique. 5 En variante, la création du dictionnaire de contenu est effectuée à partir d'un algorithme d'apprentissage automatique. Avantageusement, l'algorithme d'apprentissage automatique 10 est un algorithme K-SVD. Avantageusement, la création d'un dictionnaire de composantes d'objet comprend la génération d'un ensemble d'images ayant chacune une valeur différente du paramètre. 15 Avantageusement, le paramètre est sélectionné dans un ensemble comprenant la fréquence, la forme, la taille et l'emplacement de l'objet. 20 Avantageusement, l'étape b) de détermination conjointe est réalisée par détermination itérative d'un vecteur d'estimation (13) comprenant les différentes valeurs d'amplitude de chacune des composantes du dictionnaire de contenu et du dictionnaire d'objet présentes dans l'image, jusqu'à convergence du vecteur d'estimation (13). 25 Avantageusement, la convergence du vecteur d'estimation (13) est considérée atteinte lorsque la norme du vecteur d'estimation (13) évolue en-dessous d'un seuil prédéterminé entre deux itérations successives. 30 En variante, la convergence du vecteur d'estimation (13) est considérée atteinte au bout d'un nombre prédéterminé d'itérations. Avantageusement, le procédé comprend, à la fin de chaque itération, une étape de modification des composantes du dictionnaire de 6 composantes de contenu et/ou du dictionnaire de composantes d'objet pour minimiser une erreur de projection (e) définie par : E= iy - D1312, où y est le vecteur représentant l'image numérique, et D la matrice 5 représentant les composantes du dictionnaire à modifier. Avantageusement, l'étape b) de détermination conjointe comprend la détermination itérative d'une estimation de la variance de bruit (02), déterminée à partir du vecteur d'estimation (13) et d'une matrice N de 10 corrélation du bruit présent dans l'image numérique, jusqu'à convergence de l'estimation de ta variance de bruit (02). Avantageusement, la matrice N de corrélation du bruit est déterminée par une connaissance a priori du détecteur d'acquisition de 15 l'image numérique. Avantageusement, le procédé comprend une étape d) de reconstitution de l'image de contenu par combinaison linéaire de toutes les composantes de contenu. 20 Avantageusement, l'image numérique est acquise par imagerie à rayons X. L'invention se rapporte également à un procédé de 25 suppression d'un objet présent dans une image numérique représentant un contenu réel imagé, comprenant les étapes de : - estimation de la valeur d'au moins un paramètre de l'objet et de son amplitude conformément au procédé d'estimation paramétrique précité ; 30 - estimation de l'objet à partir de la valeur du paramètre caractérisant l'objet et de son amplitude ; soustraction pixel par pixel de l'objet estimé dans l'image numérique. 7 D'autres caractéristiques, détails et avantages de l'invention ressortiront à la lecture de la description faite en référence aux dessins annexés donnés à titre d'exemple et qui représentent, respectivement : 5 - la figure 1, déjà décrite, un exemple d'oscillation présente dans une image à rayons X, due à la présence d'une grille antidiffusante ; - la figure 2, un exemple d'objet à estimer dans une image numérique ; - la figure 3, un schéma fonctionnel du procédé selon 10 l'invention ; - la figure 4A, un exemple de dictionnaire de composantes de contenu, obtenu à partir de composantes de Fourier ; - la figure 4B, un autre exemple de dictionnaire de composantes de contenu, obtenu à partir de la réponse impulsionnelle du 15 détecteur ; - la figure 5, un exemple de dictionnaire de composantes d'objet, pour lequel le paramètre à estimer est la taille d'un carré ; - la figure 6, une illustration schématique des différentes étapes de l'algorithme de détermination conjointe des amplitudes de 20 composantes d'objet et de contenu, selon un premier mode de réalisation ; - la figure 7, une illustration schématique des différentes étapes de l'algorithme de détermination conjointe des amplitudes de composantes d'objet et de contenu, selon un deuxième mode de réalisation ; - la figure 8, un graphe illustrant l'amplitude de chaque 25 composante d'objet ; - la figure 9, l'image de la figure 2, dans laquelle l'objet a été retiré. La figure 2 illustre une image numérique IM, comprenant un 30 objet (un carré CR) qui blanchit une partie correspondante de l'image numérique. Sur la figure 2, les pointillés entourant le carré CR ne sont représentés que pour des raisons de clarté de la figure. On cherche donc à 8 estimer la valeur d'au moins un paramètre caractérisant l'objet, et l'amplitude de l'objet, pour estimer de combien le carré CR éclaircit l'image. La figure 3 illustre un schéma fonctionnel du procédé selon 5 l'invention. La première étape initiale a) consiste à créer un dictionnaire de composantes de contenu et à créer un dictionnaire de composantes d'objet. Une première sous-étape consiste donc à créer un 10 dictionnaire de composantes de contenu. Selon un premier mode de réalisation, le dictionnaire de composantes de contenu peut être créé sans a priori, en générant un ensemble d'images de signaux sinusoïdaux bidimensionnels basiques. La figure 4A illustre un tel dictionnaire. Chacune des composantes de contenu 15 (COMP_CONT_1, COMP_CONT_2, ..., COMP_CONT_K) (K supérieur ou égal à 1) représente un signal sinusoïdal bidimensionnel, à savoir à une fréquence particulière. Les composantes de contenu représentées sont appelées composantes de Fourier. Le nombre de composantes de contenu représenté sur la figure 4 n'est pas limitatif. 20 Selon un autre mode de réalisation, le dictionnaire de composantes de contenu peut être créé en générant un ensemble d'images représentatives de la fonction de transfert de modulation du capteur formant l'image numérique. La figure 4B illustre un tel dictionnaire. Ce dictionnaire est généré à partir de la réponse impulsionnelle du détecteur. Chacune des 25 composantes (COMP_CONT_1, COMP_CONT_2, ..., COMP_CONT_K) (K supérieur ou égal à 1) représente la réponse impulsionnelle du détecteur sur un pixel. Le dictionnaire de contenu est donc constitué de l'ensemble de ces images, pour chacun des pixels du détecteur. Ce dictionnaire est avantageusement proche des caractéristiques physiques du détecteur, 30 notamment en imagerie à rayons X. Selon un autre mode de réalisation, le dictionnaire de composantes de contenu peut être créé à partir d'un algorithme 9 d'apprentissage automatique (« machine learning » en terminologie anglo-saxonne), notamment à partir de l'algorithme K-SVD. L'algorithme K-SVD est particulièrement adapté à la création et à l'apprentissage de composantes de dictionnaire, ce qui le rend facilement implémentable. Le dictionnaire de 5 composantes de contenu utilisant un tel algorithme peut par exemple créer les composantes à partir d'une banque d'images de radiographies à rayons X. D'autres algorithmes d'apprentissage automatique, par exemple les réseaux de neurones, ou encore l'algorithme SVM (pour « Support Vector Machine », ou « machine à vecteurs de support »). 10 Les dimensions de chaque composante de contenu peuvent être les mêmes que les dimensions de l'image numérique dans laquelle se trouve l'objet. L'estimation paramétrique peut également être réalisée sur une partie de l'image numérique, par exemple si l'objet est très localisé. Dans ce cas, les dimensions de chaque composante de contenu sont les mêmes que 15 les dimensions de la partie de l'image numérique. Une deuxième sous-étape de la première étape initiale a) consiste à créer un dictionnaire de composantes d'objet. Le dictionnaire de composantes d'objet peut être créé en générant un ensemble d'images ayant chacune une valeur différente du paramètre. Le paramètre peut être, de façon 20 non exhaustive, la fréquence, la forme, la taille et l'emplacement de l'objet. La figure 5 un illustre un exemple de dictionnaire de composantes d'objet, pour lequel le paramètre à estimer est la taille du carré à estimer dans la figure 2. Le nombre de composantes d'objets (COMP_OBJ_1, COMP_OBJ_2, COMP_OBJ_3, COMP_OBJ_4) présentes 25 dans la figure 5, à savoir le nombre de valeurs possibles que peut prendre le paramètre à estimer, n'est pas limitatif. Il peut notamment comprendre autant de composantes qu'il y a de tailles différentes de carrés sur l'image numérique (par exemple 1x1 pixel, 2x2 pixels, 3x3 pixels, etc.). Le cas le plus fréquent est celui où il y a un paramètre à 30 estimer (par exemple la taille du carré), et l'amplitude de l'objet correspondant. Ainsi, le dictionnaire de composantes d'objet comprend autant de composantes d'objet qu'il y a de valeurs possibles du paramètre. La 10 composante d'objet ayant la plus forte amplitude correspond à la bonne valeur de paramètre. Il est également possible d'envisager le cas où il y a plusieurs paramètres à estimer (par exemple l'emplacement, la taille du carré), en plus de l'amplitude de l'objet correspondant. Si l'objet est une somme de 5 composantes connues mais d'amplitudes inconnues, il convient alors d'estimer les amplitudes de chacune des composantes. On dispose ainsi d'un dictionnaire de composantes de contenu, et d'un dictionnaire de composantes d'objet. Chaque composante est une matrice de I lignes et J colonnes que l'on écrira sous la forme d'un 10 vecteur de taille M, M étant la taille de l'image, ou de la partie de l'image dans laquelle se trouve l'objet (i.e. l'ensemble des pixels correspondants). La détermination conjointe de l'amplitude de chacune des composantes de contenu du dictionnaire de composantes de contenu et des composantes d'objet du dictionnaire de composantes d'objet présentes dans 15 l'image numérique (étape b) dans la figure 3) peut être effectuée conformément aux deux modes de réalisation décrits ci-après. 1" mode de réalisation 20 Selon ce premier mode de réalisation, le bruit présent dans l'image numérique est traité séparément du contenu et de l'objet. Ce traitement séparé peut être justifié par la nature différente des données traitées : les données du contenu et de l'objet ne sont pas aléatoires, alors que celles du bruit le sont. 25 Soit D une matrice de taille MxL contenant les différentes composantes à estimer (objet et contenu). La matrice D est donc une concaténation du dictionnaire de composantes de contenu et du dictionnaire de composantes d'objet. L est le nombre de composantes du dictionnaire de composantes de contenu et du dictionnaire de composantes d'objet. Di, qui 30 est un vecteur de taille M, est la ième composante de la matrice D. Les informations de corrélation du bruit présent dans l'image numérique se retrouvent dans une matrice N, distincte de la matrice D. Le 11 bruit présent dans l'image numérique est donc traité séparément du contenu et de l'objet. La matrice N de corrélation du bruit peut être déterminée par une connaissance a priori du détecteur d'acquisition de l'image numérique. Le bruit peut en effet être propre à chaque technique d'imagerie (rayons X, 5 infrarouges, visible). En particulier, en imagerie à rayons X, le bruit photonique n'est généralement pas un bruit blanc. La matrice N de corrélation du bruit peut être estimée d'après le théorème de Wiener-Khintchine, en prenant la transformée de Fourier inverse de la densité spectrale de puissance d'une image « plate », à savoir sans contenu. 10 En variante, si on suppose qu'il n'y a aucune corrélation entre les pixels, N peut être une matrice identité de taille MxM. Soit y un vecteur représentant l'image numérique, à savoir la valeur (également appelée luminosité) de chaque pixel. 15 Soit 13 un vecteur d'estimation de taille L contenant les amplitudes estimées de chaque composante Di. Ainsi, 13 = [131, P2, PL] Soient a, b et q des paramètres de l'algorithme, a et b dépendent de la nature et de la connaissance du bruit. S'il n'y a pas de connaissance particulière sur le bruit, alors a=0 et b=0. 20 L'algorithme de détermination conjointe de l'amplitude de chacune des composantes de contenu du dictionnaire de composantes de contenu et des composantes d'objet du dictionnaire de composantes d'objet présentes dans l'image numérique comprend une étape d'initialisation, et six étapes d'itération, comme l'illustre la figure 6. 25 Etape d'initialisation La valeur du vecteur d'estimation 13 peut être initialisée d'une manière quelconque (par exemple avec des valeurs aléatoire), ou bien en effectuant le produit scalaire entre chaque composante Di et le vecteur y 30 représentant l'image numérique : pi= Di. y, pour 1 i L 12 Par ailleurs, une estimation de la variance de bruit (a2) est initialisée, soit d'une manière quelconque (par exemple avec des valeurs aléatoire), soit avec la relation suivante : 2 = [Y - DeN-1[37 - DP] + a o- M + b 13 étant à sa valeur initiale. Etapes d'itération 1. On calcule P = 113IqZ , où q est un paramètre fixé par l'utilisateur ; il s'agit d'un réel compris entre 1 et 2. L'utilisateur peut changer la valeur de ce paramètre, en fonction du résultat final, à savoir en fonction de la qualité de la suppression de l'objet dans l'image. Z est une matrice diagonale permettant de privilégier l'estimation de certaines composantes, dans le cas général Z sera la matrice identité. Z est paramétrable par l'utilisateur. 2. On calcule une matrice R d'estimation de la matrice de covariance du signal y de l'image numérique : R = Ddiag(P)DH + o-2N 3. On calcule un terme de normalisation : Nor = 113irq 4. On calcule ensuite chaque amplitude estimée de composante pi par la formule : - D'H R-ly I3i Nor 5. L'estimation de la variance de bruit o-2 peut être calculée par exemple par la formule suivante : 2 [Y - DeN-1[37 - DP] + a o- = M + b 13 6. De façon facultative, les composantes de dictionnaire peuvent être modifiées pour minimiser une erreur de projection E définie par : £ = ly - D1312 Chaque composante Di de dictionnaire peut notamment être mise à jour avec un paramètre de mise à jour À compris entre 0 et 1 : y DI3 + D!,ew = (1- À)Dioid + À Ri Où Di," est la composante Di de dictionnaire à mettre à jour, et D!,ew la composante Di de dictionnaire mise à jour. Le paramètre de mise à jour À peut être fixé à 0, il n'y a dans ce cas pas de mise à jour des composantes du dictionnaire. La mise à jour des composantes Di de dictionnaire peut être réalisée par d'autres méthodes, notamment par dichotomie, ou par descente de gradient. Cette étape de mise à jour peut s'appliquer aux composantes de contenu du dictionnaire de composantes de contenu et/ou aux composantes d'objet du dictionnaire de composantes d'objet. 20 Ces six étapes sont répétées jusqu'à convergence de l'algorithme, à savoir lorsque le vecteur d'estimation 13 et l'estimation de la variance de bruit o-2 convergent. La variance de bruit o-2 étant directement dépendante de 13, la convergence de l'un implique la convergence de l'autre. La convergence du vecteur d'estimation 13 et la convergence 25 de l'estimation de la variance de bruit o-2 peuvent être considérées atteintes respectivement lorsque la norme du vecteur d'estimation 13 évolue en-dessous d'un seuil prédéterminé entre deux itérations successives, et lorsque la valeur 5 10 15 14 de l'estimation de la variance de bruit o-2 évolue en-dessous d'un seuil prédéterminé entre deux itérations successives. En variante, la convergence peut être considérée atteinte au bout d'un nombre prédéterminé d'itérations. Il est alors possible de déterminer 5 la valeur du ou des paramètres de l'objet (étape c) dans la figure 3). L'étape c) comprend les sous-étapes de : - détermination d'un sous-ensemble de composantes d'objet de plus forte amplitude ; - combinaison linéaire des composantes d'objet dudit sous- 10 ensemble. Le sous-ensemble de composantes d'objet de plus forte amplitude peut notamment comprendre une seule composante. A partir de l'algorithme décrit ci-dessus, il est également 15 possible de reconstituer le contenu, correspondant au contenu réel imagé (étape d) dans la figure 3), en effectuant une combinaison linéaire de toutes les composantes de contenu. Cette reconstitution permet d'obtenir une image sans bruit. 20 2ème mode de réalisation Selon ce deuxième mode de réalisation, le bruit présent dans l'image numérique est traité avec le contenu et l'objet, à savoir dans une seule et même matrice. 25 Par rapport au premier mode de réalisation, la matrice D est donc le résultat de la concaténation du dictionnaire de composantes d'objet, du dictionnaire de composantes de contenu, et des informations de corrélation du bruit présent dans l'image numérique (qui constituent la matrice N dans le premier mode de réalisation). 30 y est le vecteur représentant l'image numérique, à savoir la valeur de chaque pixel. 15 13 est le vecteur d'estimation de taille L contenant les amplitudes estimées de chaque composante Di. Ainsi, 13 = [131, P2, PL] q est un paramètre de l'algorithme. Par rapport au premier mode de réalisation, il n'y a pas de calcul de l'estimation de la variance de 5 bruit o-2 ; les paramètres a et b ne sont donc pas utilisés. L'algorithme de détermination conjointe comprend également une étape d'initialisation, et cinq étapes d'itération, comme l'illustre la figure 7. 10 Etape d'initialisation La valeur du vecteur d'estimation 13 peut être initialisée d'une manière quelconque (par exemple avec des valeurs aléatoire), ou bien en effectuant le produit scalaire entre chaque composante Di et le vecteur y représentant l'image numérique : 15 pi= Di. y, pour 1 i L Etapes d'itération 1. On calcule P = 113IqZ , où q est un paramètre fixé par l'utilisateur ; il s'agit d'un réel compris entre 1 et 2. 20 L'utilisateur peut changer la valeur de ce paramètre, en fonction du résultat final, à savoir en fonction de la qualité de la suppression de l'objet dans l'image. Z est une matrice diagonale permettant de privilégier l'estimation de certaines composantes, dans le cas général Z sera la 25 matrice identité. Z est paramétrable par l'utilisateur. 2. On calcule une matrice R d'estimation de la matrice de covariance du signal y de l'image numérique : R = Ddiag(P)DH 3. On calcule un terme de normalisation : Nor = 16 4. On calcule ensuite chaque amplitude estimée de composante pi par la formule : -H D' R-1-37 Ri = Nor 5. De façon facultative, les composantes du dictionnaire peuvent être modifiées pour minimiser une erreur de 5 projection E, de la même manière que pour le premier mode de réalisation. Ces cinq étapes sont répétées jusqu'à convergence de l'algorithme, à savoir lorsque le vecteur d'estimation 13 converge. 10 La convergence du vecteur d'estimation 13 peut être considérée atteinte lorsque la norme du vecteur d'estimation 13 évolue en-dessous d'un seuil prédéterminé entre deux itérations successives. En variante, la convergence peut être considérée atteinte au bout d'un nombre prédéterminé d'itérations. Les étapes c) et d) sont les 15 mêmes que pour le premier mode de réalisation. L'invention se rapporte aussi à un procédé de suppression d'un objet présent dans une image numérique représentant un contenu réel imagé, comprenant les étapes de : 20 - estimation de la valeur d'au moins un paramètre de l'objet et de son amplitude conformément au procédé d'estimation paramétrique précité ; - détermination de l'objet à partir du paramètre de l'objet et de son amplitude ; 25 - soustraction pixel par pixel de l'objet dans l'image numérique. En utilisant le premier mode de réalisation de l'algorithme, il est possible d'obtenir une estimation de la taille et de l'amplitude du carré. La figure 8 représente le résultat de l'exécution de l'algorithme. La valeur de 30 paramètre de chaque composante du dictionnaire d'objet se trouve en 17 abscisse, et l'amplitude estimée de chacune de ces composantes en ordonnée (en échelle logarithmique). Il ressort de cette figure que l'algorithme a détecté un carré de taille 7 avec une forte amplitude. Ainsi, le carré situé à la figure 2 est de largeur 7 et a pour amplitude 52. En utilisant ces 5 informations, il est alors possible de corriger l'image de sorte à obtenir une correction parfaite, comme l'illustre la figure 9. La méthode a été décrite pour la totalité d'une image formée par le détecteur. Elle peut également être appliquée à une partie de l'image 10 formée par le détecteur, en adaptant les dimensions des matrices et des vecteurs nécessaires à l'estimation paramétrique. La méthode a été décrite dans le cas d'un objet additif : Image = Objet + Contenu Elle peut être appliquée au cas d'un objet multiplicatif: 15 Image = Objet X Contenu Par exemple, la grille antidiffusante est un objet multiplicatif, la valeur de la fréquence et son amplitude dépendent alors du patient. Pour passer du cas additif au cas multiplicatif, il suffit de calculer le logarithme de l'image (donc du vecteur y) : 20 (Image)=log(Objet)+log(Contenu) En utilisant la méthode additive décrite précédemment pour obtenir le vecteur d'estimation 13, il est possible d'estimer log(Objet) et d'en déduire : Objet = 10 log (Objet) La présente invention peut être utilisée dans n'importe quel 25 procédé nécessitant une amélioration de la qualité de l'image. Elle s'applique entre autre au domaine de l'imagerie médicale (Radiographies à rayons X, scanner, IRM, etc...), ainsi que dans tous les autres types d'imageurs (rayons X, visible, infrarouge). 30 Les applications concrètes concernent par exemple : 18 - La suppression d'une perturbation, différente du bruit, qui dégrade la visibilité et qui peut entrainer un mauvais diagnostic dans le cas de l'imagerie médicale, comme l'illustre la figure 1. L'invention permet avantageusement de laisser le bruit dans l'image numérique, ce qui peut être 5 recherché par certains praticiens en imagerie médicale ; - La caractérisation d'objets à des fins techniques ou médicales (par exemple la taille d'une grosseur dans le corps du patient) ; - Le débruitage : si le premier mode de réalisation est utilisé, il est possible d'obtenir une estimation de l'objet, du contenu et du 10 bruit (en connaissant l'estimation de la variance de bruit a2). Il est alors possible de reconstruire le signal seulement avec les composantes de l'objet et du contenu. La même opération peut être effectuée sans objet, en estimant seulement le contenu et le bruit. 15 Le procédé selon l'invention est ainsi capable d'estimer correctement chacune des composantes des dictionnaires bien que ceux-ci comportent de nombreuses composantes. L'utilisation d'un procédé itératif permet d'estimer REVENDICATIONS1. Procédé d'estimation paramétrique d'un objet à estimer dans une image numérique représentant un contenu réel imagé, comprenant 5 au moins : a) une étape initiale comprenant la création d'un dictionnaire de composantes de contenu et la création d'un dictionnaire de composantes d'objet, les composantes de contenu et les composantes d'objet ayant les mêmes dimensions que l'image numérique ; 10 b) une étape de détermination conjointe de l'amplitude de chacune des composantes de contenu du dictionnaire de composantes de contenu et des composantes d'objet du dictionnaire de composantes d'objet présentes dans l'image numérique ; c) une étape de détermination de la valeur d'au moins un paramètre 15 caractérisant l'objet à estimer à partir de l'amplitude de chacune des composantes d'objet. 2. Procédé d'estimation paramétrique selon la revendication 1, l'étape c) de détermination de la valeur d'au moins un paramètre 20 caractérisant l'objet à estimer comprenant les sous-étapes de : détermination d'un sous-ensemble de composantes d'objet de plus forte amplitude ; détermination de la valeur dudit paramètre en fonction d'une combinaison linéaire des composantes d'objet dudit sous- 25 ensemble. 3. Procédé d'estimation paramétrique selon l'une des revendications précédentes, la création du dictionnaire de contenu comprenant la génération d'un ensemble d'images de signaux sinusoïdaux 30 bidimensionnels. 20 4. Procédé d'estimation paramétrique selon l'une des revendications 1 ou 2, la création du dictionnaire de contenu comprenant la génération d'un ensemble d'images représentatives de la fonction de transfert de modulation du capteur formant l'image numérique. 5. Procédé d'estimation paramétrique selon l'une des revendications 1 ou 2, la création du dictionnaire de contenu étant effectuée à partir d'un algorithme d'apprentissage automatique. 6. Procédé d'estimation paramétrique selon la revendication 5, l'algorithme d'apprentissage automatique étant un algorithme K-SVD. 7. Procédé d'estimation paramétrique selon l'une des revendications précédentes, la création d'un dictionnaire de composantes 15 d'objet comprenant la génération d'un ensemble d'images ayant chacune une valeur différente du paramètre. 8. Procédé d'estimation paramétrique selon l'une des revendications précédentes, le paramètre étant sélectionné dans un ensemble 20 comprenant la fréquence, la forme, la taille et l'emplacement de l'objet. 9. Procédé d'estimation paramétrique selon l'une des revendications précédentes, l'étape b) de détermination conjointe étant réalisée par détermination itérative d'un vecteur d'estimation (13) comprenant 25 les différentes valeurs d'amplitude de chacune des composantes du dictionnaire de contenu et du dictionnaire d'objet présentes dans l'image, jusqu'à convergence du vecteur d'estimation (13). 10. Procédé d'estimation paramétrique selon la revendication 30 précédente, la convergence du vecteur d'estimation (13) étant considérée atteinte lorsque la norme du vecteur d'estimation (13) évolue en-dessous d'un seuil prédéterminé entre deux itérations successives. 21 11. Procédé d'estimation paramétrique selon la revendication 9, la convergence du vecteur d'estimation (13) étant considérée atteinte au bout d'un nombre prédéterminé d'itérations. 5 12. Procédé d'estimation paramétrique selon l'une des revendications 9 à 11, comprenant, à la fin de chaque itération, une étape de modification des composantes du dictionnaire de composantes de contenu et/ou du dictionnaire de composantes d'objet pour minimiser une erreur de projection (c) définie par : 10 £ = ly - D1312, où y est le vecteur représentant l'image numérique, et D la matrice représentant les composantes du dictionnaire à modifier. 13. Procédé d'estimation paramétrique selon l'une des revendications 9 à 12, l'étape b) de détermination conjointe comprenant la détermination itérative d'une estimation de la variance de bruit (a2), déterminée à partir du vecteur d'estimation (13) et d'une matrice N de corrélation du bruit présent dans l'image numérique, jusqu'à convergence de l'estimation de la variance de bruit (a2). 14. Procédé d'estimation paramétrique selon la revendication précédente, la matrice N de corrélation du bruit étant déterminée par une connaissance a priori du détecteur d'acquisition de l'image numérique. 15. Procédé d'estimation paramétrique selon l'une des revendications précédentes, comprenant une étape d) de reconstitution de l'image de contenu par combinaison linéaire de toutes les composantes de contenu. 16. Procédé d'estimation paramétrique selon l'une des revendications précédentes, l'image numérique étant acquise par imagerie à rayons X. 22 17. Procédé de suppression d'un objet présent dans une image numérique représentant un contenu réel imagé, comprenant les étapes de : estimation de la valeur d'au moins un paramètre de l'objet et de 5 son amplitude conformément au procédé d'estimation paramétrique selon l'une des revendications précédentes ; estimation de l'objet à partir de la valeur du paramètre caractérisant l'objet et de son amplitude ; soustraction pixel par pixel de l'objet estimé dans l'image 10 numérique.