MAGNETIC RESONANCE IMAGING USING CORRECTED K-SPACE TRAJECTORIES CALCULATED FROM CURRENT SENSOR DATA

Abstract

The invention provides for a magnetic resonance imaging system (100, 300, 500) with a gradient coil system (110, 112, 113) that comprises a set of gradient coils (110) configured for generating a gradient, a gradient coil amplifier (112), and a current sensor system (113) configured for measuring current sensor data (146) descriptive of the electrical current supplied to each of the set of gradient coils. Execution of the machine executable instructions causes a processor to: control (200) the magnetic resonance imaging system with the pulse sequence commands (142) to acquire magnetic resonance imaging data; record (202) the current sensor data during the acquisition of the magnetic resonance imaging data; calculate (204) a corrected k-space trajectory (150) using the current sensor data and a gradient coil transfer function (148); and reconstruct (206) a corrected magnetic resonance image (152) using the magnetic resonance imaging data and the corrected k-space trajectory.

Claims

1. A magnetic resonance imaging system configured for acquiring magnetic resonance imaging data from an imaging zone, the magnetic resonance imaging system comprising: a magnet for generating a main magnetic field within the imaging zone; a gradient coil system for generating a gradient magnetic field within the imaging zone, wherein the gradient coil system further comprises a set of gradient coils for generating the gradient magnetic field when supplied with electrical current, wherein each of the set of gradient coils is configured for generating the gradient magnetic field along an axis, wherein the gradient coil system comprises a gradient coil amplifier configured for supplying the electrical current to each of the set of gradient coils, and wherein the gradient coil system further comprises a current sensor system configured for measuring current sensor data descriptive of the electrical current supplied to each of the set of gradient coils by the gradient coil amplifier; a memory containing machine executable instructions, wherein the memory further contains pulse sequence commands configured for controlling the magnetic resonance imaging system to acquire the magnetic resonance imaging data according to a magnetic resonance imaging protocol, wherein the memory further contains a selected gradient coil transfer function configured for mapping the current sensor data to magnetic field components within the imaging zone, wherein the magnetic field components comprise the gradient magnetic field; a processor for controlling the magnetic resonance imaging system, wherein execution of the machine executable instructions causes the processor to: control the magnetic resonance imaging system with the pulse sequence commands to acquire the magnetic resonance imaging data; record the current sensor data during the acquisition of the magnetic resonance imaging data; calculate a corrected k-space trajectory using the current sensor data and the selected gradient coil transfer function; and reconstruct a corrected magnetic resonance image using the magnetic resonance imaging data and the corrected k-space trajectory according to the magnetic resonance imaging protocol; wherein the gradient amplifier is configured to receive a gradient control signal, wherein the memory further contains a selected gradient amplifier transfer function configured for mapping the gradient control signal to the electrical current supplied by the gradient amplifier, wherein the pulse sequence commands are configured to provide the control signal during acquisition of the magnetic resonance imaging data, wherein execution of the machine executable instructions further cause the processor to calculate a corrected control signal using the control signal and the selected gradient amplifier transfer function, wherein the gradient amplifier is configured to be controlled with the corrected control signal during acquisition of the magnetic resonance imaging data.

2. The magnetic resonance imaging system of claim 1, wherein the memory further comprises a set of gradient coil transfer functions, wherein execution of the machine executable instructions further causes the processor to receive one or more acquisition parameters descriptive of the magnetic resonance imaging protocol, wherein execution of the machine executable instructions further causes the processor to choose the selected gradient coil transfer function using the acquisition parameters.

3. The magnetic resonance imaging system of claim 2, wherein the acquisition parameters comprise any one of the following: a subject height, a subject weight, a subject support position, a room temperature, a magnetic resonance imaging protocol type, a type of receive coil, and combinations thereof.

4. The magnetic resonance imaging system of claim 3, wherein the acquisition parameters further comprise any one of the following: a cryostat temperature, a cryostat state, a gradient coil coolant temperature, a gradient coil temperature, a gradient coil impedance, and combinations thereof.

5. The magnetic resonance imaging system of claim 1, wherein execution of the machine executable instructions further causes the processor to control the magnetic resonance imaging system to: measure calibration current sensor data using the current sensor system and a gradient coil system response using the magnetic resonance imaging system; and calculate the selected gradient coil transfer function using the calibration current sensor data and the gradient coil system response before acquiring the magnetic resonance imaging data.

6. The magnetic resonance imaging system of claim 1, wherein execution of the machine executable instructions further causes the processor to modify the selected gradient amplifier transfer function using the current sensor data recorded during acquisition of the magnetic resonance imaging data.

