Thermal noise variance of a receive radiofrequency coil as a respiratory motion sensor

Beat Pilot Tone (BPT): Simultaneous MRI and RF motion sensing at arbitrary frequencies

PURPOSE

To introduce a simple system exploitation with the potential to turn MRI scanners into general-purpose radiofrequency (RF) motion monitoring systems.

METHODS

Inspired by Pilot Tone (PT), this work proposes Beat Pilot Tone (BPT), in which two or more RF tones at arbitrary frequencies are transmitted continuously during the scan. These tones create motion-modulated standing wave patterns that are sensed by the receiver coil array, incidentally mixed by intermodulation in the receiver chain, and digitized simultaneously with the MRI data. BPT can operate at almost any frequency as long as the intermodulation products lie within the bandwidth of the receivers. BPT's mechanism is explained in electromagnetic simulations and validated experimentally.

RESULTS

Phantom and volunteer experiments over a range of transmit frequencies suggest that BPT may offer frequency-dependent sensitivity to motion. Using a semi-flexible anterior receiver array, BPT appears to sense cardiac-induced body vibrations at microwave frequencies (≥ $$ \ge $$ 1.2 GHz). At lower frequencies, it exhibits a similar cardiac signal shape to PT, likely due to blood volume changes. Other volunteer experiments with respiratory, bulk, and head motion show that BPT can achieve greater sensitivity to motion than PT and greater separability between motion types. Basic multiple-input multiple-output (4 × 22 $$ 4\times 22 $$ MIMO) operation with simultaneous PT and BPT in head motion is demonstrated using two transmit antennas and a 22-channel head-neck coil.

CONCLUSION

BPT may offer a rich source of motion information that is frequency-dependent, simultaneous, and complementary to PT and the MRI exam.

Current Research Status of Respiratory Motion for Thorax and Abdominal Treatment: A Systematic Review

Malignant tumors have become one of the serious public health problems in human safety and health, among which the chest and abdomen diseases account for the largest proportion. Early diagnosis and treatment can effectively improve the survival rate of patients. However, respiratory motion in the chest and abdomen can lead to uncertainty in the shape, volume, and location of the tumor, making treatment of the chest and abdomen difficult. Therefore, compensation for respiratory motion is very important in clinical treatment. The purpose of this review was to discuss the research and development of respiratory movement monitoring and prediction in thoracic and abdominal surgery, as well as introduce the current research status. The integration of modern respiratory motion compensation technology with advanced sensor detection technology, medical-image-guided therapy, and artificial intelligence technology is discussed and analyzed. The future research direction of intraoperative thoracic and abdominal respiratory motion compensation should be non-invasive, non-contact, use a low dose, and involve intelligent development. The complexity of the surgical environment, the constraints on the accuracy of existing image guidance devices, and the latency of data transmission are all present technical challenges.

Ultra-fast multi-parametric 4D-MRI image reconstruction for real-time applications using a downsampling-invariant deformable registration (D2R) model

BACKGROUND AND PURPOSE

Motion estimation from severely downsampled 4D-MRI is essential for real-time imaging and tumor tracking. This simulation study developed a novel deep learning model for simultaneous MR image reconstruction and motion estimation, named the Downsampling-Invariant Deformable Registration (D2R) model.

MATERIALS AND METHODS

Forty-three patients undergoing radiotherapy for liver tumors were recruited for model training and internal validation. Five prospective patients from another center were recruited for external validation. Patients received 4D-MRI scans and 3D MRI scans. The 4D-MRI was retrospectively down-sampled to simulate real-time acquisition. Motion estimation was performed using the proposed D2R model. The accuracy and robustness of the proposed D2R model and baseline methods, including Demons, Elastix, the parametric total variation (pTV) algorithm, and VoxelMorph, were compared. High-quality (HQ) 4D-MR images were also constructed using the D2R model for real-time imaging feasibility verification. The image quality and motion accuracy of the constructed HQ 4D-MRI were evaluated.

RESULTS

The D2R model showed significantly superior and robust registration performance than all the baseline methods at downsampling factors up to 500. HQ T1-weighted and T2-weighted 4D-MR images were also successfully constructed with significantly improved image quality, sub-voxel level motion error, and real-time efficiency. External validation demonstrated the robustness and generalizability of the technique.

CONCLUSION

In this study, we developed a novel D2R model for deformation estimation of downsampled 4D-MR images. HQ 4D-MR images were successfully constructed using the D2R model. This model may expand the clinical implementation of 4D-MRI for real-time motion management during liver cancer treatment.

Gaussian Processes for real-time 3D motion and uncertainty estimation during MR-guided radiotherapy

Respiratory motion during radiotherapy causes uncertainty in the tumor's location, which is typically addressed by an increased radiation area and a decreased dose. As a result, the treatments' efficacy is reduced. The recently proposed hybrid MR-linac scanner holds the promise to efficiently deal with such respiratory motion through real-time adaptive MR-guided radiotherapy (MRgRT). For MRgRT, motion-fields should be estimated from MR-data and the radiotherapy plan should be adapted in real-time according to the estimated motion-fields. All of this should be performed with a total latency of maximally 200 ms, including data acquisition and reconstruction. A measure of confidence in such estimated motion-fields is highly desirable, for instance to ensure the patient's safety in case of unexpected and undesirable motion. In this work, we propose a framework based on Gaussian Processes to infer 3D motion-fields and uncertainty maps in real-time from only three readouts of MR-data. We demonstrated an inference frame rate up to 69 Hz including data acquisition and reconstruction, thereby exploiting the limited amount of required MR-data. Additionally, we designed a rejection criterion based on the motion-field uncertainty maps to demonstrate the framework's potential for quality assurance. The framework was validated in silico and in vivo on healthy volunteer data (n=5) acquired using an MR-linac, thereby taking into account different breathing patterns and controlled bulk motion. Results indicate end-point-errors with a 75th percentile below 1 mm in silico, and a correct detection of erroneous motion estimates with the rejection criterion. Altogether, the results show the potential of the framework for application in real-time MR-guided radiotherapy with an MR-linac.

Pilot tone-based prospective correction of respiratory motion for free-breathing myocardial T1 mapping

OBJECTIVE

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 METHODS

During 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.

RESULTS

Non-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.

DISCUSSION

Our proposed PT-based motion correction approach allows for T1 mapping during free-breathing with the same accuracy as a corresponding breathhold T1 mapping scan.

Real-Time Non-Rigid 3D Respiratory Motion Estimation for MR-Guided Radiotherapy Using MR-MOTUS

The MR-Linac is a combination of an MR-scanner and radiotherapy linear accelerator (Linac) which holds the promise to increase the precision of radiotherapy treatments with MR-guided radiotherapy by monitoring motion during radiotherapy with MRI, and adjusting the radiotherapy plan accordingly. Optimal MR-guidance for respiratory motion during radiotherapy requires MR-based 3D motion estimation with a latency of 200-500 ms. Currently this is still challenging since typical methods rely on MR-images, and are therefore limited by the 3D MR-imaging latency. In this work, we present a method to perform non-rigid 3D respiratory motion estimation with 170 ms latency, including both acquisition and reconstruction. The proposed method called real-time low-rank MR-MOTUS reconstructs motion-fields directly from k -space data, and leverages an explicit low-rank decomposition of motion-fields to split the large scale 3D+t motion-field reconstruction problem posed in our previous work into two parts: (I) a medium-scale offline preparation phase and (II) a small-scale online inference phase which exploits the results of the offline phase for real-time computations. The method was validated on free-breathing data of five volunteers, acquired with a 1.5T Elekta Unity MR-Linac. Results show that the reconstructed 3D motion-field are anatomically plausible, highly correlated with a self-navigation motion surrogate ( R=0.975 ±0.0110 ), and can be reconstructed with a total latency of 170 ms that is sufficient for real-time MR-guided abdominal radiotherapy.

Pilot tone-based motion correction for prospective respiratory compensated cardiac cine MRI

PURPOSE

To evaluate prospective motion correction using the pilot tone (PT) as a quantitative respiratory motion signal with high temporal resolution for cardiac cine images during free breathing.

METHODS

Before cine data acquisition, a short prescan was performed, calibrating the PT to the respiratory-induced heart motion using respiratory-resolved real-time images. The calibrated PT was then applied for nearly real-time prospective motion correction of cine MRI through slice tracking (ie, updating the slice position before every readout). Additionally, in-plane motion correction was performed retrospectively also based on the calibrated PT data. The proposed method was evaluated in a moving phantom and 10 healthy volunteers.

RESULTS

The PT showed very good correlation to the phantom motion. In volunteer studies using a long-term scan over 7.96 ± 1.40 min, the mean absolute error between registered and predicted motion from the PT was 1.44 ± 0.46 mm in head-feet and 0.46 ± 0.07 mm in anterior-posterior direction. Irregular breathing could also be corrected well with the PT. The PT motion correction leads to a significant improvement of contrast-to-noise ratio by 68% (P ≤ .01) between blood pool and myocardium and sharpness of endocardium by 24% (P = .04) in comparison to uncorrected data. The image score, which refers to the cine image quality, has improved with the utilization of the proposed PT motion correction.

CONCLUSION

The proposed approach provides respiratory motion-corrected cine images of the heart with improved image quality and a high scan efficiency using the PT. The PT is independent of the MR acquisition, making this a very flexible motion-correction approach.

The noise navigator for MRI-guided radiotherapy: an independent method to detect physiological motion

Motion is problematic during radiotherapy as it could lead to potential underdosage of the tumor, and/or overdosage in organs-at-risk. A solution is adaptive radiotherapy guided by magnetic resonance imaging (MRI). MRI allows for imaging of target volumes and organs-at-risk before and during treatment delivery with superb soft tissue contrast in any desired orientation, enabling motion management by means of (real-time) adaptive radiotherapy. The noise navigator, which is independent of the MR signal, could serve as a secondary motion detection method in synergy with MR imaging. The feasibility of respiratory motion detection by means of the noise navigator was demonstrated previously. Furthermore, from electromagnetic simulations we know that the noise navigator is sensitive to tissue displacement and thus could in principle be used for the detection of various types of motion. In this study we demonstrate the detection of various types of motion for three anatomical use cases of MRI-guided radiotherapy, i.e. torso (bulk movement and variable breathing), head-and-neck (swallowing) and cardiac. Furthermore, it is shown that the noise navigator can detect bulk movement, variable breathing and swallowing on a hybrid 1.5 T MRI-linac system. Cardiac activity detection through the noise navigator seems feasible in an MRI-guided radiotherapy setting, but needs further optimization. The noise navigator is a versatile and fast (millisecond temporal resolution) motion detection method independent of MR signal that could serve as an independent verification method to detect the occurrence of motion in synergy with real-time MRI-guided radiotherapy.

The noise navigator: a surrogate for respiratory-correlated 4D-MRI for motion characterization in radiotherapy

Respiratory-correlated 4D-MRI can characterize respiratory-induced motion of tumors and organs-at-risk for radiotherapy treatment planning and is a necessity for image guidance of moving tumors treated on an MRI-linac. Essential for 4D-MRI generation is a robust respiratory surrogate signal. We investigated the feasibility of the noise navigator as respiratory surrogate signal for 4D-MRI generation. The noise navigator is based on the respiratory-induced modulation of the thermal noise variance measured by the receive coils during MR acquisition and thus is inherently present and synchronized with MRI data acquisition. Additionally, the noise navigator can be combined with any rectilinear readout strategy (e.g. radial and cartesian) and is independent of MR image contrast and imaging orientation. For radiotherapy applications, the noise navigator provides a robust respiratory signal for patients scanned with an elevated coil setup. This is particularly attractive for widely used cartesian sequences where currently a non-interfering self-navigation means is lacking for MRI-based simulation and MRI-guided radiotherapy. The feasibility of 4D-MRI generation with the noise navigator as respiratory surrogate signal was demonstrated for both cartesian and radial readout strategies in radiotherapy setup on four healthy volunteers and two radiotherapy patients on a dedicated 1.5 T MRI scanner and two radiotherapy patients on a 1.5 T MRI-linac system. Moreover, the respiratory-correlated 4D-MR images showed liver motion comparable to a reference 2D cine MRI series for the volunteers. For 2D cartesian cine MRI acquisitions, both the noise navigator and respiratory bellows were benchmarked against an image navigator. Respiratory phase detection based on the noise navigator agreed 1.4 times better with the image navigator than the respiratory bellows did. For a 3D Stack-of-Stars acquisitions, the noise navigator was compared to radial self-navigation and a 1.7 times higher respiratory phase detection agreement was observed than for the respiratory bellows compared to radial self-navigation.

MR-MOTUS: model-based non-rigid motion estimation for MR-guided radiotherapy using a reference image and minimal k-space data

Time-resolved motion estimation from MRI data has received an increasing amount of interest due to the advent of the MR-Linac. The combination of an MRI scanner and a linear accelerator enables radiation plan adaptation based on internal organ motion estimated from MRI data. However, time-resolved estimation of this motion from MRI data still remains a challenge. In light of this application, we propose MR-MOTUS, a framework to estimate non-rigid 3D motion from minimal k-space data. MR-MOTUS consists of two main components: (1) a signal model that explicitly relates the k-space signal of a deforming object to non-rigid motion-fields and a reference image, and (2) model-based reconstructions of the non-rigid motion-fields directly from k-space data. Using an a priori available reference image and the fact that internal body motion exhibits a high level of spatial correlation, we represent the motion-fields in a low-dimensional space and reconstruct them from minimal k-space data that can be acquired very rapidly. The signal model is validated through numerical experiments with a digital 3D phantom and motion-fields are reconstructed from retrospectively undersampled in vivo head and abdomen data using various undersampling strategies. A comparison is made with state-of-the-art image registration performed on images reconstructed from the same undersampled data. Results show that MR-MOTUS reconstructs in vivo 3D rigid head motion from 474-fold retrospectively downsampled k-space data, and in vivo non-rigid 3D respiratory motion from 63-fold retrospectively undersampled k-space data. Preliminary results on prospectively undersampled data acquired with a 2D golden angle acquisition during free-breathing demonstrate the practical feasibility of the method.

Scattering matrix imaging pulse design for real-time respiration and cardiac motion monitoring

Purpose

The scattering matrix (S‐matrix) of a parallel transmit (pTx) coil is sensitive to physiological motion but requires additional monitoring RF pulses to be measured. In this work, we present and evaluate pTx RF pulse designs that simultaneously excite for imaging and measure the S‐matrix to generate real‐time motion signals without prolonging the image sequence.

Theory and Methods

Three pTx waveforms for measuring the S‐matrix were identified and superimposed onto the imaging excitation RF pulses: (1) time division multiplexing, (2) frequency division multiplexing, and (3) code division multiplexing. These 3 methods were evaluated in healthy volunteers for scattering sensitivity and image artefacts. The S‐matrix and real‐time motion signals were calculated on the image calculation environment of the MR scanner. Prospective cardiac triggers were identified in early systole as a high rate of change of the cardiac motion signal. Monitoring accuracy was compared against electrocardiogram or the imaged diaphragm position.

Results

All 3 monitoring approaches measure the S‐matrix during image excitation with quality correlated to input power. No image artefacts were observed for frequency multiplexing, and low energy artefacts were observed in the other methods. The accuracy of the achieved prospective cardiac gating was 15 ± 16 ms for breath hold and 24 ± 17 ms during free breathing. The diaphragm position prediction accuracy was 1.3 ± 0.9 mm. In all volunteers, good quality cine images were acquired for breath hold scans and dual gated CINEs were demonstrated.

Conclusion

The S‐matrix can be measured during image excitation to generate real‐time cardiac and respiratory motion signals for prospective gating. No artefacts are introduced when frequency division multiplexing is used.

Understanding the physical relations governing the noise navigator

Purpose

The noise navigator is a passive way to detect physiological motion occurring in a patient through thermal noise modulations measured by standard clinical radiofrequency receive coils. The aim is to gain a deeper understanding of the potential and applications of physiologically induced thermal noise modulations.

Methods

Numerical electromagnetic simulations and MR measurements were performed to investigate the relative contribution of tissue displacement versus modulation of the dielectric lung properties over the respiratory cycle, the impact of coil diameter and position with respect to the body. Furthermore, the spatial motion sensitivity of specific noise covariance matrix elements of a receive array was investigated.

Results

The influence of dielectric lung property variations on the noise variance is negligible compared to tissue displacement. Coil size affected the thermal noise variance modulation, but the location of the coil with respect to the body had a larger impact. The modulation depth of a 15 cm diameter stationary coil approximately 3 cm away from the chest (i.e. radiotherapy setup) was 39.7% compared to 4.2% for a coil of the same size on the chest, moving along with respiratory motion. A combination of particular noise covariance matrix elements creates a specific spatial sensitivity for motion.

Conclusions

The insight gained on the physical relations governing the noise navigator will allow for optimized use and development of new applications. An optimized combination of elements from the noise covariance matrix offer new ways of performing, e.g. motion tracking.

Respiration resolved imaging with continuous stable state 2D acquisition using linear frequency SWEEP

Purpose

