Brevet FR2823968A1 — publié le 20021031

l'espace des k. - - Procédé d'étalonnage d'un système d'imagerie, support de mémoire et dispositif associé La présente invention concerne le domaine de l'étalonnage d'une châîne de traitement d'images pour corriger des distorsions, en particulier dans des images radiologiques bidimensionnelles, par exemple dans le domaine de la radiologie. De faç on connue, le s appareil s de radiol ogie comprennent un moyen d'émission d'un faisccau de rayons X ou ensemble radiogène tel qu'un tube à rayons X et un moyen de réception dudit faisccau tel qu'un détecteur à l'état solide ou encore un scintillateur et une caméra vidéo, par exemple de type CCD. Le moyen d'émission et le moyen de réception de rayons X sont généralement supportés par un système mobile à un ou plusieurs axes, permettant des prises de vue selon différentes incidences. Le moyen de réception du faisccau de rayons X est relié à des moyens de traitement d'images permettant la génération d'images tridimensionnelles à partir d'une série d'images bidimensionnelles prises par le moyen de réception. I1 s'avère important d'obtenir une bonne correspondance entre une image tridimensionnelle et une image bidimensionnelle prise lors d'une étape au cours de laquelle le patient subit une action particulière telle que la pose d'un cathéter dans le domaine de l'angiographie, pour pouvoir suivre le déplacement du cathéter dans l'image bidimensionnelle mais également dans l'image tridimensionnelle. A cet effet, un étalonnage de la machine de radiologie est effectué lors de sa mise en service puis à intervalles réguliers, pour tenir compte des imperfections de différents éléments de la machine de radiologie, tels que la déformation des supports du moyen d'émission et du i moyen de réception de rayons X, les imprécisions des capteurs, les défauts du récepteur, etc. Jusqu'à présent, un étalonnage pour certaines configurations de la machine est effectué à intervalles courts, par exemple une semaine, d'o une immobilisation de la machine et un coût de main d'_uvre importants. La présente invention propose, notamment, un procédé d'étalonnage perfectionné permettant une précision améliorce des images obtenues pour toutes les configurations de distorsion et des paramètres intrinsèques de la géométrie conique de la machine. La présente invention propose, notamment, un procédé d'étalonnage économique réduisant la durce des opérations d'étalonnage d'entretien effectuces régulièrement au cours de la vie d'une machine d'imagerie. Le procédé d'étalonnage, selon un aspect de l'invention, est destiné à un système d'imagerie, comprenant un moyen d'émission d'un faisccau d'énergie, un moyen de réception du faisccau d'énergie après qu'il a traversé un organe à étudier, un moyen de commande du moyen de réception, un moyen de support dudit organe et un moyen de mise en forme d'images. On étalonne des paramètres intrinsèques relatifs aux différents types d'images pouvant être prises dans une configuration prédétermince du système, les paramètres intrinsèques étant estimés à l'aide d'un fantôme bidimensionnel pour établir un champ de distorsion pour une configuration particulière, le champ de distorsion pour une configuration quelconque étant calculé par combinaison d'une pluralité de champs de distorsion de causes supposces indépendantes. Dans un mode de réalisation de l'invention, on détermine un champ de distorsion Do dû au champ électromagnétique terrestre à partir d'une image dudit fantôme prise dans une position centrale du système d'imagerie selon chaque axe de rotation, ledit champ de distorsion Do étant indépendant de l'angulation. Le champ de distorsion Do sera recalculé fréquemment pour tenir compte des variations du champ électromagnétique terrestre. Au contraire, les autres distorsions évoluent de façon plus lente et pourront être recalculées moins fréquemment, par exemple annuellement. Dans un mode de réalisation de l'invention, on étalonne des paramètres intrinsèques relatifs aux différents types d'images pouvant être prises dans des positions angulaires prédéterminces du moyen d'émission et du moyen de réception, les paramètres intrinsèques étant estimés à l'aide d'un fantôme bidimensionnel pour établir un champ de distorsion pour une angulation particulière, le champ de distorsion pour une angulation quelconque étant calculé par interpolation polynomiale des champs de distorsion établis chacun pour un point de contrôle, un point de contrôle étant défini par une angulation particulière par axe de rotation du système d'imagerie, devant faire l'objet d'un étalonnage. De préférence, les points de contrôle sont définis de façon telle qu'ils recouvrent totalement et régulièrement le domaine des orientations possibles de la machine considérce. Si autour d'un axe, la rotation est possible sur plus de 180 , il suffira de prendre -90 , 0 et +90 comme angulations particulières, 0 étant la position médiane par rapport aux extremums. - De préférence, la fonction d'interpolation est de type biquadratique ou bicubique. Avantageusement, on établit une relation entre des paramètres d'interpolation u, v et l'angulation, la relation étant de type u(x) = ax5+bx3 + cx, x étant une coordonnce d'un point, de façon que le champ de distorsion soit déterminé en fonction de l'angulation, la fonction d'interpolation permettant de déterminer la distorsion en fonction des paramètres d'interpolation u, v. Avantageusement, la fonction u(x) est de type u(x) = cc(A+,)x5+ %(A+)x3 + e(A+Ox, A étant une constante dépendant du type de système d'imagerie, a,, X, 8, E et p étant des constantes déterminces expérimentalement. Plus particulièrement, la fonction u(x) est de type u(x) =(A 1,5)x5 +(2,5-2A) x3 + Ax. A peut être supérieur ou égal à O et inférieur ou égal à 15/8. De préférence, on utilise un paramètre d'interpolation pour chaque angle définissant la position du système d'imagerie. A sera calculé pour chaque angle. Les points de contrôle seront au nombre de (n+1)2, avec n valant au minimum 2 pour l'interpolation biquadratique et 3 pour l'interpolation bicubique. Dans un mode de réalisation de l'invention, on détermine un champ de distorsion DI dû à une inversion d'image par inversion du sens de balayage du moyen de réception du faisccau d'énergie par comparaison d'une image dudit fantôme retournce après sa prise et d'une image dudit fantôme prise en balayage inversé, et ce pour les changements de balayage utilisés. Dans un mode de réalisation de l'invention, on détermine un champ de distorsion DR dû à un changement de résolution du moyen de réception du faisccau d'énergie par comparaison d'une image dudit fantôme prise avec une première résolution et d'une image dudit fantôme prise avec une deuxième résolution différente de la première, et ce pour les changements de résolution utilisés. Dans un mode de réalisation de l'invention, on détermine un champ de distorsion DM dû à un changement de matrice d'acquisition du moyen de réception du faisccau d'énergie par comparaison d'une image - dudit fantôme prise avec une première matrice d'acquisition et d'une image dudit fantôme prise avec une deuxième matrice d'acquisition - différente de la première, par exemple lors d'un changement d'une matrice d'acquisition de 51 2x5 12 pixcls à une matrice d'acquisition de 1 024x 1024 pixcls, et ce pour les changements de matrice d'acquisition utilisés. L'invention propose, également, un support de mémoire contenant des instructions pour mettre en _uvre un procédé d'étalonnage d'un système d'imagerie comprenant un moyen d'émission d'un faisccau d'énergie, un moyen de réception du faisccau d'énergie après qu'il a traversé un organe à étudier, un moyen de commande du moyen de réception, un moyen de support dadit organe et un moyen de mise en forme d'images, procédé dans lequel on étalonne des paramètres intrinsèques relatifs aux différents types d'images pouvant être prises dans une configuration prédétermince du système, les paramètres intrinsèques étant estimés à l'aide d'un fantôme bidimensionnel pour établir un champ de distorsion pour une configuration particulière, le champ de distorsion pour une configuration quelconque étant calculé par combinaison d'une pluralité de champs de distorsion de causes supposces indépendantes. L'invention propose, également, un système d'imagerie comprenant un moyen d'émission d'un faisccau d'énergie, un moyen de s réception du faisccau d'énergie après qu'il a traversé un organe à étudier, un moyen de commande du moyen de réception, et un moyen de mise en forme d'images. Le système comprend un moyen d'étalonnage des paramètres intrinsèques relatifs aux différents types d'images pouvant S être prises dans une configuration prédétermince du système, les paramètres intrinsèques étant estimés à l'aide d'un fantôme bidimensionnel pour établir un champ de distorsion pour une configuration particulière, le champ de distorsion pour une configuration quelconque étant calculé par combinaison d'une pluralité de champs de distorsion de causes supposces indépendantes. Dans un mode de réalisation de l'invention, le fantôme bidimensionnel comprend une grille. Dans un mode de réalisation de l'invention, on effectue un étalonnage lors de l'installation du système d'imagerie. I1 peut s'agir de 1S l'étalonnage des paramètres de la géométrie conique et de la constante A. - En outre, on peut étalonner des paramètres extrinsèques relatifs à la position dudit organe par rapport au système d'imagerie, du moyen de - réception et du moyen de support, pour déterminer la géométrie d'acquisition et le champ de distorsion quelle que soit la position du moyen d'émission et du moyen de réception, les paramètres extrinsèques étant calculés à l'aide d'un fantôme tridimensionnel pour des positions prédéterminces du moyen d'émission.Le même fantôme peut servir à l'étalonnage des deux types de paramètres. Avantageusement, les positions prédéterminces pour le calcul des paramètres intrinsèques tiennent compte de la distance entre moyen d'émission et moyen de réception, de la résolution, et de la matrice image. Dans un mode de réalisation de l'invention, au moins trois images d'un fantôme sont prises pour trois valeurs différentes de la distance entre moyen d'émission et moyen de réception. Dans un mode de réalisation de l'invention, la position du système d'imagerie est mesurce par des capteurs. Dans un mode de réalisation de l'invention, on effectue un étalonnage à intervalles de temps prédéterminés. L'invention permet de corriger des distorsions dans de nombreux types d'appareils de d'imagerie, notamment de radiologie avec des images bi ou tridimensionnelles. Une autre application serait simplement d'obtenir des images géométriquement correctes afin de pouvoir effectuer des mesures plus fiables dans les images bidimensionnelles. Ceci est également valable pour l'estimation des paramètres intrinsèques de la géométrie conique qui permettent de déduire, entre autres, le grossissement dans les images, ce qui aide au recalage 2D/3D mais également permet de traduire des mesures effectuces dans les images (donc donnces en pixcls) en des mesures métriques (donc en millimètres). On cherche à estimer le champ de distorsions présent dans les images radiographiques, quelle que soit la configuration de la machine lors de la prise d'image. L'intérêt de pouvoir estimer ce champ de distorsion est que l'on peut grâce à cette connaissance corriger les images des distorsions. Les images radiographiques sont alors géométriquement correctes, ce qui permet d'utiliser un modèle de projection conique dans un procédé de recalage 2D/3D, mais encore de pouvoir faire des mesures plus fiables dans les images 2D. Un procédé spécial est nocessaire, étant donné le grand nombre de configurations possibles de la machine. Le principe général du procédé est de retrouver la distorsion différentielle par rapport à un champ de distorsions Do de référence. Cette distorsion différentielle est due principalement à quatre causes intervenant lors d'un changement de configuration de la machine. 1) Changement de sens du balayage caméra 2) Changement de matrice d'acquisition 3) Changement de taille de champ de vue 4) Changement de l 'orientation spatiale de l' arccau portant l'émetteur et le récepteur. Une configuration de la machine peut se définir comme la donnce de l'état du balayage caméra (par exemple: normal, inverse horizontalement, inverse verticalement ou inverse horizonto verticalement), la taille de la matrice d'acquisition (par exemple: 256, 512 ou 1024 pixcls), la taille du champ de vue (par exemple: 1 1 cm, 16 cm, 22 cm ou 30 cm) et l'orientation de l'arccau dans l'espace (donnce par deux angles, par exemple les angles anatomiques). Une configuration de référence est définie en fixant: le sens du balayage caméra (par exemple normal), la taille de la matrice d'acquisition (par exemple 512 pixcls), la taille du champ de vue (par exemple 22 cm) et l'orientation de l'arccau dans l'espace (par exemple S pour le plan frontal, la position anatomique correspondant à la position de face: angles anatomiques RAO/LAO=0 et CRA/CAU=0 dogrés. Pour le plan latéral, nous prendrons de préférence la position de profil: angles anatomiques RAO/LAO=90 dogrés et CRA/CAU=0 dogrés). Le champ de distorsions Do est estimé à partir d'une image du fantôme de grille bidimensionnel acquise dans la confguration de référence. Par conception du fantome de grille, le champ de distorsion est équivalent à la donnce de la distorsion en tous les noeuds, c'est-à-dire un vecteur bidimensionnel donnant le décalage suivant les lignes et les 1S colonnes de l'image entre la position du noeud détecté dans l'image de - grille et la position théorique qu'il devrait avoir si l'image ne contenait pas de distorsion. Un champ de distorsion différentiel est la différence noeud à noeud entre deux champs de distorsions et est donc un champ de distorsion lui-même. On note "o" l'application d'un champ de distorsion à une image. Par exemple, soit I une image et D un champ de distorsions, l'image I'=D o I est l'image I distordue par le champ de distorsion D. De même, soit un autre champ D', l'image I" = D' o I' est l'image distordue par le champ de distorsion D'. Nous avons donc I" = D' o D o I = (D' o D) o I. Le champ de distorsion D" = D' o D est donc le champ de distorsion qu'il faut appliquer à I pour obtenir l'image I". On dit que D' est le champ de distorsion différentiel entre le champ D et le champ D". On suppose que quatre causes de variation de distorsion (changement, par rapport à la configuration de référence, dans le balayage caméra, dans la matrice d'acquisition, dans le champ de vue et/ou dans l'orientation de l'arccau) sont indépendantes. Nous pouvons alors écrire tout champ de distorsion D, observé dans les images acquises pour une configuration particulière de la machine comme étant composé de l'application successive du champ de référence Do' puis du champ différentiel dû au changement du balayage caméra DI' puis du champ différentiel dû au changement de la matrice d'acquisition DM, puis du champ différentiel dû au changement dans la taille du champ de vue DF' puis du champ différentiel dû au changement dans l'orientation de l'arccau DA, soit: D=DAoDF DM DI DO Le procédé permet de retrouver ces champs différentiels pour toutes les configurations possibles de la machine. Concernant les trois premières causes de variation de la distorsion (inversion caméra, matrice d'acquisition, champ de vue), les valeurs que peuvent prendre ces paramètres sont en très petit nombre (moins de 5 sur les machines actuelles). Invoquant leur indépendance, on peut estimer le champ de distorsion différentiel dû à un changement dans un seul d'entre eux, quelle que soit la valeur des autres paramètres. En particulier, on peut utiliser pour les autres paramètres les valeurs définies pour la configuration de référence. Par exemple, pour les changements dans l'inversion caméra, quatre configurations d'étalonnage peuvent être définies: I1: tous les paramètres sont les mêmes qu'en confguration de référence, sauf le balayage caméra inversé horizontalement. I2: tous les paramètres sont les mêmes qu'en configuration de référence, sauf le balayage caméra inversé verticalement. I3: tous les paramètres sont les mêmes qu'en configuration de référence, sauf le balayage caméra inversé verticalement et horizontalement. On estime les champs de distorsions Dl' D2 et D3 en acquérant une image de grille dans les configurations I1, I2 et I3 respectivement. On détermine ensuite les champs de distorsion différentiels DI1, DI2 et DI3 par comparaisons de D1, D2 et D3 respectivement avec DO. La comparaison revient à établir le vecteur de distorsion différentiel pour chaque noeud de la grille. Pour retrouver la distorsion D pour une configuration particulière, il suffira de regarder l'état du balayage caméra: S'il n'y a pas d'inversion caméra, la distorsion D1 sera nulle Si le balayage caméra se fait en inversion horizontale, DI = DII 35. Si le balayage caméra se fait en inversion verticale, DI = DI2 Si le balayage caméra se fait en inversion verticale et horizontale, DI = DI3. On procède de manière similaire pour établir des champs de distorsion différentiels correspondant aux changements de matrice d'acquisition et de taille de champ de vue. Le champ de distorsion, même lorsque l'inversion caméra, la matrice d'acquisition et le champ de vue sont gardés constants, va varier si l 'orientation de l' arccau est modifice. Ceci est dû au fait que la distorsion appara^t en partie parce que les électrons qui sont crcés dans l'amplificateur de brillance voient leur trajectoire modifiée par l'action du champ électromagnétique local, cette modification étant fonction de l'angle entre la direction du champ électromagnétique local et la trajectoire des électrons. Si l'orientation de l'amplificateur de brillance est modifiée, l'angle respectif entre le champ électromagnétique local et la trajectoire des électrons est également modifié et une variation de - distorsion s'ensuit. A la différence des trois premiers paramètres, l'orientation de l'arccau prend des valeurs continues (donc une infinité de valeurs potentielles). On ne peut donc pas a priori estimer un champ de distorsion différentiel pour chaque orientation: on le fait pour un nombre fini de positions, dites de contrôle, et on en déduit le champ différentiel pour toute autre position quand on en a besoin, et non pas pendant l'étalonnage, par interpolation biquadratique. Un mode de réalisation particulier de l'invention est illustré par les dessins annexés: - la figure 1 est une vue en perspective d'un appareil de radiologie à trois axes qui peut être utilisé pour mettre en oeuvre le procédé; - la figure 2 est une vue d'une grille d'étalonnage telle qu'observée sur un moyen de visualisation d'images associé au système de radiologie; - la figure 3 est une courbe de la fonction u(x) pour plusieurs paramétrages; et - la figure 4 est un organigramme des étapes de procédé. Comme on peut le voir sur la figure 1, l'appareil de radiologie comprend un pied 1 en forme de L, avec une base 2 sensiblement horizontale et un support 3 sensiblement vertical fixé à une extrémité 4 de la base 2. A l'extrérnité opposce 5, Ia base 2 comprend un axe de rotation parallèle au support 3 et autour duquel le pied est capable de tourner. Un bras de support 6 est fixé par une première extrémité au sommet 7 du support 3, de façon rotative selon un axe 8. Le bras de support 6 peut présenter la forme d'une baionnette. Un bras 9 circulaire en forme de C est maintenu par une autre extrémité 10 du bras de support 6. Le bras en C 9 est apte à coulisser de façon rotative autour d'un axe 13 par rapport à l'extrémité 10 du bras de support 6. Le bras en C 9 supporte un tube à rayons X 11 et un détecteur d'énergie 12 en positions diamétralement opposces se faisant face. Le détecteur 12 comprend une surface plane de détection. La direction du faisccau de rayons X est détermince par une ligne droite joignant un point focal du tube 11 au centre de la surface plane du détecteur 12. Les trois axes de rotation du pied 1, du bras de support 6 et du bras en C 9 sont sécants en un point 14. En position moyenne, ces trois axes sont matuellement perpendiculaires. Une table 15, prévue pour recevoir un patient, possède une orientation longitudinale alignce avec l'axe 8 en position de repos. A titre d'exemple, un jeu d'images d'un vaisscau sanguin est acquis pour différentes positions de l'appareil de radiologie. Une image tridimensionnelle peut être ensuite reconstruite à partir du jeu d'images bidimensionnelles. Des vues bidimensionnelles de l'image tridimensionnelle peuvent être produites de façon interactive. Sur la figure 2, est montré un exemple d'une image d'une grille d'étalonnage telle que qu'on peut l'observer sur un écran de visualisation d'un appareil de radiologie. La grille est disposoe à la place d'un objet à radiographier sur le trajet du faisccau de rayons X entre le tube à rayons X 11 et le détecteur de rayons X 12. La grille est à mailles carrces en matériau opaque ou relativement opaque aux rayons X, par exemple en plomb ou en cuivre. On constate que l'image de la grille n'en constitue pas la reproduction fidèle. Si le centre de la grille est convenablement reproduit, les bords présentent des distorsions que l'on cherche à calculer. En un point donné de la grille, la distorsion peut être caractérisce par un vecteur ayant une direction et un module. Parmi les causes de distorsion, on constate que certaines, comme les défauts du détecteur, sont indépendantes des angulations autour des axes de l'appareil de radiologie tandis que d'autres dépendent desdites angulations. Un étalonnage idéul nocessiterait de calculer la distorsion en tout point du champ image, pour toute angulation selon chacun des axes, et plus généralement pour toute configuration de la machine ce qui nocessiterait d'énormes ressources de calcul. On cherche à obtenir un étalonnage satisfaisant à coût raisonnable et relativement rapidement. Dans ce but, on effectue un calcul du champ de distorsion pour quelques angulations dites de contrôle et on en déduit un champ de distorsion pour les autres angulations. Pour simplifier les calculs, on considérera la position spatiale du tube 11 et du 1S détecteur 12, non par rapport aux axes de rotation de l'appareil de radiologie, mais par rapport à deux axes, de préférence orthogonaux entre eux et passant par l'isocentre, perpendiculaires à la direction du faisccau de rayons X, en d'autres termes définissant un plan parallèle à l'axe 13. Le troisième axe d'un repère formé à partir desdits deux axes est la direction du faisccau de rayons X. Une rotation du tube 11 et du détecteur 12 autour de ce troisième axe n'a pas d'effet sur la distorsion. Les angulations par rapport auxdits deux axes sont notées et La prise d'images peut se modéliser géométriquement comrne une projection conique perturbée par un champ de distorsion bidimensionnel dans l'image. Le détecteur d'image 12 est sensible au champ électromagnétique environnant, ce qui introduit des distorsions de l'image par rapport au modèle de projection conique. Une fonction de projection est donc vue comme la combinaison d'une projection conique et d'une fonction de distorsion qui vont dépendre toutes deux de l'état courant de l'appareil. La distorsion est une fonction du plan image dans lui-même. On suppose qu'elle est suffisamment régulière. On peut alors l'estimer de façon précise en la calculant sur un nombre de points épars dans l'image, sa valeur sur les autres points en étant déduite par interpolation. Le fantôme utilisé pour l'étalonnage des distorsions consiste en une grille. Chaque intersection de la grille est appelée noeud. La position de chaque noeud est détectée dans l'image. La distorsion à chaque noeud est définie comme l'écart entre sa position dans l'image et celle qu'il doit avoir dans le modèle idéal de la grille. Quand la distorsion est corrigée, l'hypothèse de projection conique est applicable. Cette dernière est une fonction homographique de l'espace tridimensionnel dans lequel est exprimé l'obj et vers le plan image. Elle se traduit par l 'estimation d'une matrice 3 x 4 à une constante multiplicative près. Nous avons ainsi onze paramètres à estimer, cinq paramètres intrinsèques modélisant la chaîne image et six paramètres extrinsèques donnant la position respective de l'ensemble émetteur-récepteur de rayons X et de l'objet imagé dans l'espace. Cette estimation de la projection conique est effectuce grâce à un autre fantôme consistant en un cylindre en plastique transparent aux rayons X sur lequel sont portées des billes en plomb disposoes en hélice. - Cette hélice permet d'avoir une trentaine de points visibles dans l'image. La détection de ces points et leur appariement avec ceux du modèle parfait de l'hélice permet de déterminer les paramètres de la projection conique. En connaissant à priori la configuration de la machine dans laquelle sera prise une image, on peut alors déterminer avec une grande précision la fonction de distorsion grâce à la grille et la projection conique grâce à l'hélice. Par approximation, on considère que le champ magnétique terrestre est localement uniforme et varie peu dans un intervalle de temps de l'ordre de quelques jours. On considère également que le champ magnétique terrestre est régulier, ce qui permet de considérer un faible nombre d'échantillons ou positions de contrôle, donc un faible nombre d'images de la grille. Pour chaque axe définissant l'orientation, trois angulations de contrôle, voire quatre, s'avèrent utiles. I1 en résulte donc neuf, voire seize positions de contrôle définies par une angulation et une interpolation de type biquadratique, voire bicubique. L'angulation considérce peut être une angulation "machine" par rapport à chacun des deux axes ou, préférablement, une angulation anatomique avec deux angles (3, O différents des angles "machine" mais pouvant s'exprimer en fonction d'eux. Les positions de contrôle seront choisies régulièrement espacées et symétriques par rapport à l'origine, afin d'optimiser la densité des positions dans l'espace paramètre. Dans le cas de l'interpolation biquadratique, un ensemble de positions sera donc donné par des valeurs positives 60 et p0, qui définissent des points de contrôle dans l'espace (3, O: (0,0), (-00, O), (00' O), (O. -0), (O. +0), (430, q)o), (0o, (30, pO). Un nombre plus élevé de points de contrôle peut être prévu. Pour obtenir une précision satisfaisante, il convient d'éviter l'extrapolation, c'est-à-dire la génération de valeur de distorsion pour des angulations situces en-dehors de l'intervalle compris entre - 00 et (30, - 4,O et +0. On choisira donc les valeurs de (30 et de +0 correspondant à la totalité ou à une large plage de la gamme d'angulations possible de façon que l'ensemble des estimations soit effectué par interpolation, 60 et 0 restant inférieurs ou égaux à 90 . Le principe de l'interpolation quadratique est le suivant. Trois pointsdits de contrôle Pi = (xi, Yi) dans le plan peuvent être reliés par une courbe quadratique: (P0) P(U) = (U2, U. DM(Pl) (P2) avec M une matrice caractéristique de l'interpolation quadratique: 2-4 2 M = -3 4 -1 1 0 1 et u le paramètre de l'interpolation variant entre 0 et 1. Pour l'interpolation biquadratique, une trajectoire d'un point dans le plan est paramétrce par deux paramètres d'interpolation u et v, et s'exprime: (Po O Po,1 Po,2) (u2) P(u, v) = (v2, V, l)M(Pl,oPl,l P192)M (U) (P2,0P2,l P2,2) (1) La généralisation à un nombre de points plus élevé est possible. La fonction d'interpolation relie les points de contrôle à tout point de la trajectoire par les paramètres d'interpolation u et v. Plus précisément nous relions ces points à des angles anatomique s définis en fonction d' angles de la machine. On peut prendre u et v comme coordonnces réduites des angles anatomiques, du genre u = (0-(3min)/(0max- min) Toutefois, on préférera prendre u = f((3-(3min)/(êmax-(3min)) et v = g((±min)/(+max-4)min)) f et g sont des fonctions impaires, croissantes et telle que f( 1) = 1 et f'(l) = 0, ce qui aboutit à une fonction polynomiale de type ax5 + bx3 + cx. Plus précisément, on parvient à utiliser une fonction incluant une constante A qui dépend du type d'appareil de radiologie et qui sera conservoe pour une même série d'appareils. La fonction de distorsion s'écrit alors: type f(x) = 0c(A+,B)xs + X(A+)x3 + ú(A+Ox. Une bonne approximation est fournie par la fonction suivante: f(x) = (A 1,5)x5 + (2,5 - 2A) x3 + Ax, avec A compris entre 0 et 15/8. La valeur de A = 1,25 est particulièrement intéressante dans la mesure o elle permet d'obtenir une estimation convenable de la distorsion. A titre d'exemple, pour un champ de 16, l'erreur moyenne est d'environ 0,8 pixcl avec 1,95 pixcls d'erreur maximale. Pour un champ de 22, l'erreur moyenne est d'environ 0,75 pixcl avec une erreur maximale de 1,65 pixcls. Pour un champ de 30, l'erreur moyenne est de l'ordre de 0,6 pixcl avec une erreur maximale de 1,42 pixcls. Plus généralement, on définit N points de contrôle, d'o N orientations (An)l n) coordonnces du point de contrôle numéro n. N configurations de la machine sont définies: - le balayage caméra est le même que pour la configuration de référence (par exemple balayage normal) - la matrice d'acquisition est la même que pour la configuration de référence (par exemple matrice 512) la taille du champ de vue est la même que pour la configuration de référence (par exemple champ de 22 cm) - l'orientation est celle correspondant au point de contrôle: An = ( nq>n) On acquiert une image de grille pour ces N configurations et on en daduit N champs de distorsions différentiels (DAn) l (ou 16) champs de distorsion différentiels. Un champ de distorsion différentiel est un ensemble de vecteurs bidimensionnels donnant la distorsion pour tous les noeuds de la grille. Par interpolation, connaissant l'orientation courante (é,p), on retrouve pour chaque noeud de la grille la valeur de la distorsion différentielle d. Ainsi, on détermine le champ de distorsion différentiel DA pour l'orientation courante ((3,O. Sur la figure 3, on voit les différentes allures que prend la fonction f(x, A) en faisant varier A entre 0 et 1,875. Pour parfaire l'estimation de la distorsion, on peut considérer une fonction d'interpolation différente selon chaque axe, par exemple en adaptant la constante A. La constante A peut faire l'objet d'une estimation automatique en fonction des donnces d'un étalonnage composces de neuf images. Dans un mode de réalisation particulier de l'invention, on cherche à estimer deux champs de distorsion, le champ DO lié au champ électromagnétique terrestre à l'endroit o se trouve l'appareil de radiologie et le champ DV lié aux différentes angulations de l'appareil de radiologie. On pourra ainsi effectuer un étalonnage hebdomadaire ou mensuel du champ de distorsion Do en raison de l'évolution relativement 3 0 rapi de du champ électromagnétique terre stre et n' effectuer un étal onn age du champ de distorsion DV qu'à intervalles plus larges, par exemple semestriel ou annuel en raison de la dérive plus lente de cette distorsion. On peut égal ement c alculer un ch amp de di stors ion d' inversion DI qui est lié à une inversion de l'image par inversion du balayage vidéo du détecteur. En effet, le détecteur n'étant pas parfait, l'inversion du balayage crce aussi une distorsion qu'il est utile d'appréhender, notamment dans des applications de neuroradiologie, en particulier avec des machines de radiologie biplan équipées de deux tubes et de deux détecteurs pour lesquels une inversion mécanique de la prise de vue, en passant, par exemple, le tube de la droite à la gauche et le détecteur de la gauche à la droite, n'est pas possible. L'estimation du champ de distorsion DI est effectuce en comparant une image obtenue en balayage normal puis retournce et une image obtenue en balayage inversé de la même grille. En outre, on peut améliorer l'estimation de la distorsion en tenant compte du changement de champ de l'image. On estime alors un champ de distorsion DR lié au réglage de résolution du détecteur. Ainsi, un étalonnage complet est mis en _uvre lors de la réception et de la mise en service de l'appareil. Puis un étalonnage partiel est mis en _uvre à intervalles temporels courts et un étalonnage complet est mis en _uvre à intervalles temporels longs. L'étalonnage partiel pourra consister à seulement recalculer le champ de distorsion DO lié au champ électromagnétique terrestre. L'étalonnage complet pourra comprendre les étapes suivantes, voir figure 4: étape 16: installation de la grille sur le trajet du faisccau d'énergie par un responsable de maintenance, étape 17: prise de neuf images de la grille pour les points de contrôle: (0,0), (-(30, 0), (O0, 0) , (0, -q)0), (0, Oo)' (-0o' -q)o), (30' to)' (00, 0), et (30, +0). étape 18: calcul des champs de distorsion correspondants aux points de contrôle, et étape 19: calcul des champs de distorsion Do et Dv I1 convient de noter que les étapes 17 et 18 peuvent être effectuées partiellement en même temps pour réduire la durce de l'étalonnage complet. Les autres champs de distorsion DI et DR feront, le cas échéant, l 'obj et d'étapes spécifiques, telles que pour le champ DI prise d'une image de la grille pour le point de contrôle (0,0) en balayage inversé, retournement de l'image de la grille pour le point de contrôle (0,0) en balayage normal et comparaison, et pour le champ DR prise d'images de la grille pour le point de contrôle (0,0) sous différentes résolutions et comparaison. L'étalonnage partiel pourra comprendre les étapes suivantes: a) installation de la grille sur le trajet du faisccau d'énergie par un responsable de maintenance, b) prise d'une image de la grille pour le point de contrôle (0,0). c) calcul du champ de distorsion DO correspondant. On peut ensuite effectuer une combinaison du champ de distorsion DO et des champs de distorsion Dv, DI, DR' et DF non modifiés lors de l'étalonnage partiel pour obtenir les champs de distorsion interpolés pour toute angulation (D, O lorsque l'on désire corriger une image de la distorsion qu'elle contient. Reconsidérons à présent le problème initial qui consiste à estimer le champ de distorsion D pour une configuration de la machine donnce. On connâît le champ de distorsion différentiel DI correspondant au changement de balayage caméra entre la configuration de référence et la configuration étudice ("configuration courante") On connâît le champ de distorsion différentiel DM correspondant au changement de matrice d'acquisition entre la configuration de référence et la configuration courante On connâît le champ de distorsion différentiel DF correspondant au changement de taille du champ de vue entre la configuration de référence et la configuration courante 25. On connâît le champ de distorsion différentiel DA correspondant au changement d'orientation de l'arccau entre la configuration de référence et la configuration courante (par interpolation). Ainsi, on peut en déduire le champ de distorsion pour la configuration courante de la machine: D=DAoDFoDMoDIoD0 Grâce à l'invention, on dispose d'un étalonnage économique, avec des durces d'immobilisation de l'appareil d'imagerie en vue de la maintenance réduite, et de haute précision dont une partie est réalisce lors de l'installation de l'appareil d'imagerie et n'a plus à être mise en oeuvre par la suite et qui est particulièrement adaptée pour la mise en correspondence d'images bidimensionnelles et d'images tridimensionnelles, qui permet les reconstructions tridimensionnelles à partir d'un nombre limité de vues. On pout supprimer l'utilisation de marqueurs externes, en radiologie, au profit d'un étalonnage effectué avant l'arrivoe du patient, ce qui réduit la durce de l'examen radiologique effectué sur le patient. On peut effectuer des mesures dans les images bidimensionnelles avec une précision accrue. REVENDICATIONS 1. Procédé d'étalonnage d'un système de d'imagerie, comprenant un moyen d'émission d'un faisccau d'énergie, un moyen de réception du faisccau d'énergie après qu'il a traversé un organe à étudier, un moyen de commande du moyen de réception, un moyen de support dudit organe et un moyen de mise en forme d'images, dans lequel on étalonne des paramètres intrinsèques relatifs aux différents types d'images pouvant être prises dans une configuration prédéterminée du système, les paramètres intrinsèques étant estimés à l'aide d'un fantôme bidimensionnel pour établir un champ de distorsion pour une configuration particulière, le champ de distorsion pour une configuration quelconque étant calculé par combinaison d'une pluralité de champs de distorsion de causes supposces indépendantes. 2. Procédé selon la revendication 1, dans lequel dans lequel on détermine un champ de distorsion Do dû au champ électromagnétique terrestre à partir d'une image dudit fantôme prise dans une position centrale du système d'imagerie selon chaque axe de rotation, ledit champ de distorsion Do étant indépendant de l'angulation. 3. Procédé selon la revendication 1 ou 2, dans lequel on étalonne des paramètres intrinsèques relatifs aux différents types d'images pouvant être prises dans des positions angulaires prédéterminces du moyen d'émission et du moyen de réception, les paramètres intrinsèques étant estimés à l'aide d'un fantôme bidimensionnel pour établir un champ de distorsion pour une angulation particulière, le champ de distorsion pour une angulation quelconque étant calculé par interpolation polynomiale des champs de distorsion établis chacun pour un point de contrôle, un point de contrôle étant défini par une angulation particulière par axe de rotation du système d'imagerie, devant faire l'objet d'un étalonnage. 4. Procédé selon la revendication 3, dans lequel un point de contrôle est proche de chaque extremum de position angulaire. 5. Procédé selon la revendication 3 ou 4, dans lequel la fonction d'interpolation est de type biquadratique ou bicubique. 6. Procédé selon l'une quelconque des revendications précédentes, dans lequel on établit une relation entre des paramètres d'interpolation u, v et l'angulation, la relation étant de type u(x) = ax5+ bx3 + cx, x étant une coordonnce d'un point, de façon que le champ de distorsion soit déterminé en fonction de l'angulation, la fonction d'interpolation permettant de déterminer la distorsion en fonction des paramètres d'interpolation u, v. 7. Procédé selon la revendication 6, dans lequel la fonction u(x) est de type u(x) = cc(A+,13)x5 + %(A+o)x3 + e(A+q,)x, A étant une constante dépendant du type de système d'imagerie, oc,,B, %, 8, et étant des constantes déterminces expérimentalement. 8. Procédé selon l'une quelconque des revendications précédentes, dans lequel on détermine un champ de distorsion DI dû à une inversion d'image par inversion du sens de balayage du moyen de réception du faisccau d'énergie par comparaison d'une image dudit fantôme retournce après sa prise et d'une image dudit fantôme prise en balayage inversé, et ce pour les changements de balayage utilisés. 9. Procédé selon l'une quelconque des revendications précédentes, dans lequel on détermine un champ de distorsion DR dû à un changement de résolution du moyen de réception du faisccau d'énergie par comparaison d'une image dudit fantôme prise avec une première résolution et d'une image dudit fantôme prise avec une deuxième résolution différente de la première, et ce pour les changements de résolution utilisés. 10. Procédé selon l'une quelconque des revendications précédentes, dans lequel on détermine un champ de distorsion DM dû à un changement de matrice d'acquisition du moyen de réception du faisccau d'énergie par comparaison d'une image dudit fantôme prise avec une première matrice d'acquisition et d'une image dudit fantôme prise avec une deuxième matrice d'acquisition différente de la première, et ce pour les changements de matrice d'acquisition utilisés. 11. Support de mémoire contenant des instructions pour mettre en _uvre un procédé d'étalonnage d'un système d'imagerie comprenant un moyen d'émission d'un faisccau d'énergie, un moyen de réception du faisccau d'énergie après qu'il a traversé un organe à étudier, un moyen de commande du moyen de réception, un moyen de support dudit organe et un moyen de mise en forme d'images, procédé dans lequel on étalonne des paramètres intrinsèques relatifs aux différents types d'images pouvant être prises dans une configuration prédétermince du système, les paramètres intrinsèques étant estimés à l'aide d'un fantôme bidimensionnel pour établir un champ de distorsion pour une configuration particulière, le champ de distorsion pour une configuration quelconque étant calculé par combinaison d'une pluralité de champs de distorsion de causes supposées indépendantes. 12. Système d'imagerie, comprenant un moyen d'émission (1 1) d'un faisccau d'énergie, un moyen de réception (12) du faisccau d'énergie après qu'il a traversé un organe à étudier, un moyen de commande du moyen de réception, et un moyen de mise en forme d'images, caractérisé par le fait qu 'i l comprend un moyen d'étal onn age de s paramètres intrinsèques relatifs aux différents types d'images pouvant être prises dans une configuration prédétermince du système, les paramètres intrinsèques étant estimés à l'aide d'un fantôme bidimensionnel pour établir un champ de distorsion pour une configuration particulière, le champ de distorsion pour une configuration quelconque étant calculé par combinaison d'une pluralité de champs de distorsion de causes supposces