3D registration of MR and X-ray spine images using an articulated model.

Présentation: Cet article a été publié dans le journal : Computerised medical imaging and graphics (CMIG). Le but de cet article est de recaler les vertèbres extraites à partir d’images RM avec des vertèbres extraites à partir d’images RX pour des patients scoliotiques, en tenant compte des déformations non-rigides due au changement de posture entre ces deux modalités. À ces fins, une méthode de recalage à l’aide d’un modèle articulé est proposée. Cette méthode a été comparée avec un recalage rigide en calculant l’erreur sur des points de repère, ainsi qu’en calculant la différence entre l’angle de Cobb avant et après recalage. Une validation additionelle de la méthode de recalage présentée ici se trouve dans l’annexe A. Ce travail servira de première étape dans la fusion des images RM, RX et TP du tronc complet. Donc, cet article vérifie l’hypothèse 1 décrite dans la section 3.2.1.

Recherche d’étoiles jeunes de faible masse dans le voisinage solaire

Formées lors de l’effondrement gravitationnel d’un nuage de gaz moléculaire, les étoiles naissantes auront différentes masses variant entre 0.08 et environ 100M . La majorité de la population stellaire de la Galaxie est constituée d’étoiles dont la masse est inférieure à environ 0.6 M . Le dernier évènement de formation stellaire dans le voisinage solaire s’est produit dans la bulle locale il y a au plus 100 millions d’années, vraisemblablement provoqué par le passage d’une onde de choc dans le bras local de la Galaxie. C’est ainsi que se formèrent de jeunes associations d’étoiles dont les membres se caractérisent en particulier par une vitesse spatiale et une position commune dans la Galaxie. Les associations jeunes étant peu densément peuplées et relativement proches du Soleil, leurs membres se font plutôt rares et dispersés sur toute la voûte céleste. Jusqu’à présent, surtout les étoiles les plus massives (brillantes) ont été répertoriées. Les étoiles jeunes de faible masse, constituant la majorité de la population, restent pour la plupart à être identifiées. Les étoiles jeunes de faible masse représentent une population clef pour contraindre les modèles évolutifs des étoiles M et des naines brunes. Elles sont également d’excellentes candidates pour chercher des exoplanètes via les techniques d’imagerie directe. Ce mémoire présente une nouvelle méthode utilisant un modèle cinématique enrichi d’une analyse statistique Bayesienne pour identifier des étoiles jeunes de faible masse dans les associations beta Pictoris, Tucana-Horologium et AB Doradus. À partir d’un échantillon de 1080 étoiles K et M, toutes comportant des indicateurs de jeunesse tels l’émission Halpha et une forte luminosité dans les rayons X, leurs propriétés cinématiques (mouvement propre) et photométriques sont analysées pour en extraire 98 candidates hautement probables membres d’une des trois associations. Une confirmation de leur statut comme membre nécessitera en particulier une mesure de leur vitesse radiale (prédit par notre analyse) et une mesure de la largeur équivalente du lithium à 6708 Å pour mieux contraindre leur âge.

Structure d'une tagatose-1,6-bisphosphate aldolase de classe I : étude d'une apparente perte de stéréospécificité

