A hybrid hierarchical strategy for registration of 7T TOF-MRI to 7T PC-MRI intracranial vessel data.

PURPOSE: 7T time-of-flight (TOF) MRI provides high resolution for the evaluation of cerebrovascular vessels and pathologies. In combination with 4D flow fields acquired with phase-contrast (PC) MRI, hemodynamic information can be extracted to enhance the analysis by providing direct measurements in the larger arteries or patient-specific boundary conditions. Hence, a registration between both modalities is required. METHODS: To combine TOF and PC-MRI data, we developed a hybrid registration approach. Vessels and their centerlines are segmented from the TOF data. The centerline is fit to the intensity ridges of the lower resolved PC-MRI data, which provides temporal information. We used a metric that utilizes a scaled sum of weighted intensities and gradients on the normal plane. The registration is then guided by decoupled local affine transformations. It is applied hierarchically following the branching order of the vessel tree. RESULTS: A landmark validation over Monte Carlo simulations yielded an average mean squared error of 184.73 mm and an average Hausdorff distance of 15.20 mm. The hierarchical traversal that transforms child vessels with their parents registers even small vessels not detectable in the PC-MRI. CONCLUSION: The presented work combines high-resolution tomographic information from 7T TOF-MRI and measured flow data from 4D 7T PC-MRI scan for the arteries of the brain. This enables usage of patient-specific flow parameters for realistic simulations, thus supporting research in areas such as cerebral small vessel disease. Automatization and free deformations can help address the limiting error measures in the future.

Quantitative evaluation of prospective motion correction in healthy subjects at 7T MRI.

PURPOSE: Quantitative assessment of prospective motion correction (PMC) capability at 7T MRI for compliant healthy subjects to improve high-resolution images in the absence of intentional motion. METHODS: Twenty-one healthy subjects were imaged at 7 T. They were asked not to move, to consider only unintentional motion. An in-bore optical tracking system was used to monitor head motion and consequently update the imaging volume. For all subjects, high-resolution T1 (3D-MPRAGE), T2 (2D turbo spin echo), proton density (2D turbo spin echo), and T2∗ (2D gradient echo) weighted images were acquired with and without PMC. The images were evaluated through subjective and objective analysis. RESULTS: Subjective evaluation overall has shown a statistically significant improvement (5.5%) in terms of image quality with PMC ON. In a separate evaluation of every contrast, three of the four contrasts (T1 , T2 , and proton density) have shown a statistically significant improvement (9.62%, 9.85%, and 9.26%), whereas the fourth one ( T2∗ ) has shown improvement, although not statistically significant. In the evaluation with objective metrics, average edge strength has shown an overall improvement of 6% with PMC ON, which was statistically significant; and gradient entropy has shown an overall improvement of 2%, which did not reach statistical significance. CONCLUSION: Based on subjective assessment, PMC improved image quality in high-resolution images of healthy compliant subjects in the absence of intentional motion for all contrasts except T2∗ , in which no significant differences were observed. Quantitative metrics showed an overall trend for an improvement with PMC, but not all differences were significant.

Hemodynamic Data Assimilation in a Subject-specific Circle of Willis Geometry.

PURPOSE: The anatomy of the circle of Willis (CoW), the brain's main arterial blood supply system, strongly differs between individuals, resulting in highly variable flow fields and intracranial vascularization patterns. To predict subject-specific hemodynamics with high certainty, we propose a data assimilation (DA) approach that merges fully 4D phase-contrast magnetic resonance imaging (PC-MRI) data with a numerical model in the form of computational fluid dynamics (CFD) simulations. METHODS: To the best of our knowledge, this study is the first to provide a transient state estimate for the three-dimensional velocity field in a subject-specific CoW geometry using DA. High-resolution velocity state estimates are obtained using the local ensemble transform Kalman filter (LETKF). RESULTS: Quantitative evaluation shows a considerable reduction (up to 90%) in the uncertainty of the velocity field state estimate after the data assimilation step. Velocity values in vessel areas that are below the resolution of the PC-MRI data (e.g., in posterior communicating arteries) are provided. Furthermore, the uncertainty of the analysis-based wall shear stress distribution is reduced by a factor of 2 for the data assimilation approach when compared to the CFD model alone. CONCLUSION: This study demonstrates the potential of data assimilation to provide detailed information on vascular flow, and to reduce the uncertainty in such estimates by combining various sources of data in a statistically appropriate fashion.

Phase-Contrast MRI Detection of Ventricular Shunt CSF Flow: Proof of Principle.

BACKGROUND AND PURPOSE: The evaluation of a suspected malfunction of a ventricular shunt is a common procedure in neurosurgery. The evaluation relies on either the interpretation of the ventricular width using cranial imaging or invasive techniques. Several attempts have been made to measure the flow velocity of cerebrospinal fluid (CSF) utilizing different phase-contrast magnet resonance imaging (PC MRI) techniques. In the present study, we evaluated 3 T (Tesla) MRI scanners for their effectiveness in determining of flow in the parenchymal portion of ventricular shunt systems with adjustable valves containing magnets. METHODS: At first, an MRI phantom was used to measure the phase-contrasts at different constant low flow rates. The next step was to measure the CSF flow in patients treated with ventricular shunts without suspected malfunction of the shunt under observation. RESULTS: The measurements of the phantom showed a linear correlation between the CSF flow and corresponding phase values. Despite many artifacts resulting from the magnetic valves, the ventricular catheter within the parenchymal portion of shunt was not superimposed by artifacts at each PC MRI plane and clearly distinguishable in 9 of 12 patients. Three patients suffering from obstructive hydrocephalus showed a clear flow signal. CONCLUSION: CSF flow detected within the parenchymal portion of the shunt by PC MRI may reliably provide information about the functional status of a ventricular shunt. Even in patients whose hydrocephalus was treated with magnetic adjustable valves, the CSF flow was detectable using PC MRI sequences at 3 T field strength.

Transient flow prediction in an idealized aneurysm geometry using data assimilation.

Hemodynamic simulations are restricted by modeling assumptions and uncertain initial and boundary conditions, whereas Phase-Contrast Magnetic Resonance Imaging (PC-MRI) data is affected by measurement noise and artifacts. To overcome the limitations of both techniques, the current study uses a Localization Ensemble Transform Kalman Filter (LETKF) to fully incorporate noisy, low-resolution Phase-Contrast MRI data into an ensemble of high-resolution numerical simulations. The analysis output provides an improved state estimate of the three-dimensional blood flow field in an intracranial aneurysm model. Benchmark measurements are carried out in a silicone phantom model of an idealized aneurysm under pulsatile inflow conditions. Validation is ensured with high-resolution Particle Imaging Velocimetry (PIV) obtained in the symmetry plane of the same geometry. Two data assimilation approaches are introduced, which differ in their way to propagate the ensemble members in time. In both cases the velocity noise is significantly reduced over the whole cardiac cycle. Quantitative and qualitative results indicate an improvement of the flow field prediction in comparison to the raw measurement data. Although biased measurement data reveal a systematic deviation from the truth, the LETKF is able to account for stochastically distributed errors. Through the implementation of the data assimilation step, physical constraints are introduced into the raw measurement data. The resulting, realistic high-resolution flow field can be readily used to assess further patient-specific parameters in addition to the velocity distribution, such as wall shear stress or pressure.

Wireless video transmission into the MRI magnet room: implementation and evaluation at 1.5T, 3T and 7T.

