JP6169361B2 - 超音波診断装置および輝度補正方法 - Google Patents

Description

本発明は、被検体内に超音波を送信し反射して戻ってきた超音波を受信して受信信号を得、その受信信号に基づいて被検体内の断層像を表示する超音波診断装置、およびその断層像の輝度を補正する輝度補正方法に関する。

超音波診断装置において、被検体内の深部で反射して戻ってきた超音波は概ね弱くしたがって受信信号は小さく、この受信信号をそのまま輝度信号に変換して断層像を表示すると被検体内の深い箇所ほど暗い断層像が表示されることになる。

この断層像の輝度を補正するため、特許文献１には、断層像を被検体内の深さに応じて複数の深さ領域に分割したときの各深さ領域ごとの輝度の平均値を求め、それらの平均値がほぼ同一となるように輝度補正することが提案されている。

また、深さ領域ごとの輝度が互いにほぼ同一となるように補正すると、スペックルノイズ等のノイズまで強調されてかえって見にくい断層像となるおそれがあり、特許文献２には各深さ領域ごとの信号値の分散値を算出し、分散値がある値以上の深さ領域について輝度を下げることが提案されている。

特開平６−５４８４９号公報 特開２００５−１５２４２２号公報

特許文献２にて提案された、分散値がある値以上の深さ領域の輝度を下げる提案を採用すると、分散値がある値以上の深さ領域については特許文献１で提案された、平均値をほぼ同一とすることにより強調されたノイズをある程度緩和することができる。しかしながらこの特許文献２での提案では効果は限定的であり、分散値がある値以下の深さ領域については、深さ領域ごとのノイズレベルに差があっても改善されないという問題がある。

また、単純に分散値を用いるだけでは診断部位等によってはかえってノイズを増大させるおそれもある。

本発明は、上記事情に鑑み、ノイズを有効に抑制することにより見やすい断層像が表示される超音波診断装置および輝度補正方法を提供することを目的とする。

上記目的を達成する本発明の超音波診断装置は、
被検体表面に当てがった超音波探触子から被検体内に超音波を送信し被検体内で反射して戻ってきた超音波を超音波探触子で受信して受信信号を得、その受信信号に基づいて被検体内の断層像を表示画面上に表示する超音波診断装置において、
補正モードで得られた受信信号由来の信号に基づいて、被検体内の、超音波探触子が当てがわれた被検体表面からの深さに応じて分割された複数の深さ領域それぞれについて各深さ領域内のピクセル毎の信号値の分散値を算出するとともに複数の深さ領域を統合した全領域についてのピクセル毎の信号値の分散値を算出して、全領域についての分散値に対する複数の深さ領域それぞれの分散値の比である正規化分散値を算出する正規化分散値算出部と、
正規化分散値を変数とし、正規化分散値の増大につれて輝度を下げる補正関数値から輝度を上げる補正関数値まで単調増加する連続曲線からなる補正関数を参照し、複数の深さ領域それぞれについての正規化分散値を複数の深さ領域それぞれの補正関数値に変換して、補正関数値由来の補正値を算出する補正値算出部と、
補正モードよりも後の表示モードにおける受信信号由来の信号を深さ領域ごとに深さ領域に対応する補正値で補正することにより、受信信号に基づいて表示画面上に表示される断層像の、深さ領域ごとの輝度を補正する輝度補正部と、
輝度補正部で補正された後の信号に基づいて、表示画面上に深さ領域ごとに輝度が補正された断層像を表示する断層像表示部とを備えたことを特徴とする。

本発明の超音波診断装置は、上記の正規化分散値を算出している。このため断層像全体の輝度が高い場合であっても低い場合であっても安定した指標を得ることができる。

さらに本発明の超音波診断装置は、その正規化分散値を変数とし、その正規化分散値の増大につれて輝度を下げる補正関数値から輝度を上げる補正関数値まで単調増加する連続曲線からなる補正関数を参照して補正値を算出するため、全ての深さ領域について各深さ領域のノイズレベルに応じた輝度補正が行なわれ、断層像全域の画質が向上する。