La tagatose-1,6-biphosphate aldolase de Streptococcus pyogenes est une aldolase de classe I qui fait montre d'un remarquable manque de spécificité vis à vis de ses substrats. En effet, elle catalyse le clivage réversible du tagatose-1,6-biphosphate (TBP), mais également du fructose-1,6-biphosphate (FBP), du sorbose-1,6-biphosphate et du psicose-1,6-biphosphate, quatre stéréoisomères, en dihydroxyacétone phosphate (DHAP) et en glycéraldéhyde-3-phosphate (G3P). Afin de mettre à jour les caractéristiques du mécanisme enzymatique, une étude structurale de la TBP aldolase de S. pyogenes, un pathogène humain extrêmement versatile, a été entreprise. Elle a permis la résolution de la structure native et en complexe avec le DHAP, a respectivement 1.87 et 1.92 Å de résolution. Ces mêmes structures ont permis de se représenter plus clairement le site actif de l'enzyme en général, et les résidus catalytiques en particulier. Le trempage des cristaux de TBP aldolase dans une solution saturante de DHAP a en outre permis de piéger un authentique intermédiaire iminium, ainsi que sa géométrie particulière en atteste. Des expériences d'échange de proton, entreprises afin d'évaluer le stéréoisomérisme du transfert de proton catalytique, ont également permis de faire une intéressante découverte : la TBP aldolase ne peut déprotoner le coté pro-R du C3 du DHAP, mais peut le protonner. Ce résultat, ainsi que la comparaison de la structure du complexe TBP aldolase-DHAP avec la structure du complexe FBP aldolase de muscle de lapin- DHAP, pointe vers un isomérisme cis-trans autour du lien C2-C3 de la base de Schiff formée avec le DHAP. De plus, la résolution de ces deux structures a permis de mettre en évidence trois régions très mobiles de la protéine, ce qui pourrait être relié au rôle postulé de son isozyme chez S. pyogenes dans la régulation de l’expression génétique et de la virulence de la bactérie. La cristallographie par rayons X et la cinétique enzymatique ont ainsi permis d'avancer dans l'élucidation du mécanisme et des propriétés structurales de cette enzyme aux caractéristiques particulières.

Développement et validation d’un outil clinique pour l’analyse quantitative de la posture auprès de personnes atteintes d’une scoliose idiopathique

La scoliose idiopathique (SI) est une déformation tridimensionnelle (3D) de la colonne vertébrale et de la cage thoracique à potentiel évolutif pendant la croissance. Cette déformation 3D entraîne des asymétries de la posture. La correction de la posture est un des objectifs du traitement en physiothérapie chez les jeunes atteints d’une SI afin d’éviter la progression de la scoliose, de réduire les déformations morphologiques et leurs impacts sur la qualité de vie. Les outils cliniques actuels ne permettent pas de quantifier globalement les changements de la posture attribuables à la progression de la scoliose ou à l’efficacité des interventions thérapeutiques. L’objectif de cette thèse consiste donc au développement et à la validation d’un nouvel outil clinique permettant l’analyse quantitative de la posture auprès de personnes atteintes d’une SI. Ce projet vise plus spécifiquement à déterminer la fidélité et la validité des indices de posture (IP) de ce nouvel outil clinique et à vérifier leur capacité à détecter des changements entre les positions debout et assise. Suite à une recension de la littérature, 34 IP représentant l’alignement frontal et sagittal des différents segments corporels ont été sélectionnés. L’outil quantitatif clinique d’évaluation de la posture (outil 2D) construit dans ce projet consiste en un logiciel qui permet de calculer les différents IP (mesures angulaires et linéaires). L’interface graphique de cet outil est conviviale et permet de sélectionner interactivement des marqueurs sur les photographies digitales. Afin de vérifier la fidélité et la validité des IP de cet outil, la posture debout de 70 participants âgés entre 10 et 20 ans atteints d'une SI (angle de Cobb: 15º à 60º) a été évaluée à deux occasions par deux physiothérapeutes. Des marqueurs placés sur plusieurs repères anatomiques, ainsi que des points de référence anatomique (yeux, lobes des oreilles, etc.), ont permis de mesurer les IP 2D en utilisant des photographies. Ces mêmes marqueurs et points de référence ont également servi au calcul d’IP 3D obtenus par des reconstructions du tronc avec un système de topographie de surface. Les angles de Cobb frontaux et sagittaux et le déjettement C7-S1 ont été mesurés sur des radiographies. La théorie de la généralisabilité a été utilisée pour déterminer la fidélité et l’erreur standard de la mesure (ESM) des IP de l’outil 2D. Des coefficients de Pearson ont servi à déterminer la validité concomitante des IP du tronc de l’outil 2D avec les IP 3D et les mesures radiographiques correspondantes. Cinquante participants ont été également évalués en position assise « membres inférieurs allongés » pour l’étude comparative de la posture debout et assise. Des tests de t pour échantillons appariés ont été utilisés pour détecter les différences entre les positions debout et assise. Nos résultats indiquent un bon niveau de fidélité pour la majorité des IP de l’outil 2D. La corrélation entre les IP 2D et 3D est bonne pour les épaules, les omoplates, le déjettement C7-S1, les angles de taille, la scoliose thoracique et le bassin. Elle est faible à modérée pour la cyphose thoracique, la lordose lombaire et la scoliose thoraco-lombaire ou lombaire. La corrélation entre les IP 2D et les mesures radiographiques est bonne pour le déjettement C7-S1, la scoliose et la cyphose thoracique. L’outil est suffisamment discriminant pour détecter des différences entre la posture debout et assise pour dix des treize IP. Certaines recommandations spécifiques résultents de ce projet : la hauteur de la caméra devrait être ajustée en fonction de la taille des personnes; la formation des juges est importante pour maximiser la précision de la pose des marqueurs; et des marqueurs montés sur des tiges devraient faciliter l’évaluation des courbures vertébrales sagittales. En conclusion, l’outil développé dans le cadre de cette thèse possède de bonnes propriétés psychométriques et permet une évaluation globale de la posture. Cet outil devrait contribuer à l’amélioration de la pratique clinique en facilitant l’analyse de la posture debout et assise. Cet outil s’avère une alternative clinique pour suivre l’évolution de la scoliose thoracique et diminuer la fréquence des radiographies au cours du suivi de jeunes atteints d’une SI thoracique. Cet outil pourrait aussi être utile pour vérifier l’efficacité des interventions thérapeutiques sur la posture.

