The present invention relates to a magnetic resonance imaging method (hereinafter referred to as an MRI method) and apparatus, in particular to an MRI method and apparatus capable of generating high-quality images of blood stream.
In Arterial Spin Labeling (hereinafter referred to as ASL) which is a technique for exciting in advance a desired region of an object to be examined and imaging a blood stream in downstream of the excited region, a method is disclosed in Non-Patent Document 1 for imaging by conforming direction of gradient magnetic field for slice selection and gradient magnetic field for readout.
Non-Patent Document 1: W. G. Rehwald et al.: GCFP-A New Non-Invasive Non-Contrast Cine Angiography Technique Using Selective Excitation and Global Coherent: Proc. Intl. Soc. Mag. Reson. Med. 11 (2004)
With usage of the imaging method disclosed in Non-Patent Document 1, an excited plane excited by slice selection and an imaging area imaged by accumulating NMR signals intersect orthogonally, whereby making it possible to generate an image of blood vessel extended from the excited plane. In Non-Patent Document 1, images of a blood vessel extended from an excited plane over time are illustrated particularly in FIG. 1.
However, the following problem still remains in the conventional technique disclosed in Non-Patent Document 1. That is, while the regions excited by slicing selection are only excited planes in the conventional technique disclosed in Non-Patent Document 1, the signals from blood flowed out of the excited plane attenuates its intensity over time due to relaxation phenomenon. Thus lowering performance for blood vessel description positioned apart from the excited plane remains as a problem.
Also in Non-Patent Document 1, an imaging sequence that is simulating SSFP for collecting NMR signals is used. In other words, NMR signals are collected while continuously applying a plurality of RF pulses having small flip angles. However, the imaging method by such imaging sequence has a tendency of being influenced by turbulence of signal phase. For this reason, deterioration of images occurs in regions where the static magnetic field is not uniform or in regions having high velocity of blood flow. Such problems are not taken into consideration in the method disclosed in the above-mentioned document.
The objective of the present invention is to provide an MRI method and apparatus capable of imaging a wide range of blood flow with high quality image using the ASL method.
In order to achieve the above-mentioned objective, an MRI method of the present invention comprises:
Also, the MRI apparatus of the present invention comprises:
Hereinafter, system configuration of an MRI apparatus related to the present invention will be described in detail referring to
Configuration of the MRI apparatus is classified broadly by central processing unit (hereinafter referred to as CPU) 1, sequencer 2, transmitting system 3, static magnetic field generating magnet 4, receiving system 5, gradient magnetic field generating system 21, and signal processing system 6.
CPU 1 controls sequencer 2, transmitting system 3, receiving system 5 and signal processing system 6 according to the program set in advance. Sequencer 2 is operated based on control commands from CPU 1, and transmits various commands necessary for collecting image data for generating tomographic images of object 7 to transmitting system 3, gradient magnetic field generating system 21 and receiving system 5.
Transmitting system 3 comprises devices such as high-frequency oscillator 8, modulator 9, irradiating coil 11 and RF shield, amplitude-modifies the reference high-frequency pulse from high-frequency oscillator 8 by the command of sequencer 2, and irradiates a predetermined pulsed electromagnetic waves to the object by amplifying the amplitude-modulated high-frequency pulse via high-frequency amplifier 10 and providing it to irradiating coil 11.
Static magnetic field generating magnet 4 generates a homogeneous static magnetic field around object 7 in a predetermined direction. Inside of static magnetic field generating magnet 4, irradiating coil 11, gradient magnetic field coil 13 and receiving coil 14 are disposed. Gradient magnetic field coil 13 is included in gradient magnetic field generating system 21, receives provision of current from gradient magnetic field source 12, and generates gradient magnetic field under the control of sequencer 2.
Receiving system 5 is for detecting NMR signals emitted by nuclear magnetic resonance of atomic nuclei in the body of the object, and has receiving coil 14, amplifier 15, quadrature detector 16 and A/D converter 17. NMR signals as a response of the object to the electromagnetic waves irradiated from the above-mentioned irradiating coil 11 are detected in receiving coil 14 disposed in the vicinity of the object, inputted to A/D converter via amplifier 15 and quadrature detector 16, and converted into digital quantity. Then the signals converted into digital quantity are transmitted to CPU 1.
Signal processing system 6 comprises an external memory device such as magnetic disk 20 and optical disk 19, and display 18 formed by devices such as CRT. When data from receiving system 5 is inputted, CPU 1 performs process such as signal processing and image reconstruction. Images of the desired fault plane of object 7 which are the result of the above-mentioned process are displayed on display 18, and stored in an external memory device such as magnetic disk 20.
First, an imaging sequence of an MRI method related to embodiment 1 will be described in order using
Non-Patent Document 2: H. Ochi et al.: Dual-frequency amplitude-modulated BURST Imaging, International Society for Magnetic Resonance in Medicine, 5th Scientific Meeting and Exhibition P1824(1997).
Also, 202 indicates application of a gradient magnetic field pulse for slice selection applied along with the application of burst RF pulse 201, and it is to be applied in the same direction as direction of the gradient magnetic field pulse for readout (Gr direction). 203 indicates an inverting pulse (inverting high-frequency magnetic field pulse; π pulse), 204 indicates gradient magnetic field pulse 204 in Gs direction to be applied with inverting pulse 203 and these pulses are for inverting magnetization of the selected region (the region indicated by 208 in
Next,
In
In the present embodiment, after the plurality of excited planes disposed mutually at intervals are excited simultaneously, the signals are collected, and images in each elapsed time are reconstructed from the first application of a burst RF pulse. In this way, in the present embodiment, it is possible to perform imaging of the blood flow with passage of time. The above-mentioned image data in each elapsed time obtained by the imaging sequence shown in
Next, the detail of the imaging composition step will be described referring to
By the above-mentioned diagram, while the image of the blood flowing from of the excited plane is not extracted, it can be recognized that the extracted region of the blood flow is extending from the excited plane to a location downstream as the time passes from time 1 to time 4. However, the blood flow extracted in the image (the blood vessel image) in time 1 time 4 which are generated in the regions between the respective excited planes (the region between excited planes 207-1 and 207-2 is region 1, the region between excited planes 207-2 and 207-3 is region 2 and the region between excited planes 207-3 and 207-4 is region 3) are discrete, since the blood flow is extracted in the image from each of excited plane 207-1˜excited plane 207-4. Given this factor, generation of moving images is executed in the present embodiment for making the image look like the blood is continuously flowing from a specific excited plane (from excited plane 207-1 here) by the method shown in
The moving image in
The image is created by combining the blood extraction image of the target region in each time phase and the blood extraction image at the time when the blood reaches the utmost downstream side in the upstream side region with reference to the targeted region. More specifically, first in time phase 1, a moving image is generated using an image as it is in the reference time. In time phase 2, region 1 of the moving image is generated using the blood flow of time 1 in
Next, in time phase 5, the moving image of region 1 is generated using the blood flow of time phase 3, and the moving image of region 2 is generated using the blood flow of time phase 1 in
Next, concrete procedure of the image composition step for generating the moving image in
(Step 401)
According to the time passed from the previously set reference time (elapsed time), the obtained images are rearranged. The reference time is the time when, for example, the first burst RF pulse 201 is applied, and the image generated by collecting echo signals immediately after the reference time is set as the reference image.
{Step 402}
Counter related to the time phase of the moving image is set as L, and the default value thereof is set as 1.
(Step 403)
Position and number of the excited planes which are necessary parameter in the step described below are derived based on data (imaging condition) such as how burst RF pulses or gradient magnetic field pulses are applied for collecting image data. Or, they also can be calculated based on imaging data. In concrete terms, for example, a threshold value of signal intensity is set in image data, the region having the signal intensity more than the threshold value is detected as the excited plane, and the position and number of the excited plane thereof is calculated. And the excited plane at the utmost upstream point is set as the reference excited plane.
(Step 404)
Counter related to the excited plane is set as n, the default value thereof is set as 1, and the upper limit value is set as Ns being obtained in (step 403).
(Step 405)
Attention is paid on region n which is sandwiched between excited plane n and excited plane (n+1), and time Mn which is the time that the signals of blood that has flowed from excited plane n reaches excited plane (n+1) is identified. Concretely, identification of time Mn is defined by, for example, setting a threshold value to the signal intensity with respect to the pixel which is positioned on the side of excited plane n and adjacent to excited plane (n+1), and defining the time which is more than the threshold value as the time that the blood that has flowed from excited plane n reaches excited plane (n+1) (in the case of an example illustrated in
(Step 406)
The counter related to elapsed time from the reference time for being used upon generation of a moving image of blood flow in the respective regions is set as m. Here, default value of counter m is set as 1, and maximum value of counter m upon extraction of region n is set as Mn.
(Step 407)
A blood vessel image is extracted with respect to region n being sandwiched between excited plane n and excited plane (n+1), from the image after elapsed time m from the reference time. Extraction of a blood vessel image in the present invention is carried out while imaging data is being read out from magnetic disk 20 to CPU 1, and the data of the extracted blood vessel is stored in magnetic disk 20 for the time being.
(Step 408)
The blood vessel image extracted in (step 407) and the moving image at time phase L are synthesized, and set as the moving image at time phase L+1. At that time, the blood vessel image extracted in (step 407) and image data in region 1˜region n−1 (only region 1 in the case that n=2) at time phase L are synthesized. In this regard, however, when the blood vessel image at region 1 is extracted in (step 407), the extracted image is used as it is at time phase L+1. Image composition of moving images in the present step is carried out while the blood vessel image extracted in (step 407) and stored in magnetic disk 20 and the blood vessel image at time phase L are being read out to CPU 1, and the combined result is also stored in magnetic disk 20.
(Step 409)
Counter L related to the time phase of the moving image is incremented.
(Step 410)
In (step 407), parameter m of the elapsed time from the reference time of the target image for extraction of a blood vessel image is compared with upper limit value Mn of m in region n thereof. If m is not the same as Mn step 411 is carried out, and if m is the same as Mn step 412 is to proceed.
(Step 411)
Counter m is incremented by 1, and the step moves to step 47.
(Step 412)
Parameter n related to number of the region is compared with Ns obtained in (step 403). If n is not the same as Ns step 413 is to proceed, and if n is the same as Ns the procedure is ended.
(Step 413)
Counter n is incremented, and the step moves to step 405. The moving image of the conclusively synthesized blood vessel image is displayed, for example, on display 18 and stored in magnetic disk 20.
As mentioned above, according to embodiment 1, after simultaneously exciting a plurality of excited planes at a predetermined interval, in order to obtain the signals produced from the blood flowing from the respective excited planes, it is possible to generate the blood vessel image by collecting the echo signals from the excited planes within a minute distance in a minute period of time. Also, influence of turbulence in the phases of the signals can be minimized. The example of the present embodiment also has an advantage of eliminating influence caused by nonuniformity of magnetic fields since the inverting pulses are applied.
Next, an imaging sequence of an MRI method related to embodiment 2 will be described using FIGS. 5 (a) and (b).
Embodiment 2 has the same advantage as embodiment 1 to prevent the lowering performance for blood vessel description in the downstream region, and to prevent influence due to turbulence of phases. Also, embodiment 2 has another advantage to save imaging time, since there is no application of inverting pulses.
Next, an imaging sequence of an MRI method related to embodiment 3 will be described referring to
Next, an imaging sequence of an MRI imaging method related to embodiment 4 will be described using
Next, an imaging sequence of an MRI method related to embodiment 5 will be described using
In the case that echo signals are collected while alternately changing the polar character of readout gradient magnetic field on the negative side and positive side as shown in the imaging sequence of
Next, image composition steps in an MRI method related to embodiment 6 will be described referring to
In the moving image with respect to each time phase, the state of blood flowing from an upstream region along with passage of time is generated. The image is created combining the blood extraction image of the target region in the respective time phases and the blood extraction image in the respective time phases up to that moment. More specifically, first, the moving image is generated at the reference time in time phase 1. Next, region 1 of the moving image at time phase 2 is generated using the blood flow at time 1 in
Next, an MRI method related to embodiment 7 will be described referring to FIGS. 10 (a), (b), and FIGS. 11 (a) and (b). The present embodiment is an imaging method wherein the excited plane is segmented into a plurality of groups arranged alternately to each other, and the respective groups are alternately excited. First, an image-acquisition step in the present embodiment will be described using FIGS. 10 (a) and (b).
The imaging sequence in
Next, image composition step in the present embodiment will be described using FIGS. 11 (a) and (b).
As mentioned above, according to embodiment 7 compared to embodiments 1˜6, since a plurality of excited planes are divided into a number of groups being alternately arranged and they are alternately excited, effective TR is increased upon imaging the respective imaging regions (the regions sandwiched between the excited planes). As a result, S/N ratio of blood signals are improved, since recovery time of nuclear-magnetization can be longer than the cases of embodiments 1˜5.
Next, an MRI method related to embodiment 8 will be described using FIGS. 12 (a) and (b). Embodiment 8 is an example for imaging an object while the object is being transferred, and position of the excited plane is also moved according to movement of the object.
As shown in
By using such imaging sequence, when the object is being moved with the table, the excited plane can be moved along with the movement thereof, and high-quality image of the blood vessel can be extracted.
Also, as for the image composition method, the method illustrated in the flow chart in
Next, an MRI method related to embodiment 9 will be described using FIGS. 13 (a) and (b). Embodiment 9 is an example for imaging while the object is being transferred, the excited planes are made not to move, and the same position in coordinate system viewing from the MRI apparatus is excited. First,
In order to perform imaging as illustrated in
By using such imaging sequence, it is possible to consistently excite the excited plane of the same position with respect to the MRI apparatus.
As for an image composition method, an image composition process considering the relative position with respect to the object at the excited position can be performed using the method as illustrated in
The present invention is not limited to the above-mentioned embodiments, and various changes may be made without departing from the scope of the invention. For example, methods such as spin echo method, high-speed spin echo method, and gradient echo method may be used.
Also, in the above-mentioned embodiment, while the number of excited planes were four in embodiments 1, 2, 6, 7, and three in embodiments 8 and 9, it may be less than three and more than five in the cases such as embodiments 1, 2, 6, 7, and two or more than four in the cases such as embodiments 8 and 9.
While the number for dividing the excited plane into a plurality of groups is set as two in embodiment 7, the number may be more than three. The respective divided groups of the excited plane do not have to be excited alternately, and the echo signals may be obtained by exciting a certain group for a plurality of times, and after that exciting another group for a plurality of times.
Moreover, a method shown in embodiment 7 for imaging by dividing an excited plane into a plurality of groups may be combined with a method as described in embodiments 8 and 9 for imaging while the object is being transferred.
While the blood vessel image on the excited plane is difficult to extract on the image in the above-mentioned embodiment, it is possible to perform interpolation through extracting the blood vessel on the excited plane in downstream by exciting only the excited planes in upstream.
Combination of the image acquisition step and image combination step in the above-mentioned embodiment does not have to be limited to the above-mentioned combination, and other combinations may be applied.
Also, the respective plurality of excited planes does not have to be arranged in parallel, and they may be slightly tilted.
The MRI apparatus used in the present invention includes a program for implementing the above-mentioned MRI methods stored in devices such as magnetic disk 20. The MRI apparatus used in the present invention is also provided with a memory device such as magnetic disk, with information or data generated in the respective process of the above-mentioned MRI method (parameter for executing the imaging sequence, echo signals obtained by the execution of the imaging sequence, image data reconstructed by the echo signals, time-series image data rearranged in (step 701), image data of the moving images in the respective time phases generated in (step 707), and various types of counter for carrying out the flow chart shown in
Input means is also provided for selectively displaying the generated images or moving images on a device such as display 18. By such means, it is possible for an operator to update the images or moving images displayed on display 18, and to visibly recognize the blood flow with passage of time.
The blood vessel image extracted by the above-mentioned MRI method has a tendency that the pixel value gets larger as getting closer to the excited plane on the upstream side and the pixel value gets smaller as approaching more to downstream, and interpolation may be performed to make it displayed more naturally. In other words, luminance interpolation may be performed to make the pixel value of the blood vessel image on the upstream side of the excited plane small and the pixel value of the blood vessel image on the downstream side large.
The plurality of excitation of the excited plane does not have to be performed simultaneously, and may be sequentially performed from upstream.
Also, the image or moving image from which the excited planes are deleted may be generated, stored, and displayed so that the blood vessel will be clearly visible upon being extracted.
Number | Date | Country | Kind |
---|---|---|---|
2005-038708 | Feb 2005 | JP | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/JP06/02012 | 2/7/2006 | WO | 8/2/2007 |