Understanding the influence of surgical parameters on craniofacial surgery outcomes: a computational study.

Sagittal craniosynostosis (SC) is a congenital condition whereby the newborn skull develops abnormally owing to the premature ossification of the sagittal suture. Spring-assisted cranioplasty (SAC) is a minimally invasive surgical technique to treat SC, where metallic distractors are used to reshape the newborn's head. Although safe and effective, SAC outcomes remain uncertain owing to the limited understanding of skull-distractor interaction and the limited information provided by the analysis of single surgical cases. In this work, an SC population-averaged skull model was created and used to simulate spring insertion by means of the finite-element analysis using a previously developed modelling framework. Surgical parameters were varied to assess the effect of osteotomy and spring positioning, as well as distractor combinations, on the final skull dimensions. Simulation trends were compared with retrospective measurements from clinical imaging (X-ray and three-dimensional photogrammetry scans). It was found that the on-table post-implantation head shape change is more sensitive to spring stiffness than to the other surgical parameters. However, the overall end-of-treatment head shape is more sensitive to spring positioning and osteotomy size parameters. The results of this work suggest that SAC surgical planning should be performed in view of long-term results, rather than immediate on-table reshaping outcomes.

The Esthetic Perception of Morphological Severity in Scaphocephalic Patients is Correlated With Specific Head Geometrical Features.

OBJECTIVE: To investigate the relationship between perception of craniofacial deformity, geometric head features, and 3D head shape analyzed by statistical shape modeling (SSM). PATIENTS: A total of 18 unoperated patients with scaphocephaly (age = 5.2 ± 1.1m)-6 were followed-up after spring-assisted cranioplasty (SAC) (age = 9.6 ± 1.5m)-and 6 controls (age = 6.7 ± 2.5m). MAIN OUTCOME MEASURES: 3D head shapes were retrieved from 3D scans or computed tomography (CTs). Various geometrical features were measured: anterior and posterior prominence, take-off angle, average anterior and posterior lateral and horizontal curvatures, cranial index (CI) (cranial width over length), and turricephaly index (TI) (cranial height over length). SSM and principal component analysis (PCA) described shape variability. All models were 3D printed; the perception of deformity was blindly scored by 9 surgeons and 1 radiologist in terms of frontal bossing (FB), occipital bulleting (OB), biparietal narrowing (BN), low posterior vertex (LPV), and overall head shape (OHS). RESULTS: A moderate correlation was found between FB and anterior prominence (r = 0.56, P < .01) and take-off angle (r = - 0.57, P < .01). OB correlated with average posterior lateral curvature (r = 0.43, P < 0.01) similarly to BPN (r = 0.55, P < .01) and LPV (r = 0.43, P < .01). OHS showed strong correlation with CI (r = - 0.68, P < .01) and TI (r = 0.63, P< .01). SSM Mode 1 correlated with OHS (r = 0.66, p < .01) while Mode 3 correlated with FB (r = - 0.58, P < .01). CONCLUSIONS: Esthetic cranial appearance in craniofacial patients is correlated to specific geometric parameters and could be estimated using automated methods such as SSM.

Parametrizing the genioplasty: a biomechanical virtual study on soft tissue behavior.

PURPOSE: Sliding genioplasty is used to surgically correct a retruded or misaligned chin: in this procedure, an osteotomy is performed and the bony segment is repositioned. In this study we investigate the effect of surgical parameters (bony segment movement, osteotomy design) on postop soft tissue changes in a patient cohort. METHODS: Seven patients were retrospectively recruited. Cone beam computed tomography data were obtained and soft tissue and bone shape reconstructions were performed. 3D models were created and surgical cuts were replicated according to postop scans. Each model was imported in ANSYS 2019R1 (Ansys Inc, USA) for simulation: the effect of variation in osteotomy plane as well as extent of bony segment movement were assessed by means of design of experiment: surgical parameters were varied in a surgically acceptable range and the soft tissue predictions were evaluated as displacement output of five craniometric landmarks. RESULTS: Simulation results show the overall changes of the lower third of the face are sensitive to changes in horizontal and vertical displacement of the bony segment as well as segment rotation. No significant changes in the soft tissue response were to attribute to the osteotomy design. CONCLUSIONS: Our results are consistent with experimental findings reported in the literature: when planning genioplasty in orthognathic surgery, particular focus on the segment movement (horizontal translation, vertical translation and rotation), rather than on the design of the osteotomy itself, should be considered.

A Proof of Concept of a Non-Invasive Image-Based Material Characterization Method for Enhanced Patient-Specific Computational Modeling.