Purpose To analyze the interference between a wireless high definition multimedia interface (WHDMI) and magnetic resonance imaging (MRI) image quality at 1.5T, 3T and 7T. Materials and methods A wireless video transmission system (WVTS) consisting of a WHDMI and a projector was used to transmit and display a video stream into the magnet room. MR image quality was analyzed at 1.5T, 3T and 7T. Signal-to-noise-ratio (SNR¯) $(\overline {{\rm{SNR}}} )$ and radio frequency (RF)-noise spectrum were measured at three transmitter positions (A: inside the cabin, B: in front of the waveguide and C: in the control room). WVTS system functionality tests included measurements of reliability, delay and image quality. Results With the WVTS mean SNR¯ $\overline {{\rm{SNR}}} $ values significantly decreased in comparison to the reference for all positions and fieldstrenghts, while the spectra's baseline is elevated at 1.5T and 3T. Peaks related to continuous wave interferences are apparent at all field strenghts. For WHDMI alone mean SNR¯ $\overline {{\rm{SNR}}} $ values were stable without significant differences to the reference. No elevation of the spectra's baseline could be observed. Functionality measurements confirmed high connection reliability with stable image quality and no delays for all field strengths. Conclusion We conclude that wireless transmission of video streams into the MRI magnet room is feasible at all field strengths without hampering image quality.

Comparison of intracranial aneurysm flow quantification techniques: standard PIV vs stereoscopic PIV vs tomographic PIV vs phase-contrast MRI vs CFD.

Image-based hemodynamic simulations to assess the rupture risk or improve the treatment planning of intracranial aneurysms have become popular recently. However, due to strong modeling assumptions and limitations, the acceptance of numerical approaches remains limited. Therefore, validation using experimental methods is mandatory.In this study, a unique compilation of four in-vitro flow measurements (three particle image velocimetry approaches using a standard (PIV), stereoscopic (sPIV), and tomographic (tPIV) setup, as well as a phase-contrast magnetic resonance imaging (PC-MRI) measurement) were compared with a computational fluid dynamics (CFD) simulation. This was carried out in a patient-specific silicone phantom model of an internal carotid artery aneurysm under steady flow conditions. To evaluate differences between each technique, a similarity index (SI) with respect to the velocity vectors and the average velocity magnitude differences among all involved modalities were computed.The qualitative comparison reveals that all techniques are able to provide a reasonable description of the global flow structures. High quantitative agreement in terms of SI and velocity magnitude differences was found between all PIV methods and CFD. However, quantitative differences were observed between PC-MRI and the other techniques. Deeper analysis revealed that the limited resolution of the PC-MRI technique is a major contributor to the experienced differences and leads to a systematic underestimation of overall velocity magnitude levels inside the vessel. This confirms the necessity of using highly resolving flow measurement techniques, such as PIV, in an in-vitro environment to individually verify the validity of the numerically obtained hemodynamic results.

Hippocampal vascularization patterns: A high-resolution 7 Tesla time-of-flight magnetic resonance angiography study.

Considerable evidence suggests a close relationship between vascular and degenerative pathology in the human hippocampus. Due to the intrinsic fragility of its vascular network, the hippocampus appears less able to cope with hypoperfusion and anoxia than other cortical areas. Although hippocampal blood supply is generally provided by the collateral branches of the posterior cerebral artery (PCA) and the anterior choroidal artery (AChA), different vascularization patterns have been detected postmortem. To date, a methodology that enables the classification of individual hippocampal vascularization patterns in vivo has not been established. In this study, using high-resolution 7 Tesla time-of-flight angiography data (0.3 mm isotropic resolution) in young adults, we classified individual variability in hippocampal vascularization patterns involved in medial temporal lobe blood supply in vivo. A strong concordance between our classification and previous autopsy findings was found, along with interesting anatomical observations, such as the variable contribution of the AChA to hippocampal supply, the relationships between hippocampal and PCA patterns, and the different distribution patterns of the right and left hemispheres. The approach presented here for determining hippocampal vascularization patterns in vivo may provide new insights into not only the vulnerability of the hippocampus to vascular and neurodegenerative diseases but also hippocampal vascular plasticity after exercise training.

Prospective motion correction enables highest resolution time-of-flight angiography at 7T.

