Brevet FR2954556A1 — publié le 20110624

DOMAINE TECHNIQUE GENERAL La présente invention concerne le traitement et l'affichage d'images de tomosynthèse et d'images de radiographie d'un objet d'intérêt. Plus précisément, elle concerne un procédé de traitement d'images de tomosynthèse d'un objet d'intérêt pour obtenir des coupes de reconstruction tridimensionnelle (3D) de l'objet d'intérêt ainsi qu'une image bidimensionnelle (2D) de l'objet d'intérêt. L'invention concerne également un procédé d'affichage desdites coupes et de ladite image afin de commuter rapidement de l'image 2D vers les coupes de reconstruction et inversement et des systèmes et produits de programmes d'ordinateur pour la mise en oeuvre de tels procédés. ETAT DE L'ART La radiographie est généralement utilisée pour la recherche d'anormalités dans un objet d'intérêt. Nous nous intéresserons en particulier à la mammographie qui est utilisée depuis des décennies pour le dépistage et le diagnostique du cancer du sein mais l'invention peut s'appliquer également à la radiologie conventionnelle. Une image de radiographie représente une projection d'un objet d'intérêt. Elle est généralement obtenue en plaçant l'objet entre une source émettrice de rayons X et un détecteur de rayons X, de sorte que les rayons atteignent le détecteur après avoir traversé l'objet. L'image de radiographie est ensuite construite à partir des données fournies par le détecteur, et représentent l'objet d'intérêt projeté sur le détecteur dans la direction des rayons X. Sur cette image, un praticien radiologue expérimenté peut distinguer des signes radiologiques indicateurs d'un problème potentiel, par exemple des microcalcifications ou des opacités dans le cas de la mammographie. Toutefois, dans une image de projection 2D, la superposition des tissus peut masquer des lésions, et en aucun cas leur position réelle dans l'objet d'intérêt n'est connue, le praticien n'ayant aucune information sur la position du signe radiologique dans la direction de projection. La tomosynthèse s'est récemment développée pour répondre à ces problèmes ; elle permet d'obtenir une représentation 3D d'un objet d'intérêt sous la forme d'une série de coupes successives. Ces coupes sont reconstruites à partir de projections de l'objet d'intérêt sous divers angles. Pour ce faire, l'objet d'intérêt est généralement placé entre une source émettrice de rayons X et un détecteur de rayons X. La source et/ou le détecteur sont mobiles, de sorte que la direction de projection de l'objet sur le détecteur peut varier (typiquement sur une plage angulaire de 300). On obtient ainsi plusieurs projections de l'objet d'intérêt sous différents angles, à partir desquelles on peut reconstruire une représentation tridimensionnelle de l'objet d'intérêt, généralement par une méthode de reconstruction bien connue de l'homme du métier. Pour chaque projection, les doses de radiation des rayons sont naturellement inférieures à celles utilisées pour une mammographie standard. Par exemple, en notant D la dose de radiation d'une mammographie standard, et N le nombre de projections utilisées pour une tomosynthèse, la dose de radiation utilisée pour chaque projection est de l'ordre de D/N. Ainsi, en parcourant les coupes de tomosynthèse d'un objet d'intérêt, 20 un praticien peut détecter des signes radiologiques dans l'objet d'intérêt et évaluer leur position 3D. Cependant, les praticiens ont des années d'expérience dans l'analyse d'images de mammographies, alors que l'analyse de coupes de tomosynthèse est encore à son commencement. 25 Ainsi, une période de transition s'avère préférable, période au cours de laquelle les praticiens seront amenés à analyser des coupes de tomosynthèse d'un objet d'intérêt accompagnées d'images de mammographie standard d'un même objet d'intérêt, dans le but d'acquérir une expérience et une aisance dans l'analyse des coupes de représentation 30 3D et de comparer à des acquisition 2D précédentes afin d'apprécier l'évolution temporelle des signes radiologiques. Dans cette optique, des recherches ont été menées pour réaliser des systèmes d'acquisition d'images de radiographie permettant d'obtenir des images de mammographie 2D et des coupes de tomosynthèse d'un même objet d'intérêt. Le document EP 1 795 569 propose un exemple d'un tel système, comprenant une source émettrice de rayons X à dose de radiation variable, de sorte que la dose de radiation émise puisse être adaptée soit à des projections de tomosynthèse, soit à une mammographie. Toutefois, ce système émet plus de radiations qu'un système d'acquisition de mammographie, ou qu'un système d'acquisition de tomosynthèse ; il en émet environ autant que les deux réunis. Aussi, il existe un réel besoin de disposer d'un système d'acquisition d'images de radiographie permettant d'obtenir des images de mammographie 2D et des coupes de tomosynthèse d'un même objet d'intérêt, sans que la dose de radiation émise pour ces acquisitions n'excède celle des systèmes de l'art antérieur, en particulier celle des systèmes de mammographie 2D. Deux méthodes pour reconstruire une image 2D ressemblant à une image de mammographie standard ont été testées par la Demanderesse. Ces méthodes présentent toutefois des inconvénients. Une première méthode consiste à prendre tout simplement pour image 2D une projection de tomosynthèse, préférentiellement à angle nul, c'est-à-dire lorsque le tube est orthogonal au détecteur. Cette méthode ne peut constituer une réelle solution au problème exposé, dans le sens où la dose de radiation émise pour une projection de tomosynthèse est trop faible ; l'information récoltée sur le détecteur est alors noyée dans le bruit quantique et ne permet pas de repérer les objets faiblement contrastés particulièrement les petites microcalcifications. Une deuxième méthode consiste à regrouper les coupes reconstruites en une couche épaisse (représentative de toutes les coupes successives), en effectuant pour chaque pixel de la couche épaisse une moyenne des voxels situés le long du chemin menant du pixel à la source (placée perpendiculairement au détecteur) ou en sélectionnant le voxel d'intensité maximum (méthode dite MIP pour « Maximum Intensity Projection » en terminologie anglo-saxonne). Cette méthode ne donne pas de résultats satisfaisants à cause des artefacts de reconstruction générés particulièrement par la retroprojection des gros objets. Les petits objets ne génèrent qu'une quantité limitée d'artefacts dans le volume. C'est la raison pour laquelle il est proposé d'enlever les gros objets des projections en leur appliquant un filtre passe haut avant d'appliquer la retroprojection. Les gros objets seront apportés par la projection centrale car ils sont toujours visibles même à faible dose. PRESENTATION DE L'INVENTION L'invention permet d'obtenir une reconstruction 3D d'un objet d'intérêt ainsi qu'une image 2D de l'objet d'intérêt similaire à une image de mammographie standard à partir d'acquisition de projection de tomosynthèse de l'objet d'intérêt, en s'affranchissant des inconvénients précités. Selon un premier aspect, l'invention concerne, un procédé de traitement d'images de tomosynthèse d'un objet d'intérêt, au moyen d'un système d'imagerie, le système d'imagerie comprenant une source à rayons X disposée en regard d'un détecteur sur lequel l'objet d'intérêt est disposé. Le procédé selon le premier aspect de l'invention comprend les étapes suivantes : acquisition d'une pluralité d'images de projection 2D de l'objet d'intérêt selon une pluralité d'orientations repérées par rapport à la perpendiculaire au détecteur, une orientation dite nulle étant la plus proche de la perpendiculaire ; application, d'un filtre sur les images de projections 2D acquises de manière à obtenir des images de projection filtrées de l'objet d'intérêt ; détermination de coupes de reconstruction de l'objet d'intérêt à partir d'au moins deux des projections filtrées, l'ensemble des coupes de reconstruction étant le volume reconstruit filtré de l'objet d'intérêt ; reprojection, des coupes de reconstruction selon une direction d'orientation déterminée, de préférence dans la direction d'orientation dite nulle, de manière à obtenir une image 2D intermédiaire de l'objet d'intérêt ; et obtention, d'une image 2D finale de l'objet d'intérêt par combinaison de l'image 2D intermédiaire et de l'image de projection correspondant à l'orientation déterminée. Le procédé selon le premier aspect de l'invention est avantageusement complété par les caractéristiques suivantes, prises seules ou en une quelconque de leurs combinaisons techniquement possibles : le filtre appliqué aux images de projection est un filtre passe haut, dont les caractéristiques dépendent de préférence de l'épaisseur de l'objet d'intérêt ; les coupes de reconstruction sont déterminées par un algorithme de rétroprojection, de préférence non linéaire ; la reprojection des coupes de reconstruction est une reprojection SIP dans la direction d'orientation déterminée, la reprojection SIP consistant, pour chaque pixel de l'image 2D intermédiaire, en une détermination, dans le volume reconstruit filtré, d'un voxel d'intensité calculé en utilisant un tri des valeurs de pixels le long du rayon allant de la source au pixel de l'image 2D intermédiaire et en une mémorisation, dans l'unité de mémoire du système d'imagerie, d'un identifiant de la coupe de reconstruction dans laquelle se trouve le voxel issu du tri, de manière à disposer, dans l'unité de mémoire, d'une information de profondeur reliant chaque pixel de l'image 2D intermédiaire à la coupe de reconstruction associée dont provient ce pixel ; la reprojection des de reconstruction est une reprojection MIP dans la direction d'orientation déterminée, la reprojection MIP consistant, pour chaque pixel de l'image 2D intermédiaire, en une détermination, dans le volume reconstruit filtré, du voxel d'intensité maximale le long du rayon allant de la source au pixel de l'image 2D intermédiaire et en une mémorisation, dans l'unité de mémoire du système d'imagerie, d'un identifiant de la coupe de reconstruction dans laquelle se trouve le voxel d'intensité maximale, de manière à disposer, dans l'unité de mémoire, d'une information de profondeur reliant chaque pixel de l'image 2D intermédiaire à la coupe de reconstruction associée dont provient ce pixel ; le procédé comprend une étape de détermination d'un volume, obtenu à l'aide d'un système de détection automatique de lésions, dans lequel chaque voxel représente une probabilité d'appartenir à une lésion ; la reprojection des coupes de reconstruction consiste en tout pixel de l'image 2D intermédiaire en la sélection du voxel de plus forte probabilité d'appartenir à une lésion le long du rayon allant de la source à ce pixel, et en une mémorisation dans l'unité de mémoire du système d'imagerie, d'un identifiant de la coupe de reconstruction dans laquelle se trouve le voxel de probabilité maximale, de manière à disposer, dans l'unité de mémoire, d'une information de profondeur reliant chaque pixel de l'image 2D intermédiaire à la coupe de reconstruction associée dont provient ce pixel ; le procédé comprend en outre une étape de régularisation locale, par l'unité de traitement, de l'information de profondeur, de manière à obtenir une information de profondeur régularisée ; le procédé comprend préalablement à l'étape de détermination de l'image 2D finale, une étape d'application d'un filtre à l'image de projection correspondant à l'orientation déterminée de manière à obtenir une image de projection correspondant à l'orientation déterminée, la détermination de l'image 2D finale consistant en la combinaison de l'image 2D intermédiaire et de l'image de projection correspondant à l'orientation déterminée filtrée et le filtre appliqué à l'image de projection correspondant à l'orientation déterminée de l'objet d'intérêt peut être un filtre passe-bas ; l'image 2D finale est déterminée pour chaque orientation ; chaque image ainsi obtenue étant destinée à être affichée l'une après l'autre ou bien en défilement automatique ; le procédé comprend une détermination de coupes d'affichage obtenues à partir des images de projection 2D, l'ensemble des coupes d'affichage étant le volume reconstruit de l'objet d'intérêt. Selon un second aspect, l'invention concerne un procédé d'affichage d'un ensemble d'images de radiographie d'un objet d'intérêt par une unité d'affichage d'images de radiographie, les images radiographiques étant les coupes d'affichage de l'objet et une image 2D finale de l'objet obtenues par un procédé de traitement d'images de tomosynthèse selon le premier aspect de l'invention, l'information de profondeur étant stockée dans l'unité mémoire. Le procédé d'affichage selon le second aspect de l'invention comprend les étapes suivantes : affichage de l'image 2D finale au moyen de l'unité d'affichage ; sélection d'un pixel de l'image 2D finale ; recherche, parmi, l'information de profondeur stockée dans l'unité de mémoire, d'un identifiant de la coupe de reconstruction associée au pixel sélectionné ; ajustement de la taille de la coupe d'affichage associée au pixel sélectionné à la taille de l'image 2D finale ; affichage, au moyen de l'unité d'affichage, de la coupe d'affichage ajustée et associée au pixel sélectionné. Le procédé selon le second aspect de l'invention peut comprendre à l'issu de l'affichage de la coupe d'affichage ajustée une étape de retour à un affichage de l'image 2D finale. Selon un troisième aspect, l'invention concerne un procédé d'affichage d'un ensemble d'images de radiographie d'un objet d'intérêt par une unité d'affichage d'images de radiographie, les images radiographiques étant les coupes d'affichage de l'objet et une image 2D finale de l'objet obtenues par un procédé de traitement d'images de tomosynthèse selon le premier aspect de l'invention. Le procédé d'affichage selon le troisième aspect comprend les étapes suivantes : affichage d'une coupe d'affichage au moyen de l'unité d'affichage ; ajustement de la taille de l'image 2D finale à la taille de la coupe d'affichage ; affichage sur demande de l'utilisateur, au moyen de l'unité d'affichage, de l'image 2D finale ajustée. Le procédé d'affichage selon le troisième aspect de l'invention peut comprendre à l'issu de l'affichage de l'image 2D finale ajustée une étape de 5 retour à un affichage de la coupe d'affichage. Selon un quatrième aspect, l'invention concerne un système d'imagerie comprenant un détecteur de rayons X, une source émettrice de rayons X, une unité de traitement et des moyens de mise en oeuvre d'un 10 procédé de traitement selon le premier aspect de l'invention. Selon un cinquième aspect, l'invention concerne un système d'affichage destiné à être mis en oeuvre dans une unité de traitement d'un système d'imagerie médicale selon le quatrième aspect de l'invention, 15 l'unité de traitement comprenant des moyens de mise en oeuvre d'un procédé d'affichage selon le second et/ou le troisième aspect de l'invention. Selon un sixième aspect, l'invention concerne un produit de programme d'ordinateur comprenant des instructions de code enregistrées 20 sur un support utilisable dans un ordinateur, caractérisé en ce qu'il comprend des instructions pour la mise en oeuvre du procédé de traitement selon le premier aspect de l'invention. Selon un septième aspect, l'invention concerne un produit de 25 programme d'ordinateur comprenant des instructions de code enregistrées sur un support utilisable dans un ordinateur, caractérisé en ce qu'il comprend des instructions pour la mise en oeuvre d'un procédé d'affichage selon le second et/ou le troisième aspect de l'invention. 30 L'invention permet d'obtenir, à partir d'images de projection de tomosynthèse d'un objet d'intérêt, une image 2D finale similaire à une image de mammographie standard de l'objet d'intérêt, image 2D sur laquelle les signes radiologiques de petite taille, comme les microcalcifications et les fibres, sont plus visibles qu'avec les procédés connus. Par ailleurs, le procédé d'affichage selon le deuxième ou le troisième aspect de l'invention permet à un praticien utilisateur d'accéder très rapidement à une coupe de reconstruction d'un objet d'intérêt sur laquelle est situé un signe radiologique repéré sur l'image 2D finale. Ainsi, un praticien radiologue peut examiner avec rapidité et efficacité les coupes de reconstruction d'un objet d'intérêt à l'aide de l'image similaire à une image de mammographie standard. PRESENTATION DES FIGURES D'autres caractéristiques, buts et avantages de l'invention ressortiront de la description qui suit, qui est purement illustrative et non limitative, et qui doit être lue en regard des dessins annexés sur lesquels : - la figure 1 illustre schématiquement un système d'imagerie selon l'invention ; - la figure 2 illustre schématiquement quelques fonctionnalités du système d'imagerie ; - la figure 3 illustre, sous la forme d'un schéma bloc, un procédé de 20 traitement d'images selon l'invention ; - la figure 4 illustre, sous la forme d'un schéma bloc, une étape de reprojection de coupes de reconstruction d'un objet d'intérêt selon un mode de réalisation de l'invention ; - la figure 5 illustre, sous la forme d'un schéma bloc, une étape de 25 reprojection de coupes de reconstruction selon un mode de réalisation de l'invention ; - la figure 6 illustre, sous la forme d'un schéma bloc, un procédé d'affichage selon l'invention. DESCRIPTION DETAILLEE Système d'imagerie médicale La figure 1 illustre, schématiquement, un système d'imagerie 5 médicale 10 pour l'acquisition d'images 2D de projection pour la reconstruction d'une image 3D d'une région d'intérêt. Le système d'imagerie médical 10 comprend une unité d'acquisition d'images 12, une unité de traitement d'images 13 et un système d'affichage 15. 10 La figure 2 illustre schématiquement l'unité d'acquisition 12. L'unité d'acquisition 12 permet d'acquérir une pluralité de projection 2D d'une région d'intérêt O û d'un organe û d'un patient. L'unité d'acquisition 12 est notamment constituée par un détecteur 11 situé en regard d'une source 12 de rayons X. Le détecteur 11 est par exemple une 15 caméra numérique. L'unité d'acquisition 12 est par exemple un système d'acquisition par rayons X, ce dernier comprenant tout moyen connu permettant l'émission de rayons X sur l'objet O et l'acquisition d'images résultantes. L'unité d'affichage 15 peut être intégrée dans l'unité d'acquisition 20 d'images 12 ou l'unité de traitement d'images 13, ou être séparée de l'unité d'acquisition 12 et de l'unité de traitement 13. L'unité d'affichage 15 est par exemple un écran d'ordinateur, un moniteur, un écran plat, un écran plasma ou tout type de dispositif d'affichage connu du commerce. L'unité d'affichage 15 permet à un praticien radiologue de contrôler la 25 reconstruction et/ou l'affichage des images 2D acquises. L'unité de traitement 13 est adaptée à la mise en oeuvre de procédé de traitement (par exemple de reconstruction d'une image 3D à partir d'images 2D). L'unité de traitement 13 peut être intégrée dans l'unité d'acquisition 12 ou être séparée de l'unité d'acquisition 12. L'unité de 30 traitement 13 est par exemple un/des ordinateur(s), un/des processeur(s), un/des microcontrôleur(s), un/des micro-ordinateur(s), un/des automate(s) programmable(s), un/des circuit(s) intégré(s) spécifique(s) d'application, 10 d'autres circuits programmables, ou d'autres dispositifs qui incluent un ordinateur tel qu'une station de travail. L'unité de traitement 13 est couplée à des moyens mémoires 14 qui peuvent être intégrés ou séparés de l'unité de traitement 13. Ces moyens peuvent être formés par un disque dur ou tout autre moyen de stockage amovible (un cd-rom, une disquette etc.). Ces moyens mémoires peuvent servir à stocker une image 3D de la zone de l'organe visualisée comme une image 2D acquise ou traitée. Il peut s'agir d'une mémoire ROM/RAM de l'unité de traitement 13, un CD-ROM, une clé USB, une mémoire d'un serveur central. L'unité de traitement 13 peut comprendre un dispositif de lecture (non représenté) par exemple un lecteur de disquettes ou un lecteur de CD-ROM, pour lire les instructions du procédé de traitement (qui va être décrit dans la suite) d'un support d'instructions (non montré), comme une disquette ou un CD-ROM. En variante, l'unité de traitement 13 exécute les instructions du procédé de traitement (qui va être décrit dans la suite) stockées dans des micrologiciels (non représentés). Procédé de traitement d'images Le procédé de traitement d'images que nous allons décrire permet 20 d'obtenir, à partir d'images de tomosynthèse, une image similaire à une image mammographique. On précise ici que l'on entend par image similaire à une image mammographique, une image qui ressemble visuellement à une image de mammographie 2D standard pleine dose. 25 Le procédé de traitement d'images consiste donc à traiter des images de radiographie obtenues par un système d'imagerie 10 comprenant une source à rayons X 12 disposée en regard d'un détecteur 11 sur lequel l'objet d'intérêt O est disposé. La figure 3 illustre des étapes du procédé de traitement d'images 30 obtenues par tomosynthèse. Dans une première étape S1 on acquiert une pluralité d'images de projection 2D 30 de l'objet d'intérêt O selon une pluralité d'orientations, une orientation dite nulle étant la plus proche de la perpendiculaire H. On précise ici que l'on repère chaque orientation par rapport à la perpendiculaire H au détecteur 11. Au cours de cette première étape on acquiert, en particulier, une image 2D selon une orientation déterminée 31. L'orientation déterminée est 5 de préférence l'orientation dite nulle. Le procédé comprend ensuite une étape S2 d'application d'un filtre sur les images de projections 2D acquises 30 de manière à obtenir des images de projection filtrées 40 de l'objet d'intérêt O. Ce filtre est de préférence du type passe-haut et sa fréquence de 10 coupure est de préférence fonction de l'épaisseur de l'objet d'intérêt O. Au cours d'une étape S3, on détermine des coupes 50 de reconstruction de l'objet d'intérêt O. Cette détermination S3 consiste, en particulier, en une rétroprojection des images de projection 2D filtrées 40. Cette rétroprojection peut être en particulier non linéaire de type 15 « Order Statistics Based Backprojection ». Dans la rétroprojection linéaire, chaque voxel du volume est reconstruit en utilisant N pixel d'information, chaque pixel étant déterminé par reprojection du voxel dans chacune des N projections. Dans la rétroprojection non-linéaire, le pixel d'intensité maximale parmi les N n'est pas utilisé ce qui permet de réduire nettement 20 les artéfacts de réplication provoqués par les objets les plus intenses. On note que les coupes de reconstruction de l'objet d'intérêt O représentent le volume reconstruit de l'objet d'intérêt O. Ensuite au cours d'une étape S4 on effectue une reprojection des coupes 50 de reconstruction dans la direction déterminée. Ceci permet 25 d'obtenir une image 2D intermédiaire 60 de l'objet d'intérêt O. On note que l'on reprojette selon la même direction que l'image de projection correspondant à l'orientation déterminée 31. Enfin, au cours d'une étape S5, on obtient une image 2D finale de l'objet d'intérêt par combinaison de l'image 2D intermédiaire 60 et de 30 l'image de projection correspondant à l'orientation déterminée 31. La combinaison est de préférence une combinaison linéaire pixel à pixel. L'image 2D finale étant une image similaire à une image de mammographie. Dans une alternative avantageuse, l'image 2D finale de l'objet d'intérêt est déterminée pour chaque orientation. L'étape de reprojection S4 des coupes 50 de reconstruction est, dans une implémentation préférée, une reprojection MIP (en anglais « Maximum Intensity Pixel ») dans la direction d'orientation déterminée. De façon plus générale, toute reprojection utilisant un tri des valeurs des pixels présents le long des rayons peut être utilisée (SIP, pour « Sorted Intensity Pixel »). Le tri consiste au classement selon leur intensité des pixels (tri croissant ou décroissant). La figure 4 illustre les étapes de la reprojection MIP selon un premier mode de réalisation. Ce type de reprojection MIP consiste, pour chaque pixel de l'image 2D intermédiaire 60, en une détermination S41, dans le volume reconstruit de l'objet d'intérêt, du voxel d'intensité maximale le long du rayon allant de la source à ce pixel et en une mémorisation S42, dans une unité de mémoire 14 du système d'imagerie, d'un identifiant de la coupe de reconstruction dans laquelle se trouve le voxel d'intensité maximale. De cette façon, on dispose dans l'unité de mémoire 14, d'une information argMax de profondeur reliant chaque pixel de l'image 2D intermédiaire 60 à la coupe de reconstruction associée dont provient ce pixel. La reprojection peut être mise en oeuvre d'une manière différente de celle présentée ci-dessus. La figure 5 illustre les étapes de la reprojection selon un second 25 mode de réalisation. Dans ce mode de réalisation, la reprojection S4 des coupes 50 de reconstruction consiste, en outre, en tout pixel de l'image 2D intermédiaire 60 en une sélection S43 du voxel ayant la plus forte probabilité d'appartenir à une lésion le long du rayon allant de la source à ce pixel et en une 30 mémorisation S44, dans une unité de mémoire 14 du système d'imagerie, d'un identifiant de la coupe de reconstruction dans laquelle se trouve le voxel de probabilité maximale. Ceci suppose qu'à chaque voxel soit associé une probabilité d'appartenir à une lésion. Un système de détection automatique (CAD 3D pour en anglais, « Computer Aided Detection ») permet d'obtenir un tel volume de probabilités. De cette façon, on dispose dans l'unité de mémoire 14, d'une information argMax de profondeur reliant chaque pixel de l'image 2D intermédiaire 60 à la coupe de reconstruction associée dont provient ce pixel. Le procédé de traitement peut comprendre en outre une étape de régularisation locale S6, par l'unité de traitement 13, de l'information argMax de profondeur. Cette étape de régularisation consiste à rendre localement plus uniforme l'information argMax. Le résultat est une information de profondeur notée argMax'. De manière avantageuse, le procédé comprend une étape d'application S2', par l'unité de traitement 13, d'un filtre à l'image de projection d'orientation déterminée 31 de l'objet d'intérêt O préalablement à l'élaboration S5 de l'image 2D finale 70 de manière à réduire le bruit. Le filtre appliqué est de préférence un filtre passe-bas. En outre, il est possible de mettre en oeuvre une étape de détermination S4', de coupes d'affichage 50' qui correspondent au volume reconstruit de l'objet d'intérêt. Il s'agit en d'autres termes d'un volume que l'on obtient à partir des images de projection avec les procédés de reconstruction connus de l'état de l'art dont l'objectif est la visualisation de coupes. Procédé d'affichage Le procédé de traitement ci-dessus décrit trouve une application particulière dans un procédé d'affichage d'images de radiographie d'un objet d'intérêt O. L'affichage des images de radiographie est effectué au moyen d'une unité d'affichage 15 connectée à l'unité de traitement 13 mettant en oeuvre le procédé d'affichage. La figure 5 illustre des étapes du procédé d'affichage. Le procédé d'affichage comprend une étape Dl d'affichage de l'image 2D finale 70 obtenue au moyen du procédé de traitement d'images ci-dessus décrit. Il comprend en outre une sélection d'un pixel de l'image 2D finale 70. Le procédé d'affichage comprend également une recherche D3, parmi, l'information argMax ou argMax' de profondeur stockée dans l'unité de mémoire 14, d'un identifiant de la coupe de reconstruction associée au pixel sélectionné. Suite à cette recherche, le procédé comprend une étape d'ajustement D4 de la taille de la coupe d'affichage de même identifiant à la taille de l'image 2D finale 70, le facteur d'ajustement dépendant préférentiellement de la hauteur de la coupe d'affichage. Enfin, on affiche, au moyen de l'unité d'affichage, la coupe d'affichage ainsi ajustée. Par ce biais, la coupe d'affichage et l'image 2D finale 70 ont les mêmes propriétés géométriques avec des structures placées au même endroit. Il est alors plus facile de commuter entre les coupes d'affichage et l'image 2D finale 70 et inversement. En variante, le procédé d'affichage peut consister en la succession des étapes suivantes : affichage Dl d'une coupe d'affichage au moyen de l'unité d'affichage 15; ajustement D4 de la taille de l'image 2D finale 70 à la taille de la coupe d'affichage, le facteur d'ajustement dépendant préférentiellement de la hauteur de la coupe d'affichage; affichage D5, au moyen de l'unité d'affichage 15, de l'image 2D finale 70 ajustée. Par ce biais, la coupe d'affichage et l'image 2D finale 70 ont les mêmes propriétés géométriques avec des structures placées au même endroit. Il est alors plus facile de commuter entre les coupes d'affichage et l'image 2D finale 70 et inversement. REVENDICATIONS1. Procédé de traitement d'images de tomosynthèse d'un objet d'intérêt (0), au moyen d'un système d'imagerie (10), le système d'imagerie comprenant une source à rayons X (12) disposée en regard d'un détecteur (11) sur lequel l'objet d'intérêt (0) est disposé, le procédé comprenant les étapes suivantes : - acquisition (SI) d'une pluralité d'images de projection 2D (30) de l'objet d'intérêt (0) selon une pluralité d'orientations repérées par rapport à la perpendiculaire (H) au détecteur (11), une orientation dite nulle étant la plus proche de la perpendiculaire (H) ; - application (S2), d'un filtre sur les images de projections 2D acquises (30) de manière à obtenir des images de projection filtrées (40) de l'objet d'intérêt (0) ; - détermination (S3) de coupes (50) de reconstruction de l'objet d'intérêt (0) à partir d'au moins deux des projections filtrées (40), l'ensemble des coupes (50) de reconstruction étant le volume reconstruit filtré de l'objet d'intérêt (0) ; - reprojection (S4), des coupes (50) de reconstruction selon une direction d'orientation déterminée (31), la direction déterminée étant de préférence la direction d'orientation dite nulle, de manière à obtenir une image 2D intermédiaire (60) de l'objet d'intérêt (0) ; et - obtention (S5), d'une image 2D finale (70) de l'objet d'intérêt (0) par combinaison de l'image 2D intermédiaire (60) et de l'image de projection correspondant à l'orientation déterminée (31). 2. Procédé selon la revendication 1, dans lequel le filtre appliqué aux images de projection (30) est un filtre passe haut, dont les caractéristiques dépendent de préférence de l'épaisseur de l'objet d'intérêt (0). 3. Procédé selon l'une des revendications 1 à 2, dans lequel les coupes (50) de reconstruction sont déterminées par un algorithme de rétroprojection, de préférence non linéaire.30 4. Procédé selon l'une des revendication 1 à 3, dans lequel la reprojection (S4) des coupes (50) de reconstruction est une reprojection SIP dans la direction d'orientation déterminée (31), la reprojection SIP consistant, pour chaque pixel de l'image 2D intermédiaire (60), en une détermination (S41), dans le volume reconstruit filtré, d'un voxel d'intensité calculé en utilisant un tri des valeurs de pixels le long du rayon allant de la source au pixel de l'image 2D intermédiaire (60) et en une mémorisation (S42), dans l'unité de mémoire (14) du système d'imagerie, d'un identifiant de la coupe de reconstruction dans laquelle se trouve le voxel issu du tri, de manière à disposer, dans l'unité de mémoire (14), d'une information (argMax) de profondeur reliant chaque pixel de l'image 2D intermédiaire (60) à la coupe de reconstruction associée dont provient ce pixel. 5. Procédé selon l'une des revendications 1 à 3, dans lequel la reprojection (S4) des coupes (50) de reconstruction est une reprojection MIP dans la direction d'orientation déterminée (31), la reprojection MIP consistant, pour chaque pixel de l'image 2D intermédiaire (60), en une détermination (S41), dans le volume reconstruit filtré, du voxel d'intensité maximale le long du rayon allant de la source au pixel de l'image 2D intermédiaire (60) et en une mémorisation (S42), dans l'unité de mémoire (14) du système d'imagerie, d'un identifiant de la coupe de reconstruction dans laquelle se trouve le voxel d'intensité maximale, de manière à disposer, dans l'unité de mémoire (14), d'une information (argMax) de profondeur reliant chaque pixel de l'image 2D intermédiaire (60) à la coupe de reconstruction associée dont provient ce pixel. 6. Procédé selon l'une des revendications 1 à 5, comprenant une étape de détermination (S3') d'un volume, obtenu à l'aide d'un système de détection automatique de lésions, dans lequel chaque voxel représente une probabilité d'appartenir à une lésion. 7. Procédé selon la revendication 6, dans lequel la reprojection (S4) des coupes (50) de reconstruction consiste en tout pixel de l'image 2Dintermédiaire (60) en la sélection (S43) du voxel de plus forte probabilité d'appartenir à une lésion le long du rayon allant de la source à ce pixel, et en une mémorisation (S44), dans l'unité de mémoire (14) du système d'imagerie, d'un identifiant de la coupe de reconstruction dans laquelle se trouve le voxel de probabilité maximale, de manière à disposer, dans l'unité de mémoire (14), d'une information (argMax) de profondeur reliant chaque pixel de l'image 2D intermédiaire (60) à la coupe de reconstruction associée dont provient ce pixel. 8. Procédé selon l'une des revendications 4 à 7, comprenant en outre une étape de régularisation locale (S6), par l'unité de traitement (13), de l'information (argMax) de profondeur, de manière à obtenir une information (argMax') de profondeur régularisée. 9. Procédé selon les revendications 1 à 8, comprenant préalablement à l'étape de détermination (S5) de l'image 2D finale (70), une étape d'application (S2') d'un filtre à l'image de projection correspondant à l'orientation déterminée (31) de manière à obtenir une image de projection correspondant à l'orientation déterminée filtrée, la détermination de l'image 2D finale (70) consistant en la combinaison de l'image 2D intermédiaire (60) et de l'image de projection correspondant à l'orientation déterminée filtrée. 10. Procédé selon la revendication 9 dans lequel le filtre appliqué à l'image de projection correspondant à l'orientation déterminée (31) de 25 l'objet d'intérêt (0) est un filtre passe-bas. 11. Procédé selon les revendications 1 à 10 dans lequel l'image 2D finale (70) est déterminée pour chaque orientation ; chaque image ainsi obtenue étant destinée à être affichée l'une après l'autre ou bien en 30 défilement automatique. 12. Procédé selon l'une des revendications précédentes comprenant une détermination (S4') de coupes d'affichage (50') obtenues à partir desimages de projection 2D (30), l'ensemble des coupes d'affichage (50') étant le volume reconstruit de l'objet d'intérêt. 13. Procédé selon l'une des revendications 9 à 12 comprenant en 5 outre les étapes suivantes : - affichage (Dl) de l'image 2D finale (70) au moyen de l'unité d'affichage (15), sélection (D2) d'un pixel de l'image 2D finale (70) ; - recherche (D3), parmi, l'information (argMax') de profondeur stockée 10 dans l'unité de mémoire (14), d'un identifiant de la coupe de reconstruction associée au pixel sélectionné ; ajustement (D4) de la taille de la coupe d'affichage (50') associée au pixel sélectionné à la taille de l'image 2D finale (70) ; - affichage (D5), au moyen de l'unité d'affichage (15), de la coupe 15 d'affichage (50') ajustée et associée au pixel sélectionné. 14. Procédé selon la revendication 13, comprenant à l'issu de l'affichage de la coupe d'affichage ajustée une étape de retour (D6) à un affichage de l'image 2D finale (70). 20 15. Procédé selon l'une des revendications 9 à 12 comprenant les étapes suivantes affichage (Dl) d'une coupe d'affichage (50') au moyen de l'unité d'affichage (15) ; 25 ajustement (D4) de la taille de l'image 2D finale (70) à la taille de la coupe d'affichage ; affichage (D5) sur demande de l'utilisateur, au moyen de l'unité d'affichage (15), de l'image 2D finale (70) ajustée. 30 16. Procédé selon la revendication 15, comprenant à l'issu de l'affichage de l'image 2D finale ajustée (70) une étape de retour (D6') à un affichage de la coupe d'affichage (50').5 17. Système (10) d'imagerie comprenant un détecteur (11) de rayons X, une source (12) émettrice de rayons X, une unité de traitement (13) et des moyens de mise en oeuvre d'un procédé de traitement selon l'une des revendications 1 à 16. 18. Produit (CPIO) de programme d'ordinateur comprenant des instructions de code enregistrées sur un support utilisable dans un ordinateur, caractérisé en ce qu'il comprend des instructions pour la mise en oeuvre du procédé de traitement selon l'une des revendications 1 à 12. 10