7. The magnetic resonance imaging system of claim 6, wherein the memory further comprises a set of gradient amplifier transfer functions, wherein execution of the machine executable instructions further causes the processor to receive one or more system status parameters descriptive of a status of the magnetic resonance imaging system, wherein execution of the machine executable instructions further causes the processor to choose the selected gradient amplifier transfer function from the set of gradient amplifier transfer functions using the one or more system status parameters.

8. The magnetic resonance imaging system of claim 7, wherein the one or more system status parameters comprise any one of the following: a room temperature, a prior use of the gradient amplifier, a gradient amplifier temperature, a magnetic resonance imaging protocol type, a gradient coil temperature, a gradient coil impedance, and combinations thereof.

9. The magnetic resonance imaging system of claim 7, wherein execution of the machine executable instructions further causes the processor to store the selected gradient amplifier transfer function after modifying the set of gradient amplifier transfer functions.

10. The magnetic resonance imaging system of claim 9, wherein the selected gradient amplifier transfer function is selected from the set of gradient amplifier transfer functions using a machine learning algorithm.

11. The magnetic resonance imaging system of claim 10, wherein execution of the machine executable instructions further causes the processor to train the machine learning algorithm with the system status parameters when storing the selected gradient amplifier transfer function in the set of gradient amplifier transfer functions.

12. The magnetic resonance imaging system of claim 1, wherein the corrected magnetic resonance imaging is at least partially reconstructed by regridding the magnetic resonance imaging data using the corrected k-space trajectory.

3. A method of operating a magnetic resonance imaging system configured for acquiring magnetic resonance imaging data from an imaging zone, wherein the magnetic resonance imaging system comprises a magnet for generating a main magnetic field within the imaging zone, wherein the magnetic resonance imaging system further comprises a gradient coil system for generating a gradient magnetic field within the imaging zone, wherein the gradient coil system further comprises a set of gradient coils for generating the gradient magnetic field when supplied with electrical current, wherein each of the set of gradient coils is configured for generating the gradient magnetic field along an axis, wherein the gradient coil system comprises a gradient coil amplifier configured for supplying the electrical current to each of the set of gradient coils, and wherein the gradient coil system further comprises a current sensor system configured for measuring current sensor data descriptive of the electrical current supplied to each of the set of gradient coils by the gradient coil amplifier, wherein the method comprises: controlling the magnetic resonance imaging system with pulse sequence commands to acquire the magnetic resonance imaging data according to a magnetic resonance imaging protocol; recording the current sensor data during the acquisition of the magnetic resonance imaging data; calculating a corrected k-space trajectory using the current sensor data and a selected gradient coil transfer function, wherein the selected gradient coil transfer function is configured for mapping the current sensor data to magnetic field components within the imaging zone, wherein the magnetic field components comprise the gradient magnetic field; and reconstructing a corrected magnetic resonance image using the magnetic resonance imaging data and the corrected k-space trajectory according to the magnetic resonance imaging protocol; wherein the gradient amplifier is configured to receive a gradient control signal, wherein the memory further contains a selected gradient amplifier transfer function configured for mapping the gradient control signal to the electrical current supplied by the gradient amplifier, wherein the pulse sequence commands are configured to provide the control signal during acquisition of the magnetic resonance imaging data, wherein execution of the machine executable instructions further cause the processor to calculate a corrected control signal using the control signal and the selected gradient amplifier transfer function, wherein the gradient amplifier is configured to be controlled with the corrected control signal during acquisition of the magnetic resonance imaging data.

