WO2009128013A1 - Real-time local and global sar estimation for patient safety and improved scanning performance - Google Patents

Real-time local and global sar estimation for patient safety and improved scanning performance Download PDF

Info

Publication number
WO2009128013A1
WO2009128013A1 PCT/IB2009/051531 IB2009051531W WO2009128013A1 WO 2009128013 A1 WO2009128013 A1 WO 2009128013A1 IB 2009051531 W IB2009051531 W IB 2009051531W WO 2009128013 A1 WO2009128013 A1 WO 2009128013A1
Authority
WO
WIPO (PCT)
Prior art keywords
sar
absorption rate
energy absorption
specific energy
set forth
Prior art date
Application number
PCT/IB2009/051531
Other languages
French (fr)
Inventor
Ingmar Graesslin
Sven Biederer
Ulrich Katscher
Ferdinand Schweser
Peter Boernert
Paul Royston Harvey
Wilhelmus Reinerius Maria Mens
Original Assignee
Koninklijke Philips Electronics N.V.
Philips Intellectual Property & Standards Gmbh
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips Electronics N.V., Philips Intellectual Property & Standards Gmbh filed Critical Koninklijke Philips Electronics N.V.
Priority to US12/936,540 priority Critical patent/US8941380B2/en
Priority to CN200980113091.7A priority patent/CN102007423B/en
Priority to EP09732227A priority patent/EP2269085A1/en
Priority to JP2011504587A priority patent/JP5731373B2/en
Publication of WO2009128013A1 publication Critical patent/WO2009128013A1/en

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/58Calibration of imaging systems, e.g. using test probes, Phantoms; Calibration objects or fiducial markers such as active or passive RF coils surrounding an MR active material
    • G01R33/583Calibration of signal excitation or detection systems, e.g. for optimal RF excitation power or frequency
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5611Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
    • G01R33/5612Parallel RF transmission, i.e. RF pulse transmission using a plurality of independent transmission channels
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/288Provisions 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/5659Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the RF magnetic field, e.g. spatial inhomogeneities of the RF magnetic field

