Image reconstruction in k-space from MR data encoded with ambiguous gradient fields.

PURPOSE: In this work, the limits of image reconstruction in k-space are explored when non-bijective gradient fields are used for spatial encoding. THEORY: The image space analogy between parallel imaging and imaging with non-bijective encoding fields is partially broken in k-space. As a consequence, it is hypothesized and proven that ambiguities can only be resolved partially in k-space, and not completely as is the case in image space. METHODS: Image-space and k-space based reconstruction algorithms for multi-channel radiofrequency data acquisitions are programmed and tested using numerical simulations as well as in vivo measurement data. RESULTS: The hypothesis is verified based on an analysis of reconstructed images. It is found that non-bijective gradient fields have the effect that densely sampled autocalibration data, used for k-space reconstruction, provide less information than a separate scan of the receiver coil sensitivity maps, used for image space reconstruction. Consequently, in k-space only the undersampling artifact can be unfolded, whereas in image space, it is also possible to resolve aliasing that is caused by the non-bijectivity of the gradient fields. CONCLUSION: For standard imaging, reconstruction in image space and in k-space is nearly equivalent, whereas there is a fundamental difference with practical consequences for the selection of image reconstruction algorithms when non-bijective encoding fields are involved.

Whole-body MRI-based fat quantification: a comparison to air displacement plethysmography.

PURPOSE: To demonstrate the feasibility of an algorithm for MRI whole-body quantification of internal and subcutaneous fat and quantitative comparison of total adipose tissue to air displacement plethysmography (ADP). MATERIALS AND METHODS: For comparison with ADP, whole-body MR data of 11 volunteers were obtained using a continuously moving table Dixon sequence. Resulting fat images were corrected for B1 related intensity inhomogeneities before fat segmentation. RESULTS: The performed MR measurements of the whole body provided a direct comparison to ADP measurements. The segmentation of subcutaneous and internal fat in the abdomen worked reliably with an accuracy of 98%. Depending on the underlying model for fat quantification, the resultant MR fat masses represent an upper and a lower limit for the true fat masses. In comparison to ADP, the results were in good agreement with ρ ≥ 0.97, P < 0.0001. CONCLUSION: Whole-body fat quantities derived noninvasively by using a continuously moving table Dixon acquisition were directly compared with ADP. The accuracy of the method and the high reproducibility of results indicate its potential for clinical applications.

Temporally constrained respiratory gating improves continuously moving table MRI during free breathing.

PURPOSE: To evaluate a novel breathing motion correction algorithm for continuously moving table magnetic resonance imaging (CMT-MRI) that optimizes motion consistency in a fixed time span. MATERIALS AND METHODS: In 22 patients CMT-MRI was performed during free breathing. During a preparatory phase (constant) or continuously during the scan (adaptive) gating thresholds were computed from breathing states that should allow for motion consistent k-space sampling. After data from a first k-space traversal was acquired irrespective of breathing motion, subsequently k-space lines with discordant breathing states were reacquired below the gating threshold. Time constraints of CMT-MRI were respected, because a fixed time was allocated for reacquisition. Image quality and lesion depiction were evaluated on images reconstructed from the first traversal and motion-corrected images. RESULTS: Compared to constant thresholds, gating with adaptive thresholds led to a higher number of reacquired k-space lines (60.1%/41.7%) and a larger fraction of motion consistent final k-space lines (96.6%/78.8%). Adaptive gating induced a significant increase in image quality for all regions affected by breathing motion. Only one of 22 lesions was not depicted on the adaptively corrected images, whereas 15 were readily appreciable. CONCLUSION: Temporally constrained respiratory gating with adaptive thresholds allows for fully sampled, motion-corrected CMT-MRI acquisitions during free breathing.

Tendon graft fixation sites at the coracoid process for reconstruction of the coracoclavicular ligaments: a kinematic evaluation of three different surgical techniques.

