CNN-based Landmark Detection in Cardiac CTA Scans

Fast and accurate anatomical landmark detection can benefit many medical image analysis methods. Here, we propose a method to automatically detect anatomical landmarks in medical images. Automatic landmark detection is performed with a patch-based fully convolutional neural network (FCNN) that combines regression and classification. For any given image patch, regression is used to predict the 3D displacement vector from the image patch to the landmark. Simultaneously, classification is used to identify patches that contain the landmark. Under the assumption that patches close to a landmark can determine the landmark location more precisely than patches farther from it, only those patches that contain the landmark according to classification are used to determine the landmark location. The landmark location is obtained by calculating the average landmark location using the computed 3D displacement vectors. The method is evaluated using detection of six clinically relevant landmarks in coronary CT angiography (CCTA) scans: the right and left ostium, the bifurcation of the left main coronary artery (LM) into the left anterior descending and the left circumflex artery, and the origin of the right, non-coronary, and left aortic valve commissure. The proposed method achieved an average Euclidean distance error of 2.19 mm and 2.88 mm for the right and left ostium respectively, 3.78 mm for the bifurcation of the LM, and 1.82 mm, 2.10 mm and 1.89 mm for the origin of the right, non-coronary, and left aortic valve commissure respectively, demonstrating accurate performance. The proposed combination of regression and classification can be used to accurately detect landmarks in CCTA scans.Comment: This work was submitted to MIDL 2018 Conferenc

Radiographic Assessment of Hip Disease in Children with Cerebral Palsy: Development of a Core Measurement Set and Analysis of an Artificial Intelligence System

Cerebral palsy is the most common physical disability during childhood. Cerebral palsy related hip disease is caused by an imbalance of muscle forces, resulting in progressive migration of the hip to complete dislocation. This can decrease function and quality of life. The prevention of hip dislocation is possible if detected early. Therefore, surveillance programmes have been set up to monitor children with cerebral palsy enabling clinicians to intervene early and improve outcomes. Currently, hip disease is assessed by analysing pelvic radiographs with various geometric measurements. This time-consuming task is undertaken frequently when monitoring a child with cerebral palsy. This thesis aimed to identify the key radiographic parameters used by clinicians (the core measurement set), and then build an artificial intelligence system to automate the calculation of this core measurement set. A systematic review was conducted identifying a comprehensive list of previously reported measurements from studies measuring radiographic outcomes in cerebral palsy children with hip pathologies. Fifteen measurements were identified from the systematic review, of which Reimers’ migration percentage was the most commonly reported. These measurements were used to perform a two-round Delphi study among orthopaedic surgeons and physiotherapists. Participants rated the importance of each measurement using a nine-point Likert scale (‘not important’ to critically important’). After the two rounds of the Delphi process, Reimers’ migration percentage was included in the core measurement set. Following the final consensus meeting, the femoral head-shaft angle was also included. The anteroposterior pelvic radiographs of 1650 children were then used to build an artificial intelligence system integrating the core measurement set, in collaboration with engineers from the University of Manchester. The newly developed artificial intelligence system was assessed by comparing its ability to calculate measurements and outline the pelvis and femur on a radiograph. The reliability of the dataset used to train the model was also analysed. The proposed artificial intelligence model achieved a ‘good to excellent’ inter-observer reliability across 450 radiographs when comparing its ability to calculate Reimers’ migration percentage to five clinicians. Its ability to outline the pelvis and proximal femur was ‘adequate’ with the better performance observed in the pelvis than the femur. The reliability of the training dataset used to teach the artificial intelligence model was ‘good’ to ‘very good’. Artificial intelligence systems are feasible solutions to optimise the efficiency of hip radiograph analysis in cerebral palsy. Studies are warranted to include the core measurement set as a minimum when reporting on hip disease in cerebral palsy. Future research should investigate the feasibility of implementing a risk score to predict the likelihood of hip displacement

Automated brain lesion segmentation in magnetic resonance images

