JPH06327650A - Image diagnosis apparatus - Google Patents

Image diagnosis apparatus

Info

Publication number
JPH06327650A
JPH06327650A JP5115999A JP11599993A JPH06327650A JP H06327650 A JPH06327650 A JP H06327650A JP 5115999 A JP5115999 A JP 5115999A JP 11599993 A JP11599993 A JP 11599993A JP H06327650 A JPH06327650 A JP H06327650A
Authority
JP
Japan
Prior art keywords
image
data
images
subtraction
sets
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
Application number
JP5115999A
Other languages
Japanese (ja)
Other versions
JP3408580B2 (en
Inventor
Tetsuji Tsukamoto
鉄二 塚元
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GE Healthcare Japan Corp
Original Assignee
Yokogawa Medical Systems Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yokogawa Medical Systems Ltd filed Critical Yokogawa Medical Systems Ltd
Priority to JP11599993A priority Critical patent/JP3408580B2/en
Publication of JPH06327650A publication Critical patent/JPH06327650A/en
Application granted granted Critical
Publication of JP3408580B2 publication Critical patent/JP3408580B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

PURPOSE:To provide an image diagnosis apparatus with which it is possible to display the changing position of images picked up with time clearly simultaneously with data taking. CONSTITUTION:The subtracted images of two images are calculated simultaneously with data taking and displayed. The phases of two sets of data that are the basis of two images used for calculating the subtracted images, are made the same, and the images are recomposed from the revised data and the subtracted images are calculated. Thus, subtracted images that shows the change of position of image-picking up with time clearly can be displayed simultaneously with data taking. Even if there is a position gap, the subtracted images are obtained showing correct change with time without the influence of the gap.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】この発明は、画像診断装置に関
し、さらに詳しくは、撮像部位の時間的変化を明確に表
示可能な画像診断装置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an image diagnostic apparatus, and more particularly to an image diagnostic apparatus capable of clearly displaying a temporal change of an imaged region.

【0002】[0002]

【従来の技術】図10に、従来の画像診断装置における
撮像と画像表示のタイミングの一例を示す。データ収集
は同一の撮像対象部位に対して連続して行い、複数組の
データR(1),R(2),R(3),R(4)を得る。再構成
は前記各データR(1),R(2),R(3),R(4)が得ら
れる毎に行い、画像I(1),I(2),I(3),I(4)を
得る。表示は前記各画像I(1),I(2),I(3),I
(4)が得られる毎に行う。
2. Description of the Related Art FIG. 10 shows an example of timing of image pickup and image display in a conventional image diagnostic apparatus. Data collection is continuously performed on the same imaging target site, and a plurality of sets of data R (1), R (2), R (3), and R (4) are obtained. Reconstruction is performed every time the data R (1), R (2), R (3), and R (4) are obtained, and the images I (1), I (2), I (3), and I ( 4) is obtained. The images are displayed on the images I (1), I (2), I (3), I
Do this every time (4) is obtained.

【0003】[0003]

【発明が解決しようとする課題】近年、造影ダイナミッ
クスキャンや,ファンクショナルイメージングなどによ
り、同一部位の時間的変化を観察する要求が高まってき
た。この時間的変化は、上記のような画像I(1),I
(2),…をそのまま見るよりも、図11に示すように減
算画像I(2)−I(1),…を見る方が判りやすい。しか
し、従来の画像診断装置では、画像I(1),I(2),…
をデータ収集と並行して表示するだけで、減算画像I
(2)−I(1),…はデータ収集と並行して表示できず、
一連の撮像の終了後に減算画像I(2)−I(1),…を算
出して表示させなければならない問題点があった。
In recent years, there has been an increasing demand for observing temporal changes in the same site by contrast dynamic scanning, functional imaging and the like. This temporal change is caused by the images I (1), I as described above.
It is easier to see the subtracted images I (2) -I (1), ... As shown in FIG. However, in the conventional image diagnostic apparatus, the images I (1), I (2), ...
Is displayed in parallel with the data acquisition, the subtraction image I
(2) -I (1), ... cannot be displayed in parallel with data collection,
There is a problem that the subtracted images I (2) -I (1), ... Must be calculated and displayed after completion of a series of imaging.

【0004】また、従来の画像診断装置では、画像I
(1),I(2),…の間に撮像部位の位置ずれを生じない
ことを前提にしていたが、患者の動きなどによって前記
位置ずれを生じると、減算画像に正しい時間的変化分が
表われず、図12に示すように、撮像部位のエッジ部分
などに不正な変化分が表れてしまう問題点があった。
In the conventional image diagnostic apparatus, the image I
(1), I (2), ... It was premised that the displacement of the imaged region does not occur. However, if the displacement occurs due to patient movement, etc., the correct temporal change amount will appear in the subtracted image. However, there is a problem in that an incorrect change appears in the edge portion of the imaged region as shown in FIG.

【0005】そこで、この発明の第1の目的は、減算画
像I(2)−I(1),…をデータ収集と並行して表示でき
る画像診断装置を提供することにある。また、この発明
の第2の目的は、画像I(1),I(2),…の間に撮像部
位の位置ずれを生じても、減算画像に正しい時間的変化
分が表れるようにした画像診断装置を提供することにあ
る。
Therefore, a first object of the present invention is to provide an image diagnostic apparatus capable of displaying subtracted images I (2) -I (1), ... In parallel with data acquisition. A second object of the present invention is an image in which a correct temporal change appears in the subtracted image even if the position of the imaging region is displaced between the images I (1), I (2), .... To provide a diagnostic device.

【0006】[0006]

【課題を解決するための手段】上記第1の目的を達成す
るため、この発明は、同一の撮像対象部位に対するデー
タの収集を連続して行って複数組のデータを得て、その
複数組のデータから同一の撮像対象部位の複数枚の時系
列画像を得る画像診断装置において、データ収集と並行
してその時までに得られた画像中の2枚の画像の減算画
像を算出する減算画像算出手段と、前記減算画像を表示
する減算画像表示手段とを具備したことを特徴とするの
画像診断装置を提供する。上記第2の目的を達成するた
め、この発明は、同一の撮像対象部位に対するデータの
収集を連続して行って複数組のデータを得て、その複数
組のデータから同一の撮像対象部位の複数枚の時系列画
像を得る画像診断装置において、減算画像を算出する2
枚の画像の元となる2組のデータ間の位相を揃える補正
を行う位相補正手段と、前記補正したデータから再構成
した画像の減算画像を算出する減算画像算出手段とを具
備したことを特徴とする画像診断装置を提供する。
In order to achieve the first object, the present invention continuously collects data for the same imaging target region to obtain a plurality of sets of data, and the plurality of sets of data are collected. In an image diagnostic apparatus that obtains a plurality of time-series images of the same imaging target region from data, subtraction image calculation means for calculating a subtraction image of two images among the images acquired up to that time in parallel with data collection And a subtraction image display means for displaying the subtraction image. In order to achieve the above-mentioned second object, the present invention continuously collects data for the same imaging target site to obtain a plurality of sets of data, and from the plurality of sets of data, a plurality of data of the same imaging target site is acquired. Calculating a subtraction image in an image diagnostic apparatus that obtains one time-series image 2
A phase correction unit that corrects the phase between two sets of data that are the basis of one image, and a subtraction image calculation unit that calculates a subtraction image of an image reconstructed from the corrected data are provided. The present invention provides a diagnostic imaging device.

【0007】なお、この発明の画像診断装置の具体例と
しては、MRI装置やCT装置が挙げられる。
Specific examples of the image diagnostic apparatus of the present invention include an MRI apparatus and a CT apparatus.

【0008】[0008]

【作用】上記第1の目的を達成するための画像診断装置
では、減算画像算出手段を備え、この減算画像算出手段
によりデータ収集と並行して2枚の画像の減算画像を算
出する。そして、その減算画像を減算画像表示手段によ
り表示する。従って、データ収集と並行して、撮像部位
の時間的変化を明確に表示できるようになる。
In the image diagnostic apparatus for achieving the first object, the subtraction image calculation means is provided, and the subtraction image calculation means calculates the subtraction image of the two images in parallel with the data collection. Then, the subtracted image is displayed by the subtracted image display means. Therefore, it becomes possible to clearly display the temporal change of the imaging region in parallel with the data collection.

【0009】上記第2の目的を達成するための画像診断
装置では、位相補正手段を備え、その位相補正手段によ
り減算画像を算出する2枚の画像の元となる2組のデー
タ間の位相を揃える。位置ずれは1次の位相の変化とし
てデータに表われるから、データ間の位相を揃えれば位
置ずれが補正できたことになる。そこで、この補正を行
った後のデータから画像を再構成し、減算画像算出手段
で減算画像を算出すると、位置ずれの影響を受けない減
算画像が得られる。すなわち、正しい時間的変化分が表
れた減算画像が得られる。
An image diagnostic apparatus for achieving the above-mentioned second object is provided with a phase correction means, and the phase between the two sets of data which is the source of the two images for which subtraction images are calculated by the phase correction means. Align. Since the positional deviation appears in the data as a change in the primary phase, the positional deviation can be corrected by aligning the phases between the data. Therefore, when an image is reconstructed from the data after this correction and the subtraction image is calculated by the subtraction image calculation means, a subtraction image that is not affected by the positional deviation can be obtained. That is, a subtracted image showing the correct amount of temporal change can be obtained.

【0010】[0010]

【実施例】以下、図に示す実施例に基づいてこの発明を
さらに詳細に説明する。なお、これによりこの発明が限
定されるものではない。図1は、この発明の画像診断装
置の一実施例のMRI装置1のブロック図である。計算
機2は、操作卓13からの指示に基づき、全体の作動を
制御する。シーケンスコントローラ3は、記憶している
シーケンスに基づいて、勾配磁場駆動回路4を作動さ
せ、マグネットアセンブリ5の勾配磁場コイルで勾配磁
場を発生させる。また、ゲート変調回路7を制御し、R
F発振回路6で発生したRFパルスを所定の波形に変調
して、RF電力増幅器8からマグネットアセンブリ5の
送信コイルに加える。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS The present invention will be described in more detail based on the embodiments shown in the drawings. The present invention is not limited to this. FIG. 1 is a block diagram of an MRI apparatus 1 which is an embodiment of the image diagnostic apparatus of the present invention. The computer 2 controls the overall operation based on the instruction from the console 13. The sequence controller 3 operates the gradient magnetic field driving circuit 4 based on the stored sequence, and causes the gradient magnetic field coil of the magnet assembly 5 to generate a gradient magnetic field. Also, by controlling the gate modulation circuit 7, R
The RF pulse generated by the F oscillation circuit 6 is modulated into a predetermined waveform and applied from the RF power amplifier 8 to the transmission coil of the magnet assembly 5.

【0011】マグネットアセンブリ5の受信コイルで得
られたNMR信号は、前置増幅器9を介して位相検波器
10に入力され、さらにAD変換器11を介して計算機
2に入力される。計算機2は、AD変換器11から得た
NMR信号のデータに基づき、画像を再構成し、表示装
置12で表示する。なお、この計算機2が、減算画像算
出手段,減算画像表示手段および位相補正手段として機
能する。
The NMR signal obtained by the receiving coil of the magnet assembly 5 is input to the phase detector 10 via the preamplifier 9 and further to the computer 2 via the AD converter 11. The computer 2 reconstructs an image based on the NMR signal data obtained from the AD converter 11, and displays the image on the display device 12. The computer 2 functions as a subtraction image calculation means, a subtraction image display means, and a phase correction means.

【0012】図2は、上記MRI装置1において、同一
の撮像対象部位に対するデータR(1),R(2),…の
収集を連続して行い、それら複数組のデータR(1),R
(2),…から同一の撮像対象部位の複数枚の時系列画像
I(1),I(2),…を得る動作と並列に実行され、減算
画像I(2)−I(1),…を算出して、時系列画像I
(1),I(2),…と共に表示装置12に表示する動作の
フロー図である。ステップP1では、画像数(例えば
“4”枚)を反復数Nにセットし、且つ、カウンタiに
“1”をセットする。ステップP2では、カウンタiの
カウント値が反復数Nを越えたかチェックする。越えて
いなければステップP3に進み、越えたなら処理を終了
する。ステップP3では、カウンタiのカウント値が
“1”かチェックする。“1”ならステップP4に進
み、“1”でないならステップP6に進む。ステップP
4では、画像I(i)を表示装置12に表示する。ステッ
プP5では、カウンタiのカウント値をインクリメント
し、前記ステップP2に戻る。ステップP6では、画像
I(i)から画像(1)を減算して、減算画像I(i)−I
(1)を得て、表示装置12に表示する。そして、前記ス
テップP4に進む。
FIG. 2 shows that the MRI apparatus 1 continuously collects data R (1), R (2), ... With respect to the same imaging target region, and a plurality of sets of these data R (1), R (R).
(2), ... Is executed in parallel with the operation of obtaining a plurality of time-series images I (1), I (2), ... Of the same imaging target region, and subtracted images I (2) -I (1), Is calculated and the time series image I is calculated.
(1), I (2), ... Together with FIG. In step P1, the number of images (for example, "4") is set to the number of repetitions N, and the counter i is set to "1". In step P2, it is checked whether the count value of the counter i exceeds the number of iterations N. If it does not exceed, the process proceeds to step P3, and if it does, the process ends. At step P3, it is checked whether the count value of the counter i is "1". If it is "1", the process proceeds to step P4, and if it is not "1", the process proceeds to step P6. Step P
At 4, the image I (i) is displayed on the display device 12. In step P5, the count value of the counter i is incremented, and the process returns to step P2. In Step P6, the image (1) is subtracted from the image I (i), and the subtracted image I (i) -I
(1) is obtained and displayed on the display device 12. Then, the process proceeds to step P4.

【0013】図3に、図2の動作における撮像と画像表
示のタイミングを示す。データ収集は同一の撮像対象部
位に対して連続して行い、複数組のデータR(1),R
(2),R(3),R(4)を得る。再構成は前記各データR
(1),R(2),R(3),R(4)が得られる毎に行い、画
像I(1),I(2),I(3),I(4)を得る。そして、画
像I(1)が得られるとそれを表示し、画像I(2),I
(3),I(4)が得られる毎に減算画像I(2)−I(1),
I(3)−I(1),I(4)−I(1)を算出して、それと画
像I(2),I(3),I(4)とを表示する。以上により、
データ収集と並行して、撮像部位の時間的変化を明確に
表示できることとなる。
FIG. 3 shows the timing of image pickup and image display in the operation of FIG. Data collection is performed continuously for the same imaging target region, and multiple sets of data R (1), R
(2), R (3) and R (4) are obtained. Reconstruction is for each data R
(1), R (2), R (3), R (4) is performed every time it is obtained, and images I (1), I (2), I (3), and I (4) are obtained. When the image I (1) is obtained, it is displayed, and the images I (2), I
(3), I (4) each time the subtracted image I (2) -I (1),
I (3) -I (1) and I (4) -I (1) are calculated, and the images I (2), I (3) and I (4) are displayed. From the above,
In parallel with the data collection, it is possible to clearly display the temporal change of the imaging region.

【0014】図4の(a),(b)は、データR(1),
R(2)の絶対値を表している。データR(1),R(2)の
間で撮像部位の位置ずれがあっても、絶対値には影響し
ない。図5の(a),(b)は、データR(1),R(2)
の直交軸x,y上の位相を表している。データR(1),
R(2)の間で撮像部位の位置ずれがあると、1次の位相
の変化としてデータに影響する。そこで、データ間の位
相を揃えることにより位置ずれを補正できることにな
る。そこで、次に、この位置ずれの補正を行う動作を説
明する。
FIGS. 4A and 4B show data R (1),
It represents the absolute value of R (2). Even if there is a displacement of the imaged region between the data R (1) and R (2), the absolute value is not affected. 5A and 5B show data R (1) and R (2).
Represents the phase on the orthogonal axes x and y. Data R (1),
If there is a positional deviation of the imaged region between R (2), it affects the data as a change in the primary phase. Therefore, it is possible to correct the positional deviation by aligning the phases of the data. Therefore, next, an operation of correcting the positional deviation will be described.

【0015】図6は、上記MRI装置1において、同一
の撮像対象部位に対するデータR(1),R(2),…の
収集を連続して行い、それら複数組のデータR(1),R
(2),…から同一の撮像対象部位の複数枚の時系列画像
I(1),I(2),…を得る動作と並列に実行され、画像
I(2),…について画像I(1)に合せる位置ずれの補正
を行った上で、それら補正後の画像I’(2),…より
減算画像I’(2)−I(1),…を算出して、時系列画像
I(1),I(2),…と共に表示装置12に表示する動作
のフロー図である。ステップS1では、画像数(例えば
“4”枚)を反復数Nにセットし、且つ、カウンタiに
“1”をセットする。ステップS2では、カウンタiの
カウント値が反復数Nを越えたかチェックする。越えて
いなければステップS3に進み、越えたなら処理を終了
する。ステップS3では、画像I(i)を2次元逆フーリ
エ変換し、画像I(i)の元のデータR(i)を求める。ス
テップS4では、データR(i)の位相φ{R(i)}を求
める。ステップS5では、カウンタiのカウント値が
“1”かチェックする。“1”ならステップS6に進
み、“1”でないならステップS8に進む。ステップS
6では、画像I(i)を表示装置12に表示する。ステッ
プS7では、カウンタiのカウント値をインクリメント
し、前記ステップS2に戻る。
FIG. 6 shows that the MRI apparatus 1 continuously collects data R (1), R (2), ... For the same imaging target region, and sets a plurality of sets of data R (1), R (R).
(2), ... Is performed in parallel with the operation of obtaining a plurality of time-series images I (1), I (2), ... Of the same imaging target region, and image I (1), I (2) ,. ), The subtraction images I ′ (2) −I (1), ... Are calculated from the corrected images I ′ (2) ,. 1), I (2), ... Together with FIG. In step S1, the number of images (for example, "4") is set to the number of repetitions N, and the counter i is set to "1". In step S2, it is checked whether the count value of the counter i exceeds the number of iterations N. If not exceeded, the process proceeds to step S3, and if exceeded, the process ends. In step S3, the image I (i) is two-dimensional inverse Fourier transformed to obtain the original data R (i) of the image I (i). In step S4, the phase φ {R (i)} of the data R (i) is obtained. In step S5, it is checked whether the count value of the counter i is "1". If it is "1", the process proceeds to step S6, and if it is not "1", the process proceeds to step S8. Step S
At 6, the image I (i) is displayed on the display device 12. In step S7, the count value of the counter i is incremented, and the process returns to step S2.

