JP3785128B2 - 画像診断装置、画像処理方法、画像処理装置及び記憶媒体 - Google Patents
画像診断装置、画像処理方法、画像処理装置及び記憶媒体 Download PDFInfo
- Publication number
- JP3785128B2 JP3785128B2 JP2002273537A JP2002273537A JP3785128B2 JP 3785128 B2 JP3785128 B2 JP 3785128B2 JP 2002273537 A JP2002273537 A JP 2002273537A JP 2002273537 A JP2002273537 A JP 2002273537A JP 3785128 B2 JP3785128 B2 JP 3785128B2
- Authority
- JP
- Japan
- Prior art keywords
- images
- image
- contrast
- subtraction
- degree
- 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
- 238000003672 processing method Methods 0.000 title claims description 3
- 239000002872 contrast media Substances 0.000 claims description 35
- 230000000694 effects Effects 0.000 claims description 25
- 238000002347 injection Methods 0.000 claims description 22
- 239000007924 injection Substances 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 19
- 238000002059 diagnostic imaging Methods 0.000 claims description 15
- 238000003745 diagnosis Methods 0.000 claims description 5
- 238000003709 image segmentation Methods 0.000 claims 2
- 238000000034 method Methods 0.000 description 25
- 238000010586 diagram Methods 0.000 description 12
- 210000004204 blood vessel Anatomy 0.000 description 10
- 238000011410 subtraction method Methods 0.000 description 9
- 210000000056 organ Anatomy 0.000 description 8
- 230000029058 respiratory gaseous exchange Effects 0.000 description 5
- 239000000470 constituent Substances 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 206010011224 Cough Diseases 0.000 description 2
- 238000012790 confirmation Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 208000031481 Pathologic Constriction Diseases 0.000 description 1
- 206010057469 Vascular stenosis Diseases 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000002603 single-photon emission computed tomography Methods 0.000 description 1
- 208000037804 stenosis Diseases 0.000 description 1
- 230000036262 stenosis Effects 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Description
【発明の属する技術分野】
本発明は、被検体に造影剤を注入し、その注入前後の画像間でサブトラクション処理を行う画像診断装置、画像処理方法、画像処理装置、記憶媒体に関する。
【0002】
【従来の技術】
造影剤の注入前に収集した画像(低コントラスト画像)と、造影剤の注入後に収集した血管コントラストが増強された画像(高コントラスト画像)との間でサブトラクション処理を行い、そのサブトラクション像を用いて行う診断は、例えば血管の狭窄位置や状態を観察するのに非常に有用であり、以前から使用されている。
【0003】
このサブトラクション処理は周知の通り、サブトラクション処理に供する画像間で位置や濃度の変化がない構造物は消去され、位置や濃度の変化のあった構造物は残存するものである。したがって、造影前後のタイミングに取得された各画像間でサブトラクション処理を行う際にこれらタイミングの間で被検体の体動がなければ、造影前後で濃度変化が生じる血管のみが抽出されることになるので、背景構造物に影響を受けることなしに血管狭窄等の観察を行い得る。
【0004】
しかしながら、これらタイミングの間で被検体の心拍や呼吸或いは咳などによって体動が生じた場合、被検体の構造物に位置の変化が生じてしまう。このタイミングでサブトラクション処理が行われると、造影された血管が映出されるのはもちろんであるが、位置の変化が生じた他の構造物も映出されてしまう。以降、この構造物によるものをアーチファクトと称する。実際の臨床の場面では、時々咳が止まらなくなる被検体もあり、その場合、横隔膜の移動によるアーチファクトが生じてしまい診断に支障をきたしていた。このような体動によるアーチファクトを除去する方法としては、例えば以下に示す技術が知られている。
【0005】
特開2001−45374号公報には、高コントラスト画像各々に対して、心拍や呼吸が同じ周期の低コントラスト画像をマスク像として選択する発明が記載されている。また、特開平8−317287号公報には、指標をマスク像及びコントラスト像に定位させ、画像シフトを行ってマスク像及びコントラスト像間の位置あわせをする発明が記載されている。
【0006】
しかし、前者の特開2001−45374号公報に記載の発明は、心拍又は呼吸に伴う体動による位置ズレの低減に関してはある程度、効果的であるが、それらの複合的な体動、また心拍や呼吸以外の体動による位置ズレの低減に関しては、効果的とはいえない。また、後者の特開平8−317287号公報に記載の発明は、体動による臓器等の平行移動による位置ズレの低減に関してはある程度、効果的であるが、サイズや形状の変位による位置ズレの低減に関しては、効果的とはいえない。
【0007】
【特許文献1】
特開2001−45374号公報(段落番号0046−0047、図2−図4)
【0008】
【特許文献2】
特開平8−317287号公報(段落番号0010−0011、図4−図6等)
【0009】
【発明が解決しようとする課題】
本発明の目的は、被検体に造影剤を注入し、その注入前後の画像間でサブトラクション処理を行う際に生じる、被検体の体動に起因する位置ズレアーチファクトを軽減することにある。
【0010】
【課題を解決するための手段】
本発明の第1局面による画像診断装置は、被検体への造影剤の注入の前後にわたって複数の画像を連続的に収集する手段と、前記複数の画像を、造影効果の無い又は比較的低い複数の第1画像と、造影効果の比較的高い複数の第2画像とに区分する画像区分手段と、前記複数の第1画像の中の特定の第1画像に対する前記被検体の体動に起因する画像上での位置ズレの類似性を示す一致度を部分領域に限局して求め、前記一致度に基づいて、前記複数の第1画像の中の一をマスク画像として前記第2画像ごとに選択するマスク選択手段と、前記複数の第2画像とそれぞれマスク画像として選択された第1画像とでサブトラクション処理する演算手段とを具備する。
本発明の第2局面による画像診断方法は、被検体への造影剤の注入の前後にわたって連続的に収集された複数の画像を、造影効果の無い又は比較的低い複数の第1画像と、血管造影効果の比較的高い複数の第2画像とに区分し、前記複数の第1画像の中の特定の第1画像に対する前記被検体の体動に起因する画像上での位置ズレの類似性を示す一致度を部分領域に限局して求め、前記一致度に基づいて、前記複数の第1画像の中の一をマスク画像として前記高コントラスト画像ごとに選択し、前記複数の第2画像から、それぞれマスク画像として選択された第1画像をサブトラクション演算する。
本発明の第3局面による画像診断装置は、被検体への造影剤の注入の前後にわたって連続的に収集された複数の画像を、造影効果の無い又は比較的低い複数の第1画像と、血管造影効果の比較的高い複数の第2画像とに区分する手段と、前記複数の第1画像の中の特定の第1画像に対する前記被検体の体動に起因する画像上での位置ズレの類似性を示す一致度を部分領域に限局して求め、前記一致度に基づいて、前記複数の第1画像の中の一をマスク画像として前記第2画像ごとに選択する手段と、前記複数の第2画像とそれぞれマスク画像として選択された第1画像とでサブトラクション演算する手段とを具備する。
本発明の第4局面による記憶媒体は、コンピュータに演算させるためのプログラムコードを記録したコンピュータ読み取り可能な記憶媒体であって、前記プログラムコードは、被検体への造影剤の注入の前後にわたって連続的に収集された複数の画像を、造影効果の無い又は比較的低い複数の第1画像と、血管造影効果の比較的高い複数の第2画像とに区分する手段と、前記複数の第1画像の中の特定の第1画像に対する前記被検体の体動に起因する画像上での位置ズレの類似性を示す一致度を部分領域に限局して求め、前記一致度に基づいて、前記複数の第1画像の中の一をマスク画像として前記第2画像ごとに選択する手段と、前記複数の第2画像とそれぞれマスク画像として選択された第1画像とでサブトラクション演算する手段とを具備する。
【0011】
【発明の実施の形態】
以下、図面を参照して本発明を好ましい実施形態により説明する。なお、本発明は造影画像を収集可能な画像診断装置を含んでいる。この種の画像診断装置としては、X線診断装置、超音波診断装置、磁気共鳴イメージング装置、さらにPET,SPECTが含まれ、そのいずれでも本発明に適用可能である。以下では、X線診断装置を例に説明する。
【0012】
図1は本実施形態に係るX線診断装置の構成を示している。画像収集部2は、寝台に載置された被検体Pを挟んで対向配置されたX線管1とX線検出器3とを有する。X線検出器3は、シンチレータとフォトダイオードアレイとを有するフラットパネルデテクタ(FPD)が採用される。しかし、X線検出器3は、イメージインテンシファイアとTVカメラとの組み合わせが採用されても良い。インジェクタ5は、被検体Pに造影剤を自動注入するための装置である。造影剤が被検体Pに注入された時、トリガ信号又はその注入時刻を表す信号が、インジェクタインタフェース11を介して画像処理ユニット7に供給される。X線検出器3から出力されるアナログの映像信号は、アナログディジタル変換器13から画像処理装置7に供給される。画像処理装置7は、これらインジェクタインタフェース11、アナログディジタル変換器13とともに、CPU9、画像データや各種制御データ等を記憶する記憶装置15、コンソール(操作卓)17、マスク選択部19、サブトラクション処理部21、ディジタルアナログ変換器23、ディスプレイ25から構成される。
【0013】
マスク選択部19は、本実施形態で特徴的な構成要素の一であり、その主な機能は、造影剤注入後に収集した造影効果の比較的高い複数の高コントラスト画像各々に対して、造影剤注入前に収集した造影効果の無い又は比較的低い複数の低コントラスト画像から最適な1枚をマスク画像として選択することにある。その選択基準としては、高コントラスト画像との間で、ある特定の低コントラスト画像を基準としてその画像に対する被検体Pの体動に起因する画像上での位置ズレが最も類似している低コントラスト画像をマスク画像として、複数の低コントラスト画像の中から選択することである。その基準により、サブトラクション画像上に、被検体Pの体動に起因する画像上での位置ズレに起因するアーチファクト(位置ズレによるサブトラクション演算残差によリ生じる偽像)の低減を実現し得る。
【0014】
以下、このマスク選択を中心に本実施形態の動作について詳細に説明する。図2は、当該動作の画像処理の流れを示すフローチャートである。画像処理に先立って画像収集が行われる。撮影にあたっては、インジェクタ5に造影剤が充填され、被検体Pに注入可能な状態にセットされる。その後、操作者(撮影技師)からの撮影開始指示をトリガとして撮影が開始される。撮影は、図3(a)、図3(b)に示すように、例えば、15フレーム/秒又は30フレーム/秒を実現する一定周期で連続的に繰り返され、それに伴って画像データが連続的に収集される。連続的に収集された画像データは、CPU9の制御のもとで、それぞれの撮影時刻又は撮影開始からの経過時間(タイムコード)が関連付けられて、記憶装置15に記憶される。この撮影期間中、予定したタイミング又は操作者が指示したタイミングで、インジェクタ5から被検体Pへの造影剤の注入が開始される。造影剤の注入開始に同期してインジェクタ5から画像処理装置7にトリガ信号が供給される。
【0015】
CPU9は、インジェクタ5からトリガ信号が供給された時刻に基づいて、記憶装置15に記憶された複数枚の画像を、造影効果の無い又は比較的低い複数の低コントラスト画像LCi と、造影効果の比較的高い複数の高コントラスト画像HCj とに区分する(S1)。具体的には、インジェクタ5からトリガ信号が供給された時刻、つまり造影剤の注入が開始された時刻から、所定の遅れ時間が経過した時刻を境界に、当該時刻前に収集された画像を低コントラスト画像LCi に分類し、当該時刻後に収集された画像を高コントラスト画像HCj に分類する。低コントラスト画像LCi の枚数は「n」、高コントラスト画像HCj の枚数は「m」と仮定する。上記遅れ時間は、操作者により、注入部位から撮影部位へ造影剤が流れるのに要する時間が推定され、その時間に設定される。
【0016】
なお、操作者による目視確認のもとで手動により指定された画像を境界としてその前後で区分するようにしてもよい。目視確認のために、CPU9の制御により、複数の画像が、撮影時間軸に沿った配列でディスプレイ25に一覧表示され、又は動画として再生される。
【0017】
次に、マスク選択部19によるマスク画像の選択処理S2〜S7が行われる。まず、S2では、複数の低コントラスト画像LCi の中から、他の低コントラスト画像LCi に対する位置ズレが最も大きい低コントラスト画像LCi が、ダミーマスク画像として選択される。ここでは、低コントラスト画像LCi の各々について、他の(n−1)枚の低コントラスト画像(LCi )との間のプラス極性のサブトラクション演算残差(+)を示す画素の数の総和を求め、その総和が最も多い低コントラスト画像(LCk )を、ダミーマスク画像として選択する。
【0018】
図4は、ダミーマスク画像選択処理の詳細な手順を示している。図5にはその補足説明図を示している。図4において、処理変数Nは、注目される低コントラスト画像LCi の移動を表す変数であり、Mは、上記他の低コントラスト画像LCi の移動を表す変数として定義される。まず、Nが1に初期化され(S1−1)、また、Mが1に、また総和変数PPN(N)が0にそれぞれ初期化される(S1−2)。低コントラスト画像LCM から、低コントラスト画像LCN がサブトラクション(サブトラクション演算)される(S1−3)。そのサブトラクション画像LCSM 上に対して、注目画素がフレーム内の始点に設定され(S1−4)、その注目画素の画素値がしきい値Thと比較される(S1−5)。しきい値Thは、例えば0に初期的に設定されるが、操作者の指示に従って0以上の任意の値に変更されることができる。注目画素の画素値がしきい値Thを超過しているとき、つまり注目画素の画素値がプラス極性を示しているとき、総和変数PPN(N)に1が加算される(S1−6)。注目画素の画素値がしきい値Th以下のとき、つまり注目画素の画素値がゼロ又はマイナス極性を示しているとき、注目画素が次の位置に移動される(S1ー8)。このS1−5乃至S1−8のループは、注目画素がフレーム内の終点に達するまで(S1−6)繰り返される。このS1−5乃至S1−8のループの繰り返しにより、対象画像LCM に対する注目画像LCN の位置ズレによりサブトラクション演算残差が生じ、しかも、その残差がプラス極性(+)を示している画素の数がカウントされる。
【0019】
ここで、プラス極性(+)のサブトラクション演算残差を示す画素が持つ意味について説明する。図6(a)には、低コントラスト画像LCi どうしのサブトラクションを示し、図6(b)には、高コントラスト画像HCj から低コントラスト画像LCi をサブトラクションした例を示している。周知のとおり、X線診断装置では、X線減衰の多いパス上の画素では、低い画素値を示し、X線減衰の少ないパス上の画素では、逆に高い画素値を示す。図6(a)において、呼吸等の体動により臓器や横隔膜等の構造物が移動した場合、サブトラクション処理後は、その移動前後で臓器が重なった部分では原理的にゼロ値又はその近似値を示し、ずれた部分ではその移動の向きに応じてプラス極性(+)とマイナス極性(−)を示す。一方、図6(b)に示すように、造影剤流入により、画素値が低下した部分の画素では、HCj−LCjとのサブトラクションを行った場合は、必ずマイナス極性(−)を示す。従って、サブトラクション後にプラス極性(+)を示す画素は、造影剤により濃度変化を起こした画素を除外して、臓器や組織の位置ズレを起こした画素であると認識することが可能となる。従って詳細は後述するが、各高コントラスト画像HCj に対して、複数の低コントラスト画像LCi の中から個々にマスク画像を選択するのは、位置ズレの程度の類似性を表す一致度という指標を用いて行うが、その一致度に対する造影剤の影響を極力除外して、その一致度の信頼性を向上させることが可能となる。
【0020】
図4に戻り、注目画素がフレーム内の終点に達したとき(S1−7)、変数Mが1だけインクリメントされ(S1−10)、S1−3のステップに帰還し、同じ注目画像LCN に対して次の画像LCM とサブトラクションされ(S1−3)、このサブトラクション画像LCSM に対して、S1−5乃至S1−8のループが繰り返される。それにより、当該次の画像LCM とサブトラクション画像LCSM 上のプラス極性を示す画素の数が、総和PPN(N)に累積される。
【0021】
このようにS1−3乃至S1−10のループが、M=nに達するまで繰り返され、それにより注目画像LCN に対する他の複数の低コントラスト画像LCi とのプラス極性の位置ズレを起こした画素の数の総和PPN(N)が、注目画像に対応するダミーマスク画像選択のためのパラメータとして求められる。
【0022】
S1−9でM=nに達したとき、注目画像LCN を決める変数Nが1だけインクリメントされ、注目画像が移動される(S1−12)。移動された次の注目画像に対して同様に、S1−3乃至S1−10のループにより総和PPN(N)が計算される。S1−11でN=n、つまり複数の低コントラスト画像LCi に対してそれぞれプラス画素の総和PPN(N)の計算が完了したとき、複数の低コントラスト画像LCi に対してそれぞれ求められたプラス画素の総和PPN(N)同士が比較され、最大値を示す総和PPN(N)が選択され、その選択された最大値を示す総和PPN(N)に対応する注目画像LCi が、ダミーマスク画像LCk として設定される(S1−13)。
【0023】
このような選択処理により、ダミーマスク画像LCk は、他の低コントラスト画像LCi との間での位置ズレの程度が最も大きい低コントラスト画像LCi としての性質を持つ。この性質は、後述する正規マスク画像を選択する際に用いられる「位置ズレの程度の類似性を表す一致度」のレンジを拡大して、正規マスク画像の選択性能を向上させることができる。
【0024】
なお、ダミーマスク画像LCk の選択は、上述のように正規マスク画像の選択性能の向上を実現する目的で行われる処理であるが、この正規マスク画像の選択性能は若干低下するものの、ダミーマスク画像LCk の選択処理を行わないで、デフォルトのダミーマスク画像LCi を例えば第1フレームの低コントラスト画像LC1 に固定させるようにしてもよい。
【0025】
図2に戻り、ダミーマスク画像LCk が決まった後、「位置ズレの程度の類似性を表す一致度」を計算する範囲(関心領域ROIという)を、被検体の体動が比較的大きい領域に、設定する処理が行われる(S3)。なお、計算範囲を、関心領域ROIに限定するのは、「位置ズレの程度の類似性を表す一致度」の計算処理の工数を減らして処理効率を高め、さらに一致度の信頼性を高める目的で行われるもので、必須の処理ではなく、「位置ズレの程度の類似性を表す一致度」の計算処理の工数を減らす必要性が低ければ、さらに一致度の信頼性を向上させる必要性が低いのであれば、計算範囲を関心領域ROIに限定する必要はない。
【0026】
S3では、低コントラスト画像LCi から、ダミーマスク画像LCk をサブトラクションし、そのサブトラクション画像LCSi からプラス極性のサブトラクション演算残差をを示す画素の群で関心領域ROIを設定する。
【0027】
図7には、関心領域ROIの設定処理ルーチンが示されている。図8にはその補足説明図を示している。図7において、処理変数Nは、低コントラスト画像LCi の移動を表す変数として定義される。まず、Nが1に初期化される(S2−1)。低コントラスト画像LCN から、ダミーマスク画像LCk がサブトラクション(サブトラクション演算)される(S2−2)。そのサブトラクション画像LCSN 上に対して、注目画素がフレーム内の始点に設定され(S2−3)、その注目画素の画素値がしきい値Thと比較される(S2−4)。しきい値Thは、例えば0に初期的に設定されるが、操作者の指示に従って0以上の任意の値に変更されることができる。注目画素の画素値がしきい値Thを超過しているとき、つまり注目画素の画素値がプラス極性のサブトラクション演算残差を示しているとき、その画素の画素位置がROIテーブルに追加される(S2−4)。注目画素の画素値がしきい値Th以下のとき、つまり注目画素が背景レベル若しくは血管レベルを示しているとき、注目画素が次の位置に移動される(S2−7)。このS2−4乃至S2−7のループは、注目画素がフレーム内の終点に達するまで(S2−6)繰り返される。このS2−4乃至S2−7のループの繰り返しにより、対象画像LCN に対するダミーマスク画像LCk の位置ズレによりサブトラクション演算残差が生じ、しかも、その残差がプラス極性(+)を示している複数の画素の位置情報がROIテーブルに記憶される。
【0028】
なお、一般的にサブトラクションを行うと、血管レベルはマイナス極性に、背景レベルで且つアーチファクトが現れていないところはゼロになる。しかし、ノイズの影響で背景レベルでも、ある程度のプラス極性を持つものが存在する。そこで、しきい値Thは、このノイズの影響を考慮し、ノイズレベル程度の値に設定する。これにより、しきい値Th以上のデータは、ほぼアーチファクト領域ということができる。
【0029】
S2−8でN=nに達していないとき、対象画像LCN を決める変数Nが1だけインクリメントされ、対象画像が移動され(S2−9)、S2−2のステップに帰還し、次の対象画像LCN からダミーマスク画像LCk がサブトラクションされ(S2−2)、このサブトラクション画像LCSN に対して、S2−4乃至S2−7のループが繰り返される。それにより、ダミーマスク画像LCk が他の複数の低コントラスト画像LCi に対してプラス極性の位置ズレを起こした全ての画素が関心領域ROIとして設定される。
【0030】
なお、関心領域ROIはコンソール17を介して操作者により手動で指定された範囲に設定するようにしてもよい。そのために、図9に示すように、CPU9の制御のもとで、低コントラスト画像LCi 各々からダミーマスク画像LCk がサブトラクションされ、得られた(n−1)枚の低コントラストサブトラクション画像LCSi が一覧で又はページ送りで順に表示される。それら表示画像を観察し、アーチファクトが最も顕著に現れている画像を見定め、その表示画像上に操作者の指示に従って長方形又は楕円形の1又は複数の図形が配置され、その又はそれらの図形に限定された閉じた領域に関心領域ROIが設定される。なお、図9では2つの関心領域ROI1,ROI2が設定された例を示している。
【0031】
図2に戻り、関心領域ROIが設定された後、上記「位置ズレの程度の類似性を表す一致度」を計算する処理が開始される。まず、図10に示すように、S4において、低コントラスト画像LCi と高コントラスト画像HCj 各々から、ダミーマスク画像LCk がサブトラクションされ、n枚の低コントラストサブトラクション画像LCSi と、m枚の高コントラストサブトラクション画像HCSj とが生成される。次に、高コントラストサブトラクション画像HCSj ごとに、n枚の低コントラストサブトラクション画像LCSi に対して個別に「位置ズレの程度の類似性を表す一致度」が関心領域ROIに限局して計算される(S5)。そして一致度に基づいて、位置ズレの程度の類似性が最も高い低コントラスト画像LCi のいずれか一が、高コントラスト画像HCj の各々に対して対応付けられる(S6)。一致度は、上述したように、低コントラスト画像LCi と高コントラスト画像HCj との間に存在する位置ズレの程度がどの程度の類似しているかを数量化して示す指標であり、低コントラストサブトラクション画像LCSi と、高コントラストサブトラクション画像HCSj との同じ位置の画素同士の画素値の比較情報(例えば画素値の差異の総和)に基づいて計算される。計算方法にもよるが、ここでは、一致度が低いとき、その低コントラスト画像LCi と高コントラスト画像HCj との間の位置ズレの程度は非常に類似しており、逆に、一致度が高いとき、その低コントラスト画像LCi と高コントラスト画像HCj との間の位置ズレの程度はあまり類似していないことを表すものとしている。
【0032】
図11には、S5乃至S6の処理の詳細な手順を示している。この処理では、注目画像を高コントラスト画像HCj として、注目画像に対してn枚の低コントラスト画像LCi に対する一致度を個別に計算し、最も低い一致度を示す低コントラスト画像LCi を、注目画像に対する正規のマスク画像として選択し、この選択処理を注目画像を移動しながら繰り返すものであり、図11において、処理変数Nは、注目画像の移動変数として定義され、処理変数Mは、一致度計算の対象とされる低コントラスト画像LCi の移動変数として定義される。
【0033】
まず、S3−1において、変数Nが1に初期化される。次に、変数Mが1に初期化され(S3−2)、また一致度を求める際の中間値DATA(M)が0に初期化される(S3−3)。続いて、注目画素が関心領域ROI内の始点に設定され(S3−4)、その位置で、高コントラストサブトラクション画像HCSN 上と低コントラストサブトラクション画像LCSM との間で画素値がサブトラクション演算され、その差異の二乗が中間値DATA(M)に加算される(S3−5)。そして、注目画素が次の位置に移動され(S3−7)、S3−5乃至S3−7のループが、注目画素が関心領域ROI内の終点に達するまで(S3−6)繰り返される。このS3−5乃至S3−7のループの繰り返しにより、対象画像LCSM と注目画像HCN との間の画素値の差異の二乗が関心領域ROI内で総和された値として中間値DATA(M)が求められる。注目画素がROI内の終点に達したとき(S3−6)、中間値DATA(M)がROI内の画素数Dで割り算され、その割り算値の平方根が計算される(S3−8)。この平方根値、つまりS3−5とS3−8との処理により得られるRMSE(Root Mean Square Error:差の二乗平均平方根値)が、注目画像HCN に対する対象画像LCM の位置ズレ程度の類似性を示す一致度COR(M)として保持される。
【0034】
そして、変数Mが1だけインクリメントされ(S3−10)、S3−3のステップに帰還しDATA(M)とCOR(M)とが0に初期化され、S3−4乃至S3−8の処理を経て、同じ注目画像LCN と次の対象画像LCM との間の一致度が、計算される(S3−8)。
【0035】
S3−4乃至S3−10のループが、S3−9でM=nに達するまで繰り返される。それにより注目画像HCN に対する複数の低コントラスト画像LCi おのおのに対する一致度CORが計算され、図12に示すように、その中の最小の一致度CORに対応する1枚の低コントラスト画像LCi が、当該注目画像HCN に対する正規のマスク画像として選択される(S3−11)。
【0036】
そして、注目画像HCN を決める変数Nが1だけインクリメントされ、注目画像HCN が移動される(S2−13)。移動された次の注目画像に対して同様に、S3−1乃至S3−11の処理により、移動した次の注目画像HCN に対して、正規のマスク画像が複数の低コントラスト画像LCi から選択される。この処理がS3−12で、N=mに達するまで繰り返される。それにより、それぞれ被検体の体動による画像上での臓器の位置ズレの程度が最も類似した低コントラスト画像LCi がマスク画像が、各高コントラスト画像HCj に対して個々に選択される。
【0037】
図2に戻り、S7において、高コントラスト画像HCj から、それぞれ正規のマスク画像として選択された低コントラスト画像LCi がサブトラクション処理部21でサブトラクションされる(図13参照)。
【0038】
このようにそれぞれ被検体の体動による画像上での臓器の位置ズレの程度が最も類似した低コントラスト画像LCi がマスク画像が、各高コントラスト画像HCj に対して個々に選択されるので、図14(a)に示すように、被検体の体動による画像上での臓器の位置ズレの程度が相違することにより残留していたアーチファクトが、図14(b)に示すように、低減され得る。
【0039】
なお、上述では、図2のS7で最終処理として、高コントラスト画像とそれぞれマスク画像として選択された低コントラスト画像をサブトラクションするように説明したが、高コントラストサブトラクション画像から、それぞれマスク画像として選択された低コントラストサブトラクション画像をサブトラクションするようにしてもよい。この場合、高コントラストサブトラクション画像に残留している位置ズレによるアーチファクトが、低コントラストサブトラクション画像に残留している位置ズレによるアーチファクトにより相殺されることになる。
【0040】
また、一致度CORとしては、図15のS3−14に示すように、高コントラストサブトラクション画像HCN と低コントラストサブトラクション画像LCM との差異の絶対値の総和として求めてもよい。また、一致度CORを計算する範囲を、図16のS3−15に示すように、関心領域内であって、高コントラストサブトラクション画像HCN の画素値が、臓器の位置ズレを示すプラス極性(+)を示す画素に限定しても良い。
【0041】
また、高コントラストサブトラクション画像HCSN に対してしきい値Thと比較し(S3−17)、しきい値Thより高い画素には固定値“1”を割り当て(S3−18)、しきい値Thと等価又はしきい値Thより低い画素には固定値“0”を割り当て(S3−19)、同様に、低コントラストサブトラクション画像LCSM に対してしきい値Thと比較し(S3−20)、しきい値Thより高い画素には固定値“1”を割り当て(S3−21)、しきい値Thと等価又はしきい値Thより低い画素には固定値“0”を割り当て(S3−22)、その割り当てられた画素の固定値の差異から一致度COR(M)を計算する(S3−14)ようにしてもよい。
【0042】
また、造影剤を血管に注入してコントラストの増強された血管動画像を観察するX線診断装置において、複数のサブトラクション方法を有している場合、撮影プログラム毎にサブトラクション方法を選択的に指定するようにしてもよい。複数のサブトラクション方法には、心電波形の位相が略一致したフレーム同士をサブトラクションする心電同期サブトラクション方法が含まれており、心電同期サブトラクション方法が選択されていて、且つ、心電波形の入力信号がない場合には、心電波形の入力信号がない旨のワーニングメッセージを表示する。ワーニングメッセージが表示された際、心電波形の入力信号が入ったことを確認するまで、または別の撮影プログラムが選択されるまで、若しくは別のサブトラクション方法が選択されるまでX線を曝射しないよう制御することが望ましい。また、ワーニングメッセージが表示された際、心電波形の入力信号が入ったことを確認するまで、または別の撮影プログラムが選択されるまで、若しくは別のサブトラクション方法が選択されるまでに一定時間を経過した場合、自動的にデフォルトのサブトラクション方法が選択されるようにしてもよい。また、複数のサブトラクション方法には、上述した被検体の体動が比較的大きい領域内で造影効果が比較的低い低コントラスト画像群と造影効果が比較的高い高コントラスト画像群との間で画素値を比較して最適なサブトラクションペアを決定する方法が含まれ、この方法が選択されている場合、造影効果が比較的低い画像群と造影効果が比較的高い画像群との境界フレームと関心領域との両方、又はどちらか一方をプリセットできるようにすることが望ましい。仮に、関心領域のプリセットデータがない場合、被検体の体動が比較的大きい領域が自動的に同定されるべきである。また、境界フレームのプリセットデータがない場合、造影剤注入時刻から所定の遅れ時間を経過した時刻に基づいて境界フレームを同定すべきである。さらに、境界フレーム及び造影剤からの遅れ時間に関するのプリセットデータが共にない場合、境界フレームを同定できない旨のワーニングメッセージが表示される。
【0043】
(変形例)
本発明は、上述した実施形態に限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で種々変形して実施することが可能である。さらに、上記実施形態には種々の段階が含まれており、開示される複数の構成要件における適宜な組み合わせにより種々の発明が抽出され得る。例えば、実施形態に示される全構成要件から幾つかの構成要件が削除されてもよい。
【0044】
【発明の効果】
本発明によれば、被検体に造影剤を注入し、その注入前後の画像間でサブトラクション処理を行う際に生じる被検体の体動に起因する位置ズレアーチファクトを軽減することができる。
【図面の簡単な説明】
【図1】本発明の実施形態に係る画像診断装置(X線診断装置)の構成を示す図。
【図2】本実施形態の動作を示すフローチャート。
【図3】図2のS1に関する補足説明図。
【図4】図2のS2の詳細な動作を示すフローチャート。
【図5】図4においてN=1の場合の動作の補足説明図。
【図6】図4のS1−5においてプラス極性の画素値(画素値>0)の画素が持つ特性に関する補足説明図。
【図7】図2のS4の詳細な動作を示すフローチャート。
【図8】図7の動作の補足説明図。
【図9】本実施形態において、関心領域ROIを手動設定する場合の操作画面の例を示す図。
【図10】図2のS5の補足説明図。
【図11】図2のS6〜S7の詳細な動作を示すフローチャート。
【図12】図11においてN=1の場合の動作の補足説明図。
【図13】図2のS8の補足説明図。
【図14】本実施形態によるマスク選択を通して生成したサブトラクション画像を従来のサブトラクション画像と比較して示す図。
【図15】図2のS6〜S7の他の動作を示すフローチャート。
【図16】図2のS6〜S7のさらに他の動作を示すフローチャート。
【図17】図2のS6〜S7のさらに他の動作を示すフローチャート。
【符号の説明】
1…X線管、
3…フラットパネルデテクタ、
5…インジェクタ、
7…画像処理装置、
9…CPU、
11…インジェクタインタフェース、
13…アナログディジタル変換器(ADC)、
15…記憶装置、
17…コンソール(操作卓)、
19…マスク選択部、
21…サブトラクション処理部、
23…ディジタルアナログ変換器(DAC)、
25…ディスプレイ。
Claims (20)
- 被検体への造影剤の注入の前後にわたって複数の画像を連続的に収集する手段と、
前記複数の画像を、造影効果の無い又は比較的低い複数の第1画像と、造影効果の比較的高い複数の第2画像とに区分する画像区分手段と、
前記複数の第1画像の中の特定の第1画像に対する前記被検体の体動に起因する画像上での位置ズレの類似性を示す一致度を部分領域に限局して求め、前記一致度に基づいて、前記複数の第1画像の中の一をマスク画像として前記第2画像ごとに選択するマスク選択手段と、
前記複数の第2画像とそれぞれマスク画像として選択された第1画像とでサブトラクション処理する演算手段とを具備することを特徴とする画像診断装置。 - 前記画像区分手段は、前記被検体への造影剤の注入開始時刻から所定時間経過した時刻に基づいて前記複数の画像を前記複数の第1画像と前記複数の第2画像とに区分することを特徴とする請求項1記載の画像診断装置。
- 前記画像区分手段は、操作者から指定された前記複数の画像の中の特定の画像を境界として、前記複数の画像を前記複数の第1画像と前記複数の第2画像とに区分することを特徴とする請求項1記載の画像診断装置。
- 前記特定の画像は、撮影時間軸に沿った配列で一覧表示された前記複数の画像から指定されることを特徴とする請求項3記載の画像診断装置。
- 前記特定の画像を指定するために、前記複数の画像が動画として再生されることを特徴とする請求項3記載の画像診断装置。
- 前記マスク選択手段は、前記部分領域として前記被検体の体動によるアーチファクトが発生している領域に限局して前記一致度を求めることを特徴とする請求項1記載の画像診断装置。
- 前記マスク選択手段は、前記第1画像各々から前記特定の第1画像をサブトラクション演算し、前記サブトラクション演算された第1画像の中で、造影剤の画素値から非造影剤の画素値をサブトラクション演算した値の極性とは逆の極性を示す画素から前記領域を構成する手段を有することを特徴とする請求項6記載の画像診断装置。
- 前記マスク選択手段は、前記第1画像各々を前記複数の第1画像からサブトラクション演算し、前記サブトラクション演算された第1画像の中で前記逆の極性を示す画素の数の総和が最も多い第1画像を、前記特定の第1画像として決定する手段を有することを特徴とする請求項7記載の画像診断装置。
- 前記被検体の体動によるアーチファクトが発生している領域は、前記第1画像から前記特定の第1画像をサブトラクション演算し、前記サブトラクション演算された第1画像上に操作者の指示に従って配置された1又は複数の図形に基づいて設定されることを特徴とする請求項6記載の画像診断装置。
- 前記図形は、長方形又は楕円形であることを特徴とする請求項9記載の画像診断装置。
- 前記マスク選択手段は、前記複数の第1画像及び前記複数の第2画像それぞれから前記特定の第1画像をサブトラクション演算し、前記サブトラクション演算された第2画像と前記サブトラクション演算された第1画像との間での画素値差異の総和に基づいて前記一致度を計算することを特徴とする請求項1記載の画像診断装置。
- 前記マスク選択手段は、前記複数の第1画像及び前記複数の第2画像それぞれから前記特定の第1画像をサブトラクション演算し、前記サブトラクション演算された第2画像と前記サブトラクション演算された第2画像との間での二乗誤差に基づいて前記一致度として計算することを特徴とする請求項1記載の画像診断装置
- 前記マスク選択手段は、前記複数の第1画像及び前記複数の第2画像それぞれから前記特定の第1画像をサブトラクション演算し、前記サブトラクション演算された第2画像と前記サブトラクション演算された第1画像との間の画素値の比較情報を前記一致度として計算することを特徴とする請求項12記載の画像診断装置。
- 前記マスク選択手段は、前記サブトラクション演算された第2画像上において、造影剤の画素値から非造影剤の画素値をサブトラクション演算した値の極性とは逆の極性を示す画素を対象として、前記一致度を計算することを特徴とする請求項12記載の画像診断装置。
- 前記マスク選択手段は、前記複数の第1画像及び前記複数の第2画像それぞれから前記特定の第1画像をサブトラクション演算し、前記サブトラクション演算された第2画像と前記サブトラクション演算された第1画像とに対して、造影剤の画素値から非造影剤の画素値をサブトラクション演算した値の極性とは逆の極性を示す画素に所定値を割り当て、その割り当てられた画素同士について前記一致度を計算することを特徴とする請求項1記載の画像診断装置。
- 前記第1画像どうしのサブトラクション画像内において画素値が所定のしきい値を超過している複数の画素を、前記部分領域として前記第1画像ごとに設定する部分領域設定部をさらに備えることを特徴とする請求項1記載の画像診断装置。
- 前記部分領域を操作者が設定操作するための操作部をさらに備えることを特徴とする請求項1記載の画像診断装置。
- 被検体への造影剤の注入の前後にわたって連続的に収集された複数の画像を、造影効果の無い又は比較的低い複数の第1画像と、血管造影効果の比較的高い複数の第2画像とに区分し、
前記複数の第1画像の中の特定の第1画像に対する前記被検体の体動に起因する画像上での位置ズレの類似性を示す一致度を部分領域に限局して求め、前記一致度に基づいて、前記複数の第1画像の中の一をマスク画像として前記高コントラスト画像ごとに選択し、
前記複数の第2画像から、それぞれマスク画像として選択された第1画像をサブトラクション演算することを特徴とする画像処理方法。 - 被検体への造影剤の注入の前後にわたって連続的に収集された複数の画像を、造影効果の無い又は比較的低い複数の第1画像と、血管造影効果の比較的高い複数の第2画像とに区分する手段と、
前記複数の第1画像の中の特定の第1画像に対する前記被検体の体動に起因する画像上での位置ズレの類似性を示す一致度を部分領域に限局して求め、前記一致度に基づいて、前記複数の第1画像の中の一をマスク画像として前記第2画像ごとに選択する手段と、
前記複数の第2画像とそれぞれマスク画像として選択された第1画像とでサブトラクション演算する手段とを具備することを特徴とする画像処理装置。 - コンピュータに演算させるためのプログラムコードを記録したコンピュータ読み取り可能な記憶媒体であって、前記プログラムコードは、
被検体への造影剤の注入の前後にわたって連続的に収集された複数の画像を、造影効果の無い又は比較的低い複数の第1画像と、血管造影効果の比較的高い複数の第2画像とに区分する手段と、
前記複数の第1画像の中の特定の第1画像に対する前記被検体の体動に起因する画像上での位置ズレの類似性を示す一致度を部分領域に限局して求め、前記一致度に基づいて、前記複数の第1画像の中の一をマスク画像として前記第2画像ごとに選択する手段と、
前記複数の第2画像とそれぞれマスク画像として選択された第1画像とでサブトラクション演算する手段とを具備することを特徴とする記憶媒体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002273537A JP3785128B2 (ja) | 2002-09-19 | 2002-09-19 | 画像診断装置、画像処理方法、画像処理装置及び記憶媒体 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002273537A JP3785128B2 (ja) | 2002-09-19 | 2002-09-19 | 画像診断装置、画像処理方法、画像処理装置及び記憶媒体 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2005024123A Division JP4150005B2 (ja) | 2005-01-31 | 2005-01-31 | 画像診断装置、画像処理方法、画像処理装置及び記憶媒体 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2004112469A JP2004112469A (ja) | 2004-04-08 |
JP3785128B2 true JP3785128B2 (ja) | 2006-06-14 |
Family
ID=32270268
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002273537A Expired - Lifetime JP3785128B2 (ja) | 2002-09-19 | 2002-09-19 | 画像診断装置、画像処理方法、画像処理装置及び記憶媒体 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3785128B2 (ja) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4612325B2 (ja) * | 2004-04-09 | 2011-01-12 | 株式会社東芝 | 超音波診断装置 |
GB2445322B (en) * | 2004-08-13 | 2008-08-06 | Stichting Tech Wetenschapp | Intravasular ultrasound techniques |
JP2006087631A (ja) * | 2004-09-22 | 2006-04-06 | Sangaku Renkei Kiko Kyushu:Kk | 画像診断装置、画像処理装置、及び画像処理プログラムを記録した記録媒体 |
JP4828920B2 (ja) * | 2004-11-24 | 2011-11-30 | 株式会社東芝 | 3次元画像処理装置 |
JP4745080B2 (ja) | 2006-02-20 | 2011-08-10 | 猛 中浦 | X線診断装置、画像処理装置及びプログラム |
JP5161427B2 (ja) | 2006-02-20 | 2013-03-13 | 株式会社東芝 | 画像撮影装置、画像処理装置及びプログラム |
JP5263867B2 (ja) * | 2007-10-15 | 2013-08-14 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波撮像装置 |
JP5606143B2 (ja) * | 2009-06-08 | 2014-10-15 | 株式会社東芝 | 超音波診断装置、画像処理装置、画像処理方法および画像表示方法 |
JP5523782B2 (ja) * | 2009-09-25 | 2014-06-18 | 株式会社東芝 | X線診断装置及びx線診断方法 |
JP5645399B2 (ja) * | 2009-12-18 | 2014-12-24 | キヤノン株式会社 | X線画像処理装置、x線画像処理方法、及びコンピュータプログラム |
JP5641792B2 (ja) * | 2010-06-24 | 2014-12-17 | 株式会社東芝 | 医用画像診断装置及び医用画像診断装置の制御方法 |
RU2605524C2 (ru) * | 2011-12-23 | 2016-12-20 | Конинклейке Филипс Н.В. | Магнитно-резонансная визуализация с подавлением артефактов потока |
US9384555B2 (en) * | 2012-07-02 | 2016-07-05 | Kabushiki Kaisha Toshiba | Motion correction apparatus and method |
CN106821404A (zh) * | 2017-01-20 | 2017-06-13 | 北京东软医疗设备有限公司 | 血管造影方法和系统 |
CN111913873A (zh) * | 2020-06-17 | 2020-11-10 | 浙江数链科技有限公司 | 图片校验方法、装置、系统和计算机可读存储介质 |
-
2002
- 2002-09-19 JP JP2002273537A patent/JP3785128B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JP2004112469A (ja) | 2004-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3785128B2 (ja) | 画像診断装置、画像処理方法、画像処理装置及び記憶媒体 | |
US7551721B2 (en) | X-ray diagnostic apparatus, image processing apparatus, and program | |
US7840255B2 (en) | X-ray CT apparatus and myocardial perfusion image generating system | |
CN108703764B (zh) | 血管造影方法、装置、系统、设备及存储介质 | |
US7853309B2 (en) | X-ray CT apparatus and myocardial perfusion image generating system | |
JP5161427B2 (ja) | 画像撮影装置、画像処理装置及びプログラム | |
US8111891B2 (en) | X-ray CT apparatus, myocardial perfusion information generating system, X-ray diagnostic method and myocardial perfusion information generating method | |
JP6021420B2 (ja) | 画像処理装置、画像処理方法およびプログラム | |
JP6772873B2 (ja) | 動態解析装置及び動態解析システム | |
JP4857418B2 (ja) | 検体のモーションアーチファクト補正の方法およびシステム | |
JP2007167656A (ja) | 対象の運動の解析方法および断層撮影装置 | |
JP4150005B2 (ja) | 画像診断装置、画像処理方法、画像処理装置及び記憶媒体 | |
JP5029702B2 (ja) | 画像生成装置、画像生成方法およびプログラム | |
CN113962949A (zh) | 脑部医学图像处理系统、方法及计算机可读存储介质 | |
JP2006087631A (ja) | 画像診断装置、画像処理装置、及び画像処理プログラムを記録した記録媒体 | |
CN114533096B (zh) | 一种脑血管造影中伪影去除方法及系统 | |
JP5229496B2 (ja) | 画像処理装置 | |
JPS596042A (ja) | 画像処理装置 | |
JPH06125499A (ja) | ディジタル・サブトラクション・アンギオグラフィ装置 | |
CN116725566A (zh) | 基于高速三维结构光传感器的ct成像质量提升方法 | |
Backfrieder et al. | Model based LV-reconstruction in bi-plane x-ray angiography | |
JPH04252578A (ja) | X線サブトラクションシステム | |
JPH02237278A (ja) | ディジタル・サブトラクション・アンギオグラフィ装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20040806 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20041207 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050131 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050913 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20060314 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20060316 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 Ref document number: 3785128 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100324 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100324 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110324 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120324 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130324 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140324 Year of fee payment: 8 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313115 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
EXPY | Cancellation because of completion of term |