ここで、本発明の超音波診断装置において、補正モードで得られた受信信号由来の信号に基づいて全領域についてのピクセル毎の信号値のヒストグラムを算出しそのヒストグラムに基づいて基準信号値を算出して基準信号値よりも大きな信号値を持つピクセル数と小さな信号値を持つピクセル数との間のピクセル数どうしの比率を算出する比率算出部をさらに備え、補正値算出部が、複数の深さ領域それぞれについて、補正関数値に由来するとともにさらに上記の比率にも由来する補正値であって、その比率に基づいては基準信号値よりも小さな信号値を持つピクセル数の比率が大きいほどその深さ領域の輝度を下げる向きに補正する補正値を算出するものであることが１つの好まし形態である。

上記の、補正関数値に由来する補正値を採用すると、例えば頸動脈を診断部位とする場合は血管壁からの大きな反射信号（高輝度）と血管内からの極く弱い反射信号（低輝度）とが存在する断層像の場合は、断層像上での血管壁の面積（高輝度の領域の面積）が小さくても大きな分散値を得ることができ、したがって高い輝度で見やすい表示が行なわれることになる。これに対し、上記の比率に由来する補正値では、輝度の低い大面積の領域の影響を受けて基準信号値よりも小さな信号値を持つピクセル数が優って上記の比率が下がり、低い輝度で表示されて見にくい断層像となるおそれがある。

一方、生体組織からの反射信号が大部分の面積を占める断層像の場合は、上記の比率に由来する補正値を採用することにより各深さ領域のノイズレベルが確実に反映された補正値を得ることができ、見易い断層像となる。これに対し、このような断層像の場合、上記の規格化分散値には差異が生じにくく、したがって上記の補正関数値由来の補正値を採用しても十分な補正が行なわれないおそれがある。

そこで、上記の好ましい形態の通り、正規化分散値算出部を備えるとともに比率算出部をさらに備え、複数の深さ領域それぞれについて、補正関数値に由来するとともにさらに上記の比率にも由来する補正値であって、上記の比率に基づいては基準信号値よりも小さな信号値を持つピクセル数の比率が大きいほどその深さ領域の輝度を下げる向きに補正する補正値を算出すると、補正関数値由来の補正値と上記の比率由来の補正値の双方の長所を生かすとともに互いの短所を補った最終的な補正値を算出することができ、断層像の画質をさらに向上させることができる。

ここで、この好ましい形態において、被検体の診断部位を表わす診断部位情報又は超音波探触子の種別を表わす種別情報を取得する情報取得部をさらに備え、補正値算出部が、上記の診断部位情報又は種別情報に応じて、補正関数値由来の補正値と上記の比率由来の補正値との重み付けを変更してそれらの補正値を組み合わせることにより最終的な補正値を算出するものであってもよい。この構成を採用すると、診断部位や超音波探触子の種別に応じて上記の２種類の補正方法が適応的に組み合わされた補正値が算出され、断層像の画質を一層向上させることができる。

また、本発明の超音波診断装置において、補正モードで得られた受信信号由来の信号に基づいて、複数の深さ領域それぞれについてのピクセル毎の信号値の平均値を算出する平均値算出部をさらに備え、補正値算出部が、複数の深さ領域それぞれについて、さらに平均値にも由来する補正値であって、平均値に基づいては複数の深さ領域相互の輝度が相互に近づく向きに補正する補正値を算出するものであることが好ましい。

上記の補正関数値や上記の比率に由来するとともに各深さ領域ごとの平均値にも由来する補正値を算出してその補正値を適用することにより、断層像の画質を一層向上させることができる。

また、本発明の輝度補正方法は、被検体表面に当てがった超音波探触子から被検体内に超音波を送信し被検体内で反射して戻ってきた超音波を超音波探触子で受信して受信信号を得、その受信信号に基づいて被検体内の断層像を表示画面上に表示する超音波診断装置における断層像の輝度補正方法であって、
補正モードで得られた受信信号由来の信号に基づいて、被検体内の、超音波探触子が当てがわれた被検体表面からの深さに応じて分割された複数の深さ領域それぞれについて各深さ領域内のピクセル毎の信号値の分散値を算出するとともに複数の深さ領域を統合した全領域についてのピクセル毎の信号値の分散値を算出して、全領域についての分散値に対する複数の深さ領域それぞれの分散値の比である正規化分散値を算出し、
正規化分散値を変数とし、正規化分散値の増大につれて輝度を下げる補正関数値から輝度を上げる補正関数値まで単調増加する連続曲線からなる補正関数に基づいて、複数の深さ領域それぞれについての正規化分散値を複数の深さ領域それぞれの補正関数値に変換して、補正関数値由来の補正値を算出し、
補正モードよりも後の表示モードにおける受信信号由来の信号を深さ領域ごとに深さ領域に対応する補正値で補正することにより、受信信号に基づいて表示画面上に表示される断層像の、深さ領域ごとの輝度を補正することを特徴とする。

