Simultaneous adaptive smoothing of relaxometry and quantitative magnetization transfer mapping

Attempts for in-vivo histology require a high spatial resolution that comes with the price of a decreased signal-to-noise ratio. We present a novel iterative and multi-scale smoothing method for quantitative Magnetic Resonance Imaging (MRI) data that yield proton density, apparent transverse and longitudinal relaxation, and magnetization transfer maps. The method is based on the propagation-separation approach. The adaptivity of the procedure avoids the inherent bias from blurring subtle features in the calculated maps that is common for non-adaptive smoothing approaches. The characteristics of the methods were evaluated on a high-resolution data set (500 μ isotropic) from a single subject and quantified on data from a multi-subject study. The results show that the adaptive method is able to increase the signal-to-noise ratio in the calculated quantitative maps while largely avoiding the bias that is otherwise introduced by spatially blurring values across tissue borders. As a consequence, it preserves the intensity contrast between white and gray matter and the thin cortical ribbon

Development of a tool for automatic segmentation of the cerebellum in MR images of children

The human cerebellar cortex is a highly foliated structure that supports both motor and complex cognitive functions in humans. Magnetic Resonance Imaging (MRI) is commonly used to explore structural alterations in patients with psychiatric and neurological diseases. The ability to detect regional structural differences in cerebellar lobules may provide valuable insights into disease biology, progression and response to treatment, but has been hampered by the lack of appropriate tools for performing automated structural cerebellar segmentation and morphometry. In this thesis, time intensive manual tracings by an expert neuroanatomist of 16 cerebellar regions on high-resolution T1-weighted MR images of 18 children aged 9-13 years were used to generate the Cape Town Pediatric Cerebellar Atlas (CAPCA18) in the age-appropriate National Institute of Health Pediatric Database (NIHPD) asymmetric template space. An automated pipeline was developed to process the MR images and generate lobule-wise segmentations, as well as a measure of the uncertainty of the label assignments. Validation in an independent group of children with ages similar to those of the children used in the construction of the atlas, yielded spatial overlaps with manual segmentations greater than 70% in all lobules, except lobules VIIb and X. Average spatial overlap of the whole cerebellar cortex was 86%, compared to 78% using the alternative Spatially Unbiased Infra-tentorial Template (SUIT), which was developed using adult images

Physics-based Reconstruction Methods for Magnetic Resonance Imaging

Conventional Magnetic Resonance Imaging (MRI) is hampered by long scan times and only qualitative image contrasts that prohibit a direct comparison between different systems. To address these limitations, model-based reconstructions explicitly model the physical laws that govern the MRI signal generation. By formulating image reconstruction as an inverse problem, quantitative maps of the underlying physical parameters can then be extracted directly from efficiently acquired k-space signals without intermediate image reconstruction -- addressing both shortcomings of conventional MRI at the same time. This review will discuss basic concepts of model-based reconstructions and report about our experience in developing several model-based methods over the last decade using selected examples that are provided complete with data and code.Comment: 8 figures, review accepted to Philos. Trans. R. Soc.

Steady-state anatomical and quantitative magnetic resonance imaging of the heart using RF-frequencymodulated techniques