PURPOSE: Computational models of cardiovascular structures rely on their accurate mechanical characterization. A validated method able to infer the material properties of patient-specific large vessels is currently lacking. The aim of the present study is to present a technique starting from the flow-area (QA) method to retrieve basic material properties from magnetic resonance (MR) imaging. METHODS: The proposed method was developed and tested, first, in silico and then in vitro. In silico, fluid-structure interaction (FSI) simulations of flow within a deformable pipe were run with varying elastic modules (E) between 0.5 and 32 MPa. The proposed QA-based formulation was assessed and modified based on the FSI results to retrieve E values. In vitro, a compliant phantom connected to a mock circulatory system was tested within MR scanning. Images of the phantom were acquired and post-processed according to the modified formulation to infer E of the phantom. Results of in vitro imaging assessment were verified against standard tensile test. RESULTS: In silico results from FSI simulations were used to derive the correction factor to the original formulation based on the geometrical and material characteristics. In vitro, the modified QA-based equation estimated an average E = 0.51 MPa, 2% different from the E derived from tensile tests (i.e. E = 0.50 MPa). CONCLUSION: This study presented promising results of an indirect and non-invasive method to establish elastic properties from solely MR images data, suggesting a potential image-based mechanical characterization of large blood vessels.

Developing and testing an algorithm for automatic segmentation of the fetal face from three-dimensional ultrasound images.

Fetal craniofacial abnormalities are challenging to detect and diagnose on prenatal ultrasound (US). Image segmentation and computer analysis of three-dimensional US volumes of the fetal face may provide an objective measure to quantify fetal facial features and identify abnormalities. We have developed and tested an atlas-based partially automated facial segmentation algorithm; however, the volumes require additional manual segmentation (MS), which is time and labour intensive and may preclude this method from clinical adoption. These manually refined segmentations can then be used as a reference (atlas) by the partially automated segmentation algorithm to improve algorithmic performance with the aim of eliminating the need for manual refinement and developing a fully automated system. This study assesses the inter- and intra-operator variability of MS and tests an optimized version of our automatic segmentation (AS) algorithm. The manual refinements of 15 fetal faces performed by three operators and repeated by one operator were assessed by Dice score, average symmetrical surface distance and volume difference. The performance of the partially automatic algorithm with difference size atlases was evaluated by Dice score and computational time. Assessment of the manual refinements showed low inter- and intra-operator variability demonstrating its suitability for optimizing the AS algorithm. The algorithm showed improved performance following an increase in the atlas size in turn reducing the need for manual refinement.

A novel RBF-based predictive tool for facial distraction surgery in growing children with syndromic craniosynostosis.

PURPOSE: Predicting changes in face shape from corrective surgery is challenging in growing children with syndromic craniosynostosis. A prediction tool mimicking composite bone and skin movement during facial distraction would be useful for surgical audit and planning. To model surgery, we used a radial basis function (RBF) that is smooth and continuous throughout space whilst corresponding to measured distraction at landmarks. Our aim is to showcase the pipeline for a novel landmark-based, RBF-driven simulation for facial distraction surgery in children. METHODS: An individual's dataset comprised of manually placed skin and bone landmarks on operated and unoperated regions. Surgical warps were produced for 'older' monobloc, 'older' bipartition and 'younger' bipartition groups by applying a weighted least-squares RBF fitted to the average landmarks and change vectors. A 'normalisation' warp, from fitting an RBF to craniometric landmark differences from the average, was applied to each dataset before the surgical warp. The normalisation was finally reversed to obtain the individual prediction. Predictions were compared to actual post-operative outcomes. RESULTS: The averaged change vectors for all groups showed skin and bone movements characteristic of the operations. Normalisation for shape-size removed individual asymmetry, size and proportion differences but retained typical pre-operative shape features. The surgical warps removed the average syndromic features. Reversing the normalisation reintroduced the individual's variation into the prediction. The mid-facial regions were well predicted for all groups. Forehead and brow regions were less well predicted. CONCLUSIONS: Our novel, landmark-based, weighted RBF can predict the outcome for facial distraction in younger and older children with a variety of head and face shapes. It can replicate the surgical reality of composite bone and skin movement jointly in one model. The potential applications include audit of existing patient outcomes, and predicting outcome for new patients to aid surgical planning.

Three-dimensional soft tissue prediction in orthognathic surgery: a clinical comparison of Dolphin, ProPlan CMF, and probabilistic finite element modelling.