PURPOSE: Higher magnetic field strengths enable time-of-flight (TOF) angiography with higher resolution to depict small-vessel pathologies. However, this potential is limited by the subject's ability to remain motionless. Even small-scale, involuntary motion can degrade vessel depiction, thus limiting the effective resolution. The aim of this study was to overcome this resolution limit by deploying prospectively motion-corrected (PMC) TOF. METHODS: An optical, marker-based, in-bore tracking system was used to update the imaging volume prospectively according to the subject's head motion. PMC TOF was evaluated in 12 healthy, cooperative subjects at isotropic resolution of up to 150 µm. Image quality was assessed qualitatively through reader rating and quantitatively with the average edge-strength metric. RESULTS: PMC significantly increased the average edge strength and qualitatively improved the vessel depiction in nine out of 11 cases. Image quality was never degraded by motion correction. PMC also enabled acquisition of the highest resolution human brain in vivo TOF angiography to date. CONCLUSION: With PMC enabled, high-resolution TOF is able to visualize brain vasculature beyond the effective resolution limit. Magn Reson Med 80:248-258, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Correction of B 0-induced geometric distortion variations in prospective motion correction for 7T MRI.

OBJECTIVE: Prospective motion correction can effectively fix the imaging volume of interest. For large motion, this can lead to relative motion of coil sensitivities, distortions associated with imaging gradients and B 0 field variations. This work accounts for the B 0 field change due to subject movement, and proposes a method for correcting tissue magnetic susceptibility-related distortion in prospective motion correction. MATERIALS AND METHODS: The B 0 field shifts at the different head orientations were characterized. A volunteer performed large motion with prospective motion correction enabled. The acquired data were divided into multiple groups according to the object positions. The correction of B 0-related distortion was applied to each group of data individually via augmented sensitivity encoding with additionally integrated gradient nonlinearity correction. RESULTS: The relative motion of the gradients, B 0 field and coil sensitivities in prospective motion correction results in residual spatial distortion, blurring, and coil artifacts. These errors can be mitigated by the proposed method. Moreover, iterative conjugate gradient optimization with regularization provided superior results with smaller RMSE in comparison to standard conjugate gradient. CONCLUSION: The combined correction of B 0-related distortion and gradient nonlinearity leads to a reduction of residual motion artifacts in prospective motion correction data.

Highest Resolution In Vivo Human Brain MRI Using Prospective Motion Correction.

High field MRI systems, such as 7 Tesla (T) scanners, can deliver higher signal to noise ratio (SNR) than lower field scanners and thus allow for the acquisition of data with higher spatial resolution, which is often demanded by users in the fields of clinical and neuroscientific imaging. However, high resolution scans may require long acquisition times, which in turn increase the discomfort for the subject and the risk of subject motion. Even with a cooperative and trained subject, involuntary motion due to heartbeat, swallowing, respiration and changes in muscle tone can cause image artifacts that reduce the effective resolution. In addition, scanning with higher resolution leads to increased sensitivity to even very small movements. Prospective motion correction (PMC) at 3T and 7T has proven to increase image quality in case of subject motion. Although the application of prospective motion correction is becoming more popular, previous articles focused on proof of concept studies and technical descriptions, whereas this paper briefly describes the technical aspects of the optical tracking system, marker fixation and cross calibration and focuses on the application of PMC to very high resolution imaging without intentional motion. In this study we acquired in vivo MR images at 7T using prospective motion correction during long acquisitions. As a result, we present images among the highest, if not the highest resolution of in vivo human brain MRI ever acquired.

Correction of gradient nonlinearity artifacts in prospective motion correction for 7T MRI.