In this thesis, we investigate the potential of automation in brain lesion segmentation in magnetic resonance images. We first develop a novel supervised method, which segments regions in magnetic resonance images using gated recurrent units, provided training data with pixel-wise annotations on what to segment is available. We improve on this method using the latest technical advances in the field of machine learning and insights on possible weaknesses of our method, and adapt it specifically for the task of lesion segmentation in the brain. We show the feasibility of our approach on multiple public benchmarks, consistently reaching positions at the top of the list of competing methods. Adapting our problem successfully to the problem of landmark localization, we show the generalizability of the approach. Moving away from large training cohorts with manual segmentations to data where it is only known that a certain pathology is present, we propose a weakly-supervised segmentation approach. Given a set of images with known pathology of a certain kind and a healthy reference set, our formulation can segment the difference of the two data distributions. Lastly, we show how information from already existing lesion maps can be extracted in a meaningful way by connecting lesions across time in longitudinal studies. We hence present a full tool set for the automated processing of lesions in magnetic resonance images

Novel Approaches to the Representation and Analysis of 3D Segmented Anatomical Districts

Nowadays, image processing and 3D shape analysis are an integral part of clinical practice and have the potentiality to support clinicians with advanced analysis and visualization techniques. Both approaches provide visual and quantitative information to medical practitioners, even if from different points of view. Indeed, shape analysis is aimed at studying the morphology of anatomical structures, while image processing is focused more on the tissue or functional information provided by the pixels/voxels intensities levels. Despite the progress obtained by research in both fields, a junction between these two complementary worlds is missing. When working with 3D models analyzing shape features, the information of the volume surrounding the structure is lost, since a segmentation process is needed to obtain the 3D shape model; however, the 3D nature of the anatomical structure is represented explicitly. With volume images, instead, the tissue information related to the imaged volume is the core of the analysis, while the shape and morphology of the structure are just implicitly represented, thus not clear enough. The aim of this Thesis work is the integration of these two approaches in order to increase the amount of information available for physicians, allowing a more accurate analysis of each patient. An augmented visualization tool able to provide information on both the anatomical structure shape and the surrounding volume through a hybrid representation, could reduce the gap between the two approaches and provide a more complete anatomical rendering of the subject. To this end, given a segmented anatomical district, we propose a novel mapping of volumetric data onto the segmented surface. The grey-levels of the image voxels are mapped through a volume-surface correspondence map, which defines a grey-level texture on the segmented surface. The resulting texture mapping is coherent to the local morphology of the segmented anatomical structure and provides an enhanced visual representation of the anatomical district. The integration of volume-based and surface-based information in a unique 3D representation also supports the identification and characterization of morphological landmarks and pathology evaluations. The main research contributions of the Ph.D. activities and Thesis are: \u2022 the development of a novel integration algorithm that combines surface-based (segmented 3D anatomical structure meshes) and volume-based (MRI volumes) information. The integration supports different criteria for the grey-levels mapping onto the segmented surface; \u2022 the development of methodological approaches for using the grey-levels mapping together with morphological analysis. The final goal is to solve problems in real clinical tasks, such as the identification of (patient-specific) ligament insertion sites on bones from segmented MR images, the characterization of the local morphology of bones/tissues, the early diagnosis, classification, and monitoring of muscle-skeletal pathologies; \u2022 the analysis of segmentation procedures, with a focus on the tissue classification process, in order to reduce operator dependency and to overcome the absence of a real gold standard for the evaluation of automatic segmentations; \u2022 the evaluation and comparison of (unsupervised) segmentation methods, finalized to define a novel segmentation method for low-field MR images, and for the local correction/improvement of a given segmentation. The proposed method is simple but effectively integrates information derived from medical image analysis and 3D shape analysis. Moreover, the algorithm is general enough to be applied to different anatomical districts independently of the segmentation method, imaging techniques (such as CT), or image resolution. The volume information can be integrated easily in different shape analysis applications, taking into consideration not only the morphology of the input shape but also the real context in which it is inserted, to solve clinical tasks. The results obtained by this combined analysis have been evaluated through statistical analysis