PURPOSE: The virtual graft length kinematics of 3 operative techniques were investigated and compared with kinematics of the native coracoclavicular ligaments. METHODS: Thirteen healthy volunteers underwent magnetic resonance imaging (MRI) of the shoulder in 30° increments of abduction (0° to 120°). A 3-dimensional model of the coracoid process (CP) and the clavicle (CL) was created. Footprints of the conoid and the trapezoid ligament were identified. At the CP the potential fixation sites of 3 techniques for reconstruction of the coracoclavicular ligaments (CCLs) were marked. The techniques investigated were (1) horizontal transcoracoid drilling (TH), (2) transclavicular-transcoracoid drilling (TT), and (3) tendon graft passage underneath (PU) the coracoid process. Distances between the clavicular and coracoidal footprints of the coracoclavicular ligaments and to the virtual footprints on the coracoid process were determined for each abduction increment. RESULTS: All methods investigated resulted in a significantly longer virtual trapezoidal graft (P = .001). In PU, in addition, the virtual conoidal graft was significantly longer. TT resulted in a virtual conoidal graft and conoid ligament of equal length. TH showed identical length and distance regulation of the virtual conoidal graft and the conoid ligament, but significant shortening of the virtual trapezoidal graft during abduction. PU showed isometry of the virtual trapezoidal and conoidal grafts. CONCLUSIONS: None of the described procedures for graft fixation restores the kinematics of the native coracoclavicular ligaments. Graft fixation techniques should be chosen with respect to the preoperative type of instability. Persisting isolated vertical instability might benefit from fixation of the conoidal grafts at the native clavicular footprint. For horizontal clavicular instabilities, techniques more preserving of trapezoid ligament kinematics might be favorable. CLINICAL RELEVANCE: The data suggest that the technique of fixation in conoid and trapezoid ligament reconstruction should depend on the underlying type of instability.

Compensation of breathing motion artifacts for MRI with continuously moving table.

Breathing motion is one of the main sources of artifacts in MRI acquisitions that can severely impair diagnosis. In MRI with continuously moving table, the application of common motion compensation approaches such as breath holding or the synchronization of the measurement with the breathing motion can be problematic. In this study, a technique for the reduction of breathing-motion artifacts for MRI with continuously moving table is presented, which reconstructs motion-consistent volumes from data acquired during free breathing. Axial images are acquired rapidly compared to the period of the breathing motion and consistently combined using a combination of rigid and nonrigid slice-to-volume registration. This new technique is compared to a previously reported artifact reduction method for MRI with continuously moving table that is based on the same acquisition scheme. While the latter method only suppresses ghosting artifacts, the new technique is shown to additionally reduce blurring, misregistrations, and signal cancellations in the reconstructed images.

Continuously moving table time-of-flight angiography of the peripheral veins.

Time-of-flight (TOF) MR angiography allows for noninvasive vessel imaging. To overcome the limited volumetric coverage of standard TOF techniques, the aim of this study was to investigate the combination of TOF and continuously moving table (CMT) acquisitions for peripheral vein imaging based on image subtraction. Two acquisition strategies are presented: a simple two-step method based on 2-fold CMT acquisition and an advanced one-step method requiring only one continuous scan. Image quality of both CMT TOF techniques was evaluated by semiquantitative image grading and by signal-to-noise ratio and contrast-to-noise ratio analysis for veins of the upper and lower leg in 10 healthy volunteers. Results were compared to a standard stationary two-dimensional (2D) TOF multistation acquisition. Image grading revealed good image quality for both CMT TOF methods, thereby confirming the feasibility of axial 2D CMT TOF to assess the veins of the lower extremities during a single scan. Quantitative evaluation showed no significant difference in signal-to-noise ratio and contrast-to-noise ratio compared to the stationary experiment. Additional measurements in three patients with postthrombotic changes and varicosities demonstrated the clinical applicability of the presented methods. CMT TOF provides promising results and permits the detection of various pathologic changes of the venous system.

Increasing efficiency of parallel imaging for 2D multislice acquisitions.

Parallel imaging algorithms require precise knowledge about the spatial sensitivity variation of the receiver coils to reconstruct images with full field of view (FOV) from undersampled Fourier encoded data. Sensitivity information must either be given a priori, or estimated from calibration data acquired along with the actual image data. In this study, two approaches are presented, which require very little or no additional data at all for calibration in two-dimensional multislice acquisitions. Instead of additional data, information from spatially adjacent slices is used to estimate coil sensitivity information, thereby increasing the efficiency of parallel imaging. The proposed approaches rely on the assumption that over a small range of slices, coil sensitivities vary smoothly in slice direction. Both methods are implemented as variants of the GRAPPA algorithm. For a given effective acceleration, superior image quality is achieved compared to the conventional GRAPPA method. For the latter calibration lines for coil weight computation must be acquired in addition to the undersampled k-spaces for coil weight computation, thus requiring higher k-space undersampling, that is, a higher reduction factor to achieve the same effective acceleration. The proposed methods are particularly suitable to speed up parallel imaging for clinical applications where the reduction factor is limited to two or three.

Highly k-t-space-accelerated phase-contrast MRI.