To investigate the potential of continuous radiofrequency (RF) shifting (SWEEP) as a technique for creating densely sampled data while maintaining a stable signal state for dynamic imaging.

Methods

We present a method where a continuous stable state of magnetization is swept smoothly across the anatomy of interest, creating an efficient approach to dense multiple 2D slice imaging. This is achieved by introducing a linear frequency offset to successive RF pulses shifting the excited slice by a fraction of the slice thickness with each successive repeat times (TR). Simulations and in vivo imaging were performed to assess how this affects the measured signal. Free breathing, respiration resolved 4D volumes in fetal/placental imaging is explored as potential application of this method.

Results

The SWEEP method maintained a stable signal state over a full acquisition reducing artifacts from unstable magnetization. Simulations demonstrated that the effects of SWEEP on slice profiles was of the same order as that produced by physiological motion observed with conventional methods. Respiration resolved 4D data acquired with this method shows reduced respiration artifacts and resilience to non‐rigid and non‐cyclic motion.

Conclusions

The SWEEP method is presented as a technique for improved acquisition efficiency of densely sampled short‐TR 2D sequences. Using conventional slice excitation the number of RF pulses required to enter a true steady state is excessively high when using short‐TR 2D acquisitions, SWEEP circumvents this limitation by creating a stable signal state that is preserved between slices.

Cardiac gating using scattering of an 8-channel parallel transmit coil at 7T

PURPOSE

To establish a cardiac signal from scattering matrix or scattering coefficient measurements made on a 7T 8-channel parallel transmit (pTx) system, and to evaluate its use for cardiac gating.

METHODS

Measurements of the scattering matrix and scattering coefficients were acquired using a monitoring pulse sequence and during a standard cine acquisition, respectively. Postprocessing used an independent component analysis and gating feature identification. The effect of the phase of the excitation radiofrequency (RF) field ( B1+ shim) on the cardiac signal was simulated for multiple B1+ shim configurations, and cine images were reconstructed from both the scattering coefficients and electrocardiogram (ECG).

RESULTS

The cardiac motion signal was successfully identified in all subjects with a mean signal-to-noise ratio of 33.1 and 5.7 using the scattering matrix and scattering coefficient measurements, respectively. The dominant gating feature in the cardiac signal was a peak aligned with end-systole that occurred on average at 311 and 391 ms after the ECG trigger, with a mean standard deviation of 13.4 and 18.1 ms relative to ECG when using the scattering matrix and scattering coefficients measurements, respectively. The scattering coefficients showed a dependence on B1+ shim with some shim configurations not showing any cardiac signal. Cine images were successfully reconstructed using the scattering coefficients with minimal differences compared to those using ECG.

CONCLUSION

We have shown that the scattering of a pTx RF coil can be used to estimate a cardiac signal, and that scattering matrix and coefficients can be used to cardiac gate MRI acquisitions with the scattering matrix providing a superior cardiac signal. Magn Reson Med 80:633-640, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Prospective Respiration Detection in Magnetic Resonance Imaging by a Non-Interfering Noise Navigator

Passive monitoring of the thermal noise variances of the channels of a receive array was shown to reveal respiratory motion of the underlying anatomy, a so called "noise navigator". There is, however, an inevitable trade off between the accuracy and temporal resolution of the noise navigator due to its passive nature. A temporal filter has to be added to the noise navigator to accurately reveal respiration and retain temporal resolution. For real-time applications of the noise navigator, e.g., prospective motion correction or motion tracking, the added filter must be prospective. Thus a prospective Kalman filter was designed to predict respiration from the noise navigator without a temporal delay. The performance of the noise navigator enhanced by this prospective Kalman filter was explored and the robustness of the proposed method was assessed on healthy volunteers. The respiratory signal could be measured by the noise navigator independent of magnetic resonance acquisition. The calculated respiratory signal was qualitatively compared with the respiratory bellows. In addition, a strong linear relationship was found between the prospective noise navigator and a quantitative 2-D image navigator for measurements, including free and tasked breathing.

Adaptive bulk motion exclusion for improved robustness of abdominal magnetic resonance imaging