【0016】ステップS8では、データR(i)の位相φ
{R(i)}からデータR(1)の位相φ{R(1)}を減算
し、位相差Δφ{R(i)}を求める。ステップS9で
は、位相差Δφ{R(i)}の直交軸x,yの近傍の平均
値をとり、その平均値に対するカーブフィッティングに
より平均位相差Δφ’{R(i)}を求める。ステップS
10では、データR(i)の位相φ{R(i)}から平均位
相差Δφ’{R(i)}を減算し、それをデータR(i)の
位相φ{R(i)}とする。ステップS11では、データ
R(i)を2次元フーリエ変換し、画像I’(i)を求め
る。ステップS12では、画像I’(i)から画像(1)を
減算し、減算画像I’(i)-I(1)を得て、表示装置12
に表示し、前記ステップS6に進む。
In step S8, the phase φ of the data R (i)
The phase difference Δφ {R (i)} is obtained by subtracting the phase φ {R (1)} of the data R (1) from {R (i)}. In step S9, the average value of the phase difference Δφ {R (i)} in the vicinity of the orthogonal axes x and y is calculated, and the average phase difference Δφ ′ {R (i)} is obtained by curve fitting to the average value. Step S
In 10, the average phase difference Δφ ′ {R (i)} is subtracted from the phase φ {R (i)} of the data R (i), and this is taken as the phase φ {R (i)} of the data R (i). To do. In step S11, the data R (i) is two-dimensionally Fourier transformed to obtain an image I ′ (i). In step S12, the image (1) is subtracted from the image I '(i) to obtain the subtracted image I' (i) -I (1), and the display device 12
And the process proceeds to step S6.