Definitions

  • the present application relates to the diagnostic imaging arts. It finds particular application in the context of patient safety and associated improved scanning performance (in terms of radio frequency (RF) duty cycle) and will be described with particular reference thereto. Furthermore, the estimation and suppression of local specific energy absorption rate (SAR) hot spots in connection with high field magnetic resonance imaging (MRI). It is to be appreciated, however, that it is also applicable to optimization and processing of other information, and is not necessarily limited to the aforementioned applications.
  • RF radio frequency
  • the local SAR is a limiting factor.
  • the SAR deposition increases with higher field strength and limits the RF power, duty cycle, and flip angles usable, leading to a lengthening of scan acquisition time to meet designated SAR limits.
  • the SAR was relatively easy to calculate, as all antenna elements transmitted with the same amplitude, and a fixed phase shift between them.
  • the RF pulse shapes required for the experiments are known and are stored in a shape library, associated with their SAR.
  • RF safety is a prerequisite for in vivo parallel transmission MRI scans, that is, scanning within SAR limits using multi-channel RF transmit coils must be guaranteed. Scans cannot be started unless they are "SAR safe.”
  • SAR reduced RF pulses can be calculated by incorporating electrical field information into an RF pulse design.
  • methods have been used that construct RF pulses in consideration of known SAR hotspots that are common to every individual (e.g., the eyes). This is generally not sufficient for whole body imaging, as SAR hotspots can vary in both position and magnitude from patient to patient, and from RF pulse to RF pulse.
  • an RF pulse sequence that limits SAR to acceptable levels in one patient may not be so limiting with respect to another patient.
  • RF sequences that accommodate a known static hotspot may inadvertently exacerbate unknown, patient specific hotspots at other locations.
  • RF sequences are currently not constructed on a patient-by-patient basis, is that for clinically relevant spatially RF pulses (e.g. local excitation or zoom imaging), parallel transmission systems that can accelerate these kinds of RF pulses (TxSENSE) are required.
  • TxSENSE RF pulses
  • an underlying prerequisite is the availability to efficiently estimate the SAR.
  • availability of patient related E- fields and patient position are highly desirable for an accurate SAR estimation (including global and local SAR values and optionally a SAR map.)
  • Field data obtained by simulations differ to some degree from the actual fields in the scanner. Use of bio-mesh models for E-field simulations instead of the actual patient leads to a systematic error that is difficult to characterize.
  • the RF waveforms are identical for every Tx coil element and only a phase increment (e.g. 45° for 8 elements) exists.
  • a phase increment e.g. 45° for 8 elements.
  • the calculation is more complex, as each channel may have a different but static amplitude and phase.
  • more complex scans such as for 2D/3D spatially selective pulses, each channel may have dynamically changing amplitudes and phases.
  • the system performs a high number of calculations (e.g. TeraFLOPs): as high as 10 10 calculations or more, depending on the resolution of the model and cells used for the calculation. This process would take several minutes and cannot practically be carried out in real time as the patient waits inside the scanner for the actual diagnostic scan to begin.
  • the present application provides a new and improved magnetic resonance system, which overcomes the above -referenced problems and others.
  • a magnetic resonance system In accordance with one aspect, a magnetic resonance system is provided.
  • a main magnet generates a substantially uniform main magnetic field in an examination region.
  • a radio frequency assembly induces magnetic resonance in selected dipoles of a subject in the examination region, and receives the magnetic resonance.
  • a specific energy absorption rate calculation processor calculates a specific energy absorption rate and determines local specific energy absorption rate hotspots.
  • a sequence controller designs an RF excitation pulse that accounts for the local specific energy absorption rate hotspots and keeps energy delivered to the hotspots under acceptable levels.
  • a magnetic resonance system In accordance with another aspect, a magnetic resonance system is provided.
  • a main magnet generates a substantially uniform main magnetic field in an examination region.
  • a radio frequency assembly induces magnetic resonance in selected dipoles of a subject in the examination region, and receives the magnetic resonance.
  • a specific energy absorption rate calculation processor calculates a specific energy absorption rate and determines local specific energy absorption rate hotspots.
  • a graphics card processes non- graphics information in parallel.
  • a method of magnetic resonance is provided.
  • a substantially uniform main magnetic field is generated in an examination region.
  • Magnetic resonance is induced in selected dipoles of a subject in the examination region, and the magnetic resonance is received.
  • a position of the subject within the examination region is determined.
  • a specific energy absorption rate is calculated.
  • a globally safe RF pulse waveform that accounts for the calculated specific energy absorption rate is calculated. This may be done iteratively, if e.g., the first RF pulse estimate does not meet the SAR limits, or other system parameters, like the T R (repetition time) are prolonged.
  • One advantage is the ability to efficiently verify that a parallel transmission scan does not violate existing FDA or International Electro technical Commission limits.
  • Another advantage lies in increased calculation speed of SAR values, hotspots, and spatial SAR distribution. Another advantage lies in the ability to customize the SAR calculation to individual patients. Another advantage lies in the ability to create an optimal RF pulse sequence based on a patient's SAR profile.
  • Another advantage lies in the ability to determine specific information, such as E-fields and patient positions, on a patient-by-patient basis. Another advantage lies in the ability to detect anomalies or surgical implants for SAR model adaptation.
  • the invention may take form in various components and arrangements of components, and in various steps and arrangements of steps.
  • the drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention.
  • FIGURE 1 is a diagrammatic illustration of a magnetic resonance imaging apparatus in accordance with the present application
  • FIGURE 2 includes exemplary waveforms for PUC sampling during RF excitation
  • FIGURE 3 is a graph of dependence of hot spot suppression on an empirically determined weighting factor
  • FIGURE 4 is a comparison of alternate methods that consider less than all of the available information with an embodiment that considers all of the information;
  • FIGURE 5 is a comparison of a worst-case scenario method of calculating SAR to the embodiment that considers all of the information.
  • a magnetic resonance scanner 10 is illustrated as a closed bore system that includes a solenoidal main magnet assembly 12, although open and other magnet configurations are also contemplated.
  • the main magnet assembly 12 produces a substantially constant main magnetic field B 0 oriented along a horizontal axis of an imaging region. It is to be understood that other magnet arrangements, such as vertical, and other configurations are also contemplated.
  • the main magnet 12 in a bore type system may typically have a field strength of around 0.5 T to 7.0 T or more.
  • a gradient coil assembly 14 produces magnetic field gradients in the imaging region for spatially encoding the main magnetic field.
  • the magnetic field gradient coil assembly 14 includes coil segments configured to produce magnetic field gradient pulses in three orthogonal directions, typically longitudinal or z, transverse or x, and vertical or y directions.
  • a radio frequency coil assembly 16 including n coil elements 16i, 16 2 , ... 16 n , generates radio frequency pulses for exciting resonance in dipoles of the subject.
  • the signals that the radio frequency coil assembly 16 transmits are commonly known as the B 1 field.
  • the radio frequency coil assembly 16 also serves to detect resonance signals emanating from the imaging region.
  • the illustrated radio frequency coil assembly 16 is a send/receive coil that images the entire imaging region, however, local send/receive coils, local dedicated receive coils, or dedicated transmit coils are also contemplated.
  • the radio frequency coil assembly 16 includes an 8 channel transmit/receive antenna.
  • Gradient pulse amplifiers 18 deliver controlled electrical currents to the magnetic field gradient assembly 14 to produce selected magnetic field gradients.
  • the number of coil elements and the number of transmitters are the same. However, more than one coil element can be associated with each transmit channel.
  • a radio frequency receiver array 22, including n receivers 22i, 22 2 , ... 22 n in the illustrated embodiment, is coupled to the coil assembly 16 or a separate receive coil array to receive and demodulate the induced resonance signals.
  • a sequence controller 24 communicates with the gradient amplifiers 18 and the radio frequency transmitters 20i, 2O 2 , ... 2O n to excite and manipulate magnetic resonance in the region of interest.
  • the sequence controller 24 for example, produces selected repeated echo steady-state, or other resonance sequences, spatially encodes such resonances, selectively manipulates or spoils resonances, or otherwise generates selected magnetic resonance signals characteristic of the subject.
  • the generated resonance signals are detected by the RF coil assembly 16 or local coil assembly (not shown), communicated to the radio frequency receiver 22, demodulated, and stored in a k- space memory 26.
  • the imaging data is reconstructed by a reconstruction processor 28 to produce one or more image representations that are stored in an image memory 30.
  • the reconstruction processor 28 performs an inverse Fourier transform reconstruction.
  • the resultant image representation(s) is processed by a video processor 32 and displayed on a user interface 34 equipped with a human readable display.
  • the interface 34 is preferably a personal computer or workstation. Rather than producing a video image, the image representation can be processed by a printer driver and printed, transmitted over a computer network or the Internet, or the like.
  • the user interface 34 also allows a technician or other operator to communicate with the sequence controller 24 to select magnetic resonance imaging sequences, modify imaging sequences, execute imaging sequences, and so forth.
  • the user can select a SAR model, and all or part of the remaining parameters can be determined with user interaction and feedback.
  • a specific energy absorption rate (SAR) processor 36 calculates SAR for portions of the subject within the coil assembly 16.
  • the SAR calculation processor 36 creates a SAR map of the whole body that includes regions of increased SAR or hotspots.
  • SAR can be calculated very quickly as only the phase/amplitude relation is relevant.
  • the calculation of a single RF sample is sufficient as the amplitude/phase relation does not change the pulse.
  • a standard scan could be implemented by typical MRI systems and parallel transmission systems using constant phases and amplitudes.
  • the SAR can be expressed in a quadratic form in the pulse samples b ⁇ Qb, where ⁇ denotes the conjugate transpose, b is the RF waveform sample, and Q is a Hermitian positive definite matrix resulting from the solution of Maxwell's equations and corresponding to a specific subject volume.
  • the SAR map is created by considering several inputs, including trajectory, B 1 field maps, target excitation pattern, and a global Q matrix.
  • Existing SAR optimal algorithms typically only constrain a specific known static local region, such as the eyes. As mentioned previously, this is inadequate for full body imaging because other hotspots that vary by subject may be present.
  • a memory 35 can store pre-calculated data of one or multiple patient positions to prevent recalculation of the Q-matrices when unnecessary. Additionally, the memory 35 can store SAR values so that SAR need not be calculated repetitively for identical pulses. A unique ID can be used to identify the pulses. For local SAR calculation, the SAR of each volume element of a patient model is averaged until the desired mass is reached. The SAR value for a volume element is indicative of the SAR along the edge of the volume element, and the data is interpolated to acquire the SAR value at the center of each voxel. Some of the information can be pre-calculated and stored in a look up table
  • LUT Scanner specific information such as electric fields and B 1 field maps are stored in the LUT 37.
  • An appropriate starting bio-mesh can be selected from a body model memory of the LUT 37 by knowing the patient's height, weight, sex, and position in the MR scanner. While the height, sex, and weight can be entered by an operator before the scan, the patient position is determined by a subject position processor 39.
  • One way to obtain patient position and refine the patient model is to use a moving bed approach. Images are acquired while the patient is being moved into the bore of the scanner, resulting in a low-resolution 3D volume data set. Alternately, a short pre-scan could be performed once the patient is in the final position in the bore.
  • This data can be segmented, for example, by thresholding or some other processing means.
  • the position of the patient can be obtained by correlation methods with existing models, for example, from transverse slices of the patient or the detection of landmarks. Also, the patient volume and size could be estimated. Anomalies such as implants or missing organs can also be detected. In this way, the starting body model is customized to the current patient.
  • the 36 consults the LUT 37 to get the corresponding E-field data as a function of the input parameters (weight, sex, position, etc.) using an appropriate combination thereof. If deviation from all models stored in the LUT 37 is too large, then a very conservative SAR estimation could be used. In the case of implants for a particular scan, appropriate SAR limits for the device can be retrieved from the LUT 37. After the initial position has been determined, any table movement can be monitored by the subject position processor 39 and can be used to accurately determine the new position of the patient.
  • a coarsely segmented body model can be obtained from the moving bed imaging data or pre-scan, which can be used for a fast adaptation of existing E-fields of similar patients or a fast estimation of e.g. a homogeneous model.
  • Use of a homogenous model introduces inaccuracies in conductivity and permittivity that are relatively small. The differences between using data from a homogenous model and actual data are tolerable, such that use of a model is a viable alternative.
  • patient position can be determined by using pickup coils (PUCs).
  • PUCs pickup coils
  • Each transmit element of a multi channel transmit coil is equipped with a PUC for monitoring current in each element to ensure patient safety and facilitate system adjustment.
  • the patient's presence influences the coil's properties.
  • the loading of the RF coil elements changes during the movement of the patient through the magnet bore. This movement can be detected as a phase change. This can be translated to an approximate position of the patient in the MR system. This is possible because the currents in the coil elements are sampled during the RF pulses, as shown in FIGURE 2.
  • Exemplary RF and gradient waveforms are provided. The dotted lines represent the RF excitation waveforms, the dashed lines represent the MR signal sampling waveforms, and the solid lines represent the PUC sampling waveform. Additionally, the PUCs can be used to sense abnormal currents in the RF coil channels and initiate a scan termination if safety parameters are exceeded.
  • the sequence controller 24 designs an RF pulse sequence that is tailored to the present subject's SAR map. This could also be done by a host reconstructor or a separate graphics card that calculates RF pulses.
  • the sequence controller 24 iteratively processes the SAR map to find the best weighting factors q ! to satisfy existing SAR limits and to decrease the most limiting SAR value. More specifically, the sequence controller 24 directs the gradient assembly 14 and the RF assembly 16 to apply the newly designed RF pulse sequence. The SAR calculation processor 36 then recalculates the SAR map. Positions of local hotspots are once again determined.
  • the sequence controller 24 volume averages the Q-matrices (Qcntica L region(i)) of the hotspots, and weights them.
  • the weighting factors have been determined empirically based on the distance (z) of the hotspot from the isocenter of the magnet. With reference now to FIGURE 3, the troughs of the curves represent the optimal weighting factor for that distance.
  • the weighted, volume averaged Q-matrices are added to the global Q-matrix.
  • the radius of the spatial averaging around each hotspot, the hotspot positions, the local Q- matrices, and the selected weighting factors are all taken into consideration when recalculating Q 1 , the updated global Q-matrix.
  • the sequence controller 24 designs a new, SAR optimized RF pulse sequence based on the updated Q-matrix Q 1 . As before, lower SAR values are obtained at the critical regions.
  • the sequence controller 24 and the SAR calculation processor 36 can apply one or more of the above steps iteratively until SAR converges to a minimum value at the hotspots, or alternatively, until a desired safe SAR level is reached. In some cases, it might not be necessary to apply the iterations until SAR converges to a minimum if SAR reaches a safe level before it converges.
  • the repetition time T R can also be prolonged, or the flip angle can be reduced, or a combination of the two.
  • the RF pulse can be re- optimized if the patient is moved.
  • each Q-matrix calculation accounts for the correct amplitude information and the correct phase information for each channel involved.
  • each voxel of the body is calculated separately, giving the highest resolution possible.
  • the average amount of voxels in a bio-mesh is on the order of 750,000 for a voxel size of 5 mm.
  • TeraFLOPs are required to calculate the global and local SAR and to produce a SAR map.
  • one embodiment includes an RF assembly 16 with eight channels, but it is to be understood that assemblies with more channels are possible, with any arbitrary combination of channels operating at any given time. SAR is calculated for these situations accordingly.
  • the SAR calculation processor 36 delegates the task to a sub-processor 38, such as a high-performance graphics card.
  • the sub-processor 38 can be located in the SAR calculation processor 36 itself, in a host computer, or in a spectrometer. Since the SAR calculations of the individual voxels are not dependant on one another, they do not have to be processed one after another, that is, they can be processed in parallel.
  • the sub-processor 38 such as a graphics card offers many parallel processing channels (e.g. 128, 256 etc.) to speed up the calculation of the SAR.
  • calculation of the SAR was accelerated by a factor of 100 over using a 3 GHz processor alone to calculate the SAR. Resultantly, calculation of the SAR of an RF pulse for a single biomesh can be performed in seconds instead of minutes. This allows the iterative process of converging SAR hotspots to minima described above to be performed in a practically applicable amount of time.
  • the SAR calculation processor 36 can complete the calculation so that scanning is still possible.
  • voxels can be grouped by their proximity and averaged, reducing the number of volume elements from roughly 750,000 to, for example, 100,000.
  • amplitudes are set to maximum in each corresponding channel. This method again cuts the amount of calculations down, since only the maximum amplitude for each channel is considered, but sacrifices the quality of the resultant calculations, again, erring on the side of caution.
  • FIGURES 4 and 5 illustrate some of the alternate embodiments that consider less than all the information compared to the embodiment that considers all of the available information.
  • curve 50 represents the embodiment that considers correct amplitudes, but no phase.
  • Curve 52 represents the embodiment that considers maximum amplitudes in the correct channels.
  • Curve 54 represents the worst case scenario embodiment, where not even channel information is considered.
  • the curves 50, 52, 54 depict the ratio of calculated to actual SAR as a function of the reduction factor. As can be seen, as more information is considered, the closer the estimation of SAR comes to the actual SAR.
  • FIGURE 5 illustrates calculations made using the worst case scenario method 58 compared against the actual calculations 60.
  • FIGURE 5 is the position dependence of the local trunk SAR of an 8-channel body coil emulating a standard one-channel body coil, with a phase of 45° and an amplitude of one on all channels. It is evident that the worst case scenario method greatly overestimates SAR, especially in the midsection of the patient, leading to less accurate calculations of SAR. Also evident is the dependence of SAR on position.
  • the SAR calculation processor 36, sub- processor 38 or any other components can be located on a remote server. Multiple clients can be served simultaneously by the server. When multiple requests for SAR values appear concurrently, the server can prioritize them based on the order of arrival, or based on other priorities.

