Embodiments described herein relate generally to a magnetic resonance imaging (MRI) apparatus.
An MRI apparatus is an imaging apparatus configured to excite nuclear spin of an object placed in a static magnetic field by a radio frequency (RF) pulse applied with the Larmor frequency, and reconstruct an image based on magnetic resonance (MR) signals emitted from the object due to the excitation.
Among the imaging methods using an MRI apparatus, there is known a non-contrast enhanced imaging for imaging blood flow or blood vessels without administering a contrast medium. Further, among the various non-contrast enhanced imaging methods, a Time-SLIP (time-spatial labeling inversion pulse) method is known as an imaging method in which an RF pulse (for example, an inversion pulse) is applied prior to data acquisition to control depiction of blood flow. The Time-SLIP method is an imaging method that can clearly depict or visualize the blood flow of a target blood vessel with high contrast, while suppressing the background of the imaging region. By using the Time-SLIP method, not only blood flowing through arteries or veins but also fluid such as cerebrospinal fluid (CSF) can be imaged.
Meanwhile, ECG (electrocardiogram) gating and respiratory gating are used to suppress body-motion artifacts caused by body motions due to heartbeat and/or respiration. The ECG gating is an imaging method in which data (i.e., MR signals) are acquired within a predetermined period of a desired cardiac time phase, in synchronization with an ECG signal, for example, an R-wave. On the other hand, the respiratory gating is an imaging method in which data are acquired within a predetermined period of a predetermined time phase during expiration in a cycle consisting of inspiration and expiration, in synchronization with the respiratory signal indicating the respiratory motion detected by a respiratory sensor, for example.
In a conventional MRI apparatus, the Time-SLIP method has been used in combination with the ECG gating and/or respiratory gating to suppress body motion artifacts. When the Time-SLIP method and gating imaging, such as ECG gating and respiratory gating, are used in combination, the total imaging time tend to be long.
In the accompanying drawings:
Hereinafter, an MRI apparatus 1 of the present embodiment will be described by referring to the accompanying drawings.
In one embodiment, an MRI apparatus includes: a scanner that includes a static magnetic field magnet configured to generate a static magnetic field, a gradient coil configured to generate a gradient magnetic field, and a WB (Whole Body) coil configured to apply an RF pulse to an object; and processing circuitry. The processing circuitry is configured to: set (i) a pulse sequence in which a sequence element is repeated, the sequence element including at least an inversion pulse and (ii) a data acquisition sequence executed after a delay time from the inversion pulse; and to cause the scanner to execute the pulse sequence by using virtual gating.
(Configuration and Basic Operation)
The gantry 100 includes, for example, a static magnetic field magnet 10, a gradient coil 11, a WB (whole body) coil 12, and these components are provided in a cylindrical housing. The bed 500 includes a bed body 50 and a table 51. Additionally, The MRI apparatus 1 may include an array coil 20 that is attached near an object. Further, a respiratory sensor 700 and/or an electrocardiograph 710 for respiratory gating and/or ECG (electrocardiogram) gating may be attached to the object, as described below.
The control cabinet 300 includes three gradient coil power supplies 31 (specifically, 31x for the X-axis, 31y for the Y-axis, and 31z for the Z-axis), an RF receiver 32, an RF transmitter 33, and a sequence controller 34.
The static magnetic field magnet 10 of the gantry 100 is substantially in the form of a cylinder, and generates a static magnetic field inside the bore (i.e., the space inside the cylindrical structure of the static magnetic field magnet 10) which is an imaging region of the object (for example, a patient). The static magnetic field magnet 10 includes a built-in superconducting coil, and the superconducting coil is cooled down to an extremely low temperature by liquid helium. The static magnetic field magnet 10 generates a static magnetic field by supplying the superconducting coil with electric current provided from a static magnetic field power supply (not shown) in an excitation mode. Afterward, the static magnetic field magnet 10 shifts to a permanent current mode, then the static magnetic field power supply is separated. Once it enters the permanent current mode, the static magnetic field magnet 10 continues to generate a strong static magnetic field for a long time, for example, over one year. Note that the static magnetic field magnet 10 may be configured as a permanent magnet.
The gradient coil 11 is also substantially in the form of a cylinder, and is fixed to the inside of the static magnetic field magnet 10. This gradient coil 11 applies gradient magnetic fields to the object in the respective directions of the X-axis, the Y-axis, and the Z-axis of the apparatus coordinate system, by using electric currents supplied from the gradient coil power supplies 31x, 31y, and 31z.
The bed body 50 of the bed 500 moves the table 51 in the upward and downward directions, and moves the table 51 with the object laid thereon to a predetermined height before imaging. Afterward, during the imaging, the bed body 50 moves the table 51 in the horizontal direction so as to move the object inside the bore.
The WB coil 12 is also shaped substantially in the form of a cylinder so as to surround the object, and is fixed to the inside of the gradient coil 11. The WB coil 12 applies RF pulses transmitted from the RF transmitter 33 to the object, and receives MR signals emitted from the object due to excitation of hydrogen nuclei.
The array coil 20 is a kind of RF coil, which is place adjacent to the object, and receives MR signals emitted from the object. The array coil 20 is, for example, consisted of a plurality of coil elements. Although there are various types of the array coil 20 such as a head coil, a chest coil, a spine coil, a lower-limb coil, and a whole-body coil depending on an anatomical imaging site of the object, the array coil 20 for the chest is illustrated in
The RF transmitter 33 sends RF pulses to the WB coil 12 on the basis of commands inputted from the sequence controller 34. The RF receiver 32 receives MR signals through the WB coil 12 and/or the array coil 20, and sends raw data obtained by digitizing the received MR signals to the sequence controller 34.
The sequence controller 34 performs a scan of the object by driving the gradient coil power supplies 31, the RF transmitter 33, and the RF receiver 32, under the control of the console 400. When the sequence controller 34 receives the raw data from the RF receiver 32 by performing a scan, the sequence controller 34 sends the received raw data to the console 400.
The sequence controller 34 includes processing circuitry (not shown), which is configured as hardware such as a processor for executing predetermined programs, a field programmable gate array (FPGA), and an application specific integrated circuit (ASIC).
The console 400 is configured as a computer including processing circuitry 40, a memory 41, a display 42, and an input device 43.
The memory 41 is a recording medium including a read-only memory (ROM) and a random access memory (RAM), in addition to an external memory device such as a hard disk drive (HDD) and/or an optical disc device. The memory 41 stores various programs executed by a processor of the processing circuitry 40 in addition to various data and information.
The input device 43 includes various types of devices such as a mouse, a keyboard, a trackball, and/or a touch panel, for an operator to input various data and information.
The display 42 is a display device such as a liquid crystal display panel, a plasma display panel, and an organic EL panel.
The processing circuitry 40 is, for example, a circuit provided with a CPU and/or a special-purpose or general-purpose processor. The processor realizes various functions described below by executing programs stored in the memory 41. The processing circuitry 40 may be consisted of hardware such as an FPGA and an ASIC. The various functions described below can also be realized by such hardware. Additionally, the processing circuitry 40 can achieve the various functions by combining hardware processing and software processing based on its processor and programs.
With the above components, the console 400 performs system control of the entirety of the MRI apparatus 1. Specifically, the console 400 receives various commands and information such as imaging conditions, which are inputted by an operator such as a medical imaging technologist through the mouse and/or the keyboard of the input device 43. The processing circuitry 40 causes the sequence controller 34 to perform a scan on the basis of the inputted imaging conditions, and reconstructs images by using the raw data transmitted from the sequence controller 34. The reconstructed images may be displayed on the display 42, and may be stored in the memory 41.
Further, the MRI apparatus 1 of the present embodiment generates an image of the object by reconstructing the MR signals acquired by executing a scan under pulse sequences described below. In particular, the MRI apparatus 1 of the present embodiment performs a scan using a pulse sequence that can decrease the imaging time as compared with the conventional case.
Among the components of the MRI apparatus 1 shown in
As shown in
The user-interface control function F01 performs control for users, such as doctors and/or medical imaging technologists, to input and select various data related to imaging via the user interface 401. The user interface 401 includes, for example, the input device 43 and the display 42 described above.
The user inputs or selects data regarding the imaging conditions of the imaging to be executed via the user interface 401. Here, the imaging conditions may include, for example, a type of imaging method, a type of imaging target, FOV (field of view), resolution, and so forth. The user may input or select a specific value of each of various parameters of the pulse sequence for realizing these imaging conditions.
Types of imaging methods include, for example, sampling methods in a k-space, such as two-dimensional Cartesian sampling, three-dimensional Cartesian sampling, two-dimensional radial sampling, and three-dimensional radial sampling.
In addition, the types of imaging methods may include types of pulse sequence, for example, a FLASH (fast low angel shot) method, an SPGR (spoiled gradient-recalled acquisition in the steady state) method, a balanced-SSFP (steady-state free precession) method, an FSE (fast spin echo) method, a FASE (fast advanced spin echo) method, an EPI (echo planar imaging) method, a TOF (time of flight) method, a PC (phase contrast) method, an FBI (fresh blood imaging) method, a Time-SLIP (time-spatial labelling inversion pulse) method, a DWI (diffusion weighted Imaging) method, an ASL (arterial spin labeling) method, and so forth.
The types of imaging targets selected via the user interface 401 include, for example, (i) a site of the body of the object such as the head, chest, abdomen, and inferior limbs, (ii) the type of organs such as the brain, heart, lung, liver, and kidney, and (iii) the type of tissues such as blood, cerebrospinal fluid, and bones.
Meanwhile, a standard imaging protocol and/or a recommended imaging protocol, which correspond to imaging conditions such as the type of imaging method and the type of imaging target, are stored in an imaging protocol storage region 411 of the memory 41. The imaging protocol may include, for example, the type of pulse sequence used for imaging, values of various parameters in the pulse sequence, combinations of a plurality of pulse sequences in the case of consecutively executing the plurality of pulse sequences. The values of various parameters in the pulse sequence may include, for example, values of respective parameters related to the transmission RF pulse, such as a flip angle α of the excitation pulse, a repetition period TR of the excitation pulse and their temporal change, and values of respective parameters related to the gradient magnetic field, such as the magnitude and application time of each gradient magnetic field in the X-axis, Y-axis, and Z-axis directions and their temporal changes.
The user-interface control function F01 reads out, from the imaging-protocol storage region 411 of the memory 41, the standard imaging protocol corresponding to the imaging conditions inputted or selected by the user.
The user-interface control function F01 may display the imaging protocol read out on the display 42 of the user interface 401. The user can approve and determine the displayed imaging protocol as they are. Alternatively, if necessary, the user may change the displayed imaging protocol, and then, may determine the changed imaging protocol.
The user-interface control function F01 sends the determined imaging protocol to the imaging sequence setting function F02. The imaging sequence setting function F02 sets the data related to the pulse sequence corresponding to the received imaging protocol on the sequence controller 34 of the scanner 600, as an imaging sequence.
The sequence controller 34 executes the pulse sequence according to the data related to the set pulse sequence as described above so as to acquire the MR signals.
The acquired MR signals are sent to the reconstruction function F05 of the processing circuitry 40. The reconstruction function F05 performs reconstruction processing such as Fourier transform on the MR signals so as to generate a magnetic resonance image. The generated magnetic resonance image is displayed on the display 42, for example.
The gating selection function F03 and the virtual gating period calculation function F04 are functions related to virtual gating, which is a characteristic of the MRI apparatus 1 of the present embodiment. Hereinafter, the virtual gating will be described in detail.
(Virtual Gating)
The pulse sequence targeted by the virtual gating according to the present embodiment is mainly a pulse sequence in which a sequence element is repeated, and this sequence element includes at least an inversion pulse and a data acquisition sequence executed after a predetermined delay time from the inversion pulse.
A typical example of this type of pulse sequence is the Time-SLIP. As an example of the above-described pulse sequences, the MRI apparatus 1 of the present embodiment executes a pulse sequence by applying the Time-SLIP method using virtual gating.
Prior to description of virtual gating, the Time-SLIP method will be briefly described by referring to
On the other hand,
Of the two inversion pulses, a non-selective inversion pulse A may be applied first, and a selective inversion pulse B may be applied immediately after this pulse A, for example. These two pulses, i.e., the non-selective inversion pulse A and the selective inversion pulse B may be collectively referred to as a Time-SLIP pulse in the following description.
Note that the order of application of the non-selective inversion pulse A and the selective inversion pulse B may be reversed.
The delay time from the Time-SLIP pulse to the front edge of the data acquisition sequence is called BBTI (“black blood inversion time” or “black blood traveling time”).
In the data acquisition sequence, MR signals are acquired by a plurality of excitation pulses. Various types of sequences can be used in the data acquisition sequence as described below. When the MR signals of the entire k-space, which are required for reconstructing a full image, cannot be acquired by one data acquisition sequence, the MR signals are acquired by dividing the entire k-space into a plurality of regions (i.e., into a plurality of segments). In this case, the sequence diagram shown in
The unit of repetition, i.e., the period from the front end of the Time-SLIP pulse to the front end of the next Time-SLIP pulse, is called a “sequence element” in the following description. In this case, the pulse sequence of the Time-SLIP method is composed of a plurality of sequence elements that are repeated over time by the number of segments.
The non-selective inversion pulse A is applied to the entire region (i.e., region of the square frame shown by the fine broken line in
Meanwhile, the selective inversion pulse B is applied to a desired selected region outside the imaging region. Specifically, the selective inversion pulse B is applied to a selected region including a specific blood vessel that flows into the imaging region from the outside of the imaging region on the upstream side of the imaging region. In the case of
Meanwhile, data acquisition is performed around the timing when the longitudinal magnetization of the background in the imaging region crosses the null point (i.e., around the time when BBTI elapses from the Time-SLIP pulse), and thereby, the blood of the artery α flowing from the outside of the imaging region can be depicted with high contrast with respect to the background of the imaging region. Although the blood of the artery β also flows in the imaging region, the selective inversion pulse B is not applied to the artery β, and thus, the longitudinal magnetization of the blood of the artery β flowing in the imaging region is maintained at a small value. In this way, in the Flow-Out type Time-SLIP methods, among a plurality of blood vessels in the imaging region, blood of the specified blood vessel only can be selectively depicted.
The selective inversion pulse B is not applied to the artery α, which is located outside the selected region and the imaging region. Thus, the blood in the artery α flows into the imaging region while maintaining a large longitudinal magnetization value. Meanwhile, data acquisition is performed around the timing when the longitudinal magnetization of the background in the imaging region crosses the null point (i.e., around the time when BBTI elapses from the Time-SLIP pulse), and thereby, the blood of the artery α flowing into the imaging region from the outside can be depicted with high contrast with respect to the background of the imaging region and with respect to the blood flowing into the imaging region from the vein γ in the selected region.
Data acquisition is usually done around the center of the expiration period during the respiratory cycle. If the BBTI is set to be longer than the respiratory cycle, the data acquisition may be done every two cycles of the respiratory cycle, as shown in
In the respiratory gating, the Time-SLIP pulse starts in synchronization with the respiratory-gating signal. In the case of acquiring data by dividing k-space into a plurality of segments, the sequence element specified by the period from the Time-SLIP pulse to the next Time-SLIP pulse is repeated in synchronization with the respiratory-gating signal. In this way, it is possible to acquire k-space data at the same time phase of the respiratory cycle by the data acquisition sequence of every sequence element.
Note that, in the conventional Time-SLIP method using the respiratory gating, a dead time D exists between the trailing edge of the data acquisition sequence and the leading edge of the next Time-SLIP pulse, as can be seen in
The dead time D changes depending on the length of BBTI. The shorter the BBTI becomes, the longer the dead time D becomes. As a result, the overall imaging time becomes longer.
The main aim of the present embodiment is to reduce this dead time D and decrease the overall imaging time.
The main purpose of the above-descried actual gating, such as the respiratory gating and the ECG gating, is to suppress artifacts caused by body motions such as respiratory motions and pulsations. As shown in
In the respiratory gating, data acquisition is performed in the same respiratory phase between segments, resulting in suppressing fluctuations of k-space data between segments due to body motions. Similarly, in the ECG gating, data acquisition is performed in the same cardiac time phase between segments, resulting in suppressing fluctuations of k-space data between segments due to body motions as well.
Meanwhile, it should be noted that body-motion artifacts depend on the sampling method of k-space. Specifically, as described below, it is known that body-motion artifacts are likely to occur when data acquisition is performed by using Cartesian sampling, while the body-motion artifacts are unlikely to occur when data acquisition is performed by using radial sampling.
In addition, it is also known that the degree of occurrence of the body motion artifacts varies depending on the imaging site (i.e., the target site of the object to be imaged). When the imaging site of the object is a chest or an abdomen, body-motion artifacts are likely to occur because the body motions due to respiration and heartbeat are large. Contrastively, when the imaging site is a part away from the chest and abdomen such as a head and inferior limb, body-motion artifacts are less likely to occur because such a part is less susceptible to the body motions due to respiration or heartbeat.
From the above viewpoint, the MRI apparatus 1 of the present embodiment is configured to select either of actual gating (such as respiratory gating and ECG gating) or virtual gating (without using actual gating) for executing a pulse sequence, depending on the sampling method of k-space data and/or the imaging site of the object (or, the imaging site of the patient).
In the step ST100 of
In the next step ST101, the sampling type of the data acquisition sequence in the pulse sequence of the Time-SLIP method is determined based on the imaging conditions set in the step ST100. Specifically, if the data acquisition sequence is performed by radial sampling, the processing proceeds to the step ST102, in which a pulse sequence using virtual gating is executed. On the other hand, if the data acquisition sequence is performed by Cartesian sampling instead of radial sampling, the processing proceeds to the step ST103.
As examples of Cartesian scans,
The 2D radial sampling (2) shown in
In the radial sampling of the 2D UTE method, the MR signals are sampled, for each excitation pulse, along a line extending from the center of k-space to the end of k-space (i.e., the position indicated by the tip of the arrow in
The sequence diagram of the three-dimensional UTE method using the Koosh ball type radial sampling can be obtained by adding the gradient magnetic field Gz in the Z-axis direction, which changes in the same manner as Gx and Gy, to the sequence diagram of the two-dimensional UTE shown in
The 3D radial sampling (2) shown in
Returning to
If the set imaging conditions or the selected type of data acquisition sequence uses the Cartesian sampling, not the radial sampling described above, the processing proceeds to the step ST103.
In the step ST103, based on the set imaging conditions, it is determined whether the imaging target site is a site with small body motions (for example, a head or inferior limbs) or a site with large body motions (for example, a chest or abdomen).
If the imaging target site is the head or inferior limbs, body motions due to respiration and heartbeat are small, and thus, even when using Cartesian sampling for the data acquisition sequence, body-motion artifacts are unlikely to occur. Hence, in this case, the processing proceeds to the step ST102, in which a pulse sequence using virtual gating, specifically, a pulse sequence based on the Time-SLIP method is executed. That is, when the imaging target site is the head or the inferior limbs, virtual gating can be applied to all the sampling methods shown in
If the data acquisition sequence is executed by using Cartesian sampling “and” the imaging target site is a site with large body motions due to respiration and/or heartbeat, such as the abdomen or chest, in order to suppress body-motion artifacts, the MRI apparatus 1 executes a pulse sequence based on the Time-SLIP method using actual gating such as respiratory gating and ECG gating, which are shown in
The processing of the steps ST101 and ST103 in
Next, processing related to execution of the Time-SLIP pulse sequence using virtual gating will be described in detail by referring to the flowchart of
As described above, the sequence element is defined by the behavior of the RF pulse and the gradient magnetic field during the period from the leading edge of the Time-SLIP pulse to the leading edge of the next Time-SLIP pulse. The sequence element includes at least a Time-SLIP pulse and a data acquisition sequence that is executed after a delay time BBTI from the Time-SLIP pulse.
As mentioned above, the Time-SLIP pulse may include both of the non-selective inversion pulse A and the selective inversion pulse B (Flow-Out type), or may include only the selective inversion pulse B (Flow-In type).
In execution of a pulse sequence using virtual gating, it is important to calculate the repetition period TIR-IR of virtual gating. The repetition period TIR-IR of virtual gating is the repetition period TIR-IR of the inversion pulse (i.e., Time-SLIP pulse), and is also the repetition period TIR-IR of the sequence element.
As can be seen from
The steps ST200 to ST205 in
In the step ST200, BBTI, which is set by a user via the user interface 401, is acquired. For example, for the Flow-Out type Time-SLIP method, the value of BBTI is set by a user in consideration of distance between the imaging region and the selected region (i.e., application region of the selective inversion pulse B) set on the upstream side of the imaging region, velocity of blood flowing from the selected region into the imaging region, and recovery time of the longitudinal magnetization of the background in the imaging region.
In the step ST201, sequence parameters are acquired. Specifically, in the step ST201, the values of various parameters in the Time-SLIP pulse sequence, which are set by the user-interface control function F01 and/or the imaging sequence setting function F02, are acquired.
More specifically, the values of various parameters such as (a) the duration DIR of the inversion pulse, i.e., duration DIR of the Time-SLIP pulse, (b) number NL of total lines in k-space data, (c) number Ns of segments when k-space is divided into a plurality of segments to acquire k-space data, (d) repetition period TR of the excitation pulse used in the data acquisition sequence, and (e) flip angle α of the excitation pulse, are acquired as sequence parameters.
In the step ST202, tissue parameters of the imaging target are acquired. The tissue parameters are, for example, transverse relaxation time T2 and longitudinal relaxation time T1 of the imaging target. The values of these tissue parameters can be estimated to some extent from, for example, the imaging target site of the body and/or other imaging conditions. For example, when the imaging target site is the lung field, it can be estimated that the transverse relaxation time T2 shows a very short value (for example, 1 ms or less). Alternatively, when the imaging target site is the head and the imaging target tissue is CSF, it can be estimated that the transverse relaxation time T2 shows a very long value (for example, about 1000 ms).
The flip angle of the excitation pulse acquired in the step ST201, the transverse relaxation time T2 and the longitudinal relaxation time T1 acquired (or estimated) in the step ST202 are used for determining the waiting time WT.
In the next step ST203, the period Tack of the data acquisition sequence is calculated. As can be seen from the sequence diagram of
Tack=N*TR=(NL/Ns)*TR Expression 1
In the next step ST204, the waiting time WT is determined. The waiting time WT is the time from the trailing edge of the data acquisition sequence to the leading edge of the next Time-SLIP pulse. As described above, the waiting time WT can be considered as the time during which the longitudinal magnetization tilted by the excitation pulse train of the data acquisition sequence is sufficiently recovered, or the time during which the transverse magnetization generated by the excitation pulse train of the data acquisition sequence sufficiently disappears. If the waiting time WT can be determined to be short, the overall length of the Time-SLIP pulse sequence can be decreased, and thus, the imaging time based on the Time-SLIP method can be decreased.
In the step S1, it is determined whether the transverse relaxation time T2 of the tissue of the imaging target is very short. If so, for example T2<1 ms, the waiting time WT is determined to be the predetermined minimum value WTmin, for example 5 ms.
In the next step S2, it is determined whether, by using the Time-SLIP method, blood flowing in from the outside imaging region is well depicted with the background sufficiently suppressed. If so, the contrast between blood and background can be considered to be sufficiently high. In this case as well, the waiting time WT is set to be the predetermined minimum value WTmin, for example 5 ms.
In the next step S3, it is determined whether the flip angle is very small, even when the transverse relaxation time T2 of the imaging target tissue is moderate. When the flip angle is small, for example, when the flip angle is 5° as shown in
Meanwhile, as shown in the step S4, when the transverse relaxation time T2 of the imaging target tissue is relatively long, for example, when the imaging target tissue is the synovial fluid or the CSF, the transverse relaxation time T2 shows a relatively long value such as T2=200 to 1000 ms. In such a case, the waiting time WT also must be determined to be a relatively long value. For example, the waiting time WT is about twice that of T2, and thus, in this case, the waiting time WT is determined to be a value of about 400 to 2000 ms.
Returning to
TIR-IR=DIR+BBTI+Tack+WT Expression 2
Finally, in the step ST206, virtual gating is used to perform a Time-SLIP pulse sequence as shown in
As to execution of the Time-SLIP pulse sequence, a virtual gating signal, which is shown in
Alternatively, the sequence element of the Time-SLIP sequence shown in
In the Time-SLIP pulse sequence using the conventional respiratory gating, there exists a dead time D between the trailing edge of the data acquisition sequence and the leading edge of the next Time-SLIP pulse, because the data acquisition sequence is executed at the same time phase of the respiratory cycle in order to suppress the body-motion artifacts caused by respiratory motions. The length of the dead time D depends on the length of the respiratory cycle and BBTI, and can be non-negligibly long. Thus, the presence of the dead time D is the major factor for lengthening the imaging time of the conventional Time-SLIP method using the respiratory gating.
In the Time-SLIP pulse sequence using virtual gating according to the present embodiment, the waiting time WT is used instead of the dead time D. The waiting time WT is a value determined without depending on respiration and/or heartbeat, and is determined in consideration of parameters such as the T2 value of the imaging target tissue and/or the flip angle of the excitation pulse of the data acquisition sequence. The waiting time WT is determined to be a very short value, for example, about 5 to 10 ms, except in the case for imaging a tissue having a very long transverse relaxation time T2 such as cerebrospinal fluid and synovial fluid.
Hence, the imaging time of the Time-SLIP pulse sequence using the virtual gating according to the present embodiment is significantly decreased as compared with the imaging time using the conventional respiratory gating.
Duration DIR of Time-SLIP pulse: 100 ms
Number NL of total lines in k-space: 16000 lines
Number Ns of segments: 64
Number (NL/Ns) of lines per segment: 250 lines
Repetition period TR of excitation pulse: 4 ms
In this case, the data acquisition time Tack per segment is 1000 ms (4 ms*250 lines).
Further assuming that BBTI is set to 1000 ms, the total time Ta (=DIR+BBTI+Tack) from the leading edge of the Time-SLIP pulse to the trailing edge of the data acquisition period is 2100 ms.
Assuming that the repetition period Tresp of respiratory gating (i.e., period Tresp of respiratory-cycle) is 5 seconds (5000 ms), the dead time D in the Time-SLIP pulse sequence using the conventional respiratory gating is, Tresp−Ta=5000−2100=2900 ms.
The total imaging time Ttotal of the Time-SLIP pulse sequence using the conventional respiratory gating is 320 seconds (5 minutes and 20 seconds) calculated as follows.
Ttotal=Tresp*Number of segments=5 seconds*64=320 seconds
By contrast, assuming a case where, by using the Time-SLIP method, blood flowing in from the outside imaging region is well depicted with the background sufficiently suppressed, the waiting time WT can be set to the predetermined minimum value WTmin (=5 ms).
Thus, the repetition period TIR-IR of virtual gating in the Time-SLIP pulse sequence using the virtual gating of the present embodiment (i.e., repetition period TIR-IR of the Time-SLIP pulse) is TIR-IR=Ta+WT=2100+5=2105 ms.
As a result, when using virtual gating of the present embodiment, the total imaging time Ttotal is 135 seconds (i.e., 2 minutes and 15 seconds) calculated as follows.
Ttotal=TIR-IR*number of segments=2105 ms*64=135 seconds
As described above, according to the above numerical example, the total imaging time is 320 seconds (5 minutes and 20 seconds) in the conventional respiratory gating, whereas the total imaging time in virtual gating of the present embodiment is 135 seconds (2 minutes and 15 seconds), which shows that the present embodiment can significantly reduce the imaging time.
In addition, it should be noted that, in the conventional respiratory gating, there is a possibility that the magnitude of the acquired MR signals fluctuates due to fluctuations in the length of the respiratory cycle. Contrastively, in the virtual gating of the present embodiment, the length of the virtual gating cycle is always constant. Thus, in the virtual gating of the present embodiment, the magnitude of the MR signals does not fluctuate due to the fluctuation in the length of the acquisition cycle, and occurrence of artifact due to the amplitude fluctuation can be suppressed.
As described above, the MRI apparatus of the present embodiment can decrease the imaging time of non-contrast enhanced imaging.
While certain embodiments have been described, these embodiments have been presented by way of example only, and are not intended to limit the scope of the inventions. Indeed, the novel methods and systems described herein may be embodied in a variety of other forms; furthermore, various omissions, substitutions and changes in the form of the methods and systems described herein may be made without departing from the spirit of the inventions. The accompanying claims and their equivalents are intended to cover such forms or modifications as would fall within the scope and spirit of the inventions.
Number | Name | Date | Kind |
---|---|---|---|
10488484 | Miyazaki et al. | Nov 2019 | B2 |
20090082656 | Bayram | Mar 2009 | A1 |
20150123659 | Weingartner | May 2015 | A1 |
20170219672 | Miyazaki | Aug 2017 | A1 |
Number | Date | Country |
---|---|---|
106645945 | May 2017 | CN |
1022576 | Nov 2005 | EP |
WO-2012129022 | Sep 2012 | WO |
WO-2016119054 | Aug 2016 | WO |
Number | Date | Country | |
---|---|---|---|
20230003817 A1 | Jan 2023 | US |