Medical image analysis in clinical practice is commonly carried out on 2D image data, without fully exploiting the detailed 3D anatomical information that is provided by modern non-invasive medical imaging techniques. In this paper, a statistical shape analysis method is presented, which enables the extraction of 3D anatomical shape features from cardiovascular magnetic resonance (CMR) image data, with no need for manual landmarking. The method was applied to repaired aortic coarctation arches that present complex shapes, with the aim of capturing shape features as biomarkers of potential functional relevance. The method is presented from the user-perspective and is evaluated by comparing results with traditional morphometric measurements. Steps required to set up the statistical shape modelling analyses, from pre-processing of the CMR images to parameter setting and strategies to account for size differences and outliers, are described in detail. The anatomical mean shape of 20 aortic arches post-aortic coarctation repair (CoA) was computed based on surface models reconstructed from CMR data. By analysing transformations that deform the mean shape towards each of the individual patient’s anatomy, shape patterns related to differences in body surface area (BSA) and ejection fraction (EF) were extracted. The resulting shape vectors, describing shape features in 3D, were compared with traditionally measured 2D and 3D morphometric parameters. The computed 3D mean shape was close to population mean values of geometric shape descriptors and visually integrated characteristic shape features associated with our population of CoA shapes. After removing size effects due to differences in body surface area (BSA) between patients, distinct 3D shape features of the aortic arch correlated significantly with EF (r = 0.521, p = .022) and were well in agreement with trends as shown by traditional shape descriptors. The suggested method has the potential to discover previously unknown 3D shape biomarkers from medical imaging data. Thus, it could contribute to improving diagnosis and risk stratification in complex cardiac disease.

Objective: Today's growing medical image databases call for novel processing tools to structure the bulk of data and extract clinically relevant information. Unsupervised hierarchical clustering may reveal clusters within anatomical shape data of patient populations as required for modern precision medicine strategies. Few studies have applied hierarchical clustering techniques to three-dimensional patient shape data and results depend heavily on the chosen clustering distance metrics and linkage functions. In this study, we sought to assess clustering classification performance of various distance/linkage combinations and of different types of input data to obtain clinically meaningful shape clusters. Methods: We present a processing pipeline combining automatic segmentation, statistical shape modeling, and agglomerative hierarchical clustering to automatically subdivide a set of 60 aortic arch anatomical models into healthy controls, two groups affected by congenital heart disease, and their respective subgroups as defined by clinical diagnosis. Results were compared with traditional morphometrics and principal component analysis of shape features. Results: Our pipeline achieved automatic division of input shape data according to primary clinical diagnosis with high F-score (0.902 ± 0.042) and Matthews correlation coefficient (0.851 ± 0.064) using the correlation/weighted distance/linkage combination. Meaningful subgroups within the three patient groups were obtained and benchmark scores for automatic segmentation and classification performance are reported. Conclusion: Clustering results vary depending on the distance/linkage combination used to divide the data. Yet, clinically relevant shape clusters and subgroups could be found with high specificity and low misclassification rates. Significance: Detecting disease-specific clusters within medical image data could improve image-based risk assessment, treatment planning, and medical device development in complex disease.

Computational models of congenital heart disease (CHD) have become increasingly sophisticated over the last 20 years. They can provide an insight into complex flow phenomena, allow for testing devices into patient-specific anatomies (pre-CHD or post-CHD repair) and generate predictive data. This has been applied to different CHD scenarios, including patients with single ventricle, tetralogy of Fallot, aortic coarctation and transposition of the great arteries. Patient-specific simulations have been shown to be informative for preprocedural planning in complex cases, allowing for virtual stent deployment. Novel techniques such as statistical shape modelling can further aid in the morphological assessment of CHD, risk stratification of patients and possible identification of new ‘shape biomarkers’. Cardiovascular statistical shape models can provide valuable insights into phenomena such as ventricular growth in tetralogy of Fallot, or morphological aortic arch differences in repaired coarctation. In a constant move towards more realistic simulations, models can also account for multiscale phenomena (eg, thrombus formation) and importantly include measures of uncertainty (ie, CIs around simulation results). While their potential to aid understanding of CHD, surgical/procedural decision-making and personalisation of treatments is undeniable, important elements are still lacking prior to clinical translation of computational models in the field of CHD, that is, large validation studies, cost-effectiveness evaluation and establishing possible improvements in patient outcomes.

Three-dimensional (3D) imaging is an important tool for diagnostics, surgical planning, and evaluation of surgical outcomes in craniofacial procedures. Gold standard for acquiring 3D imaging is computed tomography that entails ionizing radiations and, in young children, a general anaesthesia. Three-dimensional photographic imaging is an alternative method to assess patients who have undergone calvarial reconstructive surgery. The aim of this study was to assess the utility of 3D handheld scanning photography in a cohort of patients who underwent spring-assisted correction surgery for scaphocephaly. Pre- and postoperative 3D scans acquired in theater and at the 3-week follow-up in clinic were postprocessed for 9 patients. Cephalic index (CI), head circumference, volume, sagittal length, and coronal width over the head at pre-op, post-op, and follow-up were measured from the 3D scans. Cephalic index from 3D scans was compared with measurements from planar x-rays. Statistical shape modeling (SSM) was used to calculate the 3D mean anatomical head shape of the 9 patients at the pre-op, post-op, and follow-up. No significant differences were observed in the CI between 3D and x-ray. Cephalic index, volume, and coronal width increased significantly over time. Mean shapes from SSM visualized the overall and regional 3D changes due to the expansion of the springs in situ. Three-dimensional handheld scanning followed by SSM proved to be an efficacious and practical method to evaluate 3D shape outcomes after spring-assisted cranioplasty in individual patients and the population.

Trigonocephaly in patients with metopic synostosis is corrected by fronto-orbital remodelling (FOR). The aim of this study was to quantitatively assess aesthetic outcomes of FOR by capturing 3D forehead scans of metopic patients pre- and post-operatively and comparing them with controls. Ten single-suture metopic patients undergoing FOR and 15 age-matched non-craniosynostotic controls were recruited at Great Ormond Street Hospital for Children (UK). Scans were acquired with a three-dimensional (3D) handheld camera and post-processed combining 3D imaging software. 3D scans were first used for cephalometric measurements. Statistical shape modelling was then used to compute the 3D mean head shapes of the three groups (FOR pre-op, post-op and controls). Head shape variations were described via principal component analysis (PCA). Cephalometric measurements showed that FOR significantly increased the forehead volume and improved trigonocephaly. This improvement was supported visually by pre- and post-operative computed mean 3D shapes and numerically by PCA (p < 0.001). Compared with controls, post-operative scans showed flatter foreheads (p < 0.001). In conclusion, 3D scanning followed by 3D statistical shape modelling enabled the 3D comparison of forehead shapes of metopic patients and non-craniosynostotic controls, and demonstrated that the adopted FOR technique was successful in correcting bitemporal narrowing but overcorrected the rounding of the forehead.