Landscapes

  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

In a method and apparatus to enable increased RF duty cycle in high field MR scans, a specific energy absorption rate (SAR) calculation processor (36) calculates the local and global SAR or even a spatial SAR map. The efficient implementation by using pre- averaged data (based on E-fields) makes a fast and accurate SAR estimation possible. By incorporating additional information as e.g. patient position the SAR calculation accuracy can be increased as well as by using more patient specific precalculated information (e.g. based on different bio meshes), the so called Q-matrices. Optionally, a sequence controller (24) creates a global SAR optimal RF pulse. After the optimal RF pulse is applied, the SAR and its spatial distribution are determined. SAR hotspots are also determined. Q-matrices within an appropriate radius around the hotspots are averaged and added to a global Q-matrix in a weighted fashion. After the global Q-matrix is updated, a new optimal RF pulse is created. One or more of the steps can be repeated to converge the SAR to a minimum.

Description

Real-Time Local and Global SAR Estimation for Patient Safety and Improved Scanning Performance
FIELD OF THE INVENTION
The present application relates to the diagnostic imaging arts. It finds particular application in the context of patient safety and associated improved scanning performance (in terms of radio frequency (RF) duty cycle) and will be described with particular reference thereto. Furthermore, the estimation and suppression of local specific energy absorption rate (SAR) hot spots in connection with high field magnetic resonance imaging (MRI). It is to be appreciated, however, that it is also applicable to optimization and processing of other information, and is not necessarily limited to the aforementioned applications.
BACKGROUND OF THE INVENTION For many MR applications at higher field strengths, the local SAR is a limiting factor. The SAR deposition increases with higher field strength and limits the RF power, duty cycle, and flip angles usable, leading to a lengthening of scan acquisition time to meet designated SAR limits. For a single transmitter system, the SAR was relatively easy to calculate, as all antenna elements transmitted with the same amplitude, and a fixed phase shift between them. Furthermore, the RF pulse shapes required for the experiments are known and are stored in a shape library, associated with their SAR. With the advent of multi-transmit systems where each coil element has the potential to independently transmit its own amplitude and phase, SAR must be calculated on a per channel basis also considering the parallel RF transmission pulses, which can only be calculated based on additional information, for example, B1 maps, and are thus experiment/patient specific.
RF safety is a prerequisite for in vivo parallel transmission MRI scans, that is, scanning within SAR limits using multi-channel RF transmit coils must be guaranteed. Scans cannot be started unless they are "SAR safe." In an MR system with multiple transmit channels, SAR reduced RF pulses can be calculated by incorporating electrical field information into an RF pulse design. In the past, methods have been used that construct RF pulses in consideration of known SAR hotspots that are common to every individual (e.g., the eyes). This is generally not sufficient for whole body imaging, as SAR hotspots can vary in both position and magnitude from patient to patient, and from RF pulse to RF pulse. Thus, an RF pulse sequence that limits SAR to acceptable levels in one patient may not be so limiting with respect to another patient. Moreover, RF sequences that accommodate a known static hotspot may inadvertently exacerbate unknown, patient specific hotspots at other locations.
One possible solution is to develop a worst-case scenario estimation of SAR that would be safe for all patients. This solution, however, would limit the allowed RF duty cycle so much that the MRI system would become seriously compromised for use in conjunction with in vivo parallel transmission scans. The ability to tailor a SAR calculation to the patient would be more beneficial than using a blanket scenario or known term for all patients.
One reason in particular that RF sequences are currently not constructed on a patient-by-patient basis, is that for clinically relevant spatially RF pulses (e.g. local excitation or zoom imaging), parallel transmission systems that can accelerate these kinds of RF pulses (TxSENSE) are required. For an RF sequence, an underlying prerequisite is the availability to efficiently estimate the SAR. Moreover, availability of patient related E- fields and patient position are highly desirable for an accurate SAR estimation (including global and local SAR values and optionally a SAR map.) Field data obtained by simulations differ to some degree from the actual fields in the scanner. Use of bio-mesh models for E-field simulations instead of the actual patient leads to a systematic error that is difficult to characterize. For a standard, single channel birdcage coil RF transmit assembly, the RF waveforms are identical for every Tx coil element and only a phase increment (e.g. 45° for 8 elements) exists. For multiple Tx coil elements, the calculation is more complex, as each channel may have a different but static amplitude and phase. In more complex scans, such as for 2D/3D spatially selective pulses, each channel may have dynamically changing amplitudes and phases.
In order to calculate all SAR types as specified in the standard (local and global) and optionally a SAR map of a patient for a multi-channel RF transmit system, (e.g., eight transmit channels) the system performs a high number of calculations (e.g. TeraFLOPs): as high as 1010 calculations or more, depending on the resolution of the model and cells used for the calculation. This process would take several minutes and cannot practically be carried out in real time as the patient waits inside the scanner for the actual diagnostic scan to begin. The present application provides a new and improved magnetic resonance system, which overcomes the above -referenced problems and others.
SUMMARY OF THE INVENTION
In accordance with one aspect, a magnetic resonance system is provided. A main magnet generates a substantially uniform main magnetic field in an examination region. A radio frequency assembly induces magnetic resonance in selected dipoles of a subject in the examination region, and receives the magnetic resonance. A specific energy absorption rate calculation processor calculates a specific energy absorption rate and determines local specific energy absorption rate hotspots. A sequence controller designs an RF excitation pulse that accounts for the local specific energy absorption rate hotspots and keeps energy delivered to the hotspots under acceptable levels.
In accordance with another aspect, a magnetic resonance system is provided. A main magnet generates a substantially uniform main magnetic field in an examination region. A radio frequency assembly induces magnetic resonance in selected dipoles of a subject in the examination region, and receives the magnetic resonance. A specific energy absorption rate calculation processor calculates a specific energy absorption rate and determines local specific energy absorption rate hotspots. A graphics card processes non- graphics information in parallel.
In accordance with another aspect, a method of magnetic resonance is provided. A substantially uniform main magnetic field is generated in an examination region. Magnetic resonance is induced in selected dipoles of a subject in the examination region, and the magnetic resonance is received. A position of the subject within the examination region is determined. A specific energy absorption rate is calculated. A globally safe RF pulse waveform that accounts for the calculated specific energy absorption rate is calculated. This may be done iteratively, if e.g., the first RF pulse estimate does not meet the SAR limits, or other system parameters, like the TR (repetition time) are prolonged.
One advantage is the ability to efficiently verify that a parallel transmission scan does not violate existing FDA or International Electro technical Commission limits.
Another advantage lies in increased calculation speed of SAR values, hotspots, and spatial SAR distribution. Another advantage lies in the ability to customize the SAR calculation to individual patients. Another advantage lies in the ability to create an optimal RF pulse sequence based on a patient's SAR profile.
Another advantage lies in the ability to determine specific information, such as E-fields and patient positions, on a patient-by-patient basis. Another advantage lies in the ability to detect anomalies or surgical implants for SAR model adaptation.
Still further advantages of the present invention will be appreciated to those of ordinary skill in the art upon reading and understand the following detailed description.
BRIEF DESCRIPTION OF THE DRAWINGS
The invention may take form in various components and arrangements of components, and in various steps and arrangements of steps. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention.
FIGURE 1 is a diagrammatic illustration of a magnetic resonance imaging apparatus in accordance with the present application; FIGURE 2 includes exemplary waveforms for PUC sampling during RF excitation;
FIGURE 3 is a graph of dependence of hot spot suppression on an empirically determined weighting factor;
FIGURE 4 is a comparison of alternate methods that consider less than all of the available information with an embodiment that considers all of the information;
FIGURE 5 is a comparison of a worst-case scenario method of calculating SAR to the embodiment that considers all of the information.
DETAILED DESCRIPTION OF EMBODIMENTS
With reference to FIGURE 1, a magnetic resonance scanner 10 is illustrated as a closed bore system that includes a solenoidal main magnet assembly 12, although open and other magnet configurations are also contemplated. The main magnet assembly 12 produces a substantially constant main magnetic field B0 oriented along a horizontal axis of an imaging region. It is to be understood that other magnet arrangements, such as vertical, and other configurations are also contemplated. The main magnet 12 in a bore type system may typically have a field strength of around 0.5 T to 7.0 T or more. A gradient coil assembly 14 produces magnetic field gradients in the imaging region for spatially encoding the main magnetic field. Preferably, the magnetic field gradient coil assembly 14 includes coil segments configured to produce magnetic field gradient pulses in three orthogonal directions, typically longitudinal or z, transverse or x, and vertical or y directions.
A radio frequency coil assembly 16, including n coil elements 16i, 162, ... 16n, generates radio frequency pulses for exciting resonance in dipoles of the subject. The signals that the radio frequency coil assembly 16 transmits are commonly known as the B1 field. The radio frequency coil assembly 16 also serves to detect resonance signals emanating from the imaging region. The illustrated radio frequency coil assembly 16 is a send/receive coil that images the entire imaging region, however, local send/receive coils, local dedicated receive coils, or dedicated transmit coils are also contemplated. In one embodiment, the radio frequency coil assembly 16 includes an 8 channel transmit/receive antenna. Gradient pulse amplifiers 18 deliver controlled electrical currents to the magnetic field gradient assembly 14 to produce selected magnetic field gradients. A radio frequency transmitter array 20, including n transmitters 20i, 2O2, ... 2On, preferably digital, applies radio frequency pulses or pulse packets to the radio frequency coil assembly 16 to excite selected resonance. In the illustrated embodiment, the number of coil elements and the number of transmitters are the same. However, more than one coil element can be associated with each transmit channel. A radio frequency receiver array 22, including n receivers 22i, 222, ... 22n in the illustrated embodiment, is coupled to the coil assembly 16 or a separate receive coil array to receive and demodulate the induced resonance signals.
To acquire resonance imaging data of a subject, the subject is placed inside the imaging region. A sequence controller 24 communicates with the gradient amplifiers 18 and the radio frequency transmitters 20i, 2O2, ... 2On to excite and manipulate magnetic resonance in the region of interest. The sequence controller 24, for example, produces selected repeated echo steady-state, or other resonance sequences, spatially encodes such resonances, selectively manipulates or spoils resonances, or otherwise generates selected magnetic resonance signals characteristic of the subject. The generated resonance signals are detected by the RF coil assembly 16 or local coil assembly (not shown), communicated to the radio frequency receiver 22, demodulated, and stored in a k- space memory 26. The imaging data is reconstructed by a reconstruction processor 28 to produce one or more image representations that are stored in an image memory 30. In one suitable embodiment, the reconstruction processor 28 performs an inverse Fourier transform reconstruction.
The resultant image representation(s) is processed by a video processor 32 and displayed on a user interface 34 equipped with a human readable display. The interface 34 is preferably a personal computer or workstation. Rather than producing a video image, the image representation can be processed by a printer driver and printed, transmitted over a computer network or the Internet, or the like. Preferably, the user interface 34 also allows a technician or other operator to communicate with the sequence controller 24 to select magnetic resonance imaging sequences, modify imaging sequences, execute imaging sequences, and so forth. At the interface 34, the user can select a SAR model, and all or part of the remaining parameters can be determined with user interaction and feedback.
A specific energy absorption rate (SAR) processor 36 calculates SAR for portions of the subject within the coil assembly 16. In one embodiment, the SAR calculation processor 36 creates a SAR map of the whole body that includes regions of increased SAR or hotspots. For a standard scan in which there is a constant change of amplitude and of phase for standard RF pulses, SAR can be calculated very quickly as only the phase/amplitude relation is relevant. Thus, the calculation of a single RF sample is sufficient as the amplitude/phase relation does not change the pulse. A standard scan could be implemented by typical MRI systems and parallel transmission systems using constant phases and amplitudes.
Provided that the RF field inside the subject responds linearly to the currents driving the field, the SAR can be expressed in a quadratic form in the pulse samples b^Qb, where ^ denotes the conjugate transpose, b is the RF waveform sample, and Q is a Hermitian positive definite matrix resulting from the solution of Maxwell's equations and corresponding to a specific subject volume. The SAR map is created by considering several inputs, including trajectory, B1 field maps, target excitation pattern, and a global Q matrix. Existing SAR optimal algorithms typically only constrain a specific known static local region, such as the eyes. As mentioned previously, this is inadequate for full body imaging because other hotspots that vary by subject may be present. Statistically constraining a spatial region in which a hotspot occurs may result in new hotspots at other locations. A memory 35 can store pre-calculated data of one or multiple patient positions to prevent recalculation of the Q-matrices when unnecessary. Additionally, the memory 35 can store SAR values so that SAR need not be calculated repetitively for identical pulses. A unique ID can be used to identify the pulses. For local SAR calculation, the SAR of each volume element of a patient model is averaged until the desired mass is reached. The SAR value for a volume element is indicative of the SAR along the edge of the volume element, and the data is interpolated to acquire the SAR value at the center of each voxel. Some of the information can be pre-calculated and stored in a look up table
(LUT) 37. Scanner specific information such as electric fields and B1 field maps are stored in the LUT 37. An appropriate starting bio-mesh can be selected from a body model memory of the LUT 37 by knowing the patient's height, weight, sex, and position in the MR scanner. While the height, sex, and weight can be entered by an operator before the scan, the patient position is determined by a subject position processor 39. One way to obtain patient position and refine the patient model is to use a moving bed approach. Images are acquired while the patient is being moved into the bore of the scanner, resulting in a low-resolution 3D volume data set. Alternately, a short pre-scan could be performed once the patient is in the final position in the bore. This data can be segmented, for example, by thresholding or some other processing means. Next, the position of the patient can be obtained by correlation methods with existing models, for example, from transverse slices of the patient or the detection of landmarks. Also, the patient volume and size could be estimated. Anomalies such as implants or missing organs can also be detected. In this way, the starting body model is customized to the current patient. Once the patient position has been determined, the SAR calculation processor
36 consults the LUT 37 to get the corresponding E-field data as a function of the input parameters (weight, sex, position, etc.) using an appropriate combination thereof. If deviation from all models stored in the LUT 37 is too large, then a very conservative SAR estimation could be used. In the case of implants for a particular scan, appropriate SAR limits for the device can be retrieved from the LUT 37. After the initial position has been determined, any table movement can be monitored by the subject position processor 39 and can be used to accurately determine the new position of the patient.
In an alternate embodiment, a coarsely segmented body model can be obtained from the moving bed imaging data or pre-scan, which can be used for a fast adaptation of existing E-fields of similar patients or a fast estimation of e.g. a homogeneous model. Use of a homogenous model introduces inaccuracies in conductivity and permittivity that are relatively small. The differences between using data from a homogenous model and actual data are tolerable, such that use of a model is a viable alternative. In another alternate embodiment, patient position can be determined by using pickup coils (PUCs). Each transmit element of a multi channel transmit coil is equipped with a PUC for monitoring current in each element to ensure patient safety and facilitate system adjustment. In general, the patient's presence influences the coil's properties. Thus, the loading of the RF coil elements changes during the movement of the patient through the magnet bore. This movement can be detected as a phase change. This can be translated to an approximate position of the patient in the MR system. This is possible because the currents in the coil elements are sampled during the RF pulses, as shown in FIGURE 2. Exemplary RF and gradient waveforms are provided. The dotted lines represent the RF excitation waveforms, the dashed lines represent the MR signal sampling waveforms, and the solid lines represent the PUC sampling waveform. Additionally, the PUCs can be used to sense abnormal currents in the RF coil channels and initiate a scan termination if safety parameters are exceeded.
Once the SAR map is created, the sequence controller 24 designs an RF pulse sequence that is tailored to the present subject's SAR map. This could also be done by a host reconstructor or a separate graphics card that calculates RF pulses. The sequence controller 24 introduces weighting factors that specify a trade-off between different hotspot regions and the global SAR. For instance, depending of the spatial SAR distribution of an RF pulse that is optimal with respect to global SAR, hotspot reduction is possible via Q1 = Qgiobai + ∑ qi QcπticaLregion(i), where Q1 is the modified Q-matrix, Qgiobai is the original global Q-matrix, q! is a weighting factor, and QCπticai_region(i) is a Q-matrix of a volume immediately around the hotspot (e.g. a 3x3x3 voxel volume). The sequence controller 24 iteratively processes the SAR map to find the best weighting factors q! to satisfy existing SAR limits and to decrease the most limiting SAR value. More specifically, the sequence controller 24 directs the gradient assembly 14 and the RF assembly 16 to apply the newly designed RF pulse sequence. The SAR calculation processor 36 then recalculates the SAR map. Positions of local hotspots are once again determined. Next, the sequence controller 24 volume averages the Q-matrices (QcnticaLregion(i)) of the hotspots, and weights them. The weighting factors have been determined empirically based on the distance (z) of the hotspot from the isocenter of the magnet. With reference now to FIGURE 3, the troughs of the curves represent the optimal weighting factor for that distance. Curve 40 represents hotspot suppression with z=20 cm. Curve 42 represents hotspot suppression with z=40 cm. Curve 44 represents hotspot suppression with z=60 cm. Curve 46 represents hotspot suppression with z=80 cm. Lastly, curve 48 represents hotspot suppression with z=100 cm.
The weighted, volume averaged Q-matrices are added to the global Q-matrix. The radius of the spatial averaging around each hotspot, the hotspot positions, the local Q- matrices, and the selected weighting factors are all taken into consideration when recalculating Q1, the updated global Q-matrix. Once the Q-matrix has been updated, the sequence controller 24 designs a new, SAR optimized RF pulse sequence based on the updated Q-matrix Q1. As before, lower SAR values are obtained at the critical regions. The sequence controller 24 and the SAR calculation processor 36 can apply one or more of the above steps iteratively until SAR converges to a minimum value at the hotspots, or alternatively, until a desired safe SAR level is reached. In some cases, it might not be necessary to apply the iterations until SAR converges to a minimum if SAR reaches a safe level before it converges. Alternatively, the repetition time TR can also be prolonged, or the flip angle can be reduced, or a combination of the two. Also, the RF pulse can be re- optimized if the patient is moved.
This iterative process is computationally intense, requiring a large amount of data processing capacity each time the global Q-matrix is updated. With existing systems, each iteration could take several minutes, which is impractical with a patient waiting in the scanner. Each Q-matrix calculation accounts for the correct amplitude information and the correct phase information for each channel involved. In one embodiment, each voxel of the body is calculated separately, giving the highest resolution possible. The average amount of voxels in a bio-mesh is on the order of 750,000 for a voxel size of 5 mm. When phase, and amplitude information is processed for each RF channel's effect on each voxel, a high number of calculations (e.g. TeraFLOPs) are required to calculate the global and local SAR and to produce a SAR map. As mentioned previously, one embodiment includes an RF assembly 16 with eight channels, but it is to be understood that assemblies with more channels are possible, with any arbitrary combination of channels operating at any given time. SAR is calculated for these situations accordingly.
In the embodiment of FIGURE 1, the SAR calculation processor 36 delegates the task to a sub-processor 38, such as a high-performance graphics card. The sub-processor 38 can be located in the SAR calculation processor 36 itself, in a host computer, or in a spectrometer. Since the SAR calculations of the individual voxels are not dependant on one another, they do not have to be processed one after another, that is, they can be processed in parallel. The sub-processor 38 such as a graphics card offers many parallel processing channels (e.g. 128, 256 etc.) to speed up the calculation of the SAR. For example, by using a graphics card with 128 processing channels, calculation of the SAR was accelerated by a factor of 100 over using a 3 GHz processor alone to calculate the SAR. Resultantly, calculation of the SAR of an RF pulse for a single biomesh can be performed in seconds instead of minutes. This allows the iterative process of converging SAR hotspots to minima described above to be performed in a practically applicable amount of time. In one embodiment, if the sub-processor 38 is unavailable (for example, if the graphics card is broken) then the SAR calculation processor 36 can complete the calculation so that scanning is still possible. In an alternate embodiment, voxels can be grouped by their proximity and averaged, reducing the number of volume elements from roughly 750,000 to, for example, 100,000. This further reduces calculation time of the SAR, but sacrifices some resolution and accuracy in the calculated SAR maps. As a consequence, an extra safety margin is added to obtain the estimated SAR values for a scan. In another alternate embodiment, amplitude and channel information is considered, but phase information is not. This also speeds up the calculations, but calculates a less accurate SAR map, erring on the side of caution. The SAR values are overestimated in this embodiment.
In another alternate embodiment, amplitudes are set to maximum in each corresponding channel. This method again cuts the amount of calculations down, since only the maximum amplitude for each channel is considered, but sacrifices the quality of the resultant calculations, again, erring on the side of caution.
In another alternate embodiment, a worst case scenario embodiment, only the maximum amplitude, regardless of the channel is considered. This results in only a coarse estimation of the actual SAR map. FIGURES 4 and 5 illustrate some of the alternate embodiments that consider less than all the information compared to the embodiment that considers all of the available information. In FIGURE 4, curve 50 represents the embodiment that considers correct amplitudes, but no phase. Curve 52 represents the embodiment that considers maximum amplitudes in the correct channels. Curve 54 represents the worst case scenario embodiment, where not even channel information is considered. The curves 50, 52, 54 depict the ratio of calculated to actual SAR as a function of the reduction factor. As can be seen, as more information is considered, the closer the estimation of SAR comes to the actual SAR. If these error ratios are acceptable, however, calculation time can be saved by using one of the alternate methods. FIGURE 5 illustrates calculations made using the worst case scenario method 58 compared against the actual calculations 60. FIGURE 5 is the position dependence of the local trunk SAR of an 8-channel body coil emulating a standard one-channel body coil, with a phase of 45° and an amplitude of one on all channels. It is evident that the worst case scenario method greatly overestimates SAR, especially in the midsection of the patient, leading to less accurate calculations of SAR. Also evident is the dependence of SAR on position.
In another alternate embodiment, the SAR calculation processor 36, sub- processor 38 or any other components can be located on a remote server. Multiple clients can be served simultaneously by the server. When multiple requests for SAR values appear concurrently, the server can prioritize them based on the order of arrival, or based on other priorities.
The invention has been described with reference to the preferred embodiments. Modifications and alterations may occur to others upon reading and understanding the preceding detailed description. It is intended that the invention be construed as including all such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof.

