JP3697292B2 - 超音波診断装置 - Google Patents
超音波診断装置 Download PDFInfo
- Publication number
- JP3697292B2 JP3697292B2 JP11875895A JP11875895A JP3697292B2 JP 3697292 B2 JP3697292 B2 JP 3697292B2 JP 11875895 A JP11875895 A JP 11875895A JP 11875895 A JP11875895 A JP 11875895A JP 3697292 B2 JP3697292 B2 JP 3697292B2
- Authority
- JP
- Japan
- Prior art keywords
- parameter
- luminance
- time
- image data
- roi
- 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
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Description
【０００１】
【産業上の利用分野】
本発明は超音波診断装置に係り、とくに超音波造影剤を被検体に投与して、この超音波造影剤による超音波信号の散乱を含むエコー信号を受信し、このエコー信号に基づいて診断対象の機能情報を得るようにした超音波診断装置に係る。
【０００２】
【従来の技術】
近年、超音波造影剤を用いたコントラストエコー法が心筋画像の解析分野で注目されている。
【０００３】
このコントラストエコー法の一つとして、超音波造影剤を動脈から注入する動脈注入による心筋コントラストエコー法が研究されており、心筋分布像（perfusion)による心筋内血流の灌流域の評価に利用されている。この心筋コントラストエコー法は大動脈に留置されたカテーテルより超音波造影剤（例えば、用手的あるいはソニケータにより気泡の生成された５％ヒトアルブミン）を注入するものである。造影剤により心筋内血流の灌流域は、Ｂモード上の輝度増強領域として表示される。同様に、血流の灌流域の評価あるいは腫瘍の支配血管系を評価するために、腹部領域でも動脈注入によるコントラストエコー法が研究されている。これらのコントラストエコー法を実施する診断装置は、一般検査用の超音波診断装置あるいはさらにワークステーションが用いられる。これにより、目視によりＢモード像の輝度増強を評価したり、あるいは、メモリに記憶された画像データをワークステーション上で適当な処理をしたりした後、輝度レベルの変化を定量評価するようになっている。
【０００４】
また、近年、超音波造影剤を静脈から注入して左心系の評価が可能な超音波造影剤が開発され、これを用いた超音波コントラストエコー法が試みられている。
【０００５】
この超音波造影剤としては、塩野義製薬株式会社により輸入販売されている、『５％人血清アルブミンを超音波処理するときに生成するアルブミン膜の中に、空気を封じ込めた平均粒子径約４μｍの空気小球体』（販売名：アルブネックス注５ｍｌ）がある。
【０００６】
この静脈注入によるコントラストエコー法は、現在、試験，研究段階であり、今後、頭部・心腔・腹部などの診断で、その有用性に期待が高まっている。
【０００７】
【発明が解決しようとする課題】
このように近年、超音波造影剤を用いた診断法が盛んになってきているが、これらの診断法にあっては、超音波造影剤の超音波反射信号に対する増強作用にのみ着目したものに止まっており、造影剤の特異な性質に着目した計測は未だ行われていなかった。
【０００８】
本発明はこのような状況に鑑みてなされたもので、気体から成る泡を成分とする超音波造影剤の体積変化に因って超音波エコー信号の散乱強度が変化することを利用して、生体内の機能的情報を的確且つ簡便に得ることができる超音波診断装置を提供することを、その目的とする。
【０００９】
【課題を解決するための手段】
前記目的を達成するため、本発明の超音波診断装置は、気体から成る泡を含んだ超音波造影剤を注入した生体に超音波信号を照射する送信手段と、前記超音波信号の生体からの反射信号に基づいて断層面の画像データを得る受信処理手段と、前記画像データを表示する断層像表示手段と、前記断層像表示手段によって表示された断層像上に所望のＲＯＩを設定するＲＯＩ設定手段と、前記断層面の経時的に収集される複数枚の画像データに対して前記ＲＯＩ内の領域に含まれる画素の輝度を画像毎に演算する輝度演算手段と、この輝度演算手段により演算された時系列の輝度のデータから、前記生体の心拍由来の変化成分を表すパラメータを演算するパラメータ演算手段と、このパラメータ演算手段の演算結果を表示するパラメータ表示手段とを備え、前記パラメータ演算手段は、前記心拍由来の変化成分自体を抽出するハイパスフィルタ又はバンドパスフィルタを備たことを要部とする。
【００１０】
【作用】
生体には気体から成る泡を含んだ超音波造影剤が注入される。この状態で生体に超音波信号を照射し、その反射信号に基づいて断層面の画像データが得られ、断層像が表示される。操作者は、表示された断層像上に所望のＲＯＩを設定する。断層面の経時的に収集される複数枚の画像データに対して、設定したＲＯＩ内の領域に対応する画素の輝度が例えば平均値として画像毎に演算される。この時系列の輝度データに基づいて生体の心拍由来の変化成分を表すパラメータ（例えば、心拍に呼応して超音波造影剤の体積が変化することに伴う心拍由来の変化成分）が演算され、この演算結果が数値やグラフとして表示される。これにより、超音波造影剤の輝度増強作用による診断情報に止まらず、造影剤が気体の泡を成分とすることに伴う診断情報も得ることができる。
【００１１】
【実施例】
本発明の一実施例を図１〜図９に基づいて説明する。なお、この実施例は心機能のパラメータ演算に関して複数の実施態様を含んでいる。
【００１２】
図１に示す超音波診断装置は、フェーズド・アレイ・タイプの超音波プローブ１０（以下、「プローブ」という）と、このプローブ１０を介して被検体との間で超音波信号を送受する送受信部１１と、この送受信部１１が得た超音波エコー信号に基づいてデータ処理を行うデータ処理部１２と、このデータ処理部１２に接続された入力器１３及び表示モニタ１４とを備えている。送受信部１１及びデータ処理部１２のステアリングライン制御、動作タイミング制御などは制御回路１５によって指令される。
【００１３】
プローブ１０は超音波信号の送受を兼ねるもので、複数の振動子（圧電体）を一方向に配列させた超音波振動子群１０Ａを備えている。
【００１４】
送受信部１１は、その送信系として基準発振器２０、送信遅延回路群２１、及びパルサ回路群２２を備えており、パルサ回路群２２の出力側が超音波振動子群１０Ａにチャンネル毎に接続されている。これにより、各振動子はパルサ回路群２２から供給された駆動電圧パルスを受けて励振し、被検体に向けて超音波ビーム信号を照射する。なお、送信フォーカスを掛けない場合、送信遅延回路群２１を省くこともできる。
【００１５】
さらに送受信部１１は、受信系として、プリアンプ群２３、受信遅延回路群２４、加算回路２５、対数圧縮回路２６、検波回路２７、及びＡＤＣ（Ａ／Ｄ変換器）２８を備えており、プリアンプ群２３の入力端が超音波振動子群１０Ａにチャンネル毎に接続されている。生体内から反射してきた超音波エコー信号は、プローブ１０で受信され、振動子群１０Ａにより電気量のエコー信号に変換される。このエコー信号はプリアンプ群２３でチャンネル毎に増幅された後、受信遅延回路群２４でチャンネル毎に遅延が掛けられ、加算回路２５により加算される。
この遅延及び加算（整相加算）により受信フォーカスが掛けられる。この遅延・加算信号はその後、対数圧縮回路２７で圧縮され、検波回路２８で包絡線検波された後、ＡＤＣ２８により超音波ビームラインに沿ったデジタル量の画像データが得られる。超音波ビームを走査することで、２次元の画像データが得られる。
【００１６】
データ処理部１２はＡＤＣ２８から出力された画像データを記憶する画像メモリ（フレームメモリ）３０と、その画像データに後述するグラフデータなどを合成する合成回路３１と、この合成回路３１の出力データをアナログ量に変換するＤＡＣ（Ｄ／Ａ変換器）３２とを備え、このＤＡＣ３２の出力側が表示モニタ１４に接続されている。これにより、画像メモリ３０に一時記憶された画像データは表示モニタ１４にリアルタイムに表示される。
【００１７】
またデータ処理部１２は、ＲＯＩ設定回路３３を備えている。このＲＯＩ設定回路３３はトラックボールなどの入力器１３から出力されるＲＯＩ設定情報（ＲＯＩ位置、大きさ、形状）を受けて、これに対応するＲＯＩ設定信号（ＲＯＩのグラフィックデータを含む）を合成回路３１に出力する。このため、観察者が表示モニタ１４の画面を観察しながら入力器１３を操作することにより、例えば図２に示す如く、画面上の所望位置に所望形状及び大きさのＲＯＩ（関心領域）を設定することができる。
【００１８】
データ処理部１２はさらに、図１に示すように輝度の平均強度演算回路３５、メモリ３６、パラメータ演算回路３７、及びグラフ表示回路３８を備えている。
これらの構成は本発明の要旨を具現化するもので、以下のように構成され機能する。
【００１９】
平均強度演算回路３５はＲＯＩ設定回路３４から出力されるＲＯＩ設定信号を受けてＲＯＩの設定内容（位置、形状、大きさ）を解読するとともに、そのＲＯＩ内の領域を形成する複数の画素の各輝度値（画素値）を画像メモリ３０から読み出す。そして、平均強度演算回路３５は各輝度の平均値（又は積算値でも良い）を算出する。この算出は経時的に（フレーム毎に）連続して得られる画像データに対して行われ、各フレームの算出毎に、その算出値がグラフ表示回路３８及びメモリ３６に送られるようになっている。
【００２０】
なお、この平均強度演算回路３５の演算においては、上述の如く、対数圧縮処理がなされた画像輝度データ（すなわち超音波の散乱強度）のまま近似的に輝度平均値又は積算値を算出するとしてもよいが、以下のように補正しながら算出しても良い。この補正は
【数１】
ｙ＝ｅｘｐ（αｘ）
ここで、
α：定数 ｌｏｇ（１０）／１０（ｌｏｇは自然対数を表す）、
ｙ：散乱強度、
ｘ：画像輝度
の式を用いて対数圧縮データを元に戻す逆変換処理を行うものである。従って、この補正処理を伴う場合、演算回路３５はそのデータ（散乱強度）についてＲＯＩ内の平均値又は積算値を算出することになる。
【００２１】
なお、このように平均強度演算回路３５を作動させている間、その演算結果に誤差を生じさせる恐れのある処理（例えば画像の輪郭等を強調表示するための信号処理）は、制御回路１５によって禁止される。
【００２２】
グラフ表示回路３８は、平均強度演算回路３５から送られてくる、経時的に連続したＲＯＩ内輝度の平均強度データを入力するとともに横軸を時間、縦軸を平均強度Ｂavとするグラフデータをリアルタイムに作成し、このグラフデータを合成回路３１に出力するようになっている。合成回路３１は、リアルタイムに供給されるグラフデータとＲＯＩを設定した超音波断層像データとを例えば分割表示態様の画像データに合成する。このため、表示モニタ１４には例えば図３に示す如く、平均強度Ｂavのグラフを左側、測定している断層像を右側に置いた表示画面が現れる。平均強度Ｂavのグラフにおける時間軸のスパンは数十秒から数分のオーダーであり、造影剤注入からの時間経過とともに、このグラフは時々刻々と更新しながら表示される。なお、この平均強度グラフは超音波断層像に対して上下左右何れの位置に分割表示してもよいし、単独表示するようにしてもよい。
【００２３】
本発明で得られる平均強度曲線は、超音波造影剤が気体からなる泡でできていることに因り、通常のＸ線診断の輝度変化曲線に加えて、心拍由来の情報をも含んでいる。
【００２４】
これを詳述すると、以下のようである。超音波信号の散乱強度は造影剤の体積濃度により変化する。注入直後の濃度は高いが、時間経過とともにその濃度は低下するので、この濃度変化に対応して散乱強度も変化する。例えば、図３中の平均強度Ｂavの曲線において、点線で示した部分が造影剤自体の濃度に由来しており、時間ｔo で造影剤を注入した場合、ＲＯＩ内輝度の平均強度は時間と伴に徐々に上がって最大値に到達し、その後なだらかに低下する。この曲線部分（以下、「造影剤濃度の変化成分ＢCT」という）は、Ｘ線診断装置などで通常得られる輝度変化曲線（Time Density Curve:TDC) と同等なものである。
【００２５】
さらに、本発明の超音波造影剤は泡を成分としているため、圧力変化により個々の泡の体積が変化する。つまり、心臓などのように圧力が小刻みにかつ大きく変化する診断対象の場合、この圧力変化に伴って小刻みな体積変化が合わせて生じ、従って平均強度Ｂavも小刻みに変動する。診断対象が心臓の場合、かかる小刻みな平均強度の変動は心拍に由来する成分（以下、「心拍由来の変化成分ＢHB」という）となり、図３に示す如く、ＴＤＣに相当する造影剤濃度の変化成分ＢCT（点線の曲線）に重畳されている（実線の曲線）。この両者の散乱強度（曲線ＢCTとＢHB）の変化には、時間変化の速度に大きな差がある。造影剤濃度の変化成分ＢCTは時間的にゆっくり変化し、数十秒〜数分持続するのに対して、心拍由来の変化成分ＢHBは心拍に同期するため、その周期は１秒以下である。
【００２６】
上記造影剤濃度の変化成分及び心拍由来の変化成分を個別に抽出したり識別したりするために、前述したメモリ３６及びパラメータ演算回路３７が設けられている。メモリ３６には、平均強度演算回路３５から出力された輝度の平均強度の時系列データＢavが記憶される。
【００２７】
パラメータ演算回路３７は、メモリ３６に記憶された平均強度の時系列データを読み出して心機能に関する指令パラメータを演算する。時系列データの読出しタイミングは、入力器１３を介して演算が指令された時間でもよいし、予め決めた適宜な時間（例えば造影剤注入後の一定時間経過後）でもよい。また、演算するパラメータも入力器１３からその都度指令するようにしてもよいし、予め決めておいてもよい。
【００２８】
まず、心拍由来の変化成分そのものを抽出する場合、図４に示す如く、パラメータ演算回路３７はデジタル型のフィルタ３７ａを有する。フィルタ３７ａがハイパスフィルタの場合、そのカットオフ周波数はほぼ心拍周期の周波数未満のある値に設定される。フィルタ３７ａがバンドパスフィルタの場合、その通過帯域は心拍の周波数成分に設定される。
【００２９】
造影剤濃度の変化成分のみを抽出する場合、パラメータ演算回路３７は図５に示す如く、デジタル型のバンドパスフィルタ３７ｂと減算器３７ｃとを有する。
バンドパスフィルタ３７ｂは心拍由来の変化成分ＢHBを抽出するもので、減算器３７ｃにより、その成分ＢHBを元の時系列のデータＢavから減算することで造影剤濃度の変化成分ＢCTのみが抽出される。
【００３０】
なお、造影剤濃度の変化成分ＢCTのみを抽出するには、上記の減算する構成のほか、単にローパスフィルタを用いて心拍由来の変化成分ＢHBを平均強度Ｂavから抽出するようにしてもよい。
【００３１】
このように抽出された心拍由来の変化成分ＢHB及び造影剤濃度の変化成分ＢCTはグラフ表示回路３８に送られる。グラフ表示回路３８は前述した平均強度Ｂavのグラフデータ作成のほか、今回抽出した変化成分ＢHB、ＢCTのグラフデータ（横軸は時間軸）を作成して合成回路３１に送る。このため、表示モニタ１４には断層像などと伴に、それらの変化成分のグラフが表示され、心機能の状態を表すパラメータとして目視観察することができる。
【００３２】
さらに、パラメータ演算回路３７は心圧の最大値と最小値の比を演算することもできる。散乱体からの信号の散乱強度はその体積濃度に比例することが知られている。また気体の体積と圧力の関係はボイルシャルルの法則から互いに反比例する。このことから心拍周期内の時間における散乱強度の最大値と最小値の比をとればその間の心圧の最大値と最小値の比が得られることになる。
【００３３】
この心圧の最大値と最小値の比を演算する場合、パラメータ演算回路３７はＣＰＵ機能を備え、図６に示す処理を心周期ごとに行う。すなわち、メモリ３６に記憶されている各フレーム（時系列）のＲＯＩ内の輝度の平均強度データＢavの中から、１心周期分のデータを読み込む（同図ステップＳ１）。次いで、読み込んだ１心周期分の時系列データの中の最大値Ｂmax 及び最小値Ｂmin を演算する（同図ステップＳ２、Ｓ３）。さらに、最大値Ｂmax 及び最小値Ｂmin の比「Ｂmax ／Ｂmin 」を演算し、その演算値をグラフ表示回路３８に出力する（同図ステップＳ４）。
【００３４】
グラフ表示回路３８は上記輝度比「Ｂmax ／Ｂmin 」のデータを入力して、内部メモリに記憶していた平均強度曲線Ｂavのデータに、例えば図示しない計時信号を用いて心時相を合わせたグラフ（図７参照）のデータを作成する。作成されたグラフデータは合成回路３１に送られ、断層像などの他の必要な画像データと一緒にフレーム画像に合成される。
【００３５】
この結果、表示モニタ１４には例えば図７に示す如く、横軸を時間とし、下側に輝度の平均強度曲線Ｂav、これに心時相を合わせた状態で上側に輝度比「Ｂmax ／Ｂmin 」のグラフを配置した画像が得られる。これにより、心臓のポンプ機能の状態に関する一つの指標が好適に得られ、モニタで観察することができる。したがって、心臓のポンプ機能を判断する上で、従来のように、カテーテルを使って圧力センサを左室又は右室に留置させて心圧を測定するという困難な手法を用いなくても済む。
【００３６】
なお、上記の輝度比を演算する際、対数圧縮を掛けた数値そのものを用いる場合、輝度差「Ｂ ｍａｘ −Ｂ ｍｉｎ 」を演算することで輝度比と等価な解析を行うことになる。
【００３７】
さらに、超音波造影剤の散乱強度（エコー信号）の変化は前述したように２種の変化成分ＢCT及びＢHBの和として得られる。この内、造影剤濃度の変化成分ＢCTの時間変化は非常に緩やかであるから、その時間微分値は小さく無視できる。従って、輝度の平均強度データＢavに適宜な微小時間幅ｄｔ（例えば２秒以下）の時間微分演算を施すと、その結果は心拍由来の変化成分ＢCTの微分値にほぼ一致し、心圧の時間変化（ｄｐ／ｄｔ）が近似的に得られる。
【００３８】
この心圧の時間変化（ｄｐ／ｄｔ）を計測する場合、パラメータ演算回路３７は例えば図８に示す如く、デジタル型の微分回路３７ｄを備えている。この微分回路３７ｄで演算された心圧の時間変化（ｄｐ／ｄｔ）は同様にグラフ表示回路３８に送られ、内臓メモリに記憶されていた輝度の平均強度データＢavに心時相（時間）を合わせた状態で合成されたグラフデータに変換される。これにより、表示モニタ１４には、計測対象の断層像などと伴に、例えば図９に示すグラフが表示され、輝度の平均強度Ｂavの変化に対応する心圧の時間変化（ｄｐ／ｄｔ）を目視で観察することができる。これによって、心臓のポンプ機能に対する良い指標が得られる。
【００３９】
なお、上記の如く心圧の時間変化（ｄｐ／ｄｔ）を計測するにあたり、造影剤濃度の変化成分の影響を殆ど完全に低減するには、平均強度Ｂavのデータにハイパスフィルタやバンドパスフィルタを掛けて心拍由来の変化成分ＢHBのみを抽出し、この抽出データを時間微分するパラメータ演算回路の構成にしてもよい。
【００４０】
本発明に係る超音波診断装置の別の実施例を図１０及び図１１に基づき説明する。この実施例は、前述した心圧の時間変化（ｄｐ／ｄｔ）を計測する場合の、表示状態の改善に関する。図１と同一又は同等の構成要素には同一符号を用いて説明を省略又は簡単化する。
【００４１】
図１０に示す超音波診断装置は、図１記載の構成に加えて、ＥＣＧ（心電計）４０、ＥＣＧ用アンプ４１、ＡＤＣ４２及び表示制御回路４３を備える。ＥＣＧ４０により収集された被検体の心電図信号はＥＣＧ用アンプ４１、ＡＤＣ４２を介して画像メモリ３０に記憶される。従って、プローブ１０を介してエコー信号を収集すると同時に心電図信号も収集され、画像メモリ３０には画像データと心電図データが蓄えられる。
【００４２】
この画像データ及び心電図データは合成回路３１に送られる。これと伴に、合成回路３１にはグラフ表示回路３８から前述した如く輝度の平均強度Ｂavと心圧の時間変化（ｄｐ／ｄｔ）のグラフデータ（図９参照）が送られてくる。表示制御回路４３は入力器１３からオペレータの指令によって与えられる拡大／縮小表示などの表示指令信号を受けて、合成回路３１に、対応する表示制御信号を送るようになっている。これにより、合成回路３１は画像データやグラフデータを指令された適宜な態様で合成するとともに、拡大／縮小処理も行うようになっている。
【００４３】
この結果、オペレータが入力器１３の操作により拡大表示を指令すると、例えば図１１に示す如く、拡大された平均強度Ｂav、心圧の時間変化（ｄｐ／ｄｔ）、及び心電図のグラフ像が同時・拡大表示される。したがって、表示画面から時相評価を行い易くなるとともに、心圧変化の数値的評価のみならず、時間的変化をも容易に評価可能になる。
【００４４】
本発明のさらに別の実施例を図１２及び図１３に基づき説明する。この実施例は前述した心圧の時間変化（ｄｐ／ｄｔ）を２次元画像のまま処理／表示するようにしたものである。
【００４５】
図１２に示す超音波診断装置は、そのデータ処理部１２において、図１に示した平均強度演算回路３５、メモリ３６、パラメータ演算回路３７、及びグラフ表示回路３８に代えて、画像データ出力回路４９、空間フィルタ５０、メモリ５１、及び画像間演算回路５２を備える。
【００４６】
画像データ出力回路４９は、指令されたＲＯＩの内部領域を形成する画素の２次元画像データＢ2Dを画像メモリ３０から読み出して空間フィルタ５１に出力するようになっている。この２次元画像データＢ2Dは、空間フィルタ５０により、２次元画像のまま、スペックルノイズ除去のフィルタリングが施され、スペックルノイズの影響が著しく低減されたＲＯＩ内画像データが生成される。このフィルタリングされた画像データＢ2Dはメモリ５１に送られ、一時記憶される。画像間演算回路５２はメモリ５１に記憶された、連続する２枚の画像間の画素同士の輝度差又は強度（輝度）比を演算し、その演算値を新たな輝度値とする、ＲＯＩ内の２次元画像データＢ2D ＊を生成するようになっている。この２次元画像データＢ2D ＊の各画素は、各画素値の時間微分値であり、前述した心圧の時間変化（ｄｐ／ｄｔ）の空間分布の画素を形成する。
【００４７】
このＲＯＩ内領域に相当する２次元画像データＢ2D ＊は合成回路３１、ＤＡＣ３２を介して表示モニタ１４に、例えば図１３に示す如く表示される。これにより、心圧の時間変化（ｄｐ／ｄｔ）の空間分布を実時間で観測することができる。
【００４８】
さらに本発明の他の実施例を図１４に基づき説明する。
【００４９】
図１４に示す超音波診断装置は、データ処理部１２を外部機器として備え、送受信部１１からオンラインで画像データを送信可能にしたものである。このために、送受信部１１の出力側に画像メモリ６０を備えて、収集した画像データを一度蓄積可能になっている。この画像メモリ６０の読み出し側はインターフェース６１、６２を介してデータ処理部１２（外部機器）に接続され、制御回路１５からの指令によって画像データをオンライン伝送させる。
【００５０】
これにより、データ処理部１２の画像データ処理をバッチ処理で行うことができる一方、データ処理部１２での各コンポーネントを画像解析専用の装置に組み込むこともできる。なお、送受信部１１とデータ処理部１２間の画像データ伝送はオフライン伝送であってもよい。
【００５１】
【発明の効果】
以上説明したように、本発明の超音波診断装置によれば、断層面の経時的に収集される複数枚の画像データに対して、設定したＲＯＩ内の領域に対応する画素の輝度を画像毎に演算し、この時系列の輝度データに基づいて生体の心拍由来の変化成分を表すパラメータ（例えば、心拍に呼応して超音波造影剤の体積が変化することに伴う心拍由来の変化成分）を演算し、この演算結果を数値やグラフとして表示するようにしたので、超音波造影剤の輝度増強作用による診断情報に止まらず、造影剤が気体の泡を成分とすることに伴う心圧やその最大値、最小値の比など、生体内の機能的情報を的確且つ簡便に得ることができる。
【図面の簡単な説明】
【図１】本発明の実施例に係る超音波診断装置のブロック図。
【図２】断層像上に設定するＲＯＩの説明図。
【図３】ＲＯＩ内の輝度の平均強度グラフの表示例を示す図。
【図４】パラメータ演算回路の位置態様を示すブロック図。
【図５】パラメータ演算回路の別の態様を示すブロック図。
【図６】パラメータ演算回路の更に別の態様に係る処理を示す概略フローチャート。
【図７】図６に示す処理に係る画像例を示す図。
【図８】パラメータ演算回路の更に別の態様に係るブロック図。
【図９】図８に示す構成に係る画像例を示す図。
【図１０】本発明の他の実施例に係る超音波診断装置のブロック図。
【図１１】図１０の構成に係る装置の画像例を示す図。
【図１２】本発明の更に他の実施例を示す超音波診断装置のブロック図。
【図１３】図１２の構成に係る装置の画像例を示す図。
【図１４】本発明の更に他の実施例を示す超音波診断装置のブロック図。
【符号の説明】
１０ プローブ
１１ 送受信部
１２ データ処理部
１３ 入力器
１４ 表示モニタ
１５ 制御回路
３０ 画像メモリ
３１ 合成回路
３２ ＤＡＣ
３４ ＲＯＩ設定回路
３５ 平均強度演算回路
３６ メモリ
３７ パラメータ演算回路
３８ グラフ表示回路
３７ａ，３７ｂ フィルタ
３７ｃ 減算器
３７ｄ 微分回路
４９ 画像データ出力回路
５０ 空間フィルタ
５１ メモリ
５２ 画像間演算回路
【産業上の利用分野】
本発明は超音波診断装置に係り、とくに超音波造影剤を被検体に投与して、この超音波造影剤による超音波信号の散乱を含むエコー信号を受信し、このエコー信号に基づいて診断対象の機能情報を得るようにした超音波診断装置に係る。
【０００２】
【従来の技術】
近年、超音波造影剤を用いたコントラストエコー法が心筋画像の解析分野で注目されている。
【０００３】
このコントラストエコー法の一つとして、超音波造影剤を動脈から注入する動脈注入による心筋コントラストエコー法が研究されており、心筋分布像（perfusion)による心筋内血流の灌流域の評価に利用されている。この心筋コントラストエコー法は大動脈に留置されたカテーテルより超音波造影剤（例えば、用手的あるいはソニケータにより気泡の生成された５％ヒトアルブミン）を注入するものである。造影剤により心筋内血流の灌流域は、Ｂモード上の輝度増強領域として表示される。同様に、血流の灌流域の評価あるいは腫瘍の支配血管系を評価するために、腹部領域でも動脈注入によるコントラストエコー法が研究されている。これらのコントラストエコー法を実施する診断装置は、一般検査用の超音波診断装置あるいはさらにワークステーションが用いられる。これにより、目視によりＢモード像の輝度増強を評価したり、あるいは、メモリに記憶された画像データをワークステーション上で適当な処理をしたりした後、輝度レベルの変化を定量評価するようになっている。
【０００４】
また、近年、超音波造影剤を静脈から注入して左心系の評価が可能な超音波造影剤が開発され、これを用いた超音波コントラストエコー法が試みられている。
【０００５】
この超音波造影剤としては、塩野義製薬株式会社により輸入販売されている、『５％人血清アルブミンを超音波処理するときに生成するアルブミン膜の中に、空気を封じ込めた平均粒子径約４μｍの空気小球体』（販売名：アルブネックス注５ｍｌ）がある。
【０００６】
この静脈注入によるコントラストエコー法は、現在、試験，研究段階であり、今後、頭部・心腔・腹部などの診断で、その有用性に期待が高まっている。
【０００７】
【発明が解決しようとする課題】
このように近年、超音波造影剤を用いた診断法が盛んになってきているが、これらの診断法にあっては、超音波造影剤の超音波反射信号に対する増強作用にのみ着目したものに止まっており、造影剤の特異な性質に着目した計測は未だ行われていなかった。
【０００８】
本発明はこのような状況に鑑みてなされたもので、気体から成る泡を成分とする超音波造影剤の体積変化に因って超音波エコー信号の散乱強度が変化することを利用して、生体内の機能的情報を的確且つ簡便に得ることができる超音波診断装置を提供することを、その目的とする。
【０００９】
【課題を解決するための手段】
前記目的を達成するため、本発明の超音波診断装置は、気体から成る泡を含んだ超音波造影剤を注入した生体に超音波信号を照射する送信手段と、前記超音波信号の生体からの反射信号に基づいて断層面の画像データを得る受信処理手段と、前記画像データを表示する断層像表示手段と、前記断層像表示手段によって表示された断層像上に所望のＲＯＩを設定するＲＯＩ設定手段と、前記断層面の経時的に収集される複数枚の画像データに対して前記ＲＯＩ内の領域に含まれる画素の輝度を画像毎に演算する輝度演算手段と、この輝度演算手段により演算された時系列の輝度のデータから、前記生体の心拍由来の変化成分を表すパラメータを演算するパラメータ演算手段と、このパラメータ演算手段の演算結果を表示するパラメータ表示手段とを備え、前記パラメータ演算手段は、前記心拍由来の変化成分自体を抽出するハイパスフィルタ又はバンドパスフィルタを備たことを要部とする。
【００１０】
【作用】
生体には気体から成る泡を含んだ超音波造影剤が注入される。この状態で生体に超音波信号を照射し、その反射信号に基づいて断層面の画像データが得られ、断層像が表示される。操作者は、表示された断層像上に所望のＲＯＩを設定する。断層面の経時的に収集される複数枚の画像データに対して、設定したＲＯＩ内の領域に対応する画素の輝度が例えば平均値として画像毎に演算される。この時系列の輝度データに基づいて生体の心拍由来の変化成分を表すパラメータ（例えば、心拍に呼応して超音波造影剤の体積が変化することに伴う心拍由来の変化成分）が演算され、この演算結果が数値やグラフとして表示される。これにより、超音波造影剤の輝度増強作用による診断情報に止まらず、造影剤が気体の泡を成分とすることに伴う診断情報も得ることができる。
【００１１】
【実施例】
本発明の一実施例を図１〜図９に基づいて説明する。なお、この実施例は心機能のパラメータ演算に関して複数の実施態様を含んでいる。
【００１２】
図１に示す超音波診断装置は、フェーズド・アレイ・タイプの超音波プローブ１０（以下、「プローブ」という）と、このプローブ１０を介して被検体との間で超音波信号を送受する送受信部１１と、この送受信部１１が得た超音波エコー信号に基づいてデータ処理を行うデータ処理部１２と、このデータ処理部１２に接続された入力器１３及び表示モニタ１４とを備えている。送受信部１１及びデータ処理部１２のステアリングライン制御、動作タイミング制御などは制御回路１５によって指令される。
【００１３】
プローブ１０は超音波信号の送受を兼ねるもので、複数の振動子（圧電体）を一方向に配列させた超音波振動子群１０Ａを備えている。
【００１４】
送受信部１１は、その送信系として基準発振器２０、送信遅延回路群２１、及びパルサ回路群２２を備えており、パルサ回路群２２の出力側が超音波振動子群１０Ａにチャンネル毎に接続されている。これにより、各振動子はパルサ回路群２２から供給された駆動電圧パルスを受けて励振し、被検体に向けて超音波ビーム信号を照射する。なお、送信フォーカスを掛けない場合、送信遅延回路群２１を省くこともできる。
【００１５】
さらに送受信部１１は、受信系として、プリアンプ群２３、受信遅延回路群２４、加算回路２５、対数圧縮回路２６、検波回路２７、及びＡＤＣ（Ａ／Ｄ変換器）２８を備えており、プリアンプ群２３の入力端が超音波振動子群１０Ａにチャンネル毎に接続されている。生体内から反射してきた超音波エコー信号は、プローブ１０で受信され、振動子群１０Ａにより電気量のエコー信号に変換される。このエコー信号はプリアンプ群２３でチャンネル毎に増幅された後、受信遅延回路群２４でチャンネル毎に遅延が掛けられ、加算回路２５により加算される。
この遅延及び加算（整相加算）により受信フォーカスが掛けられる。この遅延・加算信号はその後、対数圧縮回路２７で圧縮され、検波回路２８で包絡線検波された後、ＡＤＣ２８により超音波ビームラインに沿ったデジタル量の画像データが得られる。超音波ビームを走査することで、２次元の画像データが得られる。
【００１６】
データ処理部１２はＡＤＣ２８から出力された画像データを記憶する画像メモリ（フレームメモリ）３０と、その画像データに後述するグラフデータなどを合成する合成回路３１と、この合成回路３１の出力データをアナログ量に変換するＤＡＣ（Ｄ／Ａ変換器）３２とを備え、このＤＡＣ３２の出力側が表示モニタ１４に接続されている。これにより、画像メモリ３０に一時記憶された画像データは表示モニタ１４にリアルタイムに表示される。
【００１７】
またデータ処理部１２は、ＲＯＩ設定回路３３を備えている。このＲＯＩ設定回路３３はトラックボールなどの入力器１３から出力されるＲＯＩ設定情報（ＲＯＩ位置、大きさ、形状）を受けて、これに対応するＲＯＩ設定信号（ＲＯＩのグラフィックデータを含む）を合成回路３１に出力する。このため、観察者が表示モニタ１４の画面を観察しながら入力器１３を操作することにより、例えば図２に示す如く、画面上の所望位置に所望形状及び大きさのＲＯＩ（関心領域）を設定することができる。
【００１８】
データ処理部１２はさらに、図１に示すように輝度の平均強度演算回路３５、メモリ３６、パラメータ演算回路３７、及びグラフ表示回路３８を備えている。
これらの構成は本発明の要旨を具現化するもので、以下のように構成され機能する。
【００１９】
平均強度演算回路３５はＲＯＩ設定回路３４から出力されるＲＯＩ設定信号を受けてＲＯＩの設定内容（位置、形状、大きさ）を解読するとともに、そのＲＯＩ内の領域を形成する複数の画素の各輝度値（画素値）を画像メモリ３０から読み出す。そして、平均強度演算回路３５は各輝度の平均値（又は積算値でも良い）を算出する。この算出は経時的に（フレーム毎に）連続して得られる画像データに対して行われ、各フレームの算出毎に、その算出値がグラフ表示回路３８及びメモリ３６に送られるようになっている。
【００２０】
なお、この平均強度演算回路３５の演算においては、上述の如く、対数圧縮処理がなされた画像輝度データ（すなわち超音波の散乱強度）のまま近似的に輝度平均値又は積算値を算出するとしてもよいが、以下のように補正しながら算出しても良い。この補正は
【数１】
ｙ＝ｅｘｐ（αｘ）
ここで、
α：定数 ｌｏｇ（１０）／１０（ｌｏｇは自然対数を表す）、
ｙ：散乱強度、
ｘ：画像輝度
の式を用いて対数圧縮データを元に戻す逆変換処理を行うものである。従って、この補正処理を伴う場合、演算回路３５はそのデータ（散乱強度）についてＲＯＩ内の平均値又は積算値を算出することになる。
【００２１】
なお、このように平均強度演算回路３５を作動させている間、その演算結果に誤差を生じさせる恐れのある処理（例えば画像の輪郭等を強調表示するための信号処理）は、制御回路１５によって禁止される。
【００２２】
グラフ表示回路３８は、平均強度演算回路３５から送られてくる、経時的に連続したＲＯＩ内輝度の平均強度データを入力するとともに横軸を時間、縦軸を平均強度Ｂavとするグラフデータをリアルタイムに作成し、このグラフデータを合成回路３１に出力するようになっている。合成回路３１は、リアルタイムに供給されるグラフデータとＲＯＩを設定した超音波断層像データとを例えば分割表示態様の画像データに合成する。このため、表示モニタ１４には例えば図３に示す如く、平均強度Ｂavのグラフを左側、測定している断層像を右側に置いた表示画面が現れる。平均強度Ｂavのグラフにおける時間軸のスパンは数十秒から数分のオーダーであり、造影剤注入からの時間経過とともに、このグラフは時々刻々と更新しながら表示される。なお、この平均強度グラフは超音波断層像に対して上下左右何れの位置に分割表示してもよいし、単独表示するようにしてもよい。
【００２３】
本発明で得られる平均強度曲線は、超音波造影剤が気体からなる泡でできていることに因り、通常のＸ線診断の輝度変化曲線に加えて、心拍由来の情報をも含んでいる。
【００２４】
これを詳述すると、以下のようである。超音波信号の散乱強度は造影剤の体積濃度により変化する。注入直後の濃度は高いが、時間経過とともにその濃度は低下するので、この濃度変化に対応して散乱強度も変化する。例えば、図３中の平均強度Ｂavの曲線において、点線で示した部分が造影剤自体の濃度に由来しており、時間ｔo で造影剤を注入した場合、ＲＯＩ内輝度の平均強度は時間と伴に徐々に上がって最大値に到達し、その後なだらかに低下する。この曲線部分（以下、「造影剤濃度の変化成分ＢCT」という）は、Ｘ線診断装置などで通常得られる輝度変化曲線（Time Density Curve:TDC) と同等なものである。
【００２５】
さらに、本発明の超音波造影剤は泡を成分としているため、圧力変化により個々の泡の体積が変化する。つまり、心臓などのように圧力が小刻みにかつ大きく変化する診断対象の場合、この圧力変化に伴って小刻みな体積変化が合わせて生じ、従って平均強度Ｂavも小刻みに変動する。診断対象が心臓の場合、かかる小刻みな平均強度の変動は心拍に由来する成分（以下、「心拍由来の変化成分ＢHB」という）となり、図３に示す如く、ＴＤＣに相当する造影剤濃度の変化成分ＢCT（点線の曲線）に重畳されている（実線の曲線）。この両者の散乱強度（曲線ＢCTとＢHB）の変化には、時間変化の速度に大きな差がある。造影剤濃度の変化成分ＢCTは時間的にゆっくり変化し、数十秒〜数分持続するのに対して、心拍由来の変化成分ＢHBは心拍に同期するため、その周期は１秒以下である。
【００２６】
上記造影剤濃度の変化成分及び心拍由来の変化成分を個別に抽出したり識別したりするために、前述したメモリ３６及びパラメータ演算回路３７が設けられている。メモリ３６には、平均強度演算回路３５から出力された輝度の平均強度の時系列データＢavが記憶される。
【００２７】
パラメータ演算回路３７は、メモリ３６に記憶された平均強度の時系列データを読み出して心機能に関する指令パラメータを演算する。時系列データの読出しタイミングは、入力器１３を介して演算が指令された時間でもよいし、予め決めた適宜な時間（例えば造影剤注入後の一定時間経過後）でもよい。また、演算するパラメータも入力器１３からその都度指令するようにしてもよいし、予め決めておいてもよい。
【００２８】
まず、心拍由来の変化成分そのものを抽出する場合、図４に示す如く、パラメータ演算回路３７はデジタル型のフィルタ３７ａを有する。フィルタ３７ａがハイパスフィルタの場合、そのカットオフ周波数はほぼ心拍周期の周波数未満のある値に設定される。フィルタ３７ａがバンドパスフィルタの場合、その通過帯域は心拍の周波数成分に設定される。
【００２９】
造影剤濃度の変化成分のみを抽出する場合、パラメータ演算回路３７は図５に示す如く、デジタル型のバンドパスフィルタ３７ｂと減算器３７ｃとを有する。
バンドパスフィルタ３７ｂは心拍由来の変化成分ＢHBを抽出するもので、減算器３７ｃにより、その成分ＢHBを元の時系列のデータＢavから減算することで造影剤濃度の変化成分ＢCTのみが抽出される。
【００３０】
なお、造影剤濃度の変化成分ＢCTのみを抽出するには、上記の減算する構成のほか、単にローパスフィルタを用いて心拍由来の変化成分ＢHBを平均強度Ｂavから抽出するようにしてもよい。
【００３１】
このように抽出された心拍由来の変化成分ＢHB及び造影剤濃度の変化成分ＢCTはグラフ表示回路３８に送られる。グラフ表示回路３８は前述した平均強度Ｂavのグラフデータ作成のほか、今回抽出した変化成分ＢHB、ＢCTのグラフデータ（横軸は時間軸）を作成して合成回路３１に送る。このため、表示モニタ１４には断層像などと伴に、それらの変化成分のグラフが表示され、心機能の状態を表すパラメータとして目視観察することができる。
【００３２】
さらに、パラメータ演算回路３７は心圧の最大値と最小値の比を演算することもできる。散乱体からの信号の散乱強度はその体積濃度に比例することが知られている。また気体の体積と圧力の関係はボイルシャルルの法則から互いに反比例する。このことから心拍周期内の時間における散乱強度の最大値と最小値の比をとればその間の心圧の最大値と最小値の比が得られることになる。
【００３３】
この心圧の最大値と最小値の比を演算する場合、パラメータ演算回路３７はＣＰＵ機能を備え、図６に示す処理を心周期ごとに行う。すなわち、メモリ３６に記憶されている各フレーム（時系列）のＲＯＩ内の輝度の平均強度データＢavの中から、１心周期分のデータを読み込む（同図ステップＳ１）。次いで、読み込んだ１心周期分の時系列データの中の最大値Ｂmax 及び最小値Ｂmin を演算する（同図ステップＳ２、Ｓ３）。さらに、最大値Ｂmax 及び最小値Ｂmin の比「Ｂmax ／Ｂmin 」を演算し、その演算値をグラフ表示回路３８に出力する（同図ステップＳ４）。
【００３４】
グラフ表示回路３８は上記輝度比「Ｂmax ／Ｂmin 」のデータを入力して、内部メモリに記憶していた平均強度曲線Ｂavのデータに、例えば図示しない計時信号を用いて心時相を合わせたグラフ（図７参照）のデータを作成する。作成されたグラフデータは合成回路３１に送られ、断層像などの他の必要な画像データと一緒にフレーム画像に合成される。
【００３５】
この結果、表示モニタ１４には例えば図７に示す如く、横軸を時間とし、下側に輝度の平均強度曲線Ｂav、これに心時相を合わせた状態で上側に輝度比「Ｂmax ／Ｂmin 」のグラフを配置した画像が得られる。これにより、心臓のポンプ機能の状態に関する一つの指標が好適に得られ、モニタで観察することができる。したがって、心臓のポンプ機能を判断する上で、従来のように、カテーテルを使って圧力センサを左室又は右室に留置させて心圧を測定するという困難な手法を用いなくても済む。
【００３６】
なお、上記の輝度比を演算する際、対数圧縮を掛けた数値そのものを用いる場合、輝度差「Ｂ ｍａｘ −Ｂ ｍｉｎ 」を演算することで輝度比と等価な解析を行うことになる。
【００３７】
さらに、超音波造影剤の散乱強度（エコー信号）の変化は前述したように２種の変化成分ＢCT及びＢHBの和として得られる。この内、造影剤濃度の変化成分ＢCTの時間変化は非常に緩やかであるから、その時間微分値は小さく無視できる。従って、輝度の平均強度データＢavに適宜な微小時間幅ｄｔ（例えば２秒以下）の時間微分演算を施すと、その結果は心拍由来の変化成分ＢCTの微分値にほぼ一致し、心圧の時間変化（ｄｐ／ｄｔ）が近似的に得られる。
【００３８】
この心圧の時間変化（ｄｐ／ｄｔ）を計測する場合、パラメータ演算回路３７は例えば図８に示す如く、デジタル型の微分回路３７ｄを備えている。この微分回路３７ｄで演算された心圧の時間変化（ｄｐ／ｄｔ）は同様にグラフ表示回路３８に送られ、内臓メモリに記憶されていた輝度の平均強度データＢavに心時相（時間）を合わせた状態で合成されたグラフデータに変換される。これにより、表示モニタ１４には、計測対象の断層像などと伴に、例えば図９に示すグラフが表示され、輝度の平均強度Ｂavの変化に対応する心圧の時間変化（ｄｐ／ｄｔ）を目視で観察することができる。これによって、心臓のポンプ機能に対する良い指標が得られる。
【００３９】
なお、上記の如く心圧の時間変化（ｄｐ／ｄｔ）を計測するにあたり、造影剤濃度の変化成分の影響を殆ど完全に低減するには、平均強度Ｂavのデータにハイパスフィルタやバンドパスフィルタを掛けて心拍由来の変化成分ＢHBのみを抽出し、この抽出データを時間微分するパラメータ演算回路の構成にしてもよい。
【００４０】
本発明に係る超音波診断装置の別の実施例を図１０及び図１１に基づき説明する。この実施例は、前述した心圧の時間変化（ｄｐ／ｄｔ）を計測する場合の、表示状態の改善に関する。図１と同一又は同等の構成要素には同一符号を用いて説明を省略又は簡単化する。
【００４１】
図１０に示す超音波診断装置は、図１記載の構成に加えて、ＥＣＧ（心電計）４０、ＥＣＧ用アンプ４１、ＡＤＣ４２及び表示制御回路４３を備える。ＥＣＧ４０により収集された被検体の心電図信号はＥＣＧ用アンプ４１、ＡＤＣ４２を介して画像メモリ３０に記憶される。従って、プローブ１０を介してエコー信号を収集すると同時に心電図信号も収集され、画像メモリ３０には画像データと心電図データが蓄えられる。
【００４２】
この画像データ及び心電図データは合成回路３１に送られる。これと伴に、合成回路３１にはグラフ表示回路３８から前述した如く輝度の平均強度Ｂavと心圧の時間変化（ｄｐ／ｄｔ）のグラフデータ（図９参照）が送られてくる。表示制御回路４３は入力器１３からオペレータの指令によって与えられる拡大／縮小表示などの表示指令信号を受けて、合成回路３１に、対応する表示制御信号を送るようになっている。これにより、合成回路３１は画像データやグラフデータを指令された適宜な態様で合成するとともに、拡大／縮小処理も行うようになっている。
【００４３】
この結果、オペレータが入力器１３の操作により拡大表示を指令すると、例えば図１１に示す如く、拡大された平均強度Ｂav、心圧の時間変化（ｄｐ／ｄｔ）、及び心電図のグラフ像が同時・拡大表示される。したがって、表示画面から時相評価を行い易くなるとともに、心圧変化の数値的評価のみならず、時間的変化をも容易に評価可能になる。
【００４４】
本発明のさらに別の実施例を図１２及び図１３に基づき説明する。この実施例は前述した心圧の時間変化（ｄｐ／ｄｔ）を２次元画像のまま処理／表示するようにしたものである。
【００４５】
図１２に示す超音波診断装置は、そのデータ処理部１２において、図１に示した平均強度演算回路３５、メモリ３６、パラメータ演算回路３７、及びグラフ表示回路３８に代えて、画像データ出力回路４９、空間フィルタ５０、メモリ５１、及び画像間演算回路５２を備える。
【００４６】
画像データ出力回路４９は、指令されたＲＯＩの内部領域を形成する画素の２次元画像データＢ2Dを画像メモリ３０から読み出して空間フィルタ５１に出力するようになっている。この２次元画像データＢ2Dは、空間フィルタ５０により、２次元画像のまま、スペックルノイズ除去のフィルタリングが施され、スペックルノイズの影響が著しく低減されたＲＯＩ内画像データが生成される。このフィルタリングされた画像データＢ2Dはメモリ５１に送られ、一時記憶される。画像間演算回路５２はメモリ５１に記憶された、連続する２枚の画像間の画素同士の輝度差又は強度（輝度）比を演算し、その演算値を新たな輝度値とする、ＲＯＩ内の２次元画像データＢ2D ＊を生成するようになっている。この２次元画像データＢ2D ＊の各画素は、各画素値の時間微分値であり、前述した心圧の時間変化（ｄｐ／ｄｔ）の空間分布の画素を形成する。
【００４７】
このＲＯＩ内領域に相当する２次元画像データＢ2D ＊は合成回路３１、ＤＡＣ３２を介して表示モニタ１４に、例えば図１３に示す如く表示される。これにより、心圧の時間変化（ｄｐ／ｄｔ）の空間分布を実時間で観測することができる。
【００４８】
さらに本発明の他の実施例を図１４に基づき説明する。
【００４９】
図１４に示す超音波診断装置は、データ処理部１２を外部機器として備え、送受信部１１からオンラインで画像データを送信可能にしたものである。このために、送受信部１１の出力側に画像メモリ６０を備えて、収集した画像データを一度蓄積可能になっている。この画像メモリ６０の読み出し側はインターフェース６１、６２を介してデータ処理部１２（外部機器）に接続され、制御回路１５からの指令によって画像データをオンライン伝送させる。
【００５０】
これにより、データ処理部１２の画像データ処理をバッチ処理で行うことができる一方、データ処理部１２での各コンポーネントを画像解析専用の装置に組み込むこともできる。なお、送受信部１１とデータ処理部１２間の画像データ伝送はオフライン伝送であってもよい。
【００５１】
【発明の効果】
以上説明したように、本発明の超音波診断装置によれば、断層面の経時的に収集される複数枚の画像データに対して、設定したＲＯＩ内の領域に対応する画素の輝度を画像毎に演算し、この時系列の輝度データに基づいて生体の心拍由来の変化成分を表すパラメータ（例えば、心拍に呼応して超音波造影剤の体積が変化することに伴う心拍由来の変化成分）を演算し、この演算結果を数値やグラフとして表示するようにしたので、超音波造影剤の輝度増強作用による診断情報に止まらず、造影剤が気体の泡を成分とすることに伴う心圧やその最大値、最小値の比など、生体内の機能的情報を的確且つ簡便に得ることができる。
【図面の簡単な説明】
【図１】本発明の実施例に係る超音波診断装置のブロック図。
【図２】断層像上に設定するＲＯＩの説明図。
【図３】ＲＯＩ内の輝度の平均強度グラフの表示例を示す図。
【図４】パラメータ演算回路の位置態様を示すブロック図。
【図５】パラメータ演算回路の別の態様を示すブロック図。
【図６】パラメータ演算回路の更に別の態様に係る処理を示す概略フローチャート。
【図７】図６に示す処理に係る画像例を示す図。
【図８】パラメータ演算回路の更に別の態様に係るブロック図。
【図９】図８に示す構成に係る画像例を示す図。
【図１０】本発明の他の実施例に係る超音波診断装置のブロック図。
【図１１】図１０の構成に係る装置の画像例を示す図。
【図１２】本発明の更に他の実施例を示す超音波診断装置のブロック図。
【図１３】図１２の構成に係る装置の画像例を示す図。
【図１４】本発明の更に他の実施例を示す超音波診断装置のブロック図。
【符号の説明】
１０ プローブ
１１ 送受信部
１２ データ処理部
１３ 入力器
１４ 表示モニタ
１５ 制御回路
３０ 画像メモリ
３１ 合成回路
３２ ＤＡＣ
３４ ＲＯＩ設定回路
３５ 平均強度演算回路
３６ メモリ
３７ パラメータ演算回路
３８ グラフ表示回路
３７ａ，３７ｂ フィルタ
３７ｃ 減算器
３７ｄ 微分回路
４９ 画像データ出力回路
５０ 空間フィルタ
５１ メモリ
５２ 画像間演算回路
Claims (8)
- 気体から成る泡を含んだ超音波造影剤を注入した生体に超音波信号を照射する送信手段と、
前記超音波信号の生体からの反射信号に基づいて断層面の画像データを得る受信処理手段と、
前記画像データを表示する断層像表示手段と、
前記断層像表示手段によって表示された断層像上に所望のＲＯＩを設定するＲＯＩ設定手段と、
前記断層面の経時的に収集される複数枚の画像データに対して前記ＲＯＩ内の領域に含まれる画素の輝度を画像毎に演算する輝度演算手段と、
この輝度演算手段により演算された時系列の輝度のデータから、前記生体の心拍由来の変化成分を表すパラメータを演算するパラメータ演算手段と、
このパラメータ演算手段の演算結果を表示するパラメータ表示手段とを備え、
前記パラメータ演算手段は、前記心拍由来の変化成分自体を抽出するハイパスフィルタ又はバンドパスフィルタを備えることを特徴とする超音波診断装置。 - 気体から成る泡を含んだ超音波造影剤を注入した生体に超音波信号を照射する送信手段と、
前記超音波信号の生体からの反射信号に基づいて断層面の画像データを得る受信処理手段と、
前記画像データを表示する断層像表示手段と、
前記断層像表示手段によって表示された断層像上に所望のＲＯＩを設定するＲＯＩ設定手段と、
前記断層面の経時的に収集される複数枚の画像データに対して前記ＲＯＩ内の領域に含まれる画素の輝度を画像毎に演算する輝度演算手段と、
この輝度演算手段により演算された時系列の輝度のデータから、前記生体の心拍由来の変化成分を表すパラメータを演算するパラメータ演算手段と、
このパラメータ演算手段の演算結果を表示するパラメータ表示手段とを備え、
前記パラメータ演算手段は、前記心拍由来の変化成分の一定時間内での最大値と最小値の比若しくは差分又は当該比若しくは差分の時間の関数値を前記パラメータとして演算する手段であることを特徴とする超音波診断装置。 - 気体から成る泡を含んだ超音波造影剤を注入した生体に超音波信号を照射する送信手段と、
前記超音波信号の生体からの反射信号に基づいて断層面の画像データを得る受信処理手段と、
前記画像データを表示する断層像表示手段と、
前記断層像表示手段によって表示された断層像上に所望のＲＯＩを設定するＲＯＩ設定手段と、
前記断層面の経時的に収集される複数枚の画像データに対して前記ＲＯＩ内の領域に含まれる画素の輝度を画像毎に演算する輝度演算手段と、
この輝度演算手段により演算された時系列の輝度のデータから、前記生体の心拍由来の変化成分を表すパラメータを演算するパラメータ演算手段と、
このパラメータ演算手段の演算結果を表示するパラメータ表示手段とを備え、
前記パラメータ演算手段は、前記時系列の輝度のデータに基づいて前記心拍由来の変化成分の時間微分値又は当該微分値の時間の関数値を前記パラメータとして演算する手段であることを特徴とする超音波診断装置。 - 気体から成る泡を含んだ超音波造影剤を注入した生体に超音波信号を照射する送信手段と、
前記超音波信号の生体からの反射信号に基づいて断層面の画像データを得る受信処理手段と、
前記画像データを表示する断層像表示手段と、
前記断層像表示手段によって表示された断層像上に所望のＲＯＩを設定するＲＯＩ設定手段と、
前記断層面の経時的に収集される複数枚の画像データに対して前記ＲＯＩ内の領域に含まれる画素の輝度を画像毎に演算する輝度演算手段と、
この輝度演算手段により演算された時系列の輝度のデータから、前記生体の心拍由来の変化成分をパラメータとして演算するパラメータ演算手段と、
このパラメータ演算手段の演算結果を表示するパラメータ表示手段と、
前記心拍由来の変化成分を前記時系列の輝度のデータから減じて前記超音波造影剤の体積濃度の変化成分を求める手段と、を備えることを特徴とする超音波診断装置。 - 前記輝度演算手段により演算された時系列の輝度のデータをグラフとして表示する輝度グラフ表示手段を更に備える請求項１〜４の何れか一項に記載の超音波診断装置。
- 前記パラメータは、前記超音波造影剤の濃度に関する前記生体の心拍由来の変化成分である請求項１〜４の何れか一項に記載の超音波診断装置。
- 前記パラメータ演算手段は、前記時系列の輝度のデータから前記心拍由来の変化成分を除去するローパスフィルタを備える請求項４に記載の超音波診断装置。
- 前記パラメータ表示手段は、前記パラメータ演算手段により演算されたパラメータの演算結果を数値及びグラフの少なくとも一方で表示する手段である請求項１〜４の何れか一項に記載の超音波診断装置。」
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP11875895A JP3697292B2 (ja) | 1995-05-17 | 1995-05-17 | 超音波診断装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP11875895A JP3697292B2 (ja) | 1995-05-17 | 1995-05-17 | 超音波診断装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08308831A JPH08308831A (ja) | 1996-11-26 |
| JP3697292B2 true JP3697292B2 (ja) | 2005-09-21 |
Family
ID=14744336
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP11875895A Expired - Fee Related JP3697292B2 (ja) | 1995-05-17 | 1995-05-17 | 超音波診断装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3697292B2 (ja) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4473981B2 (ja) * | 1999-07-21 | 2010-06-02 | 株式会社日立メディコ | 超音波診断装置 |
| JP4497611B2 (ja) * | 1999-12-27 | 2010-07-07 | 株式会社東芝 | 超音波診断装置 |
| EP1387637B1 (en) * | 2001-04-06 | 2007-10-31 | Bracco Research S.A. | Apparatus for measuring local physical parameters in a fluid filled cavity |
| JP4113702B2 (ja) * | 2001-12-04 | 2008-07-09 | 株式会社東芝 | 超音波診断装置及び超音波信号解析装置 |
| JP4679141B2 (ja) * | 2004-12-27 | 2011-04-27 | 株式会社東芝 | 超音波診断装置および超音波診断画像の表示方法 |
| JP4751090B2 (ja) * | 2005-04-07 | 2011-08-17 | 株式会社東芝 | 超音波診断装置 |
| JP4764521B2 (ja) * | 2010-10-12 | 2011-09-07 | 株式会社東芝 | 超音波診断装置 |
| JP2018086291A (ja) * | 2012-12-28 | 2018-06-07 | キヤノン株式会社 | 被検体情報取得装置、表示方法、プログラム、処理装置 |
| JP6292836B2 (ja) | 2012-12-28 | 2018-03-14 | キヤノン株式会社 | 被検体情報取得装置、表示方法、プログラム、処理装置 |
| JP5907185B2 (ja) * | 2014-01-27 | 2016-04-26 | コニカミノルタ株式会社 | 動態画像診断支援システム |
| WO2018023336A1 (zh) * | 2016-08-01 | 2018-02-08 | 深圳迈瑞生物医疗电子股份有限公司 | 超声弹性测量显示方法及系统 |
-
1995
- 1995-05-17 JP JP11875895A patent/JP3697292B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH08308831A (ja) | 1996-11-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3625305B2 (ja) | 超音波診断装置 | |
| US8460192B2 (en) | Ultrasound imaging apparatus, medical image processing apparatus, display apparatus, and display method | |
| EP2234543B1 (en) | Quantification analisys of immobilized contrast agent in medical imaging applications | |
| JP3863414B2 (ja) | 超音波診断装置 | |
| JP5632203B2 (ja) | 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム | |
| JP5680654B2 (ja) | 超音波診断装置及び超音波画像表示方法 | |
| JP2000333956A (ja) | 超音波診断装置及び超音波画像の表示切換方法 | |
| US20050267371A1 (en) | Ultrasonic diagnostic apparatus | |
| JP4945040B2 (ja) | 超音波診断装置 | |
| Schminke et al. | Three-dimensional power-mode ultrasound for quantification of the progression of carotid artery atherosclerosis | |
| US6612989B1 (en) | System and method for synchronized persistence with contrast agent imaging | |
| US7666141B2 (en) | Ultrasonic diagnostic apparatus | |
| EP1514516B1 (en) | Ultrasonic diagnostic equipment | |
| JPH0924047A (ja) | 超音波診断装置 | |
| JP3697292B2 (ja) | 超音波診断装置 | |
| US6607490B2 (en) | Ultrasonic diagnostic apparatus and control method thereof | |
| JPH0759781A (ja) | 超音波診断装置 | |
| JP4521204B2 (ja) | 超音波診断装置、画像処理装置及び超音波画像撮影方法 | |
| JP4113702B2 (ja) | 超音波診断装置及び超音波信号解析装置 | |
| JP3689095B2 (ja) | 超音波診断装置 | |
| JPH10118065A (ja) | 超音波診断装置 | |
| JP4497611B2 (ja) | 超音波診断装置 | |
| Khairalseed et al. | Software-based Processing for Contrast-Enhanced Ultrasound Imaging using Pulse-Inversion Spectral Deconvolution | |
| JPH08623A (ja) | 超音波診断装置 | |
| JP5454901B2 (ja) | 超音波診断装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050105 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050307 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050329 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050530 |
|
| 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: 20050621 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050704 |
|
| LAPS | Cancellation because of no payment of annual fees |