【0017】図7に、図6の動作における撮像と画像表
示のタイミングを示す。データ収集は同一の撮像対象部
位に対して連続して行い、複数組のデータR(1),R
(2),R(3),R(4)を得る。再構成は前記各データR
(1),R(2),R(3),R(4)が得られる毎に行い、画
像I(1),I(2),I(3),I(4)を得る。補正画像
I’(2),I’(3),I’(4)は、画像I(2),I
(3),I(4)が得られる毎に得る。そして、画像I(1)
が得られるとそれを表示し、補正画像I’(2),I’
(3),I’(4)が得られる毎に減算画像I’(2)−I
(1),I’(3)−I(1),I’(4)−I(1)を算出し
て、それと画像I(2),I(3),I(4)を表示する。以
上により、データ収集と並行して、撮像部位の時間的変
化を明確に且つ位置ずれの影響を受けないで表示できる
こととなる。
FIG. 7 shows the timing of image pickup and image display in the operation of FIG. Data collection is performed continuously for the same imaging target region, and multiple sets of data R (1), R
(2), R (3) and R (4) are obtained. Reconstruction is for each data R
(1), R (2), R (3), R (4) is performed every time it is obtained, and images I (1), I (2), I (3), and I (4) are obtained. The corrected images I ′ (2), I ′ (3), I ′ (4) are the images I (2), I
Obtain each time (3) and I (4) are obtained. And image I (1)
When it is obtained, it is displayed and corrected images I ′ (2), I ′
(3), I '(4) each time subtraction image I' (2) -I
(1), I '(3) -I (1), I' (4) -I (1) are calculated and the images I (2), I (3) and I (4) are displayed. As described above, in parallel with the data collection, it is possible to clearly display the temporal change of the imaged region without being affected by the positional deviation.