本発明の超音波診断装置および輝度補正方法によれば、ノイズが有効に抑制された見やすい断層像が得られる。

本発明の一実施形態としての超音波診断装置の構成を示すブロック図である。 入力部の作用説明図である。 図１に１つのブロックで示すゲイン／ＳＴＣ補正値算出部での演算処理フローを示した図である。 、図３に示すゲイン／ＳＴＣ補正値算出部での演算処理内容の説明図である。

以下、本発明の実施の形態を説明する。

図１は、本発明の一実施形態としての超音波診断装置の構成を示すブロック図である。

ここでは、この図１に示す本発明の一実施形態としての超音波診断装置について説明することで、本発明の輝度補正方法の実施形態の説明を兼ねるものとする。

この図１に示す超音波診断装置１００には、超音波探触子１が備えられている。この超音波探触子１は、着脱自在に交換され、生体３０の診断部位や診断内容等に応じて、それに適した超音波探触子が用いられる。この超音波探触子１には生体３０の体表に当てがわれる側の先端に圧電セラミックス等の振動子（図示せず）が配列されている。また、この超音波探触子１には、送信部２と受信部３が接続される。

送信部２には、送信の繰り返し周期（例えば４ＫＨｚ）を与えるレートパルスを発生するパルス発生器が含まれている。この送信部２は、例えば６４チャンネルの、パルスドライバ及び遅延回路を有する。パルスドライバはレートパルスのタイミングで送信周波部（例えば２．５ＭＨｚ）に等しい周期の振動パルスを発生し、超音波探触子１の振動子に印加する。遅延回路は超音波ビームを収束し、かつ指向性を与えるために各チャンネル毎のパルス発生タイミングに所定の遅延を与える。その結果、超音波ビームが指向性に応じた方向にパルス放射される。このようにして、超音波探触子１から、レートパルス周期で、生体３０内に延びる、例えば６４本の走査線に沿う各方向に向けて超音波探ビームが順次に送波されて、１フレーム分の走査が完了する。

一方、生体３０内の音響インピーダンスの不連続面で反射したエコーは超音波探触子１を介して受信部３でチャンネル毎に受信される。受信部３は、プリアンプ、遅延回路、および加算回路から構成される。受信信号は、プリアンプで増幅され、遅延回路により各チャンネル毎に所定の遅延が与えられて、加算回路により加算される。これにより、各走査線について順次に、超音波ビームが送波された方向からのエコーが受信される。受信部３から出力された受信信号は、対数圧縮部４で対数圧縮され、対数圧縮後の受信信号がゲイン／ＳＴＣ（ｓｅｎｓｉｔｉｖｉｔｙ ｔｉｍｅ ｃｏｎｔｒｏｌ）補正部５に入力される。ゲイン、ＳＴＣの詳細は後述する。このゲイン／ＳＴＣ補正部５では、入力部１０から入力されたゲイン／ＳＴＣ補正値、又はゲイン／ＳＴＣ補正値算出部２０で算出されたゲイン／ＳＴＣ補正値に応じて受信信号の信号値を補正するものである。このゲイン／ＳＴＣ補正部５での補正結果は、最終的に、表示部９に備えられた表示画面上に表示される断層像の輝度に反映される。

ここで、入力部１０は、オペレータによる操作を受けそのオペレータによる操作が反映されたゲイン用のエンコーダ１１およびＳＴＣ用のエンコーダ１２の値をゲイン／ＳＴＣ補正値としてゲイン／ＳＴＣ補正部５に入力するものである。

