Machine learning to automatically detect anatomical landmarks on phase-contrast enhanced magnetic resonance angiography

Abstract

Abstract Funding Acknowledgements Type of funding sources: Public grant(s) – National budget only. Main funding source(s): Guala A. received funding from the Spanish Ministry of Science, Innovation and Universities Introduction The heterogeneous characteristic of the thoracic aorta implies that all biomarkers with potential for risk stratification need to be references to a specific location. This is the case, for example, of diameter [1], stiffness [2] and wall shear stress [3]. This is normally achieved by the manual identification of a limited number of key anatomic landmarks [4], which is a time-demanding task and may impact biomarkers accuracy and reproducibility. Automatic identification of these anatomic landmarks may speed-up the analysis and allow for the creation of fully automatic image analysis pipelines. Machine learning (ML) algorithms might be suitable for this task. Purpose The aim of this study was to test the performance of a ML algorithm in localizing key thoracic anatomical landmarks on phase-contrast enhanced magnetic resonance angiograms (PC-MRA). Methods PC-MRA of 323 patients with native aorta and aortic valve and a variety of aortic conditions (141 bicuspid aortic valve patients, 60 patients with degenerative aortic aneurysms, 82 patients with genetic aortopathy and 40 healthy volunteers) were included in this study. Four anatomical landmarks were manually identified on PC-MRA by an experienced researcher: sinotubular junction, the pulmonary artery bifurcation and the first and third supra-aortic vessel braches. A reinforcement learning algorithm (DQN), combining Q-learning with deep neural networks, was trained. The algorithm was tested in a separate set of 30 PC-MRA with similar distribution of aortic conditions in which human intra-observer reproducibility was quantified. The distance between points was used as quality metric and human annotation was considered as ground-truth. Repeated-measures ANOVA was used for statistical testing. Results ML algorithm resulted in performance similar to the intra-observer variability obtained by the experienced human reader in the identification of the sinotubular junction (11.1 ± 8.6 vs 11.0 ± 8.1 mm, p = 0.949) and first (6.8 ± 5.6 vs 6.6 ± 3.9 mm, p = 0.886) and third (8.4 ± 7.4 vs 6.8 ± 4.0 mm, p = 0.161) supra-aortic vessels branches. However, the algorithm did not reach human-level performance in the localization of the pulmonary artery bifurcation (15.2 ± 13.1 vs 10.2 ± 7.0 mm, p = 0.008). The time needed to the ML algorithm to locate all points ranged between 0.8 and 1.6 seconds on a standard computer while manual annotation required around two minutes to be performed. Conclusions The rapid identification of key aortic anatomical landmarks by a reinforced learning algorithm is feasible with human-level performance. This approach may thus be used for the design of fully-automatic pipeline for 4D flow CMR analysis.