【0018】なお、ノイズが少ない場合には、上記ステ
ップS8,S9,S10を省略し、その代りに、データ
R(1)の位相φ{R(1)}でデータR(i)の位相φ{R
(i)}を置換するステップを設けてもよい。
When the noise is small, the steps S8, S9 and S10 are omitted, and instead, the phase φ of the data R (1) is replaced with the phase φ of the data R (i). {R
A step of substituting (i)} may be provided.

【0019】上記実施例の位置ずれ補正方法は、画像I
(i)を逆フーリエ変換してデータR(i)を算出し、位相
φ{R(i)}を求めて位相補正を行うため、実際に収集
したローデータR(i)がなくても実施できる。このた
め、データ収集と並行して行う場合に限定されず、一連
のデータの収集後に減算画像を得たい場合でも実施でき
る利点がある。一方、データ収集と並行して位置ずれ補
正を行うことに限定するなら、実際に収集したローデー
タR(i)を利用して位置ずれ補正を行うことが出来る。
そこで、次に、実際に収集したローデータR(i)を利用
して位置ずれ補正を行う動作を説明する。
The positional deviation correction method of the above embodiment is based on the image I
(i) is inverse-Fourier-transformed to calculate data R (i), and phase φ {R (i)} is calculated to correct the phase. it can. Therefore, the present invention is not limited to the case where the data acquisition is performed in parallel, and there is an advantage that the subtraction image can be obtained even after the series of data is acquired. On the other hand, if it is limited to performing the positional deviation correction in parallel with the data collection, the positional deviation correction can be performed using the raw data R (i) actually collected.
Therefore, next, an operation of correcting the positional deviation using the raw data R (i) actually collected will be described.