ゲイン／ＳＴＣ補正部５でゲイン／ＳＴＣ補正が行なわれた後の受信信号はダイナミックレンジ調整部６によりダイナミックレンジ調整、すなわちコントラスト調整を受け、さらに信号処理部７によりノイズ除去処理やエッジ強調処理等の各種処理を受けて、デジタルスキャンコンバータ（ＤＳＣ）８に入力される。このＤＳＣ８では、各走査線１ａごとの受信信号の、１フレーム分の集合を、表示部９の表示画面に表示するのに適したラスタ走査信号へと信号形式を変換する処理が行なわれる。このＤＳＣ８で信号形式の変換により生成された、表示部９の表示画面に表示するのに適した信号は表示部９に入力され、その表示部９の表示画面上に、ＤＳＣ８から送られてきた信号に基づく断層像が表示される。

図２は、入力部の作用説明図である。

図２（Ａ）はゲインについての説明図である。図１に示す入力部１０には、オペレータにより操作される、ゲイン調整用のスライドつまみ１１ａが備えられている。

オペレータがこのスライドつまみ１１ａを動かす操作を行なうとその操作後のスライドつまみ１１ａの位置がこの入力部１０に備えられているゲイン用のエンコーダ１１（図１参照）で読み取られてその読取値がゲイン補正値としてゲイン／ＳＴＣ補正部５に入力される。このゲイン補正値は、断層像Ｉの全域に作用し、その断層像全域の輝度が、そのゲイン補正値に応じて上下する。

また、図２（Ｂ）は，ＳＴＣについての説明図である。図１に示す入力部１０には、ゲイン調整用のスライドつまみ１１ａと同様に、オペレータにより操作される、ＳＴＣ調整用の、本実施形態では６つのスライドつまみ１２ａ，１２ｂ，・・・，１２ｆが備えられている。

オペレータがこれらのスライドつまみ１２ａ，１２ｂ，・・・，１２ｆを動かす操作を行なうとその操作後のスライドつまみ１２ａ，１２ｂ，・・・，１２ｆの各位置がこの入力部１０に備えられているＳＴＣ用のエンコーダ１２（図１参照）で読み取られて、それらの読取値がＳＴＣ補正値としてゲイン／ＳＴＣ補正部５に入力される。各スライドつまみ１２ａ，１２ｂ，・・・，１２ｆは、断層像Ｉを超音波探触子１が当てがわれた生体３０の表面からの深さに応じて複数（本実施形態では６つ）のＳＴＣ深度領域ａ，ｂ，ｃ，ｄ，ｅ，ｆに分割したときの、それら複数のＳＴＣ深度領域ａ，ｂ，ｃ，ｄ，ｅ，ｆのそれぞれに対応している。すなわちこれらのスライドつまみ１２ａ，１２ｂ，・・・，１２ｆの位置の読取値により、対応するＳＴＣ深度領域ａ，ｂ，ｃ，ｄ，ｅ，ｆの輝度がそれぞれ調整される。ただし、隣接するＳＴＣ深度領域の境界で輝度がステップ的に変化しないよう、その境界領域では一方のＳＴＣ深度領域の輝度からもう一方のＳＴＣ深度領域の輝度へと滑らかに変化するように調整される。

本実施形態では、図１に示すゲイン／ＳＴＣ補正部５では、図２を参照して説明した入力部１０からのゲイン／ＳＴＣ補正値と、以下に説明する、ゲイン／ＳＴＣ補正値算出部２０で算出されたゲイン／ＳＴＣ補正値が択一的に使用される。いずれのゲイン／ＳＴＣ補正値を使用するかは、オペレータによる、図示しない操作子の操作により切り替えられる。

図３は、図１に１つのブロックで示すゲイン／ＳＴＣ補正値算出部での演算処理フローを示した図である。また図４は、図３に示すゲイン／ＳＴＣ補正値算出部での演算処理内容の説明図である。

ここで説明している超音波診断装置１００は、補正モードと表示モードとを有する。

補正モードは、１フレーム分の受信信号を得てその１フレーム分の受信信号に基づいてゲイン補正値およびＳＴＣ補正値を算出するモードである。これらゲイン補正値やＳＴＣ補正値の算出は、図１に示すゲイン／ＳＴＣ補正値算出部２０がその役割を担っている。

