JPH0429381B2 - - Google Patents
Info
- Publication number
- JPH0429381B2 JPH0429381B2 JP59006266A JP626684A JPH0429381B2 JP H0429381 B2 JPH0429381 B2 JP H0429381B2 JP 59006266 A JP59006266 A JP 59006266A JP 626684 A JP626684 A JP 626684A JP H0429381 B2 JPH0429381 B2 JP H0429381B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- images
- contrast
- heartbeat
- ray
- 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.)
- Expired - Lifetime
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/54—Control of apparatus or devices for radiation diagnosis
- A61B6/541—Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N5/00—Details of television systems
- H04N5/30—Transforming light or analogous information into electric information
- H04N5/32—Transforming X-rays
- H04N5/3205—Transforming X-rays using subtraction imaging techniques
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physiology (AREA)
- Optics & Photonics (AREA)
- High Energy & Nuclear Physics (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Description
【発明の詳細な説明】
[発明の技術分野]
本発明は被検体を透過したX線を検出処理する
ことにより、心機能診断に有効な画像を表示する
ことができるX線を用いた心血管機能診断装置に
関するものである。
ことにより、心機能診断に有効な画像を表示する
ことができるX線を用いた心血管機能診断装置に
関するものである。
[発明の技術的背景とその問題点]
従来、X線診断装置として米国特許第4204225
号及び第4204226号に開示されているデイジタル
フルオロスコープ装置(以下「DF」という)が
知られている。
号及び第4204226号に開示されているデイジタル
フルオロスコープ装置(以下「DF」という)が
知られている。
しかしながらこのDFにおいては、動画像を目
視観察することによつてのみ心機能を診断してい
たため、その定量評価が難しかつた。
視観察することによつてのみ心機能を診断してい
たため、その定量評価が難しかつた。
また、このDFにおいてはマスク像として心拍
位相の異なる造影前画像を時間的に平均していた
ため、この造影前画像と連続した造影画像間でサ
ブトラクシヨンを実施して生成したサブトラクシ
ヨン画像には、造影血液ばかりでなく心拍動によ
る偽像も表現されてしまうという問題があつた。
位相の異なる造影前画像を時間的に平均していた
ため、この造影前画像と連続した造影画像間でサ
ブトラクシヨンを実施して生成したサブトラクシ
ヨン画像には、造影血液ばかりでなく心拍動によ
る偽像も表現されてしまうという問題があつた。
核医学診断装置の分野においては、例えば「核
医学」第19巻第5号、第765頁〜第776頁(前田寿
登ほか)に核医学像を使用して心振幅画像と心位
相画像とをフーリエ展開法でソフト的に求める手
法が開示されている。
医学」第19巻第5号、第765頁〜第776頁(前田寿
登ほか)に核医学像を使用して心振幅画像と心位
相画像とをフーリエ展開法でソフト的に求める手
法が開示されている。
しかしながら、DFによる画像と核医学診断装
置による画像とでは画像収集法及び前処理法が全
く異なつている。
置による画像とでは画像収集法及び前処理法が全
く異なつている。
[発明の目的]
本発明は上記事情に鑑みてなされたものであ
り、周期的に拍動している心血管の機能を定量的
に表現する診断画像を作成表示することができる
心血管機能診断装置の提供を目的とするものであ
る。
り、周期的に拍動している心血管の機能を定量的
に表現する診断画像を作成表示することができる
心血管機能診断装置の提供を目的とするものであ
る。
[発明の概要]
上記目的を達成するために、本発明に係るX線
診断装置は、X線を被写体に曝射するX線発生器
と、被写体を透過したX線からX線透過画像を収
集する画像収集装置と、X線透過画像を画像処理
する画像処理装置と、この画像処理装置による画
像処理後の画像を表示する画像表示装置とを有す
るX線診断装置において、前記画像処理装置は少
なくとも1心拍分の血管造影前のX線透過画像及
び血管造影後のX線造影画像から少なくとも1心
拍分の造影血管画像を得る手段と、この手段によ
つて得られた少なくとも1心拍分の造影血管画像
それぞれに所定の重み付けをして2つの画像を得
る手段と、この手段によつて得られた2つの画像
の相乗平均をとることにより振幅画像を作成する
手段と、前記2つの画像の比の逆正接をとること
により位相画像を作成する手段とを備えたことを
特徴とするものである。
診断装置は、X線を被写体に曝射するX線発生器
と、被写体を透過したX線からX線透過画像を収
集する画像収集装置と、X線透過画像を画像処理
する画像処理装置と、この画像処理装置による画
像処理後の画像を表示する画像表示装置とを有す
るX線診断装置において、前記画像処理装置は少
なくとも1心拍分の血管造影前のX線透過画像及
び血管造影後のX線造影画像から少なくとも1心
拍分の造影血管画像を得る手段と、この手段によ
つて得られた少なくとも1心拍分の造影血管画像
それぞれに所定の重み付けをして2つの画像を得
る手段と、この手段によつて得られた2つの画像
の相乗平均をとることにより振幅画像を作成する
手段と、前記2つの画像の比の逆正接をとること
により位相画像を作成する手段とを備えたことを
特徴とするものである。
[発明の実施例]
以下に本発明の実施例を詳細に説明する。
第1図はX線診断装置のブロツク図である。
同図において、1は被検体2の曝射するX線を
発生するためのX線発生器である。
発生するためのX線発生器である。
3は画像収集装置であり、被検体2を透過した
透過X線をアナログ信号に変換する例えばイメー
ジインテシフアイアと、光学系、撮像管及び撮像
管制御装置からなるX線検出装置と、該X線検出
装置から出力されるアナログ信号をデイジタル信
号に変換し出力するA/D変換器とから構成され
いる。
透過X線をアナログ信号に変換する例えばイメー
ジインテシフアイアと、光学系、撮像管及び撮像
管制御装置からなるX線検出装置と、該X線検出
装置から出力されるアナログ信号をデイジタル信
号に変換し出力するA/D変換器とから構成され
いる。
4は画像処理装置であり、第2図に具体的に示
すように、各種の制御信号を送出する画像処理コ
ントローラ20と画像収集装置3より出力される
収集画像と所定の処理を施した処理済画像を記憶
する画像記憶装置(例えばフレームメモリ)25
と、その記憶された両画像間の所定の関数領域に
ついての濃度の時間変化図を作成する濃度時間変
化図作成処理装置21と、記憶された画像の心拍
番号及び心拍位相を検出する心拍動周期検出処理
装置22と、その検出された濃度の時間変化図と
心拍位相に対応して同一心拍位相iを有する造影
前画像MIiと造影画像CIiとの間で造影剤が存在
しない背景(例えば骨、血液以外の軟組織等)を
除去し、造影血液画像SIiを作成する背景処理装
置23と、その造影血液画像SIiの一心拍分の画
像間に時間フイルター処理を実施し心血管の振幅
画像AIと位相画像PIを作成する機能画像作成装
置24と、セレクタ26とから構成されている。
すように、各種の制御信号を送出する画像処理コ
ントローラ20と画像収集装置3より出力される
収集画像と所定の処理を施した処理済画像を記憶
する画像記憶装置(例えばフレームメモリ)25
と、その記憶された両画像間の所定の関数領域に
ついての濃度の時間変化図を作成する濃度時間変
化図作成処理装置21と、記憶された画像の心拍
番号及び心拍位相を検出する心拍動周期検出処理
装置22と、その検出された濃度の時間変化図と
心拍位相に対応して同一心拍位相iを有する造影
前画像MIiと造影画像CIiとの間で造影剤が存在
しない背景(例えば骨、血液以外の軟組織等)を
除去し、造影血液画像SIiを作成する背景処理装
置23と、その造影血液画像SIiの一心拍分の画
像間に時間フイルター処理を実施し心血管の振幅
画像AIと位相画像PIを作成する機能画像作成装
置24と、セレクタ26とから構成されている。
5は画像表示装置で、画像処理装置4から出力
される画像のデイジタル信号をアナログ信号に変
換するD/A変換器とそのアナログ信号を画像と
して表示するデイスプレイとから構成されてい
る。
される画像のデイジタル信号をアナログ信号に変
換するD/A変換器とそのアナログ信号を画像と
して表示するデイスプレイとから構成されてい
る。
6は画像収集装置3、画像処理装置4、画像表
示装置5、X線発生制御器8、造影剤注入装置9
等の装置各部の動作を制御するシステムコントロ
ーラであり、中央演算処理装置を内蔵している。
示装置5、X線発生制御器8、造影剤注入装置9
等の装置各部の動作を制御するシステムコントロ
ーラであり、中央演算処理装置を内蔵している。
7はシステムコンソールであり、システムコン
トローラ6内に予めプログラムされている被検体
2内の目的部位の応じた撮影条件、画像処理デー
タ、画像処理モード及びシーケンス等を選択する
選択キーを有している。
トローラ6内に予めプログラムされている被検体
2内の目的部位の応じた撮影条件、画像処理デー
タ、画像処理モード及びシーケンス等を選択する
選択キーを有している。
X線発生制御器8はシステムコントローラ6に
より制御され、X線発生器1の曝射条件を設定す
るようになつている。
より制御され、X線発生器1の曝射条件を設定す
るようになつている。
造影剤注入装置9はシステムコントローラ6に
より制御され、所定時期に所定量の造影剤を被検
体2内へ注入するようになつている。
より制御され、所定時期に所定量の造影剤を被検
体2内へ注入するようになつている。
次に前記画像処理装置4の具体的構成例を第2
図乃至第5図を参照して詳述する。
図乃至第5図を参照して詳述する。
画像処理装置4の画像処理コントローラ20
は、システムコントローラ6から送出された制御
信号C1を受けて装置内各部の動作シーケンス等
の制御を与える各種の制御信号を出力するように
なつている。
は、システムコントローラ6から送出された制御
信号C1を受けて装置内各部の動作シーケンス等
の制御を与える各種の制御信号を出力するように
なつている。
また、画像処理コントローラ20は心拍動周期
検出処理装置22が検出した心拍番号および位相
番号等の情報を制御線uを介して入力すると共
に、背景処理装置23に対しては制御線P3を介
して、機能画像作成装置24に対しては制御線
P4を介して、画像記憶装置25に対しては制御
線Q1,Q2,……QN-1,QNを介して、セレクタ2
6に対しては制御線Rを介してそれぞれ制御信号
を送出するようになつている。
検出処理装置22が検出した心拍番号および位相
番号等の情報を制御線uを介して入力すると共
に、背景処理装置23に対しては制御線P3を介
して、機能画像作成装置24に対しては制御線
P4を介して、画像記憶装置25に対しては制御
線Q1,Q2,……QN-1,QNを介して、セレクタ2
6に対しては制御線Rを介してそれぞれ制御信号
を送出するようになつている。
濃度時間変化図作成処理装置21は、画像デー
タバスを介して画像収集装置3から伝送される
収集画像データ又は後述する画像記憶装置25か
ら画像データバス,のいずれかを介して伝送
される画像データを入力し、システムコントロー
ラ6で予め設定された所定の関心領域内の濃度値
の時間(画像番号)変化を計測し、これを画像記
憶装置25内へ記憶するようになつている。
タバスを介して画像収集装置3から伝送される
収集画像データ又は後述する画像記憶装置25か
ら画像データバス,のいずれかを介して伝送
される画像データを入力し、システムコントロー
ラ6で予め設定された所定の関心領域内の濃度値
の時間(画像番号)変化を計測し、これを画像記
憶装置25内へ記憶するようになつている。
心拍動周期検出処理装置22は、濃度時間変化
図作成処理装置21で作成された濃度時間変化図
を画像データバスを介して入力し、その濃度の
極大値、極少値を示す時間(画像番号)を求め、
かつ心拍動の心拍番号及び位相番号を制御線uを
介して画像処理コントローラ20へ伝送するよう
になつている。
図作成処理装置21で作成された濃度時間変化図
を画像データバスを介して入力し、その濃度の
極大値、極少値を示す時間(画像番号)を求め、
かつ心拍動の心拍番号及び位相番号を制御線uを
介して画像処理コントローラ20へ伝送するよう
になつている。
画像記憶装置25は、画像処理コントローラ2
0から制御線Q1〜QNを介して送られる制御信号
により制御され、画像データバスを介して入力
される画像データを前記制御信号により定められ
る所定のアドレスへ順次記憶するとともに、その
記憶された画像データを画像データバス,へ
出力するようになつている。
0から制御線Q1〜QNを介して送られる制御信号
により制御され、画像データバスを介して入力
される画像データを前記制御信号により定められ
る所定のアドレスへ順次記憶するとともに、その
記憶された画像データを画像データバス,へ
出力するようになつている。
背景処理装置23は、第3図に示すように掛算
器231,232と減算器233とから構成さ
れ、前記画像処理コントローラ20は制御線P3
を介して掛算器231,232の動作を制御する
とともに、これらに荷重係数Gi,Gjを供給する
ようになつている。
器231,232と減算器233とから構成さ
れ、前記画像処理コントローラ20は制御線P3
を介して掛算器231,232の動作を制御する
とともに、これらに荷重係数Gi,Gjを供給する
ようになつている。
画像データバスを介して背景処理装置23へ
送られた画像データは掛算器231に入力され、
荷重係数Giを掛けられてこの結果を減算器23
3へ出力する。
送られた画像データは掛算器231に入力され、
荷重係数Giを掛けられてこの結果を減算器23
3へ出力する。
他方、画像データバスを介して背景処理装置
23へ送られた画像データは掛算器232に入力
され荷重係数Gjを掛けられてこの結果を減算器
233へ出力する。減算器233は、掛算器23
1,232か送られた二つの画像データの同一空
間位置における画素間の濃度値の差を求めこの結
果を画像データバスを介して画像記憶装置25
へ送出する。
23へ送られた画像データは掛算器232に入力
され荷重係数Gjを掛けられてこの結果を減算器
233へ出力する。減算器233は、掛算器23
1,232か送られた二つの画像データの同一空
間位置における画素間の濃度値の差を求めこの結
果を画像データバスを介して画像記憶装置25
へ送出する。
機能画像作成装置24は、第4図に示すフーリ
エ級数の係数を求めるための時間フイルター装置
10と、第5図に示す振幅画像作成装置11及び
位相画像作成装置12とから構成され、画像処理
コントローラ20により制御されるようになつて
いる。
エ級数の係数を求めるための時間フイルター装置
10と、第5図に示す振幅画像作成装置11及び
位相画像作成装置12とから構成され、画像処理
コントローラ20により制御されるようになつて
いる。
時間フイルター装置10は掛算器101,10
2及び和算器103から構成され、掛算器101
には画像データバスを介して画像記憶装置25
内の所定のアドレス領域に記憶された画像データ
のゼロクリアされたものが入力される。
2及び和算器103から構成され、掛算器101
には画像データバスを介して画像記憶装置25
内の所定のアドレス領域に記憶された画像データ
のゼロクリアされたものが入力される。
また、掛算器102には画像データバスを介
して画像記憶装置25内の所定のアドレス領域に
記憶された画像データが入力される。
して画像記憶装置25内の所定のアドレス領域に
記憶された画像データが入力される。
掛算器101,102に入力されたそれぞれの
画像データには荷重係数Fi,Fjがそれぞれ掛けら
れ、この両者の結果が和算器103に入力されて
和が求められる。
画像データには荷重係数Fi,Fjがそれぞれ掛けら
れ、この両者の結果が和算器103に入力されて
和が求められる。
和算器103の出力は画像データバスを介し
て画像記憶装置25の所定のアドレス領域に順次
記憶される。
て画像記憶装置25の所定のアドレス領域に順次
記憶される。
このような動作を連続した一心拍分の画像につ
いて行なうことにより、下記(1)式で示される時間
フイルター画像が求まる。
いて行なうことにより、下記(1)式で示される時間
フイルター画像が求まる。
C=K
〓j=1
Fj・SI(j) ……(1)
(1)式において、荷重係数としてFj=2/K・
sin2π(j−1)/Kをとれば、正弦関数の係数a1
が計算され、また荷重係数としてFj=2/K・
cos2π(j−1)/Kをとれば余弦関数の係数b1が
計算される。但しここでKは1心拍周期間に収集
される画像枚数である。
sin2π(j−1)/Kをとれば、正弦関数の係数a1
が計算され、また荷重係数としてFj=2/K・
cos2π(j−1)/Kをとれば余弦関数の係数b1が
計算される。但しここでKは1心拍周期間に収集
される画像枚数である。
振幅画像作成装置11は、第5図に示すように
2つの2乗器111,112、加算器113、平
方根器114及びグレイレベル変換器115によ
り構成されている。
2つの2乗器111,112、加算器113、平
方根器114及びグレイレベル変換器115によ
り構成されている。
そして、画像記憶装置25から画像データバス
を介してa1画像(係数a1を有する画像)が2乗
器111に入力されてここで2乗され、また同様
に画像データバスを介して、b1画像(係数b1を
有する画像)が2乗器112に入力されてここで
2乗され、両者は加算器113で加算された後平
方根器114より平方根が求められ、さらにその
値がグレイレベル変換器15によりグレイレベル
変換されて振幅画像AIとなり、画像データバス
を介して画像記憶装置25の所定アドレス領域
に記憶される。
を介してa1画像(係数a1を有する画像)が2乗
器111に入力されてここで2乗され、また同様
に画像データバスを介して、b1画像(係数b1を
有する画像)が2乗器112に入力されてここで
2乗され、両者は加算器113で加算された後平
方根器114より平方根が求められ、さらにその
値がグレイレベル変換器15によりグレイレベル
変換されて振幅画像AIとなり、画像データバス
を介して画像記憶装置25の所定アドレス領域
に記憶される。
位相画像作成装置12は、第5図に示すように
逆数器121、掛算器122、逆タンジエント器
123及びグレイレベル変換器124により構成
されている。
逆数器121、掛算器122、逆タンジエント器
123及びグレイレベル変換器124により構成
されている。
そして、画像記憶装置25から画像データバス
を介してb1画像が掛算器122に、また同様に
画像データバスを介してa1画像が逆数器121
に入力されここでその逆数に変換された後掛算器
122にそれぞれ入力される。
を介してb1画像が掛算器122に、また同様に
画像データバスを介してa1画像が逆数器121
に入力されここでその逆数に変換された後掛算器
122にそれぞれ入力される。
掛算器122は両者の積をとりこの値は逆タン
ジエント器123により逆タンジエントされて位
相が求められ、次にグレイレベル変換器124に
よりグレイレベルに変換され位相画像PIとして
画像データバスを介して画像記憶装置25の所
定のアドレス領域に記憶される。
ジエント器123により逆タンジエントされて位
相が求められ、次にグレイレベル変換器124に
よりグレイレベルに変換され位相画像PIとして
画像データバスを介して画像記憶装置25の所
定のアドレス領域に記憶される。
前記セレクタ26は例えば切替スイツチで構成
され、画像データバス,を介して画像データ
が入力され、いずれか一系統の画像データを画像
処理コントローラ20から送出される制御信号R
に基づいて選択し、その選択結果を画像データバ
スを介して画像表示装置5へ送出するようにな
つている。
され、画像データバス,を介して画像データ
が入力され、いずれか一系統の画像データを画像
処理コントローラ20から送出される制御信号R
に基づいて選択し、その選択結果を画像データバ
スを介して画像表示装置5へ送出するようにな
つている。
次に上記構成のX線診断装置の作用を、造影左
心臓の振幅画像AI及び位相画像PIを作成表示す
る場合について説明する。
心臓の振幅画像AI及び位相画像PIを作成表示す
る場合について説明する。
最初に造影前画像MIiと造影画像CIiの収集手
順を説明する。
順を説明する。
システムコンソール7を操作して画像収集ルー
チンキーと検査部位キーとを押すと、システムコ
ントローラ6は所定の部位別の収集シーケンスを
選択し、X線曝射条件をX線発生制御器8へ送つ
てパラメータをセツトし、造影剤注入装置9をス
タンバイし、画像収集装置3、画像処理装置4及
び画像表示装置5にパラメータをセツトしこれら
をスタンバイ状態にする。
チンキーと検査部位キーとを押すと、システムコ
ントローラ6は所定の部位別の収集シーケンスを
選択し、X線曝射条件をX線発生制御器8へ送つ
てパラメータをセツトし、造影剤注入装置9をス
タンバイし、画像収集装置3、画像処理装置4及
び画像表示装置5にパラメータをセツトしこれら
をスタンバイ状態にする。
システムコンソール7上の収集開始キーを押す
と、上述した収集シーケンスで造影剤注入装置9
が駆動され造影剤を被検体2へ注入する。
と、上述した収集シーケンスで造影剤注入装置9
が駆動され造影剤を被検体2へ注入する。
次にX線発生制御器8を駆動してX線発生器1
からX線を発生させ、撮影部位に造影剤が出現す
る前の被検体2を透過させる。被検体2を透過し
たX線は画像収集装置3で検出され、電気信号に
変換されさらにA/D変換されて画像データバス
を介して画像処理装置4に送出される。
からX線を発生させ、撮影部位に造影剤が出現す
る前の被検体2を透過させる。被検体2を透過し
たX線は画像収集装置3で検出され、電気信号に
変換されさらにA/D変換されて画像データバス
を介して画像処理装置4に送出される。
このようにして得られた画像データはまず濃度
時間変化図作成処理装置21に入力され、所定の
関心領域内の濃度が計測された後順次画像番号順
に画像データバスを介して画像記憶装置25内
の所定のアドレス領域に格納されるとともに、画
像データバス、セレクタ26、画像データバス
を介して画像表示装置5に送られ収集画像とし
て画像表示される。
時間変化図作成処理装置21に入力され、所定の
関心領域内の濃度が計測された後順次画像番号順
に画像データバスを介して画像記憶装置25内
の所定のアドレス領域に格納されるとともに、画
像データバス、セレクタ26、画像データバス
を介して画像表示装置5に送られ収集画像とし
て画像表示される。
X線曝射が終了すると、それまでに求められて
いた関心領域内の濃度時間変化図も画像記憶装置
25内の所定のアドレス領域に記憶され、上述し
た収集画像の場合と同様画像表示装置5に送られ
て画像表示されると同時にこの濃度時間変化図は
画像データバスを介して心拍周期検出処理装置
22に入力される。
いた関心領域内の濃度時間変化図も画像記憶装置
25内の所定のアドレス領域に記憶され、上述し
た収集画像の場合と同様画像表示装置5に送られ
て画像表示されると同時にこの濃度時間変化図は
画像データバスを介して心拍周期検出処理装置
22に入力される。
心拍動周期検出処理装置22において心拍動の
心拍番号及び周期と画像番号間の対応がとられ、
この対応関係が制御線uを介して画像処理コント
ローラ20に伝送される。
心拍番号及び周期と画像番号間の対応がとられ、
この対応関係が制御線uを介して画像処理コント
ローラ20に伝送される。
画像処理コントローラ20は、心拍動周期検出
処理装置22から送られる心拍動情報と前記濃度
時間変化図作成処理装置21から送られる濃度情
報とに基づき制御線P3を介して背景処理装置2
3へ荷重係数G1,G2を送出し、画像記憶装置2
5へは制御線Q1〜QNを介して同一心拍位相を有
する造影前画像MIi及び造影画像CIiの画像番号
を送出する。
処理装置22から送られる心拍動情報と前記濃度
時間変化図作成処理装置21から送られる濃度情
報とに基づき制御線P3を介して背景処理装置2
3へ荷重係数G1,G2を送出し、画像記憶装置2
5へは制御線Q1〜QNを介して同一心拍位相を有
する造影前画像MIi及び造影画像CIiの画像番号
を送出する。
画像記憶装置25は、この2つの画像番号の画
像データを画像データバス,を介して順次背
景処理装置23へ送出する。
像データを画像データバス,を介して順次背
景処理装置23へ送出する。
背景処理装置23は、同一心拍位相iにおける
造影画像CIiと造影前画像MIiとの間で画像処理
を行ない第3図に示すように造影血液画像SIiを
作成する。
造影画像CIiと造影前画像MIiとの間で画像処理
を行ない第3図に示すように造影血液画像SIiを
作成する。
すなわち、心拍位相iにおける造影前画像
MIi、造影画像CIiに対し、背景処理装置23に
おいてSIi=MIi−CIiの演算処理が1心拍分のK
枚の画像に対して行なわれ、それぞれ造影剤が存
在しない背景が除去される。また、このようにし
て得られた造影血液画像SIiは画像データバス
を介して画像記憶装置25の所定のアドレス領域
に記憶される。
MIi、造影画像CIiに対し、背景処理装置23に
おいてSIi=MIi−CIiの演算処理が1心拍分のK
枚の画像に対して行なわれ、それぞれ造影剤が存
在しない背景が除去される。また、このようにし
て得られた造影血液画像SIiは画像データバス
を介して画像記憶装置25の所定のアドレス領域
に記憶される。
以上の処理は造影画像CIiが無くなるまで行わ
れる。
れる。
画像記憶装置25に記憶される造影血液画像
SIiは同時に画像データバス,のいずれかを
経てセレクタ26に送られ、さらに画像データバ
スを介して画像表示装置5へ転送されて画像表
示される。
SIiは同時に画像データバス,のいずれかを
経てセレクタ26に送られ、さらに画像データバ
スを介して画像表示装置5へ転送されて画像表
示される。
次に得られた造影血液画像SIiの1心拍分の画
像から心血管の振幅画像AIと位相画像PIを求め
る手順を第4図、第5図を参照して説明する。
像から心血管の振幅画像AIと位相画像PIを求め
る手順を第4図、第5図を参照して説明する。
1心拍分内の造影血液画像SIiをK枚、心拍周
期をT・K=1/f0とすれば、k枚目(但し、1
≦k≦K)の造影血液画像SIiの時刻はt=(k−
1)Tとなる。ここでTは画像のサンプル間隔時
間である。
期をT・K=1/f0とすれば、k枚目(但し、1
≦k≦K)の造影血液画像SIiの時刻はt=(k−
1)Tとなる。ここでTは画像のサンプル間隔時
間である。
また、(x,y)座標における画素濃度をPC
(x,y,t)とおき、これを下記(2)式のように
近似する。
(x,y,t)とおき、これを下記(2)式のように
近似する。
PC(x,y,t)−a0/2+a1cos2πt/kT+b1
sin2πt/kT ……(2) (2)式から係数a1,b1はそれぞれ以下の(3)、(4)式
で求めることができる。
sin2πt/kT ……(2) (2)式から係数a1,b1はそれぞれ以下の(3)、(4)式
で求めることができる。
a1=2/KK
〓K=1
PC{(x,y,(k−1)T}
・cos{2π(k−1)/K} ……(3)
b1=2/kK
〓K=1
PC{(x,y,(k−1)T}
・sin{2π(k−1)/K} ……(4)
ここで、A1=√1 2+1 2,P1=tan-1(b1/a1)
とおくと(2)式は下記(5)式となる。
とおくと(2)式は下記(5)式となる。
但しA1は振幅値を、P1は位相を示す。
PC(x,y,t)=a0/2+A1sin(2πt/KT+P1)
……(5) 次に係数a1,b1を求める手順について説明す
る。画像記憶装置25内の所定のアドレス領域の
データがゼロクリアされている。
……(5) 次に係数a1,b1を求める手順について説明す
る。画像記憶装置25内の所定のアドレス領域の
データがゼロクリアされている。
これらをa1画像、b1画像という。このa1画像は
画像データバスを介して、また造影血液画像
SIiの第1番目の画像は画像データバスを介し
てともに機能画像作成装置24の時間フイルター
装置10に入力れ、ゼロクリアされているa1画像
には荷重係数Fiとして1、第1番目の造影血液画
像SIiには荷重係数Fjとして2/Kcos2π(k−
1)/Kが掛けられた後、両者の和が求められ、
この結果が画像データバスを介して画像記憶装
置25内のa1画像を記憶するための記憶領域に記
憶される。
画像データバスを介して、また造影血液画像
SIiの第1番目の画像は画像データバスを介し
てともに機能画像作成装置24の時間フイルター
装置10に入力れ、ゼロクリアされているa1画像
には荷重係数Fiとして1、第1番目の造影血液画
像SIiには荷重係数Fjとして2/Kcos2π(k−
1)/Kが掛けられた後、両者の和が求められ、
この結果が画像データバスを介して画像記憶装
置25内のa1画像を記憶するための記憶領域に記
憶される。
このとき、a1画像として下記(6)式が記憶される
こととなる。
こととなる。
a1=2/Kcos2π(1−1)/K・PC{x,y,
(1−1)T} ……(6) 次にa1画像が画像データバスを介して第2番
目の造影血液画像SIiが画像データバスを介し
て時間フイルター装置10に入力され、上述した
場合と同様に荷重係数Fiとして1、荷重係数Fjと
して2/Kcos2π(2−1)/Kが掛けられた後両
者の和が求められ、この結果が画像データバス
を介して画像記憶装置25内の前記記憶領域に記
憶される。
(1−1)T} ……(6) 次にa1画像が画像データバスを介して第2番
目の造影血液画像SIiが画像データバスを介し
て時間フイルター装置10に入力され、上述した
場合と同様に荷重係数Fiとして1、荷重係数Fjと
して2/Kcos2π(2−1)/Kが掛けられた後両
者の和が求められ、この結果が画像データバス
を介して画像記憶装置25内の前記記憶領域に記
憶される。
このときのa1画像として下記(7)式が記憶される
ことになる。
ことになる。
a1(x,y)=K
〓K=1
2/Kcos2π(k−1)
/K・PC{x,y,(k−1)T} ……(7)
以下同様にしてK枚の造影血液画像SIiが時間
フイルター装置10により処理され、1心拍分の
a1画像として下記(8)式が求められる。
フイルター装置10により処理され、1心拍分の
a1画像として下記(8)式が求められる。
a1(x,y)=2/KK
〓K=1
cos2π(k−1)
/K・PC{x,y,(k−1)T} ……(8)
上述したa1画像の場合は同様な処理が荷重係数
Fj=2/Ksin2π(k−1)/Kを用いてb1画像に
対して行なわれ、1心拍分のb1画像として下記(9)
式が求められる。
Fj=2/Ksin2π(k−1)/Kを用いてb1画像に
対して行なわれ、1心拍分のb1画像として下記(9)
式が求められる。
b1(x,y)=2/KK
〓K=1
sin2π(k−1)
/K・PC{x,y,(k−1)T} ……(9)
次に上述のようにして求められた画像記憶装置
25内に記憶されたa1画像が画像データバスを
介して、b1画像が画像データバスを介してとも
に第5図に示す振幅画像作成装置11に入力さ
れ、ここで既述したような演算処理及びグレイレ
ベル変換が行われ1心拍分の造影血液画像SIiに
対する振幅画像AI(=√1 2+1 2)が求められる。
この結果は画像データバスを介して画像記憶装
置25の所定の記憶領域に記憶される。
25内に記憶されたa1画像が画像データバスを
介して、b1画像が画像データバスを介してとも
に第5図に示す振幅画像作成装置11に入力さ
れ、ここで既述したような演算処理及びグレイレ
ベル変換が行われ1心拍分の造影血液画像SIiに
対する振幅画像AI(=√1 2+1 2)が求められる。
この結果は画像データバスを介して画像記憶装
置25の所定の記憶領域に記憶される。
同様にa1画像、b1画像が位相画像作成装置12
に入力されここで1心拍分の位相画像PI(=tan-1
b1/a1)が求められ、この結果は画像記憶装置2
5の所定の記憶領域に記憶される。
に入力されここで1心拍分の位相画像PI(=tan-1
b1/a1)が求められ、この結果は画像記憶装置2
5の所定の記憶領域に記憶される。
この振幅画像AI及び位相画像PIは画像データ
バス,を介してセレクタ26へ送られ、さら
に画像データバスを介して画像表示装置5へ送
出され、画像表示される。
バス,を介してセレクタ26へ送られ、さら
に画像データバスを介して画像表示装置5へ送
出され、画像表示される。
本発明は上述した実施例に限定されるものでは
なく、その要旨の範囲内で種々の変形が可能であ
る。
なく、その要旨の範囲内で種々の変形が可能であ
る。
例えば上述した実施例では画像処理装置の心拍
動周期検出処理装置に対する入力データは画像デ
ータの場合について説明したが、心電図を用いる
こともできる。
動周期検出処理装置に対する入力データは画像デ
ータの場合について説明したが、心電図を用いる
こともできる。
また、振幅画像、位相画像の作成に1次フーリ
エ変換を用いる場合について説明したが、機能画
像作成装置における荷重係数を高次の正弦波、余
弦波とすることにより高次のフーリエ変換による
両画像の作成を行なうことができる。
エ変換を用いる場合について説明したが、機能画
像作成装置における荷重係数を高次の正弦波、余
弦波とすることにより高次のフーリエ変換による
両画像の作成を行なうことができる。
さらに、上述した実施例においては、画像処理
装置を濃度時間変化図作成処理装置、心拍動周期
検出処理装置、背景処理装置及び機能画像作成装
置により構成した場合について説明したが、1心
拍分のX線透過画像に対し直接機能画像作成装置
により時間フイルター処理を行ない振幅画像、位
相画像を得るようにしてもよい。
装置を濃度時間変化図作成処理装置、心拍動周期
検出処理装置、背景処理装置及び機能画像作成装
置により構成した場合について説明したが、1心
拍分のX線透過画像に対し直接機能画像作成装置
により時間フイルター処理を行ない振幅画像、位
相画像を得るようにしてもよい。
[発明の効果]
以上詳述した本発明によれば、X線透過画像の
1心拍分を利用して機能画像作成装置により心血
管の振幅画像及び位相画像を得てこれを高速に表
示することが可能となり、心血管の機能を定量的
に計測、判断することができるX線診断装置を提
供し得るものである。
1心拍分を利用して機能画像作成装置により心血
管の振幅画像及び位相画像を得てこれを高速に表
示することが可能となり、心血管の機能を定量的
に計測、判断することができるX線診断装置を提
供し得るものである。
また、被検体の心拍動に対応して収集した造影
前画像と造影画像との間で背景処理を行なうもの
であるため、心拍動の影響による偽像がない高画
質の造影血液画像を表示することができるX線診
断装置を提供し得るものである。
前画像と造影画像との間で背景処理を行なうもの
であるため、心拍動の影響による偽像がない高画
質の造影血液画像を表示することができるX線診
断装置を提供し得るものである。
第1図は本発明の実施例を示すブロツク図、第
2図は第1図に示す装置を構成する画像処理装置
のブロツク図、第3図は同上の背景処理装置の構
成及び作用を示す説明図、第4図は同上の機能画
像作成装置の一部である時間フイルター装置の構
成及び作用を示す説明図、第5図は同上の機能画
像作成装置の一部である振幅画像作成装置及び位
相画像作成装置の構成及び作用を示す説明図であ
る。 1……X線発生器、2……被検体、3……画像
収集装置、4……画像処理装置、5……画像表示
装置、6……システムコントローラ、7……シス
テムコンソール、8……X線発生制御器、9……
造影剤注入装置、10……時間フイルター装置、
20……画像処理コントローラ、21……濃度時
間変化図作成処理装置、22……心拍動周期検出
処理装置、23……背景処理装置、24……機能
画像作成装置。
2図は第1図に示す装置を構成する画像処理装置
のブロツク図、第3図は同上の背景処理装置の構
成及び作用を示す説明図、第4図は同上の機能画
像作成装置の一部である時間フイルター装置の構
成及び作用を示す説明図、第5図は同上の機能画
像作成装置の一部である振幅画像作成装置及び位
相画像作成装置の構成及び作用を示す説明図であ
る。 1……X線発生器、2……被検体、3……画像
収集装置、4……画像処理装置、5……画像表示
装置、6……システムコントローラ、7……シス
テムコンソール、8……X線発生制御器、9……
造影剤注入装置、10……時間フイルター装置、
20……画像処理コントローラ、21……濃度時
間変化図作成処理装置、22……心拍動周期検出
処理装置、23……背景処理装置、24……機能
画像作成装置。
Claims (1)
- 1 X線を被写体に曝射するX線発生器と、被写
体を透過したX線からX線透過画像を収集する画
像収集装置と、X線透過画像を画像処理する画像
処理装置と、この画像処理装置による画像処理後
の画像を表示する画像表示装置とを有するX線診
断装置において、前記画像処理装置は少なくとも
1心拍分の血管造影前のX線透過画像及び血管造
影後のX線造影画像から少なくとも1心拍分の造
影血管画像を得る手段と、この手段によつて得ら
れた少なくとも1心拍分の造影血管画像それぞれ
に所定の重み付けをして2つの画像を得る手段
と、この手段によつて得られた2つの画像の相乗
平均をとることにより振幅画像を作成する手段
と、前記2つの画像の比の逆正接をとることによ
り位相画像を作成する手段とを備えたことを特徴
とするX線診断装置。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59006266A JPS60150729A (ja) | 1984-01-19 | 1984-01-19 | X線診断装置 |
| US06/693,168 US4689670A (en) | 1984-01-19 | 1985-01-22 | Method and apparatus for measuring a biological function of a subject |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59006266A JPS60150729A (ja) | 1984-01-19 | 1984-01-19 | X線診断装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS60150729A JPS60150729A (ja) | 1985-08-08 |
| JPH0429381B2 true JPH0429381B2 (ja) | 1992-05-18 |
Family
ID=11633636
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP59006266A Granted JPS60150729A (ja) | 1984-01-19 | 1984-01-19 | X線診断装置 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US4689670A (ja) |
| JP (1) | JPS60150729A (ja) |
Families Citing this family (39)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0785B2 (ja) * | 1986-04-10 | 1995-01-11 | 株式会社東芝 | 肺換気機能診断装置 |
| DE3814089A1 (de) * | 1987-04-27 | 1988-11-17 | Toshiba Kawasaki Kk | Verfahren und anordnung zum wiedergeben von reihenbildern |
| JPS6458243A (en) * | 1987-08-28 | 1989-03-06 | Toshiba Corp | X-ray image processing apparatus |
| US5008907A (en) * | 1989-05-31 | 1991-04-16 | The Regents Of The University Of California | Therapy x-ray scanner |
| DE69425158T2 (de) * | 1993-02-26 | 2001-03-15 | E-Y Laboratories, Inc. | Optisches probe-analysesystem und verfahren |
| US5840026A (en) * | 1994-09-21 | 1998-11-24 | Medrad, Inc. | Patient specific dosing contrast delivery systems and methods |
| US6397098B1 (en) * | 1994-09-21 | 2002-05-28 | Medrad, Inc. | Data communication and control for medical imaging systems |
| US5533085A (en) * | 1995-02-27 | 1996-07-02 | University Of Washington | Automatic indexing of cine-angiograms |
| WO1999040536A1 (en) * | 1998-02-10 | 1999-08-12 | Ey Laboratories, Inc. | Reflectometry system with compensation for specimen holder topography and with lock-rejection of system noise |
| DE19824496A1 (de) * | 1998-06-02 | 1999-12-09 | Siemens Ag | Vorrichtung mit einem Betriebsparameter, dessen Wert über eine Steuereinheit vorwählbar ist |
| FR2782629B1 (fr) * | 1998-08-31 | 2001-01-12 | Ge Medical Syst Sa | Procede d'amelioration de la qualite d'une image fluoroscopique visualisee de facon permanente a l'issue d'une sequence d'acquisition d'images en fluoroscopie cardiaque pulsee |
| US6243437B1 (en) | 1998-11-25 | 2001-06-05 | General Electric Company | Coronary calcification detection using retrospective cardiac gating of imaging system |
| US6275560B1 (en) * | 1998-12-22 | 2001-08-14 | General Electric Company | Cardiac gated computed tomography system |
| US6256368B1 (en) * | 1999-10-15 | 2001-07-03 | General Electric Company | Methods and apparatus for scout-based cardiac calcification scoring |
| JP3697233B2 (ja) * | 2002-04-03 | 2005-09-21 | キヤノン株式会社 | 放射線画像処理方法及び放射線画像処理装置 |
| US6990368B2 (en) | 2002-04-04 | 2006-01-24 | Surgical Navigation Technologies, Inc. | Method and apparatus for virtual digital subtraction angiography |
| CN102488527B (zh) | 2004-11-16 | 2015-06-10 | 拜耳医疗保健公司 | 从成像过程确定病人生理参数的系统和方法 |
| EP2392379A3 (en) | 2004-11-24 | 2012-03-21 | Medrad, Inc. | Devices, systems and methods for fluid delivery |
| US7940971B2 (en) * | 2006-07-24 | 2011-05-10 | Siemens Medical Solutions Usa, Inc. | System and method for coronary digital subtraction angiography |
| EP2097835B1 (en) | 2006-12-29 | 2018-05-30 | Bayer Healthcare LLC | Patient-based parameter generation systems for medical injection procedures |
| RU2009139916A (ru) * | 2007-03-29 | 2011-05-10 | Конинклейке Филипс Электроникс Н.В. (Nl) | Способ и устройство для получения объединенных рентгеновских изображений |
| US7680240B2 (en) * | 2007-03-30 | 2010-03-16 | General Electric Company | Iterative reconstruction of tomographic image data method and system |
| US7920670B2 (en) * | 2007-03-30 | 2011-04-05 | General Electric Company | Keyhole computed tomography |
| US8428694B2 (en) | 2007-07-17 | 2013-04-23 | Medrad, Inc. | Methods for determination of parameters for a procedure, for estimation of cardiopulmonary function and for fluid delivery |
| US9421330B2 (en) | 2008-11-03 | 2016-08-23 | Bayer Healthcare Llc | Mitigation of contrast-induced nephropathy |
| JP6025718B2 (ja) | 2010-06-24 | 2016-11-16 | ベイヤー メディカル ケア インク. | 注入プロトコルのための医薬品の伝播およびパラメータ生成のモデル |
| KR101805624B1 (ko) * | 2011-08-29 | 2017-12-08 | 삼성전자주식회사 | 장기 모델 영상 생성 방법 및 장치 |
| RU2605272C2 (ru) | 2012-05-14 | 2016-12-20 | БАЙЕР ХелсКер ЛЛСи | Системы и способы определения протоколов инъекции фармацевтической жидкости исходя из напряжения на рентгеновской трубке |
| US9555379B2 (en) | 2013-03-13 | 2017-01-31 | Bayer Healthcare Llc | Fluid path set with turbulent mixing chamber, backflow compensator |
| US10898638B2 (en) | 2016-03-03 | 2021-01-26 | Bayer Healthcare Llc | System and method for improved fluid delivery in multi-fluid injector systems |
| JP7191821B2 (ja) | 2016-11-14 | 2022-12-19 | バイエル・ヘルスケア・エルエルシー | 医用流体送出に使用されるシリンジの内容物を確認するための方法およびシステム |
| WO2019046260A1 (en) | 2017-08-31 | 2019-03-07 | Bayer Healthcare Llc | METHOD FOR CONTROLLING DYNAMIC PRESSURE IN A FLUID INJECTOR SYSTEM |
| CA3066780A1 (en) | 2017-08-31 | 2019-03-07 | Bayer Healthcare Llc | Fluid injector system volume compensation system and method |
| JP7252143B2 (ja) | 2017-08-31 | 2023-04-04 | バイエル・ヘルスケア・エルエルシー | 駆動部材の位置及び流体注入器システムの機械的較正のためのシステム及び方法 |
| AU2018323442B2 (en) | 2017-08-31 | 2024-06-27 | Bayer Healthcare Llc | Fluid path impedance assessment for improving fluid delivery performance |
| JP7221885B2 (ja) | 2017-08-31 | 2023-02-14 | バイエル・ヘルスケア・エルエルシー | 注入器圧力較正システムおよび方法 |
| WO2019204605A1 (en) | 2018-04-19 | 2019-10-24 | Bayer Healthcare Llc | System and method for air detection in fluid injector |
| CA3111209A1 (en) | 2018-08-28 | 2020-03-05 | Bayer Healthcare Llc | Fluid injector system with improved ratio performance |
| CA3110859A1 (en) | 2018-08-28 | 2020-03-05 | Bayer Healthcare Llc | Fluid injector system, method of preventing fluid backflow, and computer program product |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS4972982A (ja) * | 1972-11-20 | 1974-07-15 | ||
| DE2409173A1 (de) * | 1974-02-26 | 1975-09-04 | Philips Patentverwaltung | Verfahren zur naeherungsweisen bestimmung der kontur eines pulsierenden objektes |
| US4204225A (en) * | 1978-05-16 | 1980-05-20 | Wisconsin Alumni Research Foundation | Real-time digital X-ray subtraction imaging |
| DE3124583A1 (de) * | 1981-06-23 | 1982-12-30 | Siemens AG, 1000 Berlin und 8000 München | Roentgendiagnostikeinrichtung zur darstellung von organen |
| JPS5815842A (ja) * | 1981-07-20 | 1983-01-29 | 株式会社東芝 | 画像輪郭抽出装置 |
| NL8202450A (nl) * | 1982-06-17 | 1984-01-16 | Philips Nv | Werkwijze en inrichting voor het verwerken van roentgenbeelden. |
-
1984
- 1984-01-19 JP JP59006266A patent/JPS60150729A/ja active Granted
-
1985
- 1985-01-22 US US06/693,168 patent/US4689670A/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPS60150729A (ja) | 1985-08-08 |
| US4689670A (en) | 1987-08-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH0429381B2 (ja) | ||
| US4903705A (en) | Digital radiography apparatus | |
| US6535570B2 (en) | Method for tracing organ motion and removing artifacts for computed tomography imaging systems | |
| EP1685538B1 (en) | Device and method for generating a three-dimensional vascular model | |
| US4739766A (en) | NMR blood vessel imaging method and apparatus | |
| JP2646663B2 (ja) | 動体イメージング方法およびその装置 | |
| JPH0693890B2 (ja) | 超音波診断装置 | |
| JP4086309B2 (ja) | 造影剤注入プロトコル決定方法および造影剤注入プロトコル演算装置 | |
| JPH11342132A (ja) | 血管計測方法および装置並びに医用画像装置 | |
| EP1330186B1 (en) | Method and apparatus for 3d-rotational x-ray imaging | |
| JP4316017B2 (ja) | X線ct装置 | |
| JPH0479257B2 (ja) | ||
| CN100443056C (zh) | 超声波多普勒血流测定装置 | |
| JPH119589A (ja) | X線ct装置及び画像再構成方法 | |
| JPH02164352A (ja) | 超音波診断画像処理装置 | |
| JP2007216043A (ja) | X線ct装置 | |
| JPH0432658B2 (ja) | ||
| JPS60236594A (ja) | X線診断装置 | |
| JPS60122547A (ja) | 放射線診断装置 | |
| JPH0321228A (ja) | 心筋潅流解析画像の表示方法及びその装置 | |
| CN121196597A (zh) | 脉冲多普勒填充数据确定方法及装置、电子设备 | |
| Cortesi et al. | 4D Flow MR Imaging for Blood Velocity Field Validation: A Preliminary Study in Atrial Fibrillation Patients | |
| JPH06125499A (ja) | ディジタル・サブトラクション・アンギオグラフィ装置 | |
| JP3573522B2 (ja) | 超音波診断装置 | |
| JPH0236097B2 (ja) |