【0020】図8は、上記MRI装置1において、同一
の撮像対象部位に対するデータR(1),R(2),…の
収集を連続して行い、データR(1)より画像I(1)を求
めて表示し、データR(2),…については位置ずれの補
正を行った上で、それらデータR(2),…より画像I
(2),…を算出し、且つ、減算画像I(2)−I(1),…
を算出して、両画像を表示装置12に表示する動作のフ
ロー図である。ステップR1では、画像数(例えば
“4”枚)を反復数Nにセットし、且つ、カウンタiに
“1”をセットする。ステップR2では、カウンタiの
カウント値が反復数Nを越えたかチェックする。越えて
いなければステップR3に進み、越えたなら処理を終了
する。ステップR3では、データR(i)を収集して、そ
の位相φ{R(i)}を求める。ステップR4では、カウ
ンタiのカウント値が“1”かチェックする。“1”な
らステップR5に進み、“1”でないならステップR8
に進む。ステップR5では、データR(1)を2次元フー
リエ変換し、画像I(1)を求める。ステップR6では、
画像I(i)を表示装置12に表示する。ステップR7で
は、カウンタiのカウント値をインクリメントし、前記
ステップR2に戻る。
FIG. 8 shows that the MRI apparatus 1 continuously collects data R (1), R (2), ... For the same imaging target region, and obtains an image I (1) from the data R (1). Is calculated and displayed, the positional deviation of the data R (2), ... Is corrected, and then the image I is obtained from the data R (2) ,.
(2), ... And subtracted image I (2) -I (1) ,.
FIG. 7 is a flowchart of an operation for calculating and displaying both images on the display device 12. In step R1, the number of images (for example, "4") is set to the number of repetitions N, and the counter i is set to "1". In step R2, it is checked whether the count value of the counter i exceeds the number of iterations N. If not exceeded, the process proceeds to step R3, and if exceeded, the process is terminated. In step R3, the data R (i) is collected and its phase φ {R (i)} is obtained. At step R4, it is checked whether the count value of the counter i is "1". If it is "1", proceed to step R5. If it is not "1", step R8.
Proceed to. In step R5, the data R (1) is two-dimensionally Fourier transformed to obtain the image I (1). In step R6,
The image I (i) is displayed on the display device 12. In step R7, the count value of the counter i is incremented, and the process returns to step R2.