表示モードは、順次得られる複数フレームの受信信号それぞれを、補正モードで算出されたゲイン補正値やＳＴＣ補正値で補正して、その補正後の受信信号に基づく断層像を表示部９の表示画面上に表示するモードである。補正モードは、この超音波診断装置１００の電源を立上げ超音波探触子１を生体３０に当てがって動作を開始した時点、超音波探触子１を交換して動作を開始した時点、および動作中にオペレータが図示しない操作子を操作して新たな補正を指示した時点など、様々なタイミングで実行される。

補正モードに入ると、超音波探触子１から１フレーム分の超音波の送受信が行なわれて１フレーム分の受信信号が生成され、その１フレーム分の受信信号に基づいて、先ず、各ＳＴＣ深度領域ａ〜ｆ（図２参照）のそれぞれについて、各ピクセル毎の信号値の平均値が算出される（ステップＳ０１）。図４（Ａ）は、この平均値演算により算出された各平均値（縦軸）の一例を各ＳＴＣ深度領域（横軸）ごとに示したグラフである。

図１に示すゲイン／ＳＴＣ補正値算出部２０における、図３のステップＳ０１の処理は、本発明にいう平均値算出部の一例に相当する。

この平均値演算に次いで、各ＳＴＣ深度領域ごとの深度補正値が算出される（ステップＳ０２）。

図４（Ｂ）は、各ＳＴＣ深度領域ごとの深度補正値のイメージである。

この深度補正値を受信信号に作用させると、受信信号は、断層像上の各ＳＴＣ深度領域相互の輝度が相互に近づく向きに補正される。本実施形態では、具体的には、各ＳＴＣ深度領域ａ〜ｆを互いに等しい平均輝度に補正する深度補正値が算出される。

深度補正値演算（ステップＳ０２）が終了すると、次の演算の１つとして、ヒストグラムの算出が行なわれる（ステップＳ０３）。ここで算出されるヒストグラムは、ＳＴＣ深度領域ａ〜ｆ（図２参照）を全て合わせた断層像Ｉの全域についての、各ピクセル毎の受信信号の信号値（ピクセル値）のヒストグラムである。すなわち、このヒストグラムは、横軸にピクセル値、縦軸にそのピクセル値の出現頻度を示したグラフ（図示せず）であらわされる。ただし、ここで算出されるヒストグラムは、ステップＳ０２で算出された深度補正値を用いて、この補正モードで得られた１フレーム分の受信信号を、断層像Ｉの全域が同一輝度となるように補正した後の受信信号に基づくヒストグラムである。

このように補正した後の受信信号に基づくヒストグラムを算出することにより、複数のＳＴＣ深度領域ａ〜ｆがヒストグラムに平等に反映されることになる。

この算出されたヒストグラムは、ゲイン補正値算出処理（ステップＳ０４）と、後述するＳＴＣノイズ補正値演算処理（ステップＳ０５）に反映される。

ステップＳ０４のゲイン補正値算出処理では、ステップＳ０３で算出されたヒストグラムの頻度（縦軸）のピークピクセル値（横軸）が診断部位又は超音波探触子の種別に応じた基準値となるように補正するためのゲイン補正値が算出される。

図１に示す本実施形態の超音波診断装置１００の場合、診断部位ごとに超音波探触子１の種別が定められており、超音波探触子１を装着するとその装着された超音波探触子１の種別が自動認識されるように構成されている。また、オペレータによる操作で、診断部位を手動設定することも可能である。

このゲイン補正値算出処理によりゲインの最適化処理が完了する。

ステップＳ０５のＳＴＣノイズ補正値演算処理では、ステップＳ０３で算出されたヒストグラムを使って、以下に説明するＳＴＣノイズ補正値が算出される。

ここでは、ステップＳ０３で算出されたヒストグラムに基づいて基準のピクセル値を算出して、複数のＳＴＣ深度領域ａ〜ｆ（図２参照）のそれぞれについて、その算出された基準のピクセル値よりも大きなピクセル値を持つピクセル数と小さなピクセル値を持つピクセル数との間のピクセル数どうしの比率が算出される。本実施形態では、より具体的には、上記のヒストグラムのピークに対応するピクセル値が基準のピクセル値として採用される。その基準のピクセル値よりも大きなピクセル値を持つピクセルに値‘１’を割り当てるとともにその基準のピクセル値よりも小さなピクセル値を持つピクセルに値‘０’を割り当てて、各ＳＴＣ深度領域ａ〜ｆごとにそのＳＴＣ深度領域に含まれるピクセルについて割り当てられている値‘１’又は‘０’を合計してそのＳＴＣ深度領域に含まれるピクセルの数で割り算をする。これにより、各ＳＴＣ深度領域ごとに基準のピクセル値よりも大きなピクセル値を持つピクセル数と基準のピクセル値よりも小さいピクセル値を持つピクセル数との間のピクセル数どうしの比率が算出される。