Développement de tests mesurant les habiletés de perception et d'interprétation des radiographies pulmonaires

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal

Morphologie tridimensionnelle du bassin scoliotique

La scoliose idiopathique de l’adolescence (SIA) est une déformation tridimensionnelle de la colonne vertébrale et de la cage thoracique dont la cause est inconnue. Il semble que la ceinture pelvienne soit impliquée dans la pathogénie de la SIA, car des différences géométriques des os coxaux ont été observées. Notamment, une rotation du bassin ou une inclinaison latérale dans le sens de la courbe scoliotique ont été mises en évidence en plus des distorsions osseuses. Il est difficile de dissocier la rotation du bassin de son asymétrie, car la majorité des études porte sur des données radiologiques bidimensionnelles. Une analyse tridimensionnelle de la morphologie du bassin de patientes ayant une SIA, mais n’ayant pas reçu de traitement par corset ou chirurgie permettrait d’identifier le rôle du bassin dans la progression de la scoliose. Dix-sept jeunes filles atteintes de la SIA ont participé à cette étude pour lesquelles des radiographies bi-planaires en station debout étaient disponibles au moment du diagnostic par un chirurgien orthopédiste pédiatrique et à l’instant de la prescription d'un corset. Des radiographies postéro-antérieures et latérales avaient été prises au moyen du système EOS®. Douze repères anatomiques du bassin ont été identifiés sur les paires de radiographies, alors que quatre repères ont été identifiés sur la radiographie postéro-antérieure uniquement. Ces quatre derniers n’étaient pas identifiables sur la radiographie latérale à cause de la superposition des repères droits et gauches. La reconstruction tridimensionnelle du bassin a été réalisée à partir de deux clichés radiographiques des 12 premiers repères osseux. Au total, neuf paramètres tridimensionnels ont été calculés afin de quantifier l’asymétrie et la distorsion du bassin entre les deux temps donnés. Des paramètres bidimensionnels ont également été mesurés sur les quatre derniers repères osseux afin de documenter des déformations du bassin pertinentes à la pratique clinique, telle que la rotation axiale de celui-ci. Dans le but d'évaluer une possible asymétrie entre les os coxaux du bassin, les paramètres tridimensionnels du bassin gauche ont été comparés à ceux du côté droit à chaque temps, au moyen d'un test-t pour échantillon apparié. La morphologie pelvienne a été également évaluée par l'analyse multivariée (MANOVA) à mesures répétées à deux conditions (côté*temps). En conséquence, nous avons constaté une croissance osseuse statistiquement significative du bassin dans l’intervalle de temps entre le diagnostic de la scoliose et le port du corset (p=0,033). Une asymétrie significative entre les côtés gauche et droit du bassin (p=0,013) a également été constatée. En ce qui concerne les paramètres bidimensionnels, nous avons constaté une augmentation de la version pelvienne (p=0,024) au cours de la croissance des jeunes filles. Finalement, le bassin n'a pas présenté de distorsion, puisqu'une valeur de p de 0,763 a été observée. En conclusion, la croissance des jeunes filles atteintes de la scoliose idiopathique de l'adolescence est accompagnée d'une asymétrie morphologique entre les deux os coxaux du bassin. Cette asymétrie constatée au moment du diagnostic de la scoliose des filles a évolué jusqu'à l’instant où le port du corset a été prescrit. Quant aux paramètres bidimensionnels, nous pouvons conclure que la rotation du bassin vers l'arrière a augmenté au cours de la croissance des jeunes filles, produisant ainsi une rétroversion pelvienne dans le plan sagittal. La distorsion tridimensionnelle du bassin n'a toutefois pas été observée au cours de la croissance des jeunes filles.