PURPOSE: To demonstrate the effect of gradient nonlinearity and develop a method for correction of gradient nonlinearity artifacts in prospective motion correction (Mo-Co). METHODS: Nonlinear gradients can induce geometric distortions in magnetic resonance imaging, leading to pixel shifts with errors of up to several millimeters, thereby interfering with precise localization of anatomical structures. Prospective Mo-Co has been extended by conventional gradient warp correction applied to individual phase encoding steps/groups during the reconstruction. The gradient-related displacements are approximated using spherical harmonic functions. In addition, the combination of this method with a retrospective correction of the changes in the coil sensitivity profiles relative to the object (augmented sensitivity encoding (SENSE) reconstruction) was evaluated in simulation and experimental data. RESULTS: Prospective Mo-Co under gradient fields and coils sensitivity inconsistencies results in residual blurring, spatial distortion, and coil sensitivity mismatch artifacts. These errors can be considerably mitigated by the proposed method. High image quality with very little remaining artifacts was achieved after a few iterations. The relative image errors decreased from 25.7% to below 17.3% after 10 iterations. CONCLUSION: The combined correction of gradient nonlinearity and sensitivity map variation leads to a pronounced reduction of residual motion artifacts in prospectively motion-corrected data.

Cerebral blood flow in a healthy Circle of Willis and two intracranial aneurysms: computational fluid dynamics versus four-dimensional phase-contrast magnetic resonance imaging.

Computational fluid dynamics (CFD) opens up multiple opportunities to investigate the hemodynamics of the human vascular system. However, due to numerous assumptions the acceptance of CFD among physicians is still limited in practice and validation through comparison is mandatory. Time-dependent quantitative phase-contrast magnetic resonance imaging PC-MRI measurements in a healthy volunteer and two intracranial aneurysms were carried out at 3 and 7 Tesla. Based on the acquired images, three-dimensional (3D) models of the aneurysms were reconstructed and used for the numerical simulations. Flow information from the MR measurements were applied as boundary conditions. The four-dimensional (4D) velocity fields obtained by CFD and MRI were qualitatively as well as quantitatively compared including cut planes and vector analyses. For all cases a high similarity of the velocity patterns was observed. Additionally, the quantitative analysis revealed a good agreement between CFD and MRI. Deviations were caused by minor differences between the reconstructed vessel models and the actual lumen. The comparisons between diastole and systole indicate that relative differences between MRI and CFD are intensified with increasing velocity. The findings of this study lead to the conclusion that CFD and MRI agree well in predicting intracranial velocities when realistic geometries and boundary conditions are provided. Due to the considerably higher temporal and spatial resolution of CFD compared to MRI, complex flow patterns can be further investigated in order to evaluate their role with respect to aneurysm formation or rupture. Nevertheless, special care is required regarding the vessel reconstruction since the geometry has a major impact on the subsequent numerical results.

Measurement and correction of microscopic head motion during magnetic resonance imaging of the brain.

Magnetic resonance imaging (MRI) is a widely used method for non-invasive study of the structure and function of the human brain. Increasing magnetic field strengths enable higher resolution imaging; however, long scan times and high motion sensitivity mean that image quality is often limited by the involuntary motion of the subject. Prospective motion correction is a technique that addresses this problem by tracking head motion and continuously updating the imaging pulse sequence, locking the imaging volume position and orientation relative to the moving brain. The accuracy and precision of current MR-compatible tracking systems and navigator methods allows the quantification and correction of large-scale motion, but not the correction of very small involuntary movements in six degrees of freedom. In this work, we present an MR-compatible tracking system comprising a single camera and a single 15 mm marker that provides tracking precision in the order of 10 m and 0.01 degrees. We show preliminary results, which indicate that when used for prospective motion correction, the system enables improvement in image quality at both 3 T and 7 T, even in experienced and cooperative subjects trained to remain motionless during imaging. We also report direct observation and quantification of the mechanical ballistocardiogram (BCG) during simultaneous MR imaging. This is particularly apparent in the head-feet direction, with a peak-to-peak displacement of 140 m.

Navigator accuracy requirements for prospective motion correction.

Prospective motion correction in MRI is becoming increasingly popular to prevent the image artifacts that result from subject motion. Navigator information is used to update the position of the imaging volume before every spin excitation so that lines of acquired k-space data are consistent. Errors in the navigator information, however, result in residual errors in each k-space line. This paper presents an analysis linking noise in the tracking system to the power of the resulting image artifacts. An expression is formulated for the required navigator accuracy based on the properties of the imaged object and the desired resolution. Analytical results are compared with computer simulations and experimental data.