JP2005521501A - 画像シーケンスから画像を決定する方法 - Google Patents
画像シーケンスから画像を決定する方法 Download PDFInfo
- Publication number
- JP2005521501A JP2005521501A JP2003581122A JP2003581122A JP2005521501A JP 2005521501 A JP2005521501 A JP 2005521501A JP 2003581122 A JP2003581122 A JP 2003581122A JP 2003581122 A JP2003581122 A JP 2003581122A JP 2005521501 A JP2005521501 A JP 2005521501A
- Authority
- JP
- Japan
- Prior art keywords
- image
- motion
- motion signal
- sequence
- image sequence
- 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.)
- Pending
Links
- 230000033001 locomotion Effects 0.000 claims abstract description 115
- 238000000034 method Methods 0.000 claims abstract description 55
- 238000012545 processing Methods 0.000 claims abstract description 20
- 238000004590 computer program Methods 0.000 claims abstract description 10
- 239000002872 contrast media Substances 0.000 claims description 26
- 210000004204 blood vessel Anatomy 0.000 claims description 10
- 238000002604 ultrasonography Methods 0.000 claims description 4
- 241001465754 Metazoa Species 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 2
- 230000000747 cardiac effect Effects 0.000 description 18
- 238000002565 electrocardiography Methods 0.000 description 10
- 230000001186 cumulative effect Effects 0.000 description 7
- 238000003384 imaging method Methods 0.000 description 6
- 230000002792 vascular Effects 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 4
- 229940039231 contrast media Drugs 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000029058 respiratory gaseous exchange Effects 0.000 description 3
- 238000002583 angiography Methods 0.000 description 2
- 238000010009 beating Methods 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 208000019622 heart disease Diseases 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000002107 myocardial effect Effects 0.000 description 1
- 210000004165 myocardium Anatomy 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
- 210000000689 upper leg Anatomy 0.000 description 1
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/38—Registration of image sequences
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30048—Heart; Cardiac
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Optics & Photonics (AREA)
- Heart & Thoracic Surgery (AREA)
- Theoretical Computer Science (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- General Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Physiology (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Analysis (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本発明は、動いている対象物の画像シーケンスであって、動きの状態のシーケンスとして対象物の動きを表す画像シーケンスからの基準画像に対する対応する画像を決定する方法に関する。この目的のために、関連する動きの進行を表す2つの動き信号が類似性について調べられる。得られた類似性関数を用いて、基準画像において表される少なくとも略対象物の動きの状態を表す、画像シーケンスにおける画像を決定することができる。更に、本発明は、この方法を実行するために適切なシステム、データ処理ユニットがこの方法を実行することを可能にするコンピュータプログラム及びコンピュータプログラムプロダクトに関する。
Description
本発明は、動いている対象物の画像シーケンスであって、動きの状態のシーケンスとして対象物の動きを表す画像シーケンスからの基準画像に対する対応する画像を決定する方法に関する。本発明は又、データ処理ユニットがその方法を実行することを可能にするコンピュータプログラムプロダクト及びコンピュータプログラムに対して同様にその方法が実行されることができる、システム及び検査装置に関する。
上記方法は、例えば、動きの状態の画像が画像シーケンスにおいて決定され、前記の動きの状態が又、対象物の第2の類似する状態の間に捕捉された基準画像において表されるどのような場合においても用いられる。対象物が動いている間に、各々の動きに対する動きの状態に合わせたシーケンスを表す信号が、逐次、決定される。この種の装置は、例えば、医療分野で周知である。米国特許第4,729,379号明細書において、鼓動している心臓の2つのX線画像シーケンスが捕捉された心臓検査のためのX線検査承知について開示されている。2つのシーケンスの1つの捕捉の間に、心臓の血管にコントラスト媒体が注入される。
次の段階の間に、2つのシーケンスは互いに即ち1つの画像は他の画像の後に減算され、それ故、コントラスト媒体で満たされた血管は、結果としての差分シーケンスにおけるバックグラウンドの最小量を用いて再生される。この目的のために、2つの画像シーケンスは、同じ動きの状態を表す画像が各々の画像シーケンスから減算されるような方法で、2つの画像シーケンスは互いに対して位置合わせされる必要がある。これは、心電図によりそれぞれの心電(ECG)信号の捕捉により、即ち、2つの画像シーケンスの捕捉と同時に達成される。両方のECGにおいて、2つの連続するR歪みが決定される度に、2つのECGは互いに対して位置合わせされる。2つのECGにおける2つのR歪みの間の時間経過が異なる場合、この時間差は線形補完により補償され、それ故、2つの対応する画像シーケンスのR歪みの間で捕捉される画像は互いに関連付けられることができる。
1つの基準例のみが心臓周期に対して用いられるため、画像シーケンスは1つの瞬間のみに互いに位置合わせされ、それ故、一方で、全体的な心臓周期の期間に関し、他方で、心臓の動きの個々のセグメントの圧縮又は拡張に関して、2つのECG信号の間の差が考慮される。このことは、差分シーケンスにおける不所望で厄介なアーチファクトを生じさせる。
従って、本発明の目的は、上記のようなシステムを改善することである。
この目的は、本発明に従って、第1動き信号及び第2動き信号により動いている対象物の画像シーケンスからの基準画像についての対応する画像を決定する方法であって、
− 第1動き信号及び第2動き信号は、対象物の第1動き及び第2動きの動きの状態に調子を合わせたそれぞれの変化を表し、
− 画像シーケンスは、動きの状態の画像のシーケンスとして対象物の第1動きを表し、
− 基準画像は第2対象物の動きから動きの状態を表し、対象物の第2動きの間の基準の瞬間において捕捉される、
方法であり、
a.第1動き信号及び第2動き信号の類似性比較としての類似性関数を決定する段階、
b.その類似性関数により第1動き信号における対応瞬間を決定する段階であって、その対応瞬間は第2動き信号から基準画像の捕捉瞬間に対応する、段階、並びに
c.第1動き信号を用いて、捕捉瞬間が少なくとも対応瞬間に略対応する画像シーケンスの画像を決定する段階、
を有する方法により達成される。
− 第1動き信号及び第2動き信号は、対象物の第1動き及び第2動きの動きの状態に調子を合わせたそれぞれの変化を表し、
− 画像シーケンスは、動きの状態の画像のシーケンスとして対象物の第1動きを表し、
− 基準画像は第2対象物の動きから動きの状態を表し、対象物の第2動きの間の基準の瞬間において捕捉される、
方法であり、
a.第1動き信号及び第2動き信号の類似性比較としての類似性関数を決定する段階、
b.その類似性関数により第1動き信号における対応瞬間を決定する段階であって、その対応瞬間は第2動き信号から基準画像の捕捉瞬間に対応する、段階、並びに
c.第1動き信号を用いて、捕捉瞬間が少なくとも対応瞬間に略対応する画像シーケンスの画像を決定する段階、
を有する方法により達成される。
上記の方法は、基準画像についての対応する画像を決定する役割を果たし、対応する画像は、基準画像において表される、少なくとも動いている対象物の略その動きの状態を表す。対象物の第1動きから、各々の画像が対象物の動きの状態を表す画像シーケンスをここで捕捉する。次いで、画像の連続が、対象物の動きの動き画像シーケンスを表す。対象物が第2動きを実施している間に、基準画像は、対象物の第2動きの間に生じる動きの状態について捕捉される。動きについての動きの状態に調子をあわせて変化を表し又は特徴付ける動き信号が、第2動き及び第1動きから利用可能である。この種の信号は、例えば、関連する動きが生じる間に捕捉されることができるECGである。他の動き信号は、患者の呼吸の動きの間に、呼吸センサにより生成される信号である。
この方法の第1段階の間に、2つの動き信号は類似性を調べられる。これは、対象物が2つの瞬間において少なくとも略同じ動きの状態を仮定されるような方法において、他の動き信号の対応する瞬間を1つの動き信号のいずれの瞬間に対応付けるために用いられることができる類似性関数を与える。類似性関数を用いると、2つの動きにおいて、対象物がそれぞれの動きの間に略同じ動きの状態を仮定する2つの瞬間を決定することができる。動きの動き信号が互いに対して非線型に歪んでいるような領域に対して動きが異なるとき、既知の方法と反対に、本発明に従った方法は、評価について適切である結果を生成する。
この方法の第2段階の間に、第2動き信号の基準の瞬間に対応する第1動き信号における瞬間が決定される。この方法の第3段階において、捕捉瞬間が対応する瞬間に略対応する画像シーケンスの画像は、対応する画像として決定される。このようにして選択された対応する画像は、基準画像において表される、動いている対象の少なくとも略動きの状態を表す。画像シーケンスからの基準画像及び対応する画像が互いから減算されるとき、減算画像は、異なる動きの状態の2つの画像が互いから減算される事実によることが可能である非常に少しおんアーチファクトのみを示す。
請求項2に従って、類似性の比較を、既知の“動的タイムワーピング”方法により実行することができる。この方法は、類似性の比較の非常に速く、効率的な実行を可能にする。
対象物の動きが原理的に認識される場合、それ故、付加情報が、画像シーケンスの画像の捕捉の間に対象物により仮定された動きの状態に関して提供されることができる場合、人工的な中間的画像が、補間によりこれらの動きの状態に対して生成されることができる。例えば、呼吸の間に一部の組織により実行される動きは、動きモデルにより十分正確に表されることができる。請求項3に従って、基準画像の捕捉の間に、対象物により仮定された動きの状態を同様に表す画像を生成することが可能である。これは、画像シーケンスの画像の捕捉の実際の瞬間と対応する瞬間との間の差が非常に大きいとき、例えば、減算画像において非常に多くのアーチファクトが生成され、特に有利である。中間画像の補間は又、対象物の動きの間に画像シーケンスに対してほんの僅かな画像が捕捉されるとき、有利に用いられることができる。
請求項4に従って、本発明の方法は、人間又は動物の心臓の画像又は画像シーケンスが生成され、更に、心臓の動きのECGが利用可能であるシステムにおいて、特に有利に用いられることができる。本発明に従った方法は、病気又は所定の物理的条件のためにECGが互いに対して非線形歪みを示す心臓周期を有する患者の場合に特に、心臓の動きの同じ状態を表す画像の信頼できる決定を提供する。
請求項5に従って、本発明に従った方法は、減算造影法が実行されるシステムにおいて用いられる。心臓の動きの状態の画像の捕捉に対して適切である画像化方法は、請求項6に一致して、X線システムとすること、及び、請求項7に一致して、超音波システムとすることが可能である。しかしながら、磁気共鳴断層撮影装置又はX線CT(computed tomography)装置のようなスライス画像又はボリューム画像を生成する画像化システムは又、非常に近い将来、心臓の動きの状態を記録することができるようになるであろう。この方法は、それ故、結果的に用いられることができる。
従って、請求項8に従って、医療分野において現代の画像化システムで用いられるこの種のデータ処理ユニットを有するシステムにおける方法を用いることは特に有利である。X線検査システムは、請求項9に従ったシステムを用いて提供されることが可能である。データ処理ユニットがプログラム可能であるように構成される場合、請求項10に記載されているようなコンピュータプログラム又はコンピュータプログラムプロダクトは、データ処理ユニットが本発明に従った方法を実行することを可能にする。
本発明の実施形態について、図を参照して、以下、詳細に説明する。
図1は、時間tに対する心臓周期の心電図(ECG)E1を表す図である。心電図は、人間又は動物の心臓の活動電流を記録する方法である。活性化された心筋の位置は、非活性位置に対して電荷を有する。この電荷は、心筋の残りの部分を通って所定の速度で伝播する。そのような電流は、図1に従って再生されることができ、身体に取り付けられる適切な電極により既知の方法で測定される。人間心臓周期の典型的な持続時間は約1秒である。
心臓周期の間に、心臓のそれぞれの画像I1乃至I14が、複数の瞬間t1乃至t14において捕捉される。各々の画像I1乃至I14は、非常に複雑な心臓の動きの瞬間的画像として心臓の動きの状態を表す。この種の画像は、例えば、X線画像化又は超音波画像化等の既知の画像化方法により捕捉されることができる。現代のX線システムは、最大1秒当たり30乃至60画像の捕捉が可能であり、それ故、心臓周期の画像シーケンスの画像はそれぞれ、30乃至60の心臓の動きの異なる状態を表す。しかしながら、そのような多くの画像は図に表すことは困難であり、それ故、ここでは、少ない数を用いている。
一部の心臓の検査に対しては、心臓の血管ツリーのみを画像化することが有用である。この目的のために、心臓の動きの状態について第1画像及び第2画像が捕捉される。第2画像の捕捉の間に、例えば、X線を吸収するコントラスト媒体が、例えば、カテーテルにより心臓の血管に導入され、それ故、血管はX線システムにおいて非常にうまく強調される。第1画像の捕捉は、コントラスト媒体なくして得られることはない。次の段階の間に、両方の画像は、例えば、互いから減算され、例えば、1つの画素は他の画素から減算され、それ故、理想的な場合には、コントラスト媒体で満たされた血管ツリーのみが尚、視認される。この方法は又、減算血管造影と呼ばれ、デジタル画像の場合、デジタル減算血管造影(DSA)と呼ばれる。図1は、心血管がコントラスト媒体で満たされる第2画像シーケンスの生成及び対応する心臓周期のECG E1´を図に表したものである。画像I´1乃至I´14は瞬間t´1乃至t´14において捕捉される。
減算により取り除かれる画像要素又は画素が2つの画像において実質的に同じであるとき、DSAにより生成される画像はアーチファクトを殆ど含んでいない。更に、2つの画像において強調されるべき対象物の形状は又、実質的に同じである必要がある。従って、特に、心臓の動きのように複雑な動きの場合、心臓の動きの同じ状態を表す2つの画像を互いに減算することは重要である。例えば、画像I8において表される動きの状態がDSA画像として表されるべきである場合、画像I8と同じ心臓の動きの状態を実質的に表す画像I8を第2画像シーケンスから見つけることが必要である。これは、破線で示されている。今まで用いられてきたシステムにおけるように、特に、所謂R偏向(I8の捕捉の瞬間の少し前の著しい偏向)のみが参照例として取り上げられ、この偏向がECG E1(ECG E1´の破線部分)におけるよりECG E´1において実質的に広い場合、既知の方法はR偏向の減少された部分に尚も位置されている画像I´8を生成し、それ故、それは、実際に調査されたI´8とは全く異なる動きの状態を表す。特に、心疾患を患っている患者の場合、2つの異なる心臓周期のECGの曲線が全く同様であることを単純に仮定することはできない。一般的に言って、1つのECGは、他のECGに対する非線型歪みを表し、それ故、1つのECGの捕捉の間の所定の瞬間において、心臓は、他のECGの捕捉の間の動きの状態とは異なる動きの状態にある。このことは、画像I8に対応する画像I´8の決定のための2つのECGの比較を非常に困難にする。
図2は、2つの時間的に非線形な歪み信号が類似性に関して調べられることを可能にする実行可能な方法の一般的表現である。水平方向と垂直方向にプロットされた(一次元)関数は、2つの略同じ関数に対応している。2つの関数の異なる時間的構成は、それらの時間スケールの相対的歪みにおいて明らかとなった。動きの状態のような対応するイベントの瞬間の仮定された関連性は、スケールのプロダクトグリッド(product grid)における経路として印が付けられている。
そのようなパターン比較アルゴリズムを使用することは、発話認識の分野において利用される単語認識システムにより認識される(例えば、文献、E.G.Schukat−Talamazzini ”Automatiche Spracherkennung”,ISBN 3−528−05492−1,Vieweg Verlag,1995,chapter5.1において)。図3との一致性において、基準パターンYl(又、プロトタイプ又はテンプレートと呼ばれる)により表される語彙の入力毎に語彙Y=Y1...YTが与えられ、各々の基準パターンYlは、特徴ベクトルシーケンスYl=yl1...ylsの形式で与えられる。語彙入力は、システム寸法記入フェーズの間に捕捉された単語の選択的認識を生成する。更に、短期間に解析された単語は、特徴ベクトルシーケンスX=x1...xTの形式で与えられる。個々の単語を認識するタスクは、語彙のどの単語が実際に発話された単語に最もよく対応するかを調べることにより、語彙から発話された可能性の大きい単語の同一性を決定することにある。
異なる期間S≠Tの基準シーケンスY=y1...ySと入力シーケンスX=x1...xTとの間の距離D(X,y)は、ベクトルシーケンスの間の適切なタイムワーピングに沿った局所距離の合計として決定される。局所距離関数d(・,・)は、例えば、ユークリッドメトリックにより実現される。この目的のための適切な歪み関数は、全体的な長さに亘ってYにおいてXをマッピングし、tスケール及びsスケールにおいて所定の単調性と恒常性の特性を満足し、そして最小の全体的距離を与える。この非常に複雑な離散最適化タスク(組み合わされた実現可能経路の数は、前記制限が課せられるときでさえ、プロトタイプの長さの関数として指数関数的に増加する)は、最適化の原理を満足し、それ故、一般に周知の“動的プログラミング”方法により解かれることができる。ベクトルシーケンスX、Yの初期セグメント間の累積距離Dy=D(x1...xi,y1...yj)は次の漸化式を条件とする。
この方法は、2つのECG曲線が類似性のために比較されることを可能にする。これを達成するために、ECG曲線はデジタルの方式で提供され、即ち、各々の曲線における多数の特徴的ポイント(例えば、300)により表される。これは、ECG信号をサンプリングすることによる既知の方法において実現されることができる。それ故、ECG曲線はスカラのシーケンスをなす。局所距離関数は、それ故、次の代替の基準により表される。
dij=d(xi,yi)=|xi−yj|,
dij=d(xi,yi)=(xi−yj)2,または
一般に基準として‖xi,yj‖を用いて、dij=d(xi,yi)=‖xi,yj‖である。
dij=d(xi,yi)=|xi−yj|,
dij=d(xi,yi)=(xi−yj)2,または
一般に基準として‖xi,yj‖を用いて、dij=d(xi,yi)=‖xi,yj‖である。
式‖xi,yj‖は、一般に、数学で周知であり、条件、‖x,y‖+‖y,z‖≧‖x,z‖(三角不等式)、‖x,x‖=0、‖x,y‖≧0及び‖x,y‖=‖y,x‖(対称)である。
上記のように、累積距離は前の局所距離の合計として構成され、最も好ましい歪み経路は、最小条件を用いて漸化的に決定される。一旦、歪み経路が規定されると、他のECG曲線の特性ポイントを、一のECG曲線の各々の特性ポイントに対して示すことができる。
最適な経路を決定する間に、漸化式は、漸化段階においてどのボックスに到達することができるかを規定する。上記の式に従って、3つの異なる累積距離から最小値が決定され、それら3つの距離は1つのボックスから直接隣接するボックスまで測定される。それに代えて、対応する漸化式からもたらされる次の漸化マトリクスを用いることができる。
図5及び6は、2つのECGに対する幾つかの局所距離についてのそれぞれの図を示している。簡単のために、信号値は、0と10との間の絶対値を有する、対応する整数値に四捨五入(量子化)される。時間軸tの各々の特性ポイントに対してボックスがあり、時間毎の第1特性ポイントは、x1乃至x5及びy1乃至y5により示されている。局所距離は、dij=d(xi,yi)=|xi−yj|としての時間毎に計算される。最適経路は、次の漸化式と累積距離(図示せず)により決定される。
図7は、対応する瞬間が規定されるときに、一致する、対応する画像を選択するための可能性を示している。ECG E5の捕捉の間に、心臓の血管はコントラスト媒体で満たされ、画像I51は瞬間t51において捕捉される。ECG E6及びE7および対応する画像シーケンスの捕捉の間、コントラスト媒体は使用されない。ECG E6との類似性の比較は、画像I51と同じ心臓の動きの状態を表す、対応する画像I63が捕捉された、対応する瞬間t63を与える。画像I63を減算アルゴリズムのために用いることができる。ECG E7との類似性の比較は対応する瞬間t75を与える。しかしながら、瞬間t75において破線の矢印により示されるように、その瞬間に画像は捕捉されない。対応する瞬間t75に先立ち、画像I72は捕捉され、画像I73は瞬間t75の後に捕捉される。
減算造影図に適する画像を得る第1の可能性は、捕捉の瞬間が対応する瞬間t75のできるだけ近くにある画像を選択することである。換言すれば、その画像は、画像の捕捉の瞬間と対応する瞬間t75との間の時間差が最も小さくなるように選択される。これは、画像I73である。この可能性は、コントラスト媒体を用いることなく捕捉されたECG(図示せず)を有する更なる画像が存在するときに特に、用いられる。これらの更なるECGの各々に対して、この方法は、関連する対応する瞬間と最近接の画像との間の時間差を決定することを可能にする。この時間差が最小である画像が、減算造影図に対して選択される。
減算造影図に適する画像を得るための他の可能性は、対応する瞬間の前後の各々の時間に捕捉された画像から画像を補間することである。例えば、図I72とI73に示された動きの状態の間の心臓の動きが、対応する瞬間t75における心臓の動きの状態を表す人工的画像を生成するように、補間される。他の場合、心臓の全体的動きが動きモデルで表され、それにより、画像I72の捕捉の瞬間とI73の捕捉の瞬間との間の心臓の殆ど正確な動きが補間され、対応する瞬間t75における動きの状態の画像が決定されることができる。
本発明に従った方法を使用することは、1つの基準画像に限定されない。複数の基準画像を有する画像シーケンスの場合、対応する画像は、上記の方法により各々の個々の画像に対して決定又は生成されることができる。
図8は、医療X線検査システムを図に表したものである。このシステムは、X線に対して透明であるテーブル43に配置された患者41が存在するこの場合、X線は検査される対象を横断する方式でX線42を放射するように配列されたX線源40を有する。従って、X線がテーブル43の下に配置されたX線画像検出器44により検出することができる。X線画像検出器44は、X線に対して感度を有するセンサ素子のアレイを有する。患者41を横断するX線42についての器官に特定の減衰のために、画像はX線画像検出器44において生成され、前記画像のデータはデータ処理ユニット46に供給される。データ処理ユニット46は又、ここで説明していないシステム制御タスクを実行することができ、最適画像がビューアに対して生成される方式で入力画像データ45を処理する。このようにして処理される画像データ47は、視覚化ユニット48であって、例えば、データがビューアに提供されることができるモニタに提供される。データ処理ユニットは、プログラマブルであるように任意に構成されることができる。その場合、データ処理ユニット46に結合され、データ処理ユニット46に適用するようにコンピュータプログラムプロダクトからコンピュータプログラムを読み出すことができるデータ読み出し装置52が利用可能である。コンピュータプログラムは、特に、データ処理ユニットが本発明に従った方法を実行することを可能にする。
ECG装置51に接続される電極は患者に取り付けられる。図8は、身体に取り付けられた電極を有する1つの接続リード線50のみを単なる例として示している。実際には、複数の電極が、基準電位(しばしば、接地)に対する電極(図示せず)と同様な既知の方式で使用される。データ処理ユニット46に接続されたECG装置は、データ処理ユニット46の制御下でX線画像捕捉の間に患者41のECGを生成し、本発明に従った方法を実行するようにデータ処理ユニット46にECGを供給する。
又、心臓検査の間に患者の腿の付け根の領域における血管に代表的に導かれたカテーテルを示している。このとき、医師は、身体内のナビゲーションのための支援としての役割を果たす低照射量を用いて、X線画像シーケンスを捕捉しながら、心臓までカテーテルの先端を進める。一旦、心臓に達すると、カテーテルにより、コントラスト媒体が心臓の血管に注入される。カテーテルの先端からコントラスト媒体が現れる少し前に、X線画像検出器44とX線源40は、次の捕捉のために動作の高照射量モードに切り換えられ、それ故、詳細な高照射量画像が、コントラスト媒体で満たされた心臓の血管ツリーとして生成される。
その後、X線源は動作を停止され、又は低照射量画像が再び生成され、医師は更なる動作のためのナビゲーションの支援を必要とする。コントラスト媒体により満たされた血管ツリーの高照射量画像がデータ処理ユニット46において記憶される。本発明に従った方法とECG装置41により生成されるECGを用いて、これらの高照射量画像が、低照射量画像又は更なる高照射量画像に上記の方法で重ね合わされることができ、視覚化ユニット48により医師に提供されるように、コントラスト媒体なしで血管ツリーを表すことができる。
本発明に従った方法は種々の心臓検査のために用いられることができる。
1.鼓動している心臓の複数の心臓周期がコントラスト媒体の注入を用いることなく、捕捉される。次いで、幾つかの心臓周期又はただ1つの心臓周期がコントラスト媒体を用いてが捕捉され、次いで、いくつかの心臓周期がコントラスト媒体を用いることなく捕捉される。コントラスト媒体を用いないそれぞれの画像が、本発明に従った方法により、コントラスト媒体で満たされた心臓をどの画像が表すかを、心臓周期の幾つかの又は全ての画像に対して決定された後、及びその減算の後、画像が心臓サイクルの動きの種々の段階における心臓の血管ツリーのみを実質的に全く表す心臓周期の画像シーケンスが得られる。これらの画像はデータ処理ユニットに記憶され、動いている画像シーケンスとして医師に即座に提供されることができる。コントラスト媒体を伴う画像は、それ故、基準画像を有し、コントラスト媒体を伴わない画像は対応する画像を生成する。
2.検査システムのデータ処理ユニットは、上記1項において生成されるような画像シーケンスを有し、その画像は、特に、心臓周期の動きの種々の状態における心臓の血管ツリーのみを表す。心臓の近くのカテーテルを位置付けるために、医師は、低照射量で連続的に捕捉されるX線画像シーケンスの生成において身体におけるナビゲーションのための支援を利用する。これらのシーケンスにおいて、血管ツリー及び心臓は、物理的条件のための困難のみを伴って認識されることができる。ちょうど捕捉されたX線画像に対して、対応する画像は、本発明に従った方法によりX線画像に重ね合わされるように、記憶されたシーケンスから自動的に又はマニュアルで決定される。結果として、コントラスト媒体で満たされた血管ツリーは、医師が適切なナビゲーションの支援を提供されるように、瞬間のX線画像に重ね合わされる。必要に応じて、そのような重ね合わせは又、実際には捕捉されたX線画像全てに対して連続的に実行される。それ故、コントラスト媒体を伴わない画像は基準画像を生成し、コントラスト媒体を伴う画像は対応する画像を生成する。
1.鼓動している心臓の複数の心臓周期がコントラスト媒体の注入を用いることなく、捕捉される。次いで、幾つかの心臓周期又はただ1つの心臓周期がコントラスト媒体を用いてが捕捉され、次いで、いくつかの心臓周期がコントラスト媒体を用いることなく捕捉される。コントラスト媒体を用いないそれぞれの画像が、本発明に従った方法により、コントラスト媒体で満たされた心臓をどの画像が表すかを、心臓周期の幾つかの又は全ての画像に対して決定された後、及びその減算の後、画像が心臓サイクルの動きの種々の段階における心臓の血管ツリーのみを実質的に全く表す心臓周期の画像シーケンスが得られる。これらの画像はデータ処理ユニットに記憶され、動いている画像シーケンスとして医師に即座に提供されることができる。コントラスト媒体を伴う画像は、それ故、基準画像を有し、コントラスト媒体を伴わない画像は対応する画像を生成する。
2.検査システムのデータ処理ユニットは、上記1項において生成されるような画像シーケンスを有し、その画像は、特に、心臓周期の動きの種々の状態における心臓の血管ツリーのみを表す。心臓の近くのカテーテルを位置付けるために、医師は、低照射量で連続的に捕捉されるX線画像シーケンスの生成において身体におけるナビゲーションのための支援を利用する。これらのシーケンスにおいて、血管ツリー及び心臓は、物理的条件のための困難のみを伴って認識されることができる。ちょうど捕捉されたX線画像に対して、対応する画像は、本発明に従った方法によりX線画像に重ね合わされるように、記憶されたシーケンスから自動的に又はマニュアルで決定される。結果として、コントラスト媒体で満たされた血管ツリーは、医師が適切なナビゲーションの支援を提供されるように、瞬間のX線画像に重ね合わされる。必要に応じて、そのような重ね合わせは又、実際には捕捉されたX線画像全てに対して連続的に実行される。それ故、コントラスト媒体を伴わない画像は基準画像を生成し、コントラスト媒体を伴う画像は対応する画像を生成する。
本発明に従った方法はX線検査システムに限定されるものではないことに留意されたい。本発明に従った方法は、図8に示すX線検査装置における方法と同様の方法で超音波検査システムにおいて実行されることができる。この目的のために、基準画像及び対応する画像が超音波により生成され、基準画像又は対応する画像の捕捉の間に検査される組織又は血管に対して超音波反射コントラスト媒体が与えられる。
Claims (10)
- 第1動き信号及び第2動き信号により動いている対象物の画像シーケンスからの基準画像について対応する画像を決定する方法であって:
前記第1動き信号及び前記第2動き信号は、前記対象物の第1動き及び第2動きの状態に調子を合わせたそれぞれの変化を表し;
前記画像シーケンスは、動きの状態の画像シーケンスとして前記対象物の前記第1動きを表し;
前記基準画像は第2の対象物の動きから動きの状態を表し、前記対象物の前記第2動きの間の基準の瞬間に捕捉される;
方法であり、
前記第1動き信号及び前記第2動き信号の類似性比較としての類似性関数を決定する段階;
前記類似性関数により前記第1動き信号における対応瞬間を決定する段階であって、前記対応瞬間は前記第2動き信号からの前記基準画像の捕捉瞬間に対応する、段階;並びに
前記第1動き信号を用いて、捕捉瞬間が少なくとも前記対応瞬間に略対応する画像シーケンスの画像を決定する段階;
を有することを特徴とする方法。 - 請求項1に記載の方法であって、前記類似性関数は、所謂、動的タイムワーピング方法により得られる、ことを特徴とする方法。
- 請求項1に記載の方法であって、補間画像は前記対応する画像から生成され、更なる画像は前記画像シーケンスから生成され、その補間画像は少なくとも実質的に前記対応瞬間の前記対象物の動きの状態を表す、ことを特徴とする方法。
- 請求項1に記載の方法であって、前記第1動き信号及び前記第2動き信号は心電図信号を生成し、前記画像シーケンスの前記画像及び前記基準画像は人間又は動物の心臓の動きの状態を表す、ことを特徴とする方法。
- 請求項4に記載の方法であって、前記心臓の血管は、前記画像シーケンスの画像又は前記基準画像のコントラスト媒体で少なくとも一部を満たされている、ことを特徴とする方法。
- 請求項1に記載の方法であって、前記画像シーケンスはX線画像シーケンスを生成し、及び/又は、前記基準画像はX線画像を生成する、ことを特徴とする方法。
- 請求項1に記載の方法であって、前記画像シーケンス及び/又は前記基準画像は超音波画像を生成する、ことを特徴とする方法。
- 第1動き信号及び第2動き信号により画像シーケンスから基準画像のための動いている対象物の対応する画像を決定するためのデータ処理ユニットを有するシステムであって、
請求項1に記載の方法を実行するために配置されている、ことを特徴とするデータ処理ユニット。 - 心電図信号の検出のための手段及びX線画像検出器を有する検査装置であって、請求項8に記載のシステムを有する、ことを特徴とする検査装置。
- データ処理ユニットが請求項1に記載の方法を実行することができるような方法で、前記データ処理ユニットと協働するように配置されているコンピュータプログラム又はコンピュータプログラムプロダクト。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE10214763A DE10214763A1 (de) | 2002-04-03 | 2002-04-03 | Verfahren zur Bestimmung eines Bildes aus einer Bildsequenz |
PCT/IB2003/001183 WO2003083777A2 (en) | 2002-04-03 | 2003-04-01 | Method of determining an image from an image sequence |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2005521501A true JP2005521501A (ja) | 2005-07-21 |
JP2005521501A5 JP2005521501A5 (ja) | 2006-05-25 |
Family
ID=28458565
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003581122A Pending JP2005521501A (ja) | 2002-04-03 | 2003-04-01 | 画像シーケンスから画像を決定する方法 |
Country Status (7)
Country | Link |
---|---|
US (1) | US20050207538A1 (ja) |
EP (1) | EP1500047A2 (ja) |
JP (1) | JP2005521501A (ja) |
CN (1) | CN100345162C (ja) |
AU (1) | AU2003215822A1 (ja) |
DE (1) | DE10214763A1 (ja) |
WO (1) | WO2003083777A2 (ja) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007125240A (ja) * | 2005-11-04 | 2007-05-24 | Hitachi Medical Corp | 画像診断装置 |
JP2008526420A (ja) * | 2005-01-19 | 2008-07-24 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 画像の位置合わせのための画像処理システム及び方法 |
JP2008289632A (ja) * | 2007-05-24 | 2008-12-04 | Hitachi Medical Corp | 超音波診断装置 |
JP2011161104A (ja) * | 2010-02-12 | 2011-08-25 | Fujifilm Corp | 画像生成装置、画像生成方法、及びそのプログラム |
JP2012061169A (ja) * | 2010-09-16 | 2012-03-29 | Toshiba Corp | 画像処理装置、x線ct装置、及び画像処理方法 |
JP2013505778A (ja) * | 2009-09-25 | 2013-02-21 | ジオソフト,インコーポレーテッド | 運動情報を用いた医用画像解析のためのコンピュータ可読媒体、システム、および方法 |
JP2015150371A (ja) * | 2014-02-19 | 2015-08-24 | 株式会社東芝 | X線装置 |
JP2017515604A (ja) * | 2014-05-19 | 2017-06-15 | シロナ・デンタル・システムズ・ゲゼルシャフト・ミット・ベシュレンクテル・ハフツング | 周期運動の磁気共鳴断層撮影画像を生成する方法 |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4459678B2 (ja) | 2004-03-26 | 2010-04-28 | 富士フイルム株式会社 | 放射線動態画像取得方法および装置 |
JP4939743B2 (ja) * | 2004-11-08 | 2012-05-30 | 株式会社東芝 | X線撮像装置 |
DE102004059182A1 (de) | 2004-12-08 | 2006-06-14 | Siemens Ag | Betriebsverfahren für einen Rechner und hiermit korrespondierende Einrichtungen |
DE102005036564A1 (de) * | 2005-08-03 | 2007-02-22 | Siemens Ag | Betriebsverfahren für eine bildgebende medizintechnische Anlage und hiermit korrespondierende Gegenstände |
DE102005051323B4 (de) * | 2005-10-26 | 2008-01-31 | Siemens Ag | Verfahren und Gerät zur Magnet-Resonanz-Bildgebung auf Basis einer Gradientenechosequenz |
DE102006047719A1 (de) * | 2006-10-09 | 2008-04-10 | Siemens Ag | Verfahren und Bildgebungssystem zur Kompensation von Patientenbewegungen bei Serienaufnahmen in der medizinischen Bildgebung |
DE102007005376B4 (de) | 2007-02-02 | 2017-04-13 | Siemens Healthcare Gmbh | Verfahren zum Bestimmen von EKG-getriggerten Aufnahmezeitpunkten für die Bildgebung zur Unterstützung von interventionellen und diagnostischen Eingriffen am Herzen |
RU2009139916A (ru) * | 2007-03-29 | 2011-05-10 | Конинклейке Филипс Электроникс Н.В. (Nl) | Способ и устройство для получения объединенных рентгеновских изображений |
US8165361B2 (en) * | 2008-01-14 | 2012-04-24 | General Electric Company | System and method for image based multiple-modality cardiac image alignment |
JP5685534B2 (ja) * | 2008-08-13 | 2015-03-18 | コーニンクレッカ フィリップス エヌ ヴェ | 冠状血管及び心筋灌流情報の動的視覚化 |
US8498463B2 (en) | 2008-08-13 | 2013-07-30 | Koninklijke Philips N.V. | Mask construction for cardiac subtraction |
WO2010091273A1 (en) * | 2009-02-05 | 2010-08-12 | The Research Fondation Of State University Of New York | Method and system for transfer of image data files |
US8199994B2 (en) * | 2009-03-13 | 2012-06-12 | International Business Machines Corporation | Automatic analysis of cardiac M-mode views |
CN104519798B (zh) * | 2012-08-03 | 2018-03-23 | 皇家飞利浦有限公司 | 用于路线绘图的设备位置依赖性的叠加 |
DE102016011700A1 (de) * | 2016-09-28 | 2018-03-29 | Personal Medsystems Gmbh | Überwachung von Biosignalen, insbesondere Elektrokardiogrammen |
GB201621434D0 (en) | 2016-12-16 | 2017-02-01 | Palantir Technologies Inc | Processing sensor logs |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA1264451A (en) * | 1984-05-23 | 1990-01-16 | Nanibhushan Dattagupta | Nucleic acid probe coupled to radioactive label |
JPH0657208B2 (ja) * | 1985-12-27 | 1994-08-03 | 株式会社東芝 | デイジタル・フルオログラフイ−装置 |
DE3833309A1 (de) * | 1987-09-30 | 1989-04-20 | Toshiba Kawasaki Kk | Bildverarbeitungseinrichtung |
CN1025715C (zh) * | 1988-10-27 | 1994-08-24 | 株式会社岛津制作所 | X射线摄影装置 |
JP2790778B2 (ja) * | 1994-08-09 | 1998-08-27 | 株式会社神戸製鋼所 | 放射線画像記録読取検査装置 |
JP2927209B2 (ja) * | 1995-05-22 | 1999-07-28 | 株式会社島津製作所 | ディジタルアンギオグラフィ装置 |
US5647360A (en) * | 1995-06-30 | 1997-07-15 | Siemens Corporate Research, Inc. | Digital subtraction angiography for 3D diagnostic imaging |
US5776063A (en) * | 1996-09-30 | 1998-07-07 | Molecular Biosystems, Inc. | Analysis of ultrasound images in the presence of contrast agent |
US6056691A (en) * | 1998-06-24 | 2000-05-02 | Ecton, Inc. | System for collecting ultrasound imaging data at an adjustable collection image frame rate |
US6004270A (en) * | 1998-06-24 | 1999-12-21 | Ecton, Inc. | Ultrasound system for contrast agent imaging and quantification in echocardiography using template image for image alignment |
-
2002
- 2002-04-03 DE DE10214763A patent/DE10214763A1/de not_active Ceased
-
2003
- 2003-04-01 US US10/509,457 patent/US20050207538A1/en not_active Abandoned
- 2003-04-01 AU AU2003215822A patent/AU2003215822A1/en not_active Abandoned
- 2003-04-01 WO PCT/IB2003/001183 patent/WO2003083777A2/en active Application Filing
- 2003-04-01 JP JP2003581122A patent/JP2005521501A/ja active Pending
- 2003-04-01 EP EP03745370A patent/EP1500047A2/en not_active Withdrawn
- 2003-04-01 CN CNB038076497A patent/CN100345162C/zh not_active Expired - Fee Related
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008526420A (ja) * | 2005-01-19 | 2008-07-24 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 画像の位置合わせのための画像処理システム及び方法 |
JP2007125240A (ja) * | 2005-11-04 | 2007-05-24 | Hitachi Medical Corp | 画像診断装置 |
JP2008289632A (ja) * | 2007-05-24 | 2008-12-04 | Hitachi Medical Corp | 超音波診断装置 |
JP2013505778A (ja) * | 2009-09-25 | 2013-02-21 | ジオソフト,インコーポレーテッド | 運動情報を用いた医用画像解析のためのコンピュータ可読媒体、システム、および方法 |
JP2011161104A (ja) * | 2010-02-12 | 2011-08-25 | Fujifilm Corp | 画像生成装置、画像生成方法、及びそのプログラム |
JP2012061169A (ja) * | 2010-09-16 | 2012-03-29 | Toshiba Corp | 画像処理装置、x線ct装置、及び画像処理方法 |
JP2015150371A (ja) * | 2014-02-19 | 2015-08-24 | 株式会社東芝 | X線装置 |
JP2017515604A (ja) * | 2014-05-19 | 2017-06-15 | シロナ・デンタル・システムズ・ゲゼルシャフト・ミット・ベシュレンクテル・ハフツング | 周期運動の磁気共鳴断層撮影画像を生成する方法 |
Also Published As
Publication number | Publication date |
---|---|
EP1500047A2 (en) | 2005-01-26 |
WO2003083777A2 (en) | 2003-10-09 |
AU2003215822A1 (en) | 2003-10-13 |
CN1647111A (zh) | 2005-07-27 |
AU2003215822A8 (en) | 2003-10-13 |
CN100345162C (zh) | 2007-10-24 |
DE10214763A1 (de) | 2003-10-30 |
US20050207538A1 (en) | 2005-09-22 |
WO2003083777A3 (en) | 2004-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2005521501A (ja) | 画像シーケンスから画像を決定する方法 | |
US10736593B2 (en) | X-ray diagnostic apparatus and medical image processing apparatus | |
US8094899B2 (en) | Medical image diagnostic device | |
Ohnesorge et al. | Cardiac imaging by means of electrocardiographically gated multisection spiral CT: initial experience | |
JP4717486B2 (ja) | X線装置の画像再構成装置および対象範囲の局部的3d再構成方法 | |
US7454043B2 (en) | Image processing unit and method of associating stored images with current images | |
US10052032B2 (en) | Stenosis therapy planning | |
JP5643759B2 (ja) | デジタルサブトラクションアンギオグラフィを実施するイメージングシステム及びその作動方法 | |
US20060266947A1 (en) | Method for determining positron emission measurement information in the context of positron emission tomography | |
CN105188550A (zh) | 血管数据处理和图像配准系统、方法及装置 | |
JPWO2005089651A1 (ja) | 画像データ収集制御方法及び画像データ収集装置 | |
CN112752546B (zh) | 血管内超声成像 | |
US10362943B2 (en) | Dynamic overlay of anatomy from angiography to fluoroscopy | |
US20200315584A1 (en) | Ultrasound image generation system for generating an intravascular ultrasound image | |
JP2012000135A (ja) | マルチモダリティ動画像診断装置 | |
KR20150145106A (ko) | 의료 영상 정합 방법 및 그 장치 | |
WO2021087961A1 (zh) | 测量舒张期血流速度的方法、装置、系统及存储介质 | |
JP2018057835A (ja) | 医用情報処理装置、x線ct装置及び医用情報処理プログラム | |
CA3107473A1 (en) | Ultrasound image acquisition optimization according to different respiration modes | |
KR102254365B1 (ko) | 의료영상 처리 방법 | |
JP5534704B2 (ja) | 画像処理装置およびx線診断装置 | |
RU2676435C2 (ru) | Устройство определения полости | |
CN116704068A (zh) | 一种用于双源ct设备的图像校正方法、成像方法和系统 | |
CN111202539A (zh) | 层析图像处理设备和方法以及非暂时性记录介质 | |
Martín-Leung | Navigation and positioning aids for intracoronary interventions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060329 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060329 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090428 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20091013 |