Purpose
The introduction of a linear safety factor to address peak local specific absorption rate (pSAR10g) uncertainties (eg, intersubject variation, modeling inaccuracies) bears one considerable ...drawback: It often results in over‐conservative scanning constraints. We present a more efficient approach to define a variable safety margin based on the conditional probability density function of the effectively obtained pSAR10g value, given the estimated pSAR10g value.
Methods
The conditional probability density function can be estimated from previously simulated data. A representative set of true and estimated pSAR10g samples was generated by means of our database of 23 subject‐specific models with an 8‐fractionated dipole array for prostate imaging at 7 T. The conditional probability density function was calculated for each possible estimated pSAR10g value and used to determine the corresponding safety margin with an arbitrary low probability of underestimation. This approach was applied to five state‐of‐the‐art local SAR estimation methods, namely: (1) using just the generic body model “Duke”; (2) using our model library to assess the maximum pSAR10g value over all models; (3) using the most representative “local SAR model”; (4) using the five most representative local SAR models; and (5) using a recently developed deep learning–based method.
Results
Compared with the more conventional safety factor, the conditional safety‐margin approach results in lower (up to 30%) mean overestimation for all investigated local SAR estimation methods.
Conclusion
The proposed probabilistic approach for pSAR10g correction allows more accurate local SAR assessment with much lower overestimation, while a predefined level of underestimation is accepted (eg, 0.1%).
Measuring in vivo dynamics can yield valuable information for studying the functioning of the cardiovascular or the musculoskeletal system and for the diagnosis of related diseases. MRI is a powerful ...medical imaging modality, but it shows severe limitations when dealing with motion at high spatial and temporal resolutions. In this work, a method called spectro-dynamic MRI is proposed, which can identify dynamical information directly from k-space data. It combines a measurement model, relating the measured data in k-space to the displacement fields, and a dynamical model, introducing prior knowledge about the dynamics of a system. The data sampling process is tailored to compute spatial and temporal derivatives in the spectral domain at a high temporal resolution. Preliminary results from four simple pendulum setups for which the dynamics are explicitly known show that spectro-dynamic MRI can estimate motion fields from heavily undersampled data on a millisecond timescale. Furthermore, the length of the pendula and the stiffness of the spring can be identified as the dynamical system's parameters, giving additional information about the systems under investigation.
In brain/head-and-neck radiotherapy (RT), thermoplastic immobilization masks guarantee reproducible patient positioning in treatment position between MRI, CT, and irradiation. Since immobilization ...masks do not fit in the diagnostic MR head/head-and-neck coils, flexible surface coils are used for MRI imaging in clinical practice. These coils are placed around the head/neck, in contact with the immobilization masks. However, the positioning of these flexible coils is technician dependent, thus leading to poor image reproducibility. Additionally, flexible surface coils have an inferior signal-to-noise-ratio (SNR) compared to diagnostic coils. The aim of this work was to create a new immobilization setup which fits into the diagnostic MR coils in order to enhance MR image quality and reproducibility. For this purpose, a practical immobilization setup was constructed. The performances of the standard clinical and the proposed setups were compared with four tests: SNR, image quality, motion restriction, and reproducibility of inter-fraction subject positioning. The new immobilization setup resulted in 3.4 times higher SNR values on average than the standard setup, except directly below the flexible surface coils where similar SNR was observed. Overall, the image quality was superior for brain/head-and-neck images acquired with the proposed RT setup. Comparable motion restriction in feet-head/left-right directions (maximum motion ≈1 mm) and comparable inter-fraction repositioning accuracy (mean inter-fraction movement 1 ± 0.5 mm) were observed for the standard and the new setup.
Magnetic resonance fingerprinting (MRF) is able to estimate multiple quantitative tissue parameters from a relatively short acquisition. The main characteristic of an MRF sequence is the simultaneous ...application of 1) transient states excitation and 2) highly undersampled k-space. Despite the promising empirical results obtained with MRF, no work has appeared that formally describes the combined impact of these two aspectson the reconstruction accuracy. In this paper, a mathematical model is derived that directly relates the time-varying RF excitation and the k-space sampling to the spatially dependent reconstruction errors. A subsequent in-depth analysis identifies the mechanisms by which MRF sequence properties affect accuracy, providing a formal explanation of several empirically observed or intuitively understood facts. The new insights are obtained which show how this analytical framework could be used to improve the MRF protocol.
Magnetic resonance fingerprinting (MRF) is able to estimate multiple quantitative tissue parameters from a relatively short acquisition. The main characteristic of an MRF sequence is the simultaneous ...application of 1) transient states excitation and 2) highly undersampled k -space. Despite the promising empirical results obtained with MRF, no work has appeared that formally describes the combined impact of these two aspects on the reconstruction accuracy. In this paper, a mathematical model is derived that directly relates the time-varying RF excitation and the k -space sampling to the spatially dependent reconstruction errors. A subsequent in-depth analysis identifies the mechanisms by which MRF sequence properties affect accuracy, providing a formal explanation of several empirically observed or intuitively understood facts. The new insights are obtained which show how this analytical framework could be used to improve the MRF protocol.
Purpose
Image quality obtained for brain imaging at 7T can be hampered by inhomogeneities in the static magnetic field, B0, and the RF electromagnetic field, B1. In imaging sequences such as ...fluid‐attenuated inversion recovery (FLAIR), which is used to assess neurological disorders, these inhomogeneities cause spatial variations in signal that can reduce clinical efficacy. In this work, we aim to correct for signal inhomogeneities to ensure whole‐brain coverage with 3D FLAIR at 7T.
Methods
The direct signal control (DSC) framework was used to optimize channel weightings applied to the 8 transmit channels used in this work on a pulse‐by‐pulse basis through the echo train in the FLAIR sequences. 3D FLAIR brain images were acquired on 5 different subjects and compared with imaging using a quadrature‐like mode of the transmit array. Precomputed “universal” DSC solutions calculated from a separate set of 5 subjects were also explored.
Results
DSC consistently enabled improved imaging across all subjects, with no dropouts in signal seen over the entire brain volume, which contrasted with imaging in quadrature mode. Further, the universal DSC solutions also consistently improved imaging despite not being optimized specifically for the subject being imaged.
Conclusion
3D FLAIR brain imaging at 7T is substantially improved using DSC and is able to recover regions of low signal without increasing imaging time or interecho spacing.
MR-STAT is an emerging quantitative magnetic resonance imaging technique which aims at obtaining multi-parametric tissue parameter maps from single short scans. It describes the relationship between ...the spatial-domain tissue parameters and the time-domain measured signal by using a comprehensive, volumetric forward model. The MR-STAT reconstruction solves a large-scale nonlinear problem, thus is very computationally challenging. In previous work, MR-STAT reconstruction using Cartesian readout data was accelerated by approximating the Hessian matrix with sparse, banded blocks, and can be done on high performance CPU clusters with tens of minutes. In the current work, we propose an accelerated Cartesian MR-STAT algorithm incorporating two different strategies: firstly, a neural network is trained as a fast surrogate to learn the magnetization signal not only in the full time-domain but also in the compressed low-rank domain; secondly, based on the surrogate model, the Cartesian MR-STAT problem is re-formulated and split into smaller sub-problems by the alternating direction method of multipliers. The proposed method substantially reduces the computational requirements for runtime and memory. Simulated and in-vivo balanced MR-STAT experiments show similar reconstruction results using the proposed algorithm compared to the previous sparse Hessian method, and the reconstruction times are at least 40 times shorter. Incorporating sensitivity encoding and regularization terms is straightforward, and allows for better image quality with a negligible increase in reconstruction time. The proposed algorithm could reconstruct both balanced and gradient-spoiled in-vivo data within 3 minutes on a desktop PC, and could thereby facilitate the translation of MR-STAT in clinical settings.
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 <inline-formula> <tex-math notation="LaTeX">{k} </tex-math></inline-formula>-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 (<inline-formula> <tex-math notation="LaTeX">{R}={0.975} \pm {0.0110} </tex-math></inline-formula>), and can be reconstructed with a total latency of 170 ms that is sufficient for real-time MR-guided abdominal radiotherapy.
Purpose
A purely experimental method for MRI‐based transfer function (TF) determination is presented. A TF characterizes the potential for radiofrequency heating of a linear implant by relating the ...incident tangential electric field to a scattered electric field at its tip. We utilize the previously introduced transfer matrix (TM) to determine transfer functions solely from the MR measurable quantities, that is, the
and transceive phase distributions. This technique can extend the current practice of phantom‐based TF assessment with dedicated experimental setup toward MR‐based methods that have the potential to assess the TF in more realistic situations.
Theory and Methods
An analytical description of the
magnitude and transceive phase distribution around a wire‐like implant was derived based on the TM. In this model, the background field is described using a superposition of spherical and cylindrical harmonics while the transfer matrix is parameterized using a previously introduced attenuated wave model. This analytical description can be used to estimate the transfer matrix and transfer function based on the measured
distribution.
Results
The TF was successfully determined for 2 mock‐up implants: a 20‐cm bare copper wire and a 20‐cm insulated copper wire with 10 mm of insulation stripped at both endings in respectively 4 and 3 different trajectories. The measured TFs show a strong correlation with a reference determined from simulations and between the separate experiments with correlation coefficients above 0.96 between all TFs. Compared to the simulated TF, the maximum deviation in the estimated tip field is 9.4% and 12.2% for the bare and insulated wire, respectively.
Conclusions
A method has been developed to measure the TF of medical implants using MRI experiments. Jointly fitting the incident and scattered
distributions with an analytical description based on the transfer matrix enables accurate determination of the TF of 2 test implants. The presented method no longer needs input from simulated data and can therefore, in principle, be used to measure TF's in test animals or corpses.