MRI B0 and B1 Map Correction for Conductivity Accuracy
Find Innovative SolutionsGenerate Solutions
Solution Overview
Problem
Magnetic Resonance Imaging (MRI) systems face challenges in accurately accounting for geometric distortions caused by B0 field inhomogeneities, which affect voxel size and lead to errors in calculating electrical conductivity and permittivity using differential equations.
Innovation Solution
The implementation of a method that uses Dixon pulse sequence data to estimate B0 inhomogeneity and B1 phase maps, allowing for the calculation of zero echo time B1 phase maps and subsequent electrical conductivity maps, while also segmenting fat and water images to improve boundary conditions for more accurate conductivity calculations.
Engineering Contradictions & Design Principles
Engineering Contradiction Analysis
1Measurement precision
If standard MRI imaging methods are used, then imaging speed and productivity are maintained, but geometric distortions occur due to B0 field inhomogeneities leading to inaccurate voxel size and conductivity calculations
Solution Approach 1:
The patent segments the imaging process into distinct components: acquiring B0 inhomogeneity maps, acquiring B1 phase maps, calculating geometric distortion corrections, and applying these corrections to derive accurate electrical properties. This segmentation allows each component to be optimized independently while maintaining overall system productivity.
Solution Approach 2:
The patent performs preliminary calculations of B0 inhomogeneity maps and geometric distortion corrections before conducting the main conductivity mapping. By pre-characterizing the field inhomogeneities and computing correction factors in advance, the method eliminates the need for repeated corrective measurements during the main imaging sequence, thus maintaining productivity while improving precision.
2Measurement precision
If geometric distortion correction is applied using B0 maps, then measurement precision of voxel size is improved, but the complexity of data processing and calculation increases
Solution Approach 1:
The patent calculates the geometric distortion correction factors and adjusted voxel sizes in advance using pre-acquired B0 maps. These correction parameters are stored and reused during the main imaging sequence, eliminating the need for real-time calculation of distortion corrections and their derivatives, thus reducing processing time while maintaining accuracy.
Solution Approach 2:
The system uses its own acquired B0 inhomogeneity data to generate the correction factors it needs, making the process self-sufficient. The B0 maps acquired during the imaging session are directly used to compute the geometric distortion corrections for that same session, eliminating the need for external calibration data or complex reference measurements.
3Reliability
If n-point Dixon method is used to estimate B0 and B1 maps, then reliability of field mapping is improved, but the scanning time and energy consumption increase
Solution Approach 1:
The patent employs the n-point Dixon method to simultaneously acquire multiple pieces of information: B0 inhomogeneity maps, B1 phase maps, and water-fat separation data all in a single multi-echo sequence. This multi-functionality allows reliable field mapping without requiring separate dedicated sequences for each parameter, thus maintaining productivity while improving reliability.
Solution Approach 2:
The patent maintains continuous data acquisition across multiple echoes in the Dixon sequence, where each echo contributes additional information for robust B0 and B1 estimation. By continuously acquiring data rather than using discrete separate measurements, the method achieves high reliability through multiple samples while minimizing total scanning time through efficient use of each echo.
Applied Scientific Principles
This section explains which scientific principles are used to turn an abstract innovation direction into a practical engineering solution.
Function Achieved in This Case
This approach enables the correction of geometric distortions and enhances the accuracy of electrical conductivity and permittivity mapping, reducing errors and improving the quality of MRI data for applications like hyperthermia treatment planning and SAR modeling.
Implementation Method 1
Dixon pulse sequence data, wherein the Dixon pulse sequence data comprises instructions for controlling the magnetic resonance imaging system to acquire magnetic resonance data according to an n-point Dixon method
Implementation Method 2
estimating a B0 inhomogeneity map and estimating an estimated B1 phase map by analyzing the magnetic resonance data according to the n-point Dixon method
Implementation Method 3
calculating a zero echo time B1 phase map by interpolating the estimated B1 phase map to an echo time of zero using the B0 inhomogeneity map
Data Source
AI summary
The invention provides for a magnetic resonance imaging system (100, 300, 100) for acquiring magnetic resonance data (110, 1104) from a subject (118) within an imaging zone (108). The magnetic resonance imaging system comprises a memory (136) for storing machine executable instructions (160, 162, 164, 166, 316) and pulse sequence data (140, 1102). The pulse sequence data comprises instructions for controlling the magnetic resonance imaging system to acquire magnetic resonance data according to a magnetic resonance imaging method. The magnetic resonance imaging system further comprises a processor (130) for controlling the magnetic resonance imaging system. Execution of the machine executable instructions causes the processor to: acquire (1200) the magnetic resonance data by controlling the magnetic resonance imaging system with the pulse sequence data; calculate (1202) a B0 inhomogeneity map (148) by analyzing the magnetic resonance data according to the magnetic resonance imaging method, calculate (1204) a B1 phase map (150) and/or a B1 amplitude map (1106) by analyzing the magnetic resonance data according to the magnetic resonance imaging method; and calculate (1206) a second derivative (1110) of the B1 phase map and/or a second derivative of the B1 magnitude map 1 and/or a second derivative of the B0 in homogeneity map in at least one predetermined direction. The second derivative is calculated using a corrected voxel size in the at least one predetermined direction, wherein the corrected voxel size is calculated using a correction factor calculated from the derivative of the B0 inhomogeneity map.


