SYSTEM COMPONENT IN AN IMAGING SYSTEM
20200150161 · 2020-05-14
Inventors
- Thorsten Feiweier (Poxdorf, DE)
- Michael Köhler (Nürnberg, DE)
- Helmut Lenz (Oberasbach, DE)
- Sören Grübel (Erlangen, DE)
Cpc classification
G01R33/3852
PHYSICS
G01R33/543
PHYSICS
G01R33/36
PHYSICS
International classification
Abstract
Systems and methods are provided for determining a use of a system component in an imaging system. The imaging system includes a primary side configured to provide power to the system component and a secondary side including the system component that uses the power provided by the primary side during the image sequence. The method includes determining the use of the system component during an imaging sequence, determining a time averaged power provided by the primary side during the imaging sequence, determining a maximum time averaged power that may be provided by the primary side until a temperature limit is reached on the primary side. Further, whether the time averaged power is smaller than the maximum time averaged power is determined.
Claims
1. A method for determining a utilization of a system component in an imaging system in which at least one image of an object under examination is generated during an imaging sequence, the imaging system comprising a primary side configured to provide power to the system component resulting in a thermal load on the primary side that must not exceed a predefined temperature limit, and a secondary side comprising the system component that utilizes the power provided by the primary side during the image sequence, the method comprising: determining the utilization of the system component during the imaging sequence, determining a time averaged power supplied by the primary side during the imaging sequence with the determined utilization; determining a maximum time averaged power that may be supplied by the primary side over a duration of at least one imaging sequence while not exceeding the predefined temperature limit; and determining whether the time averaged power is smaller than the maximum time averaged power, wherein when the time averaged power is not smaller than the maximum time averaged power, configuring the use of the system component during the imaging sequence until the time averaged power is smaller than the maximum time averaged power.
2. The method of claim 1 wherein the time averaged power is determined as a function of where a first time constant that describes a power induced heating of the primary side when power is supplied to the system component is longer than a time period in which the system component is continuously utilized during the imaging sequence.
3. The method of claim 1, wherein the utilization of the system component is determined without identifying a relationship of how the use of the system component during the imaging sequence influences a heating of the primary side.
4. The method of claim 1, further comprising: determining an optimized imaging sequence based on a comparison of the time averaged power to the maximum time averaged power such that a time needed to generate the at least one image, an image quality parameter of the at least one image, or the time needed to generate the at least one image and the image quality parameter of the at least one image is optimized for the imaging sequence.
5. The method of claim 1, wherein determining the time averaged power comprises determining an average of a parameter describing a heating of the primary side when power is provided to the system component, and determining the maximum time averaged power comprises determining a maximum of the average of the parameter.
6. The method according to claim 5, wherein the parameter is a square of a current provided by the primary side during the imaging sequence.
7. The method of claim 1, wherein the imaging system is an MR system configured to generate MR images, wherein the time averaged power provided by the primary side during the imaging sequence is determined when the time averaged power is determined.
8. The method of claim 7, wherein the system component comprises a gradient field generating unit used to generate magnetic field gradients applied in the MR system, wherein the utilization of the magnetic field gradients in the imaging sequence is determined and an average of a square current provided by the primary side to set up the magnetic field gradients during the imaging sequence is determined.
9. The method of claim 7, wherein the square current is determined taking into account an offset current which is flowing independent of whether magnetic field gradients are applied in the imaging sequence.
10. The method of claim 8, further comprising: determining at least one parameter a describing a relationship between the applied magnetic field gradients and a current provided by the primary side when the magnetic field gradient is applied, wherein the average of the square current is determined taking into account the determined at least one parameter a.
11. The method of claim 7, wherein the imaging sequence is a diffusion imaging sequence used to determine a diffusion property in the object under examination.
12. The method of claim 7, wherein a plurality of sequential imaging sequences are used in the MR system, wherein a utilization of magnetic field gradients applied in the plurality of sequential imaging sequences is determined taking into account a corresponding time averaged power for each of the plurality of imaging sequences and the maximum time averaged power.
13. The method of claim 1, wherein the time averaged power is determined by averaging the power provided by the primary side over an averaging period T, wherein the averaging period is larger than a time period in which the system component is continuously switched during the imaging sequence.
14. The method of claim 1, wherein the time averaged power is determined based on the utilization of the system component based on a model which translates the utilization of the system component in the power needed to use the system component.
15. An imaging system configured to generate at least one image of an object under examination during an imaging sequence, the system comprising: a system component configured to be switched on and off during the imaging sequence in order to generate the at least one image; a primary side configured to provide power to the system component, resulting in a thermal load on the primary side which must not exceed a predefined temperature limit; a secondary side comprising the system component that is configured to use the power provided by the primary side during the imaging sequence; and a control unit configured to: determine a utilization of the system component during the imaging sequence; determine a time averaged power supplied by the primary side during the imaging sequence with the determined utilization; determine a maximum time averaged power that may be supplied by the primary side over a duration of at least one imaging sequence while not exceeding the predefined temperature limit; determine whether the time averaged power is smaller than the maximum time averaged power, wherein when the time averaged power is not smaller than the maximum time averaged power, adapting the use of the system component during the imaging sequence until the time averaged power is smaller than the maximum time averaged power.
16. The imaging system of claim 15, wherein the control unit determines the time averaged power as a function of where a first time constant that describes a power induced heating of the primary side when power is supplied to the system component is longer than a time period in which the system component is continuously utilized during the imaging sequence.
17. The imaging system of claim 15, wherein the utilization of the system component is determined without identifying a relationship of how the use of the system component during the imaging sequence influences a heating of the primary side.
18. The imaging system of claim 15, wherein the control unit is further configured to determine an optimized imaging sequence based on a comparison of the time averaged power to the maximum time averaged power such that a time needed to generate the at least one image, an image quality parameter of the at least one image, or the time needed to generate the at least one image and the image quality parameter of the at least one image is optimized for the imaging sequence.
19. A non-transitory computer implemented storage medium that stores machine-readable instructions executable by at least one processor, the machine-readable instructions comprising: determining a utilization of a system component in an imaging system in which at least one image of an object under examination is generated during an imaging sequence, the imaging system comprising a primary side configured to provide power to the system component resulting in a thermal load on the primary side that must not exceed a predefined temperature limit, and a secondary side comprising the system component that utilizes the power provided by the primary side during the image sequence; determining a time averaged power supplied by the primary side during the imaging sequence with the determined utilization; determining a maximum time averaged power that may be supplied by the primary side over a duration of at least one imaging sequence while not exceeding the predefined temperature limit; and determining whether the time averaged power is smaller than the maximum time averaged power, wherein when the time averaged power is not smaller than the maximum time averaged power, configuring the use of the system component during the imaging sequence until the time averaged power is smaller than the maximum time averaged power.
20. The non-transitory computer implemented storage medium of claim 19, wherein determining the time averaged power comprises determining an average of a parameter describing a heating of the primary side when power is provided to the system component, and determining the maximum time averaged power comprises determining a maximum of the average of the parameter.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
[0031]
[0032]
[0033]
DETAILED DESCRIPTION
[0034] In the following, embodiments are described in detail with reference to the accompanying drawings. The drawings are to be regarded as being schematic representations, and elements illustrated in the drawings are not necessarily shown to scale. Rather, the various elements are represented such that their function and general purpose becomes apparent to a person skilled in the art. Any connection or coupling between functional blocks, devices, components of physical or functional units depicted in the drawings and described hereinafter may be implemented by an indirect connection or coupling. A coupling between components may be established over a wired or wireless connection. Functional blocks may be implemented in hardware, software, firmware, or a combination thereof.
[0035]
[0036] The MR system 1 includes a magnet 10 generating a polarization field B0. An object under examination 12 lying on a table 11 is moved into the center of the MR system 1 where MR signals of the RF-excitation may be detected by a receiving coil 2. One example of a gradient coil 5 is depicted that is able to generate a magnetic field gradient in one direction. Several of the gradient coils may be provided to generate the magnetic field gradients in different directions. A transmitting coil 3 is depicted that is configured to transmit RF pulses into the object under examination. By applying RF pulses and magnetic field gradients, the nuclear spins in the object 12 are excited and the currents induced by the relaxation is detected. The way how MR images are generated and how MR signals are detected using a sequence of RF pulses and a sequence of magnetic field gradients are known.
[0037] The MR system 1 includes a controller 13 that is used for controlling the MR system. The controller or control module 13 includes a gradient control unit 14 for controlling and switching the magnetic field gradients, an RF control unit 15 for controlling and generating the RF pulses for the imaging sequence. An imaging sequence control unit 16 is provided that controls the sequence of the applied RF pulses and the magnetic field gradients and thus controls the gradient control unit 14 into the RF control unit 15. In a memory 17 computer programs needed for operating the MR system and the imaging sequences necessary for generating the MR images may be stored together with the generated MR images. The generated MR images may be displayed on a display 18. An input unit 19 may be provided used by users of the MR system to control the functioning of the MR system. A processing unit 20 is provided and may coordinate the operation of the different functional units shown in
[0038] The MR system may be a spectroscopic MR system that is configured to detect MR signals of the object under examination without generating an image, but for generating spectroscopic data of the object under examination. The system component on the primary side may be the RF unit generating the RF pulses needed to excite the MR signal or the gradient unit in case of a single voxel spectroscopy.
[0039]
[0040] The protective switch 22 is configured such that it passes current to the secondary side until a certain temperature limit is reached within the protective switch 22. Thus, the protective switch 22 plays the role of a current limiter that limits the current that may be provided to the secondary side over time. The idea described below is based on that the time constant that describes the current induced heating of the protective switch 22 or the heating of the primary side itself is much longer than the time period in which the different gradient coils 31, 32, and 33 are switched during the imaging sequence to obtain the magnetic field gradients that are needed for the spatial encoding of the signals. There is no need to exactly identify the relationship between the actual use of the gradient coils and the thermal heating of the current limiter 22 provided on the primary side. As will be described below a model is used that is configured to determine a load parameter on the primary side that is needed to drive the magnetic field gradient based on the exact switching of the magnetic field gradients as identified from the used imaging sequence.
[0041] Embodiments make use of a time-averaged power as load parameter that is used by the magnetic field gradient and the maximum time averaged power that may be provided by the primary side over time and as a result determine optimized imaging sequences that may be used by the MR imaging system without exceeding the temperature limit on the primary side.
[0042] Furthermore, the knowledge may be used how a pause in the imaging sequence influences the heating of the current limiter.
[0043] The model described below and discussed above in connection with
[0044] The current needs of this part of the MR system depend on the switching of the magnetic field gradients by the imaging sequence. The switching of the current on the secondary side named as I.sub.X, I.sub.Y and I.sub.Z for the three axes X, Y, Z of the gradient system leads to a current on the primary side named I.sub.1, I.sub.2, and I.sub.3 in
[0045] The square of the current on the primary side I.sub.pri that is summed over the three different phases may be described as follows:
I.sup.2.sub.pri(t)=(I.sub.Offset+I.sub.Grad(t)).sup.2(1)
[0046] An offset current I.sub.offset is assumed that describes the current that is used by the gradient or MR system even when no actual magnetic field gradient is applied in the imaging sequence. Furthermore, a gradient related current I.sub.Grad is considered that depends on the actual switching of the magnetic field gradients and that may be mathematically describedassuming a transfer function as explained abovewith a convolution as shown in the following equation:
I.sub.Grad(t)=.sub.j=x,y,z.sub.jG.sub.j.sup.2(t)C(t)(2)
[0047] The above equations show that the load on the primary side is scaling with the fourth power of the current on the secondary side. It is also possible to include alpha into the transfer function and to consider axis-specific transfer functions.
[0048] The time constants on the primary side, here the time constants of the protective switch 22 are in the range of several 10 seconds until several minutes and describe the time constant of the current induced heating until a maximum temperature limit is obtained at which the current limiter interrupts the provision of currents. The time constant is much longer than the time period in which the magnetic field gradients are switched on which is in the range of milliseconds before they are switched off and switched on again. Accordingly, when considering the relevant thermal load on the primary side, it is sufficient to only take into account average values as shown in the following equation:
<I.sup.2.sub.pri>=1/T.sub.0.sup.TdtI.sup.2.sub.pri(t)(3)
[0049] The averaging time T used in equation 3 is much longer than the time period in which a single gradient is switched on before it is switched off again. The averaging time T is selected such that the resulting average value as determined by equation 3 does not depend on the fact which time period in a temporal evolution of the imaging sequence is selected to do the averaging. The time period T is in the same range as the time period describing the current induced heating on the primary side.
[0050] The average of the square current as determined by equation 3 is then compared to a maximum value, determined by the maximum temperature under which the protective switch 22 is operating without interrupting the current provided to the secondary side. When the average square of the current is lower than the square of the maximum current, the imaging sequence may get executed with the determined use of the magnetic field gradients. If not, an adaptation of the magnetic field gradients may be necessary, either by introducing further pauses into the imaging sequence or by reducing the amplitude of the magnetic field gradients.
[0051] In the example given above actual values of I.sub.offset and the parameter a are needed. Those values may be determined directly from the design of the gradient components. If this is not possible or too complex, it is also possible to determine both parameters by calibration measurements. Two different currents provided by the primary side are measured when different imaging sequences are executed on the secondary side with a known switching pattern of magnetic field gradients. It should be understood that the current used in the two different imaging sequences should differ by a certain amount in order to have different calibration points. With the knowledge of the used magnetic field gradients G.sub.x,y,z(t) it is possible to determine the needed parameters. Differences between the different gradient axes may either be determined by averaging based on an assumption that a similar switching pattern is used on each axis or may be determined using a relative scaling when the differences for the different gradient axes are known. Furthermore, it is possible to obtain a separate calibration for each gradient axis.
[0052] With the model, the knowledge of the switching of the magnetic field gradients in the imaging sequence and the calibration it is possible to determine the currents that have to be provided by the primary side. If the currents and thus the current induced heating are larger than the threshold provided by the current limiter, the imaging sequence may be adapted accordingly. By way of example it is possible to introduce a pause of a certain time period into the imaging sequence so that the time needed to carry out the imaging sequence is increased. As an alternative, or in addition, the amplitudes of the magnetic field gradients may be reduced.
[0053] One possible field of application is the use of the above method in diffusion imaging in which diffusion properties of a certain part of the examined body are determined. For diffusion imaging high magnetic field gradients are required that get successively switch along different directions that may lead to a large load on the primary side. However, as the time periods in which the diffusion gradients are switched on are comparatively short, the above described method may get applied. Accordingly, it is possible to determine in advance whether an image acquisition with the selected gradient switching will be possible without exceeding the temperature limit. This provides the calculation of valid parameter ranges for b-values, echo times TE or repetition times TR, and to limit the selection of parameters to the ranges. The user of the system may thus only select the parameter within the predetermined range that limits the maximum heat load and assures that the temperature limit is not exceeded. In order to stay within the temperature limit, it may be necessary to increase the echo time TE that leads to a smaller gradient for a certain b-value and thus to a smaller current, wherein the longer TE may lead to a smaller signal-to-noise ratio. As an alternative it may be necessary to increase the repetition time TR that increases the measurement duration but that also leads to a reduced current below the limit on the primary side.
[0054] The aforementioned approach may also be used when several sequential imaging sequences are applied to the object under examination. In an MR imaging of an object under examination, different imaging sequences are used, some of them have high gradient demand whereas others have lower gradient demand. As long as the time periods of the actual switching of the magnetic field gradients is much shorter than the time constant describing the current induced heating and as long as the overall time of the imaging sequences with large magnetic field gradient demand is shorter than the time constant describing the current induced heating, the above described method may increase the performance of the whole MR system.
[0055] As far as the calibration is concerned, the calibration may be determined for each individual MR imaging system before delivery or during installation, or for each type of the MR imaging system.
[0056] By way of example, the type-specific calibration relating to a certain type of MR imaging systems, e.g. having a certain magnetic field strength or certain gradient components may be carried out in the factory and pre-coded into the system so that a calibration is not necessary each time a system gets installed.
[0057] For components that contain complex electric circuitries, it might turn out that the precision of the model predictions depends on the actual type of imaging sequences. By way of example the prediction of the required currents with high precision might be possible either for sequences with strong magnetic field gradient variations or for nearly constant magnetic field gradients, and this may depend on the fact whether the calibration was obtained with one or the other type of sequences. The calibration may be determining, for example, for the class of imaging sequences for which the reaching of the temperature limit is highly likely. Different calibrations may be carried out for the different classes of imaging sequences.
[0058] If the MR systems is used in such a way that the user may only select imaging parameters within a parameter range by which the limits of the current limiter is not exceeded, it may be necessary to determine and predict the temperature induced heating very quickly. Certain assumptions for the calculation of the magnetic field gradients may be used. By way of example, a normal gradient pulse may have a trapezoidal form with an amplitude G, a ramp time TR and a constant maximum value during T.sub.D. The trapezoidal shapes may be simplified with a gradient having a square shape with a certain amplitude G and a switching time T.sub.K. The time period T.sub.K may be determined based on T.sub.R and T.sub.D such that the resulting current need corresponds to the actual use of the currents. By way of example time period T.sub.K may be determined as follows:
T.sub.K=T.sub.R*+TD(4)
[0059] Furthermore, it is possible to take into account the use of a pause within the imaging sequence. If the average load is described as an integral over the load B(t), then B is determined as follows:
B=1/T.sub.0.sup.TdtB(t)(5)
[0060] A pause of the duration TP leads to the following average load B:
B=1/(T+TP).sub.0.sup.T+TPdtB(t)=1/(T+TP).sub.0.sup.TdtB(t)=BT/(T+TP) (6)
[0061] Based on the knowledge of the average load B and the time period T, it is possible to directly determine the pause needed to stay below the limit B.sub.maxB.
[0062] The way how a pause influences the mean load may be determined numerically. It is possible to introduce longer and shorter pauses into the imaging sequence and then to determine the average load.
[0063]
[0064] Summarizing the above described method makes it possible to operate the MR system close to the temperature limits without exceeding them. Accordingly, the imaging system may be operated without interrupting the measurements due to the fact that temperature limits by the system components have been exceeded.
[0065] The use of the average values makes it possible to model very complex relationships. Furthermore, as discussed, the current provided by the primary side is determined based on a model. The model may additionally use the characteristics of the gradient amplifier by using calibration measurements in order to determine a relation between the current used by the secondary side to switch the system component and the current provided by the primary side.
[0066] It is to be understood that the elements and features recited in the appended claims may be combined in different ways to produce new claims that likewise fall within the scope of the present invention. Thus, whereas the dependent claims appended below depend from only a single independent or dependent claim, it is to be understood that these dependent claims may, alternatively, be made to depend in the alternative from any preceding or following claim, whether independent or dependent, and that such new combinations are to be understood as forming a part of the present specification.
[0067] While the present invention has been described above by reference to various embodiments, it may be understood that many changes and modifications may be made to the described embodiments. It is therefore intended that the foregoing description be regarded as illustrative rather than limiting, and that it be understood that all equivalents and/or combinations of embodiments are intended to be included in this description.