The purpose of this study was to combine a recently introduced spatiotemporal parallel imaging technique, PEAK-GRAPPA (parallel MRI with extended and averaged generalized autocalibrating partially parallel acquisition), with two-dimensional (2D) cine phase-contrast velocity mapping. Phase-contrast MRI was applied to measure the blood flow in the thoracic aorta and the myocardial motion of the left ventricle. To evaluate the performance of different reconstruction methods, fully acquired k-space data sets were used to compare conventional parallel imaging using GRAPPA with reduction factors of R = 2-6 and PEAK-GRAPPA as well as sliding window reconstruction with reduction factors R = 2-12 (net acceleration factors up to 5.2). PEAK-GRAPPA reconstruction resulted in improved image quality with considerably reduced artifacts, which was also supported by error analysis. To analyze potential blurring or low-pass filtering effects of spatiotemporal PEAK-GRAPPA, the velocity time courses of aortic flow and myocardial tissue motion were evaluated and compared with conventional image reconstructions. Quantitative comparisons of blood flow velocities and pixel-wise correlation analysis of velocities highlight the potential of PEAK-GRAPPA for highly accelerated dynamic phase-contrast velocity mapping.

Detection of pulmonary nodules with move-during-scan magnetic resonance imaging using a free-breathing turbo inversion recovery magnitude sequence.

PURPOSE: Detection of pulmonary metastases is still a challenging task for magnetic resonance imaging (MRI). It was the aim of this study to evaluate the potential of a free-breathing move-during-scan turbo inversion recovery magnitude sequence for the detection of pulmonary nodules. MATERIALS AND METHODS: The sensitivities and positive-predictive values of 2 radiologists to detect pulmonary nodules in 41 move-during-scan MRI examinations of 38 patients with different malignancies were calculated and subgroup analyses according to lesion size and localization were performed. Multidetector computed tomography served as the standard of reference. Additionally, 6 radiologists rated the confidence for the presence of nodular lesions in 212 regions-of-interest, which were randomly selected to represent lesions of various sizes as well as negative findings. Receiver-operator-characteristic was performed. RESULTS: Three hundred twenty-one nodules were found in 30 patients by multidetector computed tomography. Sensitivity and specificity of MRI to detect pulmonary nodules larger than 3 mm on a per-patient basis were 81.8% and 94.7%, respectively. On a per-lesion basis, MRI revealed a sensitivity of 79.0% to 80.7% for lesions larger than 3 mm, if high conspicuity ratings were counted as positive, and 84.6%, if medium and high conspicuity ratings were counted as positive. Sensitivity increased uniformly with lesion size, and all lesions larger than 12 mm were detected. Receiver-operator-characteristic analysis revealed a mean accuracy of 0.90 and sensitivities over 90% for lesions larger than 3 mm with a specificity of 96.1%. For lesions larger than 6 mm the accuracy was 0.99. CONCLUSION: Detection of pulmonary nodules with a move-during-scan turbo inversion recovery magnitude sequence is feasible. Excellent detection of lesions larger than 6 mm is achievable with free-breathing moving-table MRI.

Accurate characterization of SiO2 thin films using surface acoustic waves.

We have investigated the acoustic properties of silicon dioxide thin films. Therefore, we determined the phase velocity dispersion of LiNbO3 substrate covered with SiO2 deposited by a plasma enhanced chemical vapor deposition and a physical vapor deposition (PVD) process using differential delay lines and laser ultrasonic method. The density p and the elastic constants (c11 and c44) can be extracted by fitting corresponding finite element simulations to the phase velocities within an accuracy of at least +4%. Additionally, we propose two methods to improve the accuracy of the phase velocity determination by dealing with film thickness variation of the PVD process.

Accurate semiautomatic assessment of ligament length variations from MRI data.

PURPOSE: A semiautomatic method for the assessment of ligament length variations during different joint positions based on MRI data is proposed. METHODS: Ligament lengths are represented as distances between points marking characteristic locations in the ligament insertion regions on the bones. These points are defined manually for one single reference joint position and for all other joint positions they are automatically mapped with high accuracy to the correct locations using image registration methods. The methodology is validated using data from 16 volunteers depicting the coracoclavicular ligaments in the left shoulder during different arm abductions. RESULTS: The method yielded a superior reproducibility of the point locations over different joint positions compared to manual point marking. Significant ligament length variations were found for different abductions which was not possible with manual measurements. Acquisition related geometric distortions and inaccuracies during the registration and segmentation process were small. CONCLUSIONS: The proposed method provides superior accuracy for the in vivo analysis of ligament dynamics compared to manual measurements. This permits a better understanding of the ligament behavior during joint motion and offers new possibilities for presurgical planning which to date has not been possible with manual data analysis.

In vivo analysis of coracoclavicular ligament kinematics during shoulder abduction.

BACKGROUND: Anatomic reconstruction of the coracoclavicular ligaments for the treatment of acromioclavicular joint separations provides superior biomechanical stability compared with other procedures. Clavicular and coracoidal footprints of the conoid ligament (CL) and the trapezoid ligament (TL) are well described. So far, little is known about their kinematics and the changes of the coracoclavicular distance during shoulder abduction. HYPOTHESIS: The coracoclavicular distance along the coracoclavicular ligaments changes significantly with shoulder abduction and weightbearing. STUDY DESIGN: Descriptive laboratory study. METHODS: With use of an open magnetic resonance imaging scanner, the shoulders of 13 healthy volunteers were examined in supine and sitting positions. Three-dimensional magnetic resonance images of the shoulders were obtained in 30° increments of abduction (0°-120°). A manual segmentation of the scapula, the clavicle, and the coracoclavicular ligaments was performed. The insertion points of the coracoclavicular ligaments were identified, and automated measures along the ligamentous course were carried out. RESULTS: During transfer from the lying to sitting position, the coracoclavicular distance showed significant lengthening of 3 mm along the center of the CL, which significantly increased another 3 mm during shoulder abduction to a total lengthening of 6 mm. In the supine position, the coracoclavicular distance along the TL did not elongate significantly. In the sitting position, the distance along the medial portion of the TL shortened significantly, whereas the distance along the center portion did not elongate significantly during shoulder abduction. CONCLUSION: The distances between the coracoclavicular insertion points depend on both patient and shoulder positioning. To prevent overconstraining of the graft, the CL should be fixated during 90° to 120° of shoulder abduction in a sitting position. Isometric reconstruction of the TL can be achieved if precise fixation of the graft at the centers of the conoidal and clavicular footprints is performed.

k-t-Space accelerated myocardial perfusion.

PURPOSE: To investigate the performance of the recently introduced spatiotemporal parallel imaging technique called parallel MRI with extended and averaged generalized autocalibrating partially parallel acquisitions (GRAPPA) kernels (PEAK-GRAPPA) for myocardial perfusion measurements. MATERIALS AND METHODS: A study with 11 patients with myocardial infarction was performed to compare nonaccelerated perfusion imaging, i.e., fully acquired k-space data, with the results of conventional GRAPPA and PEAK-GRAPPA with a net acceleration factor of 2.4 to 3.4. Signal time courses reflecting the passage of the contrast agent bolus in different regions of the heart were evaluated for these different reconstruction methods. RESULTS: Reconstruction with PEAK-GRAPPA demonstrated considerably improved image quality compared to conventional GRAPPA. In addition, signal time courses for PEAK-GRAPPA demonstrated an excellent agreement compared to full k-space data, which is necessary for an accurate qualitative and quantitative assessment of myocardial perfusion. CONCLUSION: Qualitative and quantitative results of patient measurements illustrate that the temporal fidelity of nonperiodic processes such as myocardial perfusion are preserved with PEAK-GRAPPA up to net acceleration factors of more than 3 while showing a superior image quality compared to conventional GRAPPA and a sliding-window reconstruction.

Parallel MRI with extended and averaged GRAPPA kernels (PEAK-GRAPPA): optimized spatiotemporal dynamic imaging.

PURPOSE: To evaluate an optimized k-t-space related reconstruction method for dynamic magnetic resonance imaging (MRI), a method called PEAK-GRAPPA (Parallel MRI with Extended and Averaged GRAPPA Kernels) is presented which is based on an extended spatiotemporal GRAPPA kernel in combination with temporal averaging of coil weights. MATERIALS AND METHODS: The PEAK-GRAPPA kernel consists of a uniform geometry with several spatial and temporal source points from acquired k-space lines and several target points from missing k-space lines. In order to improve the quality of coil weight estimation sets of coil weights are averaged over the temporal dimension. RESULTS: The kernel geometry leads to strongly decreased reconstruction times compared to the recently introduced k-t-GRAPPA using different kernel geometries with only one target point per kernel to fit. Improved results were obtained in terms of the root mean square error and the signal-to-noise ratio as demonstrated by in vivo cardiac imaging. CONCLUSION: Using a uniform kernel geometry for weight estimation with the properties of uncorrelated noise of different acquired timeframes, optimized results were achieved in terms of error level, signal-to-noise ratio, and reconstruction time.