14. A computer program product comprising machine executable instructions for execution by a processor controlling a magnetic resonance imaging system configured for acquiring magnetic resonance imaging data from an imaging zone, wherein the magnetic resonance imaging system comprise a magnet for generating a main magnetic field within the imaging zone, wherein the magnetic resonance imaging system further comprise a gradient coil system for generating a gradient magnetic field within the imaging zone, wherein the gradient coil system further comprises a set of gradient coils for generating the gradient magnetic field when supplied with electrical current, wherein each of the set of gradient coils is configured for generating the gradient magnetic field along an axis, wherein the gradient coil system comprises a gradient coil amplifier configured for supplying the electrical current to each of the set of gradient coils, and wherein the gradient coil system further comprises a current sensor system configured for measuring current sensor data descriptive of the electrical current supplied to each of the set of gradient coils by the gradient coil amplifier, wherein execution of the machine executable instructions causes the processor to: control the magnetic resonance imaging system with pulse sequence commands to acquire the magnetic resonance imaging data according to a magnetic resonance imaging protocol; record the current sensor data during the acquisition of the magnetic resonance imaging data; calculate a corrected k-space trajectory using the current sensor data and a selected gradient coil transfer function, wherein the selected gradient coil transfer function is configured for mapping the current sensor data to magnetic field components within the imaging zone, wherein the magnetic field components comprise the gradient magnetic field; and reconstruct a corrected magnetic resonance image using the magnetic resonance imaging data and the corrected k-space trajectory according to the magnetic resonance imaging protocol; wherein the gradient amplifier is configured to receive a gradient control signal, wherein the memory further contains a selected gradient amplifier transfer function configured for mapping the gradient control signal to the electrical current supplied by the gradient amplifier, wherein the pulse sequence commands are configured to provide the control signal during acquisition of the magnetic resonance imaging data, wherein execution of the machine executable instructions further cause the processor to calculate a corrected control signal using the control signal and the selected gradient amplifier transfer function, wherein the gradient amplifier is configured to be controlled with the corrected control signal during acquisition of the magnetic resonance imaging data.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0057] In the following preferred embodiments of the invention will be described, by way of example only, and with reference to the drawings in which:

[0058] FIG. 1 illustrates an example of a magnetic resonance imaging system;

[0059] FIG. 2 shows a flow chart which illustrates a method of using the magnetic resonance imaging system of FIG. 1;

[0060] FIG. 3 illustrates a further example of a magnetic resonance imaging system;

[0061] FIG. 4 shows a flow chart which illustrates a method of using the magnetic resonance imaging system of FIG. 4;

[0062] FIG. 5 illustrates a further example of a magnetic resonance imaging system;

[0063] FIG. 6 shows a flow chart which illustrates a method of using the magnetic resonance imaging system of FIG. 5;

[0064] FIG. 7 illustrates an example of electrical sensors used to measure the input and output of a gradient coil amplifier;

[0065] FIG. 8 illustrates the use of the gradient coil transfer function and the gradient amplifier transfer function;

[0066] FIG. 9 illustrates an example of a gradient coil transfer function;

[0067] FIG. 10 shows a block diagram which illustrates how the input to the gradient coil amplifier and the output of the gradient coil amplifier can be used to calculate a gradient amplifier transfer function;

[0068] FIG. 11 shows a sketch of the real component of a gradient amplifier transfer function; and

[0069] FIG. 12 shows a sketch of the phase of a gradient amplifier transfer function.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0070] Like numbered elements in these figures are either equivalent elements or perform the same function. Elements which have been discussed previously will not necessarily be discussed in later figures if the function is equivalent.

[0071] FIG. 1 illustrates an example of a magnetic resonance imaging system 100 with a magnet 104. The magnet 104 is a superconducting cylindrical type magnet with a bore 106 through it. The use of different types of magnets is also possible; for instance it is also possible to use both a split cylindrical magnet and a so called open magnet. A split cylindrical magnet is similar to a standard cylindrical magnet, except that the cryostat has been split into two sections to allow access to the iso-plane of the magnet, such magnets may for instance be used in conjunction with charged particle beam therapy. An open magnet has two magnet sections, one above the other with a space in-between that is large enough to receive a subject: the arrangement of the two sections area similar to that of a Helmholtz coil. Open magnets are popular, because the subject is less confined. Inside the cryostat of the cylindrical magnet there is a collection of superconducting coils. Within the bore 106 of the cylindrical magnet 104 there is an imaging zone 108 where the magnetic field is strong and uniform enough to perform magnetic resonance imaging. A region of interest 109 is shown within the imaging zone 108. The magnetic resonance imaging data that is typically acquired for the region of interest. A subject 118 is shown as being supported by a subject support 120 such that at least a portion of the subject 118 is within the imaging zone 108 and the region of interest 109.

[0072] Within the bore 106 of the magnet there is also a set of set of gradient coils 110 which is used for acquisition of preliminary magnetic resonance imaging data to spatially encode magnetic spins within the imaging zone 108 of the magnet 104. The set of gradient coils 110 connected to a magnetic field gradient coil amplifier 112. The set of gradient coils 110 are intended to be representative. The set of gradient coils 110 contain three separate coils for spatially encoding in three orthogonal spatial directions. A magnetic field gradient power supply supplies current to the set of gradient coils. The current supplied to the set of gradient coils 110 is controlled as a function of time and may be ramped or pulsed.

