Brevet FR3071715A1 — publié le 20190405

PROCEDE D'IMAGERIE RADIOGRAPHIQUE, DISPOSITIF DE TRAITEMENT D'IMAGE RADIOGRAPHIQUE ET DISPOSITIF D'IMAGERIE RADIOGRAPHIQUE. L'invention concerne un procédé d'imagerie radiographique d'une structure interne en trois dimensions faisant partie d'un objet situé à l'intérieur d'un champ d'observation, un dispositif de traitement d'image radiographique et un dispositif d'imagerie radiographique. Notamment, le procédé d'imagerie selon l'invention est adapté pour localiser un point caractéristique d'une structure interne sur une image radiographique. Dans un mode de réalisation, le procédé d'imagerie selon l'invention est également adapté pour calculer un modèle à trois dimensions d'un objet à imager dans un champ d'observation. Des procédés d'imagerie radiographique pour la reconstruction tridimensionnelle sont connus. A titre d'exemple, la demande de brevet français FR-A-2 810 769 décrit un procédé dans lequel on calcule la forme à trois dimensions d'un modèle représentant l'objet à partir d'un modèle générique connu a priori de l'objet. L'utilisateur doit repérer manuellement sur chacune des deux images radiographiées des repères de contrôle appartenant au dit objet, qui peuvent être des points, des segments, des droites, des arcs, des contours, des arêtes ou autres, ce qui demande à la fois de sa part des connaissances précises en anatomie pour localiser précisément les repères de contrôle, du temps et une grande précision. Après ce repérage, une position géométrique de chaque repère est déterminée dans un référentiel à trois dimensions, puis utilisée pour calculer la forme du modèle à partir du modèle générique connu. Cette importante intervention humaine au niveau du repérage manuel sur les radiographies limite en particulier la reproductibilité de la reconstruction lorsque celle-ci est réalisée par un personnel non hautement spécialisé. Il existe des procédés d'imagerie cherchant à limiter l'intervention de spécialistes. C'est le cas du procédé décrit dans le document FR-A 2 856 170. Le procédé décrit utilise un volume de confinement de l'objet, estimé à partir d'un motif géométrique visible sur deux images, et de la position de la source. Il utilise un modèle géométrique comprenant des informations permettant d'établir, depuis un 2 estimateur de l'objet une caractéristique géométrique pour le modèle représentant l'objet. Le procédé décrit présente l'inconvénient de calculer la forme à trois dimensions du modèle à partir de points obtenus uniquement par 5 traitement des images radiographiques. En cas d'identification erronée des points, le modèle est faussé. La présente invention a notamment pour but de fournir un procédé fiable dans lequel l'intervention de spécialistes est limitée. En particulier, un but de l'invention est de faciliter et améliorer 10 l'identification des structures anatomiques internes sur une radiographie à deux dimensions. Un autre but est de proposer un procédé d'imagerie de reconstruction tridimensionnelle d'un objet fiable et précis. A cet effet, un premier objet de l'invention concerne un procédé 15 d'imagerie radiographique d'une structure interne en trois dimensions faisant partie d'un objet situé à l'intérieur d'un champ d'observation. Selon l'invention, le procédé comprend : a) recevoir des premières données d'image générées dans un référentiel de référence à trois dimensions, représentatives d'au moins 20 une image radiographique en deux dimensions de la structure interne de l'objet, b) recevoir des deuxièmes données d'image générées dans le référentiel de référence, représentatives d'une image en trois dimensions d'une enveloppe externe de l'objet, 25 c) estimer à partir des deuxièmes données d'image les coordonnées dans le référentiel de référence de premiers repères caractéristiques de la structure interne, et, optionnellement, la valeur d'au moins un paramètre caractéristique de la structure interne, d) calculer, dans le référentiel de référence, la forme à trois 30 dimensions d'un modèle représentant ladite structure interne à partir d'un modèle générique connu a priori de ladite structure interne, ce modèle générique comportant des points correspondants aux premiers repères, le modèle calculé étant obtenu par déformation du modèle générique de façon que ledit modèle calculé suive une forme la plus proche possible d'une isométrie du modèle générique tout en maintenant en coïncidence les points du modèle générique déformé avec 3 les premiers repères, et optionnellement en appliquant au modèle calculé le ou les paramètres caractéristiques de la structure interne, e) choisir au moins un deuxième repère caractéristique de la structure interne à localiser sur l'image radiographique, et pour chaque 5 deuxième repère : i) identifier un point appartenant au modèle calculé à l'étape d) correspondant au deuxième repère, (ii) calculer les coordonnées d'une projection de ce point sur l'image radiographique, 10 iii) déterminer à partir des premières données d'image les coordonnées du deuxième repère sur l'image radiographique en le recherchant dans une zone de recherche située autour des coordonnées de la projection. L'objet mentionné ci-dessus peut comprendre notamment le 15 rachis, le bassin, ou encore le genou d'un patient, ou plus généralement être constitué par tout ou partie du corps du patient. Dans ces différents cas, les structures internes à observer peuvent être constituées notamment par les os du patient compris dans le champ d'observation. 20 Le ou les repères mentionnés ci-dessus sont des repères anatomiques qui peuvent être des points, des segments, des droites, des arcs, des contours, des arêtes ou autres, ou plus généralement être constitué par tout point ou élément anatomique caractéristique du squelette du patient. Eventuellement, les premiers et deuxièmes repères 25 peuvent être les mêmes repères. Le ou les paramètres caractéristiques de la structure interne peuvent être des paramètres anatomiques de la structure interne, tel que par exemple des dimensions de repères anatomiques. A titre d'exemple, le tour de taille ou la hauteur séparant le tour de taille au 30 niveau sous axillaire sont des dimensions morphologiques (mesurables sur l'enveloppe externe) permettant d'estimer une dimension anatomique de la structure interne (largeur de l'os iliaque, hauteur des vertèbres). Le procédé selon l'invention permet une intervention réduite, voire 35 nulle, d'un opérateur pour localiser les points anatomiques (deuxièmes repères) sur une radiographie à deux dimensions, facilitant ainsi cette localisation et la rendant plus fiable. Ceci est obtenu grâce à 4 l'estimation de la position géométrique des premiers repères à partir des données de l'image numérique en trois dimensions de l'enveloppe externe, cette estimation permettant d'adapter un modèle générique connu a priori de la structure interne à observer, qui sert ensuite à 5 localiser les deuxièmes repères sur une image radiographiée. Ainsi, la recherche de repères sur les radiographies 2D (en deux dimensions) est initiée par la numérisation 3D (en trois dimensions) de l'enveloppe externe de l'objet. On limite de cette manière les risques de mauvaise identification des repères sur les radiographies. Ces points sont en effet 10 recherchés dans des zones de recherche limitées, définies autour des projections des repères sur un modèle calculé et ne sont pas définies uniquement par traitement des images radiographiques. Autrement dit, la localisation des deuxièmes repères sur les images radiographiques utilise une image numérisée en 3D de l'objet, 15 exploitée pour corriger un modèle générique existant. La projection de ce modèle corrigé permet de déterminer une position probable des deuxièmes repères sur une radiographie en deux dimensions. Un traitement de l'image radiographique autour de chaque position probable permet alors de localiser les deuxièmes repères de manière 20 fiable. Ainsi, ce traitement d'image nécessite moins de ressources car il est effectué dans une zone restreinte de l'image radiographique. On notera que l'étape c) est réalisée avant l'étape (d), l'étape (b) est réalisée avant l'étape (c) et que l'étape (a) peut être réalisée à tout moment avant l'étape (e). 25 Les premières et deuxièmes données d'image sont par exemple des données représentatives des coordonnées des points d'une image en 2D ou en 3D. Il peut s'agir d'un signal électrique (signal d'image). Les deuxièmes données d'image reçues à l'étape (b) peuvent être générées par un système d'imagerie comprenant au moins un capteur 30 optique. Ces deuxièmes données d'image peuvent avantageusement être générées au moyen d'un ou plusieurs capteurs optiques, notamment non irradiants ou ionisants, par exemple un ou plusieurs capteurs optiques d'image, convertissant un signal lumineux en un signal électrique. Le signal lumineux utilisé peut être visible ou non. A titre 35 d'exemple, un système d'imagerie apte à mesurer une profondeur optique est utilisable, par exemple du type décrit dans le document WO 2007/043036. 5 Dans un mode de réalisation, au cours de l'étape (c), la position géométrique estimée de la pluralité de premiers repères peut être déterminée en procédant tout d'abord à la détermination (c)(i) des coordonnées de repères externes positionnés sur l'enveloppe externe 5 dans un référentiel à trois dimensions, puis par l'estimation (c)(ii) des coordonnées des premiers repères dans le même référentiel à partir des repères externes. L'étape (c)(i) peut être obtenue par calcul au moyen d'un modèle mathématique, ou encore par traitement des deuxièmes données 10 d'image, ces repères externes étant marqués sur l'enveloppe externe. On comprendra que ce marquage est réalisé avant la génération des deuxièmes données d'image reçues à l'étape (b). Le marquage, manuel, utilisé en cas de traitement des données d'image peut être un marquage passif. Il peut consister à positionner 15 des repères externes dessinés ou collés sur la peau du patient, présentant par exemple une forme et/une couleur particulière ou présentant une surface réfléchissante pouvant réfléchir un flash de lumière. Les repères externes dessinés ou collés peuvent être repérés de 20 manière manuelle, par exemple désignés sur l'image numérisée 3D par un opérateur, ou de manière automatique par exemple par des logiciels de traitement d'image permettant de différencier des couleurs, des intensités lumineuses,.... Une intervention humaine est nécessaire dans ce cas, mais reste limitée. 25 Le marquage peut également être un marquage actif, tel que le positionnement d'émetteurs sur la peau du patient, de tels émetteurs pouvant éventuellement émettre un signal caractéristique de leur position apte à être détecté par le système d'imagerie permettant de numériser l'enveloppe externe. Ce signal, par exemple un signal 30 lumineux, dans le domaine visible ou non, peut être propre à chaque marqueur et/ou être émis ponctuellement à des intervalles de temps connus ou pendant une durée connue. L'intervention humaine est alors limitée au positionnement des marqueurs sur le patient. L'étape (c)(ii) peut notamment être mise en oeuvre par calcul au 35 moyen d'un modèle mathématique. On comprend qu'aucune intervention humaine n'est alors nécessaire lorsque l'étape (c)(i) utilise aussi un modèle mathématique, améliorant la reproductibilité de la 6 localisation. A titre d'exemple, on peut utiliser un modèle permettant de localiser des repères internes (appartenant à la structure interne) à partir de l'enveloppe externe. Le procédé d'imagerie selon l'invention permet ainsi de localiser 5 de manière fiable et quasiment automatique des repères sur une ou plusieurs images radiographiques. Les premières données d'image reçues à l'étape (a) peuvent être générées par balayage en déplaçant au moins une source radioactive dans une direction de translation non parallèle à une direction de prise 10 de vue. De manière générale, les deuxièmes repères sont de préférence choisis parmi des repères à localiser sur chacune des deux images radiographiques (aussi appelés stéréo correspondants), et éventuellement parmi des repères à localiser sur une seule image 15 radiographique (aussi appelés non stéréo correspondants). Dans un mode de réalisation particulier, au cours de l'étape (a), on reçoit des premières données d'image représentatives de deux images radiographiques enregistrées selon deux directions de prise de vue non parallèles. 20 La réception de premières données d'images représentatives de deux images radiographiques capturées selon deux directions de prise de vue non parallèle peut être utilisée pour la reproduction en trois dimensions de la structure interne d'un objet. A cet effet, au cours de l'étape (e), on choisit les deuxièmes 25 repères parmi des deuxièmes repères (stéréo correspondants) à localiser sur chacune des deux images radiographiques, et, optionnellement parmi des deuxièmes repères (non stéréo correspondants) à localiser sur une seule image radiographique. Le procédé comprend alors en outre : 30 f) calculer les coordonnées des deuxièmes repères dans le référentiel de référence à partir de leurs coordonnées sur au moins une des images radiographiques déterminées à l'étape (e), g) corriger le modèle calculé à l'étape (d) en faisant coïncider des points du modèle calculé correspondant aux deuxièmes repères avec les 35 coordonnées de ceux-ci calculés à l'étape (f). Le modèle calculé à partir d'un modèle générique connu a priori de la structure interne est ainsi établi dans un premier temps pour 7 localiser les deuxièmes repères de l'objet (étape (e) (iii)) sur une ou sur les deux radiographies, puis est simplement corrigé pour obtenir un modèle fiable à trois dimensions de la structure interne. L'utilisation de deuxièmes repères stéréo correspondants, à savoir 5 visibles sur les deux images radiographiées, permet d'améliorer la qualité du modèle corrigé. La qualité du modèle corrigé peut encore être améliorée en utilisant d'autres deuxièmes repères non stéréo correspondants. Le nombre de deuxièmes repères utilisables pour corriger le modèle pourra être choisi en fonction de la nature de la 10 structure interne à imager. A titre d'exemple, pour une vertèbre, l'utilisation d'une dizaine de repères peut suffire. En général, on utilisera de 10 à 100 repères de contrôle. Avantageusement, dans un mode de réalisation, les premières données d'image reçues à l'étape (a) sont générées simultanément, par 15 balayage, en déplaçant en synchronisme, dans une même direction de translation non parallèle aux directions de prise de vues, deux sources radioactives émettant deux faisceaux de rayons ionisants respectivement dans les deux directions de prise de vue. Ceci peut permettre de générer des données d'images radiographiques dans un 20 même référentiel à trois dimensions, par exemple après une calibration appropriée des sources radioactives. Dans un mode de réalisation, les premières données d'image et les deuxièmes données d'image sont générées simultanément. Les premières et deuxièmes données images peuvent alors être générées 25 dans un même référentiel à trois dimensions, par exemple après une calibration appropriée des dispositifs d'enregistrement des images. Ceci permet d'éviter de procéder à un recalage des différentes images dans le même référentiel à trois dimensions et de simplifier la mise en oeuvre du procédé selon l'invention. Ce référentiel peut alors servir de référentiel 30 de référence pour le procédé de l'invention. Dans un autre mode de réalisation, les premières données d'image et les deuxièmes données d'image peuvent être générées dans des référentiels à trois dimensions différents. Autrement dit, chaque image ou chaque type d'image est enregistrée dans un référentiel à trois 35 dimensions qui lui est propre. C'est par exemple le cas lorsque les images sont générées à des moments distincts. Dans ce cas, au moins une des étapes de réception (a) et (b) comprend une étape de 8 transformation des données d'image reçues vers le référentiel de référence. Cette étape de transformation peut par exemple être réalisée au moyen d'une matrice de transformation qui permet de passer d'un 5 référentiel à un autre. Cette matrice de transformation peut être déterminée de sorte que des éléments communs aux premières données d'image et aux deuxièmes données d'image coïncident dans le référentiel de référence choisi. Ces éléments communs peuvent faire partie de l'image tel que, par exemple, le contour des images ou être des 10 éléments de marquage radio-opaques de dimensions connus positionnés sur l'objet à imager avant la génération des données d'image. Le procédé d'imagerie selon l'invention peut en outre comprendre une étape de transmission de données, notamment vers un dispositif 15 d'affichage. Ces données peuvent être: - des premières et/ou des deuxièmes données d'image, permettant de visualiser sur un écran la ou les images radiographiques/ 3D ; - des données représentatives des coordonnées des premiers repères, 20 - des données représentatives des coordonnées de chaque deuxième repère déterminé à l'étape (e)iii), ce qui permet de visualiser sur un écran la position des deuxièmes repères, notamment sur la ou les images radiographiques, et, - le cas échéant des données représentatives du modèle corrigé 25 calculé à l'étape (g), ce qui permet de visualiser une reconstruction en 3D de la structure interne. La visualisation de la reconstruction en 3D pourra être séparée des autres vues. Il est également proposé un dispositif de traitement d'image 30 radiographique comprenant des moyens agencés pour effectuer les étapes du procédé d'imagerie selon l'invention. Notamment, le dispositif de traitement d'image radiographique comprend : - des moyens de réception de données d'image pour recevoir 35 des premières données d'image générées dans un référentiel à trois dimensions représentatives d'au moins une image radiographique en deux dimensions de la structure interne de l'objet, et des deuxièmes 9 données d'image générées dans un référentiel à trois dimensions représentatives d'une image en trois dimensions d'une enveloppe externe de l'objet - des moyens de traitement raccordés aux moyens de 5 réception agencés pour effectuer au moins les étapes c) à e) précédemment décrites du procédé selon l'invention, - optionnellement des moyens de transmission, notamment vers un dispositif d'affichage. Ce dispositif de traitement d'image peut comprendre ou être 10 intégré dans un ou plusieurs processeurs, par exemple des microcontrôleurs, des microprocesseurs ou autres. Les moyens de réception sont par exemple un port d'entrée, une broche d'entrée, ou autre. Les moyens de traitement peuvent par exemple comprendre un coeur de processeur ou CPU (de l'anglais 15 « Central Processing Unit »). Les moyens de transmission peuvent par exemple comprendre un port de sortie, une broche de sortie ou autre. Les moyens de traitement peuvent également être agencés pour effectuer les étapes (f) et (g) précédemment décrites. Le dispositif de traitement d'image peut également comprendre 20 des moyens de mémorisation reliés aux moyens de traitement. Il est également proposé un dispositif d'imagerie radiographique d'une structure interne en trois dimensions faisant partie d'un objet comprenant : - un premier système d'imagerie radiographique agencé pour 25 prendre au moins une image radiographique en deux dimensions de la structure interne de l'objet au moyen d'au moins une source radioactive et pour générer dans un référentiel à trois dimensions des premières données d'image représentatives de l'image radiographique prise, - un deuxième système d'imagerie comprenant au moins un 30 capteur optique, agencé pour prendre une image en trois dimensions d'une enveloppe externe de l'objet et pour générer dans un référentiel à trois dimensions des deuxièmes données d'image représentatives de cette image, - le dispositif de traitement d'image radiographique selon 35 l'invention connecté aux premier et deuxième systèmes d'imagerie pour recevoir les premières et deuxièmes données d'image. 10 Il est en outre proposé un produit programme d'ordinateur comprenant des instructions ou portions de code pour exécuter les étapes du procédé d'imagerie décrit ci-dessus lorsque ces instructions sont exécutées par un dispositif de traitement d'image radiographique, 5 notamment celui de l'invention, comprenant notamment un processeur. Ce programme peut par exemple être stocké sur un support mémoire de type disque dur, téléchargé, ou autre. D'autres caractéristiques et avantages de l'invention apparaîtront au cours de la description suivante d'une de ses formes de réalisation, 10 donnée à titre d'exemple non limitatif, en référence aux dessins annexés, dans lesquels : - la figure 1 est une vue schématique d'un dispositif d'imagerie radiographique selon une forme de réalisation de l'invention, permettant d'effectuer une prise de vue de face et une prise de vue de 15 côté par radiographie du patient et d'effectuer une numérisation 3D de l'enveloppe externe du patient, - la figure 2 est une représentation schématique d'une image numérisée en 3D de l'enveloppe externe du patient, - la figure 3 est une représentation similaire à la figure 2 sur laquelle 20 sont représentés des repères externes, - la figure 4 est une représentation schématique d'un modèle générique en 3D d'un squelette, - la figure 5 est une représentation schématique d'un modèle en 3D calculé du squelette du patient et de ses projections sur les images 25 radiographiques, - les figures 6a et 6b sont des représentations schématiques en coupe transversale montrant les déformations du modèle 3D, - la figure 7 représente les projections des points sur une radiographie et une zone de recherche, 30 - la figure 8 représente un modèle en 3D reconstruit du squelette du patient. La figure 1 représente un dispositif d'imagerie radiographique 1, notamment pour la reconstruction tridimensionnelle, comportant une armature mobile 2 déplaçable verticalement de façon motorisée sur des 35 guides verticaux 3, dans une direction de translation 4. Cette armature entoure un champ d'observation 5 dans lequel peut prendre place un patient P, par exemple debout, ce qui permet 11 d'observer la position des os du squelette de ce patient en station debout, ce qui peut être important par exemple pour les patients atteints de scoliose. Bien entendu, le patient pourrait être dans une position assise si nécessaire. 5 L'armature mobile 2 porte une première source radioactive 6a et une deuxième source radioactive 6b. Chaque source radioactive 6a, 6b est associée à un détecteur 7a, 7b, respectivement qui est disposé face à la source 6a, 6b, au-delà du champ d'observation 5. Chaque détecteur 7a, 7b comporte au moins une ligne horizontale 8a, 8b de cellules de 10 détection. Les sources 6a, 6b et les détecteurs 7a, 7b font partie d'un premier système d'imagerie radiographique 6. Les sources radioactives 6a, 6b sont adaptées pour émettre des rayons ionisants, notamment des rayons X, dans les directions de prise de vue 9a, 9b respectivement, ces rayonnements ionisants étant aptes à 15 être détectés par les détecteurs 7a, 7b. La direction de prise de vue 9a est antéropostérieure par rapport au patient P alors que la direction de prise de vue 9b est latérale par rapport au patient P. De manière classique, les rayonnements ionisants émis par les sources 6a, 6b traversent chacun une fente horizontale 10a, 10b 20 ménagée dans un réticule 11a, 1 lb tel qu'une plaque métallique, pour générer un faisceau horizontal 12a, 12b de rayonnements ionisants dans le champ d'observation 5. Bien entendu, les sources radioactives et les détecteurs pourraient être en nombre supérieur à 2, et les directions de prise de 25 vue de ces différentes sources radioactives pourraient, le cas échéant, ne pas être perpendiculaires entre elles, ni même horizontales. Le dispositif comprend en outre au moins un capteur optique 13 pour réaliser la numérisation externe du patient P. Tout type de capteur optique d'image est utilisable, tels que les capteurs CCD (Couple Charge 30 Device) et CMOS (Complementary Metal Oxyde Semiconductor) ou encore les capteurs aptes à mesurer une profondeur optique. A noter cependant qu'il est préférable d'utiliser un capteur optique présentant une précision et résolution de maillage élevée. Dans l'exemple, ce capteur optique 13 est monté sur l'un des 35 guides verticaux 3. Le patient est en en outre debout sur un plateau 14 qui peut être mis en rotation au moyen d'un moteur. Une image en 3D du patient peut ainsi être prise en faisant tourner le plateau 14 sur 12 360°. Cette image 3D pourrait être prise en positionnant plusieurs capteurs optiques 13 autour du patient P au lieu de le faire tourner. 5 ou 6 capteurs optiques répartis autour du patient pourraient par exemple être utilisés pour numériser l'enveloppe externe du patient P, 5 notamment simultanément à la prise des radiographies. Dans une autre variante, l'image 3D pourrait être prise en déplaçant un unique capteur optique 13 autour du patient, soit manuellement, soit automatiquement au moyen d'un bras motorisé, d'un rail, ou autre. Le ou les capteurs optiques 13 font partie d'un deuxième système 10 d'imagerie 15. Les deux détecteurs 7a, 7b et le capteur optique 13 sont reliés à un dispositif de traitement d'image radiographique 20, par exemple un micro-ordinateur ou autre système électronique de commande tel que, par exemple un ou plusieurs processeurs de type microprocesseur, 15 microcontrôleur ou autre. Le micro-ordinateur 20 comporte des moyens d'exécution d'un programme d'ordinateur adapté pour mettre en oeuvre le procédé décrit dans la présente invention. Il comprend ainsi des moyens de réception 21 de données d'image pour recevoir des données d'image 20 représentatives d'une ou plusieurs images radiographiques et d'une image en trois dimensions d'une enveloppe externe de l'objet. Il comprend aussi des moyens de traitement 22 raccordés aux moyens de réception 21 agencés pour exécuter les étapes (c) à (g) précédemment décrites du procédé selon l'invention. Enfin, il comprend des moyens de 25 transmission 23, notamment vers un dispositif d'affichage 24, ici un écran, raccordés aux moyens de traitement 22. Il peut également comporter des moyens de mémorisation 25 reliés aux moyens de traitement 22, qui peuvent être une mémoire vive ou une mémoire RAM (de l'anglais « Random Access Memory »), une EEPROM (de l'anglais 30 « (Electrically-Erasable Programmable Read-Only Memory »), ou autre. Ces moyens de mémorisation peuvent notamment mémoriser les différents modèles et données, et éventuellement le programme d'ordinateur. Le micro-ordinateur 20 est en outre ici équipé : d'une interface 35 d'entrée comprenant au moins un clavier 26 et généralement une souris (non représentée), et d'une interface de sortie comprenant au moins l'écran 24 et généralement une imprimante (non représentée). 13 Le micro-ordinateur 20 peut également être relié aux moyens d'entraînement motorisés (non représentés) contenus dans les guides 3 et dans le plateau tournant 14, aux sources 6a, 6b et au capteur optique 13, de façon à commander d'une part le déplacement vertical de 5 l'armature 2 et l'émission des rayonnements ionisants et d'autre part la rotation du plateau tournant 14 et le fonctionnement du capteur optique 13. Les détecteurs 7a, 7b peuvent par exemple être des détecteurs gazeux sensibles aux basses doses de rayonnements, par exemple du 10 type décrit dans les documents FR-A-2 749 402 ou FR-A-2 754 068. Bien entendu, d'autres types de détecteurs pourraient éventuellement être utilisés dans le cadre de la présente invention. La détection peut être réalisée par un détecteur purement linéaire (une ligne d'image à la fois), ou par un détecteur matriciel de n'importe quel rapport de forme. 15 Le fonctionnement du dispositif qui vient d'être décrit est maintenant détaillé en référence aux figures 2 à 8. Au moyen du micro-ordinateur 20, on fait prendre deux images radiographiques numériques I2Da, I2Db du patient P, par exemple en faisant balayer le champ d'observation 5 par les faisceaux 12a, 12b de 20 rayonnements ionisants sur une hauteur correspondant à la hauteur de la zone du patient à observer, par exemple le rachis et le bassin, voire l'ensemble du squelette. A cet effet, l'armature 2 est de préférence déplaçable sur une hauteur suffisante, par exemple de 70 cm ou plus, voire de 1 mètre ou plus. Suite à une calibration appropriée des sources 25 6a, 6b et des détecteurs 7a, 7b, le système d'imagerie radiographique 6 génère des premières données d'image des deux radiographies dans un même référentiel à trois dimensions. Avant ou après la prise des images radiographiques numériques, toujours au moyen du micro-ordinateur 20, on prend une image 30 numérique en trois dimensions I3D (fig. 2) de l'enveloppe externe du patient P en faisant tourner, notamment sur 360°, le plateau 14 supportant le patient P. Le deuxième système d'imagerie 15 génère alors des deuxièmes données d'image dans un deuxième référentiel à trois dimensions différent du précédent, les prises de vues ne pouvant être 35 simultanées. Il est alors nécessaire de transformer les premières ou les deuxièmes données d'image pour les exprimer dans un même référentiel de référence, qui peut être l'un des référentiels dans lequel sont 14 générées les premières ou deuxièmes données d'image ou un autre référentiel. Cette transformation peut être réalisée (par les moyens de traitement 22) au moyen d'une matrice de transformation tel que précédemment décrit. Pour une prise de cette image numérique en 3D 5 simultanée aux radiographies, au lieu d'utiliser un plateau tournant, on pourrait utiliser plusieurs capteurs optiques 13 disposés autour du patient P. Une calibration appropriées des sources 6a, 6b et du capteur optique 13 permettrait alors de générer les premières et deuxièmes données d'image dans un même référentiel à trois dimensions. 10 Les images numériques I2Da, I2Db, I3D de la partie examinée du patient sont transmises aux moyens de traitement 22. On peut éventuellement les enregistrer dans la mémoire du micro-ordinateur 20 et les visualiser sur l'écran 24 du micro-ordinateur. Dans l'exemple, les images radiographiques numériques I2Da, I2Db sont des images 15 antéropostérieures et latérales (fig. 5). On estime ensuite une position géométrique dans le référentiel de référence en 3D d'un ou plusieurs points d'intérêt (premiers repères) appartenant à la structure interne. Ces premiers repères sont caractéristiques de la structure interne. 20 Cette estimation passe de préférence par l'estimation de la position géométrique de repères externes, autrement dit de repères faisant partie de l'enveloppe externe. La localisation de repères externes sur l'enveloppe externe de l'objet présente l'avantage de permettre une estimation sans radiation de la structure interne de cet objet. 25 L'estimation de la position géométrique des repères externes est de préférence automatique, par exemple en utilisant des méthodes connues, telles que celles décrites par Michonski et al. (Automatic recognition of surface landmarks of anatomical structures of back and posture. J Biomed Opt. 2012 May;17(5):056015). 30 On peut aussi employer des courbures de Gauss pour identifier des repères externes qui pourront ensuite être utilisés pour estimer la position des premiers repères internes. Tel que représenté sur la figure 3, on peut identifier préférentiellement les repères externes suivants : 35 R1 : vertebra prominens de C7, R2 : sommet du pli inter-fessier, R3 et R4 : épaules gauche et droite (points immédiatement à la verticale des plis sous-axillaire), R5 et R6 : plis sous-axillaires gauche et droit, 15 R7, R8 : points gauche et droit de la taille, R9, R10 : épines iliaques postéro-supérieures gauche et droite, R11, R12 angles inférieurs de la scapula gauche et droit. Les paramètres anatomiques tels que les dimensions des 5 vertèbres et des structures osseuses peuvent être estimées par des équations de régressions liées à des mesures anthropométriques externes, tel que décrit par exemple par Bertrand et al., Estimation of external and internal human body dimensions from few external measurements, Journal of Musculoskeletal Research, Vol. 12, No. 4 10 (2009) 191-204). Une telle estimation des longueurs caractéristiques de repères anatomiques de la structure interne permet de calculer la position possible de ce repère interne (appelé ici premier repère). Plusieurs études ont également proposé des relations externes-internes pour estimer les centres des corps vertébraux à partir de la 15 surface externe (The prédiction of lumbar spine geometry: method development and validation. Campbell-Kyureghyan N, Jorgensen M, Burr D, Marras W. Clin Biomech (Bristol, Avon). 2005 Jun;20(5):455- 64). Différentes lois de transformation non-rigides peuvent aussi être 20 envisagées pour estimer la position des repères de la structure interne (Seo et al. An automatic modeling of human bodies from sizing parameters. In: SI3D '03: Proceedings of the 2003 symposium on Interactive 3D graphics, New York, NY, USA, ACM Press (2003) 19-26). Le procédé est ici présenté en faisant référence à la colonne 25 vertébrale, mais un procédé similaire pourrait également être utilisé pour une structure interne simple dont on aurait une connaissance a priori, tel qu'un os simple, tel une seule vertèbre, une structure ligamentaire, ou autre, ou un ensemble de structures anatomiques tel qu'un membre inférieur ou supérieur, ou autres. On utilisera alors 30 d'autres repères anatomiques adaptés aux structures à observer/imager. En outre, en fonction de la structure interne à observer/ imager, la numérisation de l'enveloppe externe peut être réalisée sur le corps complet du patient ou sur une partie limitée du corps (uniquement le dos du tronc par exemple). 35 Dès lors que les coordonnées (dans le référentiel de référence) des premiers repères sont estimées et que les valeurs des paramètres anatomiques ont éventuellement été estimées, on peut calculer la forme 16 à trois dimensions d'un modèle représentant la structure interne à partir d'un modèle générique de celle-ci. A cet effet, on dispose de modèles génériques des structures internes représentés sur les images radiographiques, en particulier les 5 vertèbres, mais éventuellement aussi d'autres structures anatomiques internes, tels les ligaments ou autres. Ces modèles peuvent être enregistrés dans le micro-ordinateur 20. Un modèle générique, par exemple établi à partir d'une base de données, peut être défini comme un maillage de quelques centaines à 10 quelques centaines de milliers de points d'une structure. De tels modèles génériques sont par exemple établis à partir de bases de données contenant des données relatives à des repères particuliers de la structure. Ces données peuvent comprendre des positions de points caractéristiques de la structure, des longueurs 15 caractéristiques de la structure, ou encore des segments, droites ou arcs caractéristiques de la structure, et/ou des contours et arêtes de la structure. Pour une vertèbre, la base contient par exemple la position d'une vingtaine de points caractéristiques de la vertèbre, des longueurs caractéristiques de la vertèbre. 20 La base de données peut également contenir des données relatives à l'emplacement relatif de la structure interne dans le squelette du sujet dont elle est issue. Dans le cas d'une vertèbre, il s'agit par exemple de l'orientation angulaire de la vertèbre et la courbure de la colonne au niveau de la vertèbre. 25 Il est possible d'établir des classements des objets dans différentes catégories afin de distinguer des données caractéristiques d'individus sains ou atteints d'une pathologie, et/ou caractéristiques du poids, de la taille, de l'âge ou de tout autre type de paramètre d'un individu. 30 Un modèle générique peut comporter par exemple des données statistiques (moyennes, variances, ...) pour chaque paramètre de la base de données, ou bien des équations mathématiques permettant de déterminer pour un objet donné la localisation de points caractéristiques à partir de la valeur de paramètres estimateurs de cet 35 objet. Dans un mode de réalisation, la forme à trois dimensions d'un modèle représentant la structure interne est calculée à partir d'un 17 modèle générique de celle-ci par exemple en paramétrant ce dernier au moyen des premiers repères et éventuellement des paramètres anatomiques qui ont été estimées, de préférence automatiquement, à partir de l'image numérique I3D de l'enveloppe externe. 5 La figure 4 représente schématiquement le modèle générique M de la structure osseuse du sujet, la figure 5 montrant le modèle M1 calculé obtenu par paramétrage du modèle M. La scoliose du sujet observé est ainsi visible sur le modèle calculé Ml. Ce paramétrage du modèle générique M peut être réalisé en 10 modifiant les dimensions du modèle générique conformément aux paramètres anatomiques. On peut par exemple utiliser la méthode décrite par Hwang et al. ("Rapid Development of Diverse Human Body Models for Crash Simulations through Mesh Morphing," SAE Technical Paper 2016-01- 15 1491, 2016). Pour cela, on fait correspondre les coordonnées d'un ou plusieurs points du modèle générique M correspondants à des premiers repères avec les coordonnées de ces premiers repères estimés à partir des données de l'image I3D. Ce positionnement peut être réalisé par rapport 20 aux positions des vertèbres associées à chaque paire de côtes. Puis, on déforme le modèle générique géométriquement en conservant les longueurs pour correspondre aux contraintes données par l'enveloppe externe. Un modèle mécanique peut également être utilisé à cette étape (Closkey, R. F., Schultz, A. B. and Luchies, C. W. (1992) A model for 25 studies of the deformable rib cage. Journal of Biomechanics, 25, 529- 539) Par exemple, lorsque la structure à imager est le gril costal, le modèle générique M est d'abord mis aux dimensions du patient au moyen des premiers repères et des paramètres anatomiques pour 30 obtenir un modèle M', puis, il est positionné à l'intérieur de l'image numérisée I3D de l'enveloppe externe par rapport aux positions des vertèbres associées à chaque paire de côtes, tel que représenté sur la figure 6a. Le modèle M' est alors dans le même référentiel de référence que l'image I3D de l'enveloppe externe. Sur cette figure 6a, la lettre « V » 35 désigne une vertèbre, à laquelle sont reliées les côtes Cl et C2. Le trait plein correspond à l'image de l'enveloppe externe I3D, alors que les traits pointillés correspondent au modèle M'. Les côtes Cl, C2 du 18 modèle M' sont ensuite déformées géométriquement en conservant leurs longueurs pour correspondre aux contraintes données par la forme externe du tronc, en particulier les gibbosités thoraciques, tel que représenté sur la figure 6b. On obtient ainsi un modèle 3D calculé M1 5 de la structure interne à observer/imager. Lors de cette étape, les techniques d'animation de rigging (déformation du squelette par cinématique directe ou inverse) et de skinning (association de chaque os du squelette à une portion de la peau) peuvent être utilisées (Avatar reshaping and automatic rigging using a deformable model. Feng A., 10 Casas D., Shapiro A..Proceeding MIG '15 Proceedings of the 8th ACM SIGGRAPH Conference on Motion in Games p 57-64). On procède ensuite à la projection du modèle 3D calculé M1 sur chacune des images radiographiées I2Da, I2Db (fig.5). Cette projection permet un repérage de repères anatomiques (deuxièmes repères) dans 15 les radiographies 2D sans nécessiter de traitement d'image, ce qui permet d'accélérer ce repérage et de le rendre plus fiable. Le choix des deuxièmes repères peut être automatique ou manuel avec une désignation d'un opérateur. Ces deuxièmes repères peuvent être les mêmes que les premiers repères ou non. 20 La projection utilisée est par exemple une projection linéaire. Cette projection peut être effectuée en utilisant la technique DLT (acronyme anglais de « Direct Linear Transformation »- transformation linéaire directe). Il est également possible de modifier les équations de la technique DLT en procédant à une simplification. Par projection du 25 modèle Ml, on obtient ainsi les coordonnées en 2D estimées des deuxièmes repères dans chaque image radiographiée I2Da, I2Db. Cette estimation initiale est ensuite affinée (étape (e)(iii)), en définissant une zone de recherche située autour de la projection du deuxième repère du modèle calculé dans l'image radiographiée. La figure 7 représente les 30 projections d'une pluralité de ces points Pi sur une radiographie, ou i est un nombre entier représentant un nombre de points. Pour des raisons de clarté, seuls quelques points Pl, P2, P3 sont représentés. La zone Z1 correspond à la zone de recherche du repère de projection Pl. Cette étape commence par exemple par une segmentation des ROI 35 (Region Of Interest- regions d'intérêt). On peut alors appliquer des filtres permettant de repérer les zones de fortes densités négatives et positives (filtre orienté). Un filtre de 19 Canny appliqué à ces ROI permet ensuite de déterminer des contours qui peuvent être utilisés par la technique de transformée généralisé de Hough (utilisant un modèle de vertèbre). En variante, au lieu d'utiliser les filtres orientés et de Canny, 5 l'amélioration des estimations initiales peut se faire automatiquement par une méthode de contour actif comme le flux de vecteur gradient (GVF) (Moura et al.. Fast 3D reconstruction of the spine from biplanar radiographs using a deformable articulated model. Med Eng Phys. 2011 Oct;33(8):924-33). 10 Lorsque les projections 2D du modèle 3D M1 ont été optimisées en exploitant les données radiographiques, le modèle 3D M1 peut être reconstruit. A cet effet, on corrige par exemple le modèle M1 par Transformation Linéaire Directe (DLT) en utilisant les mêmes couples de points (avec leurs coordonnées actualisées) que lors de la projection 1 5 (deuxièmes repères). Un modèle reconstruit M2 est alors obtenu (fig. 8). Cette reconstruction est réalisée en utilisant des deuxièmes repères qui sont stéréocorrespondants. La géométrie épipolaire, utilisant des repères non stéréo correspondants peut être aussi utilisée pour obtenir une reconstruction plus précise du modèle M2 (méthode 20 de type NSCP (Non Stereo Correspondent Point)). Par exemple, après la reconstruction des points stéréo-correspondants, le modèle M2 est initialisé. Les points non-stéréo-correspondants du modèle M2 seront projetés uniquement dans une des deux vues radiographiques (I2Da ou ID2b) dans laquelle on sait que la recherche sera possible. Le modèle 25 M2 est réactualisé par l'intégration des points non-stéréo- correspondants. La présente invention permet ainsi de reconstruire un modèle 3D précis d'une structure interne en limitant la dose de rayonnement émise vers le champ d'observation. REVENDICATIONS1. Procédé d'imagerie radiographique d'une structure interne en trois dimensions faisant partie d'un objet situé à l'intérieur d'un champ d'observation, le procédé comprenant : a) recevoir des premières données d'image générées dans un référentiel de référence à trois dimensions, représentatives d'au moins une image radiographique en deux dimensions de la structure interne de l'objet, b) recevoir des deuxièmes données d'image générées dans le référentiel de référence, représentatives d'une image en trois dimensions d'une enveloppe externe de l'objet, c) estimer à partir des deuxièmes données d'image les coordonnées dans le référentiel de référence de premiers repères caractéristiques de la structure interne, et, optionnellement, la valeur d'au moins un paramètre caractéristique de la structure interne, d) calculer, dans le référentiel de référence, la forme à trois dimensions d'un modèle représentant ladite structure interne à partir d'un modèle générique connu a priori de ladite structure interne, ce modèle générique comportant des points correspondants aux premiers repères, le modèle calculé étant obtenu par déformation du modèle générique de façon que ledit modèle calculé suive une forme la plus proche possible d'une isométrie du modèle générique tout en maintenant en coïncidence les points du modèle générique déformé avec les premiers repères, et optionnellement en appliquant au modèle calculé le ou les paramètres caractéristiques de la structure interne, e) choisir au moins un deuxième repère caractéristique de la structure interne à localiser sur l'image radiographique, et pour chaque deuxième repère : i) identifier un point appartenant au modèle calculé à l'étape d) correspondant au deuxième repère, (ii) calculer les coordonnées d'une projection de ce point sur l'image radiographique, iii) déterminer à partir des premières données d'image les coordonnées du deuxième repère sur l'image radiographique en le recherchant dans une zone de recherche située autour des coordonnées de la projection. 21 2. Procédé d'imagerie radiographique d'une structure interne en trois dimensions faisant partie d'un objet situé à l'intérieur d'un champ d'observation, le procédé comprenant : a) recevoir des premières données d'image générées dans un 5 référentiel à trois dimensions, représentatives d'au moins une image radiographique en deux dimensions de la structure interne de l'objet, b) recevoir -des deuxièmes données d'image générées dans un autre référentiel à trois dimensions, représentatives d'une image en trois dimensions d'une enveloppe externe de l'objet, 10 dans lequel au moins une des étapes de réception (a) et (b) comprend une étape de transformation des données d'image reçues pour exprimer les premières et les deuxièmes données d'image dans un même référentiel de référence, c) estimer à partir des deuxièmes données d'image les 15 coordonnées dans le référentiel de référence de premiers repères caractéristiques de la structure interne, et, optionnellement, la valeur d'au moins un paramètre caractéristique de la structure interne, d) calculer, dans le référentiel de référence, la forme à trois dimensions d'un modèle représentant ladite structure interne à partir 20 d'un modèle générique connu a priori de ladite structure interne, ce modèle générique comportant des points correspondants aux premiers repères, le modèle calculé étant obtenu par déformation du modèle générique de façon que ledit modèle calculé suive une forme la plus proche possible d'une isométrie du modèle générique tout en 25 maintenant en coïncidence les points du modèle générique déformé avec les premiers repères, et optionnellement en appliquant au modèle calculé le ou les paramètres caractéristiques de la structure interne, e) choisir au moins un deuxième repère caractéristique de la structure interne à localiser sur l'image radiographique, et pour chaque 30 deuxième repère : i) identifier un point appartenant au modèle calculé à l'étape d) correspondant au deuxième repère, (ii) calculer les coordonnées d'une projection de ce point sur l'image radiographique, 35 iii) déterminer à partir des premières données d'image les coordonnées du deuxième repère sur l'image radiographique en le 22 recherchant dans une zone de recherche située autour des coordonnées de la projection. 3. Procédé d'imagerie selon la revendication 1 ou 2, dans lequel, les deuxièmes données d'image reçues à l'étape b) sont générées 5 par un système d'imagerie comprenant au moins un capteur optique. 4. Procédé d'imagerie selon l'une des revendications 1 à 3, dans lequel, l'étape c) comprend : (i) la détermination des coordonnées de repères externes positionnés sur l'enveloppe externe dans le référentiel de référence, obtenue : 10 par calcul au moyen d'un modèle mathématique, ou par traitement des deuxièmes données d'image, ces repères externes étant marqués sur l'enveloppe externe, (ii) l'estimation des coordonnées des premiers repères dans le même référentiel à partir des repères externes. 15 5. Procédé d'imagerie selon l'une quelconque des revendications 1 à 4, dans lequel, les premières données d'image reçues à l'étape (a) sont générées par balayage en déplaçant au moins une source radioactive dans une direction de translation non parallèle à une direction de prise de vue. 20 6. Procédé d'imagerie selon l'une quelconque des revendications précédentes, dans lequel, au cours de l'étape a), on reçoit des premières données d'image représentatives de deux images radiographiques enregistrées selon deux directions de prise de vue non parallèles. 25 7. Procédé d'imagerie selon la revendication 6, dans lequel, au cours de l'étape e), on choisit les deuxièmes repères parmi des deuxièmes repères à localiser sur chacune des deux images radiographiques, et optionnellement, des deuxièmes repères à localiser sur une seule image radiographique, 30 le procédé comprenant en outre : f) calculer les coordonnées des deuxièmes repères dans le référentiel de référence à partir de leurs coordonnées sur au moins une des images radiographiques déterminées à l'étape e), g) corriger le modèle calculé à l'étape d) en faisant coïncider des 35 points du modèle calculé correspondants aux deuxièmes repères avec les coordonnées de ceux-ci calculés à l'étape f). 23 8. Procédé d'imagerie selon la revendication 7, dans lequel les premières données d'image reçues à l'étape a) sont générées simultanément, notamment dans un même référentiel à trois dimensions, par balayage, en déplaçant en synchronisme, dans une 5 même direction de translation non parallèle aux directions de prise de vues, deux sources radioactives émettant deux faisceaux de rayons ionisants respectivement dans les deux directions de prise de vue. 9. Procédé d'imagerie selon l'une quelconque des revendications 1 et 3 à 8 lorsqu'elles ne dépendent pas de la 10 revendication 2, dans lequel les premières données d'image et les deuxièmes données d'image sont générées simultanément, notamment dans un même référentiel à trois dimensions. 10. Procédé d'imagerie selon l'une quelconque des revendications 1 à 9 comprenant en outre : 15 une étape de transmission, notamment vers un dispositif d'affichage, des données choisies parmi les premières données d'image, les deuxièmes données d'image, des données représentatives des coordonnées des premiers repères, des données représentatives des coordonnées de chaque deuxième repère déterminé à l'étape e)iii), et, le 20 cas échéant des données représentatives du modèle corrigé calculé à l'étape g). 11. Dispositif de traitement d'image radiographique (20) comprenant : - des moyens (21) de réception de données d'image pour 25 recevoir des premières données d'image générées dans un référentiel à trois dimensions, représentatives d'au moins une image radiographique en deux dimensions de la structure interne de l'objet, et des deuxièmes données d'image générées dans un référentiel à trois dimensions représentatives d'une image en trois dimensions d'une enveloppe 30 externe de l'objet, - des moyens (22) de traitement raccordés aux moyens de réception agencés pour effectuer au moins les étapes c) à e) du procédé selon l'une quelconque des revendications 1 à 10, - optionnellement des moyens de transmission (23), 35 notamment vers un dispositif d'affichage (24). 12. Dispositif d'imagerie radiographique (1) d'une structure interne en trois dimensions faisant partie d'un objet comprenant : 24 - un premier système d'imagerie radiographique (6) agencé pour prendre au moins une image radiographique en deux dimensions de la structure interne de l'objet au moyen d'au moins une source radioactive (6a, 6b) et pour générer dans un référentiel à trois 5 dimensions des premières données d'image représentatives de l'image radiographique prise, - un deuxième système d'imagerie (15) comprenant au moins un capteur optique (13), agencé pour prendre une image en trois dimensions d'une enveloppe externe de l'objet et pour générer dans un 10 référentiel à trois dimensions des deuxièmes données d'image représentatives de cette image, - le dispositif de traitement d'image radiographique (20) selon la revendication 11 connecté aux premier et deuxième systèmes d'imagerie pour recevoir les premières et deuxièmes données d'image. 15 13. Programme d'ordinateur comprenant des portions de code de programme pour l'exécution des étapes du procédé d'imagerie selon l'une quelconque des revendications 1 à 10, lorsque ledit programme est exécuté par un dispositif de traitement d'image radiographique (20).