Classifications

• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B5/00—Measuring for diagnostic purposes; Identification of persons
  • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
  • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
  • G01R33/00—Arrangements or instruments for measuring magnetic variables
  • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
  • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
  • G01R33/443—Assessment of an electric or a magnetic field, e.g. spatial mapping, determination of a B0 drift or dosimetry
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
  • G01R33/00—Arrangements or instruments for measuring magnetic variables
  • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
  • G01R33/24—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
  • G01R33/243—Spatial mapping of the polarizing magnetic field
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
  • G01R33/00—Arrangements or instruments for measuring magnetic variables
  • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
  • G01R33/24—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
  • G01R33/246—Spatial mapping of the RF magnetic field B1
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
  • G01R33/00—Arrangements or instruments for measuring magnetic variables
  • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
  • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
  • G01R33/48—NMR imaging systems
  • G01R33/4816—NMR imaging of samples with ultrashort relaxation times such as solid samples, e.g. MRI using ultrashort TE [UTE], single point imaging, constant time imaging
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
  • G01R33/00—Arrangements or instruments for measuring magnetic variables
  • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
  • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
  • G01R33/48—NMR imaging systems
  • G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
  • G01R33/4824—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
  • G01R33/00—Arrangements or instruments for measuring magnetic variables
  • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
  • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
  • G01R33/48—NMR imaging systems
  • G01R33/4828—Resolving the MR signals of different chemical species, e.g. water-fat imaging
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
  • G01R33/00—Arrangements or instruments for measuring magnetic variables
  • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
  • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
  • G01R33/48—NMR imaging systems
  • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
  • G01R33/00—Arrangements or instruments for measuring magnetic variables
  • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
  • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
  • G01R33/48—NMR imaging systems
  • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
  • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
  • G01R33/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
  • G01R33/00—Arrangements or instruments for measuring magnetic variables
  • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
  • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
  • G01R33/48—NMR imaging systems
  • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
  • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
  • G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
  • G01R33/56563—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
  • G01R33/00—Arrangements or instruments for measuring magnetic variables
  • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
  • G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
  • G01R33/288—Provisions within MR facilities for enhancing safety during MR, e.g. reduction of the specific absorption rate [SAR], detection of ferromagnetic objects in the scanner room

Definitions