図４（Ｃ）は、各ＳＴＣ深度領域ａ〜ｆごとの比率のイメージを表した構成図である。この比率の値が小さいほど、基準のピクセル値よりも小さいピクセル値を持つピクセルの数の割合が大きいことを意味している。

図３ステップＳ０５のＳＴＣノイズ補正値演算処理では、上記の比率を算出した後、さらに、この比率に由来するＳＴＣノイズ補正値が算出される。あるいは、上記の比率自体をＳＴＣノイズ補正値としてもよい。

図４（Ｄ）は、上記の比率に基づいて算出されたＳＴＣノイズ補正値のイメージを表した構成図である。

基準のピクセル値よりも小さいピクセル値を持つピクセルの数の比率が大きい領域、すなわち、上記の‘１’‘０’を割り当てて算出した値が小さい領域には、強いノイズが含まれている場合が多い。また、これと逆に、基準のピクセル値よりも大きいピクセル値を持つピクセルの数の比率が大きい領域、すなわち、上記の‘１’‘０’を割り当てて算出した値が大きい領域はノイズ成分が小さい傾向にある。そこでここでは、ノイズ成分の大きなＳＴＣ深度領域（上記の‘１’‘０’を割り当てて算出した値が小さいＳＴＣ深度領域）ほど輝度を下げる向きに補正するＳＴＣノイズ補正値が算出される。

ただし、上記の比率に基づいて算出された補正値は、１フレーム中で、生体組織からのエコーが大部分の面積を占める断層像の場合は適正な補正が行なわれるが、一方、例えば頸動脈を診断部位とする場合など、生体組織（ここでは血管壁）からの強いエコーは存在するものの、その断層像の広い面積を血管内部の小さいエコーしか存在しない領域が占めるような断層像の場合は、輝度が不必要に下がる傾向となる。すなわち、上記のようにして算出されたＳＴＣノイズ補正値は、万全ではなく苦手な断層像も存在する。

ステップＳ０３のヒストグラム演算およびステップＳ０５のＳＴＣノイズ補正値演算における上記の比率算出処理は、それらを合わせたものが本発明にいう比率算出部の一例に相当する。ステップＳ０５のＳＴＣ補正値算出処理における、上記の比率に基づいてＳＴＣノイズ補正値を算出する処理は、ステップＳ０２における深度補正値演算処理とともに本発明にいう補正値算出部のうちの一部に相当する。

図３のステップＳ０６では、各ＳＴＣ深度領域ごとの分散値が算出される。この分散値算出にあたっても、ステップＳ０２で算出された深度補正値に基づき、複数のＳＴＣ深度領域ａ〜ｆ（図２参照）ごとの平均的な輝度が互いに等しくなるように補正された後の受信信号が採用され、その補正された受信信号に基づいて各ＳＴＣ深度領域ａ〜ｆそれぞれについてのピクセル値の分布の分散値が算出される。さらに、ここでは複数のＳＴＣ深度領域ａ〜ｆを統合した断層像Ｉの全域についての分散値も算出され、全域についての分散値に対する各ＳＴＣ深度領域ごとの分散値の比率（ここではこの比率を「正規化分散値」と称する）が算出される。

このステップＳ０６における演算処理は、本発明にいう正規化分散値算出部の一例に相当する。

図４（Ｅ）は、上記のようにして算出された各ＳＴＣ深度領域ａ〜ｆごとの正規化分散値のイメージを表した模式図である。

次に、ステップＳ０７のＳＴＣ分散補正値演算処理が行なわれる。

図４（Ｆ）は、正規化分散値を変数とし分散補正値を関数値とした分散補正関数の一例を示した模式図である。

図１に示す超音波診断装置１００のゲイン／ＳＴＣ補正値算出部２０には、例えば図４（Ｆ）に示すような分散補正関数が記憶されている。