[0073] The gradient coils 110 represent three separate sets of orthogonal gradient coils for generating a gradient magnetic field within the imaging zone 108. These are typically oriented as the axes 122, 123 and 124 show. The axis 124 is aligned with the axis of the magnet 104. This is typically referred to as the z-axis. 122 and 123 are the x and y-axes respectively. They are orthogonal to each other and also to the z-axis 124.

[0074] The magnetic field gradient coil amplifier 112 is configured for supplying current to each of the sets of gradient coils separately. The magnetic field gradient coil amplifier 112 is shown as having a current sensor system 113 for measuring the current supplied to each of the set of gradient coils 110. The current sensor system 113 could for example be part of the magnetic field gradient coil amplifier 112 or it could also be integrated into the set of gradient coils 110.

[0075] Adjacent to the imaging zone 108 is a radio-frequency coil 114 for manipulating the orientations of magnetic spins within the imaging zone 108 and for receiving radio transmissions from spins also within the imaging zone 108. The radio frequency antenna may contain multiple coil elements. The radio frequency antenna may also be referred to as a channel or antenna. The radio-frequency coil 114 is connected to a radio frequency transceiver 116. The radio-frequency coil 114 and radio frequency transceiver 116 may be replaced by separate transmit and receive coils and a separate transmitter and receiver. It is understood that the radio-frequency coil 114 and the radio frequency transceiver 116 are representative. The radio-frequency coil 114 is intended to also represent a dedicated transmit antenna and a dedicated receive antenna. Likewise the transceiver 116 may also represent a separate transmitter and receivers. The radio-frequency coil 114 may also have multiple receive/transmit elements and the radio frequency transceiver 116 may have multiple receive/transmit channels. For example if a parallel imaging technique such as SENSE is performed, the radio-frequency could 114 will have multiple coil elements.

[0076] The transceiver 116 and the gradient controller 112 are shown as being connected to a hardware interface 128 of a computer system 126. The computer system further comprises a processor 130 that is in communication with the hardware system 128, a memory 134, and a user interface 132. The memory 134 may be any combination of memory which is accessible to the processor 130. This may include such things as main memory, cached memory, and also non-volatile memory such as flash RAM, hard drives, or other storage devices. In some examples the memory 134 may be considered to be a non-transitory computer-readable medium.

[0077] The memory 134 is shown as containing machine-executable instructions 140. The machine-executable instructions 140 enable the processor 130 to control the operation and function of the magnetic resonance imaging system 100. The machine-executable instructions 140 may also enable the processor 130 to perform various data analysis and calculation functions. The computer memory 134 is further shown as containing pulse sequence commands 142. The pulse sequence commands are configured for controlling the magnetic resonance imaging system 100 to acquire magnetic resonance imaging data.

[0078] The computer memory 134 is further shown as containing the magnetic resonance imaging data 144 that was acquired when the processor 130 executed the pulse sequence commands 142. The memory 134 is also shown as containing current sensor data 146. The current sensor data 146 was acquired at the same time as the magnetic resonance imaging data 144. The two 144, 146 are correlated in time so that it is possible to reconstruct a corrected k-space trajectory 150. The corrected k-space trajectory 150 is shown as being stored in the memory 134. The memory 134 is shown as containing a selected gradient coil transfer function 148. The selected gradient coil transfer function 148 is a mapping between the current sensor data 146 and magnetic field components that are generated within the imaging zone 108. The magnetic field components include the value of the gradient magnetic field generated by the gradient coils 110. The magnetic field components may also contain higher order terms.

[0079] The computer memory 134 is further shown as containing a corrected magnetic resonance image 152 that is reconstructed using the corrected k-space trajectory 150 and the magnetic resonance imaging data 144. In some examples there may be re-sampled magnetic resonance imaging data 144 that is reconstructed as an intermediate step using the magnetic resonance imaging data and the corrected k-space trajectory 150. FIG. 2 shows a flowchart which illustrates a method of operating the magnetic resonance imaging system 100 of FIG. 1. First in step 200 the processor 130 controls the magnetic resonance imaging system 100 with the pulse sequence commands 142. This causes the magnetic resonance imaging system 100 to acquire the magnetic resonance imaging data 144. Step 202 is performed concurrently with step 200. In step 200 the processor 130 receives the current sensor data 146 from the current sensor system 113. Next in step 204 the processor 130 calculates a corrected k-space trajectory using the current sensor data 146 and the selected gradient coil transfer function 148. Finally, in step 206, the processor 130 reconstructs a corrected magnetic resonance image 152 using the magnetic resonance imaging data and the corrected k-space trajectory 150.