• the invention relates to magnetic resonance imaging, in particular to the processing of magnetic resonance images to account for geometric distortions of voxel size.
• a large static magnetic field is used by Magnetic Resonance Imaging (MRI) scanners to align the nuclear spins of atoms as part of the procedure for producing images within the body of a patient.
• This large static magnetic field is referred to as the BO field.
• Radio Frequency (RF) pulses generated by one or more transmitter coils cause perturbations to the local magnetic field, and RF signals emitted by the nuclear spins are detected by one or more receiver coils. These RF signals are used to construct the MR images.
• These coils can also be referred to as antennas.
• the transmitter and receiver coils can also be integrated into one or more transceiver coils that perform both functions. It is understood that the use of the term transceiver coil also refers to systems where separate transmitter and receiver coils are used.
• the transmitted RF field is referred to as the B 1 field.
• MRI scanners are able to construct images of either slices or volumes.
• a slice is a thin volume that is only one voxel thick.
• a voxel is a small volume over which the MR signal is averaged, and represents the resolution of the MR image.
• a voxel may also be referred to as a pixel herein.
• Dixon methods of magnetic resonance imaging include a family of techniques for producing separate water and lipid (fat) images.
• the various Dixon techniques such as, but not limited to, two-point Dixon methods, three-point Dixon methods, and multi-point Dixon methods are collectively referred to herein as Dixon techniques or methods.
• the terminology to describe the Dixon techniques is well known and has been the subject of many review articles and is present in standard texts on Magnetic Resonance Imaging. For example, the "Handbook of MRI Pulse Sequences" by Bernstein et ah, published by Elsevier Academic Press in 2004, contains a review of some Dixon techniques on pages 857 to 887.
• United States patent application US 20140184219 Al discloses calculating a zero phase estimation of a Bl phase map using a BO inhomogeneity map.
• the invention provides for a magnetic resonance imaging system, a method of operating a magnetic resonance imaging system and a computer program product in the independent claims. Embodiments are given in the dependent claims.
• aspects of the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment (including firmware, resident software, microcode, etc.) or an embodiment combining software and hardware aspects that may all generally be referred to herein as a "circuit,” “module” or “system.”
• aspects of the present invention may take the form of a computer program product embodied in one or more computer readable medium(s) having computer executable code embodied thereon.
• the computer readable medium may be a computer readable signal medium or a computer readable storage medium.
• a 'computer-readable storage medium' as used herein encompasses any tangible storage medium which may store instructions which are executable by a processor of a computing device.
• the computer-readable storage medium may be referred to as a computer-readable non-transitory storage medium.
• the computer-readable storage medium may also be referred to as a tangible computer readable medium.
• a computer-readable storage medium may also be able to store data which is able to be accessed by the processor of the computing device.
• Examples of computer- readable storage media include, but are not limited to: a floppy disk, a magnetic hard disk drive, a solid state hard disk, flash memory, a USB thumb drive, Random Access Memory (RAM), Read Only Memory (ROM), an optical disk, a magneto-optical disk, and the register file of the processor.
• Examples of optical disks include Compact Disks (CD) and Digital Versatile Disks (DVD), for example CD-ROM, CD-RW, CD-R, DVD-ROM, DVD-RW, or DVD-R disks.
• the term computer readable-storage medium also refers to various types of recording media capable of being accessed by the computer device via a network or communication link.
• a data may be retrieved over a modem, over the internet, or over a local area network.
• Computer executable code embodied on a computer readable medium may be transmitted using any appropriate medium, including but not limited to wireless, wire line, optical fiber cable, RF, etc., or any suitable combination of the foregoing.
• a computer readable signal medium may include a propagated data signal with computer executable code embodied therein, for example, in baseband or as part of a carrier wave. Such a propagated signal may take any of a variety of forms, including, but not limited to, electro-magnetic, optical, or any suitable combination thereof.
• a computer readable signal medium may be any computer readable medium that is not a computer readable storage medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device.
• 'Computer memory' or 'memory' is an example of a computer-readable storage medium.
• Computer memory is any memory which is directly accessible to a processor.
• 'Computer storage' or 'storage' is a further example of a computer-readable storage medium.
• Computer storage is any non- volatile computer-readable storage medium. In some embodiments computer storage may also be computer memory or vice versa.
• a 'processor' as used herein encompasses an electronic component which is able to execute a program or machine executable instruction or computer executable code.
• References to the computing device comprising "a processor” should be interpreted as possibly containing more than one processor or processing core.
• the processor may for instance be a multi-core processor.
• a processor may also refer to a collection of processors within a single computer system or distributed amongst multiple computer systems.
• the term computing device should also be interpreted to possibly refer to a collection or network of computing devices each comprising a processor or processors.
• the computer executable code may be executed by multiple processors that may be within the same computing device or which may even be distributed across multiple computing devices.
• Computer executable code may comprise machine executable instructions or a program which causes a processor to perform an aspect of the present invention.
• Computer executable code for carrying out operations for aspects of the present invention may be written in any combination of one or more programming languages, including an object oriented programming language such as Java, Smalltalk, C++ or the like and conventional procedural programming languages, such as the C programming language or similar programming languages and compiled into machine executable instructions.
• the computer executable code may be in the form of a high level language or in a pre-compiled form and be used in conjunction with an interpreter which generates the machine executable instructions on the fly.
• the computer executable code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server.
• the remote computer may be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet Service Provider).
• These computer program instructions may also be stored in a computer readable medium that can direct a computer, other programmable data processing apparatus, or other devices to function in a particular manner, such that the instructions stored in the computer readable medium produce an article of manufacture including instructions which implement the function/act specified in the flowchart and/or block diagram block or blocks.
• the computer program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other devices to cause a series of operational steps to be performed on the computer, other programmable apparatus or other devices to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide processes for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks.
• a 'user interface' as used herein is an interface which allows a user or operator to interact with a computer or computer system.
• a 'user interface' may also be referred to as a 'human interface device.
• a user interface may provide information or data to the operator and/or receive information or data from the operator.
• a user interface may enable input from an operator to be received by the computer and may provide output to the user from the computer.
• the user interface may allow an operator to control or manipulate a computer and the interface may allow the computer indicate the effects of the operator's control or manipulation.
• the display of data or information on a display or a graphical user interface is an example of providing information to an operator.
• the receiving of data through a keyboard, mouse, trackball, touchpad, pointing stick, graphics tablet, joystick, gamepad, webcam, headset, pedals, wired glove, remote control, and accelerometer are all examples of user interface components which enable the receiving of information or data from an operator.
• a 'hardware interface' as used herein encompasses an interface which enables the processor of a computer system to interact with and/or control an external computing device and/or apparatus.
• a hardware interface may allow a processor to send control signals or instructions to an external computing device and/or apparatus.
• a hardware interface may also enable a processor to exchange data with an external computing device and/or apparatus. Examples of a hardware interface include, but are not limited to: a universal serial bus, IEEE 1394 port, parallel port, IEEE 1284 port, serial port, RS-232 port, IEEE-488 port, bluetooth connection, wireless local area network connection, TCP/IP connection, ethernet connection, control voltage interface, MIDI interface, analog input interface, and digital input interface.
• a 'display' or 'display device' as used herein encompasses an output device or a user interface adapted for displaying images or data.
• a display may output visual, audio, and or tactile data. Examples of a display include, but are not limited to: a computer monitor, a television screen, a touch screen, tactile electronic display, Braille screen,
• Cathode ray tube (CRT), Storage tube, Bi-stable display, Electronic paper, Vector display, Flat panel display, Vacuum fluorescent display (VF), Light-emitting diode (LED) displays, Electroluminescent display (ELD), Plasma display panels (PDP), Liquid crystal display (LCD), Organic light-emitting diode displays (OLED), a projector, and Head- mounted display.
• CTR Cathode ray tube
• Storage tube Bi-stable display
• Electronic paper Electronic paper
• Vector display Flat panel display
• VF Vacuum fluorescent display
• LED Light-emitting diode
• ELD Electroluminescent display
• PDP Plasma display panels
• LCD Liquid crystal display
• OLED Organic light-emitting diode displays
• Magnetic Resonance (MR) data is defined herein as being the recorded measurements of radio frequency signals emitted by atomic spins using the antenna of a magnetic resonance apparatus during a magnetic resonance imaging scan.
• Magnetic resonance data is an example of medical image data.
• a Magnetic Resonance (MR) image is defined herein as being the reconstructed two or three dimensional visualization of anatomic data contained within the magnetic resonance imaging data. This visualization can be performed using a computer.
• a magnetic resonance imaging system for acquiring magnetic resonance data from a subject within an imaging zone.
• the magnetic resonance imaging system comprises a memory for storing machine-executable instructions and Dixon pulse sequence data.
• Pulse sequence data as used herein encompasses data which comprises instructions for controlling a magnetic resonance imaging system or data which may be readily converted into such instructions.
• a pulse sequence is commonly a timing diagram which illustrates what components of a magnetic resonance imaging system do at a particular point in time. Such a timing diagram or data descriptive of a timing diagram could be converted into instructions for controlling the magnetic resonance imaging system.
• 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, n is greater than or equal to 2.
• n is greater than or equal to 2.
• the magnetic resonance imaging system further comprises a processor for controlling the magnetic resonance imaging system.
• Execution of the machine-executable instructions causes the processor to acquire magnetic resonance data by controlling the magnetic resonance imaging system using the Dixon pulse sequence data.
• Execution of the machine-executable instructions causes the processor to estimate a BO inhomogeneity map and to estimate an estimated Bl phase map by analyzing the magnetic resonance data with an n-point Dixon method.
• n-point Dixon methods fit a signal equation, which describes the evolution of the MR signal over time, to the Dixon pulse sequence data.
• the signal equation typically includes the contributions of water and fat to the MR signal as unknowns, as well as a phase offset over time common to water and fat.
• the phase offset over time between water and fat due to the difference in chemical shift between water and fat is usually assumed to be known, as is the amplitude variation of the fat signal over time due to de- and rephasing of the spectral components of fat.
• the contributions of water and fat to the MR signal may either be considered as real variables, in which case a further unknown is introduced to characterize the phase offset common to water and fat at the time of excitation, i.e. the center of the RF pulse, or as complex variables.
• the further unknown is herein called zero echo time Bl phase map ⁇ .
• the phase of the complex variable modeling the contribution of water to the MR signal herein denoted by ⁇
• ⁇ the phase of the complex variable modeling the contribution of water to the MR signal
• ⁇ the phase offset per unit time common to water and fat
• t the first echo time, at which the Dixon pulse sequence acquired data.
• phase of the complex variable modeling the contribution of fat to the MR signal may be rewound to zero echo time, taking the additional phase offset over time between water and fat due to the difference in chemical shift between water and fat into account.
• both may be performed and the results may be suitably averaged.
• the basic equation for EPT is given by (cf. Katscher U et al, IEEE Trans Med Imag, 28 (2009) 1365)
• hyperthermia treatment planning or SAR modelling hyperthermia treatment planning or SAR modelling.
• it can be used as starting point for iterative EPT algorithms as known per se from Lee JS et al., ISMRM 2013;21 :4183 and Balidemaj E et al, ISMRM 2013;21 :4185 .
• Execution of the machine-executable instructions further causes the processor to calculate at least one calculated electrical conductivity map using the zero echo time Bl phase map.
• This example may have the benefit that the electrical conductivity map can be determined at the same time that the n-point Dixon method is performed.
• execution of the machine-executable instructions further causes the processor to calculate a fat image and a water image when analyzing the magnetic resonance data according to the n-point Dixon method.
• execution of the instructions further causes the processor to identify fat regions within the subject by segmenting the fat image.
• Execution of the instructions further causes the processor to identify water regions within the subject by segmenting the water image.
• This segmentation could for instance be done in different ways.
• a model such as a deformable-shaped body model or an anatomical atlas could be used to assist in the segmentation or a raw analysis of the voxels of the fat image and the water image could be used.
• the fat image and the water image could be thresholded to produce the data about where the fat regions and the water regions are.
• the fat regions may include regions that are partially filled with water. Likewise the water regions may partially include fatty tissue.
• execution of the instructions further causes the processor to divide the subject into multiple differential kernel regions using the fat regions and the water regions.
• Execution of the instructions further causes the processor to determine boundary conditions between the multiple different kernel regions.
• the at least one calculated electrical conductivity map is calculated using the boundary conditions between the multiple differential kernel regions. This example may be particularly beneficial because the data from the fat and water images is used to produce boundary conditions which may improve the quality or accuracy of the calculation of the calculated electrical conductivity map.
• the fat and water regions may have different electrical properties. By explicitly making boundary conditions between these regions the at least one calculated electrical conductivity map can be calculated more accurately.
• the at least one calculated electrical conductivity map is calculated using a kernel that solves differential equations.
• execution of the instructions further causes the processor to calculate an estimated electrical conductivity map and an estimated permittivity map using the fat regions and the water regions.
• the estimated electrical conductivity map and the estimated permittivity map are only calculated using knowledge of the properties of fat and water. Calculating these two maps may be useful because they can be used to directly calculate other properties or may be used as starting conditions for solving the differential equations or for other numerical methods.
• the magnetic resonance imaging system further comprises an electromagnetic tissue heating system.
• An electromagnetic tissue heating system as used herein encompasses a system which is used to heat the tissue of a subject using
• Execution of the instructions further causes the processor to estimate spatially dependent heating of the subject using the estimated electrical conductivity map and the estimated permittivity map. Using the estimated electrical conductivity map and the estimated permittivity map to get the heating of the subject may be useful for treatment planning and also to avoid overheating portions of the subject.
• the spatially dependent heating of the subject is estimated using the at least one calculated electrical conductivity map and the estimated permittivity map. This may be beneficial because with knowledge of the zero echo time Bl phase map only the electrical conductivity can be calculated.
• the tissue heating system is a microwave tissue heating system.
• the tissue heating system is a radio-frequency tissue heating system.
• execution of the instructions further causes the processor to calculate the at least one electrical conductivity map using the zero echo time Bl map with an iterative solver.
• the iterative solver may be used to solve the solution for the at least one electrical conductivity map and an iterative method for solving the differential equations.
• the iterative solver is configured for using the estimated electrical conductivity map and the estimated permittivity map at least partially to determine an initial solution for initializing the iterative solver.
• execution of the instructions further causes the processor to identify the water regions by selecting regions with a fat-to-water ratio within a second predefined range. As with the fat regions this could be done by thresholding or by plotting a fat-to-water ratio within a particular range of values. The values of the first predefined range are greater than the values of the second predefined range.
• execution of the instructions further causes the processor to calculate a combined electrical conductivity map by combining the water electrical conductivity map and the fat electrical conductivity map. In this example the two images are combined to create a combined electrical conductivity map.
• an example provides for a method of operating the magnetic resonance imaging system for acquiring magnetic resonance data from a subject within an imaging zone.
• the method comprises the step of acquiring magnetic resonance data by controlling the magnetic resonance imaging system using the Dixon pulse sequence data.
• 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, n is greater than or equal to 2.
• the method further comprises the step of estimating a BO inhomogeneity map and estimating an estimated Bl phase map by analyzing the magnetic resonance data according to the n-point Dixon method.
• the method further comprises the step of calculating a zero echo time Bl phase map by interpolating the estimated Bl phase map to an echo time of zero using the BO inhomogeneity map.
• the method further comprises the step of calculating at least one calculated electrical conductivity map using the zero echo time Bl phase map.
• an example provides for a computer program product comprising machine-executable instructions for execution by a processor controlling the magnetic resonance imaging system for acquiring magnetic resonance data from a subject within an imaging zone.
• Execution of the machine-executable instructions causes the processor to acquire magnetic resonance data by controlling the magnetic resonance imaging system using Dixon pulse sequence data.
• the Dixon pulse sequence data comprises instructions for controlling the magnetic resonance imaging system to acquire magnetic resonance data according to an /? -point Dixon method, n is equal to or greater than 2.
• Execution of the instructions further causes the processor to estimate a BO inhomogeneity map and estimate an estimated Bl phase map by analyzing the magnetic resonance data according to the n-point Dixon method. Execution of the instructions further causes the processor to calculate a zero echo time Bl phase map by interpolating the estimated Bl phase map to an echo time of zero using the BO inhomogeneity map. Execution of the machine-executable instructions further causes the processor to calculate at least one calculated electrical conductivity map using the zero echo time Bl phase map.
• the presence of metallic objects may cause magnetic field inhomogeneities which lead to geometric distortions of the imaged voxels. This geometric distortion of the voxels may therefore lead to errors in values calculated with differential equations. For example, when performing electrical properties tomography the electrical conductivity and permittivity are calculated using differential equations. However, it is possible to correct for or account for the magnetic field inhomogeneities directly within the calculation of the derivatives.
• the invention provides for a magnetic resonance imaging system for acquiring magnetic resonance data from a subject within an imaging zone.
• the magnetic resonance imaging system comprises a memory for storing machine-executable instructions and pulse sequence data.
• 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 for controlling the magnetic resonance imaging system.
• Execution of the machine-executable instructions causes the processor to acquire the magnetic resonance data by controlling the magnetic resonance imaging system with the pulse sequence data.
• the pulse sequence data may comprise more than one pulse sequence.
• the magnetic resonance data may therefore be formed from more than one dataset.
• Execution of the machine-executable instructions further cause the processor to calculate a BO inhomogeneity map by analyzing the magnetic resonance data according to the magnetic resonance imaging method.
• Execution of the instructions further cause the processor to calculate a Bl phase map and/or a Bl amplitude map by analyzing the magnetic resonance data according to the magnetic resonance imaging method.
• Execution of the machine-executable instructions further cause the processor to calculate a second derivative of the Bl phase map and/or calculate a second derivative of the Bl magnitude map and/or calculate a second derivative of the BO inhomogeneity map in at least one predetermined direction.
• the second derivative is calculated using a corrected voxel size in the at least one predetermined direction.
• the corrected voxel size is calculated using a correction factor calculated from the derivative of the BO inhomogeneity map.
• This embodiment may be beneficial because it may provide for a means of calculating the second derivative of the Bl phase map, the second derivative of the Bl magnitude map, or even the second derivative of the BO inhomogeneity map more accurately.
• execution of the machine-executable instructions further cause the processor to calculate a zero echo time Bl phase map by interpolating the Bl phase map to an echo time of zero using the BO inhomogeneity map.
• the second derivative of the Bl phase map is calculated using a zero echo time Bl phase map.
• the pulse sequence data is Dixon pulse sequence data.
• the magnetic resonance imaging method is an n-point Dixon method, n is greater than or equal to 2.
• the BO inhomogeneity map and the Bl phase map are estimated by analyzing the magnetic resonance data according to the n-point Dixon method.
• execution of the machine-executable instructions further cause the processor to calculate a fat image and a water image when analyzing the magnetic resonance data according to the n-point Dixon method.
• Execution of the machine- executable instructions further causes the processor to identify fat regions within the subject by segmenting the fat image.
• Execution of the machine-executable instructions further causes the processor to identify water regions within the subject by segmenting the water image.
• Execution of the machine-executable instructions further cause the processor to calculate an estimated electrical conductivity map and/or an estimated permittivity map using the fat regions and the water regions.
• the magnetic resonance imaging system further comprises an electromagnetic tissue heating system. Execution of the instructions further cause the processor to estimate spatially dependent heating of the subject using the estimated permittivity map. Execution of the machine-executable instructions further cause the processor to further estimate the heating of the subject using the estimated electrical conductivity map and/or the at least one calculated electrical conductivity map.
• the pulse sequence data comprises a multi-echo pulse sequence for measuring the BO map.
• the pulse sequence data comprises BO mapping pulse sequence data.
• the pulse sequence data further comprises a B 1 magnitude measuring pulse sequence for measuring the Bl magnitude map.
• the pulse sequence data further comprises Bl magnitude mapping pulse sequence data.
• the pulse sequence data further comprises a Bl phase measuring pulse sequence for measuring the Bl phase map.
• the pulse sequence data comprises Bl phase mapping pulse sequence data.
• the pulse sequence data further comprises both a Bl magnitude measuring pulse sequence and the Bl phase measuring pulse sequence.
• both the Bl magnitude and Bl phase are measured that are acquired in separate scans of the magnetic resonance data.
• the multi-echo pulse sequence comprises instructions for controlling the magnetic resonance imaging system to acquire magnetic resonance data according to any one of the following magnetic resonance imaging methods: an n-point Dixon method and a multi-echo pulse sequence method.
• the Bl magnitude measuring pulse sequence comprises instructions for controlling the magnetic resonance imaging system to acquire magnetic resonance data according to any one of the following magnetic resonance imaging methods: an actual flip-angle imaging (AFI) magnetic resonance imaging method, a dual refocusing echo acquisition mode (DREAM) magnetic resonance imaging method, and a Bloch-Siegert shift magnetic resonance imaging method.
• AFI actual flip-angle imaging
• DREAM dual refocusing echo acquisition mode
• Bloch-Siegert shift magnetic resonance imaging method Bloch-Siegert shift magnetic resonance imaging method.
• the B 1 phase measuring pulse sequence comprises instructions for controlling the magnetic resonance imaging system to acquire magnetic resonance data according to any one of the following magnetic resonance imaging methods: a spin echo based magnetic resonance imaging method and a balanced gradient echo magnetic resonance imaging method.
• execution of the machine-executable instructions further causes the processor to calculate the at least one calculated electrical conductivity map using the second derivative of the Bl phase map in the at least one predetermined direction.
• execution of the machine-executable instructions further cause the processor to calculate at least one susceptibility map using the second derivative of the BO inhomogeneity map in the at least one predetermined direction according to a quantitative susceptibility mapping method.
• the pulse sequence data specifies a readout gradient for each of the at least one predetermined direction.
• the corrected voxel size in each of the at least one predetermined direction is .
• x is one
• ⁇ is any one of the following: the Bl phase map, the Bl magnitude map, and the BO inhomogeneity map.
• the invention provides for a method of operating a magnetic resonance imaging system for acquiring magnetic resonance data from a subject within the imaging zone.
• the method comprises the step of acquiring the magnetic resonance data by controlling the magnetic resonance imaging system with pulse sequence data.
• 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 method further comprises the step of calculating a BO inhomogeneity map by analyzing the magnetic resonance data according to the magnetic resonance imaging method.
• the method further comprises the step of calculating a B 1 phase map and/or a Bl amplitude map by analyzing the magnetic resonance data according to the magnetic resonance imaging method.
• the method further comprises the step of calculating a second derivative of the Bl phase map and/or a second derivative of the Bl magnitude map and/or a second derivative of the BO inhomogeneity map in the at least one predetermined direction.
• the second derivative is calculated using a corrected voxel size in the at least one
• the corrected voxel size is calculated using a correction factor calculated from the derivative of the BO inhomogeneity map in the at least one predetermined direction.
• the invention provides for a computer program product comprising machine-executable instructions for execution by a processor controlling a magnetic resonance imaging system for acquiring magnetic resonance data from a subject within an imaging zone.
• Execution of the machine-executable instructions causes the processor to acquire the magnetic resonance data by controlling the magnetic resonance imaging system with the pulse sequence data.
• 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.
• Execution of the machine-executable instructions further cause the processor to calculate a BO inhomogeneity map by analyzing the magnetic resonance data according to the magnetic resonance imaging method.
• Execution of the machine-executable instructions further cause the processor to calculate a Bl phase map and/or a Bl amplitude map by analyzing the magnetic resonance data according to the magnetic resonance imaging method.
• Execution of the machine- executable instructions further cause the processor to calculate a second derivative of the Bl phase map and/or a second derivative of the Bl magnitude map and/or second derivative of the BO inhomogeneity map in the at least one predetermined direction.
• the second derivative is calculated using the corrected voxel size in the at least one predetermined direction.
• the corrected voxel size is calculated using a correction factor calculated from the derivative of the BO inhomogeneity map.
• Fig. 1 illustrates an example of a magnetic resonance imaging system
• Fig. 2 shows a flow chart of a method of operating the magnetic resonance imaging system of claim 1 ;
• Fig. 3 illustrates a further example of a magnetic resonance imaging system
• Fig. 4 shows the magnitude image from an mDixon image
• Fig. 5 shows the water separated image of the mDixon image of Fig. 4;
• Fig. 6 shows the fat separated image of the mDixon image of Fig. 4;
• Fig. 7 shows a segmentation of Fig. 4 using a Laplace operator
• Fig. 8 shows a segmentation of Fig. 6 using a Laplace operator
• Fig. 9 shows a segmentation of Fig. 5 using a Laplace operator
• Fig. 10 shows a combination of images 8 and 9;
• Fig. 11 shows a further example of a magnetic resonance imaging system
• Fig. 12 shows a flow chart of a method of operating the magnetic resonance imaging system of claim 11 ;
• Fig. 13 shows four images; and Fig. 14 shows a plot of the relative error of the second derivative of the phase along the readout direction which is averaged over 10 voxels along the x-direction in one slice for an experiment.
• Fig. 1 shows an example of a magnetic resonance imaging system 100 with a magnet 104.
• the magnet 104 is a superconducting cylindrical type magnet 104 with a bore 106 through it.
• the use of different types of magnets is also possible; for instance it is also possible to use both a split cylindrical magnet and a so called open magnet.
• a split cylindrical magnet is similar to a standard cylindrical magnet, except that the cryostat has been split into two sections to allow access to the iso-plane of the magnet, such magnets may for instance be used in conjunction with charged particle beam therapy.
• An open magnet has two magnet sections, one above the other with a space in-between that is large enough to receive a subject: the arrangement of the two sections area similar to that of a Helmholtz coil.
• Open magnets are popular, because the subject is less confined. Inside the cryostat of the cylindrical magnet, there is a collection of superconducting coils. Within the bore 106 of the cylindrical magnet 104 there is an imaging zone 108 where the magnetic field is strong and uniform enough to perform magnetic resonance imaging.
• the magnetic field gradient coils 110 are connected to a magnetic field gradient coil power supply 112.
• the magnetic field gradient coils 110 are intended to be representative.
• magnetic field gradient coils 110 contain three separate sets of coils for spatially encoding in three orthogonal spatial directions.
• a magnetic field gradient power supply supplies current to the magnetic field gradient coils. The current supplied to the magnetic field gradient coils 110 is controlled as a function of time and may be ramped or pulsed.
• a radio-frequency coil 114 Adjacent to the imaging zone 108 is a radio-frequency coil 114 for manipulating the orientation of magnetic spins within the imaging zone 108 and for receiving radio transmissions from spins also within the imaging zone 108.
• the radio frequency antenna may contain multiple coil elements.
• the radio frequency antenna may also be referred to as a channel or antenna.
• the radio -frequency coil 114 is connected to a radio frequency transceiver 116.
• the radio-frequency coil 114 and radio frequency transceiver 116 may be replaced by separate transmit and receive coils and a separate transmitter and receiver. It is understood that the radio -frequency coil 114 and the radio frequency transceiver 116 are representative.
• the radio-frequency coil 114 is intended to also represent a dedicated transmit antenna and a dedicated receive antenna.
• the transceiver 116 may also represent a separate transmitter and receiver.
• the radio-frequency coil 114 may also have multiple receive/transmit elements and the radio frequency transceiver 116 may have multiple receive/transmit
• a subject support 120 which is attached to an optional actuator 122 that is able to move the subject support and the subject 118 through the imaging zone 108.
• the transceiver 116, the magnetic field gradient coil power supply 112 and the actuator 122 are all seen as being connected to a hardware interface 128 of computer system 126.
• the contents of the computer storage 134 and the computer memory 136 may be interchangeable. In some examples the contents of the computer storage 134 may be duplicated in the computer memory 136.
• the computer storage 134 is shown as containing the Dixon pulse sequence data.
• the computer storage 134 is further shown as containing magnetic resonance data 142 that has been acquired by controlling the magnetic resonance imaging system 100 with the Dixon pulse sequence data 140.
• the computer storage 134 is further shown as containing a fat image 144 and a water image 146 that have been reconstructed from the magnetic resonance data 142 according to the Dixon method.
• the computer storage 134 is further shown as containing a B0 inhomogeneity map 148 and an estimated Bl phase map 150 that have also been calculated using the Dixon method.
• the computer storage 134 is further shown as containing a zero echo time Bl phase map 152 that has been calculated by interpolating the estimated B 1 phase map to an echo time of 0 using the B0 inhomogeneity map.
• the computer storage 134 is further shown as containing electrical conductivity maps 154 that have been calculated using the zero echo time Bl phase map.
• the computer storage 134 is further shown as containing fat region location data 156 and the water region location data 158 that have been determined from the fat image 144 and the water image 146 respectively. These may even be determined for instance by performing image processing on the fat image 144 and the water image 146.
• the computer memory 136 is shown as containing a control module 160.
• the control module 160 contains computer executable code which enables the processor 130 to control the magnetic resonance imaging system 100. For instance the control module 160 may enable the processor 130 to control the magnetic resonance imaging system 100 with the Dixon pulse sequence data 140 to acquire the magnetic resonance data 142.
• the computer memory 136 is further shown as containing an image reconstruction module 162.
• the image reconstruction module 162 enables the processor 130 to process the magnetic resonance data 142 into the fat image 144, the water image 146, the B0 inhomogeneity map 148, and the estimated Bl phase map 150.
• the image reconstruction module essentially enables the processor 130 to perform the data analysis aspects of the Dixon method.
• Fig. 2 shows a flowchart which illustrates a method of operating the magnetic resonance imaging system 100 of Fig. 1.
• the magnetic resonance data 142 is acquired by controlling the magnetic resonance imaging system with the Dixon pulse sequence data 140.
• the B0 inhomogeneity map 148 and the estimated Bl phase map 150 are calculated by analyzing the magnetic resonance data according to an n- point Dixon method.
• a zero echo time Bl phase map 152 is calculated by interpolating the estimated Bl phase map 150 to an echo time of 0 using the B0
• step 206 the at least one electrical conductivity map 154 is calculated using the zero echo time Bl phase map 152.
• Fig. 3 shows a further example of a magnetic resonance imaging system 300.
• the magnetic resonance imaging system 300 also includes a tissue heating system formed by an antenna 302 adjacent to the subject 118 and by a radio -frequency transmitter 304.
• the combination of the antenna 302 and the radio -frequency transmitter 304 are exemplary. For example this could be replaced with a system that generates microwave or other electromagnetic radiation for heating the subject 118.
• the subject 118 is shown as having a target zone 306 that is desired to be heated within the subject 118.
• the computer storage 134 is shown as containing a treatment plan 312.
• the treatment plan 312 may be descriptive of the internal structure of the subject 118 and contain data which enables identification or location of the target zone 306.
• Computer storage 134 further contains a set of heating system control commands 314 that have been generated using the treatment plan 312.
• the heating system control commands 314 contain commands which enable the processor 130 to control the operation and function of the heating system 302, 304.
• the computer memory 136 is further shown as containing a heating system control generation module 316.
• the heating system control generation module 316 contains computer executable code which enables the processor 130 to generate the heating system control commands 314 from the treatment plan 312 and possibly from magnetic resonance data acquired by the magnetic resonance imaging system 300.
• the radio-frequency heating system comprises an antenna 302 and a radio- frequency transmitter 304.
• the antenna 302 is in the vicinity of target zone 306.
• Radio- frequency energy generated by the transmitter 304 and radiated by the antenna 302 is used to selectively heat the target zone 306.
• the radio-frequency transmitter 304 is shown as being connected to the hardware interface 128.
• the computer storage 134 is shown as having an estimated permittivity map 310 that was calculated by knowing the permittivity properties of fat and water and then using the fat region location data 156 and the water region location data 158 to calculate the estimated electrical permittivity.
• the heating system control generation module 316 can use the estimated permittivity map 310 and the at least one electrical conductivity map 154 to estimate heating of the subject 118 by the antenna 302. This may enable a more accurate calculation or determination of the heating system control commands 314.
• EPT Electric Properties Tomography
• the human body In general the human body consists roughly of 65% water, 10% fat (sometimes more) and 20% proteins and minerals, where the latter two are often difficult to detect directly by proton MR. Water-rich tissue has obviously a significantly higher conductivity than fat, which is often seen as a kind of electrical isolator. Fat and water contributions to the received MR signal can be separated based on a chemical shift encoded acquisition using so- called Dixon methods. The chemical shift encoding is usually achieved by repeated measurements at different echo times, and the fat-water separation commonly involves an estimation of the underlying main field (BO) inhomogeneity.
• BO main field
• EPT is hampered by a couple of physical/technical drawbacks. For this invention, three of these drawbacks are considered.
• EPT requires phase data purely related to RF penetration (i.e., free of contributions from BO inhomogeneities), which usually is fulfilled sufficiently only for spin- echo based sequences. Thus, if a patient exam contains only field-echo based sequences, EPT requires additional scan time using dedicated MRI sequences.
• Iterative EPT reconstruction algorithms require a suitable starting point. This is typically realized by applying standard EPT (solving the inverse problem) or by applying literature values of the electric properties to the compartments of the individual patient. Both methods are time-consuming.
• the numerical differentiation kernel required for EPT should not contain voxels from different compartments, requiring suitable image segmentation. Image segmentation is sometimes hampered by low contrast between tissue compartments.
• Examples may have one or more of the following features
• the required RF phase purely related to RF penetration can be extracted from a Dixon scan without further scans. Furthermore, the obtained BO map can be used to apply EPT also to other field-echo based sequences.
• From the Dixon scan relative concentrations of fat Cp(r) and water c ⁇ r) are obtained. In some of the Dixon applications some care is necessary to make these numbers really quantitative (see fat fraction quantification), because the signal intensity might be colored by the sequence parameters used.
• Image segmentation is required for EPT to ensure that numerical differentiation kernels do not cross boundaries of tissue compartments with different electric properties. This case is not covered by Eqs. (1,2) and would lead to strong oscillatory artefacts in the reconstructed electric properties along the compartment boundaries.
• Image segmentation is typically based on differences in the magnitude of the
• MR signal between the compartments to be segmented.
• this difference is not always guaranteed, since electric properties have no direct impact on the MR signal magnitude.
• Fat and water images can support image segmentation, since they provide a different contrast than standard (composite) MR images. It is also possible to base segmentation on multiple images, combining areas of locally highest contrast from different images. Alternatively, image segmentation can be taken into account by applying two separate EPT reconstructions: the first reconstruction based on the water image, the second reconstruction based on the fat image.
• This procedure has the additional advantage that potential chemical shift artefacts do not deteriorate the EPT reconstruction.
• the phantom is a bottle with saline (in the lower part of bottle) and oil (in the upper part of bottle).
• the magnitude image of echo 1 is shown in Fig. 4, the water image in Fig. 5, and the fat image in Fig. 6.
• Figs. 4, 5 and 6 show results from performing a two-point Dixon method on a phantom.
• the phantom has an upper part made of a fat-like phantom labeled 402 in Fig. 4 and a lower saline portion labeled 400 in Fig. 4.
• Fig. 4 shows the magnitude image from an mDixon image.
• Fig. 5 shows the water separated image 146.
• Fig. 6 shows the corresponding fat image 144.
• Fig. 7 shows a segmentation of Fig. 4 using a Laplace operator. It can be seen in Fig. 7 that a clear delineation of the fat 402 and saline 400 regions is not shown in Fig. 7.
• Fig. 8 shows a segmentation of Fig. 6 using the Laplace operator. As 6 is a fat image the segmentation in Fig. 8 indicates a fat region 156.
• Fig. 9 shows a segmentation of Fig. 5 using the Laplace operator.
• Fig. 5 is a water image 146 the segmentation in Fig. 9 indicates a water region 158.
• Fig. 10 shows a combination of Figs. 8 and 9.
• the combination of Figs. 8 and 9 shows a region which is identified as a fat region 156 and a water region 158.
• Fig. 10 illustrates that the segmentation of images from a Dixon method can be used to identify fat and water regions within an image. This will be particularly useful in solving the differential equations for determining the electrical conductivity.
• the boundary region 1000 can be fed to the differential equation solver or kernel so that the appropriate boundary conditions between the fat region 156 and the water region 158 can be used.
• Fig. 11 shows an example of a magnetic resonance imaging system 1100 that is similar to that shown in Fig. 1 and Fig. 3.
• the magnetic resonance imaging systems shown in Figs. 1, 3 and 11 may have their features combined.
• the computer storage 134 is shown as containing a pulse sequence data 1102 for controlling the magnetic resonance imaging system 1100. In some instances the pulse sequence data 1102 may be identical with the pulse sequence data 140 of Fig. 1.
• the computer storage 134 is further shown as containing magnetic resonance data 1104. In some cases the magnetic resonance data 1104 may be identical with the magnetic resonance data 142 of Fig. 1.
• the computer storage 134 is further shown as containing a B0 inhomogeneity map 148 that was determined using the magnetic resonance data 1104.
• the computer storage 134 is further shown as containing an estimated Bl phase map 150 and a Bl amplitude map 1106. Both the Bl phase map 150 and the Bl amplitude map 1106 are also calculated or derived from the magnetic resonance data 1104. The estimated Bl phase map 150 and the Bl amplitude map 1106 may not be present in all embodiments.
• the computer memory 136 is shown as containing a control module 160, an image reconstruction module 162, and an image processing module 164. These modules are as described in Fig. 1 and/or Fig. 3.
• the image processing module 164 for instance may be used to calculate the correction factor for calculating the second derivative and also calculating the value of the second derivatives.
• the computer storage 134 is further shown as containing a corrected voxel size 1108 and a second derivative 1110.
• the second derivative may be representative of a B0 inhomogeneity map second derivative, a Bl phase map second derivative and/or a Bl amplitude map second derivative.
• the corrected voxel size 1108 and the second derivative 1110 may be calculated by numerical code which may be a separate module or may for instance be part of an image processing module 164.
• the computer storage 134 or computer memory 136 may contain additional data and/or numerical algorithms for performing such things as quantitative susceptibility mapping or electric properties tomography.
• Fig. 12 shows a flowchart which illustrates an example of a method of operating the magnetic resonance system 1100 of Fig. 11.
• the magnetic resonance data 1104 is acquired by controlling the magnetic resonance imaging system 1100 with the pulse sequence data 1102.
• a B0 inhomogeneity map 148 is calculated by analyzing the magnetic resonance data 1104 according to a magnetic resonance imaging method.
• a Bl phase map 150 and/or a Bl amplitude map 1106 is calculated by analyzing the magnetic resonance data 142 according to the magnetic resonance imaging method.
• a second derivative 1110 of the Bl phase map 150 and/or a second derivative of the Bl magnitude map 1106 and/or a second derivative of the B0 inhomogeneity map 148 is calculated in at least one predetermined direction.
• the second derivative is calculated using the corrected voxel size 1108 in the at least one predetermined direction.
• the corrected voxel size is calculated using a correction factor calculated from the derivative of the BO inhomogeneity map 148.
• Geometric distortions due to static inhomogeneities of the magnetic field adversely affect any image processing procedure which relies on spatial derivatives, e.g. Quantitative Susceptibility Mapping (QSM) or Electric Properties Tomography (EPT).
• QSM Quantitative Susceptibility Mapping
• EPT Electric Properties Tomography
• Magnetic field distortions caused by the object placed in the MR scanner lead to additional, spatially varying magnetic field gradients. These additional gradients lead to geometric distortions of the reconstructed images in which the targeted nominal voxel size (as shown in scanner GUI) is different from the true (physical) voxel size.
• the distortions of the magnetic field can be measured by acquiring a B0-map (i.e., off-resonance map) using a suitable MR sequence.
• Susceptibility Mapping or Electric Properties Tomography (EPT)
• EPT Electric Properties Tomography
• a physical model expressed as a differential equation.
• derivatives of a suitable MR image are calculated using finite differences into which the voxel sizes enters. It is therefore important that the geometric distortions are taken into account to avoid systematic errors, which can exceed 50%>.
• a correction factor can be calculated for each spatial position relating the nominal voxel size to the true voxel size (in readout direction determined solely by G R ).
• the true voxel size at a given location in space, dx ⁇ (r) is given by:
• ⁇ ⁇ ( - Nominal ) ⁇ 2 ⁇ ( ) + ⁇ ( + )
• both the error in the voxel size and the order of the derivative enter linearly into the result.
• Fig. 13 shows four images 1300, 1302, 1304, and 1306.
• the first image 1300 shows an axial slice of a phantom showing the measurement of the field map.
• Image 1302 shows the relative error of the voxel size for the voxels shown in image 1300.
• the third image 1304 shows the phase of the image 1300 at the echo time of 19.4 ms.
• the fourth image 1306 shows the relative error of the second derivative of the phase. In all four images the vertical axis is in the readout direction.
• Fig. 14 shows a plot of the relative error of the second derivative of the phase along the readout direction which is averaged over 10 voxels along the x-direction in one slice.
• a computer program may be stored/distributed on a suitable medium, such as an optical storage medium or a solid-state medium supplied together with or as part of other hardware, but may also be distributed in other forms, such as via the Internet or other wired or wireless telecommunication systems. Any reference signs in the claims should not be construed as limiting the scope.
• 166 differential equation kernel module 200 acquire magnetic resonance data by controlling the magnetic resonance imaging system using the Dixon pulse sequence data
• 1200 acquire the magnetic resonance data by controlling the magnetic resonance imaging system with pulse sequence data
• magnetic resonance data according to the magnetic resonance imaging method 1206 calculate a second derivative of the Bl phase map and/or a second derivative of the Bl magnitude map and/or a second derivative of the B0 inhomogeneity map in at least one predetermined direction 1300 magnetic resonance image of a phantom which shows an axial slice of a phantom showing the measurement of the field map.

