JP2004208807A - マイオカーディアルブラッドボリュームマップによる血液量解析・表示方法 - Google Patents
マイオカーディアルブラッドボリュームマップによる血液量解析・表示方法 Download PDFInfo
- Publication number
- JP2004208807A JP2004208807A JP2002379760A JP2002379760A JP2004208807A JP 2004208807 A JP2004208807 A JP 2004208807A JP 2002379760 A JP2002379760 A JP 2002379760A JP 2002379760 A JP2002379760 A JP 2002379760A JP 2004208807 A JP2004208807 A JP 2004208807A
- Authority
- JP
- Japan
- Prior art keywords
- line
- myocardium
- myocardial
- region
- blood volume
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/481—Diagnostic techniques involving the use of contrast agent, e.g. microbubbles introduced into the bloodstream
Abstract
【解決手段】気泡性超音波造影剤を血管に注入して超音波診断装置によって得られた心筋と心腔のコントラストエコー画像を評価する血液量解析・表示方法において、計算領域を細分化して心筋、心腔の音場が同等と思われる部位で、心筋の輝度と心腔の輝度の差分を求め、前記差分位で心筋の血液量を求め、その血液量に対してカラー表示を行う手段、を備えたことを特徴とするマイオカーディアルブラッドボリュームマップによる血液量解析・表示方法である。
【選択図】 図3
Description
【発明が属する技術分野】
本発明は、気泡性超音波造影剤を血管に注入して超音波診断装置によって得られるコントラストエコー画像に対して、心筋部と心腔部の輝度を手がかりに、心筋部と心腔部の画像を解析し、心筋の血液量を視認的、定量的に解析して、虚血性心疾患などの発見や治療の手助けとなるコントラストエコー画像解析法に関する。
【0002】
【従来の技術】
X線診断装置、X線CT装置、MRI装置、核医学診断装置などは、人体の内部組織を外科的な手術を用いずに診断できるという特徴を持っている。しかしその一方、被曝や検査時間を要するなどの問題があった。又カテーテル検査はカテーテルを血管内に挿入するという侵襲的な検査である。その点、超音波を利用した超音波診断装置は音(波長の短い音波)を用いるために、前者の装置によって生じる、核被曝や放射線被曝、あるいは磁気よる体内組織への影響がほとんどないという特徴を持っている。すなわち、超音波診断装置は無侵襲性検査装置であり、リアルタイムの画像表示が可能であり、短期間の反復使用を行っても危険性は少なく、短い時間で検査を行えるなどの特徴を有している。また装置が小型で安価なために、各病院で導入がしやすいだけでなく、医師が直接装置を操作しながら関心領域(ROI:Region Of Interest)を綿密に検査することが容易になっている。
【0003】
超音波診断装置では、超音波プローブを体表に当て、振動子から発信される超音波の反射波を受信することによって、体内組織の断面画像を得ることができる。このため、心臓、腹部、乳腺、子宮内の胎児の動きなどの観察がリアルタイムに行える。またパワードプラ法を用いることによって、血流量のイメージイングが可能となっている。
【0004】
全身を循環してきた血液は、静脈から右心房、右心室と送られ、右心室から動脈を通して肺に送られ、肺で酸素が供給される。肺で供給された酸素を含んだ血液は左心房、左心室を通り、動脈を通して全身に送られる。このとき、一部の血液は動脈から枝分かれした冠状動脈(冠動脈)へと送られる。冠動脈は心筋を蜘蛛の巣状に覆い、心臓の活動のもととなる酸素や栄養分を供給している。この冠動脈が血栓や血管の異常によって詰まったり細くなったりすると、心筋への血液量が不足し、心臓の働きに悪影響を及ぼすことになる。冠動脈が閉塞して筋細胞が壊死に陥ると、心筋細胞の収縮力が消失する。これが心筋梗塞である。閉塞までいかなくても、血管が細くなるなどの原因で心筋への血液供給量が減少した場合にも、心臓への悪影響が生じる。これが狭心症である。これらの心筋梗塞や狭心症を虚血性心疾患という。
【0005】
超音波診断装置による心臓などの検査では、造影剤を注入して血流動態評価が行われている。静脈からの造影剤の注入は侵襲性が小さいために、血流動態の評価法として使用されている。この診断では、診断部位における造影剤の空間分布の時間変化の様子を輝度変化の広がりや輝度増強の程度を観察したり、造影剤の注入から関心領域に造影剤が到達するまでの時間および関心領域の造影剤によるエコー輝度の時間変化TIC(Time Intensity Curve)などを求めることが行われている。
【0006】
生体臓器からの超音波の反射波(エコー信号)は比較的非線形挙動を起こしにくい。しかし、微小気泡(マイクロバブル)を主成分に含む超音波造影剤を使用して得られるエコー信号は非線形挙動による非基本波成分が含まれる。そのため、非基本波成分のみを分離して検出すれば、生体臓器と空洞部(血管内部や心臓の心腔部など)を良好なコントラスト比の造影剤画像として得ることができる。
【0007】
以上のようなコントラストエコー法を用いた公開特許としては、『超音波診断装置および超音波診断方法』(特許文献1)や『超音波診断装置及び超音波診断方法』(特許文献2)などがある。前者では、超音波パルスの送信音圧を最適に制御して造影剤による輝度増強を高めることによって、より効果的なコントラストエコー画像を得る方法について触れている。後者は、微小気泡を主成分とする超音波造影剤を用いたコントラストエコー法を実施する際の、操作者の労力や負担を軽減させるための装置とその方法に関するものである。微小気泡の消失情報をスピーカを通して音情報として診断者に知らせるなどの方法が提唱されている。
【0008】
【特許文献1】
特開平11−155858号公報
【特許文献2】
特開2001−178722号公報
【0009】
【発明が解決しようとする課題】
微小気泡を主成分とする超音波造影剤を用いたコントラストエコー法による超音波診断装置による心臓などの検査では、コントラストエコー画像から、心筋梗塞や狭心症等の虚血性心疾患を読みとるには熟練を要する。それは以下のような理由による。
【0010】
気泡性超音波造影剤を血管に注入し、体表に当てた超音波プローブから発信した超音波を関心領域(ROI)に当てると、造影剤に含まれる気泡で超音波が反射される。造影剤を含む血液部分は反射されるエネルギーが高く、輝度情報に変換して画像表示すると明るく映る。血液が十分に行き届いていないと、その輝度が暗くなる。このため、心筋部の輝度を比較すれば、心筋に十分血液が行き届いているかどうかの判断ができるわけである。
【0011】
ところが、プローブから発信された超音波は体表から遠くなればなるほど(体内深部ほど)、超音波信号の減衰が生じる。また、プローブはフォーカスを設定することができるが、その場合、他の場所よりも焦点周辺部がエネルギーの高い状態となる。さらに、パワードプラやBモードから得られるカラー画像の場合には、心腔部のカラーが心筋部に滲み出すなどの問題がある。このような理由から、心筋部、心腔部のコントラストエコー画像から血液の灌流状態を正しく判断することは難しい。
【0012】
以上のように、従来のコントラストエコー画像では、超音波信号の深さ方向の減衰やフォーカスの辺りでのエネルギーの最大化によって、輝度にばらつきが生じていた。ところが従来の画像処理では、この音場やフォーカスの違いによる輝度変化の補正がなされないまま計測がなされてきた。このため、従来の方法で計測される数値では、各部位における数値の差が音場の違いによって起こるものか、心腔部、心筋部の灌流状態により起こるものか、判定しづらかった。
【0013】
以上の問題点を別にしても、コントラストエコー画像から心腔部、心筋部の灌流状態を読み取ることは検者の力量(熟練度)に依存し、誰もが容易にかつ正確に把握することが難しい状態にある。なぜなら、客観的に画像処理がなされていないうえ、定量的に解析されてこなかったからである。
【0014】
以上の点に鑑み、本発明が解決しようとする課題は、超音波診断装置より得られたコントラストエコー画像から心筋の血液灌流状態を解析し、視覚的かつ定量的に検者に、前記の血液灌流状態を判断する情報を提供する方法を開発することである。
【0015】
【課題を解決するための手段】
上記の発明が解決しようとする課題を解決するために、本発明では以下のような医学的な事実を考慮する。
1)心腔(心房、心室)部には血液が充満している。
2)一般に血液の灌流障害は、心筋部の心外膜側よりも心内膜側から発症する傾向にある。
また、超音波診断装置から得られるコントラストエコー画像には以下の問題点を含んでいる。
a)超音波の特性上、心外膜方向の画像は心内膜方向の画像に比べて画像が鮮明でない場合がある。
b)パワードプラモードの場合、心腔内の色が心外膜方向に滲み出していることがある。
c)コントラストエコー画像の輝度は、深さ方向の音場の違いやフォーカス位置の違いによって、一定していない。
d)ノイズによる輝度のばらつきがある。
以上の点を考慮して、本発明では同じ音場と思われる位置で心筋部と心腔部を細分化して比較・検討することによって、上記の問題点を解決する。具体的には以下のような解析法とする。
【0016】
気泡性超音波造影剤を血管に注入して超音波診断装置によって得られた心筋部と心腔部のコントラストエコー画像を評価する心筋部の血液灌流状態を解析する際に、計算領域を細分化して心筋部、心腔部の音場が同等と思われる部位で心筋部の輝度と心腔部の輝度の差分を求め、前記差分値で心筋部をカラー表示する。これが、本発明のマイオカーディアルブラッドボリュームマップによる血液量解析法の基本である。なおマイオカーディアルブラッドボリュームマップ(Myocardial Blood Volume Map、以下“MBVM”と記述)とは、心筋の血液量分布図の意味で、心筋の血液の灌流状態を色で表示分けした画像のことである。
【0017】
本発明では、また、定量的に心筋部の血液灌流状態を以下の方法で数値化する。MBVMによる血液量解析法においては、コントラストエコー画像より得られる心筋部の輝度をI1、心腔部の輝度をI2、および定数をKとしたとき、アンプリチュードA(Amplitude A)を、
Klog10A=I1−I2
なる式から求め、現状の心筋の血液量V(単位ml/100g:Blood Volume)を
V=100×A(ml/100g)
と求める。この式によって、関心領域の心筋部における血液量が算出される。なお、定数Kはダイナミックレンジ(Dynamic Range)の値によっては10とか20などの値を取る。
【0018】
MBVMを得るために、本発明では以下の3通りの方法を用意する。まず、等間隔スライディングウィンドウ(Sliding Window)法について説明する。心筋部と心腔部との境界である心内膜、および心筋部の外側の膜をそれぞれトレースし、前記心内膜線に対してd2(1〜20mm)の間隔で内側(心腔側)および外側(心筋側)に、心内膜線に沿って線を描画する。前者を内線、後者を外線と呼ぶ。
【0019】
さらに深部方向に対して音場が同等と思われる境界線(音場線、d1<1〜20mm>)を描画し、心内膜線と内線および音場線で囲まれる領域を心腔部計算領域とし、心内膜線と外線および境界線で囲まれる領域を心筋部計算領域とする。前記の計算領域に対して、心筋部計算領域の平均輝度から心腔部計算領域の平均輝度を引き、その差分値に相当する色をカラーテーブルから選択して、音場線で区切られる心筋全体を色づけする。
【0020】
音場線をd1ずつずらして前記の差分計算を繰り返し、全心筋を色づけする。音場線と音場線との間隔をd1に取り、計算領域の輝度平均を取っているのは、ノイズの影響を極力押さえるためである。また心筋全体を、心筋部計算領域の値で代表して色づけしているのは、心外膜方向では輝度情報が正確に得られない場合があることと、血液灌流障害は心内膜の方向から起きるという事実を考慮したからである。
【0021】
重複スライディングウィンドウ法は、基本的に等間隔スライディングウィンドウ法と同じ計算領域の細分化を行う。異なる点は、より音場領域の音場差(輝度差)の影響を縮小するために、音場線を微小間隔d3で垂直軸に沿ってずらしながら心筋部計算領域の輝度と心腔部計算領域の輝度との差を計算し、計算した差分値で心筋部を色づけする点である(例えばd3<=0.5mm>)。この方法では、垂直軸にずらす間隔が短いために、計算領域が部分的に重なる形になる。色は、順次新しい計算値(差分値)で塗りつぶすことになる。音場線の間隔をd1より小さな間隔d3ずつずらしているが、計算領域はd1の範囲としているのは、ノイズの影響が現れないようにするためである(計算領域が狭ければ狭いほど、ノイズの影響を受けるため)。
【0022】
ピクセルバイピクセル(Pixel by Pixel)法も計算領域の取り方は前記のスライディングウィンドウ法と同じであるが、以下の点で異なる。すなわち、心内膜線と心外膜線および音場線で囲まれる領域の各ピクセル(画素)ごとに輝度を求め、この輝度と心腔部計算領域の平均輝度を引いて輝度の差分を求め、心筋部内の色づけを前記差分値で画素単位で行う。この方法では、心筋部の心内膜から心外膜までの範囲を元のコントラストエコー画像をそのまま評価することができる。ただし、ノイズによる不具合もそのまま表示することになるのと、心外膜付近の画像が不鮮明であるという問題点もそのまま表示することになる。
【0023】
前記のスライディングウィンドウ法やピクセルバイピクセル法による心筋の色づけ方法において、パワーモードのようなカラー画像では心腔部の色が心内膜から心筋部にかけて滲み出している場合があるので、その場合には滲み出しによる正確な心筋部の計測が行えないことがある。そこで、この滲み出しによる影響を排除するために、心内膜線の内側と外側に、ある一定の間隔でバッファ線を引き、バッファ線と外線、バッファ線と内線とに囲まれる部分を計算領域とする手段を用意する。バッファ線を引く間隔はd4(=0.5mm)単位で可変に引けるようにしておく。
【0024】
【発明の実施の形態】
本発明の実施の形態を図を用いて説明する。本発明は超音波診断装置より得たコントラストエコー画像をもとに心筋部の血液灌流状態を解析するものである。解析に使用するハードウェアはコンピュータ装置である。超音波診断装置より得られたコントラストエコー画像をCDやMOなどの大容量外部記憶媒体に保存し、その外部記憶媒体の画像をもとに画像解析する。
【0025】
図1は、超音波プローブによる走査の様子を示したものである。体表に当てた超音波プローブ2の先端部210から発信された超音波b1、b2、……bn(実際には連続的に発信)は対象物で反射され、反射波は再びプローブで受信され、デジタル化されてコントラストエコー画像として編集されてメモリあるいは記憶装置に記録される。図で示した浅部方向が体表の方向であり、深部方向が体内内部方向になる。受信した反射波(エコー)は輝度情報変換され、パワードプラ法などではカラー画像情報の形で得られる。先端部210から等距離rにある場所が同一の音場(超音波の強さが等しい場所)となる。図のl0は等音場線を表している。
【0026】
超音波のエコー画像では深さ方向(図1の深部方向)では輝度の減衰が起こる。またフォーカスの辺りでエネルギーが最大になる。このため、従来のように輝度の絶対値で画像表示したのでは、正しく画像評価ができない。そこで本発明では、音場が同等と思われる部位の心筋部と心腔部の輝度の差分で心筋部を色づけする。以下で具体的に心筋部の色づけ操作について説明する。なお以下の操作は、パワーモードのみでなく、サードハーモニック、1.5ハーモニック等の白黒画像に対しても同様に扱うことができる。なぜなら、輝度情報のみで解析を行うからである。
【0027】
図2は、心筋部と心腔部のコントラストエコー画像の例である。パワーモードではカラー画像、サードハーモニック、1.5ハーモニックでは白黒画像となっている。画像は以下のタイミングで撮影されたものである。気泡性超音波造影剤は、静脈から注入される。動脈からの注入する方法もあるが、侵襲性は静脈注射の方が小さいために、一般に静脈からの注入法が用いられる。造影剤が関心領域に到達し、十分に満たされたと思われるタイミングを計って撮影される。心臓部の場合には、心電図(ECG:Electro Cardio Gram)の波形でタイミングが計れる。
【0028】
装置によっては1回のタイミングで2枚撮影する。最初の1枚目は造影剤が満たされた状態であり、2枚目は1枚目を撮影した直後(約100msec後)に同じ部位を撮影する。気泡を含む造影剤に超音波が当たると、気泡の大半は消失する。このため、1枚目の画像と2枚目の画像を比較することによって、1枚目の画像のノイズの状態や、造影剤が予測した通り関心領域に充満していたかなどの判断がつく。例えば、1枚目と2枚目の画像が同じ場合には、造影剤が行き届いていないうちに1枚目が撮影された可能性があり、解析の対象とならないことを意味する。図の左画像は1枚目の画像、右画像は2枚目の画像である。ただし、以下でサンプルとして挙げた画像(画面)は、左画像が解析のための作業画像であり、右画像はオリジナルの画像(未処理の左の画像)である。
【0029】
図3は、図2の画像に対して心内膜と心外膜のトレース線を入力し、輪郭線と境界線(音場線)を描画した図である。輪郭線はマウスなどのデバイスで手入力する。心内膜を表す輪郭線を「心内膜線」と呼び、心外膜を表す輪郭線を「心外膜線」と呼ぶことにする。等間隔スライドウィンドウ法における心筋部の色づけは、次の手順で行う。
【0030】
1)オリジナルのコントラストエコー画像に対して、心内膜111(心内膜線)と心外膜112(心外膜線)をトレースする。この心内膜線と心外膜線の間が心筋110である。人手による入力はこのトレースだけであり、以下はプログラムが自動的に処理する。
2)心内膜の内側と外側に、心内膜に沿ってd2の間隔で内線l1と外線l2を描画する。
3)心尖部210(図1の先端部)より長軸の深さ方向(図3の深部方向)に対してd1の間隔で垂直線l0k(k=0〜n、nは分割数、図では横方向の線)を描画する。このd1間隔で引いた垂直線で挟まれる範囲をほぼ同一の音場であるとして計算するため、この垂直線l0を「等音場線」あるいは単に「音場線」と呼ぶことにする。また心内膜111、内線l1、外線l2および音場線l0で囲まれる領域が計算領域である。心内膜線111、外線l2および音場線l0kとl0k+1に囲まれる領域(心筋部計算領域)の平均輝度Ak、ならびに心内膜線111、内線l1および音場線l0kとl0k+1に囲まれる領域(心腔部計算領域)の平均輝度Bkを求める。
4)心筋部の輝度と心腔部の差分△kは
△k=Ak−Bk
と計算される。この差分の値△kに相当する色をカラーテーブルcより選択し、音場線l0kとl0k+1で挟まれる心筋部全体を、選択した色で色づけする。図4の斜線部がカラーテーブルより選択した色で色づけした部分である。
5)以下順次、計算領域をずらせて(輝度の添え字kを1ずつ増やして)差分を計算し、全心筋部を差分値に対応する色で色づけする。
【0031】
なお本発明ではd1を1〜20mm、d2を1〜20mmの幅に取る。d1を1〜20mmとしているのは、この範囲を同一音場として判断できるからである。また等音場線は、図1に示したように先端部210から等距離のl0の線上であるが、図3の例では等音場線を直線で近似している。より正確に同一音場を要求する場合には、図5のように画像の頂点を中心点にして深部方向に対しては円弧で分割すればよい。対比する計算領域(AkとBk)が狭い範囲であり、隣接していることから、直線で囲まれる領域を等音場として扱っても、それほど問題にはならない。ただし心尖部210は範囲が広くなるために、図3に示すように心尖部(図の例ではn分割目)を点Oを中心に放射状に何分割化して、さらに細分化エリアとする。図3は4分割したときの例である。このとき、差分△njは
△nj=Anj−Bnj (j=1〜4)
と求める。
【0032】
次に重複スライディングウィンドウ法を説明する。重複スライディングウィンドウ法も基本的に等間隔スライディングウィンドウ法と同じであるが、計算領域のスライド(ずらし)に違いがある。すなわち図6で示すように、計算領域Ak、Bkの取り方は等間隔スライディングウィンドウ法と同じであるが、計算領域を心内膜111に沿って微小間隔d3(例えば、1mm)ずつずらしている点である。これによって深さ方向に対する音場の違いをより細かに補正することができる。
【0033】
なお、計算領域を深さ方向にd3間隔で取っていないのは、課題を解決するための手段で述べたように、ノイズの影響を排除するためである(部分的にノイズがあっても、広範囲<といっても、最大でも20mm四方の狭い範囲ではあるが>の領域で平均化することによってノイズの影響を抑えているため)。また、図の例では計算を深部方向から行っているが、先端部(浅部)方向から行ってもよい。この点はプログラミング上の問題であって、基本的に本発明の関与する範疇ではない。
【0034】
心筋部の色を心内膜付近の心筋部計算領域の計算で代表して塗りつぶしているのは、次のような理由による。すなわち先にも述べたように、心外膜付近では画像が鮮明でない場合があることと、血液灌流異状は心内膜付近から始まる傾向にあるという医学的事実により、心筋の血液灌流状態の異変を判断するためには、上記の方法を用いても問題はないという理由による。
【0035】
もし、心内膜方向と心外膜方向とを独立して計算し色分けするには、図7に示すように心内膜部(Ak)と心外膜部(Ak')との2領域に分けて計算してもよい。さらに、もっと心筋の厚さ方向(図の縦軸方向)に細分化して計算領域を取ってもよい。いずれにしろ、差分の対象になる輝度はBkである。
【0036】
ピクセルバイピクセル法も、計算領域の取り方は等間隔スライディングウィンドウ法と同じであるが、心筋部計算領域を心内膜と心外膜の間まで拡大し、しかも差分計算ではピクセル(画素)単位で行う点で異なる。図8は、ピクセルバイピクセル法による計算領域を示した図である。心腔部計算領域は等間隔スライディングウィンドウ法と同じBkであるが、心筋部計算領域は音場線l0k、l0k+1および心内膜111と心外膜112とで囲まれる範囲であり、計算の対象はピクセルAkj(jは心筋部計算領域に含まれるすべてのピクセル)であり、差分△kjは
△kj=Akj−Bk
と計算する。差分△kjに対応するカラーテーブルの色を選択して、対応するピクセルを色づけする。これによって、心筋の精密な状態が色で判断することができる。ただし、心筋部計算領域にノイズがある場合には、そのノイズもそのまま計算して色づけすることになる。
【0037】
以上、3通りの色づけについて説明してきたが、パワーモードなどでのカラー画像では、心腔部の色が心内膜から心筋部内にかけて滲み出すこと(ブルーミング<blooming>)があるので、この影響を排除する必要がある。そこで本発明では図9に示すように、心内膜111からd4間隔を置いて、心内膜の内側と外側にそれぞれバッファ線b1、b2を引く。心腔部計算領域Bkは音場線l0k、l0k+1、バッファ線b1、内線l1に囲まれる領域とし、心筋部計算領域Akは音場線l0k、l0k+1、バッファ線b2、外線l2に囲まれる領域とすることによって、ブルーミングによる影響を排除することができる。もちろん、ピクセルバイピクセル法の場合には、心筋部計算領域はバッファ線b2から心外膜までの間となる。なおd4の値は0.5mm単位で可変で選べるようにしておけば、バッファ線は自動的に描画される。
【0038】
図10の画像が、実際に上記の方法で色づけした心臓の状態を表している。この色が心筋における血液灌流状態を表している。色が暗くなっていく(カラーテーブルcの下の方の色)ほど、灌流状態が悪いことを表す。従って、色づけされたカラー画像をみれば、一見して灌流状態が把握できる。
【0039】
以上の色づけは、視覚的に血液の灌流状態を示すものである。これに対して、実際の血液量を数値化する方法が以下の数式である。コントラストエコー画像の関心領域より得られる心筋部の輝度をI1、心腔部の輝度をI2、および定数をKとしたとき、アンプリチュードAを、
Klog10A=I1−I2 ・・・・・・・・・・・・・・ (式1)
なる式から求める。ここで、I1、I2は上記のAkとBkに対応している。またアンプリチュードAは、心腔の血液量100ml/100gとしたときに対する、現時点(コントラストエコー画像)の心筋の血液流量比を表す(実際に観測されているのは気泡の濃度比であるが、これは血液量に比例しているから、同等のものとして扱っても差し支えない)。なおKは状況によって決まる定数であり、ダイナミックレンジの値によっては、10とか20などの値を取る。
【0040】
従って正常時、関心領域の心筋部における血液量V(単位ml/100g)は
V=100×A(ml/100g)
と計算できる。例えば、
輝度差△=-10dB (△=I1−I2) および K=10
としたとき、式1より
A=1/10
となる。従って関心領域の血液量Vは
V=100・A=10ml/100g
となる。心筋内の血液量は心腔内の血液量より少ないから、この値は正常であることを示している。
【0041】
【発明の効果】
コントラストエコー画像から心筋内の血液灌流状態を判断するには経験を要し、専門の医師でも簡単には判断することが難しかった。しかもコントラストエコー画像は超音波診断装置から来る問題として、ノイズが含まれていたり、深さ方向で音場(超音波の強さ)の減衰や、フォーカス付近で輝度が強調されるなどの問題を含んでいる。ところが従来の解析法では、このような問題をなんら補正することなく放置されていたために、経験だけが頼りという状態が続いていた。
【0042】
本発明では上記の問題を、等音場での、心腔内の輝度と心筋内の輝度の差分を取ることによって解決した。すなわち、音場の減衰やフォーカス付近の強調があっても、等音場の心腔、心筋を計算領域として差分を取ることによって、場所による減衰や強調から発生する問題を解消している。さらに部分的に発生するノイズの問題は、計算領域(ほぼ等音場と思われる狭い範囲の領域)の輝度の平均値化によって解決している。またパワーモードなどのカラー画像で発生する滲み出しの問題も、バッファ線を設けることによって解決している。このような理由から、本発明の解析法で得られる、色づけされた心筋の画像をみれば、一見して心筋の血液灌流状態がわかる。というのも、カラーテーブルの色は血液の灌流状態を示しているからである。
【0043】
以上の解析を行うには、検者(一般に医師)は心筋の心内膜線と心外膜線を引くだけでよく、あとはプログラムがすべて処理するから、誰にでも簡単に利用でき、なおかつ誰にでも心筋の血液灌流状態を判断することができる。また本発明の血液量計算式を用いれば、関心領域の血液量を数値として算出されるために、定量的にも血液灌流状態を判断することができる。
【0044】
色づけに対しては、3通りの方法を提唱した。スライディングウィンドウ法では、心内膜付近の灌流状態を色で表現する。等間隔スライディングウィンドウ法の場合には計算速度が速く、重複スライディングウィンドウ法では連続的に心筋の血液灌流状態が表現される。スライディングウィンドウ法では心内膜付近の血液灌流状態で心筋全体を代表しているが、発明の実施の形態でも触れたように、計算領域を心筋の厚さ方向(心内膜から心外膜方向)に細分化して差分を取れば、心内膜から心外膜の方向にも正確な血液灌流状態を表示することもできる。
【0045】
ただし、先にも述べたように、心外膜付近は鮮明な画像が得られない場合があるので、現時点では、本発明の請求項に記述した心内膜の計算領域の差分値で心筋の全体を代表してもさほど問題はない。というのも、灌流障害は心内膜の方から起きやすいという事実により、心内膜付近の解析が正確に行われれば、医学的にそれほど問題にはならないからである。
【0046】
ピクセルバイピクセル法では、心筋全体をピクセル単位で評価するために、上記のような問題はない。すなわち、心内膜から心外膜まで、オリジナルのコントラストエコー画像を正確に評価することができる。ただし、心筋内のノイズや心外膜付近の不鮮明さもそのまま評価して色づけされるために、現時点では必ずしてもスライディングウィンドウ法より優れているとは一概にいえない。従って将来、ノイズがなく、しかも心外膜付近まで鮮明なコントラストエコー画像が得られるときには、詳細な評価法として有効活用ができることになろう。
【0047】
本発明の解析・表示方法はコントラストエコー画像さえあれば(得られれば)、ごく普通に使われているパソコンでも利用できる。このため、誰にでもコストの掛からない画像解析ソフトとして導入することができる。また操作もいたって簡単なために、パソコンに不慣れな医師にも有効なツールである。
【図面の簡単な説明】
【図1】超音波プローブによるコントラストエコー画像を撮影する方法を説明するための超音波走査線(超音波ビーム)を示す図である。
【図2】コントラストエコー画像を説明するための図である。
【図3】本発明の等間隔スライディングウィンドウ法を説明するための心臓(心筋と心腔)の断面図である(深部方向に直線で等分割した例)。
【図4】本発明の等間隔スライディングウィンドウ法での色づけ方法を説明するための心臓(心筋と心腔)の断面図である。
【図5】本発明の等間隔スライディングウィンドウ法を説明するための心臓(心筋と心腔)の断面図である(深部方向に円弧で等分割した例)。
【図6】本発明の重複スライディングウィンドウ法での色づけ方法を説明するための心臓(心筋と心腔)の断面図である。
【図7】本発明のスライディングウィンドウ法をより詳細に色分けするために、心筋を細分化して計算領域を設ける例を説明するための心臓(心筋と心腔)の断面図である。
【図8】本発明のピクセルバイピクセル法を説明するための心臓(心筋と心腔)の断面図である。
【図9】本発明の色づけ法における、計算領域から色の滲み出し部分を排除する方法を説明するための心臓(心筋と心腔)の断面図である。
【図10】本発明の色づけ法によって色づけした心臓(心筋と心腔)の画像を説明するための画面を表す図である。
【符号の説明】
1 心臓
110 心筋
111 心内膜(心内膜線)
112 心外膜(心外膜線)
120 心腔
2 超音波プローブ
210 先端部(心尖部、超音波発信部の体表に接した部分)
Ak 計算領域kの輝度(心筋部計算領域)
Bk 計算領域kの輝度(心腔部計算領域)
b1 バッファ線(心腔側)
b2 バッファ線(心筋側)
c カラーテーブル
d1 音場線間隔
d2 心内膜と内線の間隔、心内膜と外線の間隔
d3 重複スライディングウィンドウ法における計算領域のずらす間隔
d4 心内膜とバッファ線の間隔
l0 音場線(等音場線)
l1 内線
l2 外線
Claims (6)
- 気泡性超音波造影剤を血管に注入して超音波診断装置によって得られた心筋と心腔のコントラストエコー画像を評価する血液量解析・表示方法において、
計算領域を細分化して心筋、心腔の音場が同等と思われる部位で、心筋の輝度と心腔の輝度の差分を求め、前記差分位で心筋の血液量を求め、その血液量に対してカラー表示を行う手段、
を備えたことを特徴とするマイオカーディアルブラッドボリュームマップによる血液量解析・表示方法。 - 請求項1に記載のマイオカーディアルブラッドボリュームマップによる血液量解析・表示方法において、
コントラストエコー画像より得られる関心領域に対しての、心筋の輝度をI1、
心腔の輝度をI2、および定数をKとしたとき、アンプリチュードAを、
Klog10A=I1−I2
なる式から求め、関心領域における心筋の血液量V(単位ml/100g)を
V=100×A(ml/100g)
と求める手段、
を備えたことを特徴とする請求項1に記載のマイオカーディアルブラッドボリュームマップによる血液量解析・表示方法。 - 請求項1に記載のマイオカーディアルブラッドボリュームマップによる血液量解析・表示方法で使用するコントラストエコー画像において、(1)心筋と心腔との境界である心内膜、および心筋の外側の膜(心外膜)をそれぞれトレースし、心内膜の前記トレース線(心内膜線)に対してd2の間隔で内側および外側に心内膜線に沿って内線と外線を描画し、さらに深部方向に対して音場が同じと思われる境界線(音場線)をd1の等間隔に描画して、心内膜線と内線および音場線で囲まれる領域を心腔部計算領域とし、心内膜線と外線および境界線で囲まれる領域を心筋部計算領域とする手段、
(2)前記の計算領域に対して、心筋部計算領域の平均輝度から心腔部計算領域の平均輝度を差し引いて差分値を求め、前記の差分値に相当する色をカラーテーブルから選択して境界線で区切られる心筋全体を色づけする手段、
(3)音場線をd1ずつずらして前記の差分計算を繰り返し、全心筋を色づけする手段、
を備えたことを特徴とする等間隔スライディングウィンドウ法。 - 請求項3に記載のスライディングウィンドウ法において、
より音場領域の音場差(輝度差)の影響を縮小するために、計算領域を微小に深部方向(例えば微小移動間隔d3(=d1/10))にずらしながら心筋部計算領域の輝度と心腔部計算領域の輝度との差を計算し、計算した差分値で心筋を色づけする手段、
を備えたことを特徴とする重複スライディングウィンドウ法。 - 請求項3に記載のスライディングウィンドウ法において、
心内膜線と心外膜線および音場線で囲まれる領域の各ピクセル(画素)ごとに輝度を求め、この輝度と心腔部計算領域の平均輝度を引いて輝度の差分を求め、心筋内の色づけを前記差分値で画素単位で行う手段、
を備えたことを特徴とするピクセルバイピクセル法。 - 請求項3〜5に記載の心筋の色づけ方法において、
パワーモードのようなカラー画像では心腔の色が心内膜から心筋に滲み出す場合があり、この滲み出しによる影響を排除するために、心内膜線の内側と外側に一定の間隔を置いてバッファ線を引き、バッファ線と外線、バッファ線と内線とに囲まれる部分を計算領域とする手段、
を備えたことを特徴とする心筋色づけ方法。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002379760A JP3683886B2 (ja) | 2002-12-27 | 2002-12-27 | マイオカーディアルブラッドボリュームマップによる血液量解析・表示方法 |
US10/742,836 US7285095B2 (en) | 2002-12-27 | 2003-12-23 | Method of analyzing and displaying blood volume using myocardial blood volume map |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002379760A JP3683886B2 (ja) | 2002-12-27 | 2002-12-27 | マイオカーディアルブラッドボリュームマップによる血液量解析・表示方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2004208807A true JP2004208807A (ja) | 2004-07-29 |
JP3683886B2 JP3683886B2 (ja) | 2005-08-17 |
Family
ID=32708407
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002379760A Expired - Fee Related JP3683886B2 (ja) | 2002-12-27 | 2002-12-27 | マイオカーディアルブラッドボリュームマップによる血液量解析・表示方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US7285095B2 (ja) |
JP (1) | JP3683886B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011093193A1 (ja) * | 2010-01-29 | 2011-08-04 | 株式会社 日立メディコ | 超音波診断装置およびその計測点追跡方法 |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7367947B2 (en) * | 2002-11-08 | 2008-05-06 | Koninklijke Philips Electronics N.V. | Elevation beamwidth control for control for contrast imaging |
US7951083B2 (en) * | 2004-02-05 | 2011-05-31 | Siemens Medical Solutions Usa, Inc. | Motion analysis improvements for medical diagnostic ultrasound |
US7520860B2 (en) | 2005-04-13 | 2009-04-21 | Marie G. Johnson | Detection of coronary artery disease using an electronic stethoscope |
US20090234231A1 (en) * | 2008-03-13 | 2009-09-17 | Knight Jon M | Imaging Catheter With Integrated Contrast Agent Injector |
KR101138613B1 (ko) * | 2009-04-30 | 2012-04-26 | 한림대학교 산학협력단 | 컬러 도플러 모드 영상을 처리하는 초음파 시스템 및 방법 |
JP5459832B2 (ja) * | 2009-06-02 | 2014-04-02 | 東芝メディカルシステムズ株式会社 | 超音波診断装置 |
DE102010043849B3 (de) | 2010-11-12 | 2012-02-16 | Siemens Aktiengesellschaft | Vorrichtung und Computertomograph zur Bestimmung und Darstellung der Durchblutung des Herzmuskels |
KR102582540B1 (ko) * | 2015-06-16 | 2023-09-25 | 삼성메디슨 주식회사 | 초음파 장치 및 그 동작 방법 |
AU2017358321B2 (en) * | 2016-11-11 | 2022-01-06 | MedTrace Pharma A/S | Method and system for modelling a human heart and atria |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5425366A (en) * | 1988-02-05 | 1995-06-20 | Schering Aktiengesellschaft | Ultrasonic contrast agents for color Doppler imaging |
EP0758251A1 (en) * | 1994-05-03 | 1997-02-19 | Molecular Biosystems, Inc. | Composition for ultrasonically quantitating myocardial perfusion |
US5743266A (en) * | 1995-04-25 | 1998-04-28 | Molecular Biosystems, Inc. | Method for processing real-time contrast enhanced ultrasonic images |
AU6673498A (en) * | 1997-03-10 | 1998-09-29 | Robin Medical Inc. | Method and apparatus for the assessment and display of variability in mechanicalactivity of the heart, and enhancement of ultrasound contrast imaging by variab ility analysis |
GB9708246D0 (en) * | 1997-04-24 | 1997-06-18 | Nycomed Imaging As | Improvements in or relating to ultrasound imaging |
US6708055B2 (en) * | 1998-08-25 | 2004-03-16 | University Of Florida | Method for automated analysis of apical four-chamber images of the heart |
US6015384A (en) * | 1998-08-31 | 2000-01-18 | Acuson Corporation | Ultrasonic system and method for tissue viability imaging |
US6217520B1 (en) * | 1998-12-02 | 2001-04-17 | Acuson Corporation | Diagnostic medical ultrasound system and method for object of interest extraction |
US6080107A (en) * | 1999-01-26 | 2000-06-27 | Hewlett-Packard Company | Methods for the use of contrast agents in ultrasonic imaging |
JP2001212144A (ja) * | 2000-01-31 | 2001-08-07 | Toshiba Corp | 超音波診断装置及び超音波画像化方法 |
NO20010234D0 (no) * | 2001-01-12 | 2001-01-12 | Nycomed Imaging As | Perfusjon avbildningsbilde |
US6592522B2 (en) * | 2001-06-12 | 2003-07-15 | Ge Medical Systems Global Technology Company, Llc | Ultrasound display of displacement |
US6994673B2 (en) * | 2003-01-16 | 2006-02-07 | Ge Ultrasound Israel, Ltd | Method and apparatus for quantitative myocardial assessment |
-
2002
- 2002-12-27 JP JP2002379760A patent/JP3683886B2/ja not_active Expired - Fee Related
-
2003
- 2003-12-23 US US10/742,836 patent/US7285095B2/en not_active Expired - Fee Related
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011093193A1 (ja) * | 2010-01-29 | 2011-08-04 | 株式会社 日立メディコ | 超音波診断装置およびその計測点追跡方法 |
CN102724918A (zh) * | 2010-01-29 | 2012-10-10 | 株式会社日立医疗器械 | 超声波诊断装置及其测量点追踪方法 |
JP5753798B2 (ja) * | 2010-01-29 | 2015-07-22 | 株式会社日立メディコ | 超音波診断装置およびその作動方法 |
Also Published As
Publication number | Publication date |
---|---|
US7285095B2 (en) | 2007-10-23 |
US20040138567A1 (en) | 2004-07-15 |
JP3683886B2 (ja) | 2005-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11801033B2 (en) | Medical diagnostic apparatus and medical analysis method | |
US8460192B2 (en) | Ultrasound imaging apparatus, medical image processing apparatus, display apparatus, and display method | |
CN103889337B (zh) | 超声波诊断装置以及超声波诊断装置控制方法 | |
JP4473543B2 (ja) | 超音波診断装置 | |
US10743845B2 (en) | Ultrasound diagnostic apparatus and method for distinguishing a low signal/noise area in an ultrasound image | |
US6641538B2 (en) | Ultrasonic diagnostic apparatus and method of controlling a ultrasonic diagnostic apparatus | |
JP4763883B2 (ja) | 超音波診断装置 | |
JP5868067B2 (ja) | 医用画像診断装置、画像処理装置及び方法 | |
JP5680654B2 (ja) | 超音波診断装置及び超音波画像表示方法 | |
WO2014077396A1 (ja) | 超音波診断装置及び画像処理方法 | |
JP3683945B2 (ja) | 超音波診断装置 | |
JP2011067546A (ja) | 超音波診断装置及び弾性指標算出方法 | |
JP5689591B2 (ja) | 超音波診断装置及び超音波画像処理プログラム | |
JP5389118B2 (ja) | 超音波診断装置 | |
JP3683886B2 (ja) | マイオカーディアルブラッドボリュームマップによる血液量解析・表示方法 | |
JP4253494B2 (ja) | 超音波診断装置 | |
JP5300171B2 (ja) | 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム | |
JP4113702B2 (ja) | 超音波診断装置及び超音波信号解析装置 | |
JP4764520B2 (ja) | 超音波診断装置 | |
JP4764521B2 (ja) | 超音波診断装置 | |
JP2017080040A (ja) | 超音波診断装置およびプログラム | |
JP4497611B2 (ja) | 超音波診断装置 | |
Robert et al. | Spectral Analysis of Tissue Displacement for Cardiac Activation Mapping: Ex Vivo Working Heart and In Vivo Study | |
US20220330921A1 (en) | Medical image diagnosis apparatus, medical information processing apparatus, and medical image processing method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050124 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050201 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050404 |
|
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: 20050510 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050526 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110603 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110603 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120603 Year of fee payment: 7 |
|
LAPS | Cancellation because of no payment of annual fees |