JP2013244162A - 超音波診断装置 - Google Patents
超音波診断装置 Download PDFInfo
- Publication number
- JP2013244162A JP2013244162A JP2012119579A JP2012119579A JP2013244162A JP 2013244162 A JP2013244162 A JP 2013244162A JP 2012119579 A JP2012119579 A JP 2012119579A JP 2012119579 A JP2012119579 A JP 2012119579A JP 2013244162 A JP2013244162 A JP 2013244162A
- Authority
- JP
- Japan
- Prior art keywords
- region
- transmission
- unit
- reception
- sound speed
- 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
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
【課題】血管等の構造物の影響を極力抑えた音速マップの生成を短時間で行う。
【解決手段】表示部に表示されているBモード画像110上に着目領域112を設定する。ドプラ走査モードで超音波の送受信を実行して、着目領域112内のカラードプラ情報を生成する。カラードプラ情報に基づき、着目領域112内における運動体領域117を除いた領域を音速測定領域119として決定する。音速測定領域119内に複数の昇進焦点125を設定して超音波の送受信を行い、音速測定領域119に対応するRFデータを生成する。音速測定領域119に対応するRFデータに基づき局所音速値を算出することで、音速測定領域119の音速マップを生成する。
【選択図】図9
【解決手段】表示部に表示されているBモード画像110上に着目領域112を設定する。ドプラ走査モードで超音波の送受信を実行して、着目領域112内のカラードプラ情報を生成する。カラードプラ情報に基づき、着目領域112内における運動体領域117を除いた領域を音速測定領域119として決定する。音速測定領域119内に複数の昇進焦点125を設定して超音波の送受信を行い、音速測定領域119に対応するRFデータを生成する。音速測定領域119に対応するRFデータに基づき局所音速値を算出することで、音速測定領域119の音速マップを生成する。
【選択図】図9
Description
本発明は、Bモード画像内に設定された着目領域の音速マップを生成する超音波診断装置に関する。
医療分野では、被検体の内部を観察して診断を行うために超音波診断装置がよく用いられている。超音波診断装置は、超音波の送受信を行う超音波探触子(プローブ)を有しており、超音波プローブから被検体内に向けて超音波を送信し、被検体内で反射された超音波を超音波探触子で受信することで、この受信信号に基づきBモード画像を生成することができる。Bモード画像は、超音波エコーの振幅を点の明るさ(輝度)により表した画像である。このBモード画像により、被検体内に存在する構造物(例えば、内臓や病変組織等)の輪郭を観察することができる。
近年、被検体内の診断部位を精度よく診断するために、Bモード画像上で医師が指定した着目領域(診断部位)における局所音速を測定し、この着目領域における音速分布を表す音速マップを生成して、Bモード画像に重畳して表示している。このような着目領域における局所音速を測定する方法は例えば特許文献1に開示されている。
特許文献1記載の局所音速の算出方法では、着目領域の格子点XROIと超音波探触子との間に複数の格子点A1,A2,…を設定して、各格子点(XROI、A1,A2,…)における最適音速値を判定する。ここで、最適音速値とは、画像のコントラスト値、シャープネス値が最も高くなる音速値である。次いで、各格子点における最適音速値に基づいて、格子点XROIにおける局所音速を求める。具体的には、ホイヘンスの原理により、格子点XROIからの受信波と、格子点A1,A2,…からの受信波を仮想的に合成した受信波とが一致することを利用する(図11参照)。これにより、局所音速を正確に測定することができるので、高精度な音速マップが得られる。
ところで、着目領域内において、血管や動きの大きい組織(臓器)などの構造物が存在する領域は、上述の音速計算を画像のコントラスト値やシャープネス値に基づいて行う関係上、音速誤差が生じやすい。その結果、音速計算の計算結果の信頼性が十分に得られず、血管等の構造物の影響を強く受けた音速マップが生成されてしまう。
特許文献2には、着目領域から血管等の構造物を除いた領域における弾性画像を生成する超音波診断装置が開示されている。この特許文献2の超音波診断装置では、着目領域内の歪み分布を求めるとともに、超音波の受信信号から血管内の血流等に起因するドプラ信号を抽出して、このドプラ信号に基づき血流の流速等を示すドプラデータを生成する。次いで、先に求めた歪み分布の中から、ドプラーデータに基づき所定の流速以上となっている部分の歪みを除外した歪み分布に基づき、弾性画像を生成している。また、特許文献3にも、弾性画像データから血流領域を検出して除去する旨が開示されている。このため、特許文献2及び3の記載に基づき、着目領域の局所音速の分布を求めるとともに、着目領域のドプラデータを生成して、このドプラデータに基づいて着目領域から血管等の構造物を除いた領域の局所音速値に基づいて音速マップを生成する方法も提案されている。
しかしながら、着目領域内の各点の局所音速値の算出に加えてドプラデータの生成まで行う場合には膨大な演算処理を要するので、音速マップの生成に時間がかかるという問題が生じる。また、音速マップの生成時間を短縮するためには、高性能な演算処理装置(CPU)を超音波診断装置に搭載する必要があるので、超音波診断装置の製造コストが増加するという問題が生じる。
本発明の目的は、血管等の構造物の影響を極力抑えた音速マップの生成を短時間で行うことができる超音波診断装置を提供することにある。
本発明の目的を達成するための超音波診断装置は、超音波を被検体に送信し、かつ被検体によって反射された超音波を受信して超音波検出信号を出力する複数の素子を含む超音波探触子と、複数の素子をそれぞれ駆動する複数の駆動信号を複数の素子に供給し、かつ複数の素子から出力された超音波検出信号に基づき受信データを生成する送受信部であって、複数の素子に供給する複数の駆動信号の遅延量を調節することで超音波ビームを形成する送受信部と、被検体内に着目領域を設定する着目領域設定部と、送受信部を制御して、着目領域内の移動体の情報を得るための超音波の送受信を複数の素子に実行させる第1の送受信制御部と、移動体の情報を得るための超音波の送受信がなされたときに、移動体によりドプラ偏移されて反射された超音波の超音波検出信号に基づいて送受信部が生成した受信データから、移動体の動きに伴うドプラ信号を抽出して、この抽出結果に基づき着目領域内の移動体を含む移動体領域を判別する移動体領域判別部と、移動体領域判別部の判別結果に基づき、着目領域の中から移動体領域を除いた領域を、音速の分布を測定する音速測定領域として決定する音速測定領域決定部と、音速測定領域内に超音波ビームの送信焦点を複数設定し、各送信焦点に対する超音波ビームの送受信を送受信部に行わせることで、音速測定領域の音速測定用の受信データを送受信部に生成させる第2の送受信制御部と、音速測定用の受信データに基づき音速測定領域の音速の分布を示す音速マップを生成する音速マップ生成部と、を備えることを特徴とする。
本発明によれば、着目領域の中から音速誤差が生じやすい移動体領域を除いた音速測定領域の音速マップを生成することができる。
送受信部を制御して、被検体内の断層像情報を得るための超音波の送受信を複数の素子に実行させる第3の送受信制御部と、断層像情報を得るための超音波の送受信がなされたときに送受信部が生成した受信データに基づきBモード画像を生成する画像生成部と、Bモード画像を表示する表示部と、を備え、着目領域設定部は、Bモード画像上に着目領域を設定することが好ましい。Bモード画像を見ながら着目領域を設定することができる。
第1の送受信制御部による送受信部の制御と、移動体領域判別部による移動体領域の判別と、音速測定領域決定部による音速測定領域の決定と、第2の送受信制御部による送受信部の制御と、音速マップ生成部による音速マップの生成と、を含む音速マップ更新処理を実行させる更新制御部を備えることが好ましい。音速マップを再度生成する前に、移動体領域の判別や音速測定領域の決定を再度行うので、血管(血流)や動きの大きい組織等の各種移動体の影響を受けない音速マップを再生成(更新処理)することができる。
更新制御部は、予め定められたフレーム数のBモード画像が画像生成部により生成される毎に、音速マップ更新処理を実行させることが好ましい。音速マップを定期的に更新することができる。
移動体には、着目領域内の血管を流れる血液が含まれており、移動体領域判別部は、ドプラ信号の抽出結果に基づき、着目領域内での血液の流れを示す血流情報を生成して、血流情報に基づき移動体領域を判別することが好ましい。血管(血流)の影響を受けない音速マップを生成することができる。
移動体には、着目領域内で移動する組織が含まれており、移動体領域判別部は、ドプラ信号の抽出結果に基づき、着目領域内での組織の移動を示す組織移動情報を生成して、組織移動情報に基づき移動体領域を判別することが好ましい。移動する組織の影響を受けない音速マップを生成することができる。
移動体領域判別部は、送受信部が生成した受信データから、組織の中で動きの大きさが予め定めたしきい値以上となる組織に対応するドプラ信号を抽出して、ドプラ信号に基づき組織移動情報を生成することが好ましい。組織の中でも動きの大きい組織の影響を受けない音速マップを生成することができる。
表示部は、Bモード画像上に音速マップを重畳表示することが好ましい。
本発明の超音波診断装置は、Bモード画像上に設定された着目領域の中から移動体を含む移動体領域を除いた領域を音速測定領域として決定し、この音速測定領域の局所音速値を求めて音速マップを生成するので、着目領域内の全領域で局所音速値を求める必要が無くなる。その結果、音速マップの生成の際における演算処理を従来よりも低減させることができるので、血管(血流)や動きの大きい組織等の各種移動体の影響を受けない音速マップを従来よりも短時間で生成することができる。
[第1実施形態の超音波診断装置の装置構成]
図1に示すように、超音波診断装置10は、超音波探触子300から被検体OBJ(図10参照)に超音波ビームを送信して、被検体OBJによって反射された超音波ビーム(超音波エコー)を受信して得られた超音波検出信号から超音波画像を作成・表示する装置である。超音波診断装置10は、大別して、CPU100、格納部102、表示部104、表示制御部105、操作入力部(着目領域設定部)200、超音波探触子300、送受信部400、画像信号生成部(画像生成部)500、再生部600、データ解析計測部700、ドプラ信号処理部(移動体領域判別部)701、音速測定領域決定部702を備えている。
図1に示すように、超音波診断装置10は、超音波探触子300から被検体OBJ(図10参照)に超音波ビームを送信して、被検体OBJによって反射された超音波ビーム(超音波エコー)を受信して得られた超音波検出信号から超音波画像を作成・表示する装置である。超音波診断装置10は、大別して、CPU100、格納部102、表示部104、表示制御部105、操作入力部(着目領域設定部)200、超音波探触子300、送受信部400、画像信号生成部(画像生成部)500、再生部600、データ解析計測部700、ドプラ信号処理部(移動体領域判別部)701、音速測定領域決定部702を備えている。
CPU(Central Processing Unit:中央処理装置)100は、操作入力部200からの操作入力に応じて超音波診断装置10の各ブロックの制御を行う。
操作入力部200は、オペレータからの操作入力を受け付ける入力デバイスであり、操作卓202とポインティングデバイス204とを含んでいる。操作卓202は、文字情報(例えば、患者情報)の入力を受け付けるキーボードと、表示モードを切り替える表示モード切り替えボタンと、ライブモードとフリーズモードとの切り替えを指示するためのフリーズボタンと、シネメモリ再生を指示するためのシネメモリ再生ボタンと、超音波画像の解析・計測を指示するための解析・計測ボタンとを含んでいる。
表示モードには、Bモード画像を単独で表示するモードや、Bモード画像と、局所音速の分布を示す音速マップまたは血流や組織などの移動速度を示すカラーフローマッピング像とを表示するモードなどが含まれている。
ポインティングデバイス204は、表示部104の画面上における領域の指定の入力を受け付けるデバイスであり、例えば、トラックボール又はマウスである。ポインティングデバイス204は、Bモード画像上での着目領域(ROI:Region of Interest、関心領域ともいう)の設定に用いられる。なお、ポインティングデバイス204としては、タッチパネルを用いることも可能である。
格納部102は、CPU100により超音波診断装置10の各ブロックの制御を制御するための制御プログラム、パラメータを格納する記憶装置であり、例えば、ハードディスク又は半導体メモリである。
表示部104は、例えば、CRT(Cathode Ray Tube)ディスプレイ又は液晶ディスプレイであり、超音波画像(動画及び静止画)の表示、音速マップ、カラーフローマッピング像及び各種の設定画面などを表示する。表示制御部105は、画像信号生成部500から出力された画像データに基づき表示部104に超音波画像等を表示させる。
超音波探触子300は、被検体OBJに当接させて用いるプローブであり、1次元又は2次元の超音波トランスデューサアレイを構成する複数の素子302を備えている。各素子302は、送受信部400から印加される駆動信号に基づいて超音波ビームを被検体OBJに送信するとともに、被検体OBJから反射される超音波エコーを受信して検出信号を出力する。
各素子302は、圧電性を有する材料(圧電体)の両端に電極が形成されて構成された振動子を含んでいる。振動子を構成する圧電体としては、例えば、PZT(チタン酸ジルコン酸鉛:Pb (lead) zirconate titanate)のような圧電セラミック、PVDF(ポリフッ化ビニリデン:polyvinylidene difluoride)のような高分子圧電素子を用いることができる。振動子の電極に電気信号を送って電圧を印加すると圧電体が伸縮し、この圧電体の伸縮により各振動子において超音波が発生する。例えば、振動子の電極にパルス状の電気信号を送るとパルス状の超音波が発生し、振動子の電極に連続波の電気信号を送ると連続波の超音波が発生する。そして、各振動子において発生した超音波が合成されて超音波ビームが形成される。また、各振動子により超音波が受信されると、各振動子の圧電体が伸縮して電気信号を発生する。各振動子において発生した電気信号は、超音波検出信号として送受信部400に出力される。
なお、素子302としては、超音波変換方式の異なる複数種類の素子を用いることも可能である。例えば、超音波を送信する素子として上述の圧電体により構成される振動子を用いて、超音波を受信する素子として光検出方式の超音波トランスデューサを用いるようにしてもよい。ここで、光検出方式の超音波トランスデューサとは、超音波信号を光信号に変換して検出するものであり、例えば、ファブリーペロー共振器又はファイバブラッググレーティングである。
送受信部400は、超音波の送受信制御を行うものであり、送信回路402、受信回路404、及びA/D変換器406を備えている。送信回路402は、駆動信号をそれぞれ各素子302に送る。受信回路404は、各素子302から入力される検出信号を受信して音線信号を生成する。A/D変換器406は、アナログの音線信号をデジタル信号(RFデータ)に変換する。
画像信号生成部500は、送受信部400から入力されるRFデータに基づき、表示用の画像データ(例えばBモード画像データ)を生成するものであり、信号処理部502、DSC(Digital Scan Converter)504、画像処理部506、画像メモリ508、D/A変換器510を備えている。信号処理部502は、RFデータに基づきBモード画像データを生成する。DSC504は、Bモード画像データ等をテレビジョン信号の走査方式の画像データに変換する。画像処理部506は、画像データに各種の画像処理を施す。画像メモリ508は、画像処理後の画像データを一時的に格納する。D/A変換器510は、デジタルの画像データをアナログの画像信号に変換する。
再生部600は、後述のシネメモリ再生モード時に作動するものであり、シネメモリ602とシネメモリ再生部604とを有している。シネメモリ602は、A/D変換器406に接続しており、A/D変換器406から入力されたRFデータを格納する。シネメモリ再生部604は、シネメモリ602に格納されているRFデータを信号処理部502に送る。
データ解析計測部700は、RFデータを用いてオペレータ指定の解析・計測を行う。ドプラ信号処理部701は、RFデータに基づきカラードプライメージング(血流イメージング、組織ドプライメージング)用、すなわち、被検体OBJ内の血流や組織の移動状態を表すカラードプラ情報を生成する。さらに、ドプラ信号処理部701は、カラードプラ情報からカラーフローマッピング像(カラードプラ断層像ともいう、以下、「CFM像」と略す)を生成する。音速測定領域決定部702は、音速マップ生成の際に、着目領域の中から局所音速値の算出を行う音速測定領域を決定する。
超音波診断装置10は、動作モードとして、ライブモードとシネメモリ再生モードとを有している。ライブモードは、被検体OBJに超音波探触子300を当接させて超音波の送受信を行うことによって得られた超音波画像(動画)の表示、解析・計測を行うモードである。シネメモリ再生モードは、シネメモリ602に格納されているRFデータに基づいて超音波診断画像の表示、解析・計測を行うモードである。
<ライブモード>
次に、ライブモード時における超音波診断処理について説明する。超音波探触子300を被検体OBJに当接させた状態で、操作入力部200からの指示入力により超音波診断が開始される。CPU100は、送受信部400に制御信号を出力して、超音波ビームの被検体OBJへの送信、及び被検体OBJからの超音波エコーの受信を開始させる。また、CPU100は、素子302毎に超音波ビームの送信方向と超音波エコーの受信方向とを設定する。
次に、ライブモード時における超音波診断処理について説明する。超音波探触子300を被検体OBJに当接させた状態で、操作入力部200からの指示入力により超音波診断が開始される。CPU100は、送受信部400に制御信号を出力して、超音波ビームの被検体OBJへの送信、及び被検体OBJからの超音波エコーの受信を開始させる。また、CPU100は、素子302毎に超音波ビームの送信方向と超音波エコーの受信方向とを設定する。
次いで、CPU100は、超音波ビームの送信方向に応じて送信遅延パターンを選択するとともに、超音波エコーの受信方向に応じて受信遅延パターンを選択する。ここで、送信遅延パターンとは、各素子302から送信される超音波によって所望の方向に超音波ビームを形成するために駆動信号に与えられる遅延時間のパターンデータである。また、受信遅延パターンとは、複数の素子302によって受信される超音波によって所望の方向からの超音波エコーを抽出するために検出信号に与えられる遅延時間のパターンデータである。これら送信遅延パターン及び受信遅延パターンは予め格納部102に格納されている。CPU100は、格納部102の中から選択した送信遅延パターン及び受信遅延パターンに従って、送受信部400に制御信号を出力して超音波の送受信制御を行う。
送信回路402は、CPU100からの制御信号に応じて各素子302の駆動信号を生成して、各駆動信号を各素子302にそれぞれ印加する。
図2に示すように、送信回路402は、素子302毎に遅延回路τ1〜τNを有し、CPU100によって選択された送信遅延パターンに基づいて、各素子302に印加する駆動信号を遅延させる。ここで、送信回路402は、複数の素子302から送信される超音波が超音波ビームを形成するように、各素子302に駆動信号を印加するタイミングを調整(遅延)する。また、送信回路402は、超音波ビームの方向(ステア角α)を調整するように、各素子302に駆動信号を印加するタイミングを調整(遅延)する。なお、複数の素子302から一度に送信される超音波が被検体OBJの撮像領域全体に届くように、駆動信号を印加するタイミングを調節してもよい。
図1に戻って、受信回路404は、各素子302から出力される超音波検出信号を受信して増幅する。ここで、各素子302と被検体OBJ内の超音波反射源との間の距離がそれぞれ異なるため、各素子302に反射波が到達する時間が異なる。この際に受信回路404は、遅延回路を備えており、CPU100によって選択された音速(以下、「仮定音速」という)又は音速の分布に基づいて設定される受信遅延パターンに従って、反射波の到達時刻の差(遅延時間)に相当する分だけ、各検出信号を遅延させる。
さらに、受信回路404は、遅延時間を与えた検出信号を整合加算することにより受信フォーカス処理を行う。この際に、超音波反射源とは異なる位置に別の超音波反射源がある場合には、別の超音波反射源からの超音波検出信号は到達時刻が異なるので、別の超音波反射源からの超音波検出信号の位相が打ち消し合う。これにより、超音波反射源からの受信信号が最も大きくなり、フォーカスが合う。このような受信フォーカス処理によって、超音波エコーの焦点が絞り込まれた音線信号(以下、「RF信号」という)が形成される。
A/D変換器406は、受信回路404から出力されるアナログのRF信号をデジタルRF信号であるRFデータに変換する。ここで、RFデータは、受信波の位相情報を含んでいる。A/D変換器406から出力されるRFデータは、信号処理部502とシネメモリ602にそれぞれ入力される。
シネメモリ602は、A/D変換器406から入力されるRFデータを順次格納する。また、シネメモリ602は、CPU100から入力されるフレームレートに関する情報(例えば、超音波の反射位置の深度、走査線の密度、視野幅を示すパラメータ)をRFデータに関連付けて格納する。
信号処理部502は、RFデータに対して、STC(Sensitivity Time gain Control)によって、超音波の反射位置の深度に応じて距離による減衰の補正をした後、包絡線検波処理を施す。これにより、Bモード画像データ(超音波エコーの振幅を点の明るさ(輝度)により表した画像データ)が生成される。
信号処理部502によって生成されたBモード画像データは、通常のテレビジョン信号の走査方式と異なる走査方式によって得られたものである。このため、DSC504は、Bモード画像データを通常の画像データ(例えば、テレビジョン信号の走査方式(NTSC方式)の画像データ)に変換(ラスター変換)する。画像処理部506は、DSC504から入力される画像データに、各種の必要な画像処理(例えば、階調処理)を施す。
画像メモリ508は、画像処理部506から入力される画像データを格納する。D/A変換器510は、画像メモリ508から読み出された画像データをアナログの画像信号に変換して表示制御部105に出力する。これにより、超音波探触子300によって撮影された超音波画像(動画)が表示部104に表示される。
なお、本実施形態では、受信回路404において受信フォーカス処理が施された検出信号をRF信号としたが、受信フォーカス処理が施されていない検出信号をRF信号としてもよい。この場合、複数の素子302から出力される複数の超音波検出信号が、受信回路404において増幅され、増幅された検出信号、即ち、RF信号が、A/D変換器406においてA/D変換されることによってRFデータが生成される。このRFデータは、信号処理部502に供給されるとともに、シネメモリ602に格納される。受信フォーカス処理は、信号処理部502においてデジタル的に行われる。
<シネメモリ再生モード>
次に、シネメモリ再生モードについて説明する。操作卓202のシネメモリ再生ボタンが押下されると、CPU100は、超音波診断装置10の動作モードをシネメモリ再生モードに切り替える。シネメモリ再生モード時には、CPU100は、オペレータからの操作入力により指定されたRFデータの再生をシネメモリ再生部604に指令する。シネメモリ再生部604は、CPU100からの指令に従って、シネメモリ602からRFデータを読み出して信号処理部502に送信する。シネメモリ602から送信されたRFデータは、信号処理部502、DSC504及び画像処理部506において所定の処理(ライブモード時と同様の処理)が施されて画像データに変換された後、画像メモリ508及びD/A変換器510を経て表示部104に出力される。これにより、シネメモリ602に格納されたRFデータに基づく超音波画像(動画又は静止画)が表示部104に表示される。
次に、シネメモリ再生モードについて説明する。操作卓202のシネメモリ再生ボタンが押下されると、CPU100は、超音波診断装置10の動作モードをシネメモリ再生モードに切り替える。シネメモリ再生モード時には、CPU100は、オペレータからの操作入力により指定されたRFデータの再生をシネメモリ再生部604に指令する。シネメモリ再生部604は、CPU100からの指令に従って、シネメモリ602からRFデータを読み出して信号処理部502に送信する。シネメモリ602から送信されたRFデータは、信号処理部502、DSC504及び画像処理部506において所定の処理(ライブモード時と同様の処理)が施されて画像データに変換された後、画像メモリ508及びD/A変換器510を経て表示部104に出力される。これにより、シネメモリ602に格納されたRFデータに基づく超音波画像(動画又は静止画)が表示部104に表示される。
<超音波診断装置の各種機能>
ライブモード又はシネメモリ再生モード時において、超音波画像(動画)が表示されているときに操作卓202のフリーズボタンが押下されると、この押下時に表示されている超音波画像が表示部104に静止画表示される。これにより、オペレータは、着目領域の静止画を観察することができる。
ライブモード又はシネメモリ再生モード時において、超音波画像(動画)が表示されているときに操作卓202のフリーズボタンが押下されると、この押下時に表示されている超音波画像が表示部104に静止画表示される。これにより、オペレータは、着目領域の静止画を観察することができる。
操作卓202の計測ボタンが押下されると、オペレータからの操作入力により指定された解析・計測が行われる。データ解析計測部605は、各動作モード時に計測ボタンが押下された場合に、A/D変換器406又はシネメモリ602から、画像処理が施される前のRFデータを取得し、このRFデータを用いてオペレータ指定の解析・計測(例えば、組織部の歪み解析(硬さ診断)、血流の計測、組織部の動き計測、又はIMT(内膜中膜複合体厚:Intima-Media Thickness)値計測)を行う。データ解析計測部605による解析・計測結果は、DSC504に出力される。DSC504は、データ解析計測部605による解析・計測結果を超音波画像の画像データに挿入して表示制御部105に出力する。これにより、超音波画像と解析・計測結果とが表示部104に表示される。
操作卓202の表示モード切り替えボタンが押下されると、表示モードが単独表示モード、CFM像表示モード、音速マップ表示モード、並列表示モードの間で切り替わる。単独表示モードではBモード画像を単独表示する。CFM像表示モードでは、Bモード画像にCFM像(血流や組織などの移動速度の算出結果)を重畳して表示する(例えば、移動速度に応じて色分け)。音速マップ表示モードでは、Bモード画像に音速マップ(局所音速値の判定結果)を重畳して表示する(例えば、局所音速に応じて色分け又は輝度を変化させる表示、又は局所音速が等しい点を線で結ぶ表示)。並列表示モードでは、Bモード画像と、音速マップの画像またはCFM像とを並べて表示する。
音速マップを観察することで、例えば病変を発見することができる。なお、局所音速の判定結果に基づいて、送信フォーカス処理及び受信フォーカス処理の少なくとも一方を施すことにより得られたBモード画像を表示部104に表示してもよい。また、CFM像を観察することで、血流速度などの血行状態を診断することができる。
<CFM像表示モード>
CFM像表示モードが選択された場合に、送受信部400は、被検体OBJの断層像情報(Bモード画像)を得るためのBモード走査に対応した超音波の送受信と、カラードプラ情報(CFM像)を得るためのドプラモード走査に対応した超音波の送受信とを各素子302に繰り返し実行させる。Bモード走査では、ライブモード時と同様の処理によりBモード画像が生成されて、画像メモリ508及びD/A変換器510を経て表示部104に出力される。
CFM像表示モードが選択された場合に、送受信部400は、被検体OBJの断層像情報(Bモード画像)を得るためのBモード走査に対応した超音波の送受信と、カラードプラ情報(CFM像)を得るためのドプラモード走査に対応した超音波の送受信とを各素子302に繰り返し実行させる。Bモード走査では、ライブモード時と同様の処理によりBモード画像が生成されて、画像メモリ508及びD/A変換器510を経て表示部104に出力される。
ドプラモード走査では、例えば、被検体OBJ内に連続的にドプラ計測用の超音波を送信し、その超音波エコー(エコー信号)を連続的に受信する。このように超音波の送信と受信とを連続的に行うために、送信用の素子302は送信のみ、受信用の素子302は受信のみに用いられる(例えば、特許387597号参照)。なお、ドプラモード走査時の超音波の送受信は、上述の方法に限定されず、各種の方法を用いることができる。
ドプラモード走査時の超音波の送受信によって、血管を流れる血液や移動する組織などの移動体によりドプラ偏移されて反射された超音波の検出信号に基づくRFデータが送受信部400により生成される。このRFデータは、ドプラ信号処理部701に入力される。
〔ドプラ信号処理部〕
図3に示すように、ドプラ信号処理部701は、位相検波器705と、フィルタ部706と、周波数解析部707と、符号化演算部708と、を有しており、カラードプラ情報やCFM像を生成する。
図3に示すように、ドプラ信号処理部701は、位相検波器705と、フィルタ部706と、周波数解析部707と、符号化演算部708と、を有しており、カラードプラ情報やCFM像を生成する。
位相検波器705は、ミキサ及びローパスフィルタを備える。血流や移動する組織などの移動体で反射したエコー信号は、ドプラ効果によって、その周波数にドプラ偏移を受けている。位相検波器705は、そのドプラ周波数について位相検波を行い、RFデータからドプラ信号のみをフィルタ部706に出力する。
図4に示すように、フィルタ部706は、位相検波器705により位相検波されたドプラ信号から、比較的動きの小さい組織に対応するドプラ信号D1を除去し、比較的動きの大きい組織に対応するドプラ信号D2と、血流に対応するドプラ信号D3とを効率良く検出する。この場合に、フィルタ部706はハイパスフィルタとして機能する。なお、図中の符号「HF」はフィルタ部706のハイパスフィルタ特性を示す。このハイパスフィルタ特性HFは本発明のしきい値に相当するものである。従って、被検体OBJ内の組織の中で、動きの大きさがハイパスフィルタ特性HFで定められるしきい値以上となる組織に対応するドプラ信号が得られる。フィルタ部706でフィルタリングされたドプラ信号は周波数解析部707に出力される。
図3に戻って、周波数解析部707は、超音波ドプラ血流計測で用いられている周波数分析法である、FFT法及び自己相関法などを応用するものであり、スキャンされる断層面内の個々のサンプル点における観測時間(時間窓)内での平均速度や最大速度を速度データとして演算する。具体的には、例えば、FFT法又は自己相関法を用いてサンプル各点の平均ドプラ偏移周波数(即ち、その点での観測対象の移動の平均速度)や分散値(ドプラスペクトラムの乱れ度)を、さらにはFFT法を用いてドプラ偏移周波数の最大値(即ち、その点での観測対象の移動の最大速度)などをほぼリアルタイムで演算する。このドプラ周波数の解析結果は、カラードプラ情報(血流情報、組織移動情報(組織ドプラ情報ともいう))として次段の符号化演算部708に出力される。
符号化演算部708は、周波数解析部707から送られてくる、断層面の各サンプル点毎のドプラ偏移周波数fdを、図5に示すような予め指定された速度変換スケールを使って所定ビット数の速度表示データに符号化する。これにより、サンプル点毎の速度表示であるCFM像データが得られる。このCFM像データは、DSC504に出力される。
図3に戻って、DSC504は、CFM像データをテレビジョン信号の走査方式に変換する。DSC504により変換処理されたCFM像データは、画像処理部506により画像処理が施された後、D/A変換器510を経て表示制御部105に出力される。
図6(a)から図6(c)に示すように、表示制御部105は、Bモード走査により得られたBモード画像110上にドプラモード走査により得られたCFM像111を重畳した重畳画像110aを表示部104に表示させる。CFM像111(重畳画像110a)では、血管115や組織116の中で比較的に大きな(上述のしきい値以上の)動きをする組織(以下、高速移動組織という)116aが、他の箇所とは異なる色で表示される。これにより、血管115内の血流115aや高速移動組織116aなどの移動体を含む移動体領域117を判別することができる。
なお、図6に示した例では、ドプラ信号処理部701がBモード画像データの全領域に対応するCFM像データを生成しているが、例えば、後述の着目領域112に対応するCFM像データを生成してもよい。
<音速マップ表示モード>
音速マップ表示モードが選択された場合にも、送受信部400は、Bモード走査に対応した超音波の送受信と、カラードプラ情報(CFM像)を得るためのドプラモード走査に対応した超音波の送受信とを各素子302に実行させる。そして、ドプラ信号処理部701は、カラードプラ情報(またはCFM像データでも可)を音速測定領域決定部702に送る(図1、図3参照)。
音速マップ表示モードが選択された場合にも、送受信部400は、Bモード走査に対応した超音波の送受信と、カラードプラ情報(CFM像)を得るためのドプラモード走査に対応した超音波の送受信とを各素子302に実行させる。そして、ドプラ信号処理部701は、カラードプラ情報(またはCFM像データでも可)を音速測定領域決定部702に送る(図1、図3参照)。
また、図7に示すように音速マップ表示モードが選択された場合には、ポインティングデバイス204を用いて、表示部104に表示されているBモード画像110上に着目領域112が設定される。これにより、音速測定領域決定部702による着目領域112内での音速測定領域119(図9参照)の決定と、CPU100による音速測定領域119の音速マップの生成とが行われる。
図8及び図9(a)から図9(c)に示すように、音速測定領域決定部702は、カラードプラ情報またはCFM像データに基づき、音速測定領域119(図9(c)中の斜線で表示)を決定する。カラードプラ情報またはCFM像データには、局所音速算出時に誤差が生じやすい移動体領域117の位置、形状、大きさ等が含まれている(図9(b)参照)。なお、図9(b),(c)では、移動体領域117の位置、形状、大きさ等を判別し易くするため、カラードプラ情報またはCFM像データの代わりにCFM像111を図示している。
音速測定領域決定部702は、カラードプラ情報またはCFM像データに基づき、先に設定された着目領域112内に移動体領域117が含まれているか否かを判別する。そして、音速測定領域決定部702は、着目領域112内における移動体領域117を除いた領域を音速測定領域119として決定する。この決定結果(音速測定領域119の位置、形状、大きさ等)はCPU100に入力される。
<CPUの機能(音速マップの生成に係る構成)>
CPU100は、格納部102から読み出した各種プログラムを実行することにより、送受信制御部(第1の送受信制御部、第2の送受信制御部、第3の送受信制御部)120、局所音速値算出部121、音速マップ生成部122、更新制御部123として機能する(図8参照)。
CPU100は、格納部102から読み出した各種プログラムを実行することにより、送受信制御部(第1の送受信制御部、第2の送受信制御部、第3の送受信制御部)120、局所音速値算出部121、音速マップ生成部122、更新制御部123として機能する(図8参照)。
送受信制御部120は、単独表示モードが選択されたときには送受信部400を制御してBモード走査に対応した超音波の送受信を各素子302に実行させる。また、送受信制御部120は、CFM像表示モードが選択されたときには送受信部400を制御して、Bモード走査に対応した超音波の送受信と、ドプラモード走査に対応した超音波の送受信とを各素子302に繰り返し実行させる。
さらに、送受信制御部120は、音速マップ表示モードが選択されたときには送受信部400を制御して、音速マップの生成前に、ドプラモード走査に対応した超音波の送受信と、音速測定領域119内の音速測定用のRFデータ(音速測定用の受信データ)を取得するための超音波の送受信とを各素子302に実行させる。音速マップの生成前に、ドプラモード走査に対応した超音波の送受信を実行させることで、ドプラ信号処理部701によるカラードプラ情報またはCFM像データの生成が実行される。これにより、音速マップ生成前に、音速測定領域決定部702により音速測定領域119が決定される。
〔音速測定用のRFデータの取得処理〕
送受信制御部120は、音速測定領域決定部702により決定された音速測定領域119内に、超音波ビームB(図10参照)の送信焦点125を複数設定する。送信焦点125は、局所音速値算出部121による局所音速算出において格子点XROI、格子点A1、格子点A2、・・・となる点である(図11参照)。本実施形態では、送信焦点125は略格子状に設定される。
送受信制御部120は、音速測定領域決定部702により決定された音速測定領域119内に、超音波ビームB(図10参照)の送信焦点125を複数設定する。送信焦点125は、局所音速値算出部121による局所音速算出において格子点XROI、格子点A1、格子点A2、・・・となる点である(図11参照)。本実施形態では、送信焦点125は略格子状に設定される。
次いで、図10に示すように、送受信制御部120は、各送信焦点125にそれぞれ送信フォーカスを行うための送信遅延パターンと、各送信焦点125からの超音波エコーを抽出するための受信遅延パターンを選択する。そして、送受信制御部120は、選択した各送信遅延パターン及び各受信遅延パターンに従って、送受信部400に制御信号を出力して超音波の送受信制御を行う。これにより、各送信焦点125に超音波ビームの焦点を合わせつつ、各送信焦点125からの超音波エコーを各素子302で検出することで、各送信焦点125のRF信号が形成される。なお、本実施形態では、図中の各直線Lsに示すように、超音波ビームの焦点をY方向(被検体OBJの深さ方向)に変えながら1ラインずつ各送信焦点125のRF信号を形成する。各送信焦点125のRF信号は、A/D変換器406によりRFデータ(音速測定用の受信データ)に変換された後、シネメモリ602に格納される。
〔音速測定領域内の局所音速値の算出処理〕
局所音速値算出部121は、特許文献1に開示されたホイヘンスの原理を利用した局所音速算出法を用いて、着目領域112内の個々の領域の局所音速値を求める。
局所音速値算出部121は、特許文献1に開示されたホイヘンスの原理を利用した局所音速算出法を用いて、着目領域112内の個々の領域の局所音速値を求める。
図11(a),(b)は、局所音速値の算出処理を模式的に示す図である。被検体OBJ内の着目領域を代表する格子点をXROI、格子点XROIよりも浅い(即ち、超音波探触子300側の)位置にXY方向に等間隔で配置された格子点をA1,A2,…とし、少なくとも格子点XROIと各格子点A1,A2,…との間の領域RXAの音速はそれぞれ一定と仮定する。
本実施形態では、格子点A1,A2,…からの受信波(それぞれWA1,WA2,…)の(T及び遅延時間ΔT)が既知として、格子点XROIと格子点A1,A2,…の位置関係から格子点XROIにおける局所音速値を求める。具体的には、ホイヘンスの原理により、格子点XROIからの受信波WXと格子点A1,A2,…からの受信波を仮想的に合成した受信波WSUMが一致することを利用する。
格子点XROI及び格子点A1,A2,…の位置は、上述の送信焦点125の位置に基づき定められる。例えば、図10に示した各送信焦点125の配列において、図中の下から2行目の送信焦点125の1つを格子点XROIと定めた場合に、図中の下から1行目の送信焦点125の中から格子点A1,A2,…が定められる。また、図中の下から3行目の送信焦点125の1つを格子点XROIと定めた場合に、図中の下から2行目の送信焦点125の中から格子点A1,A2,…が定められる。そして、各格子点XROI,A1,A2,…(送信焦点125)について予め取得されてシネメモリ602に格納されているRFデータに基づき局所音速値の算出が行われる。
ここで、格子点XROIにおける局所音速値を求めるときの演算に使用する格子点A1,A2,…の範囲及び個数は予め決めておく。なお、局所音速値演算に使用する格子点の範囲が広いと局所音速値の誤差が大きくなり、狭いと仮想受信波との誤差が大きくなるため、格子点の範囲はこれらの兼ね合いで決める。
また、格子点A1,A2,…のX方向の間隔は、分解能と処理時間の兼ね合いで決定される。格子点A1,A2,…のX方向の間隔は一例で1mmから1cmである。なお、格子点A1,A2,…の間隔が広い場合には、合成波の演算(後述)が困難になるため、補間によって細かい格子点を生成すればよい。
図12に示すフローチャートにおいて、局所音速値算出部121は、上述のように格子点XROI、格子点A1,A2,…の設定を行った後(ステップS1)、各格子点(XROI、A1,A2,…)における最適音速値を判定する(ステップS2)。ここで、最適音速値とは、画像のコントラスト、シャープネスが最も高くなる音速値であり、各格子点における実際の局所音速値とは必ずしも一致しない。ステップS2における最適音速値の判定方法としては、例えば、画像のコントラスト、スキャン方向の空間周波数、分散などから判定する方法(例えば、特開平8−317926号公報)を適用することができる。
次いで、各格子点における最適音速値に基づいて、各格子点における局所音速値の判定が行われる(ステップS3)。以下、最適音速値に基づいて局所音速値を判定する方法を説明する。
最初に、図13に示すように、局所音速値算出部121は、格子点XROIにおける最適音速値に基づいて、格子点XROIを反射点とした時の仮想的な受信波(以下、仮想受信波という)WXの波形を算出する(ステップS4)。
次いで、局所音速値算出部121は、格子点XROIにおける仮定音速の初期値を設定する(ステップS5)。そして、局所音速値算出部121は、下記のように仮想的な合成受信波(以下、仮想合成受信波という)WSUMを算出する(ステップS6)。格子点XROIにおける局所音速値をVと仮定すると、格子点XROIから伝播した超音波が格子点A1,A2,…に到達するまでの時間はXROIA1/V,XROIA2/V,…となる。ここで、XROIA1,XROIA2,…は、それぞれ格子点A1,A2,…と格子点XROIとの間の距離である。格子点A1,A2,…における最適音速値はステップS2により既知のため、各格子点A1,A2,…からの受信波は予め求めることができる。従って、格子点A1,A2,…からそれぞれ遅延XROIA1/V,XROIA2/V,…で発した反射波(超音波エコー)を合成することにより、仮想合成受信波WSUMを求めることができる。
なお、実際には、素子データ(RF信号)上で上記処理を行うため、格子点XROIから格子点A1,A2,…に到達するまでの時間(それぞれT1,T2,…)は下記の式(1)により表される。ここで、XA1,XA2,…は、それぞれ格子点A1,A2,…と格子点Xとの間のスキャン方向(X方向)の距離である。また、Δtは格子点のY方向時間間隔である。
上記T1,T2,…に、格子点XROIと同一の走査線上の格子点Anから格子点XROIに到達するまでの時間(Δt/2)を足した遅延で各格子点A1,A2,…からの受信波を合成することにより、仮想合成受信波WSUMを求めることができる。
ここで、格子点をY方向に時間軸で等間隔(Δt)に設定する場合に、空間上での間隔は必ずしも等間隔にはならない。従って、各格子点に超音波が到達するまでの時間を計算するときに、式(1)においてΔt/2の代わりに補正したΔt/2を用いてもよい。ここで、補正したΔt/2は、例えば、格子点XROIと同音線の格子点Anに比べたA1,A2,…の深さ(Y方向の距離)の差をVで除算した値をΔt/2から加算・減算した値である。各格子点A1,A2,…の深さはそれより浅い格子点において局所音速が既知であることから求められる。
また、仮想合成受信波WSUMの算出は、実際に格子点A1,A2,…から遅延XROIA1/V,XROIA2/V,…で発した既定のパルス波(それぞれWA1,WA2,…)を重ね合わせることにより行う。
次に、局所音速値算出部121は、仮想受信波WXと仮想合成受信波WSUMとの誤差を算出する(ステップS7)。仮想受信波WXと仮想合成受信波WSUMとの誤差は、互いの相互相関をとる方法、仮想受信波WXに仮想合成受信波WSUMから得られる遅延を掛けて位相整合加算する方法、又は逆に仮想合成受信波WSUMに仮想受信波WXから得られる遅延を掛けて位相整合加算する方法により算出される。ここで、仮想受信波WXから遅延を得るには、格子点XROIを反射点とし、音速Vで伝播した超音波が各素子302に到着する時刻を遅延とすればよい。また、仮想合成受信波WSUMから遅延を得るには、隣り合う素子302間での合成受信波の位相差から等位相線を抽出し、その等位相線を遅延とするか、又は単に各素子302の合成受信波の最大(ピーク)位置の位相差を遅延としてもよい。また、各素子302からの合成受信波の相互相関ピーク位置を遅延としてもよい。位相整合加算時の誤差は、整合加算後の波形のpeak to peakとする方法、又は包絡線検波した後の振幅の最大値とする方法により求められる。
次いで、仮定音速が1ステップ変更された後(ステップS8)、ステップS6からS7が繰り返される。そして、全ての仮定音速の値での演算が終了すると、局所音速値算出部121は、格子点XROIにおける局所音速値を判定する(ステップS9)。ホイヘンスの原理を厳密に適用した場合に、上記のステップS6において求めた仮想合成受信波WSUMの波形は、格子点XROIにおける局所音速値をVと仮定した場合の仮想受信波(反射波)WXの波形と等しくなる。ステップS9では、仮想受信波WXと仮想合成受信波WSUMとの差が最小になる仮定音速の値を格子点XROIにおける局所音速値と判定する。これにより、着目領域112内の1つの格子点XROIにおける局所音速値が求められる。
以下同様にして、局所音速値算出部121は、上述の送信焦点125の位置に基づき定められる着目領域112内の他の複数の格子点XROIの局所音速値を算出する。
<音速マップ生成処理>
図8に戻って、音速マップ生成部122は、局所音速値算出部121が算出した着目領域112内の各格子点XROIの局所音速値に基づき、着目領域112内の局所音速の分布を示す音速マップ130(図16参照)を生成する。
図8に戻って、音速マップ生成部122は、局所音速値算出部121が算出した着目領域112内の各格子点XROIの局所音速値に基づき、着目領域112内の局所音速の分布を示す音速マップ130(図16参照)を生成する。
<音速マップ更新処理>
更新制御部123は、一定期間経過毎、例えば所定フレーム数(10フレームなど)のBモード画像データの生成が行われる毎に、送受信制御部120、ドプラ信号処理部701、音速測定領域決定部702、局所音速値算出部121、音速マップ生成部122を制御して、音速マップ130の更新処理を行う。具体的には、ドプラモード走査、カラードプラ情報またはCFM像データの生成と、音速測定領域119の決定と、音速測定領域119内のRFデータの取得及び局所音速値の算出と、音速マップ130の生成とが行われる。
更新制御部123は、一定期間経過毎、例えば所定フレーム数(10フレームなど)のBモード画像データの生成が行われる毎に、送受信制御部120、ドプラ信号処理部701、音速測定領域決定部702、局所音速値算出部121、音速マップ生成部122を制御して、音速マップ130の更新処理を行う。具体的には、ドプラモード走査、カラードプラ情報またはCFM像データの生成と、音速測定領域119の決定と、音速測定領域119内のRFデータの取得及び局所音速値の算出と、音速マップ130の生成とが行われる。
<第1実施形態の超音波診断装置の作用>
次に、図14に示すフローチャートを用いて上記構成の超音波診断装置10の作用について具体的に説明する。単独表示モードでの超音波診断が開始されると、CPU100は、送受信部400に制御信号を出力して、超音波ビームの被検体OBJへの送信、及び被検体OBJからの超音波エコーの受信を開始させる。各素子302で検出された超音波検出信号は送受信部400にてデジタルなRFデータに変換される。
次に、図14に示すフローチャートを用いて上記構成の超音波診断装置10の作用について具体的に説明する。単独表示モードでの超音波診断が開始されると、CPU100は、送受信部400に制御信号を出力して、超音波ビームの被検体OBJへの送信、及び被検体OBJからの超音波エコーの受信を開始させる。各素子302で検出された超音波検出信号は送受信部400にてデジタルなRFデータに変換される。
送受信部400にて生成されたRFデータは、シネメモリ602に格納されるとともに、画像信号生成部500に出力される。そして、画像信号生成部500にてRFデータに基づきBモード画像データが生成されてアナログの画像信号に変換された後に表示制御部105に出力される。この画像信号に基づき、表示制御部105は、表示部104にBモード画像110を表示させる(ステップS10)。
操作卓202の表示モード切り替えボタンが押下されて、単独表示モードから音速マップ表示モードに切り替えられると(ステップS11)、CPU100は、表示部104に着目領域112の設定を促すメッセージ等を表示させる。医師は、表示部104に表示されるBモード画像110を見ながら、ポインティングデバイス204を操作してBモード画像110上に着目領域112を設定する(ステップS12)。
着目領域112の設定がなされると、送受信制御部120は、送受信部400を制御してドプラモード走査に対応した超音波の送受信を各素子302に実行させる。この超音波の送受信により送受信部400にてRFデータが生成され(ステップS13)、このRFデータがドプラ信号処理部701に入力される。なお、ドプラモード走査に対応した超音波の送受信が行われた後は、Bモード走査に対応した超音波の送受信、Bモード画像データの生成・出力が繰り返し実行される。
ドプラ信号処理部701は、図3から図5に示したように、ドプラモード走査時に得られたRFデータからカラードプラ情報またはCFM像データを生成する(ステップS14)。このカラードプラ情報またはCFM像データは、音速測定領域決定部702に入力される。
音速測定領域決定部702は、ドプラ信号処理部701からのカラードプラ情報またはCFM像データに基づき、先に設定された着目領域112の中から移動体領域117を除いた音速測定領域119を決定する(図9(a)〜図9(c)参照、ステップS15)。音速測定領域決定部702の決定結果はCPU100に入力される。
音速測定領域決定部702からの音速測定領域119の決定結果を受けて、送受信制御部120は、図9及び図10に示したように、音速測定領域119内に送信焦点125を設定した後、送受信部400を制御して超音波の送受信(送信フォーカス、受信フォーカス)を各素子302に実行させる。これにより、各送信焦点125のRFデータが取得(生成・格納)される(ステップS16)。
各送信焦点125のRFデータの取得後、局所音速値算出部121は、図11から図13に示したように、ホイヘンスの原理を利用して、音速測定領域119内の各格子点XROIの局所音速値を算出する。次いで、音速マップ生成部122は、局所音速値算出部121による局所音速値の算出結果に基づき、音速マップ130を生成してDSC504に出力する(ステップS17)。
音速マップ130は、DSC504にてテレビジョン信号の走査方式の画像データに変換され、画像処理部506にて画像処理が施されて画像メモリ508に格納された後、D/A変換器510にてアナログの画像信号に変換されて表示制御部105に出力される。
図15に示すように、表示制御部105は、音速マップ130の画像信号と、Bモード走査で得られたBモード画像の画像信号とに基づき、Bモード画像110に音速マップ130を重畳して表示部104に表示させる(ステップS18)。音速マップ130では、音速測定領域119内が局所音速の値に応じて色分けして表示される。
図14に戻って、ステップS10で説明したように新たなBモード画像データが生成されて画像メモリ508に格納される。先の音速マップ130の生成から一定期間が経過していない場合、具体的には所定フレーム数(10フレームなど)のBモード画像データが生成されていない場合(ステップS19でNO)には、更新制御部123による音速マップ130の更新処理は実行されない。このため、新たに生成されたBモード画像上には、先に生成された音速マップ130が重畳して表示される。
一方、先の音速マップ130の生成から一定期間が経過した場合、具体的には所定フレーム数(10フレームなど)のBモード画像データが生成された場合(ステップS19でYES)に、更新制御部123は、音速マップ130の更新処理を実行させる。これにより、上述のステップS13からステップS18までの処理が繰り返し実行される。これにより、新たな音速マップ130がBモード画像110上に重畳表示される。
以下、音速マップ表示モードが終了(ステップS20でYES)するまで、上述の各処理が繰り返し実行される。
なお、図16(a)に示すように、上記第1実施形態では、図中「S」で示した音速測定領域119のRFデータ取得のための超音波の送受信の直前に、図中「D」に示したドプラモード走査に対応した超音波の送受信を行っている。なお、図中「B」は、Bモード走査に対応した超音波の送受信である。この際に、例えば、図16(b)に示すように、ドプラモード走査に対応した超音波の送受信と、音速測定領域119のRFデータ取得のための超音波の送受信との間に、Bモード走査に対応した超音波の送受信を行ってもよい。
<作用効果>
以上のように、本発明では、着目領域112内の音速測定領域119内に送信焦点125を設定し、各送信焦点125に対して超音波の送受信を行って得られたRFデータに基づき音速測定領域119内の局所音速値を求めて音速マップ130を生成しているので、着目領域112内の全領域のRFデータを取得する必要がなくなる。その結果、音速マップ130の生成の際における演算処理を従来よりも低減させることができるので、高性能な演算処理装置(CPU)を用いることなく、音速マップ130の生成並びに表示に要する時間を短縮することができる。さらに、血管(血流)や動きの大きい組織等の各種移動体の影響を受けない音速マップ130を生成することができる。
以上のように、本発明では、着目領域112内の音速測定領域119内に送信焦点125を設定し、各送信焦点125に対して超音波の送受信を行って得られたRFデータに基づき音速測定領域119内の局所音速値を求めて音速マップ130を生成しているので、着目領域112内の全領域のRFデータを取得する必要がなくなる。その結果、音速マップ130の生成の際における演算処理を従来よりも低減させることができるので、高性能な演算処理装置(CPU)を用いることなく、音速マップ130の生成並びに表示に要する時間を短縮することができる。さらに、血管(血流)や動きの大きい組織等の各種移動体の影響を受けない音速マップ130を生成することができる。
[第2実施形態の超音波診断装置の構成]
次に、図17を用いて本発明の第2実施形態の超音波診断装置20について説明を行う。上記第1実施形態では、音速測定領域119内に送信焦点125を設定して取得したRFデータに基づき、音速測定領域119内の局所音速値を算出している。これに対して、第2実施形態の超音波診断装置20では、着目領域112内に送信焦点125を設定して取得したRFデータに基づき、音速測定領域119内の局所音速値を算出する。
次に、図17を用いて本発明の第2実施形態の超音波診断装置20について説明を行う。上記第1実施形態では、音速測定領域119内に送信焦点125を設定して取得したRFデータに基づき、音速測定領域119内の局所音速値を算出している。これに対して、第2実施形態の超音波診断装置20では、着目領域112内に送信焦点125を設定して取得したRFデータに基づき、音速測定領域119内の局所音速値を算出する。
なお、超音波診断装置20は、CPU100が第1実施形態とは異なる送受信制御部120a、局所音速値算出部121aとして機能する点を除けば、第1実施形態と基本的に同じ構成であるので、第1実施形態と機能・構成上同一のものについては、同一符号を付してその説明は省略する。
送受信制御部120aは、着目領域112内の全領域に送信焦点125を設定して超音波の送受信を送受信部400(各素子302)に実行させることで、着目領域112内のRFデータを取得する。また、局所音速値算出部121aは、着目領域112内のRFデータの中から音速測定領域119内の送信焦点125に対応するRFデータを選択して、局所音速値を算出する。
<第2実施形態の超音波診断装置の作用>
図18に示すフローチャート及び図19及び図20に示す説明図を用いて、上記構成の超音波診断装置20の作用について具体的に説明する。なお、音速測定領域119の決定がなされるまでの処理(ステップS10からステップS15まで)は、第1実施形態と同じであるので説明は省略する。
図18に示すフローチャート及び図19及び図20に示す説明図を用いて、上記構成の超音波診断装置20の作用について具体的に説明する。なお、音速測定領域119の決定がなされるまでの処理(ステップS10からステップS15まで)は、第1実施形態と同じであるので説明は省略する。
図19に示すように、音速測定領域119の決定後、送受信制御部120aは、着目領域112内の全領域に送信焦点125を設定した後、送受信部400を制御して超音波の送受信(送信フォーカス、受信フォーカス)を実行させる。これにより、着目領域112の全領域内の各送信焦点125のRFデータが取得(生成・格納)される(ステップS21)。なお、本実施形態では、送信焦点125を格子状に設定しているが、送信焦点125の設定パターンは適宜変更してもよい。また、ステップS21は、ステップS13からステップS15の間で適宜実施してもよい。
送信焦点125の設定後、局所音速値算出部121aは、音速測定領域決定部702からCPU100に入力される音速測定領域119の決定結果に基づき、着目領域112内の音速測定領域119を判別する。
次いで、図20に示すように、局所音速値算出部121aは、着目領域112内の各送信焦点125に対応するRFデータの中から、音速測定領域119内の送信焦点125(図中、黒丸で表示)に対応するRFデータを選択した後、第1実施形態と同様にして局所音速値を算出する。これにより、音速測定領域119内の局所音速値が得られる。なお、図20中の符号「125’」は音速測定領域119外の送信焦点(図中で白丸で表示)である。音速マップ生成部122は、局所音速値算出部121aによる局所音速値の算出結果に基づき、音速マップ130を生成する(ステップS22)。
以下、第1実施形態と同様に、音速マップ130がBモード画像110上に重畳表示される。また、先の音速マップ130の生成から一定期間が経過した場合(例えば10フレーム毎)には、上述のステップS13〜S15,S21,S22,S18までの処理が繰り返し実行されることで、音速マップ130が更新される。
超音波診断装置20では、着目領域112の全領域に送信焦点125を設定してRFデータの取得を行うので、音速マップ130の生成のための演算処理や生成時間が増加するものの、血管(血流)等の影響を受けない音速マップ130を生成することができる。
[その他]
上記各実施形態の局所音速値算出部121,121aは、特許文献1に開示されたホイヘンスの原理を利用した局所音速算出法を用いて音速測定領域119内の局所音速値を求めているが、例えば特許文献1に開示された屈折率モデル計算を利用した局所音速算出法を用いて音速測定領域119内の局所音速値を求めてもよく、局所音速算出法は特に限定はされない。
上記各実施形態の局所音速値算出部121,121aは、特許文献1に開示されたホイヘンスの原理を利用した局所音速算出法を用いて音速測定領域119内の局所音速値を求めているが、例えば特許文献1に開示された屈折率モデル計算を利用した局所音速算出法を用いて音速測定領域119内の局所音速値を求めてもよく、局所音速算出法は特に限定はされない。
上記各実施形態では、着目領域112を矩形状に設定しているが、例えば円形状等の各種形状に着目領域を設定してもよい。
上記各実施形態では、所定フレーム数毎(例えば10フレーム)毎に音速マップ130の更新処理を行っているが、例えば先の音速マップ130の生成時から一定時間経過後などの各種一定期間が経過した後に音速マップ130の更新処理を行ってもよい。
上記各実施形態では、ドプラ信号処理部701によりカラードプラ情報やCFM像データを生成しているが、ドプラ信号処理部701の構成は、カラードプラ情報やCFM像データを生成可能であれば図3に示した構成に限定されるものではなく適宜変更してもよい。
上記各実施形態では、Bモード走査とドプラモード走査とで各素子302の駆動(超音波の送受信)を異ならせているが、カラードプラ情報やCFM像データを取得可能であれば各素子302の駆動を同じにしてもよい。この場合には、Bモード画像データの生成と、カラードプラ情報やCFM像データの生成とを同時に行うことができる。
10,20…超音波診断装置,100…CPU,110…Bモード画像,112…着目領域,119…音速測定領域,120…送受信制御部,121…局所音速値算出部,122…音速マップ生成部,123…更新制御部,125…送信焦点,130…音速マップ,200…操作入力部,300…超音波探触子,400…送受信部,701…ドプラ信号処理部,702…音速測定領域決定部
Claims (8)
- 超音波を被検体に送信し、かつ前記被検体によって反射された超音波を受信して超音波検出信号を出力する複数の素子を含む超音波探触子と、
前記複数の素子をそれぞれ駆動する複数の駆動信号を前記複数の素子に供給し、かつ当該複数の素子から出力された前記超音波検出信号に基づき受信データを生成する送受信部であって、前記複数の素子に供給する前記複数の駆動信号の遅延量を調節することで超音波ビームを形成する送受信部と、
前記被検体内に着目領域を設定する着目領域設定部と、
前記送受信部を制御して、前記着目領域内の移動体の情報を得るための超音波の送受信を前記複数の素子に実行させる第1の送受信制御部と、
前記移動体の情報を得るための超音波の送受信がなされたときに、前記移動体によりドプラ偏移されて反射された超音波の前記超音波検出信号に基づいて前記送受信部が生成した前記受信データから、前記移動体の動きに伴うドプラ信号を抽出して、この抽出結果に基づき前記着目領域内の前記移動体を含む移動体領域を判別する移動体領域判別部と、
前記移動体領域判別部の判別結果に基づき、前記着目領域の中から前記移動体領域を除いた領域を、音速の分布を測定する音速測定領域として決定する音速測定領域決定部と、
前記音速測定領域内に前記超音波ビームの送信焦点を複数設定し、各前記送信焦点に対する前記超音波ビームの送受信を前記送受信部に行わせることで、前記音速測定領域の音速測定用の受信データを前記送受信部に生成させる第2の送受信制御部と、
前記音速測定用の受信データに基づき前記音速測定領域の音速の分布を示す音速マップを生成する音速マップ生成部と、
を備える超音波診断装置。 - 前記送受信部を制御して、前記被検体内の断層像情報を得るための超音波の送受信を前記複数の素子に実行させる第3の送受信制御部と、
前記断層像情報を得るための超音波の送受信がなされたときに前記送受信部が生成した前記受信データに基づきBモード画像を生成する画像生成部と、
前記Bモード画像を表示する表示部と、を備え、
前記着目領域設定部は、前記Bモード画像上に前記着目領域を設定する請求項1記載の超音波診断装置。 - 前記第1の送受信制御部による前記送受信部の制御と、前記移動体領域判別部による前記移動体領域の判別と、前記音速測定領域決定部による前記音速測定領域の決定と、前記第2の送受信制御部による前記送受信部の制御と、前記音速マップ生成部による前記音速マップの生成と、を含む音速マップ更新処理を実行させる更新制御部を備える請求項2記載の超音波診断装置。
- 前記更新制御部は、予め定められたフレーム数の前記Bモード画像が前記画像生成部により生成される毎に、前記音速マップ更新処理を実行させる請求項3記載の超音波診断装置。
- 前記移動体には、前記着目領域内の血管を流れる血液が含まれており、
前記移動体領域判別部は、前記ドプラ信号の抽出結果に基づき、前記着目領域内での前記血液の流れを示す血流情報を生成して、前記血流情報に基づき前記移動体領域を判別する請求項1から4のいずれか1項記載の超音波診断装置。 - 前記移動体には、前記着目領域内で移動する組織が含まれており、
前記移動体領域判別部は、前記ドプラ信号の抽出結果に基づき、前記着目領域内での前記組織の移動を示す組織移動情報を生成して、前記組織移動情報に基づき前記移動体領域を判別する請求項1から5のいずれか1項記載の超音波診断装置。 - 前記移動体領域判別部は、前記送受信部が生成した前記受信データから、前記組織の中で動きの大きさが予め定めたしきい値以上となる組織に対応する前記ドプラ信号を抽出して、当該ドプラ信号に基づき前記組織移動情報を生成する請求項6記載の超音波診断装置。
- 前記表示部は、前記Bモード画像上に前記音速マップを重畳表示する請求項2から4のいずれか1項記載の超音波診断装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2012119579A JP2013244162A (ja) | 2012-05-25 | 2012-05-25 | 超音波診断装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2012119579A JP2013244162A (ja) | 2012-05-25 | 2012-05-25 | 超音波診断装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2013244162A true JP2013244162A (ja) | 2013-12-09 |
Family
ID=49844432
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2012119579A Pending JP2013244162A (ja) | 2012-05-25 | 2012-05-25 | 超音波診断装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2013244162A (ja) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019026115A1 (ja) * | 2017-07-31 | 2019-02-07 | 本多電子株式会社 | 超音波画像表示装置及び方法、並びにプログラムを格納した記録媒体 |
| KR20200110541A (ko) * | 2019-03-14 | 2020-09-24 | 삼성메디슨 주식회사 | 초음파 영상 장치 및 그 제어방법 |
| US20210038195A1 (en) * | 2018-04-28 | 2021-02-11 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Ultrasound transient elasticity measurement device and method |
-
2012
- 2012-05-25 JP JP2012119579A patent/JP2013244162A/ja active Pending
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019026115A1 (ja) * | 2017-07-31 | 2019-02-07 | 本多電子株式会社 | 超音波画像表示装置及び方法、並びにプログラムを格納した記録媒体 |
| JP6484781B1 (ja) * | 2017-07-31 | 2019-03-20 | 本多電子株式会社 | 超音波画像表示装置 |
| US20210038195A1 (en) * | 2018-04-28 | 2021-02-11 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Ultrasound transient elasticity measurement device and method |
| US12226261B2 (en) * | 2018-04-28 | 2025-02-18 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Ultrasound transient elasticity measurement device and method |
| KR20200110541A (ko) * | 2019-03-14 | 2020-09-24 | 삼성메디슨 주식회사 | 초음파 영상 장치 및 그 제어방법 |
| KR102749835B1 (ko) | 2019-03-14 | 2025-01-08 | 삼성메디슨 주식회사 | 초음파 영상 장치 및 그 제어방법 |
| US12383234B2 (en) | 2019-03-14 | 2025-08-12 | Samsung Medison Co., Ltd. | Ultrasonic imaging apparatus and method for controlling the same |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5597734B2 (ja) | 超音波イメージング装置及び超音波イメージングプログラム | |
| JP5389722B2 (ja) | 超音波診断装置及びその作動方法 | |
| JP5808325B2 (ja) | 超音波診断装置及び超音波診断装置の作動方法 | |
| JP6288996B2 (ja) | 超音波診断装置及び超音波イメージングプログラム | |
| JP5798117B2 (ja) | 超音波診断装置及び超音波診断装置の作動方法 | |
| JP5588924B2 (ja) | 超音波診断装置 | |
| JP2008168016A (ja) | 超音波診断装置、imt計測方法及びimt計測プログラム | |
| JP2012217611A (ja) | 超音波診断装置および超音波画像生成方法 | |
| JP5948411B2 (ja) | 超音波信号処理装置および超音波信号処理方法 | |
| JP5623160B2 (ja) | 超音波診断装置及びその作動方法 | |
| JP2011072566A (ja) | 超音波診断装置及びその信号処理方法 | |
| JP5623157B2 (ja) | 超音波診断装置及びその作動方法 | |
| JP5247330B2 (ja) | 超音波信号処理装置及び超音波信号処理方法 | |
| JP2010234013A (ja) | 超音波診断装置及び超音波診断方法 | |
| JP2013244159A (ja) | 超音波診断装置及び音速推定方法 | |
| JP2012071037A (ja) | 超音波画像診断装置、超音波画像生成方法およびプログラム | |
| JP2013244162A (ja) | 超音波診断装置 | |
| WO2013176112A1 (ja) | 超音波画像生成方法および超音波画像診断装置 | |
| JP5851345B2 (ja) | 超音波診断装置およびデータ処理方法 | |
| CN104363837A (zh) | 超声波图像诊断装置 | |
| JP2013102960A (ja) | 超音波診断装置及び方法 | |
| JP2013244161A (ja) | 超音波診断装置 | |
| JP2008142130A (ja) | 超音波診断装置およびその制御処理プログラム | |
| JP2026020093A (ja) | 超音波診断装置 | |
| JP2012249925A (ja) | 超音波画像生成方法および超音波画像診断装置 |
