JPH08117202A - Image reconfiguration method in mri device and mri device - Google Patents
Image reconfiguration method in mri device and mri deviceInfo
- Publication number
- JPH08117202A JPH08117202A JP6258337A JP25833794A JPH08117202A JP H08117202 A JPH08117202 A JP H08117202A JP 6258337 A JP6258337 A JP 6258337A JP 25833794 A JP25833794 A JP 25833794A JP H08117202 A JPH08117202 A JP H08117202A
- Authority
- JP
- Japan
- Prior art keywords
- group
- data
- axis direction
- echo
- echo data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Description
【０００１】[0001]
【産業上の利用分野】この発明は、ＭＲＩ（Magnetic R
esonance Imaging）装置における画像再構成方法および
ＭＲＩ装置に関する。さらに詳しくは、データ収集処理
と演算処理とを並列実行可能とすることにより、目的画
像を得られるまでの時間を短縮できるようにしたＭＲＩ
装置における画像再構成方法およびＭＲＩ装置に関す
る。BACKGROUND OF THE INVENTION This invention is applied to MRI (Magnetic R
The present invention relates to an image reconstruction method and an MRI apparatus in an esonance imaging apparatus. More specifically, by allowing the data collection processing and the arithmetic processing to be executed in parallel, the MRI that can shorten the time until the target image is obtained
The present invention relates to an image reconstruction method in an apparatus and an MRI apparatus.
【０００２】[0002]
【従来の技術】図１２は、従来のＭＲＩ装置における２
次元の画像再構成処理の動作を示すフローチャートであ
る。ステップＥ１では、操作者がパルスシーケンスのパ
ラメータを設定する。このパラメータのうちの位相軸方
向のＦＯＶ（Ｆield Of View），エンコード数Ｎ，勾配
パルスの大きさＧ，勾配パルスの幅Ｔによって、 Δθ＝２π／（Ｇ・Ｔ・ＦＯＶ・Ｎ） から位相軸方向のエンコードステップ量Δθが求められ
る。ステップＥ２では、データ番号カウンタｋ＝１とす
る。ステップＥ３では、位相軸方向のｋ番目のエンコー
ド量θ（ｋ）を、 θ（ｋ）＝−π／（Ｇ・Ｔ・ＦＯＶ）＋ｋ・Δθ により算出する。ステップＥ４では、位相軸方向のエン
コード量θ（ｋ）で、エコーデータＲ（ｋ）を収集す
る。ステップＥ５では、データ番号カウンタｋ＝Ｎか否
かをチェックする。ｋ＝ＮでないならステップＥ６に進
み、ｋ＝Ｎならデータ収集処理を終了する。ステップＥ
６では、ｋ＝ｋ＋１とし、前記ステップＥ３に戻る。上
記ステップＥ３〜Ｅ６により、図１３に示す２次元ｋ−
空間のエコーデータＲ（１）〜Ｒ（Ｎ）を順に収集する
が、一つのエコーデータを収集している間に、それより
前に収集したエコーデータについての周波数軸方向の１
次元フーリエ変換をステップＦ１で並行して行って中間
データｆ（ｋ）を算出し、その中間データｆ（ｋ）を保
存しておく。そして、前記ステップＥ５でｋ＝Ｎとなっ
たなら、ステップＦ２で、中間データｆ（１）〜ｆ
（Ｎ）に対して、位相軸方向の１次元フーリエ変換を行
う。以上のステップＥ１〜Ｅ６のデータ収集処理および
ステップＦ１〜Ｆ２の演算処理によって、目的画像が再
構成される。2. Description of the Related Art FIG.
It is a flow chart which shows operation of three-dimensional image reconstruction processing. In step E1, the operator sets the parameters of the pulse sequence. Depending on the FOV (Field Of View) in the phase axis direction, the encoding number N, the gradient pulse size G, and the gradient pulse width T among these parameters, Δθ = 2π / (G · T · FOV · N) The encoding step amount Δθ in the direction is obtained. At step E2, the data number counter k = 1. In step E3, the kth encoding amount θ (k) in the phase axis direction is calculated by θ (k) = − π / (G · T · FOV) + k · Δθ 2. At step E4, the echo data R (k) is collected with the encode amount θ (k) in the phase axis direction. At step E5, it is checked whether or not the data number counter k = N. If not k = N, the process proceeds to step E6, and if k = N, the data collection process ends. Step E
In step 6, k = k + 1 is set, and the process returns to step E3. By the above steps E3 to E6, the two-dimensional k-
Spatial echo data R (1) to R (N) are collected in order, but while collecting one echo data, 1 in the frequency axis direction for echo data collected before it is collected.
The dimensional Fourier transform is performed in parallel in step F1 to calculate the intermediate data f (k), and the intermediate data f (k) is stored. If k = N in step E5, the intermediate data f (1) to f (f) is calculated in step F2.
For (N), one-dimensional Fourier transform in the phase axis direction is performed. The target image is reconstructed by the data collection processing in steps E1 to E6 and the arithmetic processing in steps F1 to F2.
【０００３】図１４は、従来のＭＲＩ装置における３次
元の画像再構成処理の動作を示すフローチャートであ
る。ステップＧ１では、操作者がパルスシーケンスのパ
ラメータを設定する。このパラメータのうちの位相軸方
向のＦＯＶk ，エンコード数Ｎ，勾配パルスの大きさＧ
k，勾配パルスの幅Ｔkによって、 Δθ＝２π／（Ｇk・Ｔk・ＦＯＶk・Ｎ） から位相軸方向のエンコードステップ量Δθが求められ
る。又、スライス軸方向のＦＯＶs ，エンコード数Ｓ，
勾配パルスの大きさＧs，勾配パルスの幅Ｔsによって、 Δφ＝２π／（Ｇs・Ｔs・ＦＯＶs・Ｓ） からスライス軸方向のエンコードステップ量Δφが求め
られる。ステップＧ２では、位相軸方向データ番号カウ
ンタｋ＝１とする。ステップＧ３では、位相軸方向のｋ
番目のエンコード量θ（ｋ）を、 θ（ｋ）＝−π／（Ｇk・Ｔk・ＦＯＶk）＋ｋ・Δθ により算出する。ステップＧ４では、スライス軸方向デ
ータ番号カウンタｓ＝１とする。ステップＧ５では、ス
ライス軸方向のｓ番目のエンコード量φ（ｓ）を、 φ（ｓ）＝−π／（Ｇs・Ｔs・ＦＯＶs）＋ｓ・Δφ により算出する。ステップＧ６では、位相軸方向のエン
コード量θ（ｋ），スライス軸方向のエンコード量φ
（ｓ）で、エコーデータＲ{s}(k)を収集する。ステップ
Ｇ７では、スライス軸方向データ番号カウンタｓ＝Ｓか
否かをチェックする。ｓ＝ＳでないならステップＧ８に
進み、ｓ＝ＳならステップＧ９に進む。ステップＧ８で
は、ｓ＝ｓ＋１とし、前記ステップＧ５に戻る。ステッ
プＧ９では、位相軸方向データ番号カウンタｋ＝Ｎか否
かをチェックする。ｋ＝ＮでないならステップＧ１０に
進み、ｋ＝Ｎならデータ収集処理を終了する。ステップ
Ｇ１０では、ｋ＝ｋ＋１とし、前記ステップＧ３に戻
る。上記ステップＧ３〜Ｇ１０により、図１５に示す３
次元ｋ−空間のエコーデータＲ{1}(1)〜Ｒ{4}(4)を順に
収集するが、ある位相軸方向のエンコード量θ（ｋ）に
ついてスライス軸方向のエンコード量φ（ｓ）を変えて
エコーデータを収集している間に、それより前に収集し
た位相軸方向のエンコード量θ（ｋ−１）についてのエ
コーデータのスライス軸方向の１次元フーリエ変換をス
テップＨ１で並行して行って中間データを算出し、その
中間データを保存しておく。そして、前記ステップＧ９
でｋ＝Ｎとなったなら、ステップＨ２で、前記中間デー
タに対して、周波数軸および位相軸方向のフーリエ変換
を行う。以上のステップＧ１〜Ｇ１０のデータ収集処理
およびステップＨ１〜Ｈ２の演算処理によって、各スラ
イスについての目的画像が再構成される。FIG. 14 is a flowchart showing the operation of a three-dimensional image reconstruction process in a conventional MRI apparatus. In step G1, the operator sets the parameters of the pulse sequence. Of these parameters, the FOVk in the phase axis direction, the number of encodes N, the magnitude G of the gradient pulse
The encoding step amount Δθ in the phase axis direction is obtained from Δθ = 2π / (Gk · Tk · FOVk · N) according to k and the gradient pulse width Tk. Also, FOVs in the slice axis direction, the number of encodes S,
The encoding step amount Δφ in the slice axis direction is obtained from Δφ = 2π / (Gs · Ts · FOVs · S) by the gradient pulse size Gs and the gradient pulse width Ts. In step G2, the phase axis direction data number counter k = 1. In step G3, k in the phase axis direction
The th encoding amount θ (k) is calculated by θ (k) = − π / (Gk · Tk · FOVk) + k · Δθ. In step G4, the slice axis direction data number counter s = 1. In step G5, the sth encoding amount φ (s) in the slice axis direction is calculated by φ (s) = − π / (Gs · Ts · FOVs) + s · Δφ. In step G6, the encode amount θ (k) in the phase axis direction and the encode amount φ in the slice axis direction
At (s), echo data R {s} (k) is collected. In step G7, it is checked whether or not the slice axis direction data number counter s = S. If s = S is not satisfied, the process proceeds to step G8, and if s = S, the process proceeds to step G9. In step G8, s = s + 1 is set, and the process returns to step G5. In step G9, it is checked whether or not the phase axis direction data number counter k = N. If not k = N, the process proceeds to step G10, and if k = N, the data collection process ends. In step G10, k = k + 1 is set, and the process returns to step G3. By the above steps G3 to G10, 3 shown in FIG.
The echo data R {1} (1) to R {4} (4) in the dimension k-space are collected in order, and the encoding amount φ (s) in the slice axis direction for a certain encoding amount θ (k) in the phase axis direction. While collecting the echo data while changing the parameter, the one-dimensional Fourier transform in the slice axis direction of the echo data for the encoding amount θ (k−1) in the phase axis direction, which is collected earlier, is performed in parallel in step H1. Then, the intermediate data is calculated and the intermediate data is stored. Then, the step G9
If k = N, the Fourier transform in the frequency axis and phase axis directions is performed on the intermediate data in step H2. The target image for each slice is reconstructed by the data collection processing in steps G1 to G10 and the calculation processing in steps H1 to H2.
【０００４】図１６は、従来のＭＲＩ装置における３次
元の画像再構成処理の別の動作を示すフローチャートで
ある。ステップＧ１’では、操作者がパルスシーケンス
のパラメータを設定する。このパラメータにより上記と
同様に位相軸方向のエンコードステップ量Δθおよびス
ライス軸方向のエンコードステップ量Δφが求められ
る。ステップＧ２’では、スライス軸方向データ番号カ
ウンタｓ＝１とする。ステップＧ３’では、スライス軸
方向のｓ番目のエンコード量φ（ｓ）を算出する。ステ
ップＧ４’では、位相軸方向データ番号カウンタｋ＝１
とする。ステップＧ５’では、位相軸方向のｋ番目のエ
ンコード量θ（ｋ）を算出する。ステップＧ６’では、
位相軸方向のエンコード量θ（ｋ），スライス軸方向の
エンコード量φ（ｓ）で、エコーデータＲ{s}(k)を収集
する。ステップＧ７’では、位相軸方向データ番号カウ
ンタｋ＝Ｎか否かをチェックする。ｋ＝Ｎでないならス
テップＧ８’に進み、ｋ＝ＮならステップＧ９’に進
む。ステップＧ８’では、ｋ＝ｋ＋１とし、前記ステッ
プＧ５’に戻る。ステップＧ９’では、スライス軸方向
データ番号カウンタｓ＝Ｓか否かをチェックする。ｓ＝
ＳでないならステップＧ１０’に進み、ｓ＝Ｓならデー
タ収集処理を終了する。ステップＧ１０’では、ｓ＝ｓ
＋１とし、前記ステップＧ３’に戻る。上記ステップＧ
３’〜Ｇ１０’により、図１５に示す３次元ｋ−空間の
エコーデータＲ{1}(1)〜Ｒ{4}(4)を順に収集し、全エコ
ーデータの収集を完了した後、全エコーデータのスライ
ス軸方向の１次元フーリエ変換をステップＨ１’で行っ
て中間データを算出し、その中間データに対して周波数
軸および位相軸方向のフーリエ変換をステップＨ２’で
行う。以上のステップＧ１’〜Ｇ１０’のデータ収集処
理およびステップＨ１’〜Ｈ２’の演算処理によって、
各スライスについての目的画像が再構成される。FIG. 16 is a flow chart showing another operation of the three-dimensional image reconstruction processing in the conventional MRI apparatus. At step G1 ', the operator sets the parameters of the pulse sequence. With this parameter, the encode step amount Δθ in the phase axis direction and the encode step amount Δφ in the slice axis direction are obtained similarly to the above. At step G2 ′, the slice axis direction data number counter s = 1. In step G3 ′, the sth encoding amount φ (s) in the slice axis direction is calculated. At step G4 ′, the phase axis direction data number counter k = 1
And In step G5 ′, the kth encoding amount θ (k) in the phase axis direction is calculated. In step G6 ',
The echo data R {s} (k) is collected with the encode amount θ (k) in the phase axis direction and the encode amount φ (s) in the slice axis direction. In step G7 ', it is checked whether or not the phase axis direction data number counter k = N. If not k = N, the process proceeds to step G8 ', and if k = N, proceeds to step G9'. In step G8 ', k = k + 1 is set, and the process returns to step G5'. At Step G9 ', it is checked whether or not the slice axis direction data number counter s = S. s =
If it is not S, the process proceeds to step G10 ', and if s = S, the data collection process is ended. In step G10 ', s = s
The value is set to +1 and the process returns to step G3 '. Step G above
3 ′ to G10 ′, the echo data R {1} (1) to R {4} (4) in the three-dimensional k-space shown in FIG. 15 are collected in order, and after the collection of all echo data is completed, One-dimensional Fourier transform of the echo data in the slice axis direction is performed in step H1 ′ to calculate intermediate data, and Fourier transform in the frequency axis and phase axis directions is performed on the intermediate data in step H2 ′. By the data collection process of steps G1 ′ to G10 ′ and the calculation process of steps H1 ′ to H2 ′,
The target image for each slice is reconstructed.
【０００５】[0005]
【発明が解決しようとする課題】図１２に示す従来のＭ
ＲＩ装置における２次元の画像再構成処理では、図１３
の２次元ｋ−空間の全エコーデータの収集を完了した
後、ステップＦ２で位相軸方向の１次元フーリエ変換の
演算を開始している。つまり、ステップＦ２の演算処理
がデータ収集処理と並列実行されない。このため、ステ
ップＦ２の演算時間が全て待ち時間になり、画像が表示
されるまでの所要時間が長くなってしまう問題点があ
る。The conventional M shown in FIG. 12 is used.
In the two-dimensional image reconstructing process in the RI device, FIG.
After completing the collection of all the echo data in the two-dimensional k-space of, the calculation of the one-dimensional Fourier transform in the phase axis direction is started in step F2. That is, the calculation process of step F2 is not executed in parallel with the data collection process. For this reason, there is a problem that the calculation time of step F2 is all waiting time, and the time required until an image is displayed becomes long.
【０００６】一方、図１４に示す従来のＭＲＩ装置にお
ける３次元の画像再構成処理では、ある位相軸方向のエ
ンコード量を固定しておいて、スライス軸方向のエンコ
ード量を先に変化させることで、ステップＨ１の演算処
理をデータ収集処理と並列実行可能にしている。ところ
が、スピンの動的平衡状態を利用する一部のシーケンス
やプリパレーションパルスを付加したシーケンスなどで
は、シーケンスの性質から、スライス軸方向のエンコー
ド量を先に変化させることが困難であり、図１６に示す
３次元の画像再構成処理のように、あるスライス軸方向
のエンコード量を固定しておいて、位相軸方向のエンコ
ード量を先に変化させることが行なわれる。しかし、図
１６に示す従来のＭＲＩ装置における３次元の画像再構
成処理では、演算処理（ステップＨ１’，Ｈ２’）がデ
ータ収集処理と並列実行されないため、演算時間が全て
待ち時間になり、画像が表示されるまでの所要時間が長
くなってしまう問題点がある。On the other hand, in the three-dimensional image reconstruction processing in the conventional MRI apparatus shown in FIG. 14, the encoding amount in a certain phase axis direction is fixed and the encoding amount in the slice axis direction is changed first. , The calculation processing of step H1 can be executed in parallel with the data collection processing. However, it is difficult to change the encoding amount in the slice axis direction first because of the nature of the sequence in a part of the sequence that utilizes the dynamic equilibrium state of spins or a sequence to which a preparation pulse is added. Like the three-dimensional image reconstruction processing shown, the encoding amount in a certain slice axis direction is fixed and the encoding amount in the phase axis direction is changed first. However, in the three-dimensional image reconstruction process in the conventional MRI apparatus shown in FIG. 16, the calculation process (steps H1 ′ and H2 ′) is not executed in parallel with the data collection process, so the calculation time is all waiting time, There is a problem that it takes a long time to display.
【０００７】そこで、この発明の目的は、データ収集処
理と演算処理とを並列実行可能とすることにより、予め
演算できる箇所を処理し、全データ収集後に必要な演算
を減らすことで、目的画像を得られるまでの時間を短縮
できるようにしたＭＲＩ装置における画像再構成方法お
よびＭＲＩ装置を提供することにある。Therefore, an object of the present invention is to perform a data acquisition process and an operation process in parallel, thereby processing a place where an operation can be performed in advance, and reducing an operation required after collecting all data, thereby obtaining an object image. An object of the present invention is to provide an image reconstructing method in an MRI apparatus and an MRI apparatus capable of shortening the time until it is obtained.
【０００８】[0008]
【課題を解決するための手段】第１の観点では、この発
明は、あるエンコード方向にエコーデータをグループ分
けし、そのエンコード方向のエンコード量を変化させて
第１のグループから順にエコーデータを収集すると共
に、あるグループのエコーデータの収集と並行してそれ
より前にエコーデータの収集を完了したグループについ
て前記エンコード方向のフーリエ変換を行ってグループ
データを取得し、全てのグループのグループデータを加
算して、全データについての前記エンコード方向のフー
リエ変換結果とすることを特徴とするＭＲＩ装置におけ
る画像再構成方法を提供する。In a first aspect, the present invention divides echo data into groups in a certain encoding direction, changes the encoding amount in the encoding direction, and collects the echo data in order from the first group. At the same time, in parallel with the collection of echo data of a certain group, Fourier transform in the encoding direction is performed on the group for which echo data collection was completed before that to obtain group data, and the group data of all groups are added. Then, the image reconstruction method in the MRI apparatus is provided, which is the Fourier transform result of all the data in the encoding direction.
【０００９】第２の観点では、この発明は、２次元ｋ−
空間の全エコーデータを，位相軸方向の各成分が各グル
ープに同じ数だけ含まれるようにインタリーブして，位
相軸方向に第１グループから第Ｍグループまでにグルー
プ分けして第１グループのエコーデータから順にデータ
を収集するグループ分けデータ収集手段と、あるグルー
プのエコーデータの収集と並行してそれより前に収集を
完了したグループのエコーデータについて周波数軸軸方
向のフーリエ変換および位相軸方向のフーリエ変換を行
ってグループデータを取得するグループデータ取得手段
と、第１グループから第Ｍグループのグループデータを
加算するグループデータ加算手段とを具備したことを特
徴とするＭＲＩ装置を提供する。第３の観点では、この
発明は、上記構成のＭＲＩ装置において、前記グループ
データ加算手段は、グループデータを取得する毎にそれ
までに取得したグループデータを累積加算し、さらに、
前記累積加算したデータに基づいて画像を表示する画像
表示手段を具備することを特徴とするＭＲＩ装置を提供
する。In a second aspect, the invention is a two-dimensional k-
All spatial echo data are interleaved so that the same number of each component in the phase axis direction is included in each group, and are grouped in the phase axis direction from the first group to the Mth group to echo the first group. A grouping data collection means for collecting data in order from the data, and a Fourier transform in the frequency axis direction and a Fourier transform in the phase axis direction for the echo data of the group that has been completed in parallel with the collection of the echo data of a certain group There is provided an MRI apparatus characterized by comprising group data acquisition means for performing Fourier transform to acquire group data, and group data addition means for adding group data of the first group to the Mth group. According to a third aspect of the present invention, in the MRI apparatus configured as described above, the group data adding unit cumulatively adds the group data acquired up to that time each time the group data is acquired, and further,
There is provided an MRI apparatus comprising an image display means for displaying an image based on the cumulatively added data.
【００１０】第４の観点では、この発明は、３次元ｋ−
空間の全エコーデータを，スライス軸方向の各成分が各
グループに同じ数だけ含まれるようにインタリーブし
て，スライス軸方向に第１グループから第Ｑグループま
でにグループ分けし第１グループのエコーデータから順
にデータを収集するグループ分けエコーデータ収集手段
と、あるグループのエコーデータの収集と並行してそれ
より前にエコーデータの収集を完了したグループのデー
タについてスライス軸方向のフーリエ変換を行ってグル
ープデータを取得するグループデータ取得手段と、第１
グループから第Ｑグループのグループデータを加算する
グループデータ加算手段とを具備したことを特徴とする
ＭＲＩ装置を提供する。In a fourth aspect, the invention is a three-dimensional k-
All the spatial echo data are interleaved so that each group in the slice axis direction is included in the same number in each group, and are grouped in the slice axis direction from the 1st group to the Qth group to obtain the echo data of the 1st group. Grouped echo data collection means that collects data in order from 1 to 4, and in parallel with the collection of echo data of a certain group, the Fourier transform in the slice axis direction is performed on the data of the group that completed echo data collection before that Group data acquisition means for acquiring data, first
Provided is an MRI apparatus comprising group data addition means for adding group data of a Qth group from a group.
【００１１】第５の観点では、この発明は、上記構成の
ＭＲＩ装置において、各グループについてのフーリエ変
換結果を加算前に整合させるため、前記グループデータ
取得手段または前記グループデータ加算手段により、各
グループデータの位相補償演算を行うことを特徴とする
ＭＲＩ装置を提供する。According to a fifth aspect of the present invention, in the MRI apparatus configured as described above, the group data acquisition means or the group data addition means is used to match the Fourier transform results for each group before addition. Provided is an MRI apparatus characterized by performing phase compensation calculation of data.
【００１２】[0012]
【作用】上記第１の観点によるＭＲＩ装置における画像
再構成方法では、あるエンコード方向にエコーデータを
グループ分けし、そのエンコード方向のエンコード量を
変化させて第１のグループから順にエコーデータを収集
する。そして、第１のグループのエコーデータの収集が
完了すると、第２のグループのエコーデータの収集を開
始するが、それと並行して、第１のグループについて前
記エンコード方向のフーリエ変換を行ない、第１のグル
ープについてのグループデータを取得する。以下同様に
して、あるグループのエコーデータの収集と並行して、
それより前にエコーデータの収集を完了したグループに
ついて前記エンコード方向のフーリエ変換を行ってグル
ープデータを取得する。そして、全てのグループのグル
ープデータを加算して（位相補償演算し、加算して）、
全データについての前記エンコード方向のフーリエ変換
結果とする。これによれば、あるグループについてのデ
ータ収集処理と，その前にデータ収集したグループにつ
いての演算処理とを並列実行できるため、全データを収
集してから目的画像を得るまでに必要な演算数を減少で
きるようになる。In the image reconstruction method in the MRI apparatus according to the first aspect, the echo data is divided into groups in a certain encoding direction, the encoding amount in the encoding direction is changed, and the echo data is collected in order from the first group. . Then, when the collection of the echo data of the first group is completed, the collection of the echo data of the second group is started. At the same time, the Fourier transform in the encoding direction is performed for the first group, Get group data for a group of. Similarly, in parallel with the collection of echo data of a group,
Fourier transform in the encoding direction is performed for the group for which echo data collection has been completed before that, and group data is acquired. Then, add the group data of all the groups (phase compensation calculation and add),
The result is the Fourier transform result in the encoding direction for all data. According to this, since the data collection process for a certain group and the calculation process for the group that collected data before that can be executed in parallel, the number of calculations required from collecting all data to obtaining the target image can be performed. You can decrease.
【００１３】上記原理を数式で説明すると、次のように
なる。フーリエ変換前のデータをａ（ｔ）；ｉ＝１〜Ｎ
とし、全データのフーリエ変換結果をＡ（ω）とし、グ
ループ数をＭとし、各成分が各グループに同じ数だけに
含まれるようにインタリーブしてグループ分けを行なっ
たとするとき、The above principle will be described below in terms of mathematical expressions. Data before Fourier transform is a (t); i = 1 to N
When the Fourier transform result of all data is A (ω), the number of groups is M, and each component is interleaved so that each group is included in the same number, grouping is performed,
【００１４】[0014]
【数１】 [Equation 1]
【００１５】となる。上式で、Ａｍ（ω）は各グループ
についてＮ／Ｍ個のデータで部分的にフーリエ変換を行
ったグループデータであるから、各グループデータに位
相補償を行なって加算すれば、全データのフーリエ変換
結果Ａ（ω）が得られることが判る。このとき、Ａｍ
（ω）のフーリエ変換はＮ／Ｍ個のデータで演算し、各
グループの加算においてはこれを周期関数とみなしてＮ
個のデータ域に拡張して行う。各グループのフーリエ変
換Ｎ／Ｍ個のデータに対する演算でよいので、全データ
収集後に必要な演算は、最終グループのＮ／Ｍ個のデー
タに対するフーリエ変換とＮ個のデータの位相補償と累
積加算だけになる。[0015] In the above equation, Am (ω) is group data obtained by partially performing Fourier transform on N / M pieces of data for each group. Therefore, if phase compensation is performed on each group data and addition is performed, Fourier of all data is obtained. It can be seen that the conversion result A (ω) is obtained. At this time, Am
The Fourier transform of (ω) is calculated with N / M data, and in addition of each group, this is regarded as a periodic function and N
It is performed by expanding the data area. Since Fourier transform N / M data of each group is sufficient, only the Fourier transform for the final group of N / M data, phase compensation and cumulative addition of N data are necessary after collecting all data. become.
【００１６】上記第２の観点（または第５の観点）によ
るＭＲＩ装置では、２次元ｋ−空間の全エコーデータ
を，位相軸方向の各成分が各グループに同じ数だけ含ま
れるようにインタリーブして，位相軸方向に第１グルー
プから第Ｍグループまでにグループ分けし、第１グルー
プのエコーデータから順にデータを収集し、あるグルー
プのエコーデータの収集と並行してそれより前に収集を
完了したグループのエコーデータについて周波数軸軸方
向のフーリエ変換および位相軸方向のフーリエ変換を行
ってグループデータを取得し、第１グループから第Ｍグ
ループのグループデータを加算する（位相補償演算し、
加算する）。これによれば、あるグループについてのデ
ータ収集処理と，その前にデータ収集したグループにつ
いての演算処理とを並列実行できるため、全データを収
集してから目的画像を得るまでの時間を短縮できるよう
になる。上記第３の観点（または第５の観点）によるＭ
ＲＩ装置では、グループデータを取得する毎にそれまで
に取得したグループデータを累積加算する（位相補償演
算し、累積加算する）。これにより、データ記憶領域を
節約できる。In the MRI apparatus according to the second viewpoint (or the fifth viewpoint), all echo data in the two-dimensional k-space are interleaved so that each group includes the same number of components in the phase axis direction. Then, the first group to the Mth group are grouped in the phase axis direction, the data is collected in order from the echo data of the first group, and the collection of the echo data of a certain group is completed in parallel and before that. For the echo data of the group, the Fourier transform in the frequency axis direction and the Fourier transform in the phase axis direction are performed to obtain the group data, and the group data of the first group to the Mth group are added (phase compensation calculation,
to add). According to this, since the data collection process for a certain group and the calculation process for the group for which the data has been collected can be executed in parallel, it is possible to shorten the time from collecting all the data to obtaining the target image. become. M from the third viewpoint (or the fifth viewpoint)
In the RI apparatus, every time group data is acquired, the group data acquired so far is cumulatively added (phase compensation calculation is performed and cumulative addition is performed). As a result, the data storage area can be saved.
【００１７】上記第４の観点（または第５の観点）によ
るＭＲＩ装置では、３次元ｋ−空間の全エコーデータ
を，スライス軸方向の各成分が各グループに同じ数だけ
含まれるようにインタリーブして，スライス軸方向に第
１グループから第Ｑグループまでにグループ分けし、第
１グループのエコーデータから順にデータを収集し、あ
るグループのエコーデータの収集と並行してそれより前
にエコーデータの収集を完了したグループのデータにつ
いてスライス軸方向のフーリエ変換を行ってグループデ
ータを取得し、第１グループから第Ｑグループのグルー
プデータを加算する（位相補償演算し、加算する）。こ
れによれば、スライス軸方向のエンコード量を先に変化
させるか後に変化させるかに拘らず、あるグループにつ
いてのデータ収集処理と，その前にデータ収集したグル
ープについての演算処理とを並列実行できるため、全デ
ータを収集してから目的画像を得るまでの時間を常に短
縮できるようになる。In the MRI apparatus according to the fourth aspect (or the fifth aspect), all echo data in the three-dimensional k-space are interleaved so that each group in the slice axis direction is included in the same number in each group. The first group to the Qth group are grouped in the slice axis direction, the data is collected in order from the echo data of the first group, and the echo data of the echo data of a certain group is acquired in parallel with the collection of the echo data of a certain group. The Fourier transform in the slice axis direction is performed on the data of the group for which collection has been completed to obtain group data, and the group data of the first group to the Qth group are added (phase compensation calculation is performed and added). According to this, regardless of whether the encoding amount in the slice axis direction is changed first or later, the data collection processing for a certain group and the calculation processing for the group for which data collection was performed can be executed in parallel. Therefore, it is possible to always reduce the time from collecting all the data to obtaining the target image.
【００１８】[0018]
【実施例】以下、図に示す実施例によりこの発明をさら
に詳しく説明する。なお、これによりこの発明が限定さ
れるものではない。図１は、この発明のＭＲＩ装置の一
実施例のブロック図である。このＭＲＩ装置１００にお
いて、マグネットアセンブリ１は、内部に被検体を挿入
するための空間部分（孔）を有し、この空間部分を取り
まくようにして、被検体に一定の静磁場を印加する静磁
場コイルと、勾配磁場を発生するための勾配磁場コイル
（勾配コイルは、スライス軸，周波数軸，位相軸のコイ
ルを備えている）と、被検体内の原子核のスピンを励起
するためのＲＦ（高周波）パルスを与える送信コイル
と、被検体からのＮＭＲ信号を検出する受信コイル等が
配置されている。静磁場コイル，勾配磁場コイル，送信
コイルおよび受信コイルは、それぞれ主磁場電源２，勾
配磁場駆動回路３，ＲＦ電力増幅器４および前置増幅器
５に接続されている。The present invention will be described in more detail with reference to the embodiments shown in the drawings. The present invention is not limited to this. FIG. 1 is a block diagram of an embodiment of the MRI apparatus of the present invention. In this MRI apparatus 100, the magnet assembly 1 has a space portion (hole) for inserting a subject therein, and a static magnetic field for applying a constant static magnetic field to the subject so as to surround this space portion. A coil, a gradient magnetic field coil for generating a gradient magnetic field (the gradient coil has a coil of a slice axis, a frequency axis, and a phase axis), and an RF (high frequency wave) for exciting spins of atomic nuclei in a subject. ) A transmitting coil that gives a pulse, a receiving coil that detects an NMR signal from the subject, and the like are arranged. The static magnetic field coil, the gradient magnetic field coil, the transmitting coil and the receiving coil are connected to the main magnetic field power source 2, the gradient magnetic field driving circuit 3, the RF power amplifier 4 and the preamplifier 5, respectively.
【００１９】シーケンス記憶回路８は、計算機７からの
指令に従い、記憶しているパルスシーケンスに基づいて
勾配磁場駆動回路３を操作し、前記マグネットアセンブ
リ１の勾配磁場コイルから勾配磁場を発生させると共
に、ゲート変調回路９を操作し、ＲＦ発振回路１０の高
周波出力信号を所定タイミング・所定包絡線のパルス状
信号に変調し、それをＲＦパルスとしてＲＦ電力増幅器
４に加え、ＲＦ電力増幅器４でパワー増幅した後、前記
マグネットアセンブリ１の送信コイルに印加し、目的の
領域を選択励起する。The sequence storage circuit 8 operates the gradient magnetic field drive circuit 3 based on the stored pulse sequence in accordance with a command from the computer 7 to generate a gradient magnetic field from the gradient magnetic field coil of the magnet assembly 1, and The gate modulation circuit 9 is operated to modulate the high frequency output signal of the RF oscillation circuit 10 into a pulsed signal having a predetermined timing and a predetermined envelope, which is applied as an RF pulse to the RF power amplifier 4 and amplified by the RF power amplifier 4. After that, it is applied to the transmission coil of the magnet assembly 1 to selectively excite a target region.
【００２０】前置増幅器５は、マグネットアセンブリ１
の受信コイルで検出された被検体からのＮＭＲ信号を増
幅し、位相検波器１２に入力する。位相検波器１２は、
ＲＦ発振回路１０の出力を参照信号とし、前置増幅器５
からのＮＭＲ信号を位相検波して、Ａ／Ｄ変換器１１に
与える。Ａ／Ｄ変換器１１は、位相検波後のアナログ信
号をディジタル信号に変換して、計算機７に入力する。
計算機７は、Ａ／Ｄ変換器１１からのデジタル信号に対
する画像再構成処理を行い、目的画像（前記目的の領域
のＮＭＲ画像）を生成する。この目的画像は、表示装置
６にて表示される。また、計算機７は、操作卓１３から
入力された情報を受け取るなどの全体的な制御を受け持
つ。The preamplifier 5 comprises a magnet assembly 1
The NMR signal from the subject detected by the receiving coil is amplified and input to the phase detector 12. The phase detector 12 is
The output of the RF oscillation circuit 10 is used as a reference signal, and the preamplifier 5
The NMR signal from is phase-detected and given to the A / D converter 11. The A / D converter 11 converts the analog signal after the phase detection into a digital signal and inputs the digital signal to the computer 7.
The computer 7 performs image reconstruction processing on the digital signal from the A / D converter 11 to generate a target image (an NMR image of the target region). This target image is displayed on the display device 6. Further, the computer 7 is responsible for overall control such as receiving information input from the console 13.
【００２１】図２は、ＭＲＩ装置１００における２次元
の画像再構成処理を示すフローチャートである。ステッ
プＡ１では、操作者がパルスシーケンスのパラメータを
設定する。このパラメータのうちの位相軸方向のＦＯ
Ｖ，エンコード数Ｎ，勾配パルスの大きさＧ，勾配パル
スの幅Ｔによって、 Δθ＝２π／（Ｇ・Ｔ・ＦＯＶ・Ｎ） から位相軸方向のエンコードステップ量Δθが求められ
る。ステップＡ２では、グループ番号カウンタｍ＝１と
する。ステップＡ３では、グループ内データ番号カウン
タｋ’＝１とする。FIG. 2 is a flowchart showing a two-dimensional image reconstruction process in the MRI apparatus 100. In step A1, the operator sets the pulse sequence parameters. FO in the phase axis direction of these parameters
The encoding step amount Δθ in the phase axis direction is obtained from Δθ = 2π / (G · T · FOV · N) by V, the number of encodes N, the gradient pulse size G, and the gradient pulse width T. At step A2, the group number counter m = 1. In step A3, the in-group data number counter k '= 1 is set.
【００２２】ステップＡ４では、ｍ番目のグループの
ｋ’番目のエコーデータＲ（m,k'）を収集するための位
相軸方向のエンコード量θ（m,k'）を、 θ（m,k'）＝−π／（Ｇ・Ｔ・ＦＯＶ）＋ｍ・Δθ＋
（ｋ’−１）・Ｍ・Δθ により算出する。但し、Ｍは、グループ数である。ま
た、各成分が各グループに同じ数だけ含まれるようにイ
ンタリーブしてグループ分けを行なっている。ステップ
Ａ５では、位相軸方向のエンコード量θ（m,k'）で、エ
コーデータＲ（m,k'）を収集する。ステップＡ６では、
グループ内データ番号カウンタｋ’を調べて、ｍ番目の
グループの全てのエコーデータを収集したかチェックす
る。そして、ｍ番目のグループの全てのエコーデータを
未だ収集していないならステップＡ７に進み、ｍ番目の
グループの全てのエコーデータを収集したならステップ
Ａ８に進む。ステップＡ７では、ｋ’＝ｋ’＋１とし、
前記ステップＡ４に戻る。ステップＡ８では、グループ
番号カウンタｍ＝Ｍか否かチェックする。そして、ｍ＝
ＭでないならステップＡ９に進み、ｍ＝Ｍならデータ収
集処理を終了する。ステップＡ９では、ｍ＝ｍ＋１と
し、前記ステップＡ３に戻る。In step A4, the encoding amount θ (m, k ') in the phase axis direction for collecting the k'th echo data R (m, k') of the mth group is set to θ (m, k). ') =-Π / (G ・ T ・ FOV) + m ・ Δθ +
It is calculated by (k′−1) · M · Δθ. However, M is the number of groups. In addition, interleaving is performed so that each component is included in the same number in each group. In step A5, the echo data R (m, k ') is collected with the encoding amount θ (m, k') in the phase axis direction. In step A6,
The in-group data number counter k'is checked to see if all echo data for the mth group has been collected. Then, if all echo data of the m-th group have not been collected yet, the process proceeds to step A7, and if all echo data of the m-th group has been collected, the process proceeds to step A8. In step A7, k '= k' + 1 is set,
The procedure returns to step A4. At step A8, it is checked whether or not the group number counter m = M. And m =
If it is not M, the process proceeds to step A9, and if m = M, the data collection process ends. In step A9, m = m + 1 is set, and the process returns to step A3.
【００２３】上記ステップＡ３〜Ａ７により、図３に示
す２次元ｋ−空間の第ｍグループのエコーデータＲ（m,
1）〜Ｒ（m,N/M）を順に収集するが、一つのエコーデー
タを収集している間に、それより前に収集したエコーデ
ータについての周波数軸方向の１次元フーリエ変換をス
テップＢ１で並行して行って中間データｆ（m,k'）を算
出し、その中間データｆ（m,k'）を保存しておく。ま
た、上記ステップＡ８，Ａ９により、グループ毎のエコ
ーデータ収集を進めていくが、あるグループのエコーデ
ータの収集と並行して、それより前に全ての中間データ
を算出したグループの中間データについてステップＢ２
で位相軸方向の１次元フーリエ変換を行い、ステップＢ
３で位相補償を行い、グループデータを取得する。そし
て、取得したグループデータをステップＢ４でメモリ上
の画像エリアに累積加算すると共に、その画像エリアの
データに基づいて画像を表示する。By the above steps A3 to A7, the echo data R (m, m of the m-th group in the two-dimensional k-space shown in FIG. 3 is obtained.
1) to R (m, N / M) are sequentially collected, and while collecting one echo data, the one-dimensional Fourier transform in the frequency axis direction of the echo data collected before is collected in step B1. In parallel, the intermediate data f (m, k ') is calculated, and the intermediate data f (m, k') is stored. Further, although the echo data collection for each group is advanced by the above steps A8 and A9, in parallel with the collection of the echo data of a certain group, the step for the intermediate data of the group for which all the intermediate data are calculated before that is performed. B2
Performs one-dimensional Fourier transform in the phase axis direction at step B
Phase compensation is performed in 3 to obtain group data. Then, in step B4, the acquired group data is cumulatively added to the image area on the memory, and an image is displayed based on the data in the image area.
【００２４】以上の画像再構成処理によれば、あるグル
ープについてのデータ収集処理（ステップＡ３〜Ａ９）
と，その前にデータ収集したグループについての位相軸
方向の演算処理（ステップＢ２）とを並列実行できるた
め、目的画像を得られるまでの時間を短縮できる。According to the above image reconstruction processing, data collection processing for a certain group (steps A3 to A9)
Since the calculation processing in the phase axis direction (step B2) for the group for which data has been collected before that can be executed in parallel, the time until the target image is obtained can be shortened.
【００２５】ここで、フーリエ変換の演算数について説
明する。従来の画像再構成処理では、エンコード数Ｎが
２のべき乗の時には、１ステップにＮ回の複素積算が発
生し、これがlog2{Ｎ}ステップ繰り返されるので、計Ｎ
log2{Ｎ}回の複素積算が行なわれる。このＮlog2{Ｎ}回
の複素積算は、データ収集完了後に全て行なわれる（図
１２のステップＦ２）。これに対して、上記実施例の画
像再構成処理では、位相エンコード数ＮをＮ／Ｍ個ごと
に１つのグループとしＭ個のグループに分けるので、１
グループについて（Ｎ／Ｍ）log2{Ｎ／Ｍ}回の複素積算
になるが、データ収集完了までに第（Ｍ−１）グループ
までについてのフーリエ変換はデータ収集と並行して実
行しているので（図２のステップＢ２）、データ収集完
了後には、第Ｍグループについての（Ｎ／Ｍ）log2{Ｎ
／Ｍ}回の複素積算だけでよい。但し、グループデータ
を整合させるために位相補償を行なう必要があり、１グ
ループの位相補償にはＮ回の複素積算になるが、データ
収集完了までに第（Ｍ−１）グループまでについての位
相補償はデータ収集と並行して実行しているので（図２
のステップＢ３）、データ収集完了後には、第Ｍグルー
プについてのＮ回の位相補償だけでよい。従って、計
（Ｎ／Ｍ）log2{Ｎ／Ｍ}＋Ｎ回の複素積算がデータ収集
完了後に行なわれることになる。結局、上記実施例は、
従来より、Ｎlog2{Ｎ}−（Ｎ／Ｍ）log2{Ｎ／Ｍ}＋Ｎ回
の複素積算の時間だけ、目的画像を得られるまでの待ち
時間を短縮できることとなる。この式では、グループ数
Ｍを大きくすれば短縮効果が大きくなる。しかし、グル
ープ数Ｍを大きくすることにより生じるオーバヘッド
や、計算機７の演算速度やエンコード数Ｎなどによりデ
ータ収集完了までに第（Ｍ−１）グループまでについて
の演算処理を完了できない場合も考えられるので、Ｍ＝
２，３，…について試行してみて最適の分割数Ｍを決定
するのが現実的である。Here, the number of Fourier transform operations will be described. In the conventional image reconstruction processing, when the number of encodes N is a power of 2, N times of complex integrations occur in one step, and this is repeated log 2 {N} steps.
Complex integration is performed log 2 {N} times. This Nlog 2 {N} times of complex integration is all performed after the data collection is completed (step F2 in FIG. 12). On the other hand, in the image reconstruction processing of the above-described embodiment, the number N of phase encodes is set as one group for each N / M and divided into M groups.
For the group, (N / M) log 2 {N / M} times of complex integration will be performed, but by the time the data collection is completed, the Fourier transform up to the (M-1) th group is executed in parallel with the data collection. Therefore (step B2 in FIG. 2), after the data collection is completed, (N / M) log 2 {N
Only / M} complex integrations are required. However, it is necessary to perform phase compensation in order to match the group data, and the phase compensation for one group requires N times of complex integration, but the phase compensation for the (M-1) th group is not completed until the data collection is completed. Is running in parallel with data collection (see Figure 2
Step B3), after the data collection is completed, only N phase compensations for the Mth group are required. Therefore, a total of (N / M) log 2 {N / M} + N complex integrations will be performed after the completion of data collection. After all, the above example
Conventionally, the waiting time until the target image is obtained can be shortened by the time of Nlog 2 {N} − (N / M) log 2 {N / M} + N times of complex integration. In this formula, the larger the number M of groups, the greater the shortening effect. However, it is possible that the calculation processing for the (M-1) th group cannot be completed before the data collection is completed due to the overhead caused by increasing the number of groups M, the calculation speed of the computer 7, the number of encodings N, and the like. , M =
It is realistic to try 2, 3, ... To determine the optimal number of divisions M.
【００２６】図４は、ＭＲＩ装置１００における２次元
の画像再構成処理の他の例を示すフローチャートであ
る。ステップＡ１’では、操作者がパルスシーケンスの
パラメータを設定する。このパラメータのうちの位相軸
方向のＦＯＶ，エンコード数Ｎ，勾配パルスの大きさ
Ｇ，勾配パルスの幅Ｔによって、位相軸方向のエンコー
ドステップ量Δθが求められる。ステップＡ２’では、
グループ番号カウンタｍ＝１とする。ステップＡ３’で
は、グループ内データ番号カウンタｋ’＝１とする。FIG. 4 is a flowchart showing another example of the two-dimensional image reconstruction processing in the MRI apparatus 100. In step A1 ', the operator sets the parameters of the pulse sequence. Among these parameters, the encoding step amount Δθ in the phase axis direction is obtained from the FOV in the phase axis direction, the number of encodes N, the gradient pulse size G, and the gradient pulse width T. In step A2 ',
The group number counter m = 1. At step A3 ′, the in-group data number counter k ′ = 1 is set.
【００２７】ステップＡ４’では、ｍ番目のグループの
ｋ’番目のエコーデータＲ（m,k'）を収集するための位
相軸方向のエンコード量θ（m,k'）を、 θ（m,k'）＝−π／（Ｇ・Ｔ・ＦＯＶ）＋ｍ・Δθ＋
（ｋ’−１）・Ｍ・Δθ により算出する。但し、Ｍは、グループ数である。ステ
ップＡ５’では、位相軸方向のエンコード量θ（m,k'）
で、エコーデータＲ（m,k'）を収集する。ステップＡ
６’では、グループ内データ番号カウンタｋ’を調べ
て、ｍ番目のグループの全てのエコーデータを収集した
かチェックする。そして、ｍ番目のグループの全てのエ
コーデータを未だ収集していないなら、ステップＡ７’
に進む。ｍ番目のグループの全てのエコーデータを収集
したなら、ステップＡ８’に進む。ステップＡ７’で
は、ｋ’＝ｋ’＋１とし、前記ステップＡ４’に戻る。
ステップＡ８’では、グループ番号カウンタｍ＝Ｍか否
かチェックする。そして、ｍ＝ＭでないならステップＡ
９’に進み、ｍ＝Ｍならデータ収集処理を終了する。ス
テップＡ９’では、ｍ＝ｍ＋１とし、前記ステップＡ
３’に戻る。At step A4 ', the encode amount θ (m, k') in the phase axis direction for collecting the k'th echo data R (m, k ') of the mth group is set to θ (m, k ′) = − π / (G · T · FOV) + m · Δθ +
It is calculated by (k′−1) · M · Δθ. However, M is the number of groups. At step A5 ′, the encoding amount θ (m, k ′) in the phase axis direction
Then, the echo data R (m, k ') is collected. Step A
At 6 ', the in-group data number counter k'is checked to see if all echo data of the m-th group has been collected. Then, if all echo data of the m-th group have not been collected yet, step A7 '
Proceed to. When all the echo data of the m-th group have been collected, the process proceeds to step A8 '. In step A7 ', k' = k '+ 1 is set, and the process returns to step A4'.
At step A8 ', it is checked whether or not the group number counter m = M. Then, if m = M is not satisfied, step A
9 ', and if m = M, the data collection process ends. In step A9 ′, m = m + 1 is set, and the step A
Return to 3 '.
【００２８】上記ステップＡ３’〜Ａ９’により、図５
に示す２次元ｋ−空間の各グループのエコーデータを順
に収集するが、あるグループのエコーデータの収集と並
行して、それより前に全ての収集を完了したグループの
エコーデータについてステップＢ１’で周波数軸および
位相軸方向にフーリエ変換を行い、ステップＢ２’で位
相補償を行い、取得したグループデータをステップＢ
３’でメモリ上の画像エリアに累積加算する。そして、
全てのエコーデータの収集を完了すると、そのときの画
像エリアのデータに基づいて画像を表示する。By the above steps A3 'to A9', FIG.
The echo data of each group in the two-dimensional k-space shown in FIG. 2 is sequentially collected. In parallel with the collection of the echo data of a certain group, the echo data of the group for which all the acquisitions are completed before that are performed in step B1 ′. Fourier transformation is performed in the frequency axis and phase axis directions, phase compensation is performed in step B2 ', and the acquired group data is processed in step B.
At 3 ', cumulative addition is made to the image area on the memory. And
When collection of all echo data is completed, an image is displayed based on the data in the image area at that time.
【００２９】以上の画像再構成処理によれば、あるグル
ープについてのデータ収集処理（ステップＡ３’〜Ａ
９’）と，その前にデータ収集したグループについての
周波数軸および位相軸方向の演算処理（ステップＢ１’
〜Ｂ３’）とを並列実行できるため、目的画像を得られ
るまでの時間を短縮できる。According to the above image reconstruction processing, data collection processing (steps A3'-A for a certain group) is performed.
9 ') and the calculation processing in the frequency axis and phase axis directions for the group for which data was collected before (step B1'
~ B3 ') can be executed in parallel, the time required to obtain the target image can be shortened.
【００３０】図６は、ＭＲＩ装置１００における３次元
の画像再構成処理を示すフローチャートである。ステッ
プＣ１では、操作者がパルスシーケンスのパラメータを
設定する。このパラメータのうちの位相軸方向のＦＯＶ
k ，エンコード数Ｎ，勾配パルスの大きさＧk，勾配パ
ルスの幅Ｔkによって、 Δθ＝２π／（Ｇk・Ｔk・ＦＯＶk・Ｎ） から位相軸方向のエンコードステップ量Δθが求められ
る。又、スライス軸方向のＦＯＶs ，エンコード数Ｓ，
勾配パルスの大きさＧs，勾配パルスの幅Ｔsによって、 Δφ＝２π／（Ｇs・Ｔs・ＦＯＶs・Ｓ） からスライス軸方向のエンコードステップ量Δφが求め
られる。ステップＣ２では、位相軸方向グループ番号カ
ウンタｍ＝１とする。ステップＣ３では、位相軸方向グ
ループ内データ番号カウンタｋ’＝１とする。ステップ
Ｃ４では、位相軸方向のｍ番目のグループのｋ’番目の
エンコード量θ（m,k'）を、 θ（m,k'）＝−π／（Ｇk・Ｔk・ＦＯＶk）＋ｍ・Δθ
＋（ｋ’−１）・Ｍ・Δθ により算出する。但し、Ｍは、位相軸方向のグループ数
である。また、各成分が各グループに同じ数だけ含まれ
るようにインタリーブしてグループ分けを行なってい
る。FIG. 6 is a flow chart showing a three-dimensional image reconstruction process in the MRI apparatus 100. In step C1, the operator sets the pulse sequence parameters. FOV in the phase axis direction of this parameter
The encoding step amount Δθ in the phase axis direction is obtained from Δθ = 2π / (Gk · Tk · FOVk · N) by k, the encoding number N, the gradient pulse size Gk, and the gradient pulse width Tk. Also, FOVs in the slice axis direction, the number of encodes S,
The encoding step amount Δφ in the slice axis direction is obtained from Δφ = 2π / (Gs · Ts · FOVs · S) by the gradient pulse size Gs and the gradient pulse width Ts. In step C2, the phase axis direction group number counter m = 1. In step C3, the data number counter in the phase axis direction group k '= 1. In step C4, the k′-th encoding amount θ (m, k ′) of the m-th group in the phase axis direction is θ (m, k ′) = − π / (Gk · Tk · FOVk) + m · Δθ
It is calculated by + (k′−1) · M · Δθ. However, M is the number of groups in the phase axis direction. In addition, interleaving is performed so that each component is included in the same number in each group.
【００３１】ステップＣ５では、スライス軸方向グルー
プ番号カウンタｑ＝１とする。ステップＣ６では、スラ
イス軸方向データ番号カウンタｓ’＝１とする。ステッ
プＣ７では、スライス軸方向のｑ番目のグループのｓ’
番目のエンコード量φ（q,s'）を、 φ（q,s'）＝−π／（Ｇs・Ｔs・ＦＯＶs）＋ｑ・Δθ
＋（ｓ’−１）・Ｑ・Δφ により算出する。但し、Ｑは、スライス軸方向のグルー
プ数である。また、各成分が各グループに同じ数だけ含
まれるようにインタリーブしてグループ分けを行なって
いる。ステップＣ８では、位相軸方向のエンコード量θ
（m,k'），スライス軸方向のエンコード量φ（q,s'）
で、エコーデータＲ{q,s'}(m,k')を収集する。ステップ
Ｃ９では、スライス軸方向グループ内データ番号カウン
タｓ’を調べて、位相軸方向のエンコード量θ（m,k'）
に対するｑ番目のスライス軸方向グループの全てのエコ
ーデータを収集したかチェックする。そして、全てのエ
コーデータを未だ収集していないならステップＣ１０に
進み、全てのエコーデータを収集したならステップＣ１
１に進む。ステップＣ１０では、ｓ’＝ｓ’＋１とし、
前記ステップＣ７に戻る。ステップＣ１１では、スライ
ス軸方向グループ番号カウンタｑ＝Ｑか否かチェックす
る。そして、ｑ＝ＱでないならステップＣ１２に進み、
ｑ＝ＱならステップＣ１３に進む。ステップＣ１２で
は、ｑ＝ｑ＋１とし、前記ステップＣ６に戻る。At step C5, the slice axis direction group number counter q = 1. At step C6, the slice axis direction data number counter s' = 1 is set. At step C7, s'of the q-th group in the slice axis direction
The th encoding amount φ (q, s ') is φ (q, s') =-π / (Gs · Ts · FOVs) + q · Δθ
It is calculated by + (s′−1) · Q · Δφ. However, Q is the number of groups in the slice axis direction. In addition, interleaving is performed so that each component is included in the same number in each group. At step C8, the encode amount θ in the phase axis direction
(M, k '), encoding amount in slice axis direction φ (q, s')
Then, the echo data R {q, s '} (m, k') is collected. In step C9, the slice axis direction group data number counter s'is checked to determine the encoding amount θ (m, k ') in the phase axis direction.
Check whether all echo data of the q-th slice axial group with respect to are collected. Then, if all echo data have not been collected yet, the process proceeds to step C10, and if all echo data has been collected, step C1.
Go to 1. In step C10, s '= s' + 1 is set,
The procedure returns to step C7. At Step C11, it is checked whether or not the slice axis direction group number counter q = Q. If not q = Q, the process proceeds to step C12,
If q = Q, the process proceeds to step C13. In step C12, q = q + 1 is set, and the process returns to step C6.
【００３２】ステップＣ１３では、位相軸方向グループ
内データ番号カウンタｋ’を調べて、ｍ番目のグループ
の全てのエコーデータを収集したかチェックする。そし
て、ｍ番目のグループの全てのエコーデータを未だ収集
していないなら、ステップＣ１４に進む。ｍ番目のグル
ープの全てのエコーデータを収集したなら、ステップＣ
１５に進む。ステップＣ１４では、ｋ’＝ｋ’＋１と
し、前記ステップＣ４に戻る。ステップＣ１５では、位
相軸方向グループ番号カウンタｍ＝Ｍか否かチェックす
る。そして、ｍ＝ＭでないならステップＣ１６に進み、
ｍ＝Ｍならデータ収集処理を終了する。ステップＣ１６
では、ｍ＝ｍ＋１とし、前記ステップＣ３に戻る。In step C13, the in-group data number counter k'in the phase axis direction is checked to see if all echo data of the m-th group have been collected. Then, if all the echo data of the m-th group have not been collected yet, the process proceeds to step C14. Once all echo data for the mth group has been collected, step C
Proceed to 15. In step C14, k '= k' + 1 is set, and the process returns to step C4. At Step C15, it is checked whether or not the phase axis direction group number counter m = M. If not m = M, the process proceeds to step C16,
If m = M, the data collection process ends. Step C16
Then, m = m + 1 is set, and the process returns to the step C3.
【００３３】上記ステップＣ６〜Ｃ１０により、例えば
図８に示す３次元ｋ−空間のエコーデータＲ{q,1}(m,
k')〜Ｒ{q,2}(m,k')を順に収集するが、一つのエコーデ
ータを収集している間に、それより前に収集したエコー
データについての周波数軸方向の１次元フーリエ変換を
ステップＤ１で並行して行って中間データを算出し、そ
の中間データを保存しておく。また、上記ステップＣ１
１，Ｃ１２により、スライス軸方向グループ毎のエコー
データ収集を進めていくが、あるグループのエコーデー
タの収集と並行して、それより前に全ての中間データを
算出したグループの中間データについてステップＤ２で
スライス軸方向の１次元フーリエ変換を行い、ステップ
Ｄ３で位相補償を行い、取得したスライス軸中間データ
をメモリ上の中間データエリアに累積加算しておく。ま
た、上記ステップＣ１３〜Ｃ１６により、位相軸方向グ
ループ毎のエコーデータ収集を進めていくが、あるグル
ープのエコーデータの収集と並行して、それより前に全
てのスライス軸中間データを算出したグループのスライ
ス軸中間データについてステップＤ４で位相軸方向の１
次元フーリエ変換を行い、ステップＤ５で位相補償を行
い、取得したグループデータをメモリ上の画像エリアに
累積加算する。また、その画像エリアのデータに基づい
て画像を表示する。By the steps C6 to C10, the echo data R {q, 1} (m,
k ′) to R {q, 2} (m, k ′) are collected in order, but while collecting one echo data, one-dimensional in the frequency axis direction for the echo data collected before that. Fourier transform is performed in parallel in step D1 to calculate intermediate data, and the intermediate data is stored. Also, the above step C1
1 and C12, the echo data collection for each slice axis direction group is proceeded, but in parallel with the collection of the echo data of a certain group, the intermediate data of the group in which all the intermediate data are calculated before that step D2 In step D3, one-dimensional Fourier transform in the slice axis direction is performed, phase compensation is performed, and the acquired slice axis intermediate data is cumulatively added to the intermediate data area on the memory. Further, although the echo data collection for each phase axis direction group is advanced by the steps C13 to C16, in parallel with the collection of the echo data of a certain group, all slice axis intermediate data are calculated before that group. Slice slice intermediate data of 1 in the phase axis direction in step D4.
A dimensional Fourier transform is performed, phase compensation is performed in step D5, and the acquired group data is cumulatively added to the image area on the memory. Further, the image is displayed based on the data in the image area.
【００３４】以上の画像再構成処理によれば、あるスラ
イス軸方向グループについてのデータ収集処理（ステッ
プＣ７〜Ｃ１２）と，その前にデータ収集したグループ
についてのスライス軸方向の演算処理（ステップＤ２）
とを並列実行できる。また、ある位相軸方向グループに
ついてのデータ収集処理（ステップＣ１３〜Ｃ１６）
と，その前にデータ収集したグループについての位相軸
方向の演算処理（ステップＤ４）とを並列実行できる。
このため、目的画像を得られるまでの時間を短縮でき
る。According to the above image reconstruction processing, data collection processing for a certain slice axis direction group (steps C7 to C12) and slice axis direction calculation processing for the group for which data collection was performed before that (step D2).
And can be executed in parallel. Further, data collection processing for a certain phase axis direction group (steps C13 to C16)
And the arithmetic processing in the phase axis direction (step D4) for the group for which data collection was performed can be executed in parallel.
Therefore, the time until the target image is obtained can be shortened.
【００３５】図９は、ＭＲＩ装置１００における３次元
の画像再構成処理の他の例を示すフローチャートであ
る。ステップＣ１’では、操作者がパルスシーケンスの
パラメータを設定する。このパラメータのうちの位相軸
方向のＦＯＶk ，エンコード数Ｎ，勾配パルスの大きさ
Ｇk，勾配パルスの幅Ｔkによって位相軸方向のエンコー
ドステップ量Δθが求められる。又、スライス軸方向の
ＦＯＶs ，エンコード数Ｓ，勾配パルスの大きさＧs，
勾配パルスの幅Ｔsによってスライス軸方向のエンコー
ドステップ量Δφが求められる。ステップＣ２’では、
スライス軸方向グループ番号カウンタｑ＝１に設定す
る。ステップＣ３’では、スライス軸方向データ番号カ
ウンタｓ’＝１に設定する。ステップＣ４’では、スラ
イス軸方向のｑ番目のグループのｓ’番目のエンコード
量φ（q,s'）を、 φ（q,s'）＝−π／（Ｇs・Ｔs・ＦＯＶs）＋ｑ・Δθ
＋（ｓ’−１）・Ｑ・Δφ により算出する。但し、Ｑは、スライス軸方向のグルー
プ数である。ステップＣ５’では、位相軸方向データ番
号カウンタｋ＝１とする。ステップＣ６’では、位相軸
方向のｋ番目のエンコード量θ（ｋ）を、 θ（ｋ）＝−π／（Ｇk・Ｔk・ＦＯＶk）＋ｋ・Δθ により算出する。なお、位相軸方向にはグループ分けし
ていない。FIG. 9 is a flowchart showing another example of the three-dimensional image reconstruction processing in the MRI apparatus 100. In step C1 ', the operator sets the parameters of the pulse sequence. Of these parameters, the FOVk in the phase axis direction, the number of encodes N, the gradient pulse size Gk, and the gradient pulse width Tk determine the encode step amount Δθ in the phase axis direction. Further, FOVs in the slice axis direction, the number of encodes S, the magnitude Gs of the gradient pulse,
The encoding step amount Δφ in the slice axis direction is obtained from the width Ts of the gradient pulse. In step C2 ',
Slice axis direction group number counter q = 1 is set. At step C3 ', the slice axis direction data number counter s' = 1 is set. In step C4 ′, the s′-th encoding amount φ (q, s ′) of the q-th group in the slice axis direction is expressed as φ (q, s ′) = − π / (Gs · Ts · FOVs) + q · Δθ
It is calculated by + (s′−1) · Q · Δφ. However, Q is the number of groups in the slice axis direction. In step C5 ′, the phase axis direction data number counter k = 1 is set. At step C6 ′, the kth encoding amount θ (k) in the phase axis direction is calculated by θ (k) = − π / (Gk · Tk · FOVk) + k · Δθ. Note that the groups are not divided in the phase axis direction.
【００３６】ステップＣ７’では、位相軸方向のエンコ
ード量θ（ｋ），スライス軸方向のエンコード量φ（q,
s'）で、エコーデータＲ{q,s'}(ｋ)を収集する。ステッ
プＣ８’では、位相軸方向データ番号カウンタｋを調べ
て、ｋ＝Ｎかチェックする。そして、ｋ＝Ｎでないなら
ステップＣ９’に進み、ｋ＝ＮならステップＣ１０’に
進む。ステップＣ９’では、ｋ＝ｋ＋１とし、前記ステ
ップＣ６’に戻る。ステップＣ１０’では、スライス軸
方向グループ内データ番号カウンタｓ’を調べて、ｑ番
目のスライス軸方向グループの全てのエコーデータを収
集したかチェックする。そして、全てのエコーデータを
未だ収集していないなら、ステップＣ１１’に進む。ま
た、全てのエコーデータを収集したなら、ステップＣ１
２’に進む。ステップＣ１１’では、ｓ'＝ｓ'＋１と
し、前記ステップＣ４’に戻る。ステップＣ１２’で
は、スライス軸方向グループ番号カウンタｑ＝Ｑか否か
チェックする。そして、ｑ＝ＱでないならステップＣ１
３’に進み、ｑ＝Ｑならデータ収集処理を終了する。ス
テップＣ１３’では、ｑ＝ｑ＋１とし、前記ステップＣ
３’に戻る。At step C7 ', the encoding amount θ (k) in the phase axis direction and the encoding amount φ (q, q
In s '), echo data R {q, s'} (k) is collected. In step C8 ', the phase axis direction data number counter k is checked to see if k = N. If not k = N, the process proceeds to step C9 ', and if k = N, the process proceeds to step C10'. In step C9 ', k = k + 1 is set, and the process returns to step C6'. In step C10 ', the slice-axis-direction in-group data number counter s' is checked to see if all echo data of the q-th slice-axis-direction group have been collected. Then, if all the echo data have not been collected yet, the process proceeds to step C11 '. If all echo data has been collected, step C1
Go to 2 '. In step C11 ', s' = s '+ 1 is set, and the process returns to step C4'. At Step C12 ', it is checked whether or not the slice axis direction group number counter q = Q. If q = Q, step C1
The process proceeds to 3 ', and if q = Q, the data collection process ends. In step C13 ′, q = q + 1 is set, and the above step C
Return to 3 '.
【００３７】上記ステップＣ３’〜Ｃ１３’により、ス
ライス軸方向グループ毎のエコーデータ収集を進めてい
くが、あるグループのエコーデータの収集と並行して、
それより前に全ての中間データを算出したグループの中
間データについてステップＤ１’でスライス軸方向の１
次元フーリエ変換を行い、ステップＤ２’で位相補償を
行い、取得したスライス軸中間データをメモリ上の画像
エリアに累積加算しておく。そして、図１１に示す全て
のエコーデータの収集を完了したら、周波数軸および位
相軸方向のフーリエ変換をステップＤ３’で行う。以上
のステッＣ１’〜Ｃ１３’のデータ収集処理およびステ
ップＤ１’〜Ｄ３’の演算処理によって、各スライスに
ついての目的画像が再構成される。By the above steps C3 'to C13', the echo data collection for each group in the slice axis direction is proceeded, but in parallel with the collection of the echo data of a certain group,
For the intermediate data of the group for which all the intermediate data have been calculated before that, 1 in the slice axis direction is set in step D1 ′.
Dimensional Fourier transform is performed, phase compensation is performed in step D2 ', and the acquired slice axis intermediate data is cumulatively added to the image area on the memory. Then, when collection of all the echo data shown in FIG. 11 is completed, Fourier transform in the frequency axis and phase axis directions is performed in step D3 ′. The target image for each slice is reconstructed by the data collection processing of steps C1 ′ to C13 ′ and the arithmetic processing of steps D1 ′ to D3 ′.
【００３８】以上の画像再構成処理によれば、あるスラ
イス軸方向グループについてのデータ収集処理と，その
前にデータ収集したグループについてのスライス軸方向
の演算処理とを並列実行できるため、目的画像を得られ
るまでの時間を短縮することが出来る。According to the above image reconstruction processing, the data collection processing for a certain slice axis direction group and the calculation processing in the slice axis direction for the group for which data collection was performed before that can be executed in parallel. It is possible to shorten the time until it is obtained.
【００３９】[0039]
【発明の効果】この発明のＭＲＩ装置における画像再構
成方法およびＭＲＩ装置によれば、ｋ−空間のエコーデ
ータをグループに分けて収集し、あるグループについて
のデータ収集処理と，その前にデータ収集したグループ
についての演算処理とを並列実行するため、目的画像を
得るまでの時間を短縮することが出来る。According to the image reconstruction method and the MRI apparatus in the MRI apparatus of the present invention, echo data of k-space is divided into groups and collected, and data collection processing for a certain group is performed and data collection is performed before that. Since the calculation processing for the selected group is executed in parallel, the time required to obtain the target image can be shortened.
【図１】この発明のＭＲＩ装置の一実施例のブロック図
である。FIG. 1 is a block diagram of an embodiment of an MRI apparatus of the present invention.
【図２】図１のＭＲＩ装置における２次元の画像再構成
処理を示すフローチャートである。FIG. 2 is a flowchart showing a two-dimensional image reconstruction process in the MRI apparatus of FIG.
【図３】インタリーブしてグループ分けした２次元ｋ−
空間の概念図である。[Fig. 3] Two-dimensional k- interleaved and grouped
It is a conceptual diagram of space.
【図４】図１のＭＲＩ装置における２次元の画像再構成
処理の他の例を示すフローチャートである。FIG. 4 is a flowchart showing another example of a two-dimensional image reconstruction process in the MRI apparatus of FIG.
【図５】インタリーブしないでグループ分けした２次元
ｋ−空間の概念図である。FIG. 5 is a conceptual diagram of a two-dimensional k-space divided into groups without interleaving.
【図６】図１のＭＲＩ装置における３次元の画像再構成
処理を示すフローチャートである。6 is a flowchart showing a three-dimensional image reconstruction process in the MRI apparatus of FIG.
【図７】図６の続きのフローチャートである。FIG. 7 is a continuation of the flowchart of FIG.
【図８】インタリーブしてグループ分けした３次元ｋ−
空間の概念図である。FIG. 8 is a three-dimensional k- grouped by interleaving.
It is a conceptual diagram of space.
【図９】図１のＭＲＩ装置における３次元の画像再構成
処理の他の例を示すフローチャートである。9 is a flowchart showing another example of a three-dimensional image reconstruction process in the MRI apparatus in FIG.
【図１０】図９の続きのフローチャートである。FIG. 10 is a flowchart continued from FIG. 9;
【図１１】インタリーブしないでグループ分けした３次
元ｋ−空間の概念図である。FIG. 11 is a conceptual diagram of a three-dimensional k-space divided into groups without interleaving.
【図１２】従来のＭＲＩ装置における２次元の画像再構
成処理を示すフローチャートである。FIG. 12 is a flowchart showing a two-dimensional image reconstruction process in the conventional MRI apparatus.
【図１３】２次元ｋ−空間の概念図である。FIG. 13 is a conceptual diagram of a two-dimensional k-space.
【図１４】従来のＭＲＩ装置における３次元の画像再構
成処理を示すフローチャートである。FIG. 14 is a flowchart showing a three-dimensional image reconstruction process in the conventional MRI apparatus.
【図１５】３次元ｋ−空間の概念図である。FIG. 15 is a conceptual diagram of a three-dimensional k-space.
【図１６】従来のＭＲＩ装置における３次元の画像再構
成処理の他の例を示すフローチャートである。FIG. 16 is a flowchart showing another example of a three-dimensional image reconstruction process in the conventional MRI apparatus.
１００ ＭＲＩ装置 １ マグネットアセンブリ ３ 勾配磁場駆動回路 ７ 計算機 ８ シーケンス記憶回路 ９ ゲート変調回路 １０ ＲＦ発振回路 １３ 操作卓 100 MRI apparatus 1 Magnet assembly 3 Gradient magnetic field drive circuit 7 Computer 8 Sequence storage circuit 9 Gate modulation circuit 10 RF oscillation circuit 13 Operation console
Claims (5)
ループ分けし、そのエンコード方向のエンコード量を変
化させて第１のグループから順にエコーデータを収集す
ると共に、あるグループのエコーデータの収集と並行し
てそれより前にエコーデータの収集を完了したグループ
について前記エンコード方向のフーリエ変換を行ってグ
ループデータを取得し、全てのグループのグループデー
タを加算して、全データについての前記エンコード方向
のフーリエ変換結果とすることを特徴とするＭＲＩ装置
における画像再構成方法。1. Echo data is divided into groups in a certain encoding direction, the encoding amount in the encoding direction is changed to collect the echo data in order from the first group, and in parallel with the collection of the echo data of a certain group. Performing Fourier transform in the encoding direction on the group for which echo data collection is completed before that to obtain group data, adding group data of all groups, and Fourier transform result of the encoding direction for all data An image reconstructing method in an MRI apparatus, comprising:
相軸方向の各成分が各グループに同じ数だけ含まれるよ
うにインタリーブして，位相軸方向に第１グループから
第Ｍグループまでにグループ分けして第１グループのエ
コーデータから順にデータを収集するグループ分けデー
タ収集手段と、あるグループのエコーデータの収集と並
行してそれより前に収集を完了したグループのエコーデ
ータについて周波数軸軸方向のフーリエ変換および位相
軸方向のフーリエ変換を行ってグループデータを取得す
るグループデータ取得手段と、第１グループから第Ｍグ
ループのグループデータを加算するグループデータ加算
手段とを具備したことを特徴とするＭＲＩ装置。2. All echo data in two-dimensional k-space are interleaved so that each group includes the same number of components in the phase axis direction, and the first group to the Mth group in the phase axis direction are interleaved. A grouping data collecting means for grouping and collecting data in order from the first group of echo data, and a frequency axis axis for echo data of a group that has been completed in parallel with the collection of echo data of a certain group Group data acquisition means for acquiring group data by performing Fourier transform in the direction and Fourier transform in the phase axis direction, and group data addition means for adding the group data of the first group to the Mth group. MRI device.
前記グループデータ加算手段は、グループデータを取得
する毎にそれまでに取得したグループデータを累積加算
することを特徴とするＭＲＩ装置。3. The MRI apparatus according to claim 2,
The MRI apparatus characterized in that the group data adding means cumulatively adds the group data acquired so far every time the group data is acquired.
ライス軸方向の各成分が各グループに同じ数だけ含まれ
るようにインタリーブして，スライス軸方向に第１グル
ープから第Ｑグループまでにグループ分けし第１グルー
プのエコーデータから順にデータを収集するグループ分
けエコーデータ収集手段と、あるグループのエコーデー
タの収集と並行してそれより前にエコーデータの収集を
完了したグループのデータについてスライス軸方向のフ
ーリエ変換を行ってグループデータを取得するグループ
データ取得手段と、第１グループから第Ｑグループのグ
ループデータを加算するグループデータ加算手段とを具
備したことを特徴とするＭＲＩ装置。4. All echo data of three-dimensional k-space are interleaved so that each group in the slice axis direction is included in the same number in each group, and the first group to the Qth group in the slice axis direction are interleaved. Grouped echo data collection means for grouping and collecting data sequentially from the echo data of the first group, and slices for data of a group for which echo data collection was completed in parallel with collection of echo data of a certain group An MRI apparatus comprising: group data acquisition means for performing group Fourier transformation in the axial direction to obtain group data; and group data addition means for adding group data of the first to Qth groups.
のＭＲＩ装置において、前記グループデータ取得手段ま
たは前記グループデータ加算手段は、各グループデータ
を加算するための位相補償演算を行うことを特徴とする
ＭＲＩ装置。5. The MRI apparatus according to claim 2, wherein the group data acquisition unit or the group data addition unit performs a phase compensation calculation for adding each group data. Characteristic MRI device.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP25833794A JP3472635B2 (en) | 1994-10-24 | 1994-10-24 | MRI equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP25833794A JP3472635B2 (en) | 1994-10-24 | 1994-10-24 | MRI equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08117202A true JPH08117202A (en) | 1996-05-14 |
| JP3472635B2 JP3472635B2 (en) | 2003-12-02 |
Family
ID=17318848
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP25833794A Expired - Fee Related JP3472635B2 (en) | 1994-10-24 | 1994-10-24 | MRI equipment |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3472635B2 (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006021023A (en) * | 2004-06-11 | 2006-01-26 | Hitachi Medical Corp | Magnetic resonance diagnostic system |
| JP2015029677A (en) * | 2013-08-02 | 2015-02-16 | 株式会社日立メディコ | Magnetic resonance imaging apparatus and method |
-
1994
- 1994-10-24 JP JP25833794A patent/JP3472635B2/en not_active Expired - Fee Related
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006021023A (en) * | 2004-06-11 | 2006-01-26 | Hitachi Medical Corp | Magnetic resonance diagnostic system |
| JP2015029677A (en) * | 2013-08-02 | 2015-02-16 | 株式会社日立メディコ | Magnetic resonance imaging apparatus and method |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3472635B2 (en) | 2003-12-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| NL1033584C2 (en) | Method and device for multi-coil MR imaging with hybrid room calibration. | |
| JP4229487B2 (en) | How to correct Maxwell term error | |
| US4780675A (en) | Conjugate symmetry magnetic resonance imaging | |
| JP2014083445A (en) | Magnetic resonance imaging system, and magnetic resonance imaging method | |
| JP3850495B2 (en) | Method and apparatus for generating images from NMR data | |
| US20120301000A1 (en) | Coronary magnetic resonance angiography with signal separation for water and fat | |
| US6486671B1 (en) | MRI image quality improvement using matrix regularization | |
| JP2003052667A (en) | Magnetic resonance imaging system | |
| EP2080032A1 (en) | Mri rf encoding using multiple transmit coils | |
| US6611143B2 (en) | MRI apparatus generating odd and/or even echo images with sensitivity distribution of coils | |
| JP6762284B2 (en) | Magnetic resonance imaging device and noise removal method | |
| US7504828B2 (en) | Frequency synthesizer for RF pulses, MRI apparatus and RF pulse generating method | |
| JP3472635B2 (en) | MRI equipment | |
| JP3884227B2 (en) | Magnetic resonance imaging device | |
| JP3384897B2 (en) | MRI equipment | |
| JP2002165776A (en) | Measurement method in magnetic resonance imaging apparatus and magnetic resonance imaging apparatus | |
| US6181135B1 (en) | MRI system with fractional decimation of acquired data using linear phase shift | |
| EP1685419A1 (en) | Parallel mr imaging method | |
| US6600317B1 (en) | Magnetic resonance imaging method and apparatus | |
| US5072182A (en) | Fourier shifting of voxel boundaries in MR chemical shift imaging spectroscopy | |
| US5089777A (en) | Magnetic resonance imaging apparatus | |
| US6147491A (en) | MRI system with fractional decimation of acquired data using a two-dimensional fourier transformation | |
| JP3447099B2 (en) | MRI equipment | |
| JP4822379B2 (en) | Magnetic resonance imaging system | |
| US20220015661A1 (en) | Magnetic resonance imaging apparatus and method of controlling the same |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080912 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080912 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090912 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090912 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090912 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100912 Year of fee payment: 7 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100912 Year of fee payment: 7 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100912 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110912 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110912 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120912 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120912 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120912 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130912 Year of fee payment: 10 |
|
| LAPS | Cancellation because of no payment of annual fees |