【0021】ステップR8では、データR(i)の位相φ
{R(i)}からデータR(1)の位相φ{R(1)}を減算
し、位相差Δφ{R(i)}を求める。ステップR9で
は、位相差Δφ{R(i)}の直交軸x,yの近傍の平均
値をとり、その平均値に対するカーブフィッティングに
より平均位相差Δφ’{R(i)}を求める。ステップR
10では、データR(i)の位相φ{R(i)}から平均位
相差Δφ’{R(i)}を減算し、それをデータR(i)の
位相φ{R(i)}とする。ステップR11では、データ
R(i)を2次元フーリエ変換し、画像I(i)を求める。
ステップR12では、画像I(i)から画像(1)を減算
し、減算画像I(i)−I(1)を得て、表示装置12に表
示し、前記ステップR6に進む。
In step R8, the phase φ of the data R (i)
The phase difference Δφ {R (i)} is obtained by subtracting the phase φ {R (1)} of the data R (1) from {R (i)}. In step R9, the average value of the phase difference Δφ {R (i)} in the vicinity of the orthogonal axes x and y is calculated, and the average phase difference Δφ ′ {R (i)} is obtained by curve fitting to the average value. Step R
In 10, the average phase difference Δφ ′ {R (i)} is subtracted from the phase φ {R (i)} of the data R (i), and this is taken as the phase φ {R (i)} of the data R (i). To do. In step R11, the data R (i) is two-dimensionally Fourier transformed to obtain the image I (i).
In step R12, the image (1) is subtracted from the image I (i) to obtain the subtracted image I (i) -I (1), which is displayed on the display device 12, and the process proceeds to step R6.

【0022】図9に、図8の動作における撮像と画像表
示のタイミングを示す。データ収集は同一の撮像対象部
位に対して連続して行い、複数組のデータR(1),R
(2),R(3),R(4)を得る。位相算出は、各データR
(1),R(2),R(3),R(4)が得られる毎に行い、位
相φ{R(1)},φ{R(2)},φ{R(3)},φ{R
(4)}を得る。再構成は前記各位相φ{R(1)},φ
{R(2)},φ{R(3)},φ{R(4)}が得られる毎
に行い、画像I(1),I(2),I(3),I(4)を得る。
そして、画像I(1)が得られるとそれを表示し、画像I
(2),I(3),I(4)が得られる毎に減算画像I(2)−
I(1),I(3)−I(1),I(4)−I(1)を算出し、そ
れと画像I(2),I(3),I(4)を表示する。以上によ
り、データ収集と並行して、撮像部位の時間的変化を明
確に且つ位置ずれの影響を受けないで表示できることと
なる。
FIG. 9 shows the timing of image pickup and image display in the operation of FIG. Data collection is performed continuously for the same imaging target region, and multiple sets of data R (1), R
(2), R (3) and R (4) are obtained. Phase calculation is for each data R
(1), R (2), R (3), R (4) are performed every time they are obtained, and the phases φ {R (1)}, φ {R (2)}, φ {R (3)}, φ {R
(4)} is obtained. Reconstruction is performed for each phase φ {R (1)}, φ
Each time {R (2)}, φ {R (3)}, φ {R (4)} is obtained, the image I (1), I (2), I (3), I (4) is obtained. obtain.
When the image I (1) is obtained, it is displayed, and the image I (1) is displayed.
Every time (2), I (3), and I (4) are obtained, the subtraction image I (2)-
I (1), I (3) -I (1), I (4) -I (1) are calculated, and the images I (2), I (3) and I (4) are displayed. As described above, in parallel with the data collection, it is possible to clearly display the temporal change of the imaged region without being affected by the positional deviation.