Patient-specific model of a scoliotic torso for surgical planning

A method for the construction of a patient-specific model of a scoliotic torso for surgical planning via inter- patient registration is presented. Magnetic Resonance Images (MRI) of a generic model are registered to surface topography (TP) and X-ray data of a test patient. A partial model is first obtained via thin-plate spline registration between TP and X-ray data of the test patient. The MRIs from the generic model are then fit into the test patient using articulated model registration between the vertebrae of the generic model’s MRIs in prone position and the test patient’s X-rays in standing position. A non-rigid deformation of the soft tissues is performed using a modified thin-plate spline constrained to maintain bone rigidity and to fit in the space between the vertebrae and the surface of the torso. Results show average Dice values of 0.975 ± 0.012 between the MRIs following inter-patient registration and the surface topography of the test patient, which is comparable to the average value of 0.976 ± 0.009 previously obtained following intra-patient registration. The results also show a significant improvement compared to rigid inter-patient registration. Future work includes validating the method on a larger cohort of patients and incorporating soft tissue stiffness constraints. The method developed can be used to obtain a geometric model of a patient including bone structures, soft tissues and the surface of the torso which can be incorporated in a surgical simulator in order to better predict the outcome of scoliosis surgery, even if MRI data cannot be acquired for the patient.

Multimodal image registration of the scoliotic torso for surgical planning

Background This paper presents a method that registers MRIs acquired in prone position, with surface topography (TP) and X-ray reconstructions acquired in standing position, in order to obtain a 3D representation of a human torso incorporating the external surface, bone structures, and soft tissues. Methods TP and X-ray data are registered using landmarks. Bone structures are used to register each MRI slice using an articulated model, and the soft tissue is confined to the volume delimited by the trunk and bone surfaces using a constrained thin-plate spline. Results The method is tested on 3 pre-surgical patients with scoliosis and shows a significant improvement, qualitatively and using the Dice similarity coefficient, in fitting the MRI into the standing patient model when compared to rigid and articulated model registration. The determinant of the Jacobian of the registration deformation shows higher variations in the deformation in areas closer to the surface of the torso. Conclusions The novel, resulting 3D full torso model can provide a more complete representation of patient geometry to be incorporated in surgical simulators under development that aim at predicting the effect of scoliosis surgery on the external appearance of the patient’s torso.

La tagatose-1,6-bisphosphate aldolase et la fructose-1,6-bisphosphate aldolase de classe I : mécanisme et stéréospécificité

La tagatose-1,6-biphosphate aldolase de Streptococcus pyogenes est une aldolase qui fait preuve d'un remarquable manque de spécificité vis à vis de ses substrats. En effet, elle catalyse le clivage réversible du tagatose-1,6-bisphosphate (TBP), mais également du fructose-1,6-bisphosphate (FBP), du sorbose-1,6-bisphosphate et du psicose-1,6-bisphosphate, quatre stéréoisomères, en dihydroxyacétone phosphate (DHAP) et en glycéraldéhyde-3-phosphate (G3P). Aldolase de classe I, qui donc catalyse sa réaction en formant un intermédiaire covalent obligatoire, ou base de Schiff, avec son susbtrat, la TBP aldolase de S. pyogenes partage 14 % d’identité avec l’enzyme modèle de cette famille, la FBP aldolase de muscle de mammifère. Bien que le mécanime catalytique de la FBP aldolase des mammifères ait été examiné en détails et qu’il soit approprié d’en tirer des renseignements quant à celui de la TBP aldolase, le manque singulier de stéréospécificité de cette dernière tant dans le sens du clivage que celui de la condensation n’est toujours pas éclairci. Afin de mettre à jour les caractéristiques du mécanisme enzymatique, une étude structurale de la TBP aldolase de S. pyogenes, un pathogène humain extrêmement versatile, a été entreprise. Elle a permis la résolution des structures de l’enzyme native et mutée, en complexe avec des subtrats et des inhibiteurs compétitifs, à des résolutions comprises entre 1.8 Å et 2.5 Å. Le trempage des cristaux de TBP aldolase native et mutante dans une solution saturante de FBP ou TBP a en outre permis de piéger un authentique intermédiaire covalent lié à la Lys205, la lysine catalytique. La determination des profils pH de la TBP aldolase native et mutée, entreprise afin d'évaluer l’influence du pH sur la réaction de clivage du FBP et TBP et ìdentifier le(s) résidu(s) impliqué(s), en conjonction avec les données structurales apportées par la cristallographie, ont permis d’identifier sans équivoque Glu163 comme résidu responsable du clivage. En effet, le mode de liaison sensiblement différent des ligands utilisés selon la stéréochimie en leur C3 et C4 permet à Glu163, équivalent à Glu187 dans la FBP aldolase de classe I, d’abstraire le proton sur l’hydroxyle du C4 et ainsi d’amorcer le clivage du lien C3-C4. L’étude du mécanimse inverse, celui de la condensation, grâce par exemple à la structure de l’enzyme native en complexe avec ses substrats à trois carbones le DHAP et le G3P, a en outre permis d’identifier un isomérisme du substrat G3P comme possible cause de la synthèse des isomères en C4 par cette enzyme. Ce résultat, ainsi que la decouverte d’un possible isomérisme cis-trans autour du lien C2-C3 de la base de Schiff formée avec le DHAP, identifié précedemment, permet de cerner presque complètement les particularités du mécanisme de cette enzyme et d’expliquer comment elle est capable de synthétiser les quatres stéréoisomères 3(S/R), 4(S/R). De plus, la résolution de ces structures a permis de mettre en évidence trois régions très mobiles de la protéine, ce qui pourrait être relié au rôle postulé de son isozyme chez S. pyogenes dans la régulation de l’expression génétique et de la virulence de la bactérie. Enfin, la résolution de la structure du mutant Lys229→Met de la FBP aldolase de muscle en complexe avec la forme cyclique du FBP, de même que des études cristallographiques sur le mutant équivalent Lys205→Met de la TBP aldolase de S. pyogenes et des expériences de calorimétrie ont permis d’identifier deux résidus particuliers, Ala31 et Asp33 chez la FBP aldolase, comme possible cause de la discrimination de cette enzyme contre les substrats 3(R) et 4(S), et ce par encombrement stérique des substrats cycliques. La cristallographie par rayons X et la cinétique enzymatique ont ainsi permis d'avancer dans l'élucidation du mécanisme et des propriétés structurales de cette enzyme aux caractéristiques particulières.

Une nouvelle méthode pour estimer la torsion géométrique en scoliose idiopathique de l’adolescent