Non-Cartesian magnetic resonance imaging (MRI) sequences have shown great promise for abdominal examination during free breathing, but break down in the presence of bulk patient motion (i.e. voluntary or involuntary patient movement resulting in translation, rotation or elastic deformations of the body). This work describes a data-consistency-driven image stabilization technique that detects and excludes bulk movements during data acquisition. Bulk motion is identified from changes in the signal intensity distribution across different elements of a multi-channel receive coil array. A short free induction decay signal is acquired after excitation and used as a measure to determine alterations in the load distribution. The technique has been implemented on a clinical MR scanner and evaluated in the abdomen. Six volunteers were scanned and two radiologists scored the reconstructions. To show the applicability to other body areas, additional neck and knee images were acquired. Data corrupted by bulk motion were successfully detected and excluded from image reconstruction. An overall increase in image sharpness and reduction of streaking and shine-through artifacts were seen in the volunteer study, as well as in the neck and knee scans. The proposed technique enables automatic real-time detection and exclusion of bulk motion during MR examinations without user interaction. It may help to improve the reliability of pediatric MRI examinations without the use of sedation.

Diaphragm position can be accurately estimated from the scattering of a parallel transmit RF coil at 7 T

Purpose

To evaluate the use of radiofrequency scattering of a parallel transmit coil to track diaphragm motion.

Methods

Measurements made during radiofrequency excitation on an 8‐channel parallel transmit coil by the directional couplers of the radiofrequency safety monitor were combined and converted into diaphragm position. A 30‐s subject‐specific calibration with an MRI navigator was used to determine a diaphragm estimate from each directional‐coupler measure. Seven healthy volunteers were scanned at 7 T, in which images of the diaphragm were continuously acquired and directional couplers were monitored during excitation radiofrequency pulses. The ability to detect coughing was evaluated in one subject. The method was implemented on the scanner and evaluated for diaphragm gating of a free‐breathing cardiac cine.

Results

Six of the seven scans were successful. In these subjects, the root mean square difference between MRI and scattering estimation of the superior–inferior diaphragm position was 1.4 ± 0.5 mm. On the scanner, the position was calculated less than 2 ms after every radiofrequency pulse. A prospectively gated (echocardiogram and respiration) high‐resolution free‐breathing cine showed no respiratory artifact and sharp blood‐myocardium definition.

Conclusions

Transmit coil scattering is sensitive to diaphragm motion and provides rapid, quantitative, and accurate monitoring of respiration. Magn Reson Med 79:2164–2169, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Trimmed autocalibrating k-space estimation based on structured matrix completion

PURPOSE

Parallel imaging allows the reconstruction of undersampled data from multiple coils. This provides a means to reject and regenerate corrupt data (e.g. from motion artefact). The purpose of this work is to approach this problem using the SAKE parallel imaging method.

THEORY AND METHODS

Parallel imaging methods typically require calibration by fully sampling the center of k-space. This is a challenge in the presence of corrupted data, since the calibration data may be corrupted which leads to an errors-in-variables problem that cannot be solved by least squares or even iteratively reweighted least squares. The SAKE method, based on matrix completion and structured low rank approximation, was modified to detect and trim these errors from the data.

RESULTS

Simulated and actual corrupted datasets were reconstructed with SAKE, the proposed approach and a more standard reconstruction method (based on solving a linear equation) with a data rejection criterion. The proposed approach was found to reduce artefacts considerably in comparison to the other two methods.

CONCLUSION

SAKE with data trimming improves on previous methods for reconstructing images from grossly corrupted data.

Respiratory motion model based on the noise covariance matrix of a receive array

PURPOSE

Tracking of the internal anatomy by means of a motion model that uses the MR-derived motion fields and noise covariance matrix (NCM) dynamic as a surrogate signal.

METHODS

A 2D respiratory motion model was developed based on the MR-derived motion fields and the NCM of a receive array used in MRI. Temporal dynamics of the NCM were used as a motion surrogate for a linear correspondence motion model. The model performance was tested on five healthy volunteers with a liver as the target. The motion fields were calculated from the cineMR frames with an optical flow registration tool.

RESULTS

The model estimated the liver motion with an average residual error of 2.3 mm (13% of the motion amplitude). The model formation takes 3 min and the model latency was 0.5 s in the current implementation. The limiting factor for the latency is the current update time of the NCM (0.48 s), which in principle can be reduced to 0.004 s with an alternative way to determine the NCM.

CONCLUSIONS

The 2D respiratory motion of the liver can be effectively estimated with the linear motion model that uses the temporal behavior of the NCM as motion surrogate. Magn Reson Med 79:1730-1735, 2018. © 2017 International Society for Magnetic Resonance in Medicine.