【0023】[0023]

【発明の効果】この発明のMRI装置によれば、データ
収集と並行して、撮像部位の時間的変化が明確に判る減
算画像を表示できる。また、この発明のMRI装置によ
れば、位置ずれがあってもその影響を受けず、正しい時
間的変化分が表れた減算画像を得られる。
According to the MRI apparatus of the present invention, it is possible to display a subtraction image in which the temporal change of the imaged region can be clearly understood in parallel with the data acquisition. Further, according to the MRI apparatus of the present invention, even if there is a positional shift, it is not affected and a subtracted image showing a correct temporal change can be obtained.

【図面の簡単な説明】[Brief description of drawings]

【図1】この発明のMRI装置の一実施例のブロック図
である。
FIG. 1 is a block diagram of an embodiment of an MRI apparatus of the present invention.

【図2】データ収集と並行して減算画像を表示する処理
動作のフロー図である。
FIG. 2 is a flowchart of a processing operation of displaying a subtraction image in parallel with data collection.

【図3】図3の動作におけるデータ収集と表示のタイミ
ング図である。
FIG. 3 is a timing diagram of data collection and display in the operation of FIG.

【図4】データの絶対値を表わす概念図である。FIG. 4 is a conceptual diagram showing the absolute value of data.

【図5】データの位相を表わす概念図である。FIG. 5 is a conceptual diagram showing a phase of data.

【図6】データ収集と並行して位置ずれ補正し且つ減算
画像を表示する処理動作のフロー図である。
FIG. 6 is a flowchart of a processing operation of performing positional deviation correction and displaying a subtraction image in parallel with data collection.

【図7】図6の動作におけるデータ収集と表示のタイミ
ング図である。
7 is a timing diagram of data collection and display in the operation of FIG.

【図8】データ収集と並行してローデータを用いて位置
ずれ補正し且つ減算画像を表示する処理動作のフロー図
である。
FIG. 8 is a flow chart of a processing operation of performing positional deviation correction using raw data and displaying a subtraction image in parallel with data collection.

【図9】図8の動作におけるデータ収集と表示のタイミ
ング図である。
9 is a timing diagram of data collection and display in the operation of FIG.

【図10】従来のデータ収集と表示のタイミング図であ
る。
FIG. 10 is a timing diagram of conventional data collection and display.

【図11】減算画像の説明図である。FIG. 11 is an explanatory diagram of a subtraction image.

【図12】位置ずれがある場合の減算画像の説明図であ
る。
FIG. 12 is an explanatory diagram of a subtraction image when there is a position shift.

【符号の説明】[Explanation of symbols]

1 MRI装置 2 計算機 12 表示装置 1 MRI device 2 calculator 12 display device

───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.5 識別記号 庁内整理番号 FI 技術表示箇所 9219−2J G01N 24/02 J ─────────────────────────────────────────────────── ─── Continuation of the front page (51) Int.Cl. 5 Identification code Internal reference number FI technical display location 9219-2J G01N 24/02 J

Claims (3)

