JP3105520B2 - Ｘ線検査装置 - Google Patents
Ｘ線検査装置Info
- Publication number
- JP3105520B2 JP3105520B2 JP02133715A JP13371590A JP3105520B2 JP 3105520 B2 JP3105520 B2 JP 3105520B2 JP 02133715 A JP02133715 A JP 02133715A JP 13371590 A JP13371590 A JP 13371590A JP 3105520 B2 JP3105520 B2 JP 3105520B2
- Authority
- JP
- Japan
- Prior art keywords
- ray
- recording pulse
- time
- synchronization
- recording
- 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 - Fee Related
Links
- 238000007689 inspection Methods 0.000 title claims description 10
- 230000033001 locomotion Effects 0.000 claims description 21
- 238000000034 method Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 12
- 239000013598 vector Substances 0.000 claims description 5
- 230000001678 irradiating effect Effects 0.000 claims description 4
- 230000000737 periodic effect Effects 0.000 claims 5
- 230000033764 rhythmic process Effects 0.000 description 21
- 230000015654 memory Effects 0.000 description 10
- 239000002872 contrast media Substances 0.000 description 8
- 239000011159 matrix material Substances 0.000 description 6
- 230000009897 systematic effect Effects 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000005309 stochastic process Methods 0.000 description 3
- 230000002238 attenuated effect Effects 0.000 description 2
- 210000004204 blood vessel Anatomy 0.000 description 2
- 239000013256 coordination polymer Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 208000009729 Ventricular Premature Complexes Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000010247 heart contraction Effects 0.000 description 1
- XMBWDFGMSWQBCA-UHFFFAOYSA-N hydrogen iodide Chemical compound I XMBWDFGMSWQBCA-UHFFFAOYSA-N 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000011282 treatment Methods 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/48—Diagnostic techniques
- A61B6/481—Diagnostic techniques involving the use of contrast agents
-
- 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/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
-
- 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
-
- H—ELECTRICITY
- H05—ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
- H05G—X-RAY TECHNIQUE
- H05G1/00—X-ray apparatus involving X-ray tubes; Circuits therefor
- H05G1/08—Electrical details
- H05G1/60—Circuit arrangements for obtaining a series of X-ray photographs or for X-ray cinematography
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Vascular Medicine (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- X-Ray Techniques (AREA)
Description
【発明の詳細な説明】 本発明は、周期的に動く対象の時系列Ｘ線画像を作る
Ｘ線検査装置に係る。該対象はＸ線源により出射された
Ｘ線ビームにより照射される。記録装置は対象の動作周
期を測定する。記録パルスはＸ線源を活性化及び非活性
化する同期装置に印加される。
Ｘ線検査装置に係る。該対象はＸ線源により出射された
Ｘ線ビームにより照射される。記録装置は対象の動作周
期を測定する。記録パルスはＸ線源を活性化及び非活性
化する同期装置に印加される。
本発明は又Ｘ線検査装置を実施するに適したＸ線源
と、記録装置と、同期装置とからなる装置に係る。Ｘ線
検査装置は、リーフレット「フィリップス・メディカル
・システム」4522 984 09241/744、DVI−Ｖ、1985年１
月から公知である。
と、記録装置と、同期装置とからなる装置に係る。Ｘ線
検査装置は、リーフレット「フィリップス・メディカル
・システム」4522 984 09241/744、DVI−Ｖ、1985年１
月から公知である。
例えば人体、心臓又は血管の如き対象のＸ線画像の形
成中、Ｘ線ビームは対象及びその環境で減衰される。人
体では、画像されるべき対象及びその環境は主に水から
なり、これによりＸ線ビームは対象及びその環境により
実質的に同じ程度まで減衰される。これは対象の画像を
低コントラストにする。画像さるべき対象がヨー化物の
ような比較的高いＸ線吸収を有するコントラストエージ
ェントで満たし、コントラストエージェントなしで記録
された対象の像即ちいわゆるマスク像からコントラスト
エージェントで満たされた対象の像を減算することによ
り、コントラスト強調が達成される。コントラストエー
ジェントで満たされた対象の画像が対象の動きによりマ
スク画像に関してシフトされた場合、人工物がコントラ
スト強調画像中に発生する。テレビジョンモニターでの
表示する以前に雑音積分用に合計され、又は画像さるべ
き対象を通るコントラストエージェントの流れに従う一
連の連続的Ｘ線画像の形成中明らかにマスク画像からの
減算によるコントラスト強調は対象の動きにより強く低
下される。心臓の近傍にある血管の画像に対し、心臓の
動きにより該問題が生じる。記録装置、例えば心電図の
形成に対し、記録装置により対象の動作リズムを測定す
ることにより照射は動きのリズムと同期されうる。これ
は例えばフィリップスDVI−Ｖシステムで起こる。照射
が記録装置により供給された記録パルスに関して所定の
遅延の後に実行されるので、Ｘ線画像は相が合っており
マスク画像と整合する。これは例えばコントラストエー
ジェントにより起こされた画像されるべき対象の動きの
リズムの変動の結果、Ｘ線画像が異なる位置に対象を示
すという欠点を有し、これにより人工物は、このＸ線検
査装置が用いられた時更に発生する。
成中、Ｘ線ビームは対象及びその環境で減衰される。人
体では、画像されるべき対象及びその環境は主に水から
なり、これによりＸ線ビームは対象及びその環境により
実質的に同じ程度まで減衰される。これは対象の画像を
低コントラストにする。画像さるべき対象がヨー化物の
ような比較的高いＸ線吸収を有するコントラストエージ
ェントで満たし、コントラストエージェントなしで記録
された対象の像即ちいわゆるマスク像からコントラスト
エージェントで満たされた対象の像を減算することによ
り、コントラスト強調が達成される。コントラストエー
ジェントで満たされた対象の画像が対象の動きによりマ
スク画像に関してシフトされた場合、人工物がコントラ
スト強調画像中に発生する。テレビジョンモニターでの
表示する以前に雑音積分用に合計され、又は画像さるべ
き対象を通るコントラストエージェントの流れに従う一
連の連続的Ｘ線画像の形成中明らかにマスク画像からの
減算によるコントラスト強調は対象の動きにより強く低
下される。心臓の近傍にある血管の画像に対し、心臓の
動きにより該問題が生じる。記録装置、例えば心電図の
形成に対し、記録装置により対象の動作リズムを測定す
ることにより照射は動きのリズムと同期されうる。これ
は例えばフィリップスDVI−Ｖシステムで起こる。照射
が記録装置により供給された記録パルスに関して所定の
遅延の後に実行されるので、Ｘ線画像は相が合っており
マスク画像と整合する。これは例えばコントラストエー
ジェントにより起こされた画像されるべき対象の動きの
リズムの変動の結果、Ｘ線画像が異なる位置に対象を示
すという欠点を有し、これにより人工物は、このＸ線検
査装置が用いられた時更に発生する。
本発明は、なかんずくその対象に対し、動作人工物の
発生が減少するＸ線検査装置を提供する。
発生が減少するＸ線検査装置を提供する。
これを達成する為、本発明によるＸ線検査装置は次の
記録パルスの時点から同期装置の算出手段により記録パ
ルスから計算されることを特徴とする。
記録パルスの時点から同期装置の算出手段により記録パ
ルスから計算されることを特徴とする。
測定された記録パルスからの次の記録パルスの時点の
計算は２つの順次の記録パルス間の時間の一定の時点
で、即ち一定の対象位置で照射が起こるのを可能にす
る。動きのリズムの変化には計算により検出され、Ｘ線
源の同期は変化した動きのリズムに適合される。時間的
に遠い記録パルスの効果が計算でより強く又はより弱く
なるにつれ、同期はリズム変化により速く又はより遅く
応答する。対象動作が系統的増速及び減速、例えば別な
心臓収縮又は心臓の早期心臓収縮に加えて確率的動作を
示す時、記録パルスは正確に予測可能な時点で互いに続
かない。測定された記録パルスの組み合わせにより次の
パルス（これから発生するパルス）の正確な計算のため
には、同期が系統的リズム変化に無感応になることな
く、全ての変動のうちの幾つかが同期の変化を引き起こ
すよう先行パルスが正確に加重されることが必要であ
る。平均動作リズム及び変化が対象動作を示す確率過程
の特徴的パラメータである時、次の記録パルスの正確な
予測は、平均値及び偏差値が知られている時、行われう
る。通常、対象動作の拡大統計的分析に対し時間がない
が、多数のＸ線画像はコントラストエージェントの変位
により素速く形成されなければならない。計算は記録パ
ルスの平均周期内の次の記録パルス時点に対する値を発
生すべきで、これにより僅かな計算時間だけを必要とす
る。次の記録パルスの計算はＸ線画像の形成を簡素化す
る為、完全に自動的であるべきである。
計算は２つの順次の記録パルス間の時間の一定の時点
で、即ち一定の対象位置で照射が起こるのを可能にす
る。動きのリズムの変化には計算により検出され、Ｘ線
源の同期は変化した動きのリズムに適合される。時間的
に遠い記録パルスの効果が計算でより強く又はより弱く
なるにつれ、同期はリズム変化により速く又はより遅く
応答する。対象動作が系統的増速及び減速、例えば別な
心臓収縮又は心臓の早期心臓収縮に加えて確率的動作を
示す時、記録パルスは正確に予測可能な時点で互いに続
かない。測定された記録パルスの組み合わせにより次の
パルス（これから発生するパルス）の正確な計算のため
には、同期が系統的リズム変化に無感応になることな
く、全ての変動のうちの幾つかが同期の変化を引き起こ
すよう先行パルスが正確に加重されることが必要であ
る。平均動作リズム及び変化が対象動作を示す確率過程
の特徴的パラメータである時、次の記録パルスの正確な
予測は、平均値及び偏差値が知られている時、行われう
る。通常、対象動作の拡大統計的分析に対し時間がない
が、多数のＸ線画像はコントラストエージェントの変位
により素速く形成されなければならない。計算は記録パ
ルスの平均周期内の次の記録パルス時点に対する値を発
生すべきで、これにより僅かな計算時間だけを必要とす
る。次の記録パルスの計算はＸ線画像の形成を簡素化す
る為、完全に自動的であるべきである。
本発明に関するＸ線検査装置の望ましい例は反復的カ
ルマンフィルタ法が次の記録パルスを計算するのに用い
られることを特徴とする。
ルマンフィルタ法が次の記録パルスを計算するのに用い
られることを特徴とする。
カルマンフィルタ法を使用すると、動きのリズムが反
復的に示され、これはそれ自身この確率的過程の変化に
適合する。この計算方法は系統的変化から動きのリズム
の非系統的変化を適切に区別する。これは明らかに心臓
の血管検査に対し有利である。
復的に示され、これはそれ自身この確率的過程の変化に
適合する。この計算方法は系統的変化から動きのリズム
の非系統的変化を適切に区別する。これは明らかに心臓
の血管検査に対し有利である。
以下図面と共に本発明を詳細に説明する。
第１図は細部４を含む対象３をＸ線ビーム５により照
射するＸ線源１を示す。Ｘ線検出器７、この場合はＸ線
増強管を用いて、画像搬送Ｘ線ビーム６はＸ線画像増強
管の出射スクリーンに現れる光画像に変換され、この光
画像はテレビジョンカメラ管９により電気信号に変換さ
れる。この信号はアナログ・ディジタル変換器11に印加
され、次にメモリー12又は13の一つに印加される。例え
ばメモリー12では、コントラストエイジェントで満たさ
れた細部４を含む対象３の像が記憶され、メモリー13で
は対象３のマスク像が記憶される。メモリー12及び13の
信号の対数化の後、メモリーのディジタル像は減算装置
15で減算され、細部４のコントラスト強調像を発生す
る。ディジタルーアナログ変換器16によるディジタルー
アナログ変換の後、この像はモニター17上に表示され
る。メモリー12及び13に記憶されたディジタルＸ線画像
はバックグラウンドメモリー19に記憶され、そこからデ
ィジタル画像がメモリー12及び13に印加される、細部４
の動作リズムを決める為、記録装置20は細部４に接続さ
れ、該記録装置は細部４が所定位置を占有する時点でパ
ルス状信号を同期装置22に印加する。同期装置22は記録
パルスから次の記録パルスの時点を計算し、細部４の系
統的リズム変化を検出する。同期装置22は対象のリアル
タイム像を表示するデレビジョンモニター23のフレーム
周期、例えば40msと同期してＸ線源１を作動させる。
射するＸ線源１を示す。Ｘ線検出器７、この場合はＸ線
増強管を用いて、画像搬送Ｘ線ビーム６はＸ線画像増強
管の出射スクリーンに現れる光画像に変換され、この光
画像はテレビジョンカメラ管９により電気信号に変換さ
れる。この信号はアナログ・ディジタル変換器11に印加
され、次にメモリー12又は13の一つに印加される。例え
ばメモリー12では、コントラストエイジェントで満たさ
れた細部４を含む対象３の像が記憶され、メモリー13で
は対象３のマスク像が記憶される。メモリー12及び13の
信号の対数化の後、メモリーのディジタル像は減算装置
15で減算され、細部４のコントラスト強調像を発生す
る。ディジタルーアナログ変換器16によるディジタルー
アナログ変換の後、この像はモニター17上に表示され
る。メモリー12及び13に記憶されたディジタルＸ線画像
はバックグラウンドメモリー19に記憶され、そこからデ
ィジタル画像がメモリー12及び13に印加される、細部４
の動作リズムを決める為、記録装置20は細部４に接続さ
れ、該記録装置は細部４が所定位置を占有する時点でパ
ルス状信号を同期装置22に印加する。同期装置22は記録
パルスから次の記録パルスの時点を計算し、細部４の系
統的リズム変化を検出する。同期装置22は対象のリアル
タイム像を表示するデレビジョンモニター23のフレーム
周期、例えば40msと同期してＸ線源１を作動させる。
第2a図は互いに間隔Ｔで続き系統的リズム変化が時点
t3以後発生する一連の記録パルスt1−t5を示す。記録パ
ルス関して所定の時間の間隔遅延した時点x1−x4での同
期装置によるＸ線源の同期は、リズム変化の後Ｔの分数
からずれるＴの分数でのリズム変化に先行する時点で起
こる。次の時点Ｔを計算することにより、また前回の記
録パルスの後Ｔの所定の分数で同期に作用させることに
より、上記問題が避けられる。第26図は心臓リズムで起
こる確率的リズム変化：別な収縮を示す。その場合には
別な心臓収縮は間隔Ｔで互いに続く２つの心臓鼓動間で
生じる。前回の記録パルスに関して予め決められた遅延
時間の後の同期は時点x2及びx3で位相外れ照射に導く。
第2c図は、変動が計算された次の記録パルスにかかる効
果を有さない場合、位相外れ照射は時点x2でのみ起こる
ことを示す。記録パルスt3を確立的変動として認識する
ことにより、時点x2での同期はそれ（記録パルスt3）を
基に抑制されうる。心臓リズムで起こる別な変動はいわ
ゆる早期心室収縮であり、ここで、一つの記録パルスは
早く起こり、次の記録パルスは予期した時点で起こる。
位相外れ照射が公知の同期手段によって時点x1及びx2で
起こる場合位相外れ照射は、本発明による同期が用いら
れている時、専ら時点x1で起こる。時点x2での周期は確
率的変動として記録パルスt2の認識の後抑圧されうる。
t3以後発生する一連の記録パルスt1−t5を示す。記録パ
ルス関して所定の時間の間隔遅延した時点x1−x4での同
期装置によるＸ線源の同期は、リズム変化の後Ｔの分数
からずれるＴの分数でのリズム変化に先行する時点で起
こる。次の時点Ｔを計算することにより、また前回の記
録パルスの後Ｔの所定の分数で同期に作用させることに
より、上記問題が避けられる。第26図は心臓リズムで起
こる確率的リズム変化：別な収縮を示す。その場合には
別な心臓収縮は間隔Ｔで互いに続く２つの心臓鼓動間で
生じる。前回の記録パルスに関して予め決められた遅延
時間の後の同期は時点x2及びx3で位相外れ照射に導く。
第2c図は、変動が計算された次の記録パルスにかかる効
果を有さない場合、位相外れ照射は時点x2でのみ起こる
ことを示す。記録パルスt3を確立的変動として認識する
ことにより、時点x2での同期はそれ（記録パルスt3）を
基に抑制されうる。心臓リズムで起こる別な変動はいわ
ゆる早期心室収縮であり、ここで、一つの記録パルスは
早く起こり、次の記録パルスは予期した時点で起こる。
位相外れ照射が公知の同期手段によって時点x1及びx2で
起こる場合位相外れ照射は、本発明による同期が用いら
れている時、専ら時点x1で起こる。時点x2での周期は確
率的変動として記録パルスt2の認識の後抑圧されうる。
第３図はカルマンフィルタ法に対するフローチャート
を示す。画像されるべき対象の動作リズムは確率的処理
として示され、その三次元ベクトルｘ（ｔ）は処理に関
連したパラメータからなる。処理を適宜に特徴づけるx1
（tn）＝tn、 x2（tn）＝tn−tn−１及びx3（tn）＝tn−2tn−１＋t
n−２が見出された。処理は下記の式で示される： ｘ（ｋ＋１）＝Ax（ｋ）＋ｗ（ｋ） （１） ここで、バーは列ベクトルを示し大文字は行列を示す
ｗ（Ｋ）は平均ゼロを有する確率的白処理して示され
る： Ｅ（ｘ（ｋ）.wT（１）］＝Ｑ（ｋ）δk,1 （２） ここで、Ｅは予測演算子であり、ｗ（ｋ）Ｔは列ベク
トルｗ（Ｋ）の転置ベクトルである。
を示す。画像されるべき対象の動作リズムは確率的処理
として示され、その三次元ベクトルｘ（ｔ）は処理に関
連したパラメータからなる。処理を適宜に特徴づけるx1
（tn）＝tn、 x2（tn）＝tn−tn−１及びx3（tn）＝tn−2tn−１＋t
n−２が見出された。処理は下記の式で示される： ｘ（ｋ＋１）＝Ax（ｋ）＋ｗ（ｋ） （１） ここで、バーは列ベクトルを示し大文字は行列を示す
ｗ（Ｋ）は平均ゼロを有する確率的白処理して示され
る： Ｅ（ｘ（ｋ）.wT（１）］＝Ｑ（ｋ）δk,1 （２） ここで、Ｅは予測演算子であり、ｗ（ｋ）Ｔは列ベク
トルｗ（Ｋ）の転置ベクトルである。
式１に記述されるような一次自己回帰処理は一連の離
散測定置ｙ（ｋ）、すなわち、記録パルスとして観測さ
れ、これは以下の式によりシステムに関連する： ｙ（ｋ）＝Cx（ｋ）＋ｖ（ｋ） （３） ここで、Ｃ＝（1,0,0）及びｖ（Ｋ）は平均ゼロ及び
共分散行列を有する測定置を妨害する白色雑音処理であ
る： Ｅ［ｖ（ｋ）.v（１）］＝Ｒ（ｋ）δk,1 （４） 更に、ｗ（ｋ）及びｖ（ｋ）は非相関である。解決さ
るべき問題点は、動作リズムｘ（ｋ＋１）が、観測され
た記録パルスｙ（ｏ）…ｙ（ｋ）から計算されなければ
ならないことからなる。この予測子は（ｋ＋1|k）と
して示される。ここで （ｋ＋1|k）＝Ｅ［ｘ（ｋ＋1|y（ｋ）］は所与のｙ
（ｋ）に対するｘ（ｋ＋１）の条件平均値である。ｘ
（ｋ）及びｙ（ｋ）が共に正規分布を有するので、動作
リズムの統計値はこの分布の条件平均値及び分散により
完全に説明される。共通条件確率分布はカルマンフィル
タ法により反復的に与えられる。カルマンフィルタ法は
予測された次の記録パルス Ｅ［（tn+1−tn）２］の平均二乗誤差を最小にする。
（ｋ＋1|k）に対し、カルマンフィルタ法は，下式を
発生する： ｋ（ｋ＋1|k）＝Ａ（k|k−１） ＋Ｇ（ｋ）［ｙ（ｋ）−Cx（k|k−１）］ （５） 又は （ｋ＋1|k）＝Ａ（k|k） （5b） そこで、Ｇ（ｋ）は予測利得行列である。（Ｇ）ｋは
次式により与えられる。
散測定置ｙ（ｋ）、すなわち、記録パルスとして観測さ
れ、これは以下の式によりシステムに関連する： ｙ（ｋ）＝Cx（ｋ）＋ｖ（ｋ） （３） ここで、Ｃ＝（1,0,0）及びｖ（Ｋ）は平均ゼロ及び
共分散行列を有する測定置を妨害する白色雑音処理であ
る： Ｅ［ｖ（ｋ）.v（１）］＝Ｒ（ｋ）δk,1 （４） 更に、ｗ（ｋ）及びｖ（ｋ）は非相関である。解決さ
るべき問題点は、動作リズムｘ（ｋ＋１）が、観測され
た記録パルスｙ（ｏ）…ｙ（ｋ）から計算されなければ
ならないことからなる。この予測子は（ｋ＋1|k）と
して示される。ここで （ｋ＋1|k）＝Ｅ［ｘ（ｋ＋1|y（ｋ）］は所与のｙ
（ｋ）に対するｘ（ｋ＋１）の条件平均値である。ｘ
（ｋ）及びｙ（ｋ）が共に正規分布を有するので、動作
リズムの統計値はこの分布の条件平均値及び分散により
完全に説明される。共通条件確率分布はカルマンフィル
タ法により反復的に与えられる。カルマンフィルタ法は
予測された次の記録パルス Ｅ［（tn+1−tn）２］の平均二乗誤差を最小にする。
（ｋ＋1|k）に対し、カルマンフィルタ法は，下式を
発生する： ｋ（ｋ＋1|k）＝Ａ（k|k−１） ＋Ｇ（ｋ）［ｙ（ｋ）−Cx（k|k−１）］ （５） 又は （ｋ＋1|k）＝Ａ（k|k） （5b） そこで、Ｇ（ｋ）は予測利得行列である。（Ｇ）ｋは
次式により与えられる。
Ｇ（ｋ）＝AP（k|k−１）CT［CP（k|k−１）CT ＋Ｒ（ｋ）］-1 （６） カルマン利得行列Ｋ（ｋ）は AK（ｋ）＝Ｇ（ｋ） （6b） として定義され、 Ｐ（Ｋ＋1|1k）は下記の予測平均二乗誤差行列であ
る： Ｐ（ｋ＋1|k）＝（Ａ−Ｇ（ｋ）Ｃ）Ｐ（k|k −１）AT＋Ｑ（ｋ） （７） 平均二乗誤差行列Ｐ（k|1k）は次式である： Ｐ（k|k）＝Ｅ［（（ｋ）−ｘ（ｋ）） （（ｋ）−ｘ（ｋ））Ｔ］ （８） 上記適用に対し、 とすると初期値は次の通りである： １（３）＝t3 ２（３）＝t3−t2,及び ３（３）＝t3−2t2＋t1 Ｐ（3|3）に対し、初期値は次の通りである： 式（５）から、ｘの予測値に対し下式が得られる。
る： Ｐ（ｋ＋1|k）＝（Ａ−Ｇ（ｋ）Ｃ）Ｐ（k|k −１）AT＋Ｑ（ｋ） （７） 平均二乗誤差行列Ｐ（k|1k）は次式である： Ｐ（k|k）＝Ｅ［（（ｋ）−ｘ（ｋ）） （（ｋ）−ｘ（ｋ））Ｔ］ （８） 上記適用に対し、 とすると初期値は次の通りである： １（３）＝t3 ２（３）＝t3−t2,及び ３（３）＝t3−2t2＋t1 Ｐ（3|3）に対し、初期値は次の通りである： 式（５）から、ｘの予測値に対し下式が得られる。
（4|3）＝Ax（3|3）＝３（t3−t2）＋t1 式（７）からＰの予測値に対し下式が得られる。
Ｐ（4|3）＝AP（3|3）AT＋Ｑ（ｋ）のようになり、修
正が下記の次の時点で求められている。
正が下記の次の時点で求められている。
（4|4）＝（4|3）＋ｋ（ｋ）（ｙ（４） −Cx（4|3） （９） ｋを用いて、カルマン利得行列は下式になる。
ｋ（ｋ）＝Ｐ（4|3）CT（CP（4|3）CT＋0.5）-1 （1
0） Ｐ（4|4）は下式により与えられる： Ｐ（4|4）＝Ｐ（4|3）−KCP（4|3） （11） （９），（10）及び（11）式の計算の後、（5|4）
が（５）式等により計算されうる。
0） Ｐ（4|4）は下式により与えられる： Ｐ（4|4）＝Ｐ（4|3）−KCP（4|3） （11） （９），（10）及び（11）式の計算の後、（5|4）
が（５）式等により計算されうる。
第１図はＸ線検査装置を実施するに適した装置のブロッ
ク図、 第2a〜2e図は一連の記録パルスを関連した同期パルスと
共に示す図、 第３図はカルマンフィルタ法による次の記録パルスに対
する計算方法のフローチャート図である。 １……Ｘ線源、３……対象、４……細部、5,6……Ｘ線
ビーム、７……Ｘ線検出器、９……テレビジョンカメラ
管、11……アナログディジタル変換器、12,13……メモ
リー、15……減算装置、16……ディジタルアナログ変換
器、17……モニター、19……バックグラウンドメモリ
ー、20……記録装置、22……同期装置、23……テレビジ
ョンモニター、t1,t2,t3,t4,t5,t6……記録パルス、x1,
x2,x3,x4,x5……時点。
ク図、 第2a〜2e図は一連の記録パルスを関連した同期パルスと
共に示す図、 第３図はカルマンフィルタ法による次の記録パルスに対
する計算方法のフローチャート図である。 １……Ｘ線源、３……対象、４……細部、5,6……Ｘ線
ビーム、７……Ｘ線検出器、９……テレビジョンカメラ
管、11……アナログディジタル変換器、12,13……メモ
リー、15……減算装置、16……ディジタルアナログ変換
器、17……モニター、19……バックグラウンドメモリ
ー、20……記録装置、22……同期装置、23……テレビジ
ョンモニター、t1,t2,t3,t4,t5,t6……記録パルス、x1,
x2,x3,x4,x5……時点。
フロントページの続き (56)参考文献 特開 昭61−83987（ＪＰ，Ａ) 実開 昭49−112171（ＪＰ，Ｕ)
Claims (4)
- 【請求項１】Ｘ線源により出射されたＸ線ビームによ
り、対象を照射することによって、周期的に動く対象の
時系列Ｘ線画像を作成するためのＸ線検査装置におい
て、 前記Ｘ線検査装置は、記録装置と同期装置を有し、 前記記録装置は、Ｘ線源を作動させるための同期装置
に、上記周期的に動く対象の動きに応じた記録パルスを
印加し、 前記同期装置は、印加された記録パルスからこれから発
生する記録パルスの時点を計算し、 発生した記録パルスの現実の時点と前記計算された時点
を比較して、周期的動作における変化を検出し、 発生した記録パルスを統計的変動と認識することによ
り、同期を抑圧して、Ｘ線源の同期を周期的動作の変化
に適合させる ことを特徴とするＸ線検査装置。 - 【請求項２】Ｘ線源により出射されたＸ線ビームによ
り、対象を照射することによって、周期的に動く対象の
時系列Ｘ線画像を作成するためのＸ線検査装置におい
て、 前記Ｘ線検査装置は、記録装置と同期装置を有し、 前記記録装置は、Ｘ線源を作動させるための同期装置
に、上記周期的に動く対象の動きに応じた記録パルスを
印加し、 前記同期装置は、これから発生する記録パルスを計算す
る反復カルマンフィルタ法を含み、印加された記録パル
スからこれから発生する記録パルスの時点を計算し、発
生した記録パルスの現実の時点と計算された時点を比較
して周期的動作における変化を検出し、Ｘ線源の同期を
周期的動作の変化に適合させる ことを特徴とするＸ線検査装置。 - 【請求項３】前記同期装置は、３つの順次記録パルス
tn、tn-1及びtn-2から計算され、tnに等しい第１の要
素、tn−tn-1に等しい第２の要素、tn−2tn-1＋tn-2に
等しい第３の要素を計算するものであって、 第１要素、第２要素及び第３要素は、カルマンフィルタ
によって計算されるこれから発生する記録パルスの時点
を三次元ベクトルによる時間の関数として表したときの
各要素であることを特徴とする請求項２記載のＸ線検査
装置。 - 【請求項４】Ｘ線源により出射されたＸ線ビームによ
り、対象を照射することによって、周期的に動く対象の
時系列Ｘ線画像を作成するためのＸ線露光装置における
露光同期方法において、 前記同期方法は、 Ｘ線源を作動させるための同期装置に記録パルスを印加
し、 印加された記録パルスからこれから発生する記録パルス
の時点を計算し、 発生した記録パルスの現実の時点と前記計算された時点
を比較して周期的動作における変化を検出し、 発生した記録パルスを統計的変動と認識することによ
り、同期を抑圧して、Ｘ線源の同期を周期的動作の変化
に適合させる ことを特徴とするＸ線露光装置における露光同期方法。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| NL8901322A NL8901322A (nl) | 1989-05-26 | 1989-05-26 | Roentgenbeeldsynchronisatiemethode. |
| NL8901322 | 1989-05-26 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0382000A JPH0382000A (ja) | 1991-04-08 |
| JP3105520B2 true JP3105520B2 (ja) | 2000-11-06 |
Family
ID=19854717
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP02133715A Expired - Fee Related JP3105520B2 (ja) | 1989-05-26 | 1990-05-23 | Ｘ線検査装置 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US5040201A (ja) |
| EP (1) | EP0399606B1 (ja) |
| JP (1) | JP3105520B2 (ja) |
| DE (1) | DE69013571T2 (ja) |
| NL (1) | NL8901322A (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2151192A1 (en) | 2008-08-05 | 2010-02-10 | Canon Kabushiki Kaisha | X-ray imaging apparatus |
| JP6046853B1 (ja) * | 2016-07-01 | 2016-12-21 | 範子 遣水 | 食材調理用溝切りカッター |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| MX9702434A (es) | 1991-03-07 | 1998-05-31 | Masimo Corp | Aparato de procesamiento de señales. |
| US5490505A (en) | 1991-03-07 | 1996-02-13 | Masimo Corporation | Signal processing apparatus |
| US5271055A (en) * | 1992-08-19 | 1993-12-14 | General Electric Company | Methods for reducing motion induced artifacts in a projection imaging system |
| US7376453B1 (en) | 1993-10-06 | 2008-05-20 | Masimo Corporation | Signal processing apparatus |
| EP1905352B1 (en) | 1994-10-07 | 2014-07-16 | Masimo Corporation | Signal processing method |
| US8019400B2 (en) | 1994-10-07 | 2011-09-13 | Masimo Corporation | Signal processing apparatus |
| WO1998032376A1 (en) * | 1997-01-29 | 1998-07-30 | Picker Medical Systems, Ltd. | Predictive bolus tracking |
| DE102008032296B4 (de) | 2008-07-09 | 2010-06-02 | Siemens Aktiengesellschaft | Röntgengerät |
| JP6848407B2 (ja) * | 2016-12-09 | 2021-03-24 | コニカミノルタ株式会社 | 放射線画像撮影システム |
| CN111513738B (zh) * | 2020-04-10 | 2023-08-01 | 北京东软医疗设备有限公司 | 血管造影方法、装置、设备及系统 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3871360A (en) * | 1973-07-30 | 1975-03-18 | Brattle Instr Corp | Timing biological imaging, measuring, and therapeutic timing systems |
| US4245647A (en) * | 1978-12-06 | 1981-01-20 | General Electric Company | Arrhythmia rejection circuit for gated cardiac image display systems |
| EP0193712B1 (de) * | 1985-02-04 | 1990-04-25 | Siemens Aktiengesellschaft | Röntgendiagnostikeinrichtung für Subtraktionsangiographie |
-
1989
- 1989-05-26 NL NL8901322A patent/NL8901322A/nl not_active Application Discontinuation
-
1990
- 1990-05-21 DE DE69013571T patent/DE69013571T2/de not_active Expired - Fee Related
- 1990-05-21 EP EP90201274A patent/EP0399606B1/en not_active Expired - Lifetime
- 1990-05-22 US US07/527,998 patent/US5040201A/en not_active Expired - Fee Related
- 1990-05-23 JP JP02133715A patent/JP3105520B2/ja not_active Expired - Fee Related
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2151192A1 (en) | 2008-08-05 | 2010-02-10 | Canon Kabushiki Kaisha | X-ray imaging apparatus |
| US7856085B2 (en) | 2008-08-05 | 2010-12-21 | Canon Kabushiki Kaisha | X-ray imaging apparatus |
| JP6046853B1 (ja) * | 2016-07-01 | 2016-12-21 | 範子 遣水 | 食材調理用溝切りカッター |
Also Published As
| Publication number | Publication date |
|---|---|
| EP0399606B1 (en) | 1994-10-26 |
| JPH0382000A (ja) | 1991-04-08 |
| EP0399606A1 (en) | 1990-11-28 |
| DE69013571D1 (de) | 1994-12-01 |
| DE69013571T2 (de) | 1995-05-11 |
| US5040201A (en) | 1991-08-13 |
| NL8901322A (nl) | 1990-12-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US4903705A (en) | Digital radiography apparatus | |
| US5809105A (en) | Noise filter for digital x-ray imaging system | |
| US6940943B2 (en) | Continuous scan tomosynthesis system and method | |
| JP3105520B2 (ja) | Ｘ線検査装置 | |
| Reiber et al. | Accuracy and precision of quantitative digital coronary arteriography: Observer‐, short‐, and medium‐term variabilities | |
| US5293574A (en) | Digital x-ray imaging system with automatic tracking | |
| EP0362821A1 (en) | Diagnostic X-ray apparatus | |
| JP3172180B2 (ja) | Ｘ線透視法におけるダイナミックパルス制御 | |
| JPH02211131A (ja) | 周期的動きに起因する偽像を低減する装置 | |
| JPH031869B2 (ja) | ||
| JPH04212333A (ja) | Ｘ線検査装置およびその作動方法 | |
| KR860001794B1 (ko) | X선 진단 장치 | |
| US4716904A (en) | Heart output and circulatory impedance measuring method and apparatus | |
| JPH05192319A (ja) | Ｘ線診断装置 | |
| US8103078B2 (en) | Method for determining ECG-triggered recording times for imaging to support interventional and diagnostic cardiac procedures | |
| JP2004188196A (ja) | 心臓サイクルの検出の方法及び装置 | |
| US6332014B1 (en) | Method of improving the quality of a fluoroscopic image | |
| US4720843A (en) | Method for separating moving structures from a fixed background in a sequence of X-ray projection images and equipment for implementing this method | |
| US4941169A (en) | X-ray photographic equipment | |
| US4618976A (en) | X-ray diagnostic installation | |
| JP2004523329A (ja) | 患者の心血管系の画像の質を向上させるためのｘ線装置、コンピュータプログラム及び方法 | |
| EP0097355A1 (en) | Method and apparatus for performing digital intravenous subtraction angiography | |
| JPS59231985A (ja) | Ｘ線診断装置 | |
| EP0102592B2 (en) | X-ray image producing system | |
| US6404850B1 (en) | Cardioangiography apparatus |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |