Biomechanical Morphing for Personalized Fitting of Scoliotic Torso Skeleton Models

The use of patient-specific biomechanical models offers many opportunities in the treatment of adolescent idiopathic scoliosis, such as the design of personalized braces. The first step in the development of these patient-specific models is to fit the geometry of the torso skeleton to the patient’s anatomy. However, existing methods rely on high-quality imaging data. The exposure to radiation of these methods limits their applicability for regular monitoring of patients. We present a method to fit personalized models of the torso skeleton that takes as input biplanar low-dose radiographs. The method morphs a template to fit annotated points on visible portions of the spine, and it relies on a default biomechanical model of the torso for regularization and robust fitting of hardly visible parts of the torso skeleton, such as the rib cage. The proposed method provides an accurate and robust solution to obtain personalized models of the torso skeleton, which can be adopted as part of regular management of scoliosis patients. We have evaluated the method on ten young patients who participated in our study. We have analyzed and compared clinical metrics on the spine and the full torso skeleton, and we have found that the accuracy of the method is at least comparable to other methods that require more demanding imaging methods, while it offers superior robustness to artifacts such as interpenetration of ribs. Normal-dose X-rays were available for one of the patients, and for the other nine we acquired low-dose X-rays, allowing us to validate that the accuracy of the method persisted under less invasive imaging modalities

Scaled, patient-specific 3D vertebral model reconstruction based on 2D lateral fluoroscopy

Backgrounds: Accurate three-dimensional (3D) models of lumbar vertebrae are required for image-based 3D kinematics analysis. MRI or CT datasets are frequently used to derive 3D models but have the disadvantages that they are expensive, time-consuming or involving ionizing radiation (e.g., CT acquisition). An alternative method using 2D lateral fluoroscopy was developed. Materials and methods: A technique was developed to reconstruct a scaled 3D lumbar vertebral model from a single two-dimensional (2D) lateral fluoroscopic image and a statistical shape model of the lumbar vertebrae. Four cadaveric lumbar spine segments and two statistical shape models were used for testing. Reconstruction accuracy was determined by comparison of the surface models reconstructed from the single lateral fluoroscopic images to the ground truth data from 3D CT segmentation. For each case, two different surface-based registration techniques were used to recover the unknown scale factor, and the rigid transformation between the reconstructed surface model and the ground truth model before the differences between the two discrete surface models were computed. Results: Successful reconstruction of scaled surface models was achieved for all test lumbar vertebrae based on single lateral fluoroscopic images. The mean reconstruction error was between 0.7 and 1.6mm. Conclusions: A scaled, patient-specific surface model of the lumbar vertebra from a single lateral fluoroscopic image can be synthesized using the present approach. This new method for patient-specific 3D modeling has potential applications in spine kinematics analysis, surgical planning, and navigatio

Simulation of the postoperative trunk appearance in scoliosis surgery

Persistence of external trunk asymmetry after scoliosis surgical treatment is frequent and difficult to predict by clinicians. This is a significant problem considering that correction of the apparent deformity is a major factor of satisfaction for the patients. A simulation of the correction on the external appearance would allow the clinician to illustrate to the patient the potential result of the surgery and would help in deciding on a surgical strategy that could most improve his/her appearance. We describe a method to predict the scoliotic trunk shape after a spine surgical intervention. The capability of our method was evaluated using real data of scoliotic patients. Results of the qualitative evaluation were very promising and a quantitative evaluation based on the comparison of the simulated and the actual postoperative trunk surface showed an adequate accuracy for clinical assessment. The required short simulation time also makes our approach an eligible candidate for a clinical environment demanding interactive simulations.CIHR / IRS

Multimodal image fusion of anatomical structures for diagnosis, therapy planning and assistance

This paper provides an overview of work done in recent years by our research group to fuse multimodal images of the trunk of patients with Adolescent Idiopathic Scoliosis (AIS) treated at Sainte-Justine University Hospital Center (CHU). We first describe our surface acquisition system and introduce a set of clinical measurements (indices) based on the trunk's external shape, to quantify its degree of asymmetry. We then describe our 3D reconstruction system of the spine and rib cage from biplanar radiographs and present our methodology for multimodal fusion of MRI, X-ray and external surface images of the trunk We finally present a physical model of the human trunk including bone and soft tissue for the simulation of the surgical outcome on the external trunk shape in AIS.CIHR / IRS

Reconstruction 3D personnalisée de la colonne vertébrale à partir d'images radiographiques non-calibrées

Les systèmes de reconstruction stéréo-radiographique 3D -- La colonne vertébrale -- La scoliose idiopathique adolescente -- Évolution des systèmes de reconstruction 3D -- Filtres de rehaussement d'images -- Techniques de segmentation -- Les méthodes de calibrage -- Les méthodes de reconstruction 3D -- Problématique, hypothèses, objectifs et méthode générale -- Three-dimensional reconstruction of the scoliotic spine and pelvis from uncalibrated biplanar X-ray images -- A versatile 3D reconstruction system of the spine and pelvis for clinical assessment of spinal deformities -- Simulation experiments -- Clinical validation -- A three-dimensional retrospective analysis of the evolution of spinal instrumentation for the correction of adolescent idiopathic scoliosis -- Auto-calibrage d'un système à rayons-X à partir de primitives de haut niveau -- Segmentation de la colonne vertébrale -- Approche hiérarchique d'auto-calibrage d'un système d'acquisition à rayons-X -- Personalized 3D reconstruction of the scoliotic spine from hybrid statistical and X-ray image-based models -- Validation protocol

A physically based trunk soft tissue modeling for scoliosis surgery planning systems

One of the major concerns of scoliotic patients undergoing spinal correction surgery is the trunk's external appearance after the surgery. This paper presents a novel incremental approach for simulating postoperative trunk shape in scoliosis surgery. Preoperative and postoperative trunk shapes data were obtained using three-dimensional medical imaging techniques for seven patients with adolescent idiopathic scoliosis. Results of qualitative and quantitative evaluations, based on the comparison of the simulated and actual postoperative trunk surfaces, showed an adequate accuracy of the method. Our approach provides a candidate simulation tool to be used in a clinical environment for the surgery planning process.IRSC / CIH

Statistical model based 3D shape prediction of postoperative trunks for non-invasive scoliosis surgery planning

One of the major concerns of scoliosis patients undergoing surgical treatment is the aesthetic aspect of the surgery outcome. It would be useful to predict the postoperative appearance of the patient trunk in the course of a surgery planning process in order to take into account the expectations of the patient. In this paper, we propose to use least squares support vector regression for the prediction of the postoperative trunk 3D shape after spine surgery for adolescent idiopathic scoliosis. Five dimensionality reduction techniques used in conjunction with the support vector machine are compared. The methods are evaluated in terms of their accuracy, based on the leave-one-out cross-validation performed on a database of 141 cases. The results indicate that the 3D shape predictions using a dimensionality reduction obtained by simultaneous decomposition of the predictors and response variables have the best accuracy.CIHR / IRS

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.Adolescent idiopathic scoliosis (AIS) is a three dimensional (3D) deformity of the spine. For most patients, no surgical intervention is required. However, for patients with severe deformities, surgery is often recommended. Postero-anterior (PA) and antero-posterior (AP) x-rays are the most common modality for viewing and evaluating this deformity. From this imaging modality, clinical indices such as the Cobb angle can quantify and evaluate the severity of AIS. Clinical decision making is often based on this descriptor. However, the descriptors based on spinal radiographies are two- dimensional (2D) by nature. Therefore, they do not fully describe the deformity in AIS due to its three-dimensional (3D) nature. Poncet et al. have studied and presented geometric torsion as a 3D descriptor of AIS. This index measures a curve’s tendency to twist out of a plane. However, the method presented in their study is susceptible to errors from an imperfect 3D spinal reconstruction due to the local approach taken and is presented in a qualitative fashion. Hence, the objective of this study is to evaluate a new method of estimating geometric torsion in AIS employing parametric curve fitting techniques based on local arc-length approximations. The first study attempts to evaluate the sensitivity of the presented method of estimating geometric torsion against noisy data or 3D reconstruction errors. Two clinical studies will then present this descriptor as a quantitative measurement of AIS and will attempt to identify potential new sub-groups and demonstrate the clinical relevance of these new sub-groups. vi The first study evaluated the robustness of the new method in estimating geometric torsion in the presence of reconstruction errors. This study demonstrated that the new method is robust to 3D spinal reconstruction errors and achieves quantitative measures in a global fashion. The second study evaluated this new method of estimating geometric torsion in patients with Lenke type 1 deformities. This study identified two sub-groups based on torsion values, a high torsion and a low torsion group. These two sub-groups showed differences in the orientation of the planes of maximum deformity (PMC) in the thoraco-lumbar segment of the spine with the high torsion group having greater values. The last study evaluated the surgical results of patients with Lenke type 1 deformities having been sub-classified in high and low torsion groups. This study showed differences in TLL PMC with the high torsion group of patients having higher values pre and post-operatively, These studies present a novel method of estimating geometric torsion in AIS and present this 3D descriptor quantitatively. They have demonstrated the existence of new sub-groups within current classification systems that were previously undetected and have shown the clinical relevance of this new method of estimating geometric torsion in AIS. This project contributes towards the development of new 3D indices for AIS and opens the door to potential new 3D classifications