この分散補正関数は、正規化分散値の増大に伴って輝度を下げる補正関数値から輝度を上げる補正関数値まで単調増加する連続関数からなる補正関数である。これは、正規化分散値が大きい領域ほどノイズ成分が少なく、正規化分散値が小さい領域ほどノイズ成分が大きいという知見に基づくものである。

ただし、診断部位ごとに分散値とノイズレベルとの関係が異なるため、本実施形態では、分散補正関数は、診断部位又は超音波探触子の種別ごとに最適な関数が用意されている。

図３ステップＳ０７では、今回の診断部位又は超音波探触子に応じた、図４（Ｆ）に一例を示すような分散補正関数が参照されて、ステップＳ０６で算出された各ＳＴＣ深度領域ａ〜ｆごとの各正規化分散値が各分散補正値に変換される。

図４（Ｇ）は、このようにして算出された分散補正値のイメージを表した模式図である。

正規化分散値が小さいほどノイズ成分が大きい傾向にあり、したがって輝度を下げる補正値となっている。

図３ステップＳ０８では、ステップＳ０２で算出された各ＳＴＣ深度領域ａ〜ｆごとの深度補正値、ステップＳ０５で算出された、各ＳＴＣ深度領域ごとのノイズ補正値、および、ステップＳ０７で算出された、各ＳＴＣ深度領域ごとの分散補正値に基づいて、各ＳＴＣ深度領域ａ〜ｆごとの最終的なＳＴＣ補正値が算出される。ここでは、一例として、これら３種類の補正値が各ＳＴＣ深度領域ごとに互いに乗算されることにより、最終的なＳＴＣ補正値が算出される。この、各ＳＴＣ深度領域ａ〜ｆごとの最終的なＳＴＣ補正値の算出により、今回の補正モードにおけるＳＴＣ補正値算出処理が終了する。

本実施形態では、ステップＳ０２の深度補正値演算処理、およびステップＳ０５のＳＴＣノイズ補正値演算処理に、さらにステップＳ０７のＳＴＣ分散補正値演算処理およびステップＳ０８のＳＴＣ補正値算出処理を合わせたものが、本発明にいう補正値算出部の一例に相当する。

ステップＳ０４で算出されたゲイン補正値およびステップＳ０８で算出されたＳＴＣ補正値は、図１のゲイン／ＳＴＣ補正部５において、今回の補正モードに続く表示モードにおいて後述した受信信号に適用され、最適な輝度の断層像が得られるようにその受信信号が補正される。

尚、ここでは、ステップＳ０２で算出された深度補正値と、ステップＳ０５で算出されたノイズ補正値と、ステップＳ０７で算出された分散補正値とを単純に乗算して最終的なＳＴＣ補正値を算出する旨、説明したが、ステップＳ０８におけるＳＴＣ補正値算出処理は、単純な乗算に限られるものではない。例えば、診断部位又は超音波探触子の種別に応じた重みを付して重み付け加算により最終的なＳＴＣ補正値を算出してもよい。

あるいは、診断部位又は超音波探触子の種別に応じて、ステップＳ０５で算出されたノイズ補正値とステップＳ０７で算出された分散補正値とを切り替えることにより、それら２種類の補正値のうちの一方の補正値を使って最終的なＳＴＣ補正値を算出してもよい。

さらには、ここで説明した超音波診断装置１００は、ステップＳ０５のノイズ補正値演算処理機能を備えたものであるが、ノイズ補正値と分散補正値のうちのノイズ補正値は採用せずに常に分散補正値を採用する構成としてもよい。

ａ〜ｆ ＳＴＣ深度領域
１ 超音波探触子
２ 送信部
３ 受信部
４ 対数圧縮部
５ ゲイン／ＳＴＣ補正部
６ ダイナミックレンジ調整部
７ 信号処理部
８ デジタルスキャンコンバータ（ＤＳＣ）
９ 表示部
１０ 入力部
１１，１２ エンコーダ
１１ａ，１２ａ〜１２ｆ スライドつまみ
２０ ゲイン／ＳＴＣ補正値算出部
３０ 生体
１００ 超音波診断装置

Claims (4)