Three-dimensional surgical planning is used widely in orthognathic surgery. Although numerous computer programs exist, the accuracy of soft tissue prediction remains uncertain. The purpose of this study was to compare the prediction accuracy of Dolphin, ProPlan CMF, and a probabilistic finite element method (PFEM). Seven patients (mean age 18years; five female) who had undergone Le Fort I osteotomy with preoperative and 1-year postoperative cone beam computed tomography (CBCT) were included. The three programs were used for soft tissue prediction using planned and postoperative maxillary position, and these were compared to postoperative CBCT. Accurate predictions were obtained with each program, indicated by root mean square distances: RMSDolphin=1.8±0.8mm, RMSProPlan=1.2±0.4mm, and RMSPFEM=1.3±0.4mm. Dolphin utilizes a landmark-based algorithm allowing for patient-specific bone-to-soft tissue ratios, which works well for cephalometric radiographs but has limited three-dimensional accuracy, whilst ProPlan and PFEM provide better three-dimensional predictions with continuous displacements. Patient or population-specific material properties can be defined in PFEM, while no soft tissue parameters are adjustable in ProPlan. Important clinical considerations are the topological differences between predictions due to the three algorithms, the non-negligible influence of the mismatch between planned and postoperative maxillary position, and the learning curve associated with sophisticated programs like PFEM.

Comparison and calibration of a real-time virtual stenting algorithm using Finite Element Analysis and Genetic Algorithms.

In this paper, we perform a comparative analysis between two computational methods for virtual stent deployment: a novel fast virtual stenting method, which is based on a spring-mass model, is compared with detailed finite element analysis in a sequence of in silico experiments. Given the results of the initial comparison, we present a way to optimise the fast method by calibrating a set of parameters with the help of a genetic algorithm, which utilises the outcomes of the finite element analysis as a learning reference. As a result of the calibration phase, we were able to substantially reduce the force measure discrepancy between the two methods and validate the fast stenting method by assessing the differences in the final device configurations.

Printed three-dimensional airway model assists planning of single-lung ventilation in a small child.

BACKGROUND: Single-lung ventilation in infants and small children is challenging because suitable sizes of double-lumen cuffed tracheal tubes are not available. A 6-yr-old child required pulmonary saline washout for primary alveolar proteinosis, and therefore needed sequential single-lung ventilation in order to achieve safe oxygenation. Before undertaking this potentially hazardous procedure, we practised bronchial intubation on an anatomical model of her airway constructed from computed tomography (CT) data. METHODS: We created a full-scale, anatomically accurate, transparent plastic model of the trachea and main bronchi on a three-dimensional printer using data from a CT scan. We then performed several different airway approaches to identify those likely to be most suitable, ex vivo, before the clinical procedure was carried out on the patient. RESULTS: The model helped us to choose the type and size of bronchial tubes and to practise their insertion beforehand. Subsequently, during anaesthesia, the chosen technique was successful. CONCLUSIONS: Three-dimensional printing of a model of the airway of a small child aided planning of bronchial intubation and single-lung ventilation. Three-dimensional printing of airway structures may have wider application in anaesthesia practice.

Patient-specific simulations of transcatheter aortic valve stent implantation.

Transcatheter aortic valve implantation (TAVI) enables treatment of aortic stenosis with no need for open heart surgery. According to current guidelines, only patients considered at high surgical risk can be treated with TAVI. In this study, patient-specific analyses were performed to explore the feasibility of TAVI in morphologies, which are currently borderline cases for a percutaneous approach. Five patients were recruited: four patients with failed bioprosthetic aortic valves (stenosis) and one patient with an incompetent, native aortic valve. Three-dimensional models of the implantation sites were reconstructed from computed tomography images. Within these realistic geometries, TAVI with an Edwards Sapien stent was simulated using finite element (FE) modelling. Engineering and clinical outcomes were assessed. In all patients, FE analysis proved that TAVI was morphologically feasible. After the implantation, stress distribution showed no risks of immediate device failure and geometric orifice areas increased with low risk of obstruction of the coronary arteries. Maximum principal stresses in the arterial walls were higher in the model with native outflow tract. FE analyses can both refine patient selection and characterise device mechanical performance in TAVI, overall impacting on procedural safety in the early introduction of percutaneous heart valve devices in new patient populations.

Use of rapid prototyping models in the planning of percutaneous pulmonary valved stent implantation.

Percutaneous replacement of the pulmonary valve is a recently developed interventional technique which involves the implantation of a valved stent in the pulmonary trunk. It relies upon careful consideration of patient anatomy for both stent design and detailed procedure planning. Medical imaging data in the form of two-dimensional scans and three-dimensional interactive graphics offer only limited support for these tasks. The paper reports the results of an experimental investigation on the use of arterial models built by rapid prototyping techniques. An analysis of clinical needs has helped to specify proper requirements for such model properties as cost, strength, accuracy, elastic compliance, and optical transparency. Two different process chains, based on the fused deposition modelling technique and on the vacuum casting of thermoset resins in rubber moulds, have been tested for prototype fabrication. The use of anatomical models has allowed the cardiologist's confidence in patient selection, prosthesis fabrication, and final implantation to be significantly improved.