Landscapes

• Physics & Mathematics (AREA)
• Condensed Matter Physics & Semiconductors (AREA)
• General Physics & Mathematics (AREA)
• High Energy & Nuclear Physics (AREA)
• Health & Medical Sciences (AREA)
• Engineering & Computer Science (AREA)
• Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
• Signal Processing (AREA)
• General Health & Medical Sciences (AREA)
• Radiology & Medical Imaging (AREA)
• Life Sciences & Earth Sciences (AREA)
• Computer Vision & Pattern Recognition (AREA)
• Artificial Intelligence (AREA)
• Biophysics (AREA)
• Pathology (AREA)
• Biomedical Technology (AREA)
• Heart & Thoracic Surgery (AREA)
• Medical Informatics (AREA)
• Molecular Biology (AREA)
• Surgery (AREA)
• Animal Behavior & Ethology (AREA)
• Public Health (AREA)
• Veterinary Medicine (AREA)
• Magnetic Resonance Imaging Apparatus (AREA)

Priority Applications (4)

Applications Claiming Priority (4)

Publications (1)

Family

ID=52358665

Family Applications (1)

Country Status (5)

Cited By (5)

Families Citing this family (18)

Citations (1)

Family Cites Families (9)

• 2016
  • 2016-01-21 US US15/544,397 patent/US10330757B2/en active Active
  • 2016-01-21 EP EP16701150.1A patent/EP3248023A1/en not_active Withdrawn
  • 2016-01-21 CN CN201680017078.1A patent/CN107407714B/zh not_active Expired - Fee Related
  • 2016-01-21 JP JP2017538307A patent/JP6640859B2/ja not_active Expired - Fee Related
  • 2016-01-21 WO PCT/EP2016/051199 patent/WO2016116545A1/en not_active Ceased

Patent Citations (1)

Non-Patent Citations (11)

Also Published As

Similar Documents

Legal Events