(57)【要約】【課題】  超音波イメージングにおいて、事前に予測することが困難なダイナミックレンジを最適な設定に自動調整することが可能な超音波診断装等を提供すること。【解決手段】  複数フレームに対応する画像データの各輝度を用いて統計的処理を行い、被検体からのエコー信号に対応する輝度に基づいてダイナミックレンジの上限を自動的に決定し、また、ノイズレベルに対応する輝度に基づいて同のダイナミックレンジの下限値を自動的に決定するダイナミックレンジの自動調整機能を有する。本機能により、マニュアル操作等により最初に設定したダイナミックレンジが広すぎたり狭すぎたりする場合、画像上の輝度分布が動的に変動する場合等であっても、ダイナミックレンジを最適な設定に自動調整することができる。【選択図】      図３
2008264530.tif 000002
【特許請求の範囲】
【請求項１】
複数フレームに対応する複数の第１の画像データを取得する取得ユニットと、
前記複数の第１の画像データのうち、被検体からのエコー信号に対応する輝度に基づいて上限値を決定する決定ユニットと、
前記複数の画像データと前記上限値を少なくとも用いて定義されるダイナミックレンジとを用いて、前記複数フレームのうちの少なくとも一つに対応する第２の画像データを生成する画像生成ユニットと、
を具備する超音波診断装置。
【請求項２】
前記決定ユニットは、前記複数の第１の画像データから取得されるノイズレベルに対応する輝度に基づいて下限値を決定し、
前記ダイナミックレンジは、前記上限値と前記下限値とを用いて定義されるものである請求項１記載の超音波診断装置。
【請求項３】
前記決定ユニットは、
前記被検体に注入された造影剤に起因するエコー信号に対応する輝度に基づいて前記上限値を決定する請求項１又は２記載の超音波診断装置。
【請求項４】
前記決定ユニットは、
前記複数の第１の画像データを用いて階調に関するヒストグラムを生成し、
前記ヒストグラムを用いて、前記上限値及び前記下限値のうちの少なくとも一方を決定する請求項２又は３記載の超音波診断装置。
【請求項５】
前記複数フレームのうちの少なくとも一つは、超音波受信のみが実行されたフレームであり、
前記決定ユニットは、前記超音波受信のみが実行されたフレームに対応する前記第１の画像データに含まれるホワイトノイズに対応する輝度に基づいて前記下限値を決定する請求項２乃至４のうちいずれか一項記載の超音波診断装置。
【請求項６】
前記複数の第１の画像データの収集期間に関する開始時刻及び終了時刻のうちの少なくとも一方を設定する設定ユニットをさらに具備する請求項１乃至５のうちいずれか一項記載の超音波診断装置。
【請求項７】
前記設定ユニットは、前記第１の画像データの表示をフリーズさせるフリーズ指示の操作タイミングに基づいて、前記終了時刻を設定する請求項６記載の超音波診断装置。
【請求項８】
前記設定ユニットは、被検体への造影剤注入タイミングを基準として、前記開始時刻を設定する請求項６記載の超音波診断装置。
【請求項９】
前記設定ユニットは、フラッシュ送信のタイミングを基準として、前記開始時刻を設定する請求項６記載の超音波診断装置。
【請求項１０】
前記画像生成ユニットは、前記第１の画像データの表示をフリーズさせるフリーズ指示に応答して、前記複数フレームのうちの少なくとも一つに対応する前記第２の画像データを生成する請求項１乃至９のうちいずれか一項記載の超音波診断装置。
【請求項１１】
前記画像生成ユニットは、前記ダイナミックレンジを用いて複数の前記第２の画像データを生成し、
前記複数の第２の画像データをループ再生する再生ユニットをさらに具備する請求項１乃至１０のうちいずれか一項記載の超音波診断装置。
【請求項１２】
前記決定ユニットは、前記第２の画像データの最大の輝度に基づいて前記上限値を決定する請求項１乃至１１のうちいずれか一項記載の超音波診断装置。
【請求項１３】
前記ダイナミックレンジに含まれる前記複数の画像データの輝度値の平均値が、所定の値と同じ又は前記所定の値を規準とする範囲に含まれるように、前記決定された前記上限値を調整する調整ユニットをさらに具備する請求項１乃至１２のうちいずれか一項記載の超音波診断装置。
【請求項１４】
複数フレームに対応する複数の第１の画像データを記憶する記憶ユニットと、
前記複数の第１の画像データのうち、被検体からのエコー信号に対応する輝度に基づいて上限値を決定する決定ユニットと、
前記複数の画像データと前記上限値を少なくとも用いて定義されるダイナミックレンジとを用いて、前記複数フレームのうちの少なくとも一つに対応する第２の画像データを生成する画像生成ユニットと、
を具備する超音波画像処理装置。
【請求項１５】
前記決定ユニットは、前記複数の第１の画像データから取得されるノイズレベルに対応する輝度に基づいて下限値を決定し、
前記ダイナミックレンジは、前記上限値と前記下限値とを用いて定義されるものである請求項１４記載の超音波画像処理装置。
【請求項１６】
前記決定ユニットは、
前記被検体に注入された造影剤に起因するエコー信号に対応する輝度に基づいて前記上限値を決定する請求項１４又は１５記載の超音波画像処理装置。
【請求項１７】
前記決定ユニットは、
前記複数の第１の画像データを用いて階調に関するヒストグラムを生成し、
前記ヒストグラムを用いて、前記上限値及び前記下限値を決定する請求項１５記載の超音波画像処理装置。
【請求項１８】
前記複数フレームのうちの少なくとも一つは、超音波受信のみが実行されたフレームであり、
前記決定ユニットは、前記超音波受信のみが実行されたフレームに対応する前記第１の画像データに含まれるホワイトノイズに対応する輝度に基づいて前記下限値を決定する請求項１５記載の超音波画像処理装置。
【請求項１９】
前記複数の第１の画像データは、前記第１の画像データの表示をフリーズさせるフリーズ指示の操作タイミングを基準として収集されたものである請求項１４乃至１８のうちいずれか一項記載の超音波画像処理装置。
【請求項２０】
前記複数の第１の画像データは、被検体への造影剤注入タイミングを基準として収集されたものである請求項１４乃至１８のうちいずれか一項記載の超音波画像処理装置。
【請求項２１】
前記複数の第１の画像データは、フラッシュ送信のタイミングを基準として収集されたものである請求項１４乃至１８のうちいずれか一項記載の超音波画像処理装置。
【請求項２２】
前記画像生成ユニットは、前記第１の画像データの表示をフリーズさせるフリーズ指示に応答して、前記複数フレームのうちの少なくとも一つに対応する前記第２の画像データを生成する請求項１４乃至１８のうちいずれか一項記載の超音波画像処理装置。
【請求項２３】
前記画像生成ユニットは、前記ダイナミックレンジを用いて複数の前記第２の画像データを生成し、
前記複数の第２の画像データをループ再生する再生ユニットをさらに具備する請求項１４乃至２２うちいずれか一項記載の超音波画像処理装置。
【請求項２４】
前記決定ユニットは、前記第２の画像データの最大の輝度に基づいて前記上限値を決定する請求項１４乃至２３のうちいずれか一項記載の超音波画像処理装置。
【請求項２５】
コンピュータに、
複数フレームに対応する複数の第１の画像データのうち、被検体からのエコー信号に対応する輝度に基づいて上限値を決定させる決定機能と、
前記複数の画像データと前記上限値を少なくとも用いて定義されるダイナミックレンジとを用いて、前記複数フレームのうちの少なくとも一つに対応する第２の画像データを生成させる生成機能と、
を実現させる超音波画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波造影剤を用いて行う造影エコー法において、事前に予測することが困難な、GAINおよびダイナミックレンジの設定を、視認性が良くなるように自動的に調整することが可能な超音波診断装置、超音波画像処理装置、超音波画像処理方法等に関する。
【背景技術】
【０００２】
超音波診断装置は生体内情報の画像を表示する診断装置であり、Ｘ線診断装置やＸ線コンピュータ断層撮影装置などの他の画像診断装置に比べ、安価で被爆が無く、非侵襲性に実時間で観測するための有用な装置として利用されている。係る特性から、超音波診断装置の適用範囲は広く、心臓などの循環器から肝臓、腎臓などの腹部、抹消血管、産婦人科、脳血管などの診断に利用されている。
【０００３】
近年、静脈投与型の超音波造影剤が製品化され、造影エコー法が行われるようになってきている。この手法は、例えば、心臓および肝臓などの検査で静脈から超音波造影剤を注入して血流信号を増強し、血流動態の評価を行うことを目的としている。造影剤の多くは微小気泡（マイクロバブル）が反射源として機能するものである。超音波診断装置メーカーでは、造影剤からの反射信号を感度良く映像化するために、送受信方法に工夫を凝らし、バブルからの非線形な信号成分を感度良く受信することが可能となってきている。
【０００４】
ところで、超音波診断装置では、画像診断を行う場合に調節しなければならないパラメータが非常に多いため、その操作は煩雑であると言える。この問題に関して、各メーカーは、ユーザサポート機能を開発しており、種々のゲイン、ＳＴＣ（Sensitivity Time Control）の自動調節機能等を装置に搭載することで対処している。
【０００５】
なお、本願に関連する公知文献としては、例えば次のようなものがある。
【特許文献１】特表２００４−５００９１５号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
しかしながら、従来の超音波診断装置を用いて造影エコー法による撮影を行う場合、依然として例えば次の様な問題がある。
【０００７】
すなわち、造影剤を投与したときの、信号輝度レベルは、実際に投与してみないと分からないことが多い。患者や周波数等のパラメータやプローブの当て方により、同じ肝臓を見る場合でも輝度が異なる。つまり、造影剤投与前に適当に設定したダイナミックレンジでは、広すぎてしまったり、逆に狭すぎて輝度が飽和してしまったりする危険性がある。
【０００８】
また、造影剤エコー法では、造影剤の流入に従って（時間経過に従って）画像が暗→明→暗と変化するという特徴がある。例えば、肝臓の血流状態を造影剤で観察するとき、最初に動脈が勢い良く染まり、その後、少し時間をおいて、門脈が染まるといったように、輝度が単調増加といったような簡単な変化でなく、複雑な時間変化を示す。従って、造影剤投与後でも、あるタイミングに調節したダイナミックレンジが、他のタイミングでは不適切な設定となってしまう場合がある。また、ユーザが血管を流れる血流を見たいのか、パフュージョンを見たいのかという目的によっても、ダイナミックレンジの調整法は異なってくる。
【０００９】
本発明は、上記事情を鑑みてなされたもので、超音波イメージングにおいて、事前に予測することが困難なダイナミックレンジを最適な設定に自動調整することが可能な超音波診断装置又は超音波画像処理装置及び超音波画像処理プログラムを提供することを目的としている。
【課題を解決するための手段】
【００１０】
本発明は、上記目的を達成するため、次のような手段を講じている。
【００１１】
請求項１に記載の発明は、複数フレームに対応する複数の第１の画像データを取得する取得ユニットと、前記複数の第１の画像データのうち、被検体からのエコー信号に対応する輝度に基づいて上限値を決定する決定ユニットと、前記複数の画像データと前記上限値を少なくとも用いて定義されるダイナミックレンジとを用いて、前記複数フレームのうちの少なくとも一つに対応する第２の画像データを生成する画像生成ユニットと、を具備する超音波診断装置である。
【００１２】
請求項１４に記載の発明は、複数フレームに対応する複数の第１の画像データを記憶する記憶ユニットと、前記複数の第１の画像データのうち、被検体からのエコー信号に対応する輝度に基づいて上限値を決定する決定ユニットと、前記複数の画像データと前記上限値を少なくとも用いて定義されるダイナミックレンジとを用いて、前記複数フレームのうちの少なくとも一つに対応する第２の画像データを生成する画像生成ユニットと、を具備する超音波画像処理装置である。
【００１３】
請求項２５に記載の発明は、コンピュータに、複数フレームに対応する複数の第１の画像データのうち、被検体からのエコー信号に対応する輝度に基づいて上限値を決定させる決定機能と、前記複数の画像データと前記上限値を少なくとも用いて定義されるダイナミックレンジとを用いて、前記複数フレームのうちの少なくとも一つに対応する第２の画像データを生成させる生成機能と、を実現させる超音波画像処理プログラムである。
【発明の効果】
【００１４】
以上本発明によれば、超音波イメージングにおいて、事前に予測することが困難なダイナミックレンジを最適な設定に自動調整することが可能な超音波診断装置又は超音波画像処理装置及び超音波画像処理プログラムを実現することができる。
【発明を実施するための最良の形態】
【００１５】
以下、本発明の実施形態を図面に従って説明する。なお、以下の説明において、略同一の機能及び構成を有する構成要素については、同一符号を付し、重複説明は必要な場合にのみ行う。
【００１６】
図１は、本実施形態に係る超音波診断装置１０のブロック構成図を示している。同図に示すように、本超音波診断装置１０は、超音波プローブ１２、入力装置１３、モニター１４、送受信ユニット２１、Ｂモード処理ユニット２２、ドプラ処理ユニット２３、画像生成ユニット２４、制御プロセッサ（ＣＰＵ）２５、記憶ユニット２６、画像メモリ２７、ソフトウェア格納ユニット２８、インターフェースユニット２９を具備している。
【００１７】
超音波プローブ１２は、圧電セラミック等の音響／電気可逆的変換素子としての圧電振動子を有する。複数の圧電振動子は並列され、プローブ１２の先端に装備される。
【００１８】
なお、本超音波プローブ１２は、被検体の３次元領域を超音波走査可能なものであってもよい。係る場合、超音波プローブ１２は、振動子をその配列方向の直交方向に沿って機械的に揺動させ、３次元領域を超音波走査する構成、又は二次元的に配列された二次元振動素子を用いて電気的制御により３次元領域を超音波走査する構成等を有する。前者の構成を採用する場合、被検体の３次元的走査は前記揺動回路によって行われるため、検査者はプローブ本体を被検体に接触させるだけで、自動的に複数の二次元断層像を取得することができる。制御された揺動速度から断面間の正確な距離も検知できる。また、後者の構成を採用する場合には、原理的には、従来の二次元断層像を取得するのと同じ時間で、３次元領域を超音波走査することができる。
【００１９】
入力装置１３は、各種指示・命令・情報を入力するためのスイッチ、ボタン、マウス、キーボード、トラックボールを有している。
【００２０】
なお、入力装置１３は、図２に示すように、後述するダイナミックレンジ自動調整機能の動作モードをＯＮ／ＯＦＦ指示するためのＯＮ／ＯＦＦボタン１３０、当該ダイナミックレンジ自動調整機能を適用する対象期間を指定するためのスタート／ストップボタン１３１を有している。
【００２１】
モニター１４は、画像生成ユニット２４からのビデオ信号に基づいて、生体内の形態学的情報（Ｂモード画像）、血流情報（平均速度画像、分散画像、パワー画像等）、後述するダイナミックレンジ調整機能により最適化されたダイナミックレンジに従う超音波画像等を所定の形態で表示する。
【００２２】
送受信ユニット２１は、図示しないトリガ発生回路、遅延回路およびパルサ回路等を有している。パルサ回路では、所定のレート周波数ｆｒ  Ｈｚ（周期；１／ｆｒ秒）で、送信超音波を形成するためのレートパルスが繰り返し発生される。また、遅延回路では、チャンネル毎に超音波をビーム状に集束し且つ送信指向性を決定するのに必要な遅延時間が、各レートパルスに与えられる。この遅延情報を変化させることで、プローブ振動子面からの送信方向を任意に調整することが可能となる。トリガ発生回路は、このレートパルスに基づくタイミングで、プローブ１２に駆動パルスを印加する。
【００２３】
なお、送受信ユニット２１は、制御プロセッサ２５の指示に従って、遅延情報、送信周波数、送信駆動電圧等を瞬時に変更可能な機能を有している。特に送信駆動電圧の変更については、瞬間にその値を切り替え可能なリニアアンプ型の発信回路、又は複数の電源ユニットを電気的に切り替える機構によって実現される。
【００２４】
また、送受信ユニット２１は、図示していないアンプ回路、Ａ／Ｄ変換器、加算器等を有している。アンプ回路では、プローブ１２を介して取り込まれたエコー信号をチャンネル毎に増幅する。Ａ／Ｄ変換器では、増幅されたエコー信号に対し受信指向性を決定するのに必要な遅延時間を与え、その後加算器において加算処理を行う。この加算により、エコー信号の受信指向性に応じた方向からの反射成分が強調され、受信指向性と送信指向性とにより超音波送受信の総合的なビームが形成される。
【００２５】
Ｂモード処理ユニット２２は、送受信ユニット２１からエコー信号を受け取り、対数増幅、包絡線検波処理などを施し、信号強度が輝度の明るさで表現されるデータを生成する。
【００２６】
ドプラ処理ユニット２３は、送受信ユニット２１から受け取ったエコー信号から速度情報を周波数解析し、ドプラ効果による血流や組織、造影剤エコー成分を抽出し、平均速度、分散、パワー等の血流情報を多点について求める。
【００２７】
画像生成ユニット２４は、超音波スキャンの走査線信号列を、テレビなどに代表される一般的なビデオフォーマットの走査線信号列に変換し、表示画像としての超音波診断画像を生成する。なお、当該画像生成ユニット２３に入る以前のデータは、「生データ」と呼ばれることがある。
【００２８】
制御プロセッサ２５は、情報処理装置（計算機）としての機能を持ち、本超音波診断装置全体の動作を制御する。制御プロセッサ２５は、記憶ユニット２６からダイナミックレンジの自動調整機能を実現するための専用プログラム、所定のスキャンシーケンス、画像生成・表示等を実行するための制御プログラムを読み出してソフトウェア格納ユニット２８上に展開し、各種処理に関する演算・制御等を実行する。
【００２９】
記憶ユニット２６は、磁気ディスク（フロッピー（登録商標）ディスク、ハードディスクなど）、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤなど）、半導体メモリなどの記録媒体、及びこれらの媒体に記録された情報を読み出す装置である。この記憶ユニット２６は、種々のスキャンシーケンス、後述するダイナミックレンジの自動調整機能を実現するための専用プログラム、画像生成、表示処理を実行するための制御プログラムや、診断情報（患者ＩＤ、医師の所見等）、診断プロトコル、送受信条件、その他のデータ群を格納している。また、記憶ユニット２６は、必要に応じて、画像メモリ２７中の画像の保管などにも使用される。記憶ユニット２６内のデータは、インターフェースユニット２９を経由して外部周辺装置へ転送することも可能となっている。
【００３０】
画像メモリ２７は、画像生成ユニット２４において生成された画像データを記憶する。この画像メモリ２７に記憶される画像データは、モニター１４に表示される画像データに比して２倍程度のダイナミックレンジを有する。当該画像メモリ２７に記憶された画像データは、例えば診断の後に操作者が呼び出すことが可能となっており、静止画的に、あるいは複数枚を使って動画的に再生することが可能でなる。また、画像メモリ２７は、超音波送受信ユニット２１直後の出力信号（ｒａｄｉｏ  ｆｒｅｑｕｅｎｃｙ(ＲＦ)信号と呼ばれる）、送受信ユニット２１通過後の画像輝度信号、その他の生データ、ネットワークを介して取得した画像データ等を必要に応じて記憶する。
【００３１】
インターフェースユニット２９は、入力装置１３、ネットワーク、新たな外部記憶装置（図示せず）に関するインターフェースである。当該装置によって得られた超音波画像等のデータや解析結果等は、インターフェースユニット２９よって、ネットワークを介して他の装置に転送可能である。
【００３２】
（ダイナミックレンジの自動調整機能）
次に、本超音波診断装置１０が有するダイナミックレンジの自動調整機能について説明する。この機能は、例えば造影剤エコー法の様に画像の輝度が経時的に変動する場合等において、ダイナミックレンジを自動調整するものである。
【００３３】
なお、本実施形態では、本ダイナミックレンジの自動調整機能が超音波診断装置１０によって実現される場合例として説明する。しかしながら、これに拘泥されず、例えば当該ダイナミックレンジの自動調整機能を実現するプログラムがインストールされたワークステーション等の超音波画像処理装置によって実現するようにしてもよい。
【００３４】
図３は、本ダイナミックレンジの自動調整機能に従う処理（ダイナミックレンジの自動調整処理）の流れを示したフローチャートである。
【００３５】
同図において、まず、制御プロセッサ２５は、造影剤が注入された被検体を超音波によって走査し、第１のダイナミックレンジを用いて複数フレームに対応する複数の第１の画像データを取得する（ステップＳａ）。
【００３６】
ここで、第１のダイナミックレンジとは、初期設定、マニュアル設定等によって、本ダイナミックレンジの自動調整処理の前段において設定されるものである。また、本ステップにおけるスキャンシーケンスには、造影剤エコー法に準ずるものであれば、特に限定はない。
【００３７】
次に、入力装置１３を介して本自動調整処理の開始指示を受けると、制御プロセッサ２５は、ステップＳ１において取得された複数フレームに対応する第１の画像データを用いて、輝度に関するヒストグラムを生成する（ステップＳｂ）。
【００３８】
図４は、ステップＳ２において生成されるヒストグラムの一例を示している。なお、同図においては、第１のダイナミックレンジの上限がＶＵ１で、下限がＶＬ１でそれぞれ示されている。
【００３９】
次に、画像生成ユニット２４は、生成された輝度に関するヒストグラムを用いて、第２のダイナミックレンジの上限ＶＵ２及び下限ＶＬ２を決定する（ステップＳｃ）。すなわち、画像生成ユニット２４は、ステップＳ２において生成されたヒストグラムにおいて、所定の閾値ＶＴＨを超える輝度のうちの最大値を、第２のダイナミックレンジの上限ＶＵ２として決定する。また、画像生成ユニット２４は、ヒストグラムにおいて、ホワイトノイズレベルに対応する輝度を、第２のダイナミックレンジの下限ＶＬ２として決定する。なお、図５に、ヒストグラムに設定された第２のダイナミックレンジの上限ＶＵ２、下限ＶＬ２をそれぞれ示した。
【００４０】
なお、第２のダイナミックレンジの上限ＶＵ２の決定手法は、上記例に拘泥されない。例えば、下限ＶＬ２から所定の最大値までの範囲の所定の割合（例えば１０％）だけ、所定の最大値から下方に低い値を、第２のダイナミックレンジの上限ＶＵ２として決定するようにしてもよい。また、上記ステップＳｂにおいて、輝度値に関するヒストグラムを各画像毎に作成し、ステップＳｃにおいて、各フレーム毎の最大輝度値を計算し、これを用いて全体の最大値を計算するようにしてもよい。さらに、対象期間Ｔに含まれる複数の画像のうち、最大輝度値を有する画像を選択し、当該画像を用いて第２のダイナミックレンジの上限ＶＵ２を決定するようにしてもよい。
【００４１】
第２のダイナミックレンジの下限ＶＬ２の決定をより好適にするために、ステップＳ１において取得する画像データに対応する複数フレームのうちの少なくとも一つは、空打ち後の超音波受信が実行された（超音波送信を行わず超音波受信のみが実行された）フレームとすることが好ましい。特に、造影剤投与直後に対応するフレーム、造影剤投与後の所定期間経過後（例えば、造影剤投与開始ボタン操作後１０秒後）に対応するフレーム、高音圧送信直後のフレーム等とすることがさらに好ましい。これにより、ステップＳ２において生成されるヒストグラムに、ホワイトノイズレベルに対応する輝度を正確に含めることができる。
【００４２】
次に、画像生成ユニット２４は、ステップＳ１において生成された第１の画像データのうち、第２のダイナミックレンジに属する輝度を有する信号を用いて、同じ複数フレームに対応する複数の第２の画像データを取得する（ステップＳｄ）。
【００４３】
以上述べたように、本ダイナミックレンジの自動調整処理により、通常の手法によって設定された第１のダイナミックレンジ（図４参照）は、造影剤エコー法についてより好適な第２のダイナミックレンジ（図５参照）に自動的に調整される。従って、通常の手法によって設定された第１のダイナミックレンジが広すぎたため、表示画像が例えば図６に示すように暗くなってしまった場合であっても、本自動調整機能によりダイナミックレンジを最適化することで、例えば図７に示すように観察に好適な階調を設定することができる。
【００４４】
なお、上記説明では、広すぎた第１のダイナミックレンジを、既述の自動調整機能によって、より狭い第２のダイナミックレンジに変更する場合を例に説明した。当該例に拘泥されず、例えば狭すぎた第１のダイナミックレンジを、本自動調整機能によって、より広い第２のダイナミックレンジに変更する場合についても、本自動調整機能は適用可能である。この場合においても、第１のダイナミックレンジが狭すぎたため、例えば図８に示すように信号が飽和気味となった画像を、図７に示すように観察に好適な階調を設定することができる。
【００４５】
本ダイナミックレンジの自動調整機能は、動画像（ＡＶＩファイルやＲａｗデータ）を対象とすることも可能である。
【００４６】
図９は、各フレームにおける最大輝度（階調）Ｘｍａｘと最小輝度（階調）Ｘｍｉｎ（すなわち、最大輝度Ｘｍａｘと最小輝度Ｘｍｉｎの時間変化）を示した図である。対象期間Ｔに属する複数のフレームに対して本自動調整機能を適用した場合には、上限をＶＵ２とし下限をＶＬ２とする第２のダイナミックレンジ内の輝度のみを対象として、対象期間Ｔ内の各フレームに対応する超音波画像が生成される。従って、図１０から図１４のヒストグラムに示すように、経時的に輝度が大きく変化する場合であっても、図１０から図１４の超音波画像に示すように、ダイナミックレンジが最適化された超音波動画像を生成・表示することができる。
【００４７】
一般に、造影剤エコー法では、映像化対象に造影剤が流入するまでに時間が掛かる。また、例えばフラッシュ送信とモニタリング送信とを繰り返し、造影剤が映像化対象に再間環流する様子を経時的に観察する場合等もある。従って、画像診断に有益な動画像対して本ダイナミックレンジの自動調整機能を適用する場合には、対象期間Ｔの開始時刻Ｔ１及び終了時刻Ｔ２を好適に決定することが重要である。
【００４８】
そのため、本超音波診断装置１０では、例えば以下の各手法によって対象期間Ｔ、開始時刻Ｔ１、終了時刻Ｔ２を決定することができる。
【００４９】
［実施例１］
対象期間Ｔの開始時刻Ｔ１、終了時刻Ｔ２を、図２に示したスタート／ストップスイッチ１３１がマニュアル操作されたタイミングに基づいて決定することができる。すなわち、制御プロセッサ２５は、ダイナミックレンジの自動調整モード下においてスタート／ストップスイッチ１３１が最初に押されたタイミングを開示時刻Ｔ１とし、二回目に押されたタイミングを終了時刻Ｔ２とする。このＴ１〜Ｔ２によって定義される対象期間Ｔ内のフレームに対応する画像データに対して、本ダイナミックレンジの自動調整処理が実行される。
【００５０】
［実施例２］
対象期間Ｔの終了時刻Ｔ２を、現在撮影中の超音波画像の取り込み動作（フリーズボタン操作等）タイミングに基づいて決定することができる。すなわち、制御プロセッサ２５は、ダイナミックレンジの自動調整モード下においてフリーボタンが押されると、これに連動して、フリーズ対象となったフレームに対応する時刻を終了時刻Ｔ２を決定する。このとき、開示時刻Ｔ１は、どの様な手法によって設定してもよい。
【００５１】
［実施例３］
対象期間Ｔの開始時刻Ｔ１を、造影剤注入開始におけるタイマーＯＮ動作タイミングを基準として決定することができる。すなわち、制御プロセッサ２５は、造影剤注入開始におけるタイマーＯＮと連動して、当該タイマーＯＮ時刻の所定期間後の時刻を開始時刻Ｔ１を決定する。このとき、終了時刻Ｔ２は、どの様な手法によって設定してもよい。この手法は、造影剤による輝度の時間的変化を観察する場合等に好適である。
【００５２】
［実施例４］
対象期間Ｔの開始時刻Ｔ１及び終了時刻Ｔ２の少なくとも一方を、フラッシュ送信タイミングを基準として決定することができる。すなわち、制御プロセッサ２５は、例えば第ｎ回目（ｎは自然数）のフラッシュ送信と連動して開始時刻Ｔ１を決定する。また、制御プロセッサ２５は、例えば第ｎ＋１回目のフラッシュ送信と連動して終了開始時刻Ｔ１を決定する。この手法は、造影剤再環流の様子を観察する場合等に好適である。
【００５３】
［実施例５］
時間的に連続する画像データに対して、開示時刻Ｔ１及び終了時刻Ｔ２を事後的に決定することができる。すなわち、制御プロセッサ２５は、予め取得された時間的に連続する画像データに対して、所定の手法により設定された開示時刻Ｔ１から終了時刻Ｔ２までに属するフレームを切り出し、これに対して本ダイナミックレンジの自動調整処理を実行する。
【００５４】
（動作）
次に、ダイナミックレンジの自動調整処理を含む典型的なスキャンシーケンスについて説明する。
【００５５】
図１５は、ダイナミックスタディ（血流の動的観察）を行う場合において、本ダイナミックレンジの自動調整機能を適用する場合の処理の流れを示したフローチャートである。同図において、まず、造影剤投与が開始されると、制御プロセッサ２５は、これに連動して時間経過を把握するためのタイマーをＯＮにし（ステップＳ１）、当該タイマーＯＮ時刻の所定期間後の時刻を開始時刻Ｔ１を決定する（ステップＳ２）。
【００５６】
次に、制御プロセッサ２５は、造影剤が注入された被検体を超音波によって走査し、第１のダイナミックレンジを用いて複数フレームに対応する複数の第１の画像データを取得する（ステップＳ３）。このとき、ユーザは、関心領域（例えば疑腫瘍）の染影の時間変化を第１のダイナミックレンジによって、所定時間（例えば、ｗａｓｈ  ｉｎを見るのであれば６０秒程度）観察することができる。
【００５７】
次に、ユーザから所望のタイミングによるフリーズボタンの操作を受けると（ステップＳ４）、制御プロセッサ２５は、当該フリーズボタン操作と連動して、対象期間Ｔの終了時刻Ｔ２を決定する（ステップＳ５）。この終了時刻Ｔ２及びステップＳ２における開始時刻Ｔ１の決定により、第２のダイナミックレンジが定義されることになる。なお、フリーズボタンの操作後であっても、必要に応じて、被検体の臓器に残っている造影剤を観察するために、スキャンを実行するようにしてもよい。
【００５８】
次に、制御プロセッサ２５は、既述のダイナミック自動調整処理を実行して第２のダイナミックレンジに基づいて輝度変換を行った画像データを取得し（ステップＳ６）、各画像を自動ループ再生した後、ユーザからの指定に基づいて動画データ或いは静止画データを記憶ユニット２６に記憶する（ステップＳ７）。
【００５９】
また、図１６は、再環流の観察を行う場合において、本ダイナミックレンジの自動調整機能を適用する場合の処理の流れを示したフローチャートである。同図において、まず、造影剤投与が開始されると、制御プロセッサ２５は、これに連動して時間経過を把握するためのタイマーをＯＮにする（ステップＳ１１）。
【００６０】
次に、制御プロセッサ２５は、造影剤を使って再還流の様子を見る為に、モニタリング送信（造影剤バブルを破壊しない音圧の超音波を連続的に送信し、造影剤の流入状況をモニタリングするための送信）及びフラッシュ送信（造影剤バブルを破壊可能な音圧の超音波を送信し、走査面から一時的に造影剤の流入をリセットするための送信）を含むスキャンシーケンスに従って、造影剤が注入された被検体を超音波によって走査し、第１のダイナミックレンジを用いて複数フレームに対応する複数の第１の画像データを取得する（ステップＳ１２）。また、制御プロセッサ２５は、フラッシュ送信と連動して開始時刻Ｔ１を決定する（ステップＳ１３）。
【００６１】
次に、例えば輝度がプラトーに達した後のストップボタン１３１の操作を受けると（ステップＳ１４）、制御プロセッサ２５は、当該フリーズボタン操作と連動して、対象期間Ｔの終了時刻Ｔ２を決定する（ステップＳ１５）。この終了時刻Ｔ２及びステップＳ１３における開始時刻Ｔ１の決定により、第２のダイナミックレンジが定義されることになる。
【００６２】
次に、制御プロセッサ２５は、既述のダイナミック自動調整処理を実行して第２のダイナミックレンジに基づく画像データを取得し（ステップＳ１６）、各画像を自動ループ再生した後、ユーザからの指定に基づいて動画データ或いは静止画データを記憶ユニット２６に記憶する（ステップＳ１７）。
【００６３】
なお、必要に応じて、造影剤一回の投与で、図１５に示すシーケンスと図１６に示すシーケンスを組み合わせて実行するようにしてもよい。
【００６４】
（効果）
以上述べた構成によれば、次の効果を実現することができる。
【００６５】
本超音波診断装置は、複数フレームに対応する画像データの各輝度を用いて統計的処理を行い、被検体からのエコー信号に対応する輝度に基づいてダイナミックレンジの上限を自動的に決定し、また、ノイズレベルに対応する輝度に基づいて同のダイナミックレンジの下限値を自動的に決定するダイナミックレンジの自動調整機能を有する。従って、信号輝度レベルが解りづらい造影剤エコー法において、マニュアル操作等により最初に設定したダイナミックレンジが広すぎたり狭すぎたりする場合であっても、本機能を用いてダイナミックレンジを最適な設定に自動調整することができる。
【００６６】
また、特に造影剤エコー法では、造影剤の流入に従って画像が暗→明→暗と変化する。このように画像上の輝度分布が動的に変動する場合であっても、上記ダイナミックレンジの自動調整機能を用いることで、画像の変化に追従してダイナミックレンジを動的に調整することができる。
【００６７】
（第２の実施形態）
次に、本発明の第２の実施形態について説明する。一般に、あるダイナミックレンジを設定した場合には、その設定したダイナミックレンジで輝度変調が発生し、画像コントラストは、ダイナミックレンジ設定前とは異なるものとなる。第１の実施形態では、造影剤エコー法において、既述の自動調整機能により第２のダイナミックレンジを設定して輝度変調させ、好適なコントラストを有する画像を提供する例を示した。
【００６８】
しかしながら、患者、観察部位の個体差のため、画像内の輝度値の分布にバラツキがある場合がある。また、診断目的によっては、種々のダイナミックレンジで観察することが有効である場合がある。これらの場合には、第１の実施形態に係る自動調整機能により設定されたダイナミックレンジを微調整できることが好ましい。
【００６９】
そこで、本実施形態は、第２のダイナミックレンジ設定後において、コントラストが希望する値となるように微調整することができる超音波診断装置について説明する。
【００７０】
図１７は、本実施形態に係るダイナミックレンジの微調整機能に従う処理（ダイナミックレンジの微調整処理）の流れを示したフローチャートである。以下、同図に従ってダイナミックレンジの微調整処理について説明する。なお、このダイナミックレンジの微調整処理は、例えば図１５のステップＳ６、図１６のステップＳ１６においてダイナミックレンジが設定された後に実行される。
【００７１】
まず、制御プロセッサ２５は、第２のダイナミックレンジの上限ＶＵ２と下限ＶＬ２とを設定し（ステップＳ６１）、これを用いて対象期間に含まれる全ての画像を用いて、輝度に関するヒストグラムを作成する（ステップＳ６２）。
【００７２】
次に、制御プロセッサ２５は、作成されたヒストグラムの平均値（平均輝度値）を計算し（ステップＳ６３）、得られた平均輝度値が予め設定された希望値と同じか否か（或いは、平均輝度値が希望値を規準とする所定の範囲にあるか否か）を判定する（ステップＳ６４）。
【００７３】
ステップＳ６４において、平均輝度値が予め設定された希望値と異なる（或いは、平均輝度値が希望値を規準とする所定の範囲内にない）と判定された場合には、制御プロセッサ２５は、第２のダイナミックレンジの上限ＶＵ２を変更し（ステップＳ６５）、新たな第２のダイナミックレンジの上限ＶＵ２を用いて、ステップＳ６１〜Ｓ６４の処理を繰り返す。
【００７４】
一方、平均輝度値が予め設定された希望値と同じ（或いは、平均輝度値が希望値を規準とする所定の範囲内にある）と判定された場合には、現在の上限ＶＵ２を用いて第２のダイナミックレンジを決定し、ステップＳ７、又はステップＳ１７の処理を実行する。なお、希望値は、装置の初期設定値（推奨値）、マニュアル操作によって設定された値等、その設定形態に限定はない。
【００７５】
なお、上記ステップＳ６２において、輝度値に関するヒストグラムを各画像毎に作成し、ステップＳ６４において、各フレーム毎の平均輝度値を計算し、これを用いて全体の平均値を計算するようにしてもよい。また、上記ダイナミックレンジの微調整処理は、フレーム毎に行うことも可能である。さらに、上記ダイナミックレンジの微調整処理を、対象期間Ｔに含まれる画像のうち最大の輝度値を有する画像についてのみ実行するようにしてもよい。
【００７６】
また、上記ダイナミックレンジの微調整処理においては、平均輝度値が希望値と一致等するように、上限値を変更しながら調整する場合を例とした。しかしながら、これに拘泥されず、例えば所定の計算により、平均輝度値と希望値とを一致等させる上限値を取得するようにしてもよい。
【００７７】
以上述べた構成によれば、ダイナミックレンジを微調整することができ、患者、観察部位の個体差等に関わらず、ユーザにとって最適なダイナミックレンジを設定することができる。また、極端に高い値を有するノイズが発生した場合であっても、安定して高画質の超音波画像を提供することができる。
【００７８】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。具体的な変形例としては、例えば次のようなものがある。
【００７９】
（１）本実施形態に係る各機能は、当該処理を実行するプログラムをワークステーション等のコンピュータにインストールし、これらをメモリ上で展開することによっても実現することができる。このとき、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスク（フロッピー（登録商標）ディスク、ハードディスクなど）、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤなど）、半導体メモリなどの記録媒体に格納して頒布することも可能である。
【００８０】
（２）また、本実施形態に係るダイナミックレンジの自動調整機能を超音波画像処理装置で実現する場合には、ダイナミックレンジの下限ＶＬ２に関する情報として、当該下限ＶＬ２の値のみならず、例えば下限ＶＬ２の決定に用いられた画像データを記憶するようにし、これを用いて事後的に下限ＶＬ２を決定するようにしてもよい。また、対象期間Ｔ、開始時刻Ｔ１、終了時刻Ｔ２に関する情報、上限ＶＵ２に関する情報、下限ＶＬ２に関する情報は、当該ダイナミックレンジの自動調整機能の対象となる画像の付帯情報として管理され、記憶されることが好ましい。さらに、これらの情報は、必要に応じて事後的に変更し再設定することも可能である。
【００８１】
（３）本ダイナミックレンジの自動調整機能は、造影剤エコー法において特に効果を発揮するものであるが、その適用は、これに拘泥されない。例えば、造影剤を用いないで超音波動画像を観察する場合等においても、本自動調整機能によってダイナミックレンジを最適化することによっても、好適な超音波画像を提供することができる。
【００８２】
また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【産業上の利用可能性】
【００８３】
以上本発明によれば、超音波イメージングにおいて、事前に予測することが困難なダイナミックレンジを最適な設定に自動調整することが可能な超音波診断装置又は超音波画像処理装置及び超音波画像処理プログラムを実現することができる。
【図面の簡単な説明】
【００８４】
【図１】図１は、本実施形態に係る超音波診断装置１０のブロック構成図を示している。
【図２】図２は、ダイナミックレンジ自動調整機能のＯＮ／ＯＦＦ指示を入力するための専用ボタン等を示した図である。
【図３】図３は、ダイナミックレンジの自動調整処理の流れを示したフローチャートである。
【図４】図４は、第１のダイナミックレンジに従う輝度に関するヒストグラムの一例を示した図である。
【図５】図５は、本ダイナミックレンジの自動調整処理により調整される第２のダイナミックレンジを説明するための図である。
【図６】図６は、第１のダイナミックレンジに従って取得された超音波画像の一例を示した図である。
【図７】図７は、本ダイナミックレンジの自動調整処理により調整された第２のダイナミックレンジに従って取得された超音波画像の一例を示した図である。
【図８】図８は、第１のダイナミックレンジに従って取得された超音波画像の他の例を示した図である。
【図９】図９は、各フレームにおける最大輝度（階調）Ｘｍａｘと最小輝度（階調）Ｘｍｉｎ（すなわち、最大輝度Ｘｍａｘと最小輝度Ｘｍｉｎの時間変化）を示した図である。
【図１０】図１０は、輝度が時間変化する超音波動画像に対して、本ダイナミックレンジの自動調整処理を適用した場合の効果を説明するための図である。
【図１１】図１１は、輝度が時間変化する超音波動画像に対して、本ダイナミックレンジの自動調整処理を適用した場合の効果を説明するための図である。
【図１２】図１２は、輝度が時間変化する超音波動画像に対して、本ダイナミックレンジの自動調整処理を適用した場合の効果を説明するための図である。
【図１３】図１３は、輝度が時間変化する超音波動画像に対して、本ダイナミックレンジの自動調整処理を適用した場合の効果を説明するための図である。
【図１４】図１４は、輝度が時間変化する超音波動画像に対して、本ダイナミックレンジの自動調整処理を適用した場合の効果を説明するための図である。
【図１５】図１５は、ダイナミックスタディ（血流の動的観察）を行う場合において、本ダイナミックレンジの自動調整機能を適用する場合の処理の流れを示したフローチャートである。
【図１６】図１６は、再環流の観察を行う場合において、本ダイナミックレンジの自動調整機能を適用する場合の処理の流れを示したフローチャートである。
【図１７】図１７は、第２の実施形態に係るダイナミックレンジの微調整機能に従う処理（ダイナミックレンジの微調整処理）の流れを示したフローチャートである。
【符号の説明】
【００８５】
１０…超音波診断装置、１２…超音波プローブ、１３…入力装置１３、１４…モニター、２１…送受信ユニット、２２…Ｂモード処理ユニット、２３…ドプラ処理ユニット、２４…画像生成ユニット、２５…制御プロセッサ（ＣＰＵ）、２６…記憶ユニット、２７…画像メモリ、２８…ソフトウェア格納ユニット、２９…インターフェースユニット、３０…インタフェース部
【図１】
2008264530.tif 000003
【図２】
2008264530.tif 000004
【図３】
2008264530.tif 000005
【図４】
2008264530.tif 000006
【図５】
2008264530.tif 000007
【図６】
2008264530.tif 000008
【図７】
2008264530.tif 000009
【図８】
2008264530.tif 000010
【図９】
2008264530.tif 000011
【図１０】
2008264530.tif 000012
【図１１】
2008264530.tif 000013
【図１２】
2008264530.tif 000014
【図１３】
2008264530.tif 000015
【図１４】
2008264530.tif 000016
【図１５】
2008264530.tif 000017
【図１６】
2008264530.tif 000018
【図１７】
2008264530.tif 000019