[0080] FIG. 3 illustrates a further example of a magnetic resonance imaging system 300. The magnetic resonance imaging system 300 of FIG. 3 is similar to that as is depicted in FIG. 1. The magnetic resonance imaging system 300 has several different additional elements or components within its memory 134. The memory 134 is shown as additionally containing a set of gradient coil transfer functions 302. The set of gradient coil transfer functions 302 is a number of gradient coil transfer functions which can be selected on the basis of the number of acquisition parameters. The selected gradient coil transfer function 148 in this example was selected from the set of gradient coil transfer functions 302. The acquisition parameters 304 are a number of parameters which describe the acquisition or conditions of the acquisition when the magnetic resonance imaging data 144 was acquired.

[0081] Examples of what the acquisition parameters could be could be the subject 118 height, a weight of the subject 118, a cryostat temperature of the magnet 104, a state or operational condition of the cryostat for the magnetic 104, a position of the subject support 120 within the bore of the magnet 106, a temperature of the room which contains the magnet 104, a gradient coil 110 coolant temperature, a gradient coil temperature, a gradient coil impedance, a type of the magnetic resonance imaging protocol used for implementing the pulse sequence commands 142, a type of the receive coil 114 or combinations thereof. Some of the various parameters above may have an effect on the gradient coil transfer function. By using the acquisition parameters to select one it may facilitate the selection of a correct gradient coil transfer function 148. For example, the memory 134 could contain a machine learning algorithm 306 that is used to select the selected gradient coil transfer function 148 from the set of gradient coil transfer functions 302 using the acquisition parameters 304. This could be implemented in any number of ways such as a neural net or deep learning algorithm and also an algorithm which measures the distance such as a nearest neighbors' algorithm could be used for the selection also.

[0082] As an alternative the selected gradient coil transfer function could also be measured before the acquisition of the magnetic resonance imaging data 144. For example, the memory 134 may contain magnetic field measuring pulse sequence commands 310 that contain pulse sequence commands which enable a magnetic resonance imaging protocol that can measure the magnetic field within the imaging zone 108. Execution of these pulse sequence commands 310 may enable the processor 130 to acquire magnetic field magnetic resonance imaging data 312 and the calibration current sensor data 308. From the magnetic field magnetic resonance imaging data 312 a gradient system response 314 may be calculated. The calibration current sensor data 308 and the gradient system response 314 may then be used to calculate the selected gradient coil transfer function 148.

[0083] FIG. 4 shows a flowchart which illustrates a method of operating the magnetic resonance imaging system 300 of FIG. 3. First in step 400 the processor 130 receives one or more acquisition parameters 304. Next in step 402 the processor 130 chooses the selected gradient coil transfer function 148 using the acquisition parameters 304. The selected gradient coil transfer function 148 is selected from the set of gradient coil transfer functions 302. For example, the machine learning algorithm 306 may be used for this. After step 402 the method proceeds to step 200 of the method illustrated in FIG. 2.

[0084] FIG. 5 illustrates a further example of a magnetic resonance imaging system 500. The magnetic resonance imaging system 500 is similar to that depicted in FIG. 1 with the addition of several additional elements or components in its memory 134. The features of the magnetic resonance imaging system which is illustrated in FIGS. 1, 3 and 5 may be freely combined.

[0085] The memory 134 is further shown as containing a set of gradient amplifier transfer functions. A gradient amplifier transfer function 502 is a transfer function which interpolates an input into the gradient amplifier 112 to an output current. The memory 134 is further shown as containing a set of system status parameters 504. The system status parameters are parameters which may affect the operation of the gradient coil amplifier 112. These for example may comprise a room temperature of the room in which the gradient coil amplifier 112 is placed in, it may contain a prior use of the gradient amplifier 112 for prior executions of pulse sequence commands, it may contain a temperature 112 measured within the gradient amplifier itself, a magnetic resonance imaging protocol type used for execution of the pulse sequence commands 142 and combinations thereof. The memory 134 is further shown as containing a selected gradient amplifier transfer function 506 that was selected from the set of selected gradient amplifier transfer functions 502 using the system status parameters 504. This may for example be done using a machine learning algorithm 508. The machine learning algorithm 508 may for example have a neural network, or another algorithm used to determine a closest neighbor based on the system status parameters. The set of gradient amplifier transfer functions may contain entries for the system status parameters 504 for each member of the set of gradient amplifier transfer functions. The machine learning algorithm 508 may then use the system status parameters to select the selected gradient amplifier transfer function 506.

