(57)【要約】      （修正有）【課題】画像の明るさをリアルタイムで調整することが可能な画像処理装置および超音波診断装置を提供する。【解決手段】画像処理装置は、多重解像度分解部１１と、調整部１２と、復元部１４と、を有する。多重解像度分解部１１は、画素値が２次元または３次元の方向に配列された原画像を高周波成分と低周波成分とに分解することを、画素値が各方向に所定個になるまで繰り返すことにより、低周波成分における最下層の画像を生成する。調整部１２は、最下層の画像における各画素の画素値を変更することにより、最下層の画像の明るさを調整する。復元部１４は、明るさが調整された最下層の画像を用いて、画素値が各方向に前記原画像と同じ個数になるように画像を復元する。【選択図】図１
2015083056.tif 000002
【特許請求の範囲】
【請求項１】
画素が２次元または３次元の方向に配列された原画像を高周波成分と低周波成分とに分解することを、前記画素が前記各方向に所定個になるまで繰り返すことにより、前記低周波成分における最下層の画像を生成する多重解像度分解部と、
前記最下層の画像における各画素の画素値を変更することにより、前記最下層の画像の明るさを調整する調整部と、
前記明るさが調整された前記最下層の画像を用いて、前記画素が前記各方向に前記原画像と同じ個数になるように画像を復元する復元部と、
を有する
ことを特徴とする画像処理装置。
【請求項２】
前記調整部は、前記各画素間の画素値の変化を少なくすることにより、前記最下層の画像の明るさを調整するフィルタをさらに有することを特徴とする請求項１に記載の画像処理装置。
【請求項３】
前記フィルタは、各画素値をそれに隣接した画素値の平均値とすることにより、前記明るさを調整するように構成されることを特徴とする請求項２に記載の画像処理装置。
【請求項４】
前記明るさが調整された前記最下層の画像を受けて、ゲイン調整を行うシステムゲイン設定部をさらに有することを特徴とする請求項１から３のいずれかに記載の画像処理装置。
【請求項５】
前記調整部は、前記変更した画素値を調整可能に構成されることを特徴とする請求項１に記載の画像処理装置。
【請求項６】
ユーザーおよび／または検査対象に対応づけられて係数が予め記憶部に記憶され、
前記変更した画素値を調整するとき、前記ユーザーおよび／または検査対象に対応づけられた前記係数に基づいて行われるように構成されることを特徴とする請求項５に記載の画像処理装置。
【請求項７】
前記復元部は、前記明るさが調整された前記最下層の画像に基づいて、前記各方向に前記原画像と同じ個数の画素値が配列された第１の画像を生成し、前記画像を復元するとき、前記第１の画像を用いることを特徴とする請求項１に記載の画像処理装置。
【請求項８】
前記復元部は、前記分解された前記高周波成分の画像に基づいて、前記各方向に前記原画像と同じ個数の画素値が配列された第２の画像を生成し、前記画像を復元するとき、前記第１の画像に加えて、前記第２の画像を用いることを特徴とする請求項７に記載の画像処理装置。
【請求項９】
請求項１から８のいずれかに記載の画像処理装置を有し、
検査対象に超音波を送信し、検査対象からの反射エコー信号を受信し、前記反射エコー信号に基づいて画像を取得し、前記取得した画像を前記原画像として前記画像処理装置に出力することを特徴とする超音波診断装置。
【請求項１０】
前記超音波診断装置によりフリーズ前に取得された画像、および／または、前記超音波診断装置により過去に取得された画像が、前記原画像として前記画像処理装置に出力されるように構成されることを特徴とする請求項９に記載の超音波診断装置。
【請求項１１】
Ｍモード画像、Ｂモード画像から再構成されるＭモード画像、パルスドプラ方式により取得された画像、および、連続波ドプラ方式により取得された画像の少なくとも一つを、前記原画像として前記画像処理装置に出力することを特徴とする請求項９または１０に記載の超音波診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、画像処理装置および超音波診断装置に関する。
【背景技術】
【０００２】
超音波診断装置、Ｘ線ＣＴ装置、およびＭＲＩ装置などのモダリティにより、検査対象が診断され、診断時における検査対象の画像が取得される。
【０００３】
例えば、超音波診断装置は、非侵襲でリアルタイムに検査ができることから多くの医療機関に普及している。超音波診断装置は、超音波パルスを体内に照射する電気音響変換器および画像処理装置を有する。生体に照射された超音波パルスは生体組織の固有音響インピーダンスの差に応じた反射をする。たとえば、インピーダンスの差が大きな境界面でより強い反射信号が得られ、インピーダンスの差が小さな境界面でより弱い反射信号が得られる。反射信号は電気音響変換器によって電気信号に変換される。電気信号は増幅され、断層画像などに再構築されてディスプレイに表示される。ディスプレイに画像が表示されるとき、反射信号が大きなものを明るく、小さなものを暗くなるように、信号振幅に応じた輝度変調が行われる。
【０００４】
画像処理装置は、電気信号を増幅する調整部を有する。調整部は操作部を有する。操作部による指示を受けて、調整部が電気信号を増幅する（ゲイン調整する）ことにより、画像の明るさが調整される。
画像に、軟部組織や筋肉などの低周波成分の構造から骨などの高周波成分の構造までが描画されるときがある。このような画像に対し、点や線状の陰影を強調するための周波数処理が行われる。
【０００５】
ゲイン調整が暗めに働くと、画像において必要な情報が失われるおそれがある。また、画像の明るさをリアルタイムで調整したいという要請がある。また、ユーザーや検査対象に応じて、画像の明るさを簡易に調整したいという要請もある。さらに、周波数処理において、低周波数成分における画像のみの明るさを調整したいという要請もある。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２０１２−１０５９７３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、空間的に２次元でゲイン調整をすると、計算量が膨大なため、ゲイン調整に多くの時間を要し、画像の明るさをリアルタイムで調整することが困難となる。さらに、周波数処理において、低周波成分における画像のみの明るさを調整することは困難を伴う。
【０００８】
この実施形態は、上記の問題を解決するものであり、画像の明るさをリアルタイムで調整することが可能な画像処理装置および超音波診断装置を提供することを目的とする。
【課題を解決するための手段】
【０００９】
上記課題を解決するために、実施形態の画像処理装置は、多重解像度分解部と、調整部と、復元部と、を有する。多重解像度分解部は、画素が２次元または３次元の方向に配列された原画像を高周波成分と低周波成分とに分解することを、画素が各方向に所定個になるまで繰り返すことにより、低周波成分における最下層の画像を生成する。調整部は、最下層の画像における各画素の画素値を変更することにより、最下層の画像の明るさを調整する。復元部は、明るさが調整された最下層の画像を用いて、画素が各方向に前記原画像と同じ個数になるように画像を復元することを特徴とする。
【図面の簡単な説明】
【００１０】
【図１】第１実施形態に係る画像処理装置の構成ブロック図。
【図２】画像処理装置の一連の動作を示すフローチャート。
【図３】離散ウェーブレット変換の一連の動作を示すフローチャート。
【図４】３階層までの離散ウェーブレット変換の一例を示す図。
【図５】１次の高周波成分を示す図。
【図６】２次の高周波成分を示す図。
【図７】３次の高周波成分を示す図。
【図８】最下層の画像における低周波成分を示す図。
【図９】ユーザーおよび検査対象と係数との対応関係を示す図。
【図１０】第２実施形態において、離散ウェーブレット変換から拡大処理までの動作を示すフローチャート。
【図１１】実施形態の種類と動作との対応関係を示す図。
【図１２】第３実施形態において、離散ウェーブレット変換から拡大処理までの動作を示すフローチャート。
【発明を実施するための形態】
【００１１】
（原画像）
モダリティにより検査対象が診断され、診断時における検査対象の原画像ｇが取得される。なお、原画像ｇは、画素が２次元または３次元の方向に配列されたものである。
【００１２】
２次元の方向にｍ個配列された原画像ｇは、次のように表される。
【００１３】
【数１】
2015083056.tif 000003
原画像ｇにおいて、画素の画素値（輝度値または濃度値ともいう）をｐｉｊで表す。ｐの添え字である“ｉ”は行番号（ｉ＝１，２，．．．，ｍ）、“ｊ”は列番号を表す（ｊ＝１，２，．．．，ｍ）。
【００１４】
以下の説明では、原画像ｇを、画素が２次元の方向に６４個（ｍ＝６４）配列されたものとして説明する。
【００１５】
［第１実施形態］
画像処理装置の第１実施形態について各図を参照して説明する。
先ず、画像処理装置の基本的な構成について、図１および図２を参照して簡単に説明する。図１は、画像処理装置１０の構成ブロック図である。
【００１６】
図１に示すように、画像処理装置１０は、多重解像度分解部１１と、調整部１２と、フィルタ１３と、復元部１４と、記憶部１５と、ユーザーインターフェース（ＵＩＦ）１７と、システムゲイン設定部１９と、を有する。
【００１７】
〔多重解像度分解部１１〕
多重解像度解析にはウェーブレット（Ｗａｖｅｌｅｔ）やガウシアンピラミット等が用いられる。ここでは、一例として、ウェーブレットを用いて多重解像度解析を行う多重解像度分解部１１について説明する。
多重解像度分解部１１は、原画像ｇを、離散ウェーブレット変換により高周波成分と低周波成分とに分解する。第１実施形態では、各方向の画素が所定個ｒになるまで、離散ウェーブレット変換が繰り返される。
【００１８】
第１実施形態では、所定個ｒは８である（ｒ＝８）。離散ウェーブレット変換が繰り返される都度、各方向の画素の個数が半分（後述するダウンサンプリング）になるので、各方向の画素の個数が６４個から８個になるためには、離散ウェーブレット変換が階層数Ｎ（＝３）回繰り返される必要がある。
【００１９】
図２は、画像処理装置１０の一連の動作を示すフローチャートである。図２に、画像処理装置１０の一連の動作をステップＳ１０１〜Ｓ１１０で示す。なお、画像処理装置１０の一連の動作は自動的に実行されるものとする。
図２に示すように、この例では、画像処理装置１０の一連の動作において、３階層の離散ウェーブレット変換が行われる（ステップＳ１０１〜Ｓ１０３）。図２に、各階層の離散ウェーブレット変換を、破線により囲まれた部分で示す。
【００２０】
〔１階層目の離散ウェーブレット変換：ステップＳ１０１〕
１階層目の離散ウェーブレット変換について、図３を参照して説明する。図３は、離散ウェーブレットの一連の動作を示すフローチャートである（ステップＳ２０１〜Ｓ２０４）。
【００２１】
（水平方向の離散ウェーブレット変換：ステップＳ２０１）
図３に示すように、最初に、原画像ｇに対して、水平方向（行方向）に離散ウェーブレット変換が行われる（ステップＳ２０１）。
【００２２】
水平方向（行方向）に離散ウェーブレット変換が行われたときの変換後のＦは、次の式で表される。
【００２３】
【数２】
2015083056.tif 000004
ここで、ｇｔは、ｇの転置行列である。また、Ｃは、基底関数（正方行列）である。
【００２４】
基底関数Ｃは、次の式で表される。
【００２５】
【数３】
2015083056.tif 000005
【００２６】
水平方向（行方向）に行われる離散ウェーブレット変換において、原画像ｇが水平方向の低周波成分と高周波成分とに分解される。
【００２７】
（垂直方向の離散ウェーブレット変換：ステップＳ２０２）
次に、変換後のＦに対し、垂直方向（列方向）に離散ウェーブレット変換が行われる（ステップＳ２０２）。
【００２８】
垂直方向（列方向）に離散ウェーブレット変換が行われたときの変換後のＳは、次の式で表される。
【００２９】
【数４】
2015083056.tif 000006
【００３０】
変換後のＳは、各周波数成分（水平方向・垂直方向ともに低周波成分；水平方向が低周波成分でかつ垂直方向が高周波成分；水平方向が高周波成分でかつ垂直方向が低周波成分；水平方向・垂直方向ともに高周波成分）から構成される。
【００３１】
（離散ウェーブレット逆変換：Ｓ２０３）
次に、変換後のＳについて、離散ウェーブレット逆変換が行われる（Ｓ２０３）。
各周波数成分は、離散ウェーブレット逆変換により、画像ｇ（１）に変換される。
【００３２】
（ダウンサンプリング：Ｓ２０４）
次に、画像ｇ（１）が、半分にダウンサンプリングされる（Ｓ２０４）。ダウンサンプリングにより、画素がｍ×ｍ（６４×６４）から３２×３２に配列される。
【００３３】
図４は、３階層まで行われた離散ウェーブレット変換の一例を示す図である。図４に、画像ｇ（１）における、水平方向・垂直方向ともに低周波成分を“ＬＬ１”、水平方向が低周波成分でかつ垂直方向が高周波成分を“ＬＨ１”、水平方向が高周波成分でかつ垂直方向が低周波成分を“ＨＬ１”、水平方向・垂直方向ともに高周波成分を“ＨＨ１”で示す。また、Ｈ、Ｌの添え字は、変換されたときの階層数を示す。
【００３４】
画像ｇ（１）における各周波数成分は、次の式で表される。
【００３５】
【数５】
2015083056.tif 000007
【００３６】
１階層目の離散ウェーブレット変換において、単に「低周波成分」というときは、水平方向・垂直方向ともに低周波成分のことをいい、それ以外の成分を、単に「高周波成分」という（２階層目以降の離散ウェーブレット変換においても同様とする）。
【００３７】
なお、画像ｇ（１）における高周波成分ＬＨ１、ＨＬ１、および、ＨＨ１は、記憶部１５に記憶される。
図５は１次の高周波成分を示す図である。画像ｇ（１）における高周波成分の一例を図５に示す。
【００３８】
〔２階層目の離散ウェーブレット変換：ステップＳ１０２〕
２階層目の離散ウェーブレット変換も１階層目の離散ウェーブレット変換と同様に行われるため、その説明を省略する。
【００３９】
２階層目の離散ウェーブレット変換により、画像ｇ（１）における低周波成分ＬＬ１は、各周波数成分（水平方向・垂直方向ともに低周波成分；水平方向が低周波成分でかつ垂直方向が高周波成分；水平方向が高周波成分でかつ垂直方向が低周波成分；水平方向・垂直方向ともに高周波成分）に分解される。
図４に、画像ｇ（２）における、水平方向・垂直方向ともに低周波成分を“ＬＬ２”、水平方向が低周波成分でかつ垂直方向が高周波成分を“ＬＨ２”、水平方向が高周波成分でかつ垂直方向が低周波成分を“ＨＬ２”、水平方向・垂直方向ともに高周波成分を“ＨＨ２”で示す。
【００４０】
画像ｇ（２）における各周波数成分は、次の式で表される。
【００４１】
【数６】
2015083056.tif 000008
【００４２】
２階層目の離散ウェーブレット変換において、ダウンサンプリングにより、画素が３２×３２から１６×１６に配列される。
なお、高周波成分ＬＨ２、ＨＬ２、および、ＨＨ２は、記憶部１５に記憶される。
図６は２次の高周波成分を示す図である。画像ｇ（２）における高周波成分の一例を図６に示す。
【００４３】
〔３階層目の離散ウェーブレット変換：ステップＳ１０３〕
３階層目の離散ウェーブレット変換も１階層目の離散ウェーブレット変換と同様に行われるため、その説明を省略する。
【００４４】
３階層目の離散ウェーブレット変換により、画像ｇ（２）における低周波成分ＬＬ２は、各周波数成分（水平方向・垂直方向ともに低周波成分；水平方向が低周波成分でかつ垂直方向が高周波成分；水平方向が高周波成分でかつ垂直方向が低周波成分；水平方向・垂直方向ともに高周波成分）に分解される。
【００４５】
図４に、画像ｇ（３）における、水平方向・垂直方向ともに低周波成分を“ＬＬ３”、水平方向が低周波成分でかつ垂直方向が高周波成分を“ＬＨ３”、水平方向が高周波成分でかつ垂直方向が低周波成分を“ＨＬ３”、水平方向・垂直方向ともに高周波成分を“ＨＨ３”で示す。
【００４６】
画像ｇ（３）における各周波数成分は、次の式で表される。
【００４７】
【数７】
2015083056.tif 000009
【００４８】
３階層目の離散ウェーブレット変換において、ダウンサンプリングにより、画像値が１６×１６から８×８の配列となる。すなわち、各方向に配列された画素の個数が所定個である８個になる（ｒ＝８）。それにより、離散ウェーブレット変換の繰り返しが３階層目で終了する。画像ｇ（３）が最下層の画像となる。この画像ｇ（３）における低周波成分ＬＬ３は、次の調整部１２による調整対象となる。ここで、「最下層の画像」とは、離散ウェーブレット変換が予め定められた階層数繰り返されたとき、最後の繰り返しが終了したときの画像をいう。
【００４９】
図７は３次の高周波成分を示す図である。画像ｇ（３）における高周波成分の一例を図７に示す。また、図８は、最下層の画像における低周波成分を示す図である。画像ｇ（３）における低周波成分の一例を図８に示す。
【００５０】
〔調整部１２による明るさ調整１〕
最下層の画像ｇ（３）における低周波成分をＬＬ３で示した。ここで、最下層の画像ｇ（３）における画像の画素値も“ＬＬ３”で表すものとする。最下層の画像ｇ（３）における画像の画素値ＬＬ３は、８×８に配列される。
調整部１２は、画像ｇ（３）における画像の画素値ＬＬ３に対し係数αを乗算することにより、画素値ＬＬ３を変更する（図２に示すステップＳ１０４；明るさ調整１）。変更された画素値ＬＬ３を所定の値ｖ０とすると、所定の値ｖ０は、次の式で表される。
【００５１】
【数８】
2015083056.tif 000010
所定の値ｖ０も、画素が８×８に配列されたものとなる。
【００５２】
〔フィルタ１３〕
フィルタ１３は、ぼかし処理をすることにより、所定の値ｖ０（８×８に配列された画素の画素値）の明るさを均一にする（図２に示すステップＳ１０５；明るさの均一化）。ここで、「明るさを均一にする」あるいは「明るさの均一化」とは、明るさの変化を少なくすることを意味する。
ぼかし処理においては、８×８に配列された画素の画素値に対して所定の値ｖ０を求めればよいため、計算量が膨大にならない。
【００５３】
ぼかし処理の一例としては、平滑化フィルタにおける係数の重みをいう。平滑化フィルタとして、例えば、移動平均フィルタ（平均化フィルタとも呼ばれる）、ガウシアンフィルタ、メディアンフィルタが用いられる。
【００５４】
移動平均フィルタは、対象となる画素（８×８に配列された画素）の画素値を計算するために、それに隣接した画素を用いるものである。
フィルタサイズ３×３の移動平均フィルタを次に示す。
2015083056.tif 000011
ここで、係数の重みを“ｗ”で示し、ｗの添え字は、行番号および列番号を示す。
【００５５】
フィルタサイズ３×３の移動平均フィルタの一例としては、対象となる画素の画素値を計算するために、それに隣接した８つの画素の画素値を用いるものであって、重みｗを合計した値は１とし、全ての重みｗを１／９とする。
【００５６】
明るさの均一化において、対象となる所定の値ｖ０（８×８に配列された画素の画素値）を次に示す。
2015083056.tif 000012
所定の値ｖ０において、８×８に配列された画素の画素値をｑｉｊで表す。ｑの添え字である“ｉ”は行番号（ｉ＝１，２，．．．，８）、“ｊ”は、列番号を表す（ｊ＝１，２，．．．，８）。
【００５７】
明るさの均一化には、重みｗ＝１／９であって、フィルタサイズが３×３である移動平均フィルタが用いられる。このとき、例えば、対象となる所定の値ｖ０（８×８に配列された画素の画素値）のうち、行番号３で列番号３の位置に配置された画素の画素値ｑ３３は、次のように表される。
2015083056.tif 000013
ここで、ｐ´３３は、画素値ｑ３３がフィルタ１３のぼかし処理により、変更された値を示す。
【００５８】
以上のようなフィルタ１３のぼかし処理により、所定の値ｖ０がｖ１に変更され、最下層の画像の明るさが均一になる。
変換された値ｖ１は次の式で表される。
【００５９】
【数９】
2015083056.tif 000014
βは、フィルタサイズおよび重みを示す数値の並びに相当するものである。画素の画素値が８×８に配列された所定の値ｖ０に対し値ｖ１を求めればよいので、計算量が膨大にならずに済む。
【００６０】
なお、移動平均フィルタは、３×３のフィルタサイズのものに限らない。例えば、ｎ(≧４）×ｎのフィルタサイズのものを用いることで、明るさをさらに均一にするものであってもよい。さらに、フィルタ１３は、移動平均フィルタに限らず、他の平滑化フィルタであってもよい。
【００６１】
なお、次のように構成してもよい。これらのフィルタ１３を、ユーザーおよび検査対象に対応づけて予め記憶部１５に記憶する。画像処理装置１０は、ユーザーおよび検査対象の情報を受けて、ユーザー等に対応するフィルタを用いることにより、ぼかし処理するようにしてもよい。ぼかし処理は、ＵＩＦ１７による操作を受けて、フィルタ調整部１８により行われる。
【００６２】
〔調整部１２による明るさ調整２〕
調整部１２は、フィルタ１３により変換された値ｖ１（８×８に配列された画素の画素値）にオフセット値や補正係数を乗算／加算することにより、ｖ２に調整する（図２に示すステップＳ１０６；明るさ調整２）。調整された値ｖ２は、次の式で表される。
【００６３】
【数１０】
2015083056.tif 000015
ここで、ｋ０は補正係数の初期値、ａ０はオフセット値の初期値である。
【００６４】
初期値ｋ０、ａ０は、ユーザーおよび検査対象に対応づけられて予め記憶部１５に記憶される。
【００６５】
図９は、ユーザーおよび検査対象と係数との対応関係を示す図である。図９に、初期値ｋ０、ａ０に乗算すべき係数を、ユーザーおよび検査対象に対応づけて示す。初期値ｋ０、ａ０、に乗算すべき係数が記憶部１５に記憶される。調整部１２は、記憶部１５から読み出した初期値ｋ０、ａ０および、記憶部１５から読み出した、初期値に乗算すべき係数に基づき、ユーザーおよび検査対象に対応づけられた係数を参照して、値ｖ１を値ｖ２に調整する。
【００６６】
このとき、値ｖ２は次の式で表される。
【００６７】
【数１１】
2015083056.tif 000016
ここで、γ、εは、初期値ｋ０、ａ０に乗算される、ユーザーに対応づけられた係数である。また、δ、ηは、初期値ｋ０、ａ０に乗算される、検査対象に対応づけられた係数である。
【００６８】
図９に示すように、初期値ｋ０、ａ０、に乗算すべき係数は、ユーザーの識別番号（ＵＩＤ）が“ＡＡＡＡＡＡＡＡ”のとき、“１．０”、“１．０”である。さらに、ＵＩＤが“ＢＢＢＢＢＢＢＢ”のとき、係数は“１．１”、“１．４”である。さらに、検査対象が“ＣＣＣ”のとき、係数は“１．０”、“１．０”である。さらに、検査対象が“ＤＤＤ”のとき、係数は“１．２”、“１．５”である。
【００６９】
例えば、ＵＩＤが“ＢＢＢＢＢＢＢＢ”、かつ、検査対象が“ＤＤＤ”のとき、調整された値ｖ２は、次の式で表される。
【００７０】
【数１２】
2015083056.tif 000017
以上のように、調整部１２は、式（１１）に基づき、記憶部１５に記憶された初期値ｋ０、ａ０およびそれらに乗算すべき係数を参照して、フィルタ１３により変換された値ｖ１（８×８に配列された画素の画素値）をユーザーおよび検査対象に対応づけられた値ｖ２に調整する。
【００７１】
さらに、調整部１２は、値ｖ２を調整可能に構成される。値ｖ２の調整は、調整された値ｖ２に対し直接的に行ってもよく、初期値ｋ０、ａ０やそれらに乗算すべき係数を変更することにより間接的に行ってもよい。値ｖ２、初期値ｋ０、ａ０、係数の変更はユーザーインターフェース（ＵＩＦ）１７による操作を受けて、数値調整部１６により行われる。
【００７２】
なお、第１実施形態では、明るさの均一化（ステップＳ１０５）の後に、明るさ調整２（ステップＳ１０６）をしたが、動作の順番を逆にしてもよい。
【００７３】
〔復元部１４〕
復元部１４は、明るさが調整された、各方向に８個の画素が配列された低周波成分における最下層の画像（８×８の画像）を、各方向に６４個の画素が配列された原画像と同じサイズ（６４×６４の画像）にアップサンプリングする（図２に示すステップＳ１０７；拡大処理）。アップサンプリングされた画像をｇ（４）とする。このアップサンプリングでは、水平、垂直の２次元方向に画素を８倍すればよいため、短時間に計算することが可能となる。ここで、「拡大処理」とは、対象の画像における各方向の画素の個数を原画像と同じ個数に戻すことをいう。以下、ステップＳ１０８、Ｓ１０９で行われる拡大処理において同様である。
【００７４】
画像の復元には、高周波成分の画像が用いられる。ここで、高周波成分の画像とは、前述したように、水平方向・垂直方向ともに低周波成分（ＬＬ１、および、ＬＬ２）以外の周波数成分（ＬＨ１、ＨＬ１、および、ＨＨ１、並びに、ＬＨ２、ＨＬ２、および、ＨＨ２）の画像をいう。
【００７５】
復元部１４は、記憶部１５に記憶された、高周波成分ＬＨ２、ＨＬ２、および、ＨＨ２の画像（１６×１６の画像）をそれぞれ原画像と同じサイズ（６４×６４の画像）にアップサンプリングする（図２に示すステップＳ１０８；拡大処理）。アップサンプリングされた画像をｇ（５）とする。このアップサンプリングでは、水平、垂直の２次元方向に画素の個数を４倍すればよいため、短時間に計算することが可能となる。
【００７６】
さらに、復元部１４は、記憶部１５に記憶された、高周波成分ＬＨ１、ＨＬ１、および、ＨＨ１の画像（３２×３２の画像）をそれぞれ原画像と同じサイズ（６４×６４の画像）にアップサンプリングする（図２に示すステップＳ１０９；拡大処理）。アップサンプリングされた画像をｇ（６）とする。このアップサンプリングでは、水平、垂直の２次元方向に画素の個数を２倍すればよいため、短時間に計算することが可能となる。
【００７７】
復元部１４は、アップサンプリングされた画像ｇ（４）、ｇ（５）、ｇ（６）を加算することにより、画像ｇ（ｎ）を復元する（図２に示すステップＳ１１０；加算処理）。
【００７８】
復元された画像ｇ（ｎ）は、次の式で表される。
【００７９】
【数１３】
2015083056.tif 000018
画像ｇ（ｎ）を復元するためには、６４×６４の画像である、画像ｇ（４）、ｇ（５）、ｇ（６）を加算すればよいので、算出に多くの時間を要しない。
なお、第１実施形態では、復元部１４により画像ｇ（４）に画像ｇ（５）、ｇ（６）の両方が加算されたものを示したが、いずれか一方が加算されてもよい。
【００８０】
さらに、第１実施形態では、フィルタ１３による明るさの均一化をし（ステップＳ１０５）、その後、明るさ調整２をした（ステップＳ１０６）。しかし、明るさ調整２をし、その後、明るさの均一化をしてもよい。
【００８１】
前記第１実施形態によれば、骨などの高周波成分の構造や軟部組織などの低周波成分の構造が混在していても、構造はそのままで、骨の明るさと軟部組織の明るさを近いレベルに自動でかつリアルタイムで調整でき、ゲイン調整の手間を省くことができる。さらに、心臓の検査においては、心筋部以外のノイズ成分（高周波成分）を暗く表示することも可能となる。
【００８２】
［第２実施形態］
次に、画像処理装置１０の第２実施形態について図１０および図１１を参照して説明する。図１０は、第２実施形態において、離散ウェーブレット変換から拡大処理までの動作を示すフローチャート、図１１は、実施形態の種類と動作との対応関係を示す図である。なお、第２実施形態において、第１実施形態と同じ構成については同一番号を付してその説明を省略し、異なる構成について主に説明する。
【００８３】
第１実施形態では、図２に示すように、離散ウェーブレット変換による低周波成分ＬＬ３の抽出（ステップＳ１０３）、その後、低周波成分ＬＬ３を所定の値ｖ０とする明るさ調整１（ステップＳ１０４）、その後、フィルタ１３による明るさの均一化（ステップＳ１０５）、その後、明るさ調整２（ステップＳ１０６）、その後、拡大処理（ステップＳ１０７）を行うように構成した。
【００８４】
第２実施形態では、図１１に示すように、明るさ調整１（ステップＳ１０４）、明るさの均一化（ステップＳ１０５）、および、明るさ調整２（ステップＳ１０６）の３つの処理のうち、一つの処理を省略する。それにより、画像処理の迅速化を図ることができる。
【００８５】
一つの形態は、図１０に示すように、明るさ調整２（ステップＳ１０６）を省略する。すなわち、明るさ調整１（ステップＳ１０４）の後に、明るさの均一化（ステップＳ１０５）をするように構成する。
【００８６】
他の形態は、明るさの均一化（ステップＳ１０５）を省略する。すなわち、明るさ調整１（ステップＳ１０４）の後に、明るさ調整２（ステップＳ１０６）をするように構成してもよい。
【００８７】
さらに、他の形態は、明るさ調整１（ステップＳ１０４）を省略する。このとき、処理の順番としては、明るさの均一化（ステップＳ１０５）の後に、明るさ調整２（ステップＳ１０６）をするように構成してもよく、明るさ調整２（ステップＳ１０６）の後に、明るさの均一化（ステップＳ１０５）をするように構成してもよい。
【００８８】
これらの形態における各処理が終了した後に、拡大処理（ステップＳ１０７、Ｓ１０８、Ｓ１０９）が行われる。
【００８９】
前記第２実施形態によれば、高周波成分の構造や低周波成分の構造が混在していても、構造はそのままで、高周波成分の構造と低周波成分の構造の明るさを近いレベルに自動でかつリアルタイムで調整できる。
【００９０】
［第３実施形態］
次に、画像処理装置１０の第３実施形態について図１１および図１２を参照して説明する。図１２は、第３実施形態において、離散ウェーブレット変換から拡大処理までの動作を示すフローチャートである。なお、第３実施形態において、第１実施形態と同じ構成については同一番号を付してその説明を省略し、異なる構成について主に説明する。
【００９１】
第３実施形態では、図１１に示すように、明るさ調整１（ステップＳ１０４）、明るさの均一化（ステップＳ１０５）、および、明るさ調整２（ステップＳ１０６）の３つの処理のうち、二つの処理を省略する。それにより、画像処理のさらなる迅速化を図ることができる。
【００９２】
一つの形態は、図１２に示すように、明るさの均一化（ステップＳ１０５）、および、明るさ調整２（ステップＳ１０６）を省略する。すなわち、明るさ調整１（ステップＳ１０４）のみを行うように構成してもよい。
【００９３】
他の形態は、図１１に示すように、明るさ調整１（ステップＳ１０４）、および、明るさ調整２（ステップＳ１０６）を省略する。すなわち、明るさの均一化（ステップＳ１０５）のみを行うように構成してもよい。この明るさの均一化（ステップＳ１０５）の後に、ある任意の係数を乗算して、輝度変化を強調または弱めてもよい。輝度変化を強調すれば、輝度の低いノイズを低減することができる。
【００９４】
さらに、他の形態は、図１１に示すように、明るさ調整１（ステップＳ１０４）、および、明るさの均一化（ステップＳ１０５）を省略する。すなわち、明るさ調整２（ステップＳ１０６）のみを行うように構成してもよい。
【００９５】
これらの形態における各処理が終了した後に、拡大処理（ステップＳ１０７、Ｓ１０８、Ｓ１０９）が行われる。
【００９６】
前記第３実施形態によれば、高周波成分の構造や低周波成分の構造が混在していても、構造はそのままで、高周波成分の構造と低周波成分の構造の明るさを近いレベルに自動でかつリアルタイムで調整できる。
【００９７】
［第４実施形態］
次に、画像処理装置１０の第４実施形態について図１を参照して説明する。なお、第４実施形態において、第１実施形態と同じ構成については同一番号を付してその説明を省略し、異なる構成について主に説明する。
【００９８】
第４実施形態では、図１に示すように、フィルタ１３のぼかし処理により計算された値ｖ１（８×８に配列された画素の画素値）が、システムゲイン設定部１９に出力されるように構成する。システムゲイン設定部１９は、値ｖ１に基づいて、ゲイン調整することにより、画像処理装置１０の性能を最適化することができる。
システムゲイン設定部１９によるゲイン調整は、値ｖ１（８×８に配列された画素の画素値）に基づいて行われるので、計算量が膨大にならないで済む。
【００９９】
［超音波診断装置］
次に、前記各実施形態に係る画像処理装置１０を備えた超音波診断装置について、説明する。
【０１００】
前述したように、超音波診断装置は、超音波パルスを体内に照射する電気音響変換器および画像処理装置１０を有する。生体から反射して戻ってきた反射信号は電気音響変換器によって電気信号に変換される。電気信号は増幅器で増幅される。増幅された信号に基づいて、画像が構築され、画像がディスプレイに出力される。構築される画像としては、Ｂモード画像、Ｍモード画像、パルスドプラ方式により取得された画像、連続波ドプラ方式により取得された画像がある。
【０１０１】
Ｂモード画像は、反射信号が発生した時間と反射信号の強度から、生体内部組織の音響的反射強度分布を断層画像としたものである。Ｍモード画像は、縦軸には反射強度を輝度に変換したものを、横軸には縦軸で輝度表示したものを、時系列で表示したものである。パルスドプラ方式により取得された画像は、超音波パルスを繰り返し照射し、反射信号の周波数偏移から移動組織（血流）の速度を求め、電気音響変換器（振動子）に向かってくる速度の方向は赤色成分で、反対に、遠ざかる方向は青色成分で表示されるものである。連続波ドプラ方式により取得された画像は、音波パルスを連続的に照射し、反射信号を周波数偏移してから移動組織（血流）の速度を求め、電気音響変換器（振動子）に向かってくる速度の方向は赤色成分で、反対に、遠ざかる方向は青色成分で表示されるものである。
【０１０２】
超音波診断装置は、前述する画像処理装置１０を備えている。超音波診断装置は、Ｂモード画像；Ｍモード画像；パルスドプラ方式により取得された画像；連続波ドプラ方式により取得された画像の少なくとも一つを原画像ｇとして画像処理装置１０に出力する。
【０１０３】
さらに、超音波診断装置によりフリーズ前に取得された画像、および／または、超音波診断装置により過去に取得された画像が、原画像ｇとして画像処理装置１０に出力されるように構成される。
【０１０４】
ここで、「フリーズ」とは、Ｂモード画像等が表示されているとき、例えば、ユーザーの操作を受けて、画像が静止されることをいう。フリーズされたときの画像は、自動的に、または、ユーザーの操作により、画像処理装置１０の記憶部１５に記憶される。
【０１０５】
なお、前記各実施形態では、１階層目の離散ウェーブレット変換により、原画像ｇを低周波成分ＬＬ１と高周波成分ＬＨ１、ＨＬ１、および、ＨＨ１とに分解し、高周波成分ＬＨ１、ＨＬ１、および、ＨＨ１を、画像ｇ（６）として画像の復元に用いたが、これに限らない。高周波成分を変更せず、そのまま画像の復元に用いるように構成してもよい。それにより、診断に必要となる検査対象に関する情報を失うことがない。
【０１０６】
また、前記各実施形態では、明るさ調整１（ステップＳ１０４）、明るさの均一化（ステップＳ１０５）、および、明るさ調整２（ステップＳ１０６）の３つの処理のうち、少なくとも一つの処理を行うように構成したが、全部の処理を省略してもよい。それにより、画像処理の迅速化をさらに図ることができる。
【０１０７】
さらに、前記各実施形態では、離散ウェーブレット変換が２次元データとしての原画像ｇに対し行われるものを示したが、３次元データとしての原画像に対し行われるように構成されてもよく、このように構成されても、２次元データのときと同様に、原画像の明るさ調整をリアルタイムに行うことができる。
【０１０８】
なお、前記各実施形態では、画素が６４×６４に配列された原画像ｇに対し、離散ウェーブレット変換を３階層繰り返すことにより、画素が８×８に配列された最下層の画像を生成したが、画素が４×４に配列された最下層の画像を生成するには、離散ウェーブレット変換を４階層繰り返す必要があることはいうまでもない。さらに、画素の個数がｍ（≧１２８）×ｍに配列された原画像ｇに対し、画素の個数がｒ（≦８）×ｒに配列された最下層の画像を生成するには、離散ウェーブレット変換を４階層以上繰り返す必要があることはいうまでもない。
【０１０９】
また、前記各実施形態では、離散ウェーブレット変換の対象としての原画像ｇを画素の個数が２次元または３次元の方向にｍ個配列されたデータとして説明したが、座標変換前のデータ形式である生データ（Ｒａｗ  ｄａｔａ）であってもよい。
【０１１０】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるととともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
【符号の説明】
【０１１１】
１０  画像処理装置
１１  多重解像度分解部
１２  調整部
１３  フィルタ
１４  復元部
１５  記憶部
１６  数値調整部
１７  ユーザーインターフェース（ＵＩＦ）
１８  フィルタ調整部
１９  システムゲイン設定部
【図１】
2015083056.tif 000019
【図２】
2015083056.tif 000020
【図３】
2015083056.tif 000021
【図４】
2015083056.tif 000022
【図５】
2015083056.tif 000023
【図６】
2015083056.tif 000024
【図７】
2015083056.tif 000025
【図８】
2015083056.tif 000026
【図９】
2015083056.tif 000027
【図１０】
2015083056.tif 000028
【図１１】
2015083056.tif 000029
【図１２】
2015083056.tif 000030