(57)【要約】【課題】超音波探触子に正対していない体内組織を含む超音波画像を好適に生成することができる超音波診断装置および画像合成方法を提供する。【解決手段】複数の異なるステアリング角度で超音波を送信し、被検体の反射面で反射された超音波を受信する超音波探触子２が受信した超音波に基づいて生成した受信信号に基づいて、複数の異なるステアリング角度に対応する受信信号からそれぞれ異なる複数の超音波画像データを生成する画像生成部１４と、複数の異なるステアリング角度と、複数の異なるステアリング角度に対応する複数の超音波画像データと、に基づいて回帰分析を行い、回帰分析の結果に基づいて所定の重み付け値である回帰推定値を生成する回帰推定値生成部１６と、複数の超音波画像データに対して、回帰推定値に基づく重み付けを行って合成を行い、合成画像データを生成する画像合成部１７と、を備える。【選択図】図２
2018051040.tif 000002
【特許請求の範囲】
【請求項１】
複数の異なるステアリング角度で超音波を送信し、被検体の反射面で反射された反射超音波を受信する超音波探触子が前記受信した反射超音波に基づいて生成した受信信号に基づいて、前記複数の異なるステアリング角度に対応する前記受信信号からそれぞれ異なる複数の超音波画像データを生成する画像生成部と、
前記複数の異なるステアリング角度と、前記複数の異なるステアリング角度に対応する複数の超音波画像データと、に基づいて回帰分析を行い、前記回帰分析の結果に基づいて所定の重み付け値である回帰推定値を生成する回帰推定値生成部と、
前記複数の超音波画像データに対して、前記回帰推定値に基づく重み付けを行って合成を行い、合成画像データを生成する画像合成部と、
を備える超音波診断装置。
【請求項２】
前記回帰推定値生成部は、前記回帰分析を用いて、前記超音波画像データにおいて、前記超音波の送信方向に対して反射面が傾いている前記被検体内の組織に対応する領域に対しては重みが相対的に大きくなるように、それ以外の領域に対しては重みが相対的に小さくなるように、前記回帰推定値を生成する、
請求項１に記載の超音波診断装置。
【請求項３】
前記回帰推定値は、前記回帰分析により得られた回帰式の係数、残差平方和、決定係数のうちの少なくともいずれかに基づいて生成される値である、
請求項１または２に記載の超音波診断装置。
【請求項４】
前記回帰推定値は、単回帰分析により得られた１次回帰式の傾きである、
請求項１または２に記載の超音波診断装置。
【請求項５】
前記画像合成部は、前記複数のステアリング角度に対応する複数の超音波データの前記所定の領域毎に、画素値を第１の関数に入力して取得した第１の出力値と、前記画素値を第２の関数に入力して取得した第２の出力値と、前記回帰推定値と、を用いてαブレンドにより合成画像データを生成する、
請求項１から４のいずれか一項に記載の超音波診断装置。
【請求項６】
前記画像合成部は、前記複数のステアリング角度に対応する複数の超音波データの前記所定の領域毎に、画素値を第１の関数に入力して取得した第１の出力値と、前記画素値を第２の関数に入力して取得した第２の出力値と、前記回帰推定値と、所定の固定値である固定係数と、を用いて合成画像データを生成する、
請求項１から４のいずれか一項に記載の超音波診断装置。
【請求項７】
前記第１の関数および第２の関数は、入力値に対して、平均値、最大値、最小値、中央値、分散、標準偏差、のうちのいずれかの値を取得する関数である、
請求項５または６に記載の超音波診断装置。
【請求項８】
超音波診断装置の画像合成方法であって、
複数の異なるステアリング角度で送信された超音波が被検体の反射面で反射された反射超音波を受信する超音波探触子が前記受信した反射超音波に基づいて生成した受信信号に基づいて、前記複数の異なるステアリング角度に対応する前記受信信号からそれぞれ異なる複数の超音波画像データを生成し、
前記複数の異なるステアリング角度と、前記複数の異なるステアリング角度に対応する複数の超音波画像データと、に基づいて回帰分析を行い、
前記回帰分析の結果に基づいて所定の重み付け値である回帰推定値を生成し、
前記複数の超音波画像データに対して、前記回帰推定値に基づく重み付けを行って合成を行い、合成画像データを生成する、
画像合成方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波診断装置および画像合成方法に関する。
【背景技術】
【０００２】
従来、多数の振動子を配列して備える超音波探触子を有し、生体等の被検体に対して超音波の送受信を行い、受信した超音波から得られた信号に基づいて超音波画像データを生成し、これに基づく超音波画像を画像表示装置に表示する超音波診断装置が知られている。このような装置による超音波画像診断は、超音波探触子を被検体の体表に当てるだけの簡単な操作で心臓の拍動や胎児の動き等の様子がリアルタイムで得られ、かつ非侵襲で安全性が高いため、繰り返して実施することができる。
【０００３】
しかしながら、このような超音波診断装置によって得られる画像には、被検体内の組織に関する情報以外にも、各種のノイズや、超音波探触子で受信した超音波から得られた受信信号の干渉現象により発生するスペックルが存在し、これらは被検体内の組織の境界の位置や形状を正確に把握する場合においてしばしば妨げとなっている。
【０００４】
近年では、このようなノイズやスペックルを低減する処理方法として、例えば、空間コンパウンド法を用いた超音波診断装置が普及している。空間コンパウンド法は、被検体の同一部位に対して同時期に複数の異なる方向で超音波の送受信を行い、取得された複数の超音波画像データの平均的重畳を行う方法である。これにより、ノイズやスペックルは、例えば、Ｎ枚の超音波画像データが得られた場合には、これらを合成した合成画像データにおいて、Ｎの平方根で低減される。
【０００５】
また、この空間コンパウンド法によれば、異方性部位の抽出性能を向上させることができる。異方性部位とは、超音波が当たった際の散乱、反射等の受信信号強度が角度により異なる部位であり、具体的には、例えば被検体内の骨格筋における腱や靭帯のように、繊維状で反射強度は骨表面ほど強くないが鏡面反射特性を示す軟部組織における部位である。
【０００６】
このような空間コンパウンド法を用いた超音波診断装置であって、診断検査のタイプに応じて選択された制御信号によって、複数方向からの反射信号によって得られた超音波画像の各画素値の平均値、最大値、最小値、中央値等を合成することで、より品質が高い超音波画像を生成することができる技術が、例えば特許文献１に開示されている。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特許第３９３５４３３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
探触子（超音波探触子）が送信する超音波が被検体内の体内組織に反射されるとき、一般に、超音波探触子に正対した（超音波の送信方向と反射面が垂直な）体内組織からの反射信号強度は、超音波探触子と正対していない（超音波の送信方向と反射面が垂直ではない）体内組織からの反射強度よりも強くなる。このため、超音波の反射信号に基づいて生成された超音波診断画像は、超音波探触子に正対していない体内組織の方が、超音波探触子に正対した体内組織と比較して不明瞭に描写されてしまう場合がある。
【０００９】
例えば特許文献１に開示された技術において、画素値の最大値を合成することにより画素値を増幅すれば、超音波探触子と正対していない体内組織からの反射信号を増幅することができるため、超音波探触子に正対していない体内組織を比較的明瞭に描写することができる。しかしながら、超音波探触子が受信した他の受信信号、すなわち、例えば超音波探触子と正対した体内組織からの反射信号や、スペックルノイズを含む信号等も増幅されるため、不自然な画像となってしまう恐れがある。
【００１０】
本発明は、超音波探触子に正対していない体内組織を含む超音波画像を好適に生成することができる超音波診断装置および画像合成方法を提供することを目的とする。
【課題を解決するための手段】
【００１１】
本発明の超音波診断装置は、複数の異なるステアリング角度で超音波を送信し、被検体の反射面で反射された反射超音波を受信する超音波探触子が前記受信した反射超音波に基づいて生成した受信信号に基づいて、前記複数の異なるステアリング角度に対応する前記受信信号からそれぞれ異なる複数の超音波画像データを生成する画像生成部と、前記複数の異なるステアリング角度と、前記複数の異なるステアリング角度に対応する複数の超音波画像データと、に基づいて回帰分析を行い、前記回帰分析の結果に基づいて所定の重み付け値である回帰推定値を生成する回帰推定値生成部と、前記複数の超音波画像データに対して、前記回帰推定値に基づく重み付けを行って合成を行い、合成画像データを生成する画像合成部と、を備える。
【００１２】
本発明の画像合成方法は、超音波診断装置の画像合成方法であって、複数の異なるステアリング角度で送信された超音波が被検体の反射面で反射された反射超音波を受信する超音波探触子が前記受信した反射超音波に基づいて生成した受信信号に基づいて、前記複数の異なるステアリング角度に対応する前記受信信号からそれぞれ異なる複数の超音波画像データを生成し、前記複数の異なるステアリング角度と、前記複数の異なるステアリング角度に対応する複数の超音波画像データと、に基づいて回帰分析を行い、前記回帰分析の結果に基づいて所定の重み付け値である回帰推定値を生成し、前記複数の超音波画像データに対して、前記回帰推定値に基づく重み付けを行って合成を行い、合成画像データを生成する。
【発明の効果】
【００１３】
本発明によれば、超音波探触子に正対していない体内組織を含む超音波画像を好適に生成することができる超音波診断装置および画像合成方法を提供することができる。
【図面の簡単な説明】
【００１４】
【図１】超音波診断装置の構成の一例を示す図
【図２】超音波診断装置本体の構成の一例を示すブロック図
【図３】超音波探触子から送信される送信ビームの角度と操作領域との関係を示す図
【図４Ａ】超音波探触子に対して正対したターゲットに対してステア角度０°で超音波ビームを送信する様子を例示した図
【図４Ｂ】超音波探触子に対して正対したターゲットに対してステア角度θ（θ≠０）で超音波ビームを送信する様子を例示した図
【図５Ａ】超音波探触子２に対して正対したターゲットに対してステア角度０°で超音波ビームを送信する様子を例示した図
【図５Ｂ】超音波探触子２に対して正対したターゲットに対してステア角度θ（θ≠０）で超音波ビームを送信する様子を例示した図
【図６Ａ】正対ターゲットからの反射信号強度とステア角度との関係を示す図
【図６Ｂ】非正対ターゲットからの反射信号強度とステア角度との関係を示す図
【図６Ｃ】スペックルノイズを含む反射信号強度とステア角度との関係を示す図
【図７Ａ】正対ターゲットの場合の１次回帰式を例示した図
【図７Ｂ】非正対ターゲットの場合の１次回帰式を例示した図
【図７Ｃ】スペックルノイズの場合の１次回帰式を例示した図
【図８Ａ】正対ターゲットの場合の２次回帰式を例示した図
【図８Ｂ】非正対ターゲットの場合の２次回帰式を例示した図
【図８Ｃ】スペックルノイズの場合の２次回帰式を例示した図
【発明を実施するための形態】
【００１５】
以下、本発明の実施の形態に係る超音波診断装置について、図面を参照して説明する。ただし、発明の範囲は図示した例に限定されない。なお、以下の説明において、同一の機能および構成を有するものについては、同一の符号を付し、その説明を省略する。
【００１６】
図１は、超音波診断装置１００の構成の一例を示す図である。超音波診断装置１００は、図１に示すように、超音波診断装置本体１と、超音波探触子２と、を有する。超音波探触子２は、図示しない生体等の被検体内に対して超音波（送信超音波）を送信するとともに、この被検体内で反射した超音波の反射波（反射超音波：エコー）を受信する。
【００１７】
超音波診断装置本体１は、超音波探触子２とケーブル３を介して接続され、超音波探触子２に電気信号の駆動信号を送信することによって超音波探触子２に被検体に対して送信超音波を送信させる。そして、超音波探触子２にて受信した被検体内からの反射超音波に応じて超音波探触子２で生成された電気信号である受信信号に基づいて被検体内の内部状態を超音波画像として画像化する。
【００１８】
超音波探触子２は、複数の圧電素子からなる振動子２ａ（図２参照）を有し、この振動子２ａは、例えば、方位方向（走査方向）に一次元アレイ状に複数配列されている。振動子２ａの個数は、任意に設定することができる。
【００１９】
図２は、超音波診断装置本体１の構成の一例を示すブロック図である。超音波診断装置本体１は、図２に示すように、例えば、操作入力部１１と、送信部１２と、受信部１３と、画像生成部１４と、ＤＳＣ（Digital Scan Converter）１５と、回帰推定値生成部１６と、画像合成部１７と、表示部１８と、制御部１９と、を有する。
【００２０】
操作入力部１１は、例えば、診断開始を指示するコマンド、被検体の個人情報等のデータ、および、超音波画像を表示部１８に表示するための各種パラメータの入力などを行うための操作デバイスであり、具体的には、例えば各種スイッチ、ボタン、トラックボール、マウス、キーボード等である。操作入力部１１は、入力された操作に基づく操作信号を制御部１９に出力する。
【００２１】
送信部１２は、制御部１９の制御に従って、超音波探触子２にケーブル３を介して電気信号である駆動信号を供給し、超音波探触子２に送信超音波を発生させる回路である。送信部１２は、例えば、図示しないクロック発生回路、遅延回路、パルス発生回路を有する。クロック発生回路は、駆動信号の送信タイミングや送信周波数を決定するクロック信号を発生させる回路である。遅延回路は、振動子２ａ毎に対応した個別経路毎に遅延時間を設定し、設定された遅延時間だけ駆動信号の送信を遅延させて送信超音波によって構成される送信ビームの集束（送信ビームフォーミング）や、送信ビームの角度の設定（ステアリング）を行うための回路である。パルス発生回路は、所定の周期で駆動信号としてのパルス信号を発生させるための回路である。
【００２２】
上述のように構成された送信部１２は、例えば、超音波探触子２に配列された複数（例えば、百数十個〜二百数十個）の振動子２ａのうちの連続する一部（例えば、数十個）を駆動して送信超音波を発生させる。そして、送信部１２は、送信超音波を発生させる毎に駆動する振動子２ａを方位方向にずらすことで走査（スキャン）を行う。また、送信部１２は、送信ビームの角度を変更しながら走査を行うことで、角度の異なる複数の反射信号を受信することができる。
【００２３】
受信部１３は、制御部１９の制御に従って、超音波探触子２からケーブル３を介して電気信号である受信信号を受信する回路である。受信部１３は、例えば、増幅器、Ａ／Ｄ変換回路、整相加算回路を有する。増幅器は、受信信号を、振動子２ａ毎に対応した個別経路毎に、予め設定された増幅率で増幅させるための回路である。Ａ／Ｄ変換回路は、増幅された受信信号をアナログ／デジタル変換（Ａ／Ｄ変換）するための回路である。整相加算回路は、Ａ／Ｄ変換された受信信号に対して、振動子２ａ毎に対応した個別経路毎に遅延時間を与えて時相を整え、これらを加算（整相加算）して音線データを生成するための回路である。すなわち、整相加算回路は、振動子２ａ毎の受信信号に対して受信ビームフォーミングを行って音線データを生成する。
【００２４】
画像生成部１４は、制御部１９の制御に従って、受信部１３から入力された音線データに対して包絡線検波処理や対数圧縮などを実施し、ダイナミックレンジやゲインの調整を行って輝度変換することにより、Ｂモードの超音波画像データを生成する。Ｂモードの超音波画像データは、受信信号の強さを輝度によって表したものである。なお、本実施の形態において、画像生成部１４は、Ｂモード画像データの他、Ａモード画像（振幅画像）データおよびＭモード画像（動き画像）データや、ドップラ法による超音波画像データを生成できるようにしてもよい。
【００２５】
なお、本実施の形態において、上記したように超音波探触子２が角度をずらして走査を行った場合、画像生成部１４は、角度の異なる複数の反射信号に基づいて、複数の超音波画像データを生成する。このように生成された複数の超音波画像データは、走査領域の一部または全部がそれぞれ重複している。
【００２６】
ＤＳＣ１５は、制御部１９の制御に従って、画像生成部１４が生成した複数の超音波画像データに対して走査周波数変換を行い、回帰推定値生成部１６に対して出力する。
【００２７】
回帰推定値生成部１６は、制御部１９の制御に従って、ＤＳＣ１５が出力した複数の超音波画像データの画素毎、あるいは所定の領域毎に、所定の回帰分析を行って回帰推定値を生成する回帰分析処理を行い、回帰推定値を算出して画像合成部１７に出力する。なお、回帰推定値とは、回帰推定値生成部１６が行う回帰分析により得られるパラメータにより算出された値である。回帰推定値生成部１６が行う回帰分析処理の詳細は後述する。
【００２８】
画像合成部１７は、制御部１９の制御に従って、回帰推定値生成部１６が算出した回帰推定値と、ＤＳＣ１５が出力した複数の超音波画像データと、に基づいて、走査領域の重複した部分で合成して合成画像データを生成する画像合成処理を行って合成画像データを生成する。画像合成部１７が行う画像合成処理の詳細は後述する。
【００２９】
表示部１８は、制御部１９の制御に従って、ＬＥＤ（Light-Emitting Diode）、ＬＣＤ（Liquid Crystal Display）、ＣＲＴ（Cathode-Ray Tube）ディスプレイ、有機ＥＬ（Electronic Luminescence）ディスプレイ、無機ＥＬディスプレイ、およびプラズマディスプレイ等の表示装置である。表示部１８は、画像合成部１７から出力された合成画像データに基づいて超音波画像の表示を行う。
【００３０】
制御部１９は、例えば、ＣＰＵ（Central Processing Unit）、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）を有し、ＲＯＭに記憶されているシステムプログラム等の各種処理プログラムを読み出してＲＡＭに展開し、展開したプログラムに従って超音波診断装置１００の各部の動作を集中制御する。
【００３１】
ＲＯＭは、半導体等の不揮発メモリ等により構成され、超音波診断装置１００に対応するシステムプログラム、およびシステムプログラム上で実行可能な各種処理プログラムや、ガンマテーブル等の各種データ等を記憶する。これらのプログラムは、コンピュータが読み取り可能なプログラムコードの形態で格納され、ＣＰＵは、当該プログラムコードに従った動作を逐次実行する。ＲＡＭは、ＣＰＵにより実行される各種プログラムおよびこれらプログラムに係るデータを一時的に記憶するワークエリアを形成する。
【００３２】
次に、回帰推定値生成部１６が行う回帰分析処理について詳細に説明する。まず、本実施の形態の超音波診断装置１００における超音波探触子２の送信ビームの角度と被検体内の体内組織の反射面の方向との関係について説明する。
【００３３】
図３は、超音波探触子２から送信される送信ビームの角度と操作領域との関係を示す図である。例えば、図３に示すように、超音波診断装置１００は、超音波探触子２から出力する送信ビームの角度を、方位方向に対して垂直な方向にして走査を行うことで走査領域Ｑ１のコンポーネント画像データを取得する。なお、本実施の形態においては、超音波探触子２から出力される送信ビームの角度をステア角度（ステアリング角度）と称する。より具体的には、ステア角度は、超音波探触子２の中心軸と超音波ビーム送信方向とのなす角である。
【００３４】
そして、超音波診断装置１００は、ステアリング角度を左右にそれぞれ所定角度（例えば、１０°）つけてそれぞれ走査を行う場合、走査領域Ｑ２のコンポーネント画像データおよび走査領域Ｑ３のコンポーネント画像データを取得することができる。走査領域Ｑ１〜Ｑ３は、少なくとも一部が重複している。そして、走査領域Ｑ１〜Ｑ３の各コンポーネント画像データを合成し、合成画像データを得る。このようにして得られた合成画像データは、上述したように合成した画像の枚数に応じてノイズやスペックルが低減された超音波画像データである。
【００３５】
なお、図３に示す例では、ステアリング角度が０°の場合と±１０°の場合の３つのコンポーネント画像データを合成して合成画像を作成しているが、実際には超音波診断装置１００は、より多くのステアリング角度でより多くのコンポーネント画像データを取得し、これらを合成して超音波画像データを生成することが望ましい。
【００３６】
図４Ａおよび図４Ｂは、超音波探触子２に正対した体内組織に対する超音波の送信ビームと反射ビームとを例示した図である。超音波探触子２に正対した体内組織とは、超音波探触子２からの超音波ビーム送信方向と体内組織の反射面とのなす角が９０°である体内組織を意味する。なお、以下では、被検体内の走査対象の体内組織を、ターゲットと記載する。
【００３７】
図４Ａは、超音波探触子２に対して正対したターゲットに対してステア角度０°で超音波ビームを送信する様子を例示した図であり、図４Ｂは、超音波探触子２に対して正対したターゲットに対してステア角度θ（θ≠０）で超音波ビームを送信する様子を例示した図である。
【００３８】
図４Ａに示すように、超音波探触子２に正対したターゲットから反射される超音波は、ステア角度０°の場合、その多くが超音波探触子２に受信される。一方、ステア角度がθである場合、超音波探触子２が受信できる反射ビームの量は、図４Ｂに示すように、図４Ａに示すステア角度０°の場合と比較して少なくなる。換言すれば、ターゲットが超音波探触子２に正対しているとき、ステア角度がθである場合に超音波探触子２が受信できる反射信号の強度は、ステア角度０°の場合と比較して弱くなる。
【００３９】
一方、図５Ａおよび図５Ｂは、超音波探触子２に正対していないターゲットに対する超音波の送信ビームと反射ビームとを例示した図である。図５Ａは、超音波探触子２に対して正対したターゲットに対してステア角度０°で超音波ビームを送信する様子を例示した図であり、図５Ｂは、超音波探触子２に対して正対したターゲットに対してステア角度θ（θ≠０）で超音波ビームを送信する様子を例示した図である。
【００４０】
図５Ｂに示すように、送信ビーム方向に対してターゲットの反射面が傾いている角度と、ステア角度θとがほぼ一致する場合、超音波探触子２が受信できる反射ビームの量は、図４に例示した、正対したターゲットに対するステア角度０°で超音波ビームを送信した場合と同等となる。すなわち、送信ビーム方向に対してターゲットの反射面が傾いている角度と、ステア角度θとがほぼ一致する場合、超音波探触子２が受信できる反射ビームの量は、図５Ａに示す、超音波探触子２に正対していないターゲットに対してステア角度０°で超音波ビームを送信すると比較して多くなる。換言すれば、ターゲットが超音波探触子２に正対していないとき、送信ビーム方向に対してターゲットの反射面が傾いている角度とほぼ一致するステア角度で超音波ビームを送信した場合に超音波探触子２が受信できる反射信号の強度は、ステア角度０°の場合と比較して強くなる。
【００４１】
なお、以下では、超音波探触子２に正対したターゲットを正対ターゲット、超音波探触子２に正対していないターゲットを非正対ターゲットと記載する。
【００４２】
図６Ａから図６Ｃは、反射信号の強度とステア角度との関係を例示した図である。図６Ａは、正対ターゲットからの反射信号強度とステア角度との関係を、図６Ｂは、非正対ターゲットからの反射信号強度とステア角度との関係を、図６Ｃは、スペックルノイズを含む反射信号強度とステア角度との関係を、それぞれ示している。
【００４３】
正対ターゲットの場合、図４Ａおよび図４Ｂにて説明したように、ステア角度が０°である場合に反射信号強度は最も大きくなり、ステア角度の絶対値が大きくなるにつれて反射信号強度が小さくなる。図６Ａは、正対ターゲットからの反射信号強度とステア角度とを同一平面にプロットした図である。
【００４４】
一方、非正対ターゲットの場合、図５Ａおよび図５Ｂにて説明したように、送信ビーム方向に対してターゲットの反射面が傾いている角度が＋θであるとすると、ステア角度を＋θより小さい角度から＋θまで次第に大きくしていくと、ステア角度が＋θになるまで反射信号強度は大きくなる。図６Ｂは、非正対ターゲットからの反射信号強度とステア角度とを同一平面にプロットした図である。なお、図６Ｂの例では、角度θは１５°より大きい。
【００４５】
なお、図６Ｂは送信ビーム方向に対してターゲットの反射面が傾いている角度が＋θである場合について例示しているが、例えば送信ビーム方向に対してターゲットの反射面が傾いている角度が−θである場合は、上記と反対に、ステア角度を−θより大きい角度から−θまで次第に小さくしていった場合に、ステア角度が−θになるまで反射信号強度は徐々に大きくなる。
【００４６】
スペックルノイズは、超音波ビームが被検体内の境界面等において散乱され、散乱ビームが干渉して生じたノイズである。このため、スペックルノイズを含む反射信号の信号強度は、ステア角度に依存せず、図６Ｃに示すようになる。図６Ｃは、スペックルノイズを含む反射信号の信号強度とステア角度とを同一平面にプロットした図である。
【００４７】
回帰推定値生成部１６は、このようなステア角度と反射信号強度との関係に基づいて、回帰分析を行う。具体的には、回帰推定値生成部１６は、画像生成部１４が生成した、複数のステア角度に対応する複数の超音波データにおいて、画素毎、あるいは所定の領域毎に、反射信号強度に対応する画素値（輝度値）に基づく回帰分析を行う。所定の領域とは、例えば、互いに隣接する所定数の画素数で構成される領域や、画素値（輝度値）の大きさに基づいて分類された互いに隣接する領域である。所定の領域は、例えば１つの画素で構成される領域で会ってもよい。
【００４８】
以下では、回帰推定値生成部１６が、単回帰分析を使用する場合の回帰推定値算出方法について説明する。
【００４９】
例えば上記図６Ａから図６Ｃに示した、正対ターゲット、非正対ターゲット、およびスペックルノイズに対応する反射信号の強度とステア角度との関係に基づいて、単回帰分析により得られる１次回帰式を図７Ａから図７Ｃに例示する。このような回帰分析により得られるパラメータとして、例えば以下の表１に示すようなパラメータがある。表１は、単回帰分析により得られる各種パラメータの大きさの関係を示した表である。
【００５０】
【表１】
2018051040.tif 000003
【００５１】
表１において、ａ１は回帰分析により算出された１次回帰式の傾き、ａ０は回帰分析により算出された１次回帰式の切片、ＲＳＳは残差平方和（residual sum of squares）であり、Ｒ２は決定係数である。すなわち、１次回帰式はｙ＝ａ１ｘ＋ａ０で表される。ＲＳＳは、実測値と推定モデルの各点の差の二乗和であり、実測値と推定モデルとの間の不一致を評価する尺度である。また、Ｒ２は回帰分析の精度に関する尺度である。
【００５２】
図７Ａおよび表１に示すように、正対ターゲットでは、１次回帰式の傾き｜ａ１｜の大きさは非正対ターゲットと比較して小さくなり、切片ａ０＊の大きさは大きくなる。図７Ｂに示すように、非正対ターゲットでは、１次回帰式の傾き｜ａ１｜の大きさは他のターゲットあるいはスペックルノイズの場合と比較して大きくなる。図７Ｃに示すように、スペックルノイズは信号強度に対してランダムな値をとるので、１次回帰式の傾き｜ａ１｜の大きさは非正対ターゲットと比較して小さくなる。
【００５３】
回帰推定値生成部１６は、このような回帰分析により得られたパラメータ（１次回帰式の傾き｜ａ１｜、切片ａ０、ＲＳＳ、Ｒ２等）に基づいて、以下のような算出方法により回帰推定値αを算出し、画像合成部１７に出力する。
【００５４】
なお、本発明において、回帰推定値αは、非正対ターゲットに対応する領域の画素値に対しては重みが相対的に重く、それ以外の領域の画素値に対しては重みが軽くなるように、重み付けを行う重み付け値である。従って、回帰推定値生成部１６は、非正対ターゲットに対応する領域の画素値に対してはαが大きくなるように、それ以外の領域の画素値に対してはαが小さくなるように、αを決定する。
【００５５】
回帰推定値生成部１６が行う回帰推定値αの導出は、例えば以下例示する数式のうちのいずれかを用いて行われればよい。
【００５６】
（回帰推定値αの導出方法１）
α＝｜ａ１｜        （１）
【００５７】
（回帰推定値αの導出方法２）
α＝｜ａ１｜・γ        （２）
但し、γは所定の係数である。
【００５８】
（回帰推定値αの導出方法３）
【数１】
2018051040.tif 000004
【００５９】
（回帰推定値αの導出方法４）
α＝｜ａ１｜・（ＲＳＳｍａｘ−ＲＳＳ）        （４）
【００６０】
（回帰推定値αの導出方法５）
α＝｜ａ１｜・Ｒ２        （５）
【００６１】
なお、本実施の形態では回帰推定値生成部１６による回帰推定値αの導出方法として、式（１）から式（５）の５通りの導出方法を例示したが、本発明はこれに限定されず、他の導出方法を用いてもよい。上述したように、回帰推定値生成部１６は、非正対ターゲットに対応する領域に対してはαが大きくなるように、それ以外の領域に対してはαが小さくなるように、回帰推定値αを決定すればよい。また、回帰推定値生成部１６は、回帰推定値αが０以上１以下の範囲となるように、正規化処理やクリップ処理等を行ってもよい。
【００６２】
次に、回帰推定値生成部１６が算出した回帰推定値αを用いて、画像合成部１７が行う画像合成処理について説明する。
【００６３】
画像合成部１７は、互いにステア角度が異なる超音波ビームの反射信号に基づいて画像生成部１４が生成した複数の超音波画像データの画素毎、あるいは所定の領域毎の画素値（輝度値）を以下のように合成することで、複数の超音波画像データを合成する。
【００６４】
（画像合成方法１）
１つ目の画像合成方法１において、画像合成部１７は、回帰推定値αの絶対値｜α｜を用いて、複数の超音波画像データの所定の領域毎の画素値に対してαブレンドを行う。例えば５枚の超音波画像データＡ〜Ｅを用いる場合の合成式は以下の式（６）のようになる。
【数２】
2018051040.tif 000005
【００６５】
ここで、式（６）において、ｆ１およびｆ２は所定の関数である。関数ｆ１およびｆ２としては、例えばパラメータの単純平均を取るＭＥＡＮ関数、最大値を抽出するＭＡＸ関数、最小値を抽出するＭＩＮ関数、中央値を抽出するＭＥＤＩＡＮ関数、分散値を取るＶＡＲＩＡＮＣＥ関数、標準偏差を取るＳＤ関数、加重平均を取るＷｅｉｇｈｔｅｄ−ＭＥＡＮ関数等のうちのいずれかが使用されればよい。
【００６６】
画像合成部１７がどの関数を使用するかは、例えば回帰推定値生成部１６が生成した回帰推定値αの導出方法や、最終的に超音波診断装置１００が生成する超音波画像の使用目的等に合わせて適宜選択されればよい。なお、この選択は、例えば操作入力部１１を用いた超音波診断装置１００の使用者（医師等）により行われてもよいし、例えば超音波診断装置１００の制御部１９によって自動的に行われるようにしてもよい。関数ｆ１とｆ２とは同じでもよいし、異なる関数であってもよい。
【００６７】
（画像合成方法２）
２つ目の画像合成方法２において、画像合成部１７は、回帰推定値生成部１６が生成した回帰推定値αと、所定の係数βと、を用いて、以下の式（７）に基づいて複数の超音波画像データの所定の領域毎に画素値の合成を行う。
【数３】
2018051040.tif 000006
【００６８】
式（７）において、所定の係数βは、例えば操作入力部１１を用いた超音波診断装置１００の使用者（医師等）により任意に設定されてもよいし、制御部１９が自動的に決定してもよいし、予め設定された固定値であってもよい。関数ｆ１、ｆ２、および所定の係数βを制御部１９が自動的に決定する場合、超音波診断装置１００の使用者は、関数および係数を決定する手間と操作の手間なしで好適な合成画像データを得ることができる。関数ｆ１、ｆ２、および所定の係数βを使用者が任意で決定する場合、目的に合わせて使用者が画像の表示具合を好適に調整することができる。
【００６９】
以上説明したような画像合成方法によって、画像合成部１７は、非正対ターゲットに対応する領域に対しては重みが相対的に重くなるように、それ以外の領域に対しては重みが軽くなるように、重み付けをして複数の超音波画像データの合成を行うことができる。
【００７０】
以上説明したように、本発明の超音波診断装置１００は、複数の異なるステアリング角度で超音波を送信し、被検体の反射面で反射された反射超音波を受信する超音波探触子２が受信した反射超音波に基づいて生成した受信信号に基づいて、複数の異なるステアリング角度に対応する受信信号からそれぞれ異なる複数の超音波画像データを生成する画像生成部１４と、複数の異なるステアリング角度と、複数の異なるステアリング角度に対応する複数の超音波画像データと、に基づいて回帰分析を行い、回帰分析の結果に基づいて所定の重み付け値である回帰推定値を生成する回帰推定値生成部１６と、複数の超音波画像データに対して、回帰推定値に基づく重み付けを行って合成を行い、合成画像データを生成する画像合成部１７と、を備える。
【００７１】
回帰推定値生成部１６は、超音波画像データの画素毎あるいは所定領域毎に、画素値（輝度値）に基づく回帰分析により、送信ビームに対して反射面が傾いたターゲット（非正対ターゲット）に対して比重が重くなるような回帰推定値αを算出する。このため、画像合成部１７は、非正対ターゲットに対応する領域に対しては比重が重くなるように、それ以外の領域に対しては比重が相対的に軽くなるように、重み付けを行って合成画像データを生成することができる。このため、本発明の超音波診断装置１００は、このような重み付けを行わない場合と比較して、非正対ターゲットが比較的明瞭に描写された合成画像データを生成することができる。
【００７２】
以上、図面を参照しながら本発明の実施の形態について説明したが、本発明はかかる例に限定されない。特許請求の範囲の記載範囲内において、当業者が想到できる各種の変更例または修正例についても、本発明の技術的範囲に含まれる。また、開示の趣旨を逸脱しない範囲において、上記実施の形態における各構成要素を任意に組み合わせてもよい。
【００７３】
上記実施の形態では、超音波探触子２の振動子２ａは、方位方向（走査方向）に一次元アレイ状に複数配列されているとしたが、二次元アレイ状に配列されたものであってもよい。また、上記実施の形態では、超音波探触子２としてリニア電子スキャンプローブを用いて、リニア走査方式による超音波の走査を行うものとしたが、セクタ走査方式あるいはコンベックス走査方式等、他の方式を採用してもよい。
【００７４】
上記実施の形態では、回帰推定値生成部１６は、単回帰分析により回帰推定値を生成していたが、本発明はこれに限定されない。回帰推定値生成部１６は、例えば、２次回帰式やそれ以外の回帰分析を行ってもよい。
【００７５】
以下、回帰推定値生成部１６が２次回帰式を生成する場合について説明する。図８Ａから図８Ｃは、上記図６Ａから図６Ｃに示した、正対ターゲット、非正対ターゲット、およびスペックルノイズに対応する反射信号の強度とステア角度との関係に基づいて、回帰分析により得られる２次回帰式ｙ＝ａ２ｘ２＋ａ１ｘ＋ａ０を例示した図である。
【００７６】
図８Ａおよび図８Ｃに示すように、正対ターゲットおよびスペックルノイズに対応する反射信号の強度とステア角度との関係式は、２次曲線に近い関係を示す。しかし、図８Ｂに示すように、非正対ターゲットに対応する反射信号の強度とステア角度との関係式は、１次回帰式の場合と同様に、ほぼ直線状の１次直線として表れる。このため、２次線形回帰分析により得られるパラメータは、以下の表２のような大小関係となる。
【００７７】
【表２】
2018051040.tif 000007
【００７８】
このように、２次線形回帰分析により得られるパラメータの場合でも、非正対ターゲットに対応する独立変数の係数（パラメータ｜ａ２｜および｜ａ１｜）は、正対ターゲットおよびスペックルノイズに対応する係数と比較して大きくなる。このため、回帰推定値生成部１６は、パラメータ｜ａ２｜および｜ａ１｜に基づいて、非正対ターゲットに対応する領域に対してはαが大きくなるように、それ以外の領域に対してはαが小さくなるように、回帰推定値αを決定することができる。
【産業上の利用可能性】
【００７９】
本発明は、複数のステア角度の超音波送信ビームに対する反射信号に基づいて画像を合成する超音波診断装置に好適である。
【符号の説明】
【００８０】
１００  超音波診断装置
１  超音波診断装置本体
２  超音波探触子
２ａ  振動子
１１  操作入力部
１２  送信部
１３  受信部
１４  画像生成部
１５  ＤＳＣ
１６  回帰推定値生成部
１７  画像合成部
１８  表示部
１９  制御部
【図１】
2018051040.tif 000008
【図２】
2018051040.tif 000009
【図３】
2018051040.tif 000010
【図４Ａ】
2018051040.tif 000011
【図４Ｂ】
2018051040.tif 000012
【図５Ａ】
2018051040.tif 000013
【図５Ｂ】
2018051040.tif 000014
【図６Ａ】
2018051040.tif 000015
【図６Ｂ】
2018051040.tif 000016
【図６Ｃ】
2018051040.tif 000017
【図７Ａ】
2018051040.tif 000018
【図７Ｂ】
2018051040.tif 000019
【図７Ｃ】
2018051040.tif 000020
【図８Ａ】
2018051040.tif 000021
【図８Ｂ】
2018051040.tif 000022
【図８Ｃ】
2018051040.tif 000023