[0086] The memory is further shown as containing a corrected control signal 510 which is used to control the gradient coil amplifier 112 during the execution of the pulse sequence commands 142. The selected gradient amplifier transfer function 506 is essentially a curve which measures the output of the amplifier 112 in relation to the input. Knowing this, the input can be corrected so that the output is more to the desired value. During the acquisition of the magnetic resonance imaging data 144 the current sensor data 146 is measured. The current sensor data 146 can be used to correct the selected gradient amplifier transfer function 506. The memory 134 shows a modified selected gradient amplifier transfer function 512 that was calculated using the current sensor data 146. The modified selected gradient amplifier transfer function 512 can for example be stored in the set of gradient amplifier transfer functions 502 along with the system status parameters 504. This can then be used to further train the machine learning algorithm 508.

[0087] FIG. 6 shows a flowchart which illustrates a method of operating the magnetic resonance imaging system 500 of FIG. 5. First in step 600, the processor 130 receives one or more system status parameters system status parameters 504 descriptive of a status of the magnetic resonance imaging system 500. Next in step 602 the processor 130 chooses the selected gradient amplifier transfer function using the system status parameters 504 from the set of gradient amplifier transfer functions 502. This for example may be one using the machine learning algorithm 504. Then in step 604 the processor 130 calculates a corrected control signal 510 using the control signal from the pulse sequence commands 142 and the selected gradient amplifier transfer function 506. Next the method proceeds to step 200 as is illustrated in FIG. 2.

[0088] The gradient chain (the magnetic field gradient coil amplifier and the set of gradient coils) is an essential part of any MRI system. Its proper function is essential for the correct spatial encoding (of the magnetic resonance imaging data). Since real world hardware is prone to imperfections, any method to measure, characterize or predict deviations of the desired from the actual gradient trajectory (the corrected k-space trajectory) is of major interest. In this way the gradient waveform can either be pre-compensated or the measured data (magnetic resonance imaging data) can be accordingly corrected during reconstruction. Especially non-Cartesian trajectories are sensitive to trajectory deviations and thus would benefit from these improvements.

[0089] Since the deviations can evolve during scanning, e.g. caused by heating of the gradient coil (the set of gradient coils), characterizations of the gradient chain or output waveform “in realtime”, i.e. simultaneously with the imaging may be beneficial. One way to accurately monitor the actual gradient waveform is the implementation of field probes, which quasi-continuously measure the MR-frequency at different locations. Although this is a very accurate method, it requires substantial additions of hardware, which can be very expensive. Furthermore, this additional hardware needs to be interfaced with the existing MR hard- and software, which adds further complexity to the entire MRI system.

[0090] Examples may use electrical sensors (current sensor system) to measure the currents and/or the voltages of the gradient coil and to deduce “on the fly” transfer functions of the gradient system (the gradient coil transfer functions and/or the gradient amplifier transfer functions), respectively their temporal changes. These sensors may already be present in the MR system and are used for the gradient amplifier control circuit, thus no additional hardware has to be implemented. It has turned out, that the accuracy of these measures is close to MR based methods in for characterization of the gradient chain. Additionally, the time signal of the current sensor can be directly used to determine the gradient trajectory applying a known (i.e. measured) transfer function from the current to the actual gradient (the selected gradient coil transfer function).

[0091] The actually applied gradient waveform is of major importance for the image quality, especially for non-Cartesian sampling. If the applied waveform significantly deviates from the desired one, substantial artefacts can occur. These deviations might be caused by:

[0092] Frequency dependent amplitude variations and phase shifts (delays) occurring in the gradient chain (e.g. low pass behavior)

[0093] Variations thereof due to temperature induced changes of the gradient chain/coil

[0094] Non-linearities in elements of the gradient chain, e.g. the amplifier

[0095] Etc.

[0096] Examples may provide for a means of identifying and characterizing one or more of the above effects. Using the measured output current of the gradient amplifier as a basis for trajectory calculation, all non-idealities occurring earlier in the gradient chain can be reduced or removed. In addition, examples may allow in general for a comparison of the measured gradient chain characteristics with the expected performance, thus enabling to identify potential hardware issues or defects (e.g. via the frequency dependent

[0097] impedance of the gradient coil).

[0098] A time dependent eddy current behavior may not necessarily directly measured by examples, however, this can be linked to the applied gradient currents and their history. Since these can be traced over time by examples, it is possible to apply predefined models that include also eddy current related effects.

[0099] Examples may use current and voltage sensors as well as the known gradient waveform input signal of the gradient amplifier as first order “field camera”. Based on this one-dimensional information (Gx, Gy, Gz):

[0100] The actual gradient system output and the resulting k-space trajectory can be estimated almost in real-time in a sense of a low budget field camera;