1. 被検体表面に当てがった超音波探触子から該被検体内に超音波を送信し該被検体内で反射して戻ってきた超音波を該超音波探触子で受信して受信信号を得、該受信信号に基づいて該被検体内の断層像を表示画面上に表示する超音波診断装置において、
   補正モードで得られた前記受信信号由来の信号に基づいて、前記被検体内の、前記超音波探触子が当てがわれた該被検体表面からの深さに応じて分割された複数の深さ領域それぞれについて各深さ領域内のピクセル毎の信号値の分散値を算出するとともに該複数の深さ領域を統合した全領域についてのピクセル毎の信号値の分散値を算出して、該全領域についての分散値に対する該複数の深さ領域それぞれの分散値の比である正規化分散値を算出する正規化分散値算出部と、
   前記正規化分散値を変数とし、該正規化分散値の増大につれて輝度を下げる補正関数値から輝度を上げる補正関数値まで単調増加する連続曲線からなる補正関数を参照し、前記複数の深さ領域それぞれについての前記正規化分散値を該複数の深さ領域それぞれの補正関数値に変換して、該補正関数値由来の補正値を算出する補正値算出部と、
   前記補正モードよりも後の表示モードにおける前記受信信号由来の信号を前記深さ領域ごとに該深さ領域に対応する補正値で補正することにより、該受信信号に基づいて前記表示画面上に表示される断層像の、該深さ領域ごとの輝度を補正する輝度補正部と、
   前記輝度補正部で補正された後の信号に基づいて、前記表示画面上に前記深さ領域ごとに輝度が補正された断層像を表示する断層像表示部とを備えたことを特徴とする超音波診断装置。
2. 補正モードで得られた前記受信信号由来の信号に基づいて前記全領域についてのピクセル毎の信号値のヒストグラムを算出し該ヒストグラムに基づいて基準信号値を算出して該基準信号値よりも大きな信号値を持つピクセル数と小さな信号値を持つピクセル数との間のピクセル数どうしの比率を算出する比率算出部をさらに備え、
   前記補正値算出部が、前記複数の深さ領域それぞれについて、前記補正関数値に由来するとともにさらに前記比率にも由来する補正値であって、該比率に基づいては前記基準信号値よりも小さな信号値を持つピクセル数の比率が大きいほど当該深さ領域の輝度を下げる向きに補正する補正値を算出するものであることを特徴とする請求項１記載の超音波診断装置。
3. 補正モードで得られた前記受信信号由来の信号に基づいて、前記複数の深さ領域それぞれについてのピクセル毎の信号値の平均値を算出する平均値算出部をさらに備え、
   前記補正値算出部が、前記複数の深さ領域それぞれについて、さらに前記平均値にも由来する補正値であって、該平均値に基づいては該複数の深さ領域相互の輝度が相互に近づく向きに補正する補正値を算出するものであることを特徴とする請求項１または２に記載の超音波診断装置。
4. 被検体表面に当てがった超音波探触子から該被検体内に超音波を送信し該被検体内で反射して戻ってきた超音波を該超音波探触子で受信して受信信号を得、該受信信号に基づいて該被検体内の断層像を表示画面上に表示する超音波診断装置における断層像の輝度補正方法であって、
   補正モードで得られた前記受信信号由来の信号に基づいて、前記被検体内の、前記超音波探触子が当てがわれた該被検体表面からの深さに応じて分割された複数の深さ領域それぞれについて各深さ領域内のピクセル毎の信号値の分散値を算出するとともに該複数の深さ領域を統合した全領域についてのピクセル毎の信号値の分散値を算出して、該全領域についての分散値に対する該複数の深さ領域それぞれの分散値の比である正規化分散値を算出し、
   前記正規化分散値を変数とし、該正規化分散値の増大につれて輝度を下げる補正関数値から輝度を上げる補正関数値まで単調増加する連続曲線からなる補正関数に基づいて、前記複数の深さ領域それぞれについての前記正規化分散値を該複数の深さ領域それぞれの補正関数値に変換して、該補正関数値由来の補正値を算出し、
   前記補正モードよりも後の表示モードにおける前記受信信号由来の信号を前記深さ領域ごとに該深さ領域に対応する補正値で補正することにより、該受信信号に基づいて前記表示画面上に表示される断層像の、該深さ領域ごとの輝度を補正することを特徴とする輝度補正方法。

Images