Claims

CLAIMS:Having thus described the preferred embodiments, the invention is now claimed to be:
1. A magnetic resonance system comprising: a main magnet (12) which generates a substantially uniform main magnetic field in an examination region; a multi-channel radio frequency coil assembly (16) including a plurality of coil elements (16i, 162, ... 16n) and a transmitter array (20) including a plurality of transmitters (20i, 2O2, ... 2On) which induces magnetic resonance in selected dipoles of a subject in the examination region; a receive array (22) which receives magnetic resonance signals; a specific energy absorption rate calculation processor (36) that calculates a specific energy absorption rate for the multi-channel radio frequency coil assembly (16); and a sequence controller (24) that reconfigures RF excitation pulses that account for the local specific energy absorption rate hotspots.
2. The magnetic resonance system as set forth in claim 1, further including a graphics processing unit (38) that includes a plurality of parallel processing channels that aids in processing the calculation of the specific energy absorption rate and energy absorption rate hotspots.
3. The magnetic resonance system as set forth in claim 1, further including: a memory (37) that includes pre-determined values for at least one of electric fields, Q-matrices, and model patient data.
4. The magnetic resonance system as set forth in claim 1, further including: a subject position calculation processor (39) that determines a position of the subject within the main magnet (12) and provides the position to the specific energy absorption rate calculation processor (36).
5. A method of magnetic resonance comprising: generating a substantially uniform main magnetic field in an examination region; inducing magnetic resonance in selected dipoles of a subject in the examination region, and receiving the magnetic resonance; determining a position of the subject within the examination region; calculating a specific energy absorption rate; and calculating a globally safe RF pulse waveform that accounts for the calculated specific energy absorption rate.
6. The method as set forth in 5, further including: determining a position of at least one local specific energy absorption rate hotspot in the subject.
7. The method as set forth in claim 6, further including: applying the globally safe RF pulse; and calculating a spatial specific energy absorption rate distribution induced by the globally safe RF pulse in the subject;
8. The method as set forth in claim 6, further including: updating a global Q-matrix by adding in a weighted volume-averaged Q- matrix of the at least one hotspot.
9. The method as set forth in claim 8, wherein the weighted volume-averaged Q- matrix is calculated by considering a radius of spatial averaging around the at least one hotspot, a position of the at least one hotspot, at least one local Q-matrix, and weighting factors.
10. The method as set forth in claim 8, further including: calculating an optimized RF pulse waveform from the updated global Q- matrix; applying the optimized RF pulse waveform.
11. The method as set forth in claim 10, further including: repeating at least the steps of updating the global Q-matrix and calculating an optimized RF pulse waveform until the specific energy absorption rate converges to a desired value.
12. The method as set forth in claim 5 wherein the step of calculating the globally safe RF pulse waveform includes accounting for a trajectory of the RF pulse, Bl field maps, a target excitation pattern, and a global Q-matrix.
13. The method as set forth in claim 5, further including: updating the globally safe RF pulse based on the position of the subject.
14. A computer readable medium storing a program which controls a magnetic resonance system to perform the method of claim 5.
15. The method as set forth in claim 5, wherein the step of calculating a specific energy absorption rate includes calculating the specific energy absorption rate based on at least one of a standard scan that includes RF pulses with constant amplitudes and phases and a spatially selective scan that uses 3D RF pulses for at least one of patient safety and hot spot suppression.
PCT/IB2009/051531 2008-04-16 2009-04-13 Real-time local and global sar estimation for patient safety and improved scanning performance WO2009128013A1 (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
US12/936,540 US8941380B2 (en) 2008-04-16 2009-04-13 Real-time local and global SAR estimation for patient safety and improved scanning performance
CN200980113091.7A CN102007423B (en) 2008-04-16 2009-04-13 Magnetic resonance system and method
EP09732227A EP2269085A1 (en) 2008-04-16 2009-04-13 Real-time local and global sar estimation for patient safety and improved scanning performance
JP2011504587A JP5731373B2 (en) 2008-04-16 2009-04-13 Real-time local and global SAR estimation for improved patient safety and scanning performance

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP08103561.0 2008-04-16
EP08103561 2008-04-16

Publications (1)

Publication Number Publication Date
WO2009128013A1 true WO2009128013A1 (en) 2009-10-22

Family

ID=40765597

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/IB2009/051531 WO2009128013A1 (en) 2008-04-16 2009-04-13 Real-time local and global sar estimation for patient safety and improved scanning performance

Country Status (5)

Country Link
US (1) US8941380B2 (en)
EP (1) EP2269085A1 (en)
JP (1) JP5731373B2 (en)
CN (1) CN102007423B (en)
WO (1) WO2009128013A1 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103154761A (en) * 2010-08-20 2013-06-12 皇家飞利浦电子股份有限公司 Virtual coil emulation in parallel transmission mri
DE102010015066B4 (en) * 2010-04-15 2014-12-11 Siemens Aktiengesellschaft Method and device for determining a magnetic resonance system drive sequence and method for operating a magnetic resonance system
CN105093144A (en) * 2014-05-22 2015-11-25 西门子公司 Method and apparatus for generating magnetic resonance image using combined excitation by whole body coil and local coil

Families Citing this family (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011515179A (en) * 2008-03-26 2011-05-19 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Determination of in vivo local SAR and conductivity mapping
DE102009030721B4 (en) * 2009-06-26 2013-04-04 Siemens Aktiengesellschaft SAR calculation for multi-channel MR transmission systems
EP2296000B1 (en) * 2009-09-10 2014-03-12 Commissariat à l'Énergie Atomique et aux Énergies Alternatives Parallel excitation of nuclear spins with local SAR control
JP5542591B2 (en) * 2009-11-12 2014-07-09 株式会社東芝 Magnetic resonance imaging apparatus and magnetic resonance imaging method
DE102010011588B4 (en) * 2010-03-12 2012-04-05 Friedrich-Alexander-Universität Erlangen-Nürnberg Method for determining sensitivity matrices for critical hotspots
DE102010015044A1 (en) * 2010-04-15 2011-10-20 Siemens Aktiengesellschaft Method and device for determining a magnetic resonance system drive sequence and method for operating a magnetic resonance system
BR112012029855A8 (en) * 2010-05-27 2016-09-13 Koninklijke Philips Electronics Nv MAGNETIC RESONANCE EXAMINATION SYSTEM FOR THE EXAMINATION OF AN OBJECT, METHOD AND COMPUTER PROGRAM
US8653818B2 (en) * 2011-04-08 2014-02-18 Siemens Aktiengesellschaft Parallel transmission RF pulse design with local SAR constraints
WO2013011406A2 (en) * 2011-07-20 2013-01-24 Koninklijke Philips Electronics N.V. Wireless local transmit coils and array with controllable load
EP2745133B1 (en) 2011-08-17 2022-03-09 Koninklijke Philips N.V. Reducing the radio-frequency transmit field in a predetermined volume during magnetic resonance imaging
CN104023627A (en) * 2012-08-29 2014-09-03 株式会社东芝 Magnetic resonance imaging device and method for computing specific absorption ratio thereof
CN103705237B (en) * 2012-09-28 2016-05-11 上海联影医疗科技有限公司 Guard method, the device of magnetic resonance imaging system and certainly turn-off method of testing
WO2014141109A1 (en) * 2013-03-13 2014-09-18 Koninklijke Philips N.V. Multi-element rf transmit coil for magnetic resonance imaging
JP6395388B2 (en) * 2013-06-26 2018-09-26 株式会社日立製作所 Magnetic resonance imaging apparatus and control method thereof
US20160299202A1 (en) * 2013-12-18 2016-10-13 Hitachi, Ltd. Magnetic resonance imaging apparatus and control method thereof
US10661091B2 (en) * 2014-03-26 2020-05-26 The General Hospital Corporation System and method for hyperthermia treatment using radiofrequency phased arrays
US11428761B2 (en) * 2014-06-04 2022-08-30 Koninklijke Philips N.V. Patient proximity-modulated specific absorption rate
DE102014211574A1 (en) * 2014-06-17 2015-12-17 Siemens Aktiengesellschaft Adjusting the table position in MR imaging
US20160091583A1 (en) * 2014-09-30 2016-03-31 Siemens Aktiengesellschaft Patient-Specific Estimation of Specific Absorption Rate
US9939501B2 (en) * 2014-11-12 2018-04-10 The Board Of Trustees Of The Leland Stanford Junior University Iterative minimization procedure with uncompressed local SAR estimate
DE112016001879T5 (en) 2015-04-24 2018-01-04 Koninklijke Philips N.V. A multi-channel transmit / receive radio frequency (RF) system
DE102016200611A1 (en) * 2016-01-19 2017-07-20 Siemens Healthcare Gmbh Storage evaluation of an examination object
CN105740534A (en) * 2016-01-28 2016-07-06 珠海和佳医疗设备股份有限公司 Calculation method of axisymmetric physical field
JP2017164211A (en) * 2016-03-15 2017-09-21 株式会社日立製作所 Magnetic resonance imaging device, q value calculation method and specific absorption rate control method
US10222472B2 (en) * 2016-09-30 2019-03-05 Veoneer Us, Inc. System and method for detecting heading and velocity of a target object
US10444305B2 (en) * 2017-06-21 2019-10-15 Canon Medical Systems Corporation Magnetic resonance imaging apparatus and specific absorption rate calculation method
CN107495967B (en) * 2017-08-24 2020-06-19 上海联影医疗科技有限公司 Method, device and system for predicting and controlling radio frequency energy deposition and storage medium
DE102018212091A1 (en) * 2018-07-19 2020-01-23 Siemens Healthcare Gmbh Determination of a load-optimized MR sequence
DE102021208408A1 (en) 2021-08-03 2023-02-09 Siemens Healthcare Gmbh Absorption monitoring method, magnetic resonance apparatus and computer program product

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005001502A1 (en) * 2003-06-30 2005-01-06 Koninklijke Philips Electronics N.V. Control of specific absorption rate (asr) in mri
US20070108976A1 (en) * 2003-10-31 2007-05-17 Koninkijke Hpilips Electronics Nv B1 field control in magnetic resonance imaging
US20070241753A1 (en) * 2006-02-21 2007-10-18 Beth Israel Deaconess Medical Center, Inc. Magnetic resonance imaging and radio frequency impedance mapping methods and apparatus

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6249121B1 (en) * 1999-05-17 2001-06-19 General Electric Company RF body coil
DE10150137B4 (en) * 2001-10-11 2006-08-17 Siemens Ag Method and apparatus for magnetic resonance imaging
DE10150138B4 (en) * 2001-10-11 2009-10-08 Siemens Ag Method for magnetic resonance imaging
US6989673B2 (en) * 2003-11-26 2006-01-24 General Electric Company Method and apparatus to reduce RF power deposition during MR data acquisition
DE102005019859A1 (en) * 2004-05-11 2005-12-08 Siemens Ag Method for carrying out a magnetic resonance examination of a patient
US8099059B2 (en) * 2005-04-29 2012-01-17 Koninklijke Philips Electronics N.V. Method and circuit arrangement for operating multi-channel transmit/recieve antenna devices
US8238678B2 (en) * 2006-08-30 2012-08-07 Siemens Medical Solutions Usa, Inc. Providing representative image information
DE102006052711A1 (en) * 2006-11-08 2008-05-15 Siemens Ag Method and device for surface detection of an examination person
DE102007013060A1 (en) * 2007-03-19 2008-09-25 Siemens Ag Method for setting a B1 field

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005001502A1 (en) * 2003-06-30 2005-01-06 Koninklijke Philips Electronics N.V. Control of specific absorption rate (asr) in mri
US20070108976A1 (en) * 2003-10-31 2007-05-17 Koninkijke Hpilips Electronics Nv B1 field control in magnetic resonance imaging
US20070241753A1 (en) * 2006-02-21 2007-10-18 Beth Israel Deaconess Medical Center, Inc. Magnetic resonance imaging and radio frequency impedance mapping methods and apparatus

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
A.C.ZELINSKI ET AL.: "Designing RF Pulses With Optimal Specific Absorption Rate (SAR) Characteristics And Exploring Excitation Fidelity, SAR And Pulse Duration Tradeoffs", PROC.INTL.SOC.MAG.RESON.MED. 15, 2007, pages 1699, XP002533633 *
A.C.ZELINSKI ET AL.: "Pulse Design Methods For Reduction Of Specific Absorption Rate In Parallel RF Excitation", PROC.INTL.SOC.MAG.RESON.MED. 15, 2007, pages 1698, XP002533632 *
C.A.T.VAN DEN BERG ET AL.: "Simultaneous B1+ Homogenization and Specific Absorption Rate Hotspot Suppression Using a Magnetic Resonance Phased Array Transmit Coil", MAGNETIC RESONANCE IN MEDICINE, vol. 57, 2007, pages 577 - 586, XP002533629 *
D.O.BRUNNER ET AL.: "Enforcing strict constraints in multiple-channel RF pulse optimization", PROC.INTL.SOC.MAG.RESON.MED. 15, 2007, pages 1690, XP002533631 *
F.SEIFERT ET AL.: "Patient Safety Concept for Multichannel Transmit Coils", JOURNAL OF MAGNETIC RESONANCE IMAGING, vol. 26, 2007, pages 1315 - 1321, XP002533634 *
I.GRAESSLIN ET AL.: "A Minimum SAR RF Pulse Design Approach for Parallel Tx with Local Hot Spot Suppression and Exact Fidelity Constraint", PROC.INTL.SOC.MAG.RESON.MED. 16, May 2008 (2008-05-01), pages 621, XP002533628 *
T.S.IBRAHIM ET AL.: "A Whole-Body 7 Tesla RF Excitation Scheme with Much Improved B1+ Field Homogeneity and Local/Global SARs over Quadrature Excitation", PROC.INTL.SOC.MAG.RESON.MED. 15, 2007, pages 1013, XP002533630 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102010015066B4 (en) * 2010-04-15 2014-12-11 Siemens Aktiengesellschaft Method and device for determining a magnetic resonance system drive sequence and method for operating a magnetic resonance system
CN103154761A (en) * 2010-08-20 2013-06-12 皇家飞利浦电子股份有限公司 Virtual coil emulation in parallel transmission mri
US9689945B2 (en) 2010-08-20 2017-06-27 Koninklijke Philips Electronics N.V. Virtual coil emulation in parallel transmission MRI
CN105093144A (en) * 2014-05-22 2015-11-25 西门子公司 Method and apparatus for generating magnetic resonance image using combined excitation by whole body coil and local coil
US10048344B2 (en) 2014-05-22 2018-08-14 Siemens Aktiengesellschaft Method and apparatus for generating a magnetic resonance image using combined excitation by a whole body coil and a local coil

Also Published As

Publication number Publication date
JP2011517983A (en) 2011-06-23
US20110043205A1 (en) 2011-02-24
US8941380B2 (en) 2015-01-27
CN102007423B (en) 2015-11-25
JP5731373B2 (en) 2015-06-10
EP2269085A1 (en) 2011-01-05
CN102007423A (en) 2011-04-06

Similar Documents

Publication Publication Date Title
US8941380B2 (en) Real-time local and global SAR estimation for patient safety and improved scanning performance
US9835700B2 (en) Magnetic resonance imaging apparatus
US8030920B2 (en) Method and system for modifying pulse sequences
RU2592039C2 (en) Formation of magnetic resonance images using multiple-point dixon method
US10393840B2 (en) Magnetic resonance apparatus and method for acquiring measurement data during simultaneous manipulation of spatially separate sub-volumes
US8938281B2 (en) MR imaging using multi-channel RF excitation
US9689945B2 (en) Virtual coil emulation in parallel transmission MRI
WO2013165571A1 (en) System and method for quiet magnetic resonance imaging
EP3191862B1 (en) Zero echo time mr imaging
US7358731B2 (en) Determination of spatial sensitivity profiles of RF coils in magnetic resonance imaging
EP1445623B1 (en) Eddy current correction method and magnetic resonance imaging apparatus
US8456163B2 (en) Variable image resolution MR imaging system
EP2503348A1 (en) MRI using a multi-point Dixon technique and a low resolution calibration scan
CN105093144B (en) Magnetic resonance image creation by means of a combination of a whole-body coil and a local transmission coil
US10928469B2 (en) Magnetic resonance imaging apparatus
JP3802891B2 (en) Gain adjusting method and magnetic resonance imaging apparatus
US20090099443A1 (en) Magnetic resonance device and method
JP3753668B2 (en) RF pulse tuning device
JP2023516474A (en) Phase encoding using frequency-swept pulses for magnetic resonance imaging in inhomogeneous magnetic fields
US20220206097A1 (en) Magnetic resonance image processing method and device, and computer readable storage medium
CN111090067B (en) Method and system for sampling K-space data in magnetic resonance imaging

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 200980113091.7

Country of ref document: CN

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 09732227

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2009732227

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2011504587

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 7172/CHENP/2010

Country of ref document: IN

WWE Wipo information: entry into national phase

Ref document number: 12936540

Country of ref document: US