[0101] Furthermore, this information can be used to determine corresponding transfer functions. These can be updated regularly and thus changes of the gradient chain can be detected. (Defining the transfer function in the frequency domain allows to characterize the gradient chain and possible changes thereof spectrally, thus providing additional information.)

[0102] The ability to continuously monitor the gradient system under different load/duty cycle/working conditions allows to derive system specific generalized transfer functions using deep learning algorithms. This permits also to predict individual system performance learned from actual history to better steer reconstruction. The transfer functions, which can be regularly computed, are e.g. those relating the output currents to the desired input currents and those relating the actual output voltage to the output current (=the impedance). The continuous measurement of these functions allows to identify temperature drifts, changing delays, nonlinearities and potential defects.

[0103] In one example the time series of the current sensor is directly used to predict the actual gradient by applying a known transfer function between current and gradient (e.g. measured by an MR measurement)

[0104] Some examples use sensors at least for the gradient amplifier input signal, the output current and the output voltage for each axis of the gradient coils. These sensors may be connected to an ADC in order to allow for subsequent processing of the measured data. The current sensor may be of high accuracy, such sensors (e.g. current transducers) may be already used in some gradient amplifiers for the current control loop. Consequently, these sensors can also be used to implement examples and the additional hardware effort is minimal. Furthermore, depending of the hardware configuration the “demand input” information could already be present as digital information making explicit sampling obsolete.

[0105] The same holds for voltage and input current sensors. An overview of the proposed setup is shown in FIG. 7. FIG. 7 illustrates an example of a current sensor system 710 which is an electrical sensor to measure and record input signals and output 702 of the gradient coil amplifier 112. In this example the input signal is an analogue signal which is then measured as the demand current 704 by an analogue-to-digital converter. In other examples the input into the gradient amplifier 112 could itself be a digital signal in which case the demand current 704 is already known. The current sensor system 710 may include a current 704 and optionally a voltage 706 measurement. Both the current 708 and the voltage 706 measurements could be measured by an analogue-to-digital converter. References herein to current sensor data may also include the digitized voltage 706 measurements.

[0106] FIG. 8 shows a block diagram which illustrates the relationship of the selected gradient coil transfer function 148 to the selected gradient amplifier transfer function 506. The gradient amplifier 112 is shown as supplying current to the gradient coils 110. The transfer function for the gradient coil amplifier 112 relays a demand current or signal to resulting current. The current supplied to the gradient coil is then modeled using the selected gradient coil transfer function 148. This then results in the magnetic resonance imaging data 144 being acquired. The current sensor data can then be used to calculate a corrected k-space trajectory. FIG. 7 illustrates how to predicting the actual gradient trajectory from the measured currents applying a previously measured transfer function such as is shown in FIG. 8.

[0107] FIG. 9 illustrates an example of a gradient coil transfer function 148. This for example may relate the current supplied to one of the sets of gradient coils to the field within a particular voxel or location within the imaging zone.

[0108] The processing of the recorded (or available) signals performed in different ways. In a simple approach, the measured current is used in conjunction with a previously measured transfer (current to field gradient) to predict the actually applied gradient trajectory (i.e. based on the actually applied current). This approach is depicted in FIG. 8. The required transfer function can e.g. be deduced from a calibration measurement applying chirped or other gradients to determine a gradient impulse response function (GIRF). The calculation of the gradient trajectory can either be performed in frequency or time domain (convolution).

[0109] The determined gradient trajectory can finally be used for calculation of a highly accurate k-space trajectory during reconstruction.

[0110] In another example, transfer functions (gradient amplifier transfer function) of the sensor signals in frequency domain are continuously calculated, as shown in FIG. 10. FIG. 10 shows a block diagram which illustrates how the input to the gradient coil amplifier 702 and the output of the gradient coil amplifier 702 can be used to calculate a gradient amplifier transfer function 512.

[0111] The transfer functions of interest are those relating input signal to output

[0112] current and those relating output voltage to output current (=impedance). Knowing these functions and also their evolution over time, the gradient trajectory can be deduced (as described above) but also temporal changes of the gradient chain can be identified. As an example a transfer function relating input signal to output current is shown in FIGS. 11 and 12. FIGS. 11 and 12 show sketches of the real 1100 and the phase 1200 of a gradient amplifier transfer function. The units in these two figures are arbitrary. These functions shown in FIGS. 11 and 12 and their temporal behavior provide substantial information on the performance of the gradient chain.

[0113] Similarly, the frequency dependent impedance provides important information on the gradient coil, the resistance can serve as a measure of gradient coil temperature and the resonance-like features in the inductance are due to mechanical resonances of the coil, which can provide insights in the proper function of the coil.