【特許請求の範囲】[Claims] 【請求項1】 同一の撮像対象部位に対するデータの収
集を連続して行って複数組のデータを得て、その複数組
のデータから同一の撮像対象部位の複数枚の時系列画像
を得る画像診断装置において、 データ収集と並行してその時までに得られた画像中の2
枚の画像の減算画像を算出する減算画像算出手段と、前
記減算画像を表示する減算画像表示手段とを具備したこ
とを特徴とする画像診断装置。
1. An image diagnostic in which data is continuously collected for the same imaging target region to obtain a plurality of sets of data, and a plurality of time-series images of the same imaging target region are obtained from the plurality of sets of data. 2 of the images acquired by then in parallel with the data acquisition
An image diagnostic apparatus comprising: a subtraction image calculation means for calculating a subtraction image of a single image; and a subtraction image display means for displaying the subtraction image.
【請求項2】 請求項1に記載の画像診断装置におい
て、減算画像を算出する2枚の画像の元となる2組のデ
ータ間の位相を揃える補正を行う位相補正手段をさらに
具備すると共に、前記減算画像算出手段は、前記補正し
たデータから再構成した画像の減算画像を算出すること
を特徴とする画像診断装置。
2. The image diagnostic apparatus according to claim 1, further comprising a phase correction unit that corrects a phase between two sets of data that are the bases of the two images for calculating the subtracted image, and The subtraction image calculation means calculates a subtraction image of an image reconstructed from the corrected data, the image diagnostic apparatus.
【請求項3】 同一の撮像対象部位に対するデータの収
集を連続して行って複数組のデータを得て、その複数組
のデータから同一の撮像対象部位の複数枚の時系列画像
を得る画像診断装置において、 減算画像を算出する2枚の画像の元となる2組のデータ
間の位相を揃える補正を行う位相補正手段と、前記補正
したデータから再構成した画像の減算画像を算出する減
算画像算出手段とを具備したことを特徴とする画像診断
装置。
3. An image diagnostic for obtaining data of a plurality of sets by continuously collecting data for the same image-capturing target site and obtaining a plurality of time-series images of the same image-capturing target site from the plurality of sets of data. In the apparatus, a phase correction unit that corrects a phase between two sets of data that are the bases of two images for which a subtraction image is calculated, and a subtraction image that calculates a subtraction image of an image reconstructed from the corrected data An image diagnostic apparatus comprising: a calculating unit.
JP11599993A 1993-05-18 1993-05-18 Diagnostic imaging device Expired - Fee Related JP3408580B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP11599993A JP3408580B2 (en) 1993-05-18 1993-05-18 Diagnostic imaging device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP11599993A JP3408580B2 (en) 1993-05-18 1993-05-18 Diagnostic imaging device

Publications (2)

Publication Number Publication Date
JPH06327650A true JPH06327650A (en) 1994-11-29
JP3408580B2 JP3408580B2 (en) 2003-05-19

Family

ID=14676345

Family Applications (1)

Application Number Title Priority Date Filing Date
JP11599993A Expired - Fee Related JP3408580B2 (en) 1993-05-18 1993-05-18 Diagnostic imaging device

Country Status (1)

Country Link
JP (1) JP3408580B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1014900A (en) * 1996-07-02 1998-01-20 Hitachi Medical Corp Medical image diagnostic device
JP4678916B2 (en) * 2000-06-06 2011-04-27 株式会社日立メディコ Magnetic resonance imaging system

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1014900A (en) * 1996-07-02 1998-01-20 Hitachi Medical Corp Medical image diagnostic device
JP4678916B2 (en) * 2000-06-06 2011-04-27 株式会社日立メディコ Magnetic resonance imaging system

Also Published As

Publication number Publication date
JP3408580B2 (en) 2003-05-19

Similar Documents

Publication Publication Date Title
Szumowski et al. Phase unwrapping in the three-point Dixon method for fat suppression MR imaging.
US7432710B2 (en) Apparatus and method for reducing image artifact
US8368401B2 (en) Techniques for correcting measurement artifacts in magnetic resonance thermometry
JP4121174B2 (en) Method for correcting Maxwell term errors generated by a nuclear magnetic resonance system
JPH11113878A (en) Magnetic resonance imaging method
JPH0731606A (en) Magnetic resonance tomographic camera
JP2007075387A (en) Breath-holding mr imaging method, mri apparatus, breath-holding tomographic imaging method, and tomographic imaging apparatus
JP6986659B2 (en) Magnetic resonance imaging device, Nyquist ghost correction method, and Nyquist ghost correction program
JP4052721B2 (en) System for correcting Maxwell term phase error generated by nuclear magnetic resonance system and system for forming phase contrast angiogram by nuclear magnetic resonance system
US20080218167A1 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
JPH05123312A (en) Magnetic resonance imaging device
JP3408580B2 (en) Diagnostic imaging device
JP2641486B2 (en) NMR imaging equipment
US5785042A (en) Magnetic resonance imaging method providing for correction of striation artifacts
JP3472620B2 (en) MRI equipment
JP5105823B2 (en) Magnetic resonance imaging apparatus, medical image processing method, program, and medical image display system
JP2001187039A5 (en) MRI equipment
JP3547552B2 (en) Inspection device using nuclear magnetic resonance
US7415142B2 (en) Method of visualizing the perfusion of an organ while utilizing a perfusion measurement
JPH06327649A (en) Method for recomposing images in nuclear magnetic resonance imaging apparatus
JP3152684B2 (en) MR device
JPH08182661A (en) Magnetic resonance photographing method
US7365541B2 (en) MRI apparatus which automatically determines and displays operating instructions
JPH07108288B2 (en) NMR imaging method
JPH06296600A (en) Magnetic resonance angiography and apparatus therefor by using two dimensional phase contrast

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

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: 20080314

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090314

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100314

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100314

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100314

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: 20100314

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: 20100314

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110314

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110314

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120314

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120314

Year of fee payment: 9

LAPS Cancellation because of no payment of annual fees