La scoliose idiopathique de l’adolescent (SIA) est une déformation tridimensionnelle (3D) de la colonne vertébrale. Pour la plupart des patients atteints de SIA, aucun traitement chirurgical n’est nécessaire. Lorsque la déformation devient sévère, un traitement chirurgical visant à réduire la déformation est recommandé. Pour déterminer la sévérité de la SIA, l’imagerie la plus utilisée est une radiographie postéroantérieure (PA) ou antéro-postérieure (AP) du rachis. Plusieurs indices sont disponibles à partir de cette modalité d’imagerie afin de quantifier la déformation de la SIA, dont l’angle de Cobb. La conduite thérapeutique est généralement basée sur cet indice. Cependant, les indices disponibles à cette modalité d’imagerie sont de nature bidimensionnelle (2D). Celles-ci ne décrivent donc pas entièrement la déformation dans la SIA dû à sa nature tridimensionnelle (3D). Conséquemment, les classifications basées sur les indices 2D souffrent des mêmes limitations. Dans le but décrire la SIA en 3D, la torsion géométrique a été étudiée et proposée par Poncet et al. Celle-ci mesure la tendance d’une courbe tridimensionnelle à changer de direction. Cependant, la méthode proposée est susceptible aux erreurs de reconstructions 3D et elle est calculée localement au niveau vertébral. L’objectif de cette étude est d’évaluer une nouvelle méthode d’estimation de la torsion géométrique par l’approximation de longueurs d’arcs locaux et par paramétrisation de courbes dans la SIA. Une première étude visera à étudier la sensibilité de la nouvelle méthode présentée face aux erreurs de reconstructions 3D du rachis. Par la suite, deux études cliniques vont présenter la iv torsion géométrique comme indice global et viseront à démontrer l’existence de sous-groupes non-identifiés dans les classifications actuelles et que ceux-ci ont une pertinence clinique. La première étude a évalué la robustesse de la nouvelle méthode d’estimation de la torsion géométrique chez un groupe de patient atteint de la SIA. Elle a démontré que la nouvelle technique est robuste face aux erreurs de reconstructions 3D du rachis. La deuxième étude a évalué la torsion géométrique utilisant cette nouvelle méthode dans une cohorte de patient avec des déformations de type Lenke 1. Elle a démontré qu’il existe deux sous-groupes, une avec des valeurs de torsion élevées et l’autre avec des valeurs basses. Ces deux sous-groupes possèdent des différences statistiquement significatives, notamment au niveau du rachis lombaire avec le groupe de torsion élevée ayant des valeurs d’orientation des plans de déformation maximales (PMC) en thoraco-lombaire (TLL) plus élevées. La dernière étude a évalué les résultats chirurgicaux de patients ayant une déformation Lenke 1 sous-classifiées selon les valeurs de torsion préalablement. Cette étude a pu démontrer des différences au niveau du PMC au niveau thoraco-lombaire avec des valeurs plus élevées en postopératoire chez les patients ayant une haute torsion. Ces études présentent une nouvelle méthode d’estimation de la torsion géométrique et présentent cet indice quantitativement. Elles ont démontré l’existence de sous-groupes 3D basés sur cet indice ayant une pertinence clinique dans la SIA, qui n’étaient pas identifiés auparavant. Ce projet contribue dans la tendance actuelle vers le développement d’indices 3D et de classifications 3D pour la scoliose idiopathique de l’adolescent.

Scoliosis curve type classification from 3D trunk image

Adolescent idiopathic scoliosis (AIS) is a deformity of the spine manifested by asymmetry and deformities of the external surface of the trunk. Classification of scoliosis deformities according to curve type is used to plan management of scoliosis patients. Currently, scoliosis curve type is determined based on X-ray exam. However, cumulative exposure to X-rays radiation significantly increases the risk for certain cancer. In this paper, we propose a robust system that can classify the scoliosis curve type from non invasive acquisition of 3D trunk surface of the patients. The 3D image of the trunk is divided into patches and local geometric descriptors characterizing the surface of the back are computed from each patch and forming the features. We perform the reduction of the dimensionality by using Principal Component Analysis and 53 components were retained. In this work a multi-class classifier is built with Least-squares support vector machine (LS-SVM) which is a kernel classifier. For this study, a new kernel was designed in order to achieve a robust classifier in comparison with polynomial and Gaussian kernel. The proposed system was validated using data of 103 patients with different scoliosis curve types diagnosed and classified by an orthopedic surgeon from the X-ray images. The average rate of successful classification was 93.3% with a better rate of prediction for the major thoracic and lumbar/thoracolumbar types.

Scoliosis Follow-Up Using Noninvasive Trunk Surface Acquisition