[0114] The frequency dependent impedance of the gradient coil as it can be determined using a transfer function approach. Heating can be observed as a gradual increase of the resistance near DC in the real part of the impedance. The inductance may show resonance-like structures at several frequencies, which are due to mechanical resonances. These features provide valuable information on the proper function of the gradient coil.

[0115] In a further application example, the sensor outputs are used, employing appropriate transfer function postprocessing to continuously monitor the gradient system under different load/duty-cycle/working conditions. Deep learning algorithms are employed to distil the key features of the transfer functions depending on the different system working conditions. This will result in system specific generalized transfer functions permitting to predict individual system performance learned from actual history to better steer reconstruction.

[0116] While the invention has been illustrated and described in detail in the drawings and foregoing description, such illustration and description are to be considered illustrative or exemplary and not restrictive; the invention is not limited to the disclosed embodiments.

[0117] Other variations to the disclosed embodiments can be understood and effected by those skilled in the art in practicing the claimed invention, from a study of the drawings, the disclosure, and the appended claims. In the claims, the word “comprising” does not exclude other elements or steps, and the indefinite article “a” or “an” does not exclude a plurality. A single processor or other unit may fulfill the functions of several items recited in the claims. The mere fact that certain measures are recited in mutually different dependent claims does not indicate that a combination of these measured cannot be used to advantage. A computer program may be stored/distributed on a suitable medium, such as an optical storage medium or a solid-state medium supplied together with or as part of other hardware, but may also be distributed in other forms, such as via the Internet or other wired or wireless telecommunication systems. Any reference signs in the claims should not be construed as limiting the scope.

LIST OF REFERENCE NUMERALS

[0118] 100 magnetic resonance imaging system [0119] 104 magnet [0120] 106 bore of magnet [0121] 108 imaging zone [0122] 109 region of interest [0123] 110 set of gradient coils [0124] 112 magnetic field gradient coil amplifier [0125] 113 current sensor system [0126] 114 radio-frequency coil [0127] 116 transceiver [0128] 118 subject [0129] 120 subject support [0130] 122 x-axis [0131] 123 y-axis [0132] 124 z-axis [0133] 126 computer system [0134] 128 hardware interface [0135] 130 processor [0136] 132 user interface [0137] 134 computer memory [0138] 140 machine executable instructions [0139] 142 pulse sequence commands [0140] 144 magnetic resonance imaging data [0141] 146 current sensor data [0142] 148 selected gradient coil transfer function [0143] 150 corrected k-space trajectory [0144] 152 corrected magnetic resonance image [0145] 200 control the magnetic resonance imaging system with the pulse sequence commands to acquire the magnetic resonance imaging data [0146] 202 record the current sensor data during the acquisition of the magnetic resonance imaging data [0147] 204 calculate a corrected k-space trajectory using the current sensor data and the selected gradient coil transfer function [0148] 206 reconstruct a corrected magnetic resonance image using the magnetic resonance imaging data and the corrected k-space trajectory according to the magnetic resonance imaging protocol [0149] 300 magnetic resonance imaging system [0150] 302 set of gradient coil transfer functions [0151] 304 acquisition parameters [0152] 306 machine learning algorithm [0153] 308 calibration current sensor data [0154] 310 magnetic field measuring pulse sequence commands [0155] 312 magnetic field magnetic resonance data [0156] 314 gradient system response [0157] 400 receive one or more acquisition parameters descriptive of the magnetic resonance imaging protocol [0158] 402 choose the selected gradient coil transfer function using the acquisition parameters [0159] 500 magnetic resonance imaging system [0160] 502 set of gradient amplifier transfer functions [0161] 504 system status parameters [0162] 506 selected gradient amplifier transfer function [0163] 508 machine learning algorithm [0164] 510 corrected control signal [0165] 512 modified selected gradient amplifier transfer function [0166] 600 receive one or more system status parameters descriptive of a status of the magnetic resonance imaging system [0167] 602 choose the selected gradient amplifier transfer function from the set of gradient amplifier transfer functions using the one or more system status parameters [0168] 604 calculate a corrected control signal using the control signal and the selected gradient amplifier transfer function, wherein the gradient amplifier is controlled with the corrected control signal during acquisition of the magnetic resonance imaging data [0169] 700 input to gradient coil amplifier [0170] 702 output of gradient coil amplifier [0171] 704 demand current [0172] 706 voltage [0173] 708 current [0174] 710 current sensor system [0175] 1100 sketch of magnitude of gradient amplifier transfer function