Brevet FR2863749A1 — publié le 20050617

L'invention concerne le domaine de l'imagerie radiologique et plus particulièrement celui de la mammographie. La mammographie est aujourd'hui largement utilisée pour la détection de lésions et la prévention de cancers du sein. Les premiers signes que les radiologues cherchent à détecter sur les images mammographiques sont des dépôts de calcium appelés microcalcifications qui constituent généralement des éléments plus opaques aux rayons X que les tissus environnants. Les radiologues scrutent les images à la recherche des signes de ces microcalcifications. Grâce à leur expérience, les radiologues sont capables dans une certaine mesure de faire la différence entre des microcalcifications malignes (souvent groupées en clusters) et des dépôts calciques bénins. A cet effet, ils interprètent la forme et la brillance de chaque microcalcification individuelle ainsi que la forme et l'étendue des clusters qu'ils observent sur les images radiographiques. Toutefois, il arrive que les radiologues ne repèrent pas certains dépôts. Ce phénomène a plusieurs causes. Notamment, les images mammographiques étant le résultat de projections, elles représentent des structures superposées qui perturbent la visibilité des microcalcifications. En outre, les radiologues doivent couramment analyser de nombreuses images mammographiques à la suite, il s'ensuit une fatigue visuelle qui peut affecter leur capacité à repérer les clusters. C'est pourquoi, il existe aujourd'hui des outils de Diagnostic Assisté par Ordinateur (DAO) qui permettent d'aider les radiologues dans la détection et la caractérisation des lésions. Un tel système est par exemple décrit dans le document US 6 137 898. Ces systèmes de DAO sont capables de traiter des images film préalablement numérisées par un scanner et de détecter les clusters de microcalcifications malignes. Ces systèmes indiquent les positions des clusters sur l'image par des marqueurs tels que des rectangles ou des cercles disposés autour des clusters ou encore des flèches. Le radiologue peut sélectionner la zone ainsi marquée pour la visualiser en gros plan. Ces systèmes comprennent des algorithmes basés sur des paramètres caractéristiques leur permettant de distinguer de manière automatique les éléments malins des éléments bénins. L'inconvénient de ces systèmes est qu'ils tendent à se substituer au radiologue dans l'interprétation de l'image, notamment en ce qui concerne la distinction entre éléments malins et éléments bénins. Or, ces systèmes ne sont pas suffisamment fiables pour se substituer totalement au radiologue. Le but de l'invention est donc de proposer un système de DAO permettant d'aider les radiologues dans leur travail de repérage des microcalcifications 'sans toutefois se substituer à eux dans l'étape d'interprétation de l'image. A cet effet, l'invention propose un procédé de traitement d'image 15 radiologique comprenant les étapes consistant à : - localiser sur une image radiologique fournie par un appareil d'imagerie des éléments susceptibles de constituer des signes de microcalcification, - afficher une image sur laquelle l'intensité des éléments ainsi 20 localisés a été rehaussée. Chaque élément localisé sur l'image radiologique est constitué par un amas de pixels (ou points) connexes. Le procédé de l'invention permet de mettre en évidence tous les signes de microcalcification qui apparaissent sur l'image sans distinction, que ces signes soient groupés ou isolés, bénins ou malins. Par conséquent, le radiologue obtient une vue d'ensemble immédiate de la cartographie de répartition des microcalcifications. Il peut directement apprécier la forme et l'étendue des clusters sans nécessairement commander un zoom sur les zones concernées. Dans tous les cas, l'appréciation de l'intérêt clinique des signes rehaussés est laissée au radiologue. Dans une mise en oeuvre de l'invention, le rehaussement de l'intensité des éléments susceptibles de constituer des microcalcifications est réalisé sur une image dite de présentation visualisable par le radiologue, distincte de l'image brute à partir de laquelle sont détectés les éléments. Cette image de présentation peut être obtenue par d'autres 5 procédés de traitements destinés à améliorer le confort de lecture du radiologue. Bien entendu, il existe plusieurs manières possibles de rehausser l'intensité des éléments. Dans une mise en oeuvre de l'invention, le rehaussement de l'intensité des éléments comprend une étape consistant à appliquer à chaque point de l'image radiologique appartenant à l'un de ces éléments un facteur d'amplification de la différence de niveau de gris entre ledit point et un ensemble de points situés dans un voisinage du point n'appartenant pas à l'élément considéré. Dans une mise en oeuvre de l'invention, le facteur d'amplification 15 appliqué aux points appartenant à un même élément est fonction d'un paramètre de similarité ii dudit élément avec un signal calcique. Dans une mise en oeuvre de l'invention, l'étape de localisation des éléments comprend: - une étape consistant à corréler, pour chacun des points d'un ensemble de points de l'image radiologique, l'intensité de l'image dans un voisinage du point avec chacune des fonctions d'une série de fonctions filtres, pour produire des images filtrées, - une étape consistant à sélectionner à partir des images filtrées des éléments susceptibles de constituer des signes de microcalcification. De manière avantageuse, les fonctions filtres sont des fonctions de type j3-spline. Ces fonctions présentent l'avantage de conduire à un filtrage rapide et autorisent une convolution dans le domaine spatial. Elles évitent donc des conversions notamment vers le domaine fréquentiel. L'invention concerne également un dispositif de traitement d'image radiologique caractérisé en ce qu'il comprend des moyens pour localiser sur une image radiologique fournie par un appareil d'imagerie des éléments susceptibles de constituer des signes de microcalcification et des moyens pour afficher une image sur laquelle l'intensité des éléments ainsi localisés a été rehaussée. D'autres caractéristiques et avantages ressortiront encore de la 5 description qui suit, laquelle est purement illustrative et non limitative et doit être lue en regard des figures annexées parmi lesquelles: - la figure 1 représente, sous forme d'un diagramme, les différentes phases et étapes d'une mise en oeuvre possible de l'invention - les figures 2A et 2B représentent des exemples de fonctions de 10 type (3-spline, - la figure 3 représente un exemple de fonction d'appartenance Fabre d'un pixel à la classe Fibre en fonction de l'élongation du segment auquel il appartient, - les figures 4A et 4B représentent respectivement une image de 15 présentation avant rehaussement et une image de présentation finale après rehaussement des signes de calcium. Dans la mise en oeuvre de l'invention représentée sur la figure 1, le procédé de traitement d'image est appliqué à une image brute R, c'est-àdire qui est fournie directement par le détecteur numérique d'un système de radiologie sans avoir été traitée au préalable. C'est à partir de cette image brute R que le procédé de traitement permet de localiser des éléments susceptibles de constituer des signes de microcalcification. Le procédé de traitement fait également intervenir une image dite de présentation sur laquelle ces éléments seront finalement rehaussés. Cette image de présentation S est obtenue par un autre procédé de traitement à partir de l'image brute R. Ce procédé permet de fournir une image directement visualisable par le radiologue. Le document EP 1 113 392 décrit par exemple un procédé de compensation en épaisseur permettant d'adapter les variations de niveaux de gris pour afficher une image dont l'apparence permet une lecture sur toute l'étendue du sein de la patiente. Sur la figure 1, le procédé de traitement d'image comprend cinq phases principales: 1) une phase 100 de prétraitement de l'image brute R, 2) une phase 200 d'analyse d'image, 3) une phase 300 de segmentation d'image, 4) une phase 400 de mesure de similarité, 5) une phase 500 de synthèse d'une image finale destinée à être visualisée par le radiologue. La première phase 100 de prétraitement consiste à convertir l'image 10 numérique brute R directement issue du dispositif de détection numérique en une image exploitable pour les traitements ultérieurs. Selon une première étape 110, l'image brute R est recadrée sur la région d'intérêt (ROI), c'est à dire la région de l'image contenant le sein de la patiente. Cette étape est destinée à appliquer la détection d'éléments susceptibles de constituer des signes de microcalcification uniquement à la région d'intérêt et ainsi à réduire les temps de traitement ultérieurs. Selon une deuxième étape 120, l'image brute recadrée obtenue est convertie en une image des contrastes C définie pour chaque pixel de l'image par: C(x,y) = K R(x,Y) M(x, Y) [1] M(x, y) où R(x,y) est l'intensité de niveau de gris de l'image brute R au pixel (x,y), M(x,y) est l'intensité moyenne de niveau de gris calculée au voisinage du pixel (x,y) sur l'image brute R, et K est une constante. Le voisinage du pixel (x,y) est par exemple délimité par une fenêtre de côté (2W+1)2 pixels, où W=4. La deuxième phase 200 d'analyse d'image consiste à filtrer l'image des contrastes C en utilisant un banc de filtres. Chaque filtre du banc est conçu pour réagir fortement à la présence de structures ayant une dimension et une orientation données. Selon cette phase 200, pour chaque pixel (xo, yo) de l'image, on définit une zone de voisinage f' de côté (2W+1)2 centré sur le pixel (xo,yo). On considère la fonction f(x,y) représentant l'intensité de niveau de gris de l'image des contrastes C dans la zone de voisinage du point (xo, yo) f'(x, y) = N(xo +x, yo + y) avec (x, y) e [-W,W[2 [2] On corrèle la fonction f' successivement avec chacun des filtres du banc de filtres. De manière avantageuse, ces filtres sont des fonctions de type ondelettes [ispline. Des exemples de fonctions [3-spline sont représentées aux figures 2A et 2B. Ces fonctions sont notamment décrites dans les publications suivantes: Michael UNSER, Akram ALDOURI & Murray EDEN, On the Asymptotic Convergence of [i-Spline Wavelets to Gabor Functions , IEEE Transactions on Information Theory, vol. 38, no. 2, pp. 864-872, Mars 1992, - Michael UNSER, Fast Gabor-Like Windowed Fourier and Continuous Wavelet Transforms , IEEE Signal Processing Letters. Vol. 1, no. 5, pp. 76-79, May 1994. Ce type de fonction permet un filtrage rapide et autorise une convolution dans le domaine spatial. Le banc comprend par exemple 16 ondelettes présentant 8 orientations différentes et 2 échelles différentes. Pour une ondelette gs,, d'échelle s et d'orientation a, on obtient le facteur de corrélation suivant: (l V;.f'(x,y).g(x,y).dx dy ll gs,a Il li +cc +GO 2 gs,a (t) dt, Çps,a e [0,2[ et W = 2s. _ce _.0 La troisième phase 300 de segmentation consiste à localiser à partir des images filtrées, des éléments P constitués par un amas de pixels connexes, susceptibles de constituer des signes de microcalcification. A cet effet, on génère à partir des images filtrées les deux images suivantes: - une image filtrée isotrope B, - une image filtrée directionnellement O. [3] avec llgs,a Selon une étape 310, on calcule l'image filtrée isotrope B comme la somme de toutes les images filtrées générées par la phase 200 d'analyse. On a donc: B(x, y) = N E Asa (x, y) É cos(w cos(a) x + co sin(a) É y ço a (x, y)) [4] s,a avec co = 7-c et N est le nombre de filtres composant le banc, N =16. 2s Selon une étape 311, on applique à l'image B un seuil TI de manière à sélectionner les pixels présentant une intensité supérieure à ce seuil. Ceci permet de ne conserver que les éléments d'intensité suffisante dans l'image filtrée isotrope. Selon une étape 312, on applique à l'image binaire ainsi obtenue un algorithme d'étiquetage des composantes connexes afin d'extraire des amas de pixels connectés entre eux (ces amas seront aussi appelés particules) dont l'intensité est supérieure au seuil TI dans l'image B. Les amas ainsi sélectionnés sont susceptibles de constituer des signes de microcalcification. Selon l'étape 320, on calcule l'image filtrée directionnellement O en appliquant à l'ensemble de l'image des contrastes C, le filtre qui présente la réponse la plus importante. On a donc: O(x, y) = A6 e (x, y) . cos(w É cos(0) x + w É sin(9) É y cp(x, y)) [5] avec (0%0) = Arg max As a (x, y) . s,a Selon une étape 321, on applique à l'image O un seuil T2 de manière à sélectionner les pixels présentant une intensité supérieure à ce seuil, ce qui permet de ne conserver que les éléments d'intensité suffisante. Selon une étape 322, on applique à l'image binaire ainsi obtenue un algorithme d'étiquetage des composantes connexes afin d'extraire des amas de pixels connectés entre eux (ces amas seront aussi appelés segments) dont l'intensité est supérieure à T2 dans l'image O. Selon une étape 323, on détermine pour chaque pixel (x, y) d'un segment son degré d'appartenance à la classe Fibre IlFibre(x, y) en fonction de l'élongation du segment auquel il appartient. La figure 3 représente un exemple de fonction d'appartenance ltFibre en fonction de l'élongation du segment auquel il appartient. Certaines des particules P détectées à l'étape 312 sont certainement des signes de microcalcification, d'autres n'en sont probablement pas. Pour la majorité de ces particules, il est difficile de décider si elles constituent ou non des dépôts calciques. La quatrième phase 400 consiste à déterminer pour chaque particule P un paramètre t ici m(P) qui caractérise la similarité entre la particule et un signal calcique. Ce paramètre icium(P) représente le degré d'appartenance de la particule P à la classe Signal Calcique . A cet effet, on génère à partir des images filtrées les deux images suivantes: - une image des amplitudes maximales Amax, - une image des amplitudes maximales en direction orthogonale Amax1 Selon une étape 410, on calcule l'image des amplitudes maximales Amax en appliquant à chaque pixel (x,y) de l'image des contrastes C, le filtre qui présente la réponse maximale. On a donc: Amax (x, y) = 4,0 (x, y) [6] avec (6, 0) = Arg max AS a (x, y) s,a Selon une étape 420, on calcule l'image des amplitudes maximales en direction orthogonale Amax' en appliquant à chaque pixel (x,y) de l'image des contrastes C, le filtre présentant une direction 01 orthogonale à la direction 0 du filtre présentant la réponse maximale. On a donc: Amax (x, y) = Acre,. (x, y) [7I Selon une étape 430, à partir de l'image des amplitudes maximales Amax obtenue à l'étape 410 et pour chaque particule P constituée par un nombre de pixels n, on calcule un paramètre Ap comme la moyenne des amplitudes maximales: A p = 1 A. (x, y) = 1 max As (x, y) [81 n (x,y)EP n (x,y)EP \ s'a Le paramètre Ap indique quelle est la brillance de la particule P par rapport à son environnement. Selon une étape 440, on réalise le même calcul à partir de l'image des amplitudes maximales dans la direction orthogonale Amaxl. On obtient 5 un paramètre Ap' : 1 1 1 AP 1 Amax (x, y) n (x,Y)EP Le paramètre Api caractérise à la fois la brillance et la forme de la particule P. Plus la particule P apparaît brillante et ronde, plus le paramètre Ap' sera élevé. Selon une étape 450, on calcule pour chaque particule P un paramètre 1. 1caicium(P) qui caractérise le degré d'appartenance de la particule P à la classe Signal Calcique . Ce paramètre caicium(P) est calculé en fonction des paramètres Ap et Ap' ainsi que du degré d'appartenance 1Fibre(xc, Yc) du pixel central ()cc, yc) de la particule P à la classe Fibre déterminée à l'étape 322. Ce paramètre caicium(P) tient ainsi compte à la fois de critères de taille, de forme et de brillance. Plus la particule P ressemble (respectivement ne ressemble pas) à un signal calcique, plus le paramètre caicium(P) obtenu sera grand et proche de 1 (respectivement de 0). Ce paramètre calcium(P) peut être obtenu par apprentissage sur une base d'images dans lesquelles les microcalcifications avérées ont été marquées (en utilisant par exemple un arbre de décision, un réseau de neurones ou bien les machines à vecteurs de support). La phase 500 de synthèse consiste à rehausser sur l'image de présentation S l'intensité des pixels appartenant aux particules P 25 sélectionnées à l'étape 312. Pour chaque pixel (x,y) de l'image de présentation S appartenant à une particule P, on augmente l'intensité du pixel par rapport à l'intensité de son environnement immédiat. De manière avantageuse, l'augmentation de l'intensité dépend de la similarité de la particule P avec un signal calcique. On calcule la nouvelle intensité FBF des pixels de la manière suivante: [9l 2863749 10 FBF(x, y) = Ms (x, y) + Pcaki, n (P)(SP (x, y) Ms (x, y)) [10] où Sp(x,y) est le niveau de gris d'un pixel (x,y) appartenant à la particule P, Ms(x,y) est le niveau de gris moyen des pixels situés dans un voisinage de 5 la particule P et n'appartenant pas à celle-ci. On a: Ms(x,y) # - S(x+k,y+l) [11] A où A est l'ensemble des points de l'image situés à une distance inférieure à N pixels (N typiquement égal à 5) et non nulle de la particule P et #A est le 10 cardinal de cet ensemble. Les particules P sont donc plus ou moins rehaussées en fonction de leur similitude avec un signal calcique. Ceci permet de mettre plus particulièrement en valeur les particules pour lesquelles la probabilité de constituer des microcalcifications est la plus forte. Les figures 4A et 4B représentent respectivement une image de présentation S avant rehaussement et une image de présentation finale FBF après rehaussement. L'image de présentation finale sur laquelle les particules P ont été rehaussées est affichée sur un dispositif d'affichage en vue d'être visualisée par un radiologue. Sur cette image finale, tous les signes de calcium ont été mis en évidence sans distinction quant à leur caractère suspect ou quant à leur caractère groupé ou isolé. Par conséquent, le radiologue obtient une vue d'ensemble immédiate de la cartographie de répartition des microcalcifications. II peut directement apprécier l'étendue de ces microcalcifications sans nécessairement zoomer sur les zones concernées. Le procédé de traitement comprend en outre une phase 600 optionnelle consistant à détecter les clusters de microcalcifications et à disposer des marqueurs à proximité ou autour de ces clusters. A cet effet, le procédé de traitement est apte à regrouper les particules présentant un paramètre I-icalcium(P) élevé caractérisant leur degré d'appartenance à la classe Signal Calcique , en clusters d'au moins 3 particules. Avantageusement, le radiologue peut choisir d'activer ou de désactiver cette option de marquage. Comme on l'aura compris, le rehaussement de l'intensité des éléments susceptibles de constituer des signes de microcalcification comprend la phase 500 consistant à appliquer à chaque point (x,y) de l'image radiologique appartenant à un élément P un facteur d'amplification de la différence de niveau de gris entre ledit point (x,y) et un ensemble de points (x+k,j+l) situés dans un voisinage du point (x,y) et n'appartenant pas à l'élément P. Comme on l'aura noté, l'étape 100 consiste à convertir l'image brute R en une image des contrastes C par une normalisation locale du contraste. L'invention se rapporte également à un procédé dans lequel, pour chaque particule P, le paramètre de similarité calcium(P) est déterminé selon les étapes suivantes: a) pour chacun des points (x,y) formant un élément P: - déterminer une amplitude maximale Amax en corrélant, l'intensité de l'image radiologique dans un voisinage W du point (x,y) avec une fonction filtre g6,0 conduisant à une réponse maximale, - déterminer une amplitude maximale en direction orthogonale Arr.-1 en corrélant, l'intensité de l'image radiologique dans un voisinage (W) du point (x,y) avec la fonction filtre g6,01 de direction orthogonale à la fonction filtre g6,o conduisant à une réponse maximale, b) pour chaque élément P: - déterminer un paramètre de brillance Ap comme une moyenne des amplitudes maximales Amax sur l'ensemble des points (x,y) de l'élément P, - déterminer un paramètre de brillance et de rondeur Ap' comme étant une moyenne des amplitudes maximales en direction orthogonale Amax1 sur l'ensemble des points (x,y) de l'élément P, - en déduire un paramètre de similarité calclum(P) de l'élément P à un signal calcique en fonction des paramètres de brillance et/ou de rondeur Ap, Api ainsi que du degré d'appartenance du point central ()cc, yc) de l'élément P à une classe Fibre . REVENDICATIONS 1. Procédé de traitement d'image radiologique comprenant les étapes consistant à : - localiser sur une image radiologique (R) fournie par un appareil d'imagerie des éléments (P) susceptibles de constituer des signes de microcalcification, - afficher une image (S) sur laquelle l'intensité des éléments (P) ainsi localisés a été rehaussée. 2. Procédé selon la revendication 1, caractérisé en ce que l'étape de localisation des éléments (P) est réalisée à partir d'une image brute numérique (R) telle qu'elle a été enregistrée par un détecteur numérique de l'appareil d'imagerie. 3. Procédé selon la revendication 2, caractérisé en ce que l'étape de localisation des éléments (P) comprend une étape (100) consistant à convertir l'image brute (R) en une image des contrastes (C) par une normalisation locale du contraste. 4. Procédé selon l'une des revendications 1 à 4, caractérisé en ce que l'étape de localisation des éléments (P) comprend: - une étape (200) consistant à corréler, pour chacun des points (x,y) d'un ensemble de points de l'image radiologique, l'intensité de l'image dans un voisinage du point avec chacune des fonctions (gs,a) d'une série de fonctions filtres, pour produire des image filtrées, - une étape (300) consistant à sélectionner à partir des images 25 filtrées des éléments (P) susceptibles de constituer des signes de microcalcification. 6. Procédé selon la revendication 4, caractérisé en ce que les fonctions filtres (gs,a) sont des fonctions de type 3-spline. 7. Procédé selon l'une des revendications 4 ou 5, caractérisé en ce 30 que l'étape (300) de sélection des particules comprend: - une étape (310) consistant à déterminer une image filtrée isotrope (B) définie par un ensemble de points (x,y), l'intensité de chacun de ces points étant la moyenne des intensités des points correspondants appartenant aux images filtrées, - une étape (312) consistant à en déduire par application à l'image filtrée isotrope (B) obtenue d'un seuil et d'un algorithme d'étiquetage des composantes connexes, un ensemble d'éléments (P) susceptibles de constituer des signes de microcalcification. 7. Procédé selon l'une des revendications qui précèdent, caractérisé en ce qu'il comprend: - une étape (320) consistant à déterminer une image filtrée 10 directionnellement (0), en corrélant l'ensemble de l'image radiologique avec une fonction filtre conduisant à une réponse maximale, une étape (322) consistant à en déduire par application à l'image filtrée directionnellement (0) d'un seuil et d'un algorithme d'étiquetage des composantes connexes, un ensemble de structures (F) susceptibles de constituer des fibres, - une étape (323) consistant à déterminer pour chaque point d'une structure (F) son degré d'appartenance à une classe Fibre . 8. Procédé selon l'une des revendications 1 à 7, caractérisé en ce que le rehaussement de l'intensité des éléments (P) est réalisé sur une 20 image de présentation (S) visualisable par le radiologue. 9. Procédé selon l'une des revendications 1 à 8, caractérisé en ce que le rehaussement de l'intensité des éléments comprend l'étape consistant à, pour chaque point (x,y) de l'image radiologique appartenant à une particule (P) présentant une intensité donnée, appliquer à cette intensité un facteur d'amplification. 10. Procédé selon la revendication 9, caractérisé en ce que le rehaussement de l'intensité des éléments susceptibles de constituer des signes de microcalcification comprend l'étape (500) consistant à appliquer à chaque point (x,y) de l'image radiologique appartenant à un élément (P) un facteur d'amplification de la différence de niveau de gris entre ledit point (x,y) et un ensemble de points (x+k,j+l) situés dans un voisinage du point (x,y) et n'appartenant pas à l'élément (P). 11. Procédé selon l'une des revendications 9 ou 10, caractérisé en ce que le facteur d'amplification appliqué aux points appartenant à une particule (P) est fonction d'un paramètre de similarité (p.calcium(P)) de ladite particule (P) avec un signal calcique. 12. Procédé selon la revendication 11, caractérisé en ce que pour chaque particule, le paramètre de similarité (p.calcium(P)) est déterminé selon les étapes suivantes: a) pour chacun des points (x,y) formant un élément (P): - déterminer une amplitude maximale (Amax) en corrélant, l'intensité 10 de l'image radiologique dans un voisinage (W) du point (x,y) avec une fonction filtre (g6,8) conduisant à une réponse maximale, - déterminer une amplitude maximale en direction orthogonale (Amax') en corrélant, l'intensité de l'image radiologique dans un voisinage (W) du point (x,y) avec la fonction filtre (g,,01) de direction orthogonale à la fonction filtre (g6,0) conduisant à une réponse maximale, b) pour chaque élément (P): - déterminer un paramètre de brillance (Ap) comme une moyenne des amplitudes maximales (Amax) sur l'ensemble des points (x,y) de l'élément (P), - déterminer un paramètre de brillance et de rondeur (Ap') comme étant une moyenne des amplitudes maximales en direction orthogonale (Amaxl) sur l'ensemble des points (x,y) de l'élément (P), - en déduire un paramètre de similarité ( ca,c;um(P)) de l'élément (P) à un signal calcique en fonction des paramètres de brillance et/ou de rondeur (Ap, Ap') ainsi que du degré d'appartenance du point central ()cc, yc) de l'élément (P) à une classe Fibre . 13. Procédé selon l'une des revendications 1 à 14, caractérisé en ce qu'il comprend une étape (600) consistant à détecter des clusters de microcalcifications et à disposer des marqueurs à proximité ou autour de ces clusters. 14. Dispositif de traitement d'image radiologique caractérisé en ce qu'il comprend des moyens pour localiser sur une image radiologique (R) fournie par un appareil d'imagerie des éléments (P) susceptibles de constituer des signes de microcalcification et des moyens pour afficher une image (S) sur laquelle l'intensité des éléments (P) ainsi localisées a été rehaussée.