Adolescent idiopathic scoliosis (AIS) is a musculoskeletal pathology. It is a complex spinal curvature in a 3-D space that also affects the appearance of the trunk. The clinical follow-up of AIS is decisive for its management. Currently, the Cobb angle, which is measured from full spine radiography, is the most common indicator of the scoliosis progression. However, cumulative exposure to X-rays radiation increases the risk for certain cancers. Thus, a noninvasive method for the identification of the scoliosis progression from trunk shape analysis would be helpful. In this study, a statistical model is built from a set of healthy subjects using independent component analysis and genetic algorithm. Based on this model, a representation of each scoliotic trunk from a set of AIS patients is computed and the difference between two successive acquisitions is used to determine if the scoliosis has progressed or not. This study was conducted on 58 subjects comprising 28 healthy subjects and 30 AIS patients who had trunk surface acquisitions in upright standing posture. The model detects 93% of the progressive cases and 80% of the nonprogressive cases. Thus, the rate of false negatives, representing the proportion of undetected progressions, is very low, only 7%. This study shows that it is possible to perform a scoliotic patient's follow-up using 3-D trunk image analysis, which is based on a noninvasive acquisition technique.

Personalized 3D reconstruction of the rib cage for clinical assessment of trunk deformities

Scoliosis is a 3D deformity of the spine and rib cage. Extensive validation of 3D reconstruction methods of the spine from biplanar radiography has already been published. In this article, we propose a novel method to reconstruct the rib cage, using the same biplanar views as for the 3D reconstruction of the spine, to allow clinical assessment of whole trunk deformities. This technique uses a semi-automatic segmentation of the ribs in the postero-anterior X-ray view and an interactive segmentation of partial rib edges in the lateral view. The rib midlines are automatically extracted in 2D and reconstructed in 3D using the epipolar geometry. For the ribs not visible in the lateral view, the method predicts their 3D shape. The accuracy of the proposed method has been assessed using data obtained from a synthetic bone model as a gold standard and has also been evaluated using data of real patients with scoliotic deformities. Results show that the reconstructed ribs enable a reliable evaluation of the rib axial rotation, which will allow a 3D clinical assessment of the spine and rib cage deformities.

Semiautomatic Detection of Scoliotic Rib Borders From Posteroanterior Chest Radiographs

3-D assessment of scoliotic deformities relies on an accurate 3-D reconstruction of bone structures from biplanar X-rays, which requires a precise detection and matching of anatomical structures in both views. In this paper, we propose a novel semiautomated technique for detecting complete scoliotic rib borders from PA-0° and PA-20° chest radiographs, by using an edge-following approach with multiple-path branching and oriented filtering. Edge-following processes are initiated from user starting points along upper and lower rib edges and the final rib border is obtained by finding the most parallel pair among detected edges. The method is based on a perceptual analysis leading to the assumption that no matter how bent a scoliotic rib is, it will always present relatively parallel upper and lower edges. The proposed method was tested on 44 chest radiographs of scoliotic patients and was validated by comparing pixels from all detected rib borders against their reference locations taken from the associated manually delineated rib borders. The overall 2-D detection accuracy was 2.64 ± 1.21 pixels. Comparing this accuracy level to reported results in the literature shows that the proposed method is very well suited for precisely detecting borders of scoliotic ribs from PA-0° and PA-20° chest radiographs.

A Novel Method for the 3-D Reconstruction of Scoliotic Ribs From Frontal and Lateral Radiographs

Among the external manifestations of scoliosis, the rib hump, which is associated with the ribs' deformities and rotations, constitutes the most disturbing aspect of the scoliotic deformity for patients. A personalized 3-D model of the rib cage is important for a better evaluation of the deformity, and hence, a better treatment planning. A novel method for the 3-D reconstruction of the rib cage, based only on two standard radiographs, is proposed in this paper. For each rib, two points are extrapolated from the reconstructed spine, and three points are reconstructed by stereo radiography. The reconstruction is then refined using a surface approximation. The method was evaluated using clinical data of 13 patients with scoliosis. A comparison was conducted between the reconstructions obtained with the proposed method and those obtained by using a previous reconstruction method based on two frontal radiographs. A first comparison criterion was the distances between the reconstructed ribs and the surface topography of the trunk, considered as the reference modality. The correlation between ribs axial rotation and back surface rotation was also evaluated. The proposed method successfully reconstructed the ribs of the 6th-12th thoracic levels. The evaluation results showed that the 3-D configuration of the new rib reconstructions is more consistent with the surface topography and provides more accurate measurements of ribs axial rotation.