Abstract:
A nuclear magnetic resonance imaging system divides a measurement space, which is a time band in which measurement signals are generated, into a plurality of partial regions, and measures data under the application of the corresponding gradient magnetic fields to the partial regions so that signals for a blood stream part have their phase consistent at the measurement time point for each partial region. The system synthesizes the data in each partial region, thereby producing image data consisting of only the data with a consistent phase over the entire measurement space, and extracts only the image data for the blood stream part by cancelling image data for a quiescent part from the image data synthesized.

Description:
BACKGROUND OF THE INVENTION 
     This invention relates to an apparatus for nuclear magnetic resonance imaging, and particularly to a method and apparatus for imaging moving liquid. 
     The technique related to this invention is discussed in publication: IEEE TRANSACTIONS ON MEDICAL IMAGAING, Vol. MI-5, No. 3, September 1986, pp. 140-151, and it will be briefed in the following. The operational procedures of imaging blood vessels including a moving part by the nuclear magnetic resonance imaging apparatus will be explained in connection with FIG. 3. In the operation, a pair of pictures called &#34;sensitive picture&#34; and &#34;insensitive picture&#34; are taken, and the latter is subtracted from the former to produce a picture of a blood vessel running system. 
     Initially, a sequence (a), (b), (c) and (d) is followed to produce a sensitive picture. On assumption that a blood stream has a dominant direction coincident with the x direction, the imaging procedures are as follows. 
     (1) A magnetic resonating RF pulse 301 is applied and, at the same time, a z-direction gradient magnetic field 302 is applied so that magnetization on a specific slice plane is excited. Another z-direction gradient magnetic field 303 is applied so as to align the phase of the excited magnetization. 
     (2) Next, a y-direction gradient magnetic field 304 is applied for implementing phase encoding. 
     (b 3) At the same time, x-direction gradient magnetic fields 305 and 306 are applied sequentially in the direction coincident with the dominant blood flow direction, and the MR signal 307 is measured during the application of the magnetic field 306. 
     (4) The measured MR signal 307 is rendered the Fourier transformation process, and a picture is reproduced. 
     The x, y and z directions are the axes of a 3-dimensional orthogonal coordinate system established in the space where a body under test is placed. 
     Next, a sequence (a), (b), (c) and (f) of FIG. 3 is followed to produce an insensitive picture. The above items (1) and (2) are equally applied. 
     Next, as item (b 3), a series of x-direction gradient magnetic fields 308, 309 and 310 are applied, and the MR signal 307 is measured during the application of the magnetic field 301. The measured MR signal 307 is operated by the Fourier transformation process, and a picture is reproduced. 
     The following explains the principle of producing a picture of the blood vessel running system by subtracting the sensitive picture from the insensitive picture obtained as described above. 
     In FIG. 3, phase rotations of magnetization in the body under test caused by the applications of the x-direction gradient magnetic field 305 and magnetic field 306 in its former half period T are called &#34;first phase&#34; and &#34;second phase&#34;, respectively. In a quiescent part, the sum of the first phase and second phase is zero, and the phases are consistent at the center time of echo. This is a signal for the quiescent part. In a blood stream part, magnetization moves at a velocity within the gradient magnetic field, and in this case the sum of the first phase and second phase does not become zero, but results in a variable value depending on the blood flow velocity v. Assuming that the v is constant at all times and it is in the x direction, the sum φ of the first phase and second phase is formulated as follows. ##EQU1## where γ is a constant, G is the gradient of the x-direction magnetic field, and T is the application time length of the x-direction gradient magnetic field. 
     Generally, the blood flow velocity v is variable along the radial direction of blood vessel, i.e., slicing direction, and therefore the phase φ, which is dependent on the blood flow velocity v, also varies along the slicing direction. Consequently, measured signals cancel out each other, resulting in a small amplitude of signals. From this viewpoint, the blood vessel running system is not imaged on the basis of only a sensitive picture. 
     In the insensitive sequence, the phase rotation of magnetization caused by the application of the serial x-direction gradient magnetic field 308 and magnetic field 309 in its former half period T is called &#34;third phase&#34;, and the phase rotation of magnetization caused by the application of the x-direction gradient magnetic field 309 in its latter half period T and magnetic field 310 in its former half period T is called &#34;fourth phase&#34;. The sum of the third and fourth phases in a quiescent part is zero, as in the case of the sensitive sequence mentioned above. In a part of blood stream, if the blood flow velocity v is constant at all times and is in the x direction, the sum φ&#39; of the third and fourth phases is also zero. The reason is that the third phase φ 3  and fourth phase φ 4  are expressed as follows through the same mathematical treatment as for the above case of φ: 
     
         φ.sub.3 =-2πγT.sup.2 Gv 
    
     
         φ.sub.4 =2πγT.sup.2 Gv 
    
     
         And, φ&#39;=φ.sub.3 +φ.sub.4 =0 
    
     Accordingly, the variability of phase φ&#39; due to flow velocities does not occur, and the signal for a blood stream part is produced when phases are consistent again at the center time of echo. Consequently, blood vessels appear in the insensitive picture. By subtracting the sensitive picture from the insensitive picture, the quiescent part is cancelled out and only the blood vessel running system is imaged. 
     Another related technique includes the one in which the angle of magnetization (of a quiescent part) fallen by the magnetic resonating RF pulse is made 360°. Although magnetization of quiescent part is not excited, magnetization of a blood stream part is not fallen completely by 360° due to the flow, and therefore it has a lateral magnetization. After that, by measuring the signal through the induction of echo by following the sequence shown by (a), (b), (c) and (f) in FIG. 3, only the signal for the blood stream part can be produced. The sequence is identical to the case of producing a sensitive picture, except that the magnetic resonating RF pulse 301 of (a) in FIG. 3 is replaced with a 360°-falling RF pulse. As a result, the blood vessel running system is imaged from a single insensitive picture. 
     However, the foregoing related techniques of blood vessel imaging do not comprehend the nullification of phase difference, at the center time of echo at which the measured signal exhibits a maximum value, between the resonance signal for a part having a flow velocity and an expected resonance signal for the same part with the assumption of quiescence (this phase difference will be termed &#34;flow velocity phase&#34;). Therefore, due to various blood flow velocities in reality, there is a great disparty in flow velocity phases as the time point is moved forward and backward from the center time of echo, and therefore even an insensitive picture cannot extract blood vessels clearly. 
     SUMMARY OF THE INVENTION 
     For the achievement of high-fedelity imaging, it is necessary to have signals before and after the center time of echo. The reason is that if signals are measured, as time passes on, always in a state of consistent phase, the rotational speed of phase of each spin, i.e., frequency, is made known. Once the frequency is known, it becomes possible to determine accurately the positions of spins having various frequencies through fast Fourier Transformation (FFT). 
     The present invention is based on this situation, and its prime object is to provide a method and apparatus for nuclear magnetic resonance imaging, in which the flow velocity phase is prevented from fluctuating in the measured data thereby to produce a satisfactory picture of a blood vessel running system. 
     In order to achieve the above objective, the invention recognizes the phase difference between a resonance signal produced for a part of body having a flow velocity and a resonance signal which would be produced for the same part when the body is assumed to be quiescent, and sets a predetermined value for the phase difference (this value is called &#34;target phase value&#34;). 
     The invention includes means which divides a time band in which signals are generated, i.e., measurement space, into a plurality of partial regions, and takes pictures sequentially for the partial regions so that the flow velocity phase of signal in each specific partial region is consistent with the target phase value, and synthesizes image data of the whole measurement space from all image data of pictures, and extracts the signal of the part having a flow velocity from the synthesized data. 
     Namely, the above-mentioned means performs the following. 
     1. It divides a measurement space X into a plurality of partial regions X 1 , X 2 , . . . , X n , and establishes a target phase value φ j  of flow velocity phase in each area X j . 
     2. It takes pictures corresponding to the partial regions X j  as follows. 
     (i) A magnetic resonating RF pulse is applied so that the magnetization of part of a body under test is excited. The time point of pulse application is recognized to be time 0. 
     (ii) A gradient magnetic field G j  which meets the following equation is applied and resonance signals S j  are measured in the whole measurement space. ##EQU2##  where t j  is a representative value selected from the time points of measuring signals in partial areas X j  (this value is called &#34;fetching time point&#34;), T o  is the center time of echo, and γ&#39; is a constant. Values of φ j  and t j  are set to solve the equation (1) for G j , and the G j  to be applied is evaluated. 
     (iii) For the imaging based on the Fourier transformation method using phase encoding, the above items (i) and (ii) are repeated as many times as the number of the times of phase encoding. 
     3. The measured data S j  are weighted so that values of S j  in the partial regions X j  are extracted, and synthesized data S of a piece of picture consisting solely of the data which are obtained, when the flow velocity phase is consistent with the target phase value. 
     4. The synthesized data S is operated by a computation process such that signals from the quiescent part are cancelled, and only a signal from the blood stream part is extracted. 
     The theoretical base of the inventive imaging technique will be explained mathematically in the following. It is assumed that the blood flow velocity in a body under test is constant at all times. It will be appreciated from the following explanation that the flow velocity phase is proportional to the left side of the equation (1). Namely, assuming that, at an arbitrary point Q, protones in blood, i.e., atoms of hydrogen represented by {p k  } (k=1, 2, . . . , n) which flow into close to the point Q at time t=t j  have individual velocities v k . And positions Q k  of P k  at each time point are given as follows. 
     
         Q.sup.k (t)=(Q-t.sub.j ·V.sup.k)+t·V.sup.k ( 2) 
    
     There is the following fundamental expression for the phase rotation speed: 
     
         Phase rotation speed=γ·(strength of magnetic field) 
    
     The strength of magnetic field at a position Q k  is expressed as follows. 
     
         G.sub.j (t)·Q.sup.k (t) 
    
     Accordingly, the degree of rotation of phase of protons p k  at time t j  is indicated by φ j   k  as follows. ##EQU3## 
     At time t j  a quiescent part around the point Q has the phase φ j  given by the following equation. ##EQU4## Accordingly, the flow velocity phase φ around the point Q at time t j  is given by the following equation. ##EQU5## Substituting equation (2) into (3) gives: ##EQU6## It in now appreciated from the above equation that the flow velocity phase φ is proportional to the left side of the equation (1). 
     In the case where {V k  } is virtually constant irrespective of k, with the right side of equation (1) being γ&#39; obtained by the multiplication between the mean value V k  and γ, the right side of the equation (b 4) results as follows. ##EQU7## And, it is appreciated that the actual flow velocity phase φ is equal to θ j . 
     With a prior assumption of θ j  =0 being assumed in the above explanation, when the solution G j  of equation (1) with θ j  =0 is applied, the result of equation (4) indicates that the actual flow velocity phase φ is zero. 
     Based on the foregoing, by application of a gradient magnetic field so that the equation (1) is met, the flow velocity phase can be made consistent with the target phase value at time t j . Accordingly, by extracting values of measured resonance signals S j  for each value of j at time t j  and implementing synthesis, synthesized data, with flow velocity phases being virtually consistent with respective target phase values at all times in the measurement space, can be obtained. Particularly, with a setting being θ j  =0, the right side of the equation (1) becomes zero, and the application timing of G j  is determined independently of {V k  }, and moreover, in this case, the assumption that {V k  } is constant irrespective of each proton k is no more required. 
     In case of θ j  =C·V k , where C is a constant being independent of the value of V k , V k  is the mean value of {V k  }, the equation (4) gives: γ&#39;=V k  ·γ, and therefore the right side of the equation (1) is reduced to as follows. ##EQU8## According, the application timing of G j  is determined independently of {V k  γ. 
     Furthermore, by calculating θ j  from the phase of the picture reproduced from the synthesized data, the mean value of {V k  } for θ j  /C can be evaluated. 
     The nuclear magnetic resonance imaging apparatus according to this invention enables the measurement of phase variations (flow velocity phases) of signals created by flow velocities over the entire measurement space of MR signals so that they are virtually consistent with predetermined values. Consequently, a high quality picture of a blood vessel running system, with little deterioration of picture quality due to flow velocities, can be obtained. 
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS 
     FIGS. 1A and 1B are flowcharts of an embodiment of this invention; 
     FIGS. 2 and 6 are diagrams showing the order of producing the measurement space, partial regions and synthesized data; 
     FIGS. 3a, 3b, 3c, 3d, 3e, and 3f show a sequence chart used in the related technique; 
     FIG. 4 is a block diagram showing an embodiment of the nuclear magnetic resonance imaging apparatus according to this invention; and 
     FIGS. 5a, 5b, 5c, 5d, 5e, and 5f along with FIGS. 7a, 7b, 7c, 7d and 7e are diagrams showing the sequence for RF and gradient fields used in the embodiment of this invention. 
    
    
     DESCRIPTION OF THE PREFERRED EMBODIMENTS 
     (Embodiment 1) 
     FIG. 4 is a systematic block diagram showing an embodiment of the nuclear magnetic resonance imaging apparatus according to this invention. In the figure, a sequence controller 401 is designed to control the generation of various pulses for detecting MR signals from a body under test and control the strength of magnetic field. The sequence controller 401 operates on a transmitter 402 to generate a high-frequency pulse for causing the resonation of a specific nuclide in the body under test. 
     The sequence controller 401 incorporates a computation unit including a CPU, memory and input/output interface, and upon receiving phase value θ j  and fetching time t j  shown in the equation (1) through an input device 409 such as a keyboard, it solves the equation (1) using the internal computation processor to determine the gradient magnetic field G j  suitable for the selected imaging condition, and issues a relevant command to a magnetic field controller 403. It executes the control of gradient magnetic field in steps 101, 102 and 105 in FIG. 1A and steps 101 and 102 in FIG. 1B. 
     The magnetic field controller 403 is designed to issue a signal for controlling arbitrarily the static magnetic field which determines the resonance frequency of NMR signal and the strength and direction of the static magnetic field on the basis of the command from the sequence controller 401. This signal is entered to a magnetic field driver 404, which generates a magnetic field necessary for the measurement on the basis of the signal from the magnetic field controller 403. 
     A receiver 405 receives the NMR signal generated by the body under test, and, on the basis of the measured signal taken from the receiver 405, the processing unit 406 performs picture formation and various computations, so that the reformed picture is displayed on a CRT display unit 407. 
     The processing unit 405 incorporates a computation unit including a CPU, memory and input/output interface, and it functions as follows. The unit 405 stores data at each time of measurement. This operation corresponds to steps 102 and 105 in FIG. 1A and step 102 in FIG. 1B. The unit 405 implements the weighting process for the measured data. This operation corresponds to step 103 in FIG. 1A and step 111 in FIG. 1B. The unit 405 synthesizes necessary imaging data. This operation corresponds to step 104 in FIG. 1A and step 112 in FIG. 1B. The unit 405 extracts only data for a blood stream part. This operation corresponds to step 106 in FIG. 1A and step 113 in FIG. 1B. 
     An ECG (electrocardiogram) 408 detects the heart beat of the body under test so that imaging is timed to the heart beat. 
     The control procedures of this embodiment based on the foregoing arrangement will be described in connection with FIG. 1A. 
     (1) Step 101 divides a measurement space X into X 1 , X 2 , X 3 , X j . Specifically, the space is divided as shown by 205, 206 and 207 in FIG. 2 to establish partial regions of the measurement space X. Indicated by 201, 202 and 203 are total imaging data, and weights are applied to the total imaging data so that data in the partial regions 205, 206 and 207 are taken out. Namely, with the center time of echo being T 0 , for a point t j  which equally divides by n a section t indicative of the range of total measurement space: 
     
         T.sub.0 -T≦t≦T.sub.0 +T 
    
     X j  is determined as follows. 
     
         X.sub.j ={(t, k)εX|t.sub.j ≦t≦t.sub.j +1 } 
    
     where j=1, 2, . . . , n and t 1  =T 0  -T, t n+1  =T 0  +T. In this embodiment, T 0  =22 ms, T=4 ms, and θ j  =0. 
     (2) Step 102 applies gradient magnetic fields G j  for X j  to control the flow velocity phase of signals in X j . This sequence is shown by (a), (b), (c) and (d) in FIG. 5. In the figure, a j  is altered by j=1, j=2, . . . , j=n, and imaging is carried out once for each a j . 
     In this case, a j  and c j  are as follows. ##EQU9## At this time, the equation (1) is met by the following reason. In this embodiment, the equation (1) can be written as follows. ##EQU10## Since θ j  is determined to be zero, the equation is solved for a j , with the left side of the equation (1) being zero, as follows. ##EQU11## On the horizontal axis of measurement time t j , the a j  draws a parabolic curve having a peak of T at the center time t j  =T 0  of the measurement space. 
     Accordingly, time points with the same time length front and back from t=T 0  have the same value of a j . This fact implies that the measured signal has a consistent phase at symmetrical time points with respect to t=T 0 . By increasing the value of j from 1 so that the value of a j  increases, the time point of consistent phase of measurement signal moves from T 0  -T toward T 0 . 
     The procedure of imaging at each time is as follows. The value of k is incremented by 1 from -127 to 128 so that the following items (i)-(iv) are repeated 256 times, and measured signals S j  (t,k) are produced. During the operation, the body under test is assumed to be placed in a static magnetic field of 5000 G. 
     (i) The ECG unit 408 produces a specific wait time from the R-wave of heart beat. 
     (ii) A magnetic resonating RF pulse 501 for falling the magnetization by 90° and a z-direction gradient magnetic field Gz (0.3 G/cm) 502 are applied for a duration of 4 ms so that magnetization of a specific slice plane is excited, and a further z-direction gradient magnetic field -Gz 503 is applied for a duration of 2 ms so that the magnetization has a consistent phase. 
     (iii) A y-direction gradient magnetic field Gy·k/128 (Gy=0.3 G/cm) 504 is applied for a duration of 4 ms. 
     (iv) At the same time, an x-direction gradient magnetic field Gx (Gx=0.3 G/cm) 505 is applied for a duration of a j  ms, and next an x-direction gradient magnetic field -Gx 506 is applied for a duration of T+a j  ms, and next an x-direction gradient magnetic field Gx 507 is applied for a duration of 2T ms and the MR signal 508 generated on the specific slice plane excited by item (ii) is measured. 
     (b 3) Step 103 implements such a weighting process that only data of partial regions X 1 , X 2 , . . . , X n  are taken out from the first to n-th total imaging data S j  (t,k). The weighting factors are given as a function of time. 
     (4) Step 104 uses the weighted measurement data S j  (t,k) provided by the step 103 to produce synthesized data S(t,k): 
     
         S(t,k)=S.sub.j (t,k) if t.sub.j ≦t≦t.sub.j+1 (j=1, 2, . . . , n) 
    
     (5) Step 105 takes a picture in the aforementioned usual sensitive sequence to obtain measured data S&#39;(t, k). In FIG. 5, the gradient magnetic field Gx shown by (f) for a j  =0 in FIG. 5 is used for the measurement. 
     The procedure implements the items (i)-(iv) of the step 102 by incrementing k by 1 from -127 to 128 in exactly the same manner as the conventional sensitive sequence, but without dividing the measurement space. The body under test is placed in a static magnetic field of 5000 G. 
     (6) Step 106 implements the second-order discrete FFT process for two data S(t,k) and S&#39;(t,k) to produce pictures A and A&#39;: 
     
         A=FS (where F signifies the FFT process) 
    
     
         A&#39;=FS&#39; 
    
     Finally, the two pieces of pictures are differentiated to extract a picture B of the blood vessel part: 
     
         B=A-A&#39; 
    
     Based on the extracted signal B of blood vessel part, the picture of the blood vessel running system is displayed. 
     In this manner, the synthesized data S(t,k) provides only signals at consistent flow velocity phases in the total measurement space X, and therefore the picture A has the image of blood stream part, while the picture A&#39; does not, as explained previously for the sensitive picture based on the related technique. The picture B which is the difference of the two pictures has the image of blood stream and does not have the image of quiescent part because of cancellation. Accordingly, a satisfactory picture of blood vessel running system can be produced. 
     Next, another embodiment of this invention will be described in comparison with the preceding embodiment. 
     (Embodiment 2) 
     This embodiment is designed to cancel the quiescent part, only through the weighting process, in the imaging data S j  (t,k) produced by gradient magnetic field G j  corresponding to divided regions. FIG. 1B shows the flowchart of the operation. 
     (1) Steps 101 and 102 are identical to those in Embodiment 1. 
     (2) Step 111 implements the weighting process, provided as a function of time, for the total imaging data. This weighting process is designed such that data of the quiescent part are cancelled each other in obtaining synthesized data S&#34; over the total measurement space X. Accordingly, the synthesized data S&#34; is equivalent to the data B in Embodiment 1. 
     (3) Step 112 produces the above-mentioned synthesized data S&#34;. 
     (4) Step 113 implements the second-order discrete FFT process for the synthesized data S&#34; to produce a picture A&#34;=FS&#34;, and displays it as a picture of the blood vessel running system. 
     (Embodiment 3) 
     This embodiment is applied to the measurement method in which the excitation of magnetization by the RF pulse is fallen by 360°, data for a quiescent part is not detected as a picture from the beginning, and only data for a blood stream part is detected. 
     (1) Step 101 is identical to step 101 of Embodiment 1. 
     (2) Step 102 is identical to step 102 of Embodiment 1, except that the magnetic resonating RF pulse 501 for falling the magnetization by 90° in item (ii) is altered to fall the magnetization by 360°. Magnetization of the quiescent part is returned the original state by the 360° fall, whereas the blood stream part does not fall by 360°, resulting in the generation of a lateral magnetization, and the phase of this lateral magnetization is going to be controlled. 
     (b 3) Step 103 is identical to step 103 of Embodiment 1. 
     (4) Step 104 implements the second-order discrete FFT for the data S(t,k), and a picture A for displaying only the blood vessel running system is produced as follows. 
     
         A=FS 
    
     (Embodiment 4) 
     This embodiment is applied to the case in which the angle of fall of magnetization by the RF pulse is made shallow and the time interval from one imaging to the next imaging, which is called &#34;TR&#34;, is made shorter with the intention of high-speed imaging. 
     In this embodiment, signals are measured as: ##EQU12## As in Embodiment 1. The a j  draws a parabolic curve having a peak T at t=T 0  on the time axis, as has been explained previously. Accordingly, for one value of a j , the measured signal has its phase consistent on both sides of t=T 0 , providing data of two partial regions simultaneously. 
     (1) Step 101 establishes partial regions X j  on both sides of t=T 0  of the measurement space X, as shown by 605, 606 and 607 in FIG. 6. Namely, with a center time T 0  of echo, X j  is established as shown in the following formula for points t j  which equally divide by n the section of T: 
     
         T.sub.0 -T≦t≦T.sub.0 
    
     
         X.sub.j ={(t,k)εX|t.sub.j-1 ≦t≦t.sub.j+1 or t*.sub.j+1 ≦t≦t*.sub.j-1 } 
    
     (where j=1, 2, . . . n, t 0  =T 0  -T, t n  =T 0 , t n+1  =T 0  +(T/n), t* j  =2T 0  -t j ) 
     In this embodiment, measurement is carried out for the setting of T 0  =22 ms, T=4 ms and θ j  =0. 
     (2) Step 102 implements imaging by alternation of j=1, 2, 3, . . . n for a j  in items (a), (b), (c) and (d) in FIG. 5. 
     Here, a j  is as follows. ##EQU13## At this time, the above equation (1) is met. The reason is identical to the case of step 102 in Embodiment 1. 
     The imaging procedures at each occasion are as follows. 
     (i) A body under test is placed in a static magnetic field of 5000 G. 
     (ii) The magnetic resonating RF pulse which falls the magnetization by 20° is applied for a duration of 4 ms, which is followed by a wait time of 26 ms. These operations are repeated 40 times. 
     (iii) k is incremented by 1 from -127 to 128, and the following items (a)-(c) are repeated 256 times to obtain measured signals S j  (t,k). 
     (a) In the sequence of FIG. 5, a magnetic resonating RF pulse 501 which falls the magnetization by 20° and a z-direction gradient magnetic field Gz (0.3 G/cm) 502 are applied for a duration of 4 ms so that the magnetization of a specific slice plane is excited. Furthermore, a z-direction gradient magnetic field -Gz 503 is applied for a duration of 2 ms so that the phase of magnetization is made consistent. 
     (b) After that, a y-direction gradient magnetic field Gy·(k/128) (Gy=0.3 G/cm) 504 is applied for a duration of 4 ms. 
     (c) At the same time, an x-direction gradient magnetic field Gx (Gx=0.3 g/cm) 505 is applied for a duration of a j  ms, and next an x-direction gradient magnetic field -Gx 506 is applied for a duration of T+a j  ms, and next an x-direction gradient magnetic field Gx 507 is applied for a duration of 2T ms and at the same time the MR signal 508 is measured. The operation is followed by a wait time of 4 ms. 
     (b 3) Step 103 implements such a weighting process as taking out only data in partial regions X 1 , X 2 , . . . , X n  from the first to n-th total imaging data S j  (t,k). Weighting factors are given as a function of time. 
     (4) Step 104 uses the weighted measurement data S j  (t,k) obtained in the previous step 103 to produce a piece of synthesized data S(t,k) over the entire measurement space X, as shown by the following formula. ##EQU14## where function μ(τ) is as follows. ##EQU15## 
     (5) Step 105 takes pictures in the usual sensitive sequence shown by (a), (b), (c) and (f) in FIG. 5 to obtain measured data S&#39;(t,k). 
     The procedure does not divide the measurement space, and implements items (i)-(iii) of step 102 of this embodiment. 
     (6) Step 106 is identical to step 106 of Embodiment 1. Data S(t,k) and S&#39;(t,k) obtained in steps 104 and 105 are rendered the second-order FFT process to obtain pictures as follows. 
     
         A=FS 
    
     
         A&#39;=FS&#39; 
    
     By obtaining the difference B=A-A&#39; of both image data, image data of the blood vessel running system is extracted. 
     (Embodiment 5) 
     This embodiment is applied to the imaging method in which the RF pulse is applied twice and the phase of magnetization is rotated by 180° by the second RF pulse so that unevennes of the static magnetic field is cancelled. 
     (1) Step 101 is identical to step 101 of Embodiment 4, and it sets up partial regions X j  of the measurement space X on both sides of t=T 0 . T 0  is 28 ms, T is 4 ms, and the target phase θ j  is defined to be θ j  =γ·(average flow velocity v k ) 
     (2) Step 102 implements imaging by alternation of j=1, 2, 3, . . . , n for b j  in items (a), (b), (c) and (d) in FIG. 7. b j  is as follows. ##EQU16## 
     In this case the equation (1) is met, the reason of which will be explained in the following. 
     As in step 102 of Embodiment 1, the left side of the equation (1) can be written as follows. 
     
         G/2{2T.sup.2 -(T.sub.0 -t.sub.j).sup.2 -2b.sub.j.sup.2 } 
    
     The right side of the equation (1) can be written as follows. ##EQU17## And therefore: 
     
         G/2 {2T.sup.2 -(T.sub.0 -t.sub.j).sup.2 -2b.sub.j.sup.2 }=1 
    
     Solving the equation for b j  gives: ##EQU18## 
     In this case, the gradient magnetic fields (707, 708), which has been applied before the pulse for falling the magnetization by 180° in 702 of FIG. 7, has its polarity reversed before calculating the equation (1). 
     The imaging procedures at each occasion are as follows. k is incremented by 1 from -127 to 128 so that the following items (i)-(vi) are repeated 256 times to obtain measured signals S j  (t,k). The body under test is placed in a static magnetic field of 5000 G. 
     (i) The ECG operates to wait for a specific time length from the R-wave of heart beat. 
     (ii) After that, a magnetic resonating RF pulse 701 which falls the magnetization by 90° and a z-direction gradient magnetic field Gz (0.3 G/cm) 703 are applied for a duration of 4 ms so that the magnetization of a specific slice plane is excited, and furthermore a z-direction gradient magnetic field -Gz 704 is applied for a duration of 2 ms so that the phase of magnetization is made consistent. 
     (iii) After that, an x-direction gradient magnetic field -G (0.3 G/cm) 707 is applied for a duration of b j  ms, and next an x-direction gradient magnetic field Gx 708 is applied for a duration of b j  ms. 
     (iv) After that, a magnetic resonating RF pulse 702 which falls the magnetization by 180° and a z-direction gradient magnetic field Gz (0.3 G/cm) 705 are applied for a duration of 4 ms so that the magnetization of the slice plane is reversed. 
     (v) After that, a y-direction gradient magnetic field Gy·k/128 (Gy=0.3 G/cm) 706 is applied for a duration of 4 ms. 
     (vi) At the same time, an x-direction gradient magnetic field -Gx 709 is applied for a duration of T ms, and next an X-direction gradient magnetic field Gx 710 is applied for a duration of 2T ms and simultaneously the MR signal 711 is measured. 
     (b 3) Step 103 is identical to step 103 of Embodiment 4, and it implemets such a weighting process as for taking out only data of partial regions X 1 , X 2 , . . . , X n  for the whole first through n-th imaging data S j  (t,k). Weighting factors are given as a function of time. 
     (4) Step 104 is identical to step 104 of Embodiment 4, and it produces synthesized data S(t,k) over the entire measurement space X. 
     (5) Step 105 implements imaging in the sequence of (a), (b), (c) and (f) of FIG. 7 to obtain measured data S&#39;(t,k). This operation is equivalent to the usual sensitive sequence, and it is an imaging method without the detection of signal from a blood stream part. The procedures are identical to the case of step 102 in this embodiment, except that Gx of (f) in FIG. 7 is used as Gy. 
     (6) Step 106 is identical to step 106 of Embodiment 6, and it takes differences among image data obtained to extract only image data of the blood vessel running system. 
     Moreover, in this embodiment, the phase of signals on the picture of blood vessel running system is proportional to the average blood flow velocity v k . Accordingly, by calculating θ j  /γ, the average flow velocity v k  can be evaluated.