Cardiovascular disease (CVD) is the leading cause of death in the United States and Europe and generates healthcare costs of hundreds of billions of dollars annually. Conventional methods of diagnosing CVD are often invasive and carry risks for the patient. For example, the gold standard for diagnosing coronary artery disease, a major class of CVD, is x-ray coronary angiography, which has the disadvantages of being invasive, being expensive, using ionizing radiation, and having a ris k of complications. Conversely, coronary MR angiography (MRA) does not use ionizing radiation, can effectively visualize tissues without the need for exogenous contrast agents, and benefits from an adaptable temporal resolution. However, the acquisition time of cardiac MRI is far longer than the temporal scales of cardiac and respiratory motion, necessitating some method of compensating for this motion. The free-running framework is a novel development in our lab, benefitting from advances over the past three decades, that attempts to address disadvantages of previous cardiac MRI approaches: it provides fully self-gated 5D cardiac MRI with a simplified workflow, improved ease-of-use, reduced operator dependence, and automatic patient-specific motion detection. Free-running imaging increases the amount of information available to the clinician and is flexible enough to be translated to different app lications within cardiac MRI. Moreover, the self-gating of the free-running framework decoupled the acquisition from the motion compensation and thereby opened up cardiac MRI to the wider class of steady-state-based techniques utilizing balanced steady-state free precession (bSSFP) sequences, which have the benefits of practical simplicity and high signal-to-noise ratio. The focus of this thesis was therefore on the application of steady- state techniques to cardiac MRI. The first part addressed the long acquisition time of the current free-running framework and focused on anatomical coronary imaging. The published protocol of the free- running framework used an interrupted bSSFP acquisition where CHESS fat saturation modules were inserted to provide blood-fat contrast, as they suppress the signal of fat tissue surrounding the coronary arteries, and were followed by ramp-up pulses to reduce artefacts arising from the return to steady-state. This interrupted acquisition, however, suffered from an interrupted steady-state, reduced time efficiency, and higher specific absorption rate (SAR). Using novel lipid-insensitive binomial off-resonant RF excitation (LIBRE) pulses developed in our lab, the first project showed that LIBRE pulses incorporated into an uninterrupted free-running bSSFP sequence could be successfully used for 5D cardiac MRI at 1.5T. The free-running LIBRE approach reduced the acquisition time and SAR relative to the previous interrupted approach while maintaining image quality and vessel conspicuity. Furthermore, this had been the first successful use of a fat-suppressing RF excitation pulse in an uninterrupted bSSFP sequence for cardiac imaging, demonstrating that uninterrupted bSSFP can be used for cardiac MRI and addressing the problem of clinical sequence availability. Inspired by the feasibility of uninterrupted bSSFP for cardiac MRI, the second part investigated the potential of PLANET, a novel 3D multiparametric mapping technique, for free-running 5D myocardial mapping. PLANET utilizes a phase-cycled bSSFP acquisition and a direct ellipse-fitting algorithm to calculate T1 and T2 relaxation times, which suggested that it could be readily integrated into the free-running framework without interrupting the steady-state. After initially calibrating the acquisition, the possibility of accelerating the static PLANET acquisition was explored prior to applying it to the moving heart. It was shown that PLANET accuracy and precision could be maintained with two-fold acceleration with a 3D Cartesian spiral trajectory, suggesting that PLANET for myocardial mapping with the free-running 5D radial acquisition is feasible. Further work should investigate optimizing the reconstruction scheme, improving the coil sensitivity estimate, and examining the use of the radial trajectory with a view to implementing free-running 5D myocardial T1 and T2 mapping. This thesis presents two approaches utilizing RF-frequency-modulated steady-state techniques for cardiac MRI. The first approach involved the novel application of an uninterrupted bSSFP acquisition with off-resonant RF excitation for anatomical coronary imaging. The second approach investigated the use of phase-cycled bSSFP for free-running 5D myocardial T1 and T2 mapping. Both methods addressed the challenge of clinical availability of sequences in cardiac MRI, by showing that a common and simple sequence like bSSFP can be used for acquisition while the steps of motion compensation and reconstruction can be handled offline, and thus have the potential to improve adoption of cardiac MRI. -- Les maladies cardiovasculaires (MCV) représentent la principale cause de décès aux États-Unis et en Europe et génèrent des coûts de santé de plusieurs centaines de milliards de dollars par an. Les méthodes conventionnelles de diagnostic des MCV sont souvent invasives et comportent des risques pour le patient. Par exemple, la méthode de référence pour le diagnostic de la maladie coronarienne, une catégorie majeure de MCV, est la coronarographie par rayons X qui a comme inconvénients son caractère invasif, son coût, l’utilisation de rayonnements ionisants et le risque de complications. A l’inverse, l'angiographie coronarienne par résonance magnétique (ARM) n'utilise pas de rayonnements ionisants, permet de visualiser efficacement les tissus sans avoir recours à des agents de contraste exogènes et bénéficie d'une résolution temporelle ajustable. Cependant, le temps d'acquisition en IRM cardiaque est bien plus long que les échelles temporelles des mouvements cardiaques et respiratoires en jeu, ce qui rend la compensation de ces mouvements indispensable. Le cadre dit de « free -running » est un nouveau développement de notre laboratoire qui bénéficie des progrès réalisés au cours des trois dernières décennies et tente de remédier aux inconvénients des approches précédentes pour l'IRM cardiaque : il fournit une IRM cardiaque en cinq dimensions (5D) complètement « self-gated » , c’est-à-dire capable de détecter les mouvements cardiaques et respiratoires, forte d’une implémentation simplifiée, d’une plus grande facilité d'utilisation, d’une dépendance réduite vis-à-vis de l'opérateur et d’une détection automatique des mouvements spécifiques du patient. L'imagerie « free- running » augmente la quantité d'informations à disposition du clinicien et est suffisamment flexible pour être appliquée à différents domaines de l'IRM cardiaque. De plus, le « self-gating » du cadre « free-running » a découplé l'acquisition de la compensation de mouvement et a ainsi ouvert l'IRM cardiaque à la classe plus large des techniques basées sur l'état stationnaire utilisant des séquences de précession libre équilibrée en état stationnaire (bSSFP), qui se distinguent par leur simplicité d’utilisation et leur rapport signal sur bruit élevé. Le thème de cette thèse est donc l'application des techniques basées sur l'état stationnaire à l'IRM cardiaque. La première partie porte sur le long temps d'acquisition de l'actuel cadre « free-running» et se concentre sur l'imagerie anatomique coronaire. Le protocole publié utilise une acquisition bSSFP interrompue où des modules de saturation de graisse (CHESS) sont insérés de façon à fournir un contraste sang-graisse puisqu’ils suppriment le signal du tissu graisseux entourant les artères coronaires, et sont suivis par des impulsions en rampe pour réduire les artefacts résultant du retour à l'état stable. Cette acquisition interrompue souffre cependant d'un état d'équilibre interrompu, d'une efficacité temporelle réduite et d'un débit d'absorption spécifique (DAS) plus élevé. En utilisant les nouvelles impulsions d'excitation radiofréquence (RF) binomiales hors -résonance insensibles aux lipides (LIBRE) développées dans notre laboratoi re, ce premier projet montre que les impulsions LIBRE incorporées dans une séquence bSSFP ininterrompue et « free-running » peuvent être utilisées avec succès pour l'IRM cardiaque 5D à 1,5 T. L'approche « free-running LIBRE » permet de réduire le temps d'acquisition et le DAS par rapport à l'approche interrompue précédente, tout en maintenant la perceptibilité des artères coronariennes. En outre, il s'agit de la première utilisation réussie d'une impulsion d'excitation RF supprimant la graisse dans une séquence bSSFP ininterrompue pour l'imagerie cardiaque, ce qui démontre le potentiel d’utilisation de la séquence bSSFP ininterrompue pour l'IRM cardiaque et résout le problème de la disponibilité de la séquence en clinique. Inspirée par la faisabilité d’utilisation de la séquence bSSFP ininterrompue pour l'IRM cardiaque, la deuxième partie étudie le potentiel de PLANET, une nouvelle technique de cartographie 3D multiparamétrique, pour la cartographie 5D du myocarde via l’imagerie « free-running ». PLANET utilise une acquisition bSSFP à cycle de phase et un algorithme d'ajustement d'ellipse direct pour calculer les temps de relaxation T1 et T2, ce qui suggère que cette méthode pourrait être facilement intégrée au cadre « free - running » sans interruption de l’état d'équilibre. Après calibration de l'acquisition, nous explorons la possibilité d'accélérer l'acquisition statique de PLANET pour l'appliquer au cœur. Nous démontrons que l'exactitude et la précision de PLANET peuvent être maintenues pour une accélération double avec une trajectoire 3D cartésienne en spirale, ce qui suggère que PLANET est réalisable pour la cartographie du myocarde avec une acquisition radiale 5D « free-running ». D'autres travaux devraient porter sur l'optimisation du schéma de reconstruction, l'amélioration de l'estimation de la sensibilité de l’antenne et l'examen de l'utilisation de la trajectoire radiale en vue de la mise en œuvre de la cartographie 5D « free-running » T1 et T2 du myocarde. Cette thèse présente deux approches utilisant des techniques de modulation de fréquence radio en état stationnaire pour l'IRM cardiaque. La première approche implique l'application nouvelle d'une acquisition bSSFP ininterrompue avec une excitation RF hors résonance pour l'imagerie anatomique coronaire. La seconde approche porte sur l'utilisation d’une séquence bSSFP à cycle de phase pour la cartographie 5D T1 et T2 du myocarde. Ces deux méthodes permettent de répondre au défi posé par la disponibilité des séquences en IRM cardiaque en montrant qu'une séquence commune et simple comme la bSSFP peut être utilisée pour l'acquisition, tandis que les étapes de compensation du mouvement et de reconstruction peuvent être traitées hors ligne. Ainsi, ces méthodes ont le potentiel de favoriser l'adoption de l'IRM cardiaque

Multi-modality cardiac image computing: a survey

Multi-modality cardiac imaging plays a key role in the management of patients with cardiovascular diseases. It allows a combination of complementary anatomical, morphological and functional information, increases diagnosis accuracy, and improves the efficacy of cardiovascular interventions and clinical outcomes. Fully-automated processing and quantitative analysis of multi-modality cardiac images could have a direct impact on clinical research and evidence-based patient management. However, these require overcoming significant challenges including inter-modality misalignment and finding optimal methods to integrate information from different modalities. This paper aims to provide a comprehensive review of multi-modality imaging in cardiology, the computing methods, the validation strategies, the related clinical workflows and future perspectives. For the computing methodologies, we have a favored focus on the three tasks, i.e., registration, fusion and segmentation, which generally involve multi-modality imaging data, either combining information from different modalities or transferring information across modalities. The review highlights that multi-modality cardiac imaging data has the potential of wide applicability in the clinic, such as trans-aortic valve implantation guidance, myocardial viability assessment, and catheter ablation therapy and its patient selection. Nevertheless, many challenges remain unsolved, such as missing modality, modality selection, combination of imaging and non-imaging data, and uniform analysis and representation of different modalities. There is also work to do in defining how the well-developed techniques fit in clinical workflows and how much additional and relevant information they introduce. These problems are likely to continue to be an active field of research and the questions to be answered in the future