Brevet FR3160098A1 — publié le 20250919

Ce procédé de reconstruction de plan de tomosynthèse d’un organe cible (14) est mis en œuvre par ordinateur et comprend le chargement (302) d’images de projection acquises par un détecteur de rayons X (10) d’un appareil de tomosynthèse (2), et des itérations d’étapes de calcul (303, 304, 305) d’images reconstruites à partir des images de projection acquises et d’un mouvement tridimensionnel présélectionné ou actualisé de l’organe cible (14), d’images de reprojection, d’images de résidu de projection, et d’actualisation (306) du mouvement tridimensionnel en fonction des images de résidu de projection. Chaque itération est mise en œuvre avec le mouvement tridimensionnel actualisé par l’itération précédente. Figure pour l’abrégé : Fig 3 Procédé de reconstruction de plan de tomosynthèse d’un organe cible lors de mesures de tomosynthèse, et appareil de tomosynthèse associé L’invention a pour domaine technique la radiographie et, plus précisément, la tomosynthèse qui permet de reconstruire le volume d’un organe cible, par exemple le sein d’une patiente, à partir d’un nombre restreint de projections radiographiques de l’organe cible. En particulier, la présente invention concerne un procédé de reconstruction de plan de tomosynthèse de l’organe cible et un appareil de tomosynthèse. Techniques antérieures La tomosynthèse mammaire est une technologie d'imagerie tridimensionnelle consistant à acquérir des images bidimensionnelles, dites de projection, d'un sein comprimé d’une patiente. Chaque image de projection est acquise depuis un détecteur de rayons X fixe pour un angle de mesure particulier au cours du mouvement d’une source de rayons X. Le mouvement de la source de rayons X est un mouvement angulaire restreint qui s’accompagne d’un nombre limité d’acquisitions, par exemple moins de dix images de projection. La tomosynthèse consiste à reconstruire des vues en coupe du sein comprimé dans des plans parallèles au détecteur de rayons X à partir des images de projection, permettant ainsi de réduire l’ambiguïté due à la superposition des tissus des images de mammographie conventionnelle plein champ. Le mouvement de la patiente au cours de la tomosynthèse, par exemple sa respiration, peut créer des artefacts, par exemple des zones floues, des réplications d’objets ou des bavures verticales, dans les vues en coupe reconstruites. Il est donc nécessaire de corriger le mouvement pour détecter efficacement d’éventuelles caractéristiques radiologiques, telles que des masses ou microcalcifications. Le nombre limité d’acquisitions et le mouvement angulaire restreint de la source de rayons X rendent plus délicate la correction du mouvement. La tomosynthèse permet de recueillir une quantité de données très inférieure à une méthode de tomographie, la correction du mouvement de l’organe cible en tomosynthèse étant donc plus complexe qu’en tomographie. L’opération de reconstruction en tomosynthèse n’est pas un problème complet du fait du nombre limité d’acquisitions et donc de données qui caractérise cette opération. Le brevet US 2021/0204899 A1 présente une méthode permettant de quantifier le mouvement du sein comprimé pendant le mouvement de la source de rayons X en comparant un emplacement réel d’un détail d’intérêt dans des images de projection avec un emplacement prévu de l’objet, mais ne donne pas entière satisfaction pour la correction de ce mouvement. La présente invention a donc pour objectif de pallier tout ou partie des inconvénients précités et, en particulier, d’estimer le mouvement tridimensionnel d’un organe cible lors de mesures de tomosynthèse, de déterminer précisément ce mouvement, de le corriger lors de la reconstruction d’un volume de l’organe cible, et ainsi de faciliter la détection et la reconnaissance de détails d’intérêt dans le volume reconstruit de l’organe cible, et de réduire ou au moins de ne pas augmenter l’exposition aux rayons X de l’organe cible. La présente invention a pour objet un procédé d’estimation du mouvement d’un organe cible, en particulier un sein, lors de mesures de tomosynthèse de l’organe cible acquises pour des positions de mesure d’une source de rayons X d’un appareil de tomosynthèse muni d’un détecteur de rayons X, le procédé étant mis en œuvre par ordinateur et comprenant le chargement d’images de projection acquises par le détecteur de rayons X, le procédé comprenant en outre des itérations des étapes : de calcul d’une pluralité d’images reconstruites à partir des images de projection acquises et d’un mouvement tridimensionnel présélectionné ou actualisé de l’organe cible, la pluralité d’images reconstruites formant un volume reconstruit de l’organe cible ; de calcul d’une image de reprojection pour chaque position de mesure de la source de rayons X, chaque image de reprojection étant calculée à partir du volume reconstruit et du mouvement tridimensionnel présélectionné ou actualisé ; de calcul d’une image de résidu de projection pour chaque position de mesure de la source de rayon X, chaque image de résidu de projection étant calculée comme la différence entre l’image de projection acquise considérée et l’image de reprojection correspondante ; et d’actualisation du mouvement tridimensionnel en fonction des images de résidu de projection ; Chaque itération desdites étapes de calcul est mise en œuvre avec le mouvement tridimensionnel actualisé par l’itération précédente. Chaque itération effectuée permet d’augmenter la précision d’estimation du mouvement tridimensionnel. On utilise à chaque itération un volume reconstruit corrigé du mouvement tridimensionnel de l’organe cible pour estimer plus précisément ce mouvement tridimensionnel. Avantageusement, les itérations sont effectuées jusqu’à atteindre un critère de convergence de l’étape d’actualisation du mouvement tridimensionnel. Chaque image de reprojection est calculée comme une image de projection du volume reconstruit depuis une position de mesure. Chaque image de résidu de projection est calculée comme la différence entre les images de projection acquise et de reprojection qui correspondent à une même position de mesure. Le mouvement tridimensionnel présélectionné est librement sélectionné à l’initialisation du procédé. Selon un premier mode de réalisation, le mouvement tridimensionnel présélectionné est nul à l’initialisation du procédé. Le premier mode de réalisation correspond à une mise en œuvre du procédé sans connaissance a priori du mouvement de l’organe cible. Selon un deuxième mode de réalisation le mouvement tridimensionnel présélectionnée est estimé à l’initialisation du procédé par une caméra visible ou par un capteur de mouvement. Le deuxième mode de réalisation permet une mise en œuvre plus rapide du procédé en ayant une connaissance a priori du mouvement de l’organe cible. Avantageusement, le mouvement tridimensionnel présélectionné est décomposé sur une base cinématique présélectionnée définissant les mouvements de l’organe cible autorisés pour les étapes de calculs. La base cinématique présélectionnée est librement sélectionnée à l’initialisation du procédé. La base cinématique définit un espace de solutions pour l’estimation du mouvement tridimensionnel effectué par l’organe cible au cours des mesures de tomosynthèse. Alternativement, la base cinématique pourrait être actualisée lors des itérations du procédé. De préférence, la base cinématique est une base cinématique représentative des mouvements possibles de l’organe cible lors des mesures de tomosynthèse. Optionnellement, le procédé comprend la sélection par un opérateur d’une base cinématique adaptée aux mesures de tomosynthèse de l’organe cible effectuées. Selon un exemple de mise en œuvre, l’étape d’actualisation du mouvement tridimensionnel comprend, pour chaque image de résidu de projection, le calcul d’un champ de déplacement de l’organe cible permettant de réduire la norme de l’image de résidu de projection considérée, en particulier en appliquant une procédure de minimisation notamment de Gauss-Newton. Avantageusement, le chargement des images de projection acquises comprend une sous-étape de sélection d’une zone restreinte de l’organe cible, chaque itération desdites étapes de calcul étant mise en œuvre de sorte que chaque image reconstruite, chaque image de reprojection et chaque image de résidu de projection correspond à la zone restreinte de l’organe cible. En appliquant le procédé à la zone restreinte, on réduit le temps de calcul, on simplifie la cinématique de l’organe cible et on estime localement son mouvement tridimensionnel. Avantageusement, la zone restreinte est sélectionnée de manière à comprendre un détail d’intérêt. Ladite zone restreinte peut être choisie comme comportant un détail d’intérêt soit pour son interprétation radiologique clinique, soit car montrant un artefact d’acquisition ou de reconstruction, soit encore car indiquant la présence d’un mouvement nuisant à l’interprétation. Optionnellement, la zone restreinte est sélectionnée dans au moins une image de projection acquise, la zone restreinte étant sélectionnée de manière à comprendre un détail d’intérêt. Optionnellement, la zone restreinte est sélectionnée dans au moins une image reconstruite, en particulier calculée à partir des images de projection acquises et d’un mouvement tridimensionnel présélectionné notamment nul, la zone restreinte étant sélectionnée de manière à comprendre un détail d’intérêt. Optionnellement, la zone restreinte est sélectionnée dans au moins une image de résidu de projection, en particulier calculée à partir des images de projection acquises et d’un mouvement tridimensionnel présélectionné notamment nul, la zone restreinte étant sélectionnée de manière à comprendre un détail d’intérêt. Avantageusement, le détail d’intérêt est une caractéristique radiologique de l’organe cible. Optionnellement, la zone restreinte peut être sélectionnée par un opérateur au moyen d’une interface graphique. L’opérateur peut effectuer la sélection de la zone restreinte en fonction d’une observation visuelle de ladite image de projection acquise, de ladite image reconstruite ou de ladite image de résidu de projection. Alternativement, la sélection de la zone restreinte peut être effectuée automatiquement par une méthode d’apprentissage artificiel, notamment profond, appliquée à ladite image de projection acquise, de ladite image reconstruite ou de ladite image de résidu de projection. La méthode de sélection de la zone restreinte dans ladite image est susceptible de diverses variantes et modifications qui apparaîtront à l’homme de l’art. Optionnellement le procédé est mis en œuvre pour une pluralité de zones restreintes de l’organe cible et comprend en outre une étape d’assemblage d’images reconstruites restreintes pour former des images reconstruites étendues de l’organe cible. Les images reconstruites étendues de l’organe cible peuvent représenter l’ensemble du volume de l’organe cible acquis lors des mesures de tomosynthèse, ou une partie de celui-ci. Avantageusement, le chargement des images de projection acquises comprend une sous-étape de modification de la résolution des images de projection acquises utilisées pour mettre en œuvre le procédé, le procédé comprenant des répétitions du chargement et desdites itérations de sorte que chaque répétition est effectuée avec une résolution de chaque image reconstruite, de chaque image de reprojection et de chaque image de résidu de projection strictement supérieure à la résolution utilisée pour la répétition précédente, le mouvement tridimensionnel présélectionné utilisé pour chaque répétition étant égal au mouvement tridimensionnel actualisé par la répétition précédente. En effectuant les répétitions pour plusieurs résolutions, on augmente la précision d’estimation de la cinématique de l’organe cible. En commençant par une répétition effectuée pour une résolution grossière, on actualise rapidement et grossièrement la valeur du mouvement tridimensionnel. En effectuant ensuite des répétitions pour des résolutions de plus en plus fines, on actualise de manière incrémentale le mouvement tridimensionnel pour augmenter la précision de son estimation. Autrement dit, le procédé peut comprendre des répétitions du chargement et desdites itérations, lesdites répétitions étant effectuées en chargeant des images de projection acquises dont la résolution augmente strictement d’une répétition à la suivante, le mouvement tridimensionnel actualisé lors d’une répétition considérée étant utilisé pour initialiser la répétition suivante, les répétitions étant effectuées jusqu’à la résolution de mesure. De préférence, le procédé comprend en outre une étape d’affichage d’au moins une image reconstruite, en particulier à la résolution de mesure du détecteur de rayons X, calculée à partir des images de projection acquises et du mouvement tridimensionnel actualisé par l’ensemble des itérations du procédé. Plus préférentiellement, le procédé comprend en outre une étape d’affichage du volume reconstruit de l’organe cible, en particulier à la résolution de mesure du détecteur de rayons X, calculé à partir des images de projection acquises et du mouvement tridimensionnel actualisé à l’issue du procédé. Avantageusement, les mesures de tomosynthèse sont effectuées de sorte qu’un mouvement de la source de rayons X a une amplitude angulaire par rapport à l’organe cible inférieure à 90°, en particulier inférieure à 60°. La présente invention a également pour objet un appareil de tomosynthèse d’un organe cible, en particulier d’un sein, comprenant un détecteur de rayons X et une source de rayons X mobile entre des positions de mesure, l’appareil comprenant en outre un processeur apte à charger des images de projection acquises par le détecteur de rayons X aux positions de mesure et étant muni : d’un module de calcul d’une pluralité d’images reconstruites à partir des images de projection acquises et d’un mouvement tridimensionnel présélectionné ou actualisé de l’organe cible, la pluralité d’images reconstruites formant un volume reconstruit de l’organe cible ; d’un module de calcul d’une image de reprojection pour chaque position de mesure de la source de rayons X, chaque image de reprojection étant calculée à partir du volume reconstruit et du mouvement tridimensionnel présélectionné ou actualisé ; d’un module de calcul d’une image de résidu de projection pour chaque position de mesure de la source de rayon X, chaque image de résidu de projection étant calculée comme la différence entre l’image de projection acquise considérée et l’image de reprojection correspondante ; et d’un module d’actualisation du mouvement tridimensionnel en fonction des images de résidu de projection ; Le processeur est configuré pour mettre en œuvre des itérations desdits modules de sorte que chaque itération desdits modules de calcul est effectuée avec le mouvement tridimensionnel actualisé par le module d’actualisation du mouvement tridimensionnel lors de l’itération précédente. Le module d’actualisation du mouvement tridimensionnel peut être configuré pour effectuer le calcul, pour chaque image de résidu de projection, d’un champ de déplacement de l’organe cible permettant de réduire la norme de l’image de résidu de projection considérée, en particulier en appliquant une procédure de minimisation notamment de Gauss-Newton. Le processeur peut être apte à effectuer une sélection d’une zone restreinte de l’organe cible et à mettre en œuvre chaque itération desdits modules de calcul de sorte que chaque image reconstruite, chaque image de reprojection et chaque image de résidu de projection correspond à la zone restreinte de l’organe cible. Le processeur peut être apte à effectuer plusieurs sélections de zones restreintes de l’organe cible et à effectuer un assemblage d’images reconstruites restreintes pour former des images reconstruites étendues de l’organe cible. Le processeur peut être apte à effectuer la modification de la résolution des images de projection acquises utilisées pour mettre en œuvre lesdits modules. Le processeur peut être configuré pour mettre en œuvre des répétitions desdits modules de sorte que chaque répétition est effectuée avec une résolution de chaque image reconstruite, de chaque image de reprojection et de chaque image de résidu de projection strictement supérieure à la résolution utilisée pour la répétition précédente, le mouvement tridimensionnel présélectionné utilisé pour chaque répétition étant égal au mouvement tridimensionnel actualisé à la résolution précédente. En particulier, la résolution utilisée peut être augmentée jusqu’à la résolution de mesure. L’appareil de tomosynthèse peut comprendre un écran. Le processeur peut être apte à effectuer l’affichage sur l’écran d’au moins une image reconstruite, en particulier à une résolution de mesure du détecteur de rayons X, calculée à partir des images de projection acquises et du mouvement tridimensionnel actualisé par l’ensemble des itérations. Le processeur peut comprendre un module de sélection du mouvement tridimensionnel présélectionné de l’organe cible. L’appareil de tomosynthèse peut en outre incorporer tous les moyens matériels, modules et logiciels pour la mise en œuvre du procédé précédemment défini. D’autres buts, caractéristiques et avantages de l’invention apparaîtront à la lecture de la description suivante, donnée uniquement à titre d’exemple non limitatif et faite en référence aux dessins annexés sur lesquels : FIG. 1 est une vue schématique d’un appareil de tomosynthèse mammaire selon un exemple de réalisation de l’invention ; FIG. 2 est une vue schématique d’une mesure de tomosynthèse mammaire selon un exemple de réalisation de l’invention ; FIG. 3 , FIG. 4 , FIG. 5 , FIG. 6 et FIG. 7 illustrent schématiquement des étapes d’un procédé de reconstruction de plan de tomosynthèse d’un organe cible selon des exemples de réalisation de l’invention. Description détaillée La FIG. 1 représente schématiquement un appareil de tomosynthèse 2 comprenant un système d’imagerie 4 et un système de traitement 6. Le système d’imagerie 4 est apte à acquérir des mesures de tomosynthèse et est muni d’une source de rayons X 8 et d’un détecteur de rayons X 10. La source de rayons X 8 est mobile, ici sur un arc de cercle 12. Alternativement, la source de rayons X 8 pourrait être mobile sur un segment. La source de rayons X est apte à émettre un faisceau de rayons X collimaté vers un organe cible 14 depuis un ensemble de positions de mesure 16, comme illustré sur la FIG. 2 . Ces positions de mesure 16 sont réparties sur le trajet de la source 12 qui s’étend selon une amplitude angulaire limitée. Le faisceau de rayons X irradie l’organe cible 14, par exemple le sein d’une patiente comprimé par une pelote de compression 18. Le détecteur de rayons X 10 est agencé du côté opposé à la source de rayons X 8 par rapport à l’organe cible 14 et reçoit le rayonnement ayant traversé l’organe cible 14. Le détecteur de rayons X 10 comprend une grille de pixels, par exemple une grille de N lignes et M colonnes de pixels, chaque pixel étant apte à mesurer l’intensité de radiation X reçue pour chaque position de mesure 16 de la source de rayons X 8. Une mesure de tomosynthèse permet d’acquérir, pour chaque position de mesure 16, une image dite de projection ayant une définition de mesure égale au nombre de pixels du détecteur de rayons X 10, par exemple une définition de M*N pixels. La taille physique des pixels du détecteur de rayons X 10 définit la résolution de mesure. De préférence, chaque image de projection acquise subit un traitement préliminaire connu de l’homme de l’art comprenant, par exemple, une correction par une image de blanc de l’intensité reçue par le détecteur de rayons X 10. Le système de traitement 6 comprend un processeur 20 pour traiter les mesures de tomosynthèse et un écran 22 pour interagir avec un opérateur. Le processeur 20 comprend une unité de traitement et une mémoire. La mémoire est apte à stocker des instructions qui, lorsqu’elles sont exécutées par l’unité de traitement, entraînent la mise en œuvre d’un procédé de reconstruction de plan de tomosynthèse de l’organe cible 14, ce procédé permettant en particulier d’estimer le mouvement tridimensionnel de l’organe cible 14. Ce procédé est mis en œuvre par ordinateur à partir des mesures de tomosynthèse et permet d’estimer quel a été le mouvement tridimensionnel de l’organe cible 14 lors de ces mesures. Selon un premier exemple de réalisation illustré à la FIG. 3 , le procédé comprend des étapes 300. On commence, lors d’une étape 301, par présélectionner un mouvement tridimensionnel de l’organe cible 14. Par exemple, on initialise le mouvement tridimensionnel avec une valeur nulle. Cette valeur nulle peut être une valeur par défaut enregistrée par la mémoire du processeur 20. Alternativement, l’appareil de tomosynthèse pourrait comprendre une caméra visible, ou un capteur de mouvement, apte à mesurer un mouvement de l’organe cible. Dans ce cas, on sélectionne le mouvement tridimensionnel de l’organe cible comme le mouvement mesuré par la caméra visible ou par le capteur de mouvement. Ce mouvement tridimensionnel présélectionné sera ensuite actualisé par les étapes du procédé pour estimer la cinématique de l’organe cible 14 au cours du déplacement de la source de rayons X 8. Lors de l’étape 301, on sélectionne également une base cinématique sur laquelle on décompose le mouvement tridimensionnel présélectionné ou actualisé. Cette base cinématique définit les mouvements tridimensionnels de l’organe cible autorisés lors d’étapes de calculs numériques du procédé. La recherche de la cinématique de l’organe cible est donc effectuée dans la base cinématique présélectionnée. De préférence, cette base cinématique permet de décrire les mouvements tridimensionnels de l’organe cible 14 lors de la tomosynthèse. Par exemple, l’opérateur sélectionne une base cinématique apte à décrire des translations dans le plan du détecteur de rayons X 10 et une rotation autour de la direction normale au détecteur de rayons X 10. Le mouvement tridimensionnel et la base cinématique peuvent être sélectionnés via l’écran 22 par l’opérateur, notamment dans des listes enregistrées sur la mémoire du processeur 20. Par exemple, l’opérateur sélectionne le mouvement tridimensionnel en fonction d’une observation visuelle d’images de projection à la résolution de mesure. Si l’opérateur observe, dans une image de projection, une trainée selon une direction particulière, celui-ci peut, par exemple, sélectionner le mouvement tridimensionnel comme comprenant au moins un mouvement de translation selon la direction particulière. La base cinématique est indépendante des données des mesures de tomosynthèse. Ainsi, l’opérateur peut modifier la base cinématique considérée jusqu’à obtenir une meilleure exécution du procédé. Puis, lors d’une étape 302, on charge des images de projection de l’organe cible 14 acquises par le détecteur de rayons X 10. On calcule ensuite, lors d’une étape 303, des images reconstruites de l’organe cible 14 à partir des images de projection acquises et du mouvement tridimensionnel présélectionné ou actualisé, par exemple par une méthode de reconstruction algébrique. Ces images reconstruites constituent un volume numérique reconstruit de l’organe cible 14 que l’opérateur peut consulter tranche par tranche selon la direction allant du détecteur de rayons X 10 à la source 8. Le nombre de tranches du volume numérique est conventionnellement choisi de sorte que la distance entre deux tranches successives correspond à environ 1 millimètre de l’organe cible. Bien entendu, un autre nombre de tranches peut être utilisé. On calcule ensuite, lors d’une étape 304, une image de reprojection pour chaque position de mesure 16 de la source de rayons X 8, une image de projection ayant été acquise à chaque position de mesure 16. On calcule chaque image de reprojection à partir du volume reconstruit et du mouvement tridimensionnel présélectionné ou actualisé associé aux positions de mesure 16 de la source de rayons X 8. Chaque image de reprojection est calculée en appliquant une opération de projection au volume numérique, l’opération de projection décrivant le cumul de l’atténuation des rayons X 8 lors de l’acquisition de données des mesures de tomosynthèse depuis la position de mesure 16 considérée. Puis, lors d’une étape 305, on calcule une image de résidu de projection pour chaque position de mesure 16 de la source de rayons X 8. Chaque image de résidu de projection est calculée comme la différence entre les images de reprojection et de projection acquise qui correspondent à une même position de mesure 16 de la source de rayons X 8. Enfin, lors d’une étape 306, on calcule, à partir de chaque image de résidu de projection, un mouvement tridimensionnel de l’organe cible 14 associé aux positions de mesure 16 de la source de rayons X 8. Par exemple, pour chaque image de résidu de projection, on calcule un champ de déplacement de l’organe cible 14 pour réduire la norme de l’image de résidu de projection considérée. On décompose le champ de déplacement sur la base cinématique présélectionnée puis on résout un problème des moindres carrés en appliquant une procédure de minimisation à la norme de l’image du résidu de la projection considérée, par exemple une procédure de minimisation de Gauss-Newton. Lors de l’étape 306, on actualise également le mouvement tridimensionnel avec la valeur calculée lors de cette étape. On effectue des itérations des étapes de calcul 303, 304, 305 et 306 en remplaçant pour chaque itération de ces étapes de calcul 303, 304, 305 le mouvement tridimensionnel concerné par le mouvement tridimensionnel actualisé lors de l’itération précédente de l’étape d’actualisation 306. Ainsi, on calcule, à chaque itération, un volume reconstruit de l’organe cible 14 à partir des images de projection acquises et du mouvement tridimensionnel actualisé lors de l’étape d’actualisation 306 de l’itération précédente, une image de reprojection de ce volume reconstruit pour chaque position de mesure 16, et une image de résidu de projection pour chaque position de mesure 16. A partir de ces résidus de projection, on actualise le mouvement tridimensionnel de l’organe cible 14. On effectue ces itérations jusqu’à atteindre un critère de convergence fondé, par exemple, sur la norme de l’écart relatif entre deux mouvements de l’organe cible 14 calculés pour deux itérations successives, notamment lorsque cet écart devient inférieur à 10%. Lorsque le critère de convergence est atteint, la mémoire du processeur 20 effectue l’enregistrement, lors d’une étape 307, du mouvement tridimensionnel de l’organe cible 14 calculé au cours des itérations des étapes de calculs 303, 304, 305 et 306. Selon un deuxième exemple de réalisation illustré à la FIG. 4 , on effectue les étapes 300 illustrées à la FIG. 3 de manière répétée en modifiant à chaque itération la résolution des images de projection acquises. Plus précisément, on réduit la résolution, puis on l’augmente à chaque itération jusqu’à atteindre la résolution de mesure du détecteur. On initialise chaque itération des étapes 300 par le mouvement tridimensionnel de l’organe cible 14 enregistré lors de la répétition précédente. Ainsi, on effectue des répétitions des étapes 300 du procédé. Le terme itération est employé pour décrire chaque itération des étapes de calculs et d’actualisation 303, 304, 305 et 306 effectuée à une résolution donnée. Le terme répétition est employé pour décrire chaque répétition des étapes 300 du procédé effectuée avec des images de projection acquises de résolution modifiée, chaque répétition comprenant donc plusieurs itérations. Lors de chaque répétition de l’étape de chargement 302, on effectue une sous-étape de modification de la résolution des images de projections acquises utilisées pour la répétition concernée. La sous-étape de modification de la résolution de chaque répétition est effectuée de manière à augmenter strictement la résolution des images de projections acquises par rapport à la répétition précédente. Chaque image de projection acquise de résolution modifiée est, par exemple, obtenue en appliquant un filtre de Gauss à l’image de projection acquise à la résolution de mesure qui correspond à la position de mesure 16 considérée, puis en effectuant un sous-échantillonnage de l’image filtrée. Par exemple, on applique un filtre de Gauss de variance L à l’image acquise à la résolution de mesure, puis on remplace une zone de L*L pixels par un unique pixel ayant pour intensité la moyenne de l’intensité des L*L pixels de la zone considérée. On obtient, lors des répétitions des étapes 300 du procédé, une approximation du mouvement de l’organe cible 14 comme un développement de Taylor à un nombre réduit d’ordres. Ainsi, lors de chaque répétition de l’étape de calcul 303, on calcule chaque image reconstruite à la résolution modifiée lors de la sous-étape de modification de la résolution de la répétition considérée. De même, lors de chaque répétition de l’étape de calcul 304, on calcule chaque image de reprojection à la résolution modifiée lors de la sous-étape de modification de la résolution de la répétition considérée. Enfin, lors de chaque répétition de l’étape de calcul 305, on calcule chaque image de résidu de projection comme la différence entre des images de projection acquise et de reprojection chacune à la résolution modifiée lors de la sous-étape de modification de la résolution de la répétition considérée. Le mouvement tridimensionnel présélectionné utilisé pour initialiser les étapes 300 de chaque répétition est égal au mouvement tridimensionnel actualisé par la répétition précédente, c’est-à-dire le mouvement tridimensionnel enregistré lors de l’étape 307 de la répétition précédente. Dans l’exemple représenté, on effectue des répétitions pour des résolutions égales au dixième, au cinquième, à la moitié de la résolution de mesure, et à la résolution de mesure. Avantageusement, la résolution utilisée pour chaque répétition des étapes 300 peut être sélectionnée par l’opérateur via l’écran 22. L’opérateur sélectionne notamment une diminution de la résolution d’au moins un facteur égal à un mouvement estimé de l’organe cible 14 divisé par la taille d’un pixel. Par exemple, l’opérateur estime un mouvement de l’organe cible 14 de 1 millimètre. Si chaque pixel du détecteur de rayon X mesure 100 microns, l’opérateur sélectionne une résolution égale au dixième de la résolution de mesure de sorte que le mouvement de 1 millimètre soit détectable. Ainsi, le mouvement tridimensionnel de l’organe cible 14 sera estimé précisément. Optionnellement, à la suite de la répétition effectuée à la résolution de mesure, on effectue une étape de calcul du volume reconstruit à la résolution de mesure à partir des images de projection acquises à la résolution de mesure et du mouvement tridimensionnel enregistré lors de l’étape 307 qui est effectuée pendant la répétition effectuée à la résolution de mesure. On effectue ensuite une étape d’affichage de ce volume reconstruit à la résolution de mesure et corrigé du mouvement tridimensionnel actualisés par les étapes du procédé, c’est-à-dire par les itérations et les répétitions. Selon un troisième exemple de réalisation illustré aux figures 5 à 7, on effectue les étapes 300 illustrées à la FIG. 3 pour une zone restreinte de l’organe cible 14 et donc pour une partie des données des mesures de tomosynthèse. Plus précisément, lors de l’étape de chargement 302, on effectue une sous-étape de sélection d’une zone restreinte comprenant un détail d’intérêt de l’organe cible 14, par exemple une zone de strictement moins de M*N pixels. Puis on effectue chaque étape de calcul 303, 304, 305 pour la zone restreinte sélectionnée. Ainsi, lors de chaque étape de calcul 303, on calcule des images reconstruites de la zone restreinte pour former un volume reconstruit de la zone restreinte. De même, lors de chaque étape de calcul 304, on calcule chaque image de reprojection à partir du volume reconstruit de la zone restreinte. Enfin, lors de chaque étape de calcul 305, on calcule chaque image de résidu de projection comme la différence entre des images de projection acquise et de reprojection correspondant à la zone restreinte. Dans l’exemple illustré à la FIG. 5 , on sélectionne la zone restreinte de l’organe cible 14 à partir d’une image reconstruite de l’ensemble de l’organe cible 14, en particulier reconstruite pour un mouvement tridimensionnel de l’organe cible 14 nul. Par exemple, l’opérateur repère visuellement une caractéristique radiologique 501 dans l’organe cible 14 et sélectionne via l’écran 22 la zone restreinte sur laquelle effectuer les étapes 300. Dans l’exemple illustré à la FIG. 6 , on sélectionne la zone restreinte de l’organe cible 14 à partir d’une image de résidu de projection de l’ensemble de l’organe cible 14. Par exemple, l’opérateur repère visuellement un résidu de projection d’intérêt 601, l’assimile à une caractéristique radiologique dans l’organe cible 14, et sélectionne via l’écran 22 la zone restreinte sur laquelle effectuer les étapes 300. Dans l’exemple illustré à la FIG. 7 , on décompose l’organe cible 14 en une pluralité de zones restreintes et on met en œuvre les étapes 300 pour chaque zone restreinte. On effectue ensuite, lors d’une étape 700, l’assemblage d’images reconstruites restreintes correspondant à chaque zone restreinte pour former des images reconstruites étendues de l’organe cible 14. Dans chacun des exemples illustrés aux figures 5 à 7, en restreignant la zone sur laquelle on effectue les étapes 300, on simplifie la description du mouvement de l’organe cible. Dans chacun des exemples illustrés aux figures 5 à 7, on peut effectuer les étapes 300 pour chaque zone restreinte à la résolution de mesure ou on peut effectuer des répétitions des étapes 300 pour chaque zone restreinte en augmentant la résolution à chaque répétition et ce jusqu’à la résolution de mesure. Ainsi, le deuxième exemple de réalisation illustré à la FIG. 4 et le troisième exemple de réalisation illustré aux figures 5 à 7 sont compatibles. Procédé de reconstruction de plan de tomosynthèse d’un organe cible (14), en particulier un sein, lors de mesures de tomosynthèse de l’organe cible (14) acquises pour des positions de mesure (16) d’une source de rayons X (8) d’un appareil de tomosynthèse (2) muni d’un détecteur de rayons X (10), le procédé étant mis en œuvre par ordinateur et comprenant le chargement (302) d’images de projection acquises par le détecteur de rayons X (10), caractérisé en ce que le procédé comprend en outre des itérations des étapes : de calcul (303) d’une pluralité d’images reconstruites à partir des images de projection acquises et d’un mouvement tridimensionnel présélectionné ou actualisé de l’organe cible (14), la pluralité d’images reconstruites formant un volume reconstruit de l’organe cible (14) ; de calcul (304) de l’image de reprojection pour chaque position de mesure (16) de la source de rayons X (8), chaque image de reprojection étant calculée à partir du volume reconstruit et du mouvement tridimensionnel présélectionné ou actualisé ; de calcul (305) d’une image de résidu de projection pour chaque position de mesure (16) de la source de rayons X (8), chaque image de résidu de projection étant calculée comme la différence entre l’image de projection acquise considérée et l’image de reprojection correspondante ; et d’actualisation (306) du mouvement tridimensionnel en fonction des images de résidu de projection ; chaque itération desdites étapes de calcul (303, 304, 305) étant mise en œuvre avec le mouvement tridimensionnel actualisé par l’itération précédente (306). Procédé selon la revendication 1, dans lequel le chargement (302) des images de projection acquises comprend une sous-étape de sélection d’une zone restreinte de l’organe cible (14), chaque itération desdites étapes de calcul (303, 304, 305) étant mise en œuvre de sorte que chaque image reconstruite, chaque image de reprojection et chaque image de résidu de projection corresponde à la zone restreinte de l’organe cible (14). Procédé selon la revendication 2, dans lequel la zone restreinte est sélectionnée dans au moins une image de projection acquise, la zone restreinte étant sélectionnée de manière à comprendre un détail d’intérêt. Procédé selon la revendication 2, dans lequel la zone restreinte est sélectionnée dans au moins une image reconstruite, en particulier calculée à partir des images de projection acquises et d’un mouvement tridimensionnel présélectionné notamment nul, la zone restreinte étant sélectionnée de manière à comprendre un détail d’intérêt. Procédé selon la revendication 2, dans lequel la zone restreinte est sélectionnée dans au moins une image de résidu de projection, en particulier calculée à partir des images de projection acquises et d’un mouvement tridimensionnel présélectionné notamment nul, la zone restreinte étant sélectionnée de manière à comprendre un détail d’intérêt. Procédé selon l’une quelconque des revendications 2 à 5, mis en œuvre pour une pluralité de zones restreintes de l’organe cible (14) et comprenant en outre une étape d’assemblage d’images reconstruites restreintes (700) pour former des images reconstruites étendues de l’organe cible (14). Procédé selon l’une quelconque des revendications 1 à 6, dans lequel le chargement (302) des images de projection acquises comprend une sous-étape de modification de la résolution des images de projection acquises utilisées pour mettre en œuvre le procédé, le procédé comprenant des répétitions du chargement (302) et desdites itérations de sorte que chaque répétition est effectuée avec une résolution de chaque image reconstruite, de chaque image de reprojection et de chaque image de résidu de projection strictement supérieure à la résolution utilisée pour la répétition précédente, le mouvement tridimensionnel présélectionné utilisé pour chaque répétition étant égal au mouvement tridimensionnel actualisé par la répétition précédente. Procédé selon l’une quelconque des revendications 1 à 7, comprenant en outre une étape d’affichage d’au moins une image reconstruite, en particulier à une résolution de mesure du détecteur de rayons X (10), calculée à partir des images de projection acquises et du mouvement tridimensionnel actualisé par l’ensemble des itérations du procédé. Procédé selon l’une quelconque des revendications 1 à 8, dans lequel l’étape d’actualisation (306) du mouvement tridimensionnel comprend, pour chaque image de résidu de projection, le calcul d’un champ de déplacement de l’organe cible (14) permettant de réduire la norme de l’image de résidu de projection considérée, en particulier en appliquant une procédure de minimisation notamment de Gauss-Newton. Procédé selon l’une quelconque des revendications 1 à 9, dans lequel le mouvement tridimensionnel présélectionné est décomposé sur une base cinématique présélectionnée définissant les mouvements de l’organe cible autorisés pour les étapes de calculs. Procédé selon l’une quelconque des revendications 1 à 10, dans lequel le mouvement tridimensionnel présélectionné est nul à l’initialisation du procédé. Procédé selon l’une quelconque des revendications 1 à 10, dans lequel le mouvement tridimensionnel présélectionnée est estimé à l’initialisation du procédé par une caméra visible ou par un capteur de mouvement. Procédé selon l’une quelconque des revendications 1 à 12, dans lequel les mesures de tomosynthèse sont effectuées de sorte qu’un mouvement de la source de rayons X (8) a une amplitude angulaire par rapport à l’organe cible inférieure à 90°, en particulier inférieure à 60°. Appareil de tomosynthèse (2) d’un organe cible (14), en particulier d’un sein, comprenant un détecteur de rayons X (8) et une source de rayons X (8) mobile entre des positions de mesure (16), l’appareil (2) comprenant en outre un processeur (20) apte à charger des images de projection acquises par le détecteur de rayons X (8) aux positions de mesure (16) et étant muni : d’un module de calcul d’une pluralité d’images reconstruites à partir des images de projection acquises et d’un mouvement tridimensionnel présélectionné ou actualisé de l’organe cible (14), la pluralité d’images reconstruites formant un volume reconstruit de l’organe cible (14) ; d’un module de calcul d’une image de reprojection pour chaque position de mesure (16) de la source de rayons X (8), chaque image de reprojection étant calculée à partir du volume reconstruit et du mouvement tridimensionnel présélectionné ou actualisé ; d’un module de calcul d’une image de résidu de projection pour chaque position de mesure (16) de la source de rayon X (8), chaque image de résidu de projection étant calculée comme la différence entre l’image de projection acquise considérée et l’image de reprojection correspondante ; et d’un module d’actualisation du mouvement tridimensionnel en fonction des images de résidu de projection ; le processeur étant configuré pour mettre en œuvre des itérations desdits modules de sorte que chaque itération desdits modules de calcul est effectuée avec le mouvement tridimensionnel actualisé par le module d’actualisation du mouvement tridimensionnel lors de l’itération précédente.