Diffusion-weighted MRI can assist preoperative planning by reconstructing the trajectory of eloquent fiber pathways, such as the corticospinal tract (CST). However, accurate reconstruction of the full extent of the CST remains challenging with existing tractography methods. We suggest a novel tractography algorithm exploiting unused fiber orientations to produce more complete and reliable results.
MethodsOur novel approach, referred to as multi-level fiber tractography (MLFT), reconstructs fiber pathways by progressively considering previously unused fiber orientations at multiple levels of tract propagation. Anatomical priors are used to minimize the number of false-positive pathways. The MLFT method was evaluated on synthetic data and in vivo data by reconstructing the CST while compared to conventional tractography approaches.
ResultsThe radial extent of MLFT reconstructions is comparable to that of probabilistic reconstruction: \(p=0.21\) for the left and \(p=0.53\) for the right hemisphere according to Wilcoxon test, while achieving significantly higher topography preservation compared to probabilistic tractography: \(p<0.01\).
DiscussionMLFT provides a novel way to reconstruct fiber pathways by adding the capability of including branching pathways in fiber tractography. Thanks to its robustness, feasible reconstruction extent and topography preservation, our approach may assist in clinical practice as well as in virtual dissection studies.
相似文献Amide proton transfer (APT) weighted chemical exchange saturation transfer (CEST) imaging is increasingly used to investigate high-grade, enhancing brain tumours. Non-enhancing glioma is currently less studied, but shows heterogeneous pathophysiology with subtypes having equally poor prognosis as enhancing glioma. Here, we investigate the use of CEST MRI to best differentiate non-enhancing glioma from healthy tissue and image tumour heterogeneity.
Materials & MethodsA 3D pulsed CEST sequence was applied at 3 Tesla with whole tumour coverage and 31 off-resonance frequencies (+6 to -6 ppm) in 18 patients with non-enhancing glioma. Magnetisation transfer ratio asymmetry (MTRasym) and Lorentzian difference (LD) maps at 3.5 ppm were compared for differentiation of tumour versus normal appearing white matter. Heterogeneity was mapped by calculating volume percentages of the tumour showing hyperintense APT-weighted signal.
ResultsLDamide gave greater effect sizes than MTRasym to differentiate non-enhancing glioma from normal appearing white matter. On average, 17.9 % ± 13.3 % (min–max: 2.4 %–54.5 %) of the tumour volume showed hyperintense LDamide in non-enhancing glioma.
ConclusionThis works illustrates the need for whole tumour coverage to investigate heterogeneity in increased APT-weighted CEST signal in non-enhancing glioma. Future work should investigate whether targeting hyperintense LDamide regions for biopsies improves diagnosis of non-enhancing glioma.
相似文献A method for Orthogonal Phase Encoding Reduction of Artifact (OPERA) was developed and tested.
Materials and methodsBecause the position of ghosts and aliasing artifacts is predictable along columns or rows, OPERA combines the intensity values of two images acquired using the same parameters, but with swapped phase-encoding directions, to correct the artifacts. Simulations and phantom experiments were conducted to define the efficacy, robustness, and reproducibility. Clinical validation was performed on a total of 1003 images by comparing the OPERA-corrected images and the corresponding image standard in terms of Signal-to-Noise Ratio (SNR) and Contrast-to-Noise Ratio (CNR). The method efficacy was also rated using a Likert-type scale response by two experienced independent radiologists using a single-blinded procedure.
ResultsSimulations and phantom experiments demonstrated the robustness and effectiveness of OPERA in reducing artifacts strength. OPERA application did not significantly change the SNR [+?4.16%; inter-quartile range (IQR): 2.72–5.01%] and CNR (+?4.30%; IQR: 2.86–6.04%) values. The two radiologists observed a total of 893 original images with artifacts (89.03% of the total images), a reduction in the perceived artifacts of 82.0% and 83.9% (p?<?0.0001), and an improvement in the perceived SNR (82.8% and 88.5%; K?=?0.714) and perceived CNR (86.9–88.9%; K?=?0.722).
DiscussionThe study demonstrated that OPERA reduces MR artifacts and improves the perceived image quality.
相似文献To provide respiratory motion correction for free-breathing myocardial T1 mapping using a pilot tone (PT) and a continuous golden-angle radial acquisition.
Materials and methodsDuring a 45 s prescan the PT is acquired together with a dynamic sagittal image covering multiple respiratory cycles. From these images, the respiratory heart motion in head-feet and anterior–posterior direction is estimated and two linear models are derived between the PT and heart motion. In the following scan through-plane motion is corrected prospectively with slice tracking based on the PT. In-plane motion is corrected for retrospectively. Our method was evaluated on a motion phantom and 11 healthy subjects.
ResultsNon-motion corrected measurements using a moving phantom showed T1 errors of 14 ± 4% (p < 0.05) compared to a reference measurement. The proposed motion correction approach reduced this error to 3 ± 4% (p < 0.05). In vivo the respiratory motion led to an overestimation of T1 values by 26 ± 31% compared to breathhold T1 maps, which was successfully corrected to an average difference of 3 ± 2% (p < 0.05) between our free-breathing approach and breathhold data.
DiscussionOur proposed PT-based motion correction approach allows for T1 mapping during free-breathing with the same accuracy as a corresponding breathhold T1 mapping scan.
相似文献The aim of this study was to develop and validate an MRI protocol based on a variable echo time (vTE) sensitive to the short T2* components of the sciatic nerve.
Materials and methods15 healthy subjects (M/F: 9/6; age: 21–62) were scanned at 3T targeting the sciatic nerve at the thigh bilaterally, using a dual echo variable echo time (vTE) sequence (based on a spoiled gradient echo acquisition) with echo times of 0.98/5.37 ms. Apparent T2* (aT2*) values of the sciatic nerves were calculated with a mono-exponential fit and used for data comparison.
ResultsThere were no significant differences in aT2* related to side, sex, age, and BMI, even though small differences for side were reported. Good-to-excellent repeatability and reproducibility were found for geometry of ROIs (Dice indices: intra-rater 0.68–0.7; inter-rater 0.70–0.72) and the related aT2* measures (intra-inter reader ICC 0.95–0.97; 0.66–0.85) from two different operators. Side-related signal-to-noise-ratio non-significant differences were reported, while contrast-to-noise-ratio measures were excellent both for side and echo.
DiscussionOur study introduces a novel MR sequence sensitive to the short T2* components of the sciatic nerve and may be used for the study of peripheral nerve disorders.
相似文献In this study we address the automatic segmentation of selected muscles of the thigh and leg through a supervised deep learning approach.
Material and methodsThe application of quantitative imaging in neuromuscular diseases requires the availability of regions of interest (ROI) drawn on muscles to extract quantitative parameters. Up to now, manual drawing of ROIs has been considered the gold standard in clinical studies, with no clear and universally accepted standardized procedure for segmentation. Several automatic methods, based mainly on machine learning and deep learning algorithms, have recently been proposed to discriminate between skeletal muscle, bone, subcutaneous and intermuscular adipose tissue. We develop a supervised deep learning approach based on a unified framework for ROI segmentation.
ResultsThe proposed network generates segmentation maps with high accuracy, consisting in Dice Scores ranging from 0.89 to 0.95, with respect to “ground truth” manually segmented labelled images, also showing high average performance in both mild and severe cases of disease involvement (i.e. entity of fatty replacement).
DiscussionThe presented results are promising and potentially translatable to different skeletal muscle groups and other MRI sequences with different contrast and resolution.
相似文献To improve the precision of a free-breathing 3D saturation-recovery-based myocardial T1 mapping sequence using a post-processing 3D denoising technique.
MethodsA T1 phantom and 15 healthy subjects were scanned on a 1.5 T MRI scanner using 3D saturation-recovery single-shot acquisition (SASHA) for myocardial T1 mapping. A 3D denoising technique was applied to the native T1-weighted images before pixel-wise T1 fitting. The denoising technique imposes edge-preserving regularity and exploits the co-occurrence of 3D spatial gradients in the native T1-weighted images by incorporating a multi-contrast Beltrami regularization. Additionally, 2D modified Look-Locker inversion recovery (MOLLI) acquisitions were performed for comparison purposes. Accuracy and precision were measured in the myocardial septum of 2D MOLLI and 3D SASHA T1 maps and then compared. Furthermore, the accuracy and precision of the proposed approach were evaluated in a standardized phantom in comparison to an inversion-recovery spin-echo sequence (IRSE).
ResultsFor the phantom study, Bland–Altman plots showed good agreement in terms of accuracy between IRSE and 3D SASHA, both on non-denoised and denoised T1 maps (mean difference −1.4 ± 18.9 ms and −4.4 ± 21.2 ms, respectively), while 2D MOLLI generally underestimated the T1 values (69.4 ± 48.4 ms). For the in vivo study, there was a statistical difference between the precision measured on 2D MOLLI and on non-denoised 3D SASHA T1 maps (P = 0.005), while there was no statistical difference after denoising (P = 0.95).
ConclusionThe precision of 3D SASHA myocardial T1 mapping was substantially improved using a 3D Beltrami regularization based denoising technique and was similar to that of 2D MOLLI T1 mapping, while preserving the higher accuracy and whole-heart coverage of 3D SASHA.
相似文献Object
We propose a new tracking method based on time-of-arrival (TOA) maps derived from simulated diffusion processes.Materials and methods
The proposed diffusion simulation-based tracking consists of three steps that are successively evaluated on small overlapping sub-regions in a diffusion tensor field. First, the diffusion process is simulated for several time steps. Second, a TOA map is created to store simulation results for the individual time steps that are required for the tract reconstruction. Third, the fiber pathway is reconstructed on the TOA map and concatenated between neighboring sub-regions. This new approach is compared with probabilistic and streamline tracking. All methods are applied to synthetic phantom data for an easier evaluation of their fiber reconstruction quality.Results
The comparison of the tracking results did show severe problems for the streamline approach in the reconstruction of crossing fibers, for example. The probabilistic method was able to resolve the crossing, but could not handle strong curvature. The new diffusion simulation-based tracking could reconstruct all problematic fiber constellations.Conclusion
The proposed diffusion simulation-based tracking method used the whole tensor information of a neighborhood of voxels and is, therefore, able to handle problematic tracking situations better than established methods. 相似文献With the hypothesis that 4D flow can be used in evaluation of cardiac shunts, we seek to evaluate the multilevel and interreader reproducibility of measurements of the blood flow, shunt fraction and shunt volume in patients with atrial septum defect (ASD) in practice at multiple clinical sites.
Materials and methodsFour-dimensional flow MRI examinations were performed at four institutions across Europe and the US. Twenty-nine patients (mean age, 43 years; 11 male) were included in the study. Flow measurements were performed at three levels (valve, main artery and periphery) in both the pulmonary and systemic circulation by two independent readers and compared against stroke volumes from 4D flow anatomic data. Further, the shunt ratio (Qp/Qs) was calculated. Additionally, shunt volume was quantified at the atrial level by tracking the atrial septum.
ResultsMeasurements of the pulmonary blood flow at multiple levels correlate well whether measuring at the valve, main pulmonary artery or branch pulmonary arteries (r = 0.885–0.886). Measurements of the systemic blood flow show excellent correlation, whether measuring at the valve, ascending aorta or sum of flow from the superior vena cava (SVC) and descending aorta (r = 0.974–0.991). Intraclass agreement between the two observers for the flow measurements varies between 0.96 and 0.99. Compared with stroke volume, pulmonic flow is underestimated with 0.26 l/min at the main pulmonary artery level, and systemic flow is overestimated with 0.16 l/min at the ascending aorta level. Direct measurements of ASD flow are feasible in 20 of 29 (69%) patients.
ConclusionBlood flow and shunt quantification measured at multiple levels and performed by different readers are reproducible and consistent with 4D flow MRI.
相似文献To implement magnetic resonance fingerprinting (MRF) on a permanent magnet 50 mT low-field system deployable as a future point-of-care (POC) unit and explore the quality of the parameter maps.
Materials and methods3D MRF was implemented on a custom-built Halbach array using a slab-selective spoiled steady-state free precession sequence with 3D Cartesian readout. Undersampled scans were acquired with different MRF flip angle patterns and reconstructed using matrix completion and matched to the simulated dictionary, taking excitation profile and coil ringing into account. MRF relaxation times were compared to that of inversion recovery (IR) and multi-echo spin echo (MESE) experiments in phantom and in vivo. Furthermore, B0 inhomogeneities were encoded in the MRF sequence using an alternating TE pattern, and the estimated map was used to correct for image distortions in the MRF images using a model-based reconstruction.
ResultsPhantom relaxation times measured with an optimized MRF sequence for low field were in better agreement with reference techniques than for a standard MRF sequence. In vivo muscle relaxation times measured with MRF were longer than those obtained with an IR sequence (T1: 182 ± 21.5 vs 168 ± 9.89 ms) and with an MESE sequence (T2: 69.8 ± 19.7 vs 46.1 ± 9.65 ms). In vivo lipid MRF relaxation times were also longer compared with IR (T1: 165 ± 15.1 ms vs 127 ± 8.28 ms) and with MESE (T2: 160 ± 15.0 ms vs 124 ± 4.27 ms). Integrated ΔB0 estimation and correction resulted in parameter maps with reduced distortions.
DiscussionIt is possible to measure volumetric relaxation times with MRF at 2.5 × 2.5 × 3.0 mm3 resolution in a 13 min scan time on a 50 mT permanent magnet system. The measured MRF relaxation times are longer compared to those measured with reference techniques, especially for T2. This discrepancy can potentially be addressed by hardware, reconstruction and sequence design, but long-term reproducibility needs to be further improved.
相似文献To improve accelerated MRI reconstruction through a densely connected cascading deep learning reconstruction framework.
Materials and methodsA cascading deep learning reconstruction framework (reference model) was modified by applying three architectural modifications: input-level dense connections between cascade inputs and outputs, an improved deep learning sub-network, and long-range skip-connections between subsequent deep learning networks. An ablation study was performed, where five model configurations were trained on the NYU fastMRI neuro dataset with an end-to-end scheme conjunct on four- and eightfold acceleration. The trained models were evaluated by comparing their respective structural similarity index measure (SSIM), normalized mean square error (NMSE), and peak signal to noise ratio (PSNR).
ResultsThe proposed densely interconnected residual cascading network (DIRCN), utilizing all three suggested modifications achieved a SSIM improvement of 8% and 11%, a NMSE improvement of 14% and 23%, and a PSNR improvement of 2% and 3% for four- and eightfold acceleration, respectively. In an ablation study, the individual architectural modifications all contributed to this improvement for both acceleration factors, by improving the SSIM, NMSE, and PSNR with approximately 2–4%, 4–9%, and 0.5–1%, respectively.
ConclusionThe proposed architectural modifications allow for simple adjustments on an already existing cascading framework to further improve the resulting reconstructions.
相似文献Fetal brain diffusion tensor imaging (DTI) offers quantitative analysis of the developing brain. The objective was to 1) quantify DTI measures across gestation in a cohort of fetuses without brain abnormalities using full retrospective correction for fetal head motion 2) compare results obtained in utero to those in preterm infants.
Materials and methodsMotion-corrected DTI analysis was performed on data sets obtained at 1.5T from 32 fetuses scanned between 21.29 and 37.57 (median 31.86) weeks. Results were compared to 32 preterm infants scanned at 3T between 27.43 and 37.14 (median 33.07) weeks. Apparent diffusion coefficient (ADC) and fractional anisotropy (FA) were quantified by region of interest measurements and tractography was performed.
ResultsFetal DTI was successful in 84% of fetuses for whom there was sufficient data for DTI estimation, and at least one tract could be obtained in 25 cases. Fetal FA values increased and ADC values decreased with age at scan (PLIC FA: p = 0.001; R2 = 0.469; slope = 0.011; splenium FA: p < 0.001; R2 = 0.597; slope = 0.019; thalamus ADC: p = 0.001; R2 = 0.420; slope = − 0.023); similar trends were found in preterm infants.
ConclusionThis study demonstrates that stable DTI is feasible on fetuses and provides evidence for normative values of diffusion properties that are consistent with aged matched preterm infants.
相似文献To perform a systematic review of the literature exploring magnetic resonance imaging (MRI) methods for measuring natural brain tissue pulsations (BTPs) in humans.
MethodsA prospective systematic search of MEDLINE, SCOPUS and OpenGrey databases was conducted by two independent reviewers using a pre-determined strategy. The search focused on identifying reported measurements of naturally occurring BTP motion in humans. Studies involving non-human participants, MRI in combination with other modalities, MRI during invasive procedures and MRI studies involving externally applied tests were excluded. Data from the retrieved records were combined to create Forest plots comparing brain tissue displacement between Chiari-malformation type 1 (CM-I) patients and healthy controls using an independent samples t-test.
ResultsThe search retrieved 22 eligible articles. Articles described 5 main MRI techniques for visualisation or quantification of intrinsic brain motion. MRI techniques generally agreed that the amplitude of BTPs varies regionally from 0.04 mm to ~ 0.80 mm, with larger tissue displacements occurring closer to the centre and base of the brain compared to peripheral regions. Studies of brain pathology using MRI BTP measurements are currently limited to tumour characterisation, idiopathic intracranial hypertension (IIH), and CM-I. A pooled analysis confirmed that displacement of tissue in the cerebellar tonsillar region of CM-I patients was + 0.31 mm [95% CI 0.23, 0.38, p < 0.0001] higher than in healthy controls.
DiscussionMRI techniques used for measurements of brain motion are at an early stage of development with high heterogeneity across the methods used. Further work is required to provide normative data to support systematic BTPs characterisation in health and disease.
相似文献MRI of excised hearts at ultra-high field strengths (\({\mathrm{B}}_{0}\)≥7 T) can provide high-resolution, high-fidelity ground truth data for biomedical studies, imaging science, and artificial intelligence. In this study, we demonstrate the capabilities of a custom-built, multiple-element transceiver array customized for high-resolution imaging of excised hearts.
MethodA dedicated 16-element transceiver loop array was implemented for operation in parallel transmit (pTx) mode (8Tx/16Rx) of a clinical whole-body 7 T MRI system. The initial adjustment of the array was performed using full-wave 3D-electromagnetic simulation with subsequent final fine-tuning on the bench.
ResultsWe report the results of testing the implemented array in tissue-mimicking liquid phantoms and excised porcine hearts. The array demonstrated high efficiency of parallel transmits characteristics enabling efficient pTX-based B1+-shimming.
ConclusionThe receive sensitivity and parallel imaging capability of the dedicated coil were superior to that of a commercial 1Tx/32Rx head coil in both SNR and T2*-mapping. The array was successfully tested to acquire ultra-high-resolution (0.1 × 0.1 × 0.8 mm voxel) images of post-infarction scar tissue. High-resolution (isotropic 1.6 mm3 voxel) diffusion tensor imaging-based tractography provided high-resolution information about normal myocardial fiber orientation.
相似文献The aim of this study was to compare the use of high-resolution phase and QSM images acquired at ultra-high field in the investigation of multiple sclerosis (MS) lesions with peripheral rings, and to discuss their usefulness for drawing inferences about underlying tissue composition.
Materials and methodsThirty-nine Subjects were scanned at 7 T, using 3D T 2*-weighted and T 1-weighted sequences. Phase images were then unwrapped and filtered, and quantitative susceptibility maps were generated using a thresholded k-space division method. Lesions were compared visually and using a 1D profiling algorithm.
ResultsLesions displaying peripheral rings in the phase images were identified in 10 of the 39 subjects. Dipolar projections were apparent in the phase images outside of the extent of several of these lesions; however, QSM images showed peripheral rings without such projections. These projections appeared ring-like in a small number of phase images where no ring was observed in QSM. 1D profiles of six well-isolated example lesions showed that QSM contrast corresponds more closely to the magnitude images than phase contrast.
ConclusionsPhase images contain dipolar projections, which confounds their use in the investigation of tissue composition in MS lesions. Quantitative susceptibility maps correct these projections, providing insight into the composition of MS lesions showing peripheral rings.
相似文献To evaluate the feasibility of in-vivo quantitative susceptibility mapping (QSM) of the human kidney.
MethodsAn axial single-breath-hold 3D multi-echo sequence (acquisition time 33 s) was completed on a 3 T-MRI-scanner (Magnetom Prisma, Siemens Healthineers, Erlangen, Germany) in 19 healthy volunteers. Graph-cut-based unwrapping combined with the T2*-IDEAL approach was performed to remove the chemical shift of fat and to quantify QSM of the upper abdomen. Mean susceptibility values of the entire, renal cortex and medulla in both kidneys and the liver were determined and compared. Five subjects were measured twice to examine the reproducibility. One patient with severe renal fibrosis was included in the study to evaluate the potential clinical relevance of QSM.
ResultsQSM was successful in 17 volunteers and the patient with renal fibrosis. Anatomical structures in the abdomen were clearly distinguishable by QSM and the susceptibility values obtained in the liver were comparable to those found in the literature. The results showed a good reproducibility. Besides, the mean renal QSM values obtained in healthy volunteers (0.04?±?0.07 ppm for the right and ? 0.06?±?0.19 ppm for the left kidney) were substantially higher than that measured in the investigated fibrotic kidney (? 0.43?±?? 0.02 ppm).
ConclusionQSM of the human kidney could be a promising approach for the assessment of information about microscopic renal tissue structure. Therefore, it might further improve functional renal MR imaging.
相似文献Neonatal brain and cardiac imaging would benefit from the increased signal-to-noise ratio levels at 7 T compared to lower field. Optimal performance might be achieved using purpose designed RF coil arrays. In this study, we introduce an 8-channel dipole array and investigate, using simulations, its RF performances for neonatal applications at 7 T.
MethodsThe 8-channel dipole array was designed and evaluated for neonatal brain/cardiac configurations in terms of SAR efficiency (ratio between transmit-field and maximum specific-absorption-rate level) using adjusted dielectric properties for neonate. A birdcage coil operating in circularly polarized mode was simulated for comparison. Validation of the simulation model was performed on phantom for the coil array.
ResultsThe 8-channel dipole array demonstrated up to 46% higher SAR efficiency levels compared to the birdcage coil in neonatal configurations, as the specific-absorption-rate levels were alleviated. An averaged normalized root-mean-square-error of 6.7% was found between measured and simulated transmit field maps on phantom.
ConclusionThe 8-channel dipole array design integrated for neonatal brain and cardiac MR was successfully demonstrated, in simulation with coverage of the baby and increased SAR efficiency levels compared to the birdcage. We conclude that the 8Tx-dipole array promises safe operating procedures for MR imaging of neonatal brain and heart at 7 T.
相似文献To refine a new technique to measure respiratory-resolved left ventricular end-diastolic volume (LVEDV) in mid-inspiration and mid-expiration using a respiratory self-gating technique and demonstrate clinical feasibility in patients.
Materials and methodsTen consecutive patients were imaged at 1.5 T during 10 min of free breathing using a 3D golden-angle radial trajectory. Two respiratory self-gating signals were extracted and compared: from the k-space center of all acquired spokes, and from a superior–inferior projection spoke repeated every 64 ms. Data were binned into end-diastole and two respiratory phases of 15% respiratory cycle duration in mid-inspiration and mid-expiration. LVED volume and septal–lateral diameter were measured from manual segmentation of the endocardial border.
ResultsRespiratory-induced variation in LVED size expressed as mid-inspiration relative to mid-expiration was, for volume, 1 ± 8% with k-space-based self-gating and 8 ± 2% with projection-based self-gating (P = 0.04), and for septal–lateral diameter, 2 ± 2% with k-space-based self-gating and 10 ± 1% with projection-based self-gating (P = 0.002).
DiscussionMeasuring respiratory variation in LVED size was possible in clinical patients with projection-based respiratory self-gating, and the measured respiratory variation was consistent with previous studies on healthy volunteers. Projection-based self-gating detected a higher variation in LVED volume and diameter during respiration, compared to k-space-based self-gating.
相似文献