Exploring abnormal arch shape patterns using CMR-based hierarchical 3D shape clustering: Application to a generic imaging population of repaired coarctation of the aorta
Jan L Bruse1, Abbas Khushnood1, Tain-Yen Hsia1, Andrew M Taylor1, Vivek Muthurangu1, and Silvia Schievano1

1Centre for Cardiovascular Imaging, UCL Institute of Cardiovascular Science & Great Ormond Street Hospital for Children, London, United Kingdom

Synopsis

We present a novel method for hierarchical 3D shape clustering of aortic arch shape models segmented and reconstructed from CMR imaging data. We apply the method to a cohort of 45 patients post aortic coarctation repair in order to explore previously unknown arch shape patterns that may relate to clinical outcome. Exploring a pathologic shape population using data mining and statistical shape modeling techniques can provide novel insight for improved diagnosis and treatment strategies and can thereby assisst in clinical decision making when analysing complex cases.

PURPOSE

Despite improving survival rates, patients diagnosed with aortic coarctation (CoA) suffer from late complications post CoA repair and regular screening via cardiac imaging is warranted [1]. Cardiovascular magnetic resonance (CMR) imaging being a cost-effective and accurate tool to provide detailed 3D anatomical shape information, imaging databases are growing and novel approaches to structure such data are required. We hypothesised that combining statistical shape modelling (SSM) with hierarchical clustering techniques from the field of data mining may aid in discovering previously unknown shape patterns in a population of aortic arches.

METHODS

A retrospective cohort of 45 CoA patients post repair (aged 15-38 years; BSA 1.47–2.2m2, Fig. 1) and 20 healthy Control subjects (aged 10-19 years; BSA 1.3-2.3m²) who underwent routine CMR imaging (1.5T Avanto scanner, Siemens Medical Solutions, Erlangen, Germany; 3D balanced steady-state free precession sequence, SSFP) were included in the study. Aortas were segmented semi-automatically [2] from the CMR data and reconstructed [3] to obtain 3D surface models of the arch. The arches were cut at the level of the sub-annular plane and the level of the diaphragm; head and neck vessels and coronary arteries were cut as close as possible to the arch. The meshed 3D surface models constituted the input for the SSM framework Deformetrica [4], which essentially computes the mean anatomical shape (atlas) and its deformations towards each subject shape within a population of surface shapes without requiring any landmarking. Our approach was two-fold: First we computed the atlas of the Control population as an ideal reference shape (Fig. 2). All CoA arch models were then rigidly registered [5] to the Control atlas and the surface deformations of the atlas towards each CoA subject were calculated in Deformetrica (Fig. 3). The set of all subject-specific deformation vectors constituted the input for a hierarchical clustering algorithm (MATLAB, Natick, MA) that does not require specifying an expected number of clusters prior to calculation. Results were visualised as a dendrogram where objects are grouped together based on their arch shape similarity. Cutting the dendrogram horizontally at a particular height partitions the data into shape sub-groups [6]. The obtained arch clusters were compared with regard to traditional morphometric parameters (surface to volume ratio SVol; centreline tortuosity CLT; minimum to maximum diameter ratio Dmin/Dmax; ascending aortic diameter to proximal descending diameter ratio Dasc/Ddesc) and functional parameters derived from CMR (left ventricular ejection fraction LVEF; indexed end-diastolic volume iLVEDV; left ventricular mass LVmass) using Kruskal-Wallis tests (IBM SPSS Statistics, Chicago, IL).

RESULTS

The clustering yielded two main aortic arch shape clusters, with one splitting into two large subgroups when cut at the same height of the dendrogram (Fig. 4). Distributions of the measured morphometric parameters SVol (p < .001), Dmin/Dmax (p = .006) and Dasc/Ddesc (p = .005) were significantly different between the three main shape clusters. Interestingly, distributions of iLVEDV were significantly different (p = .032) between the three groups, with larger values occurring in group three, which clustered subjects with a dilated aortic root together.

DISCUSSION

We present the first application of data mining techniques and 3D statistical shape modelling to a population of aortic arch shapes post CoA repair. Visually, similar arch shapes were clustered together correctly, which was corroborated by traditional shape parameters being different between the groups. Found differences between clusters in iLVEDV suggest that clinical outcomes can be associated with arch shape features detected via the proposed deformation-based shape clustering. Studies in a larger cohort and including more clinically relevant data are necessary to support our findings. Yet, particularly the high degree of automation makes our method an attractive tool for the retrieval of similarly shaped subjects from medical imaging databases for clinical decision support [7].

CONCLUSION

The combination of SSFP CMR data acquisition and advanced 3D shape modelling with data mining techniques provides a sensitive analytical platform, which could improve and facilitate risk assessment, follow-up strategies and treatment planning in complex cardiac lesions.

Acknowledgements

This report incorporates independent research from the National Institute for Health ResearchBiomedical Research Centre Funding Scheme. The views expressed in this publication arethose of the author(s) and not necessarily those of the NHS, the National Institute for HealthResearch or the Department of Health.

The authors gratefully acknowledge support from Fondation Leducq.

References

1. Brown, M.L., Burkhart, H.M., Connolly, H.M., Dearani, J.A., Cetta, F., Li, Z., Oliver, W.C., Warnes, C.A., Schaff, H.V.: Coarctation of the Aorta: Lifelong Surveillance Is Mandatory Following Surgical Repair. Journal of the American College of Cardiology. 62, 1020–1025 (2013).

2. Zuluaga, M.A., Cardoso, M.J., Modat, M., Ourselin, S.: Multi-atlas Propagation Whole Heart Segmentation from MRI and CTA Using a Local Normalised Correlation Coefficient Criterion. In: Ourselin, S., Rueckert, D., and Smith, N. (eds.) Functional Imaging and Modeling of the Heart. pp. 174–181. Springer Berlin Heidelberg (2013).

3. Yushkevich, P.A., Piven, J., Hazlett, H.C., Smith, R.G., Ho, S., Gee, J.C., Gerig, G.: User-guided 3D active contour segmentation of anatomical structures: Significantly improved efficiency and reliability. NeuroImage. 31, 1116–1128 (2006).

4. Durrleman, S., Prastawa, M., Charon, N., Korenberg, J.R., Joshi, S., Gerig, G., Trouvé, A.: Morphometry of anatomical shape complexes with dense deformations and sparse parameters. NeuroImage. 101, 35–49 (2014).

5. Besl, P.J., McKay, N.D.: A method for registration of 3-D shapes. IEEE Transactions on Pattern Analysis and Machine Intelligence. 14, 239–256 (1992).

6. Hastie, T., Tibshirani, R., Friedman, J.: The Elements of Statistical Learning. Springer New York, New York, NY (2009).

7. Tsymbal, A., Meissner, E., Kelm, M., Kramer, M.: Towards cloud-based image-integrated similarity search in big data. In: 2014 IEEE-EMBS International Conference on Biomedical and Health Informatics (BHI). pp. 593–596 (2014).

Figures

Control cohort of healthy subjects and 3D mean shape (atlas) computed using Deformetrica

Individual subject shapes are characterised by their unique set of deformation vectors φ, registering the Control atlas towards each subject shape rather than their actual point coordinates

Final clustering of the CoA arch shapes in three subgroups obtained by cutting the computed dendrogram

3D aortic arch surface models reconstructed from CMR data of 45 CoA patients included in the study (random order)



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
4319