(57)【要約】【課題】点広がり関数を、少ないリソースの使用で迅速かつ適切に決定するようにすることによって、高解像度の超音波画像を生成できるようにする画像処理モジュールと画像生成方法を提供する。【解決手段】画像処理モジュールは、対象体の目標部位に超音波を照射し、対象体の目標部位で反射されるエコー超音波を用いて画像を取得する画像処理モジュールであって、入力信号をビームフォーミングして、ビームフォーミングされた信号を出力するビームフォーミング部と、ビームフォーミングされた信号に対して少なくとも一つの状況変数に基づいて予め取得された少なくとも一つの２次元の点広がり関数で構築された点広がり関数データベースと、点広がり関数データベースから少なくとも一つの２次元の点広がり関数を選択し、ビームフォーミングされた信号及び選択された点広がり関数を用いてデコンボリューションを行って画像を生成する画像生成部とを有する。【選択図】      図１２
2014023934.tif 000002
【特許請求の範囲】
【請求項１】
対象体の目標部位に超音波を照射し、前記対象体の目標部位で反射されるエコー超音波を用いて画像を取得する画像処理モジュールであって、
入力信号をビームフォーミングして、ビームフォーミングされた信号を出力するビームフォーミング部と、
前記ビームフォーミングされた信号に対して少なくとも一つの状況変数に基づいて予め取得された少なくとも一つの２次元の点広がり関数（ＰＳＦ、ｐｏｉｎｔ  ｓｐｒｅａｄ  ｆｕｎｃｔｉｏｎ）で構築された点広がり関数データベースと、
前記点広がり関数データベースから少なくとも一つの２次元の点広がり関数を選択し、前記ビームフォーミングされた信号及び前記選択された点広がり関数を用いてデコンボリューション（ｄｅｃｏｎｖｏｌｕｔｉｏｎ）を行って画像を生成する画像生成部とを有することを特徴とする画像処理モジュール。
【請求項２】
前記画像生成部は、前記ビームフォーミングされた信号を前記少なくとも一つの状況変数によって分類し、前記分類した結果に基づいて、前記点広がり関数データベースから少なくとも一つの点広がり関数を抽出することを特徴とする請求項１に記載の画像処理モジュール。
【請求項３】
対象体の目標部位に超音波を照射し、前記対象体の目標部位で反射されるエコー超音波を用いて画像を取得する画像処理モジュールであって、
入力信号をビームフォーミングして、ビームフォーミングされた信号を出力するビームフォーミング部と、
前記ビームフォーミングされた信号に基づいて少なくとも一つの第１点広がり関数を推定し、前記推定された少なくとも一つの第１点広がり関数を用いて第２点広がり関数を推定した後、前記ビームフォーミングされた信号及び前記推定された第２点広がり関数を用いてデコンボリューションを行って画像を生成する画像生成部とを有することを特徴とする画像処理モジュール。
【請求項４】
前記ビームフォーミングされた信号に対する少なくとも一つの点広がり関数で構築された点広がり関数データベースをさらに含み、
前記少なくとも一つの点広がり関数は、少なくとも一つの状況変数に基づいて予め取得された少なくとも一つの点広がり関数であることを特徴とする請求項３に記載の画像処理モジュール。
【請求項５】
前記少なくとも一つの状況変数は、前記超音波の音速、又は前記対象体内部の目標部位との距離であることを特徴とする請求項１又は４に記載の画像処理モジュール。
【請求項６】
前記画像生成部は、前記点広がり関数データベースから少なくとも一つの点広がり関数を抽出して、少なくとも一つの第１点広がり関数を推定することを特徴とする請求項４に記載の画像処理モジュール。
【請求項７】
前記画像生成部は、前記ビームフォーミングされた信号を前記少なくとも一つの状況変数によって分類し、前記分類した結果に基づいて、前記点広がり関数データベースから少なくとも一つの点広がり関数を抽出して、少なくとも一つの第１点広がり関数を推定することを特徴とする請求項６に記載の画像処理モジュール。
【請求項８】
前記画像生成部は、前記ビームフォーミングされた信号及び前記推定された少なくとも一つの第１点広がり関数を用いて第２点広がり関数を推定することを特徴とする請求項３に記載の画像処理モジュール。
【請求項９】
前記画像生成部は、少なくとも一つのフィルターを用いて第１点広がり関数を推定することを特徴とする請求項３に記載の画像処理モジュール。
【請求項１０】
前記少なくとも一つのフィルターは、最小二乗フィルター（ＬＳＦ、ｌｅａｓｔ  ｓｑｕａｒｅ  ｆｉｌｔｅｒ）又はケプストラムフィルター（ｃｅｐｓｔｒｕｍ  ｆｉｌｔｅｒ）であることを特徴とする請求項９に記載の画像処理モジュール。
【請求項１１】
前記画像生成部は、前記第１点広がり関数と少なくとも一つの第２点広がり関数推定変数とを用いて第２点広がり関数を推定することを特徴とする請求項９に記載の画像処理モジュール。
【請求項１２】
前記少なくとも一つの第２点広がり関数推定変数は、ビームフォーミング係数であることを特徴とする請求項１１に記載の画像処理モジュール。
【請求項１３】
前記少なくとも一つの第２点広がり関数推定変数は、前記超音波の音速、前記音速の変化量、前記対象体内部の目標部位との距離、入力信号生成モジュールの配置形態、及び入力信号のチャネル別減衰率の内の少なくとも一つであることを特徴とする請求項１１に記載の画像処理モジュール。
【請求項１４】
前記第１点広がり関数は、１次元の点広がり関数であり、前記第２点広がり関数は２次元の点広がり関数であることを特徴とする請求項３に記載の画像処理モジュール。
【請求項１５】
対象体の目標部位に超音波を照射した後、前記対象体の目標部位で反射されるエコー超音波を受信し、前記受信されたエコー超音波を超音波信号に変換するステップと、
前記超音波信号に基づいてビームフォーミングされた超音波信号を出力するステップと、
前記ビームフォーミングされた超音波信号に基づいて、点広がり関数データベースから少なくとも一つの点広がり関数を抽出するステップと、
前記少なくとも一つの点広がり関数を用いてデコンボリューションを行って、前記ビームフォーミンされた超音波信号から画像を生成するステップとを有し、
前記点広がり関数データベースは、前記出力されたビームフォーミングされた超音波信号に対して少なくとも一つの状況変数に基づいて予め取得された少なくとも一つの点広がり関数で構築されることを特徴とする画像生成方法。
【請求項１６】
前記ビームフォーミングされた超音波信号に基づいて、点広がり関数データベースから少なくとも一つの点広がり関数を抽出するステップは、前記ビームフォーミングされた超音波信号を前記少なくとも一つの状況変数によって分類し、前記分類した結果に基づいて、前記点広がり関数データベースから少なくとも一つの点広がり関数を抽出するステップを含むことを特徴とする請求項１５に記載の画像生成方法。
【請求項１７】
前記少なくとも一つの点広がり関数を用いてデコンボリューションを行って、前記ビームフォーミングされた超音波信号から画像を生成するステップは、前記抽出された少なくとも一つの点広がり関数を用いて他の点広がり関数を推定するステップと、
前記推定された他の点広がり関数を用いてデコンボリューションを行って、前記ビームフォーミングされた信号から画像を生成するステップとを含むことを特徴とする請求項１５に記載の画像生成方法。
【請求項１８】
前記少なくとも一つの状況変数は、前記超音波の音速、又は前記対象体内部の目標部位との距離であることを特徴とする請求項１５又は１７に記載の画像生成方法。
【発明の詳細な説明】
【技術分野】
【０００１】
画像処理モジュールと画像生成方法に関し、特に、複雑な点広がり関数の演算を簡単にする画像処理モジュールと画像生成方法に関する。
【背景技術】
【０００２】
超音波画像装置は、対象体の外面、例えば、人体の表皮から内部の所定の目標地点に向かって超音波を照射し、内部の目標地点で反射された超音波を収集した後、収集された超音波情報を用いて対象体内部の各種組織や構造などに対する断層画像、例えば、各種臓器、軟部組織などの断層画像を取得する装置である。
このような超音波画像装置は、低廉で且つ他の装置に比べて小型であり、Ｘ線のような放射線被曝の危険性がない上に、リアルタイムで対象体の内部に対する画像の再生が可能であるので、医療業界で広く用いられている。
【０００３】
超音波画像装置が超音波画像を得る方法では、超音波画像装置の超音波プローブの超音波発生手段、例えば、トランスデューサ（ｔｒａｎｓｄｕｃｅｒ）を用いて電気的信号を超音波に変換した後、変換された超音波を目標地点に向かって照射する。
そして、目標地点で反射された超音波を、超音波センサ、例えば、トランスデューサで受信し、これを電気的信号に変換して、複数のチャネルの超音波信号を取得する。
次いで、超音波信号に基づいてビームフォーミング（ｂｅａｍ  ｆｏｒｍｉｎｇ）が行われる。ビームフォーミングは、超音波センサが収集した超音波信号の時間差を補正し、各超音波信号毎に所定の加重値を付加して特定位置の信号を強調し、他の位置の信号は相対的に減殺させて超音波信号を集束するようにする。
その後、超音波画像装置は、ビームフォーミングされた超音波信号に基づいて超音波画像を取得する。
【０００４】
この時、超音波信号に基づいた画像生成において、適切な点広がり関数を迅速に決定しなければならないという問題がある。
【発明の概要】
【発明が解決しようとする課題】
【０００５】
本発明の目的は、ビームフォーミングの後に、超音波画像を生成する過程において、複雑な点広がり関数の演算を簡単にして、点広がり関数を、少ないリソースの使用で迅速に決定して、画像生成過程をより迅速に行うことができるようにし、さらに、迅速に画像を生成しながらも点広がり関数を適切に決定するようにすることによって、高解像度の超音波画像を生成できるようにする画像処理モジュールと画像生成方法を提供することにある。
【課題を解決するための手段】
【０００６】
上記目的を達成するためになされた本発明による画像処理モジュールは、対象体の目標部位に超音波を照射し、前記対象体の目標部位で反射されるエコー超音波を用いて画像を取得する画像処理モジュールであって、入力信号をビームフォーミングして、ビームフォーミングされた信号を出力するビームフォーミング部と、前記ビームフォーミングされた信号に対して少なくとも一つの状況変数に基づいて予め取得された少なくとも一つの２次元の点広がり関数（ＰＳＦ、ｐｏｉｎｔ  ｓｐｒｅａｄ  ｆｕｎｃｔｉｏｎ）で構築された点広がり関数データベースと、前記点広がり関数データベースから少なくとも一つの２次元の点広がり関数を選択し、前記ビームフォーミングされた信号及び前記選択された点広がり関数を用いてデコンボリューション（ｄｅｃｏｎｖｏｌｕｔｉｏｎ）を行って画像を生成する画像生成部とを有することを特徴とする。
【０００７】
前記少なくとも一つの状況変数は、前記超音波の音速、又は前記対象体内部の目標部位との距離であることが好ましい。
記画像生成部は、前記ビームフォーミングされた信号を前記少なくとも一つの状況変数によって分類し、前記分類した結果に基づいて、前記点広がり関数データベースから少なくとも一つの点広がり関数を抽出することが好ましい。
【０００８】
また、上記目的を達成するためになされた本発明による画像処理モジュールは、対象体の目標部位に超音波を照射し、前記対象体の目標部位で反射されるエコー超音波を用いて画像を取得する画像処理モジュールであって、入力信号をビームフォーミングして、ビームフォーミングされた信号を出力するビームフォーミング部と、前記ビームフォーミングされた信号に基づいて少なくとも一つの第１点広がり関数を推定し、前記推定された少なくとも一つの第１点広がり関数を用いて第２点広がり関数を推定した後、前記ビームフォーミングされた信号及び前記推定された第２点広がり関数を用いてデコンボリューションを行って画像を生成する画像生成部とを有することを特徴とする。
【０００９】
前記ビームフォーミングされた信号に対する少なくとも一つの点広がり関数で構築された点広がり関数データベースをさらに含み、前記少なくとも一つの点広がり関数は、少なくとも一つの状況変数に基づいて予め取得された少なくとも一つの点広がり関数であることが好ましい。
前記少なくとも一つの状況変数は、前記超音波の音速、又は前記対象体内部の目標部位との距離であることが好ましい。
前記画像生成部は、前記点広がり関数データベースから少なくとも一つの点広がり関数を抽出して、少なくとも一つの第１点広がり関数を推定することが好ましい。
前記画像生成部は、前記ビームフォーミングされた信号を前記少なくとも一つの状況変数によって分類し、前記分類した結果に基づいて、前記点広がり関数データベースから少なくとも一つの点広がり関数を抽出して、少なくとも一つの第１点広がり関数を推定することが好ましい。
前記画像生成部は、前記ビームフォーミングされた信号及び前記推定された少なくとも一つの第１点広がり関数を用いて第２点広がり関数を推定することが好ましい。
【００１０】
前記画像生成部は、少なくとも一つのフィルターを用いて第１点広がり関数を推定することが好ましい。
前記少なくとも一つのフィルターは、最小二乗フィルター（ＬＳＦ、ｌｅａｓｔ  ｓｑｕａｒｅ  ｆｉｌｔｅｒ）又はケプストラムフィルター（ｃｅｐｓｔｒｕｍ  ｆｉｌｔｅｒ）であることが好ましい。
前記画像生成部は、前記第１点広がり関数と少なくとも一つの第２点広がり関数推定変数とを用いて第２点広がり関数を推定することが好ましい。
前記少なくとも一つの第２点広がり関数推定変数は、ビームフォーミング係数であることが好ましい。
前記少なくとも一つの第２点広がり関数推定変数は、前記超音波の音速、前記音速の変化量、前記対象体内部の目標部位との距離、入力信号生成モジュールの配置形態、及び入力信号のチャネル別減衰率の内の少なくとも一つであることが好ましい。
前記第１点広がり関数は、１次元の点広がり関数であり、前記第２点広がり関数は２次元の点広がり関数であることが好ましい。
【００１１】
上記目的を達成するためになされた本発明による画像生成方法は、対象体の目標部位に超音波を照射した後、前記対象体の目標部位で反射されるエコー超音波を受信し、前記受信されたエコー超音波を超音波信号に変換するステップと、前記超音波信号に基づいてビームフォーミングされた超音波信号を出力するステップと、前記ビームフォーミングされた超音波信号に基づいて、点広がり関数データベースから少なくとも一つの点広がり関数を抽出するステップと、前記少なくとも一つの点広がり関数を用いてデコンボリューションを行って、前記ビームフォーミンされた超音波信号から画像を生成するステップとを有し、前記点広がり関数データベースは、前記出力されたビームフォーミングされた超音波信号に対して少なくとも一つの状況変数に基づいて予め取得された少なくとも一つの点広がり関数で構築されることを特徴とする。
【００１２】
前記ビームフォーミングされた超音波信号に基づいて、点広がり関数データベースから少なくとも一つの点広がり関数を抽出するステップは、前記ビームフォーミングされた超音波信号を前記少なくとも一つの状況変数によって分類し、前記分類した結果に基づいて、前記点広がり関数データベースから少なくとも一つの点広がり関数を抽出するステップを含むことが好ましい。
前記少なくとも一つの点広がり関数を用いてデコンボリューションを行って、前記ビームフォーミングされた超音波信号から画像を生成するステップは、前記抽出された少なくとも一つの点広がり関数を用いて他の点広がり関数を推定するステップと、前記推定された他の点広がり関数を用いてデコンボリューションを行って、前記ビームフォーミングされた信号から画像を生成するステップとを含むことが好ましい。
前記少なくとも一つの状況変数は、前記超音波の音速、又は前記対象体内部の目標部位との距離であることが好ましい。
【発明の効果】
【００１３】
本発明に係る画像処理モジュール及び画像生成方法によれば、画像の復元、例えば、超音波信号に基づいた画像生成において、適切な点広がり関数を決定できるようになる。
従って、超音波画像装置による超音波画像生成過程において、迅速に点広がり関数を決定できるようになることで、超音波画像の復元速度が改善されて、迅速に超音波画像を生成することができる。
さらに、超音波画像の迅速な生成と共に、適切な点広がり関数を決定できるようになることで、高解像度の超音波画像を生成できるようになる効果もある。
また、１次元点広がり関数に基づいて高速に２次元点広がり関数を取得できるので、高解像度及び良質の超音波画像を容易に得ることができる。
従って、使用者、例えば、医師や超音波検査者が、超音波イメージング装置によって取得した迅速で且つ正確な高解像度の超音波画像を用いて、患者をより正確に診断できるようになる効果も得ることができる。
また、さらに、点広がり関数の決定方法は、レーダーや音響信号処理にも適用できるので、その応用実施形態によってレーダーや音響信号処理装置の動作を改善できるという効果も得ることができる。
【図面の簡単な説明】
【００１４】
【図１】本発明の実施形態による画像生成部を説明するためのブロック図である。
【図２】点広がり関数を説明するための図である。
【図３】原画像とＲＦ画像間の関係及びデコンボリューションを説明するための図である。
【図４ａ】原画像とＲＦ画像間の関係を説明するための図である。
【図４ｂ】深さによる目標部位（ｏｂ１）に対するＲＦ信号ベースの超音波画像の一例を示す図である。
【図４ｃ】目標部位の深さを説明するための超音波画像を示す図である。
【図５】本発明の実施形態による画像生成部の一実施例のブロック図である。
【図６】点広がり関数データベースの一実施例を説明するための図である。
【図７】本発明の実施形態による画像生成部の他の実施例のブロック図である。
【図８】音速による点広がり関数を説明するための図である。
【図９】音速による点広がり関数を説明するための図である。
【図１０】本発明の実施形態による画像生成部の更に他の実施例のブロック図である。
【図１１】本発明の実施形態による超音波画像装置の一実施例の斜視図である。
【図１２】超音波画像装置の一実施例の構成を示すブロック図である。
【図１３】超音波探針部の一実施例に係る平面図である。
【図１４】ビームフォーミング部の一実施例に係る構成図である。
【図１５】本発明の実施形態による超音波画像装置の画像生成部の構成の一実施例を示すブロック図である。
【図１６】第２点広がり関数を説明するための図である。
【図１７】入力画像に基づく原画像復元方法の様々な実施例を説明するためのフローチャートである。
【図１８】入力画像に基づく原画像復元方法の様々な実施例を説明するためのフローチャートである。
【図１９】入力画像に基づく原画像復元方法の様々な実施例を説明するためのフローチャートである。
【図２０】超音波画像装置を制御する方法の様々な実施例を説明するためのフローチャートである。
【図２１】超音波画像装置を制御する方法の様々な実施例を説明するためのフローチャートである。
【図２２】超音波画像装置を制御する方法の様々な実施例を説明するためのフローチャートである。
【発明を実施するための形態】
【００１５】
次に、本発明に係る画像処理モジュール及び画像生成方法を実施するための形態の具体例を図面を参照しながら説明する。
【００１６】
以下、図１〜図１０を参照して、画像生成部の一実施形態について説明する。
図１は、本発明の実施形態による画像生成部を説明するためのブロック図である。
図１に示すように、画像生成部１０は、入力信号ｄに基づいて画像信号ｍを生成して出力する。ここで、入力信号ｄは、音波や電磁波のような波動に起因して発生した信号であってもよい。
【００１７】
例えば、音波による信号である場合には、人間が聞くことができる周波数、すなわち、可聴周波数の音波から取得した信号であってもよく、可聴周波数２０ｋＨｚより大きい音波である超音波から取得した信号であってもよい。
また、電磁波による信号である場合、入力信号（ｄ）は、レーダーなどで用いられる極超短波（マイクロ波、波長１０ｃｍ〜１００ｃｍ）などから取得した信号であってもよい。
【００１８】
画像生成部１０は、入力信号（ｄ）から画像信号（ｍ）を生成するために、少なくとも一つの点広がり関数（ＰＳＦ、ｐｏｉｎｔ  ｓｐｒｅａｄ  ｆｕｎｃｔｉｏｎ、以下、図面上ではＰＳＦと表記）を推定し、推定された結果を用いてデコンボリューション（ｄｅｃｏｎｖｏｌｕｔｉｏｎ：逆畳み込み演算）を行うことによって、元の撮影しようとする対象体と同一、又はほぼ類似した画像を生成及び取得する。
画像生成部１０は、生成及び取得した画像を画像信号（ｍ）の形態で出力する。
【００１９】
点広がり関数は、画像撮影装置の撮影を通じて取得した画像データと組み合わされて、最終的な画像データを生成するための関数であって、理想的な画像データを復元するのに主に用いられる。
【００２０】
図２は、点広がり関数を説明するための図である。
図２に示すように、画像撮影装置が対象体に対する画像を取得する過程において、画像撮影装置の技術的性質や物理的特性、又は雑音（η）（ｎｏｉｓｅ）によって、元の画像である原画像（Ｏ）（ｉｄｅａｌ  ｉｍａｇｅ）とは異なる信号、例えば、超音波画像装置での超音波信号などのようなＲＦ信号（ｄ）（ｒａｄｉｏ  ｆｒｅｑｕｅｎｃｙ  ｓｉｇｎａｌ）を出力する。
言い換えると、画像撮影装置によって取得されるＲＦ信号（ｄ）は、原画像（Ｏ）に画像撮影装置の技術的性質や物理的特性によって変形され、ここに雑音（η）が付加されて出力される信号である。
【００２１】
図３は、一例として、超音波画像装置での原画像とＲＦ画像間の関係及びデコンボリューションを説明するための図である。
図３の最左側には、人体内部の組織に対する理想的な形状を示す。
図３に示すように、原画像（ｉｄｅａｌ）がｆＲのように与えられると、超音波画像装置において超音波プローブによって収集され、ビームフォーミングされた超音波画像は、中段のｇＲのように表現される。
すなわち、原画像とＲＦ信号による画像は互いに異なる。
これを、図４ａ〜図４ｃを参照して、より具体的に説明する。
【００２２】
図４ａは、原画像とＲＦ画像間の関係を説明するための図である。
図４ａに示す入力信号ベースの画像は、超音波画像装置を用いて取得した超音波画像の一例を図式的に表現したものである。
図４ａに示すように、仮に、理想的な状況での目標部位（ｏｂ１）に対する原画像（Ｏ）が、図４ａの左側にある図のように表示されるとすれば、目標部位（ｏｂ１）に対する入力信号（ｄ）、例えば、ＲＦ信号ベースの画像は、右側の図のように表示される。
【００２３】
具体的には、入力信号ベースの画像での目標部位（ｏｂ１）は、原画像での目標部位（ｏｂ１）が上下及び左右に伸びて表示される。
すなわち、入力信号ベースの画像は、原画像と非常に異なるものとなるため、もし、入力信号（ｄ）、すなわち、ＲＦ信号に基づいて画像をそのまま復元する場合には、目標部位が元の姿と異なるものとなる。
【００２４】
このような原画像と入力信号ベースの画像とは、深さなどによってさらに異なるものとなる。
図４ｂは、深さによる目標部位（ｏｂ１）に対するＲＦ信号ベースの超音波画像の一例を示したもので、図４ｃは、目標部位の深さを説明するための超音波画像を示す図である。
【００２５】
図４ｂに示すように、目標部位（ｏｂ１）と画像データ収集手段、例えば、超音波探針部間の距離が短いと、例えば、図４ｃに示すように、人体内部の病変が第１深さ（Ｄｅｐｔｈ＃１）に位置すると、目標部位（ｏｂ１）に対する入力信号ベースの画像は、同部位に対する原画像と同一、又は類似する。
【００２６】
しかし、仮に、画像データ収集手段と目標部位（ｏｂ１）との間の距離が離れると、例えば、人体内部の病変が、図４ｃの第４深さ（Ｄｅｐｔｈ＃４）や第５深さ（Ｄｅｐｔｈ＃５）などに位置すると、目標部位（ｏｂ１）の入力信号ベースの画像は側面方向に伸びて表示されて、目標部位（ｏｂ１）に対する原画像とかなりの相違が発生するようになる。すなわち、原画像（Ｏ）での目標部位（ｏｂ１）と、入力信号ベースの画像の目標部位（ｏｂ１）の形態は、データ収集手段と目標部位（ｏｂ１）との間の距離によってさらに異なるものとなる。
【００２７】
したがって、ＲＦ信号（ｄ）を用いて原画像（Ｏ）を復元する場合には、目標部位（ｏｂ１）に対する正確な画像を得るために、このような原画像（Ｏ）とＲＦ信号（ｄ）による画像との間の差を補正しなければならない。
この場合、原画像（Ｏ）と取得したＲＦ信号（ｄ）との間には所定の関係が成立するということを前提として、所定の関係に対応する所定の関数を用いてＲＦ信号（ｄ）を補正して画像を復元する。
このときに用いられる所定の関数が、点広がり関数（Ｈ）である。
【００２８】
図２に示した原画像（Ｏ）、点広がり関数（Ｈ）、雑音（η）、及び入力信号（ｄ）、すなわち、ＲＦ信号間の関係を数式で表現すると、下記の数式１のように表すことができる。
（数１）
ｄ＝Ｈｆ＋η
ここで、ｄは出力されるＲＦ信号、Ｈは点広がり関数、ｆは原画像に対する信号であり、ηは雑音を意味する。
【００２９】
仮に雑音がないと仮定すると、ＲＦ信号（ｄ）は、原画像（Ｏ）と点広がり関数（Ｈ）間の積で表現できる。
したがって、測定されたＲＦ信号ｄに対して適切な点広がり関数（Ｈ）を知れば、測定されたＲＦ信号（ｄ）に対応する原画像（Ｏ）を得ることができる。言い換えると、点広がり関数（Ｈ）及びＲＦ信号（ｄ）を知ると、対象体と同一、又は最も類似した画像を復元できるようになる。
【００３０】
画像生成部１０は、図２と反対方向に、ＲＦ信号、すなわち入力信号（ｄ）と、入力信号（ｄ）に適切な点広がり関数（Ｈ）を用いて原画像（Ｏ）と同一又は非常に類似した画像を生成し、生成された画像に対応する画像信号（ｍ）を出力する。言い換えると、画像生成部１０は、図３の中段に示すように、ＲＦ信号（ｇＲ）に適切な点広がり関数（ｈＲ）を付加及び合成してデコンボリューションを行うことによって、原画像（Ｏ、ｆＲ）と同一又はほぼ類似した復元画像（ｒｅｓｔｏｒｅｄ  ｉｍａｇｅ）を生成する。
【００３１】
画像生成部１０は、入力信号（ｄ）に基づいて対象体と同一又は最も類似した画像を復元するために、１次元又は２次元の点広がり関数で構築された点広がり関数データベース１２（ＰＳＦ  ｄａｔａｂａｓｅ）（図５参照）を用いたり、又は１次元の点広がり関数に基づいて２次元の点広がり関数を推定する等の方法を用いて、複数の方向での解像度低下などの問題点なしに、点広がり関数を推定するようにする。
【００３２】
そして、推定された点広がり関数（Ｈ）を用いて、対象体を撮影して取得した入力信号（ｄ）、例えば、ＲＦ信号を変換して、元の対象体の形態や形状と同一又は類似した形態や形状に見えるように入力信号（ｄ）を変換する。
すなわち、画像生成部１０は、点広がり関数を用いて、取得した入力信号（ｄ）に対する画像を復元して、適切に復元された画像に対する画像信号（ｍ）を出力する。
【００３３】
図５は、本発明の実施形態による画像生成部の一実施例のブロック図である。
図５に示すように、画像生成部１０の一実施例では、画像生成部１０は、点広がり関数（ＰＳＦ）選択部１１及びデコンボリューション部１４を含む。
点広がり関数選択部１１は、別途の点広がり関数データベース１２を閲覧し、点広がり関数データベース１２から少なくとも一つの点広がり関数を選択して呼び出す。
【００３４】
図６は、点広がり関数データベース１２の一実施例を説明するための図である。
点広がり関数データベース１２は、図６に示すように、少なくとも一つの点広がり関数の集合で構築される。
点広がり関数データベース１２の一実施例によれば、点広がり関数データベース１２は、１次元の点広がり関数のみを含んでもよく、２次元の点広がり関数のみを含むこともできる。また、１次元の点広がり関数及び２次元の点広がり関数を全て含んでいてもよい。実施例によっては、３次元や４次元などのさらに高い次元の点広がり関数を含んでいてもよい。
【００３５】
点広がり関数データベース１２に格納された点広がり関数は、少なくとも一つの状況変数によって予め演算された点広がり関数であってもよい。
ここで、一実施例によれば、状況変数は波動の速度（図６の水平方向）であってもよい。波動の速度は、例えば、超音波のような音波の速度（音速）であってもよい。
また、他の実施例によれば、状況変数は、データ収集手段、一例として、超音波探針部と目標部位との間の距離（図６の垂直方向）であってもよい。例えば、データ収集手段と目標部位間の距離は、人体内部の特定の目標部位、一例として、病変の深さ（ｄｅｐｔｈ）であってもよい。
【００３６】
もちろん、点広がり関数データベース１２に格納された点広がり関数は、二つの状況変数によって予め演算された点広がり関数であってもよい。
この場合、二つの状況変数は、図６に示すように、音速と目標部位の深さであってもよい。
言い換えると、図６に示すように、点広がり関数データベース１２は、音速及び深さという状況変数によって各々別途に取得された点広がり関数の集合であってもよい。
【００３７】
点広がり関数データベース１２は、２次元の点広がり関数を迅速に推定できるように、軸方向の１次元の点広がり関数はほぼ同一であると仮定し、様々な状況変数、例えば、上述した超音波音速や目標部位（ｏｂ１）の深さによって実測して取得した少なくとも一つの２次元の点広がり関数を含んでいてもよい。これについては後述する。
【００３８】
点広がり関数選択部１１は、上述したように構築することができる点広がり関数データベース１２から少なくとも一つの点広がり関数を選択する。
特に、点広がり関数選択部１１は、一実施例において、点広がり関数データベース１２から２次元の点広がり関数（２Ｄ  ｐｏｉｎｔ  ｓｐｒｅａｄ  ｆｕｎｃｔｉｏｎ）を選択することができる。
【００３９】
この場合、点広がり関数選択部１１は、入力信号（ｄ）に応じて少なくとも一つの点広がり関数を選択することができる。
仮に、目標部位（ｏｂ１）が、人体などの内部で深くない所、例えば、図４ｃの第１深さに位置していれば、点広がり関数選択部１１は、図５に示した点広がり関数データベース１２において第１深さに該当する点広がり関数を選択する。
仮に、図４ｃの第４深さや第５深さの目標部位（ｏｂ１）を撮影する場合には、点広がり関数データベース１２において第４深さや第５深さに該当する点広がり関数を選択するようにする。
【００４０】
点広がり関数選択部１１は、選択された点広がり関数、又はこれに関する情報をデコンボリューション部１４に伝達する。
デコンボリューション部１４は、点広がり関数選択部１１で選択された点広がり関数を用いて、入力信号（ｄ）をデコンボリューションする。
それによって、入力信号（ｄ）、すなわち、図２に示したＲＦ信号（ｄ）から、原画像（Ｏ）と同一又はほぼ類似した画像信号（ｍ）を取得する。
【００４１】
仮に、点広がり関数選択部１１が、点広がり関数データベース１２から２次元の点広がり関数を選択する場合であれば、デコンボリューション部１４は、選択された２次元の点広がり関数に基づいて、入力信号（ｄ）に対して２次元デコンボリューションを行うことができる。
もちろん、仮に、点広がり関数選択部１１が、点広がり関数データベース１２から１次元の点広がり関数を選択した場合であれば、デコンボリューション部１４は、入力信号（ｄ）に対して１次元デコンボリューションを行うようにすることも可能である。
【００４２】
これらを通じて、画像生成部１０は、入力信号（ｄ）から、原画像（Ｏ）と同一又は類似した画像を生成できるようになる。生成された画像は画像信号（ｍ）で出力される。
上述したように、画像生成部１０は、入力信号（ｄ）に応じて少なくとも一つの２次元の点広がり関数を点広がり関数データベース１２から選択し、選択した２次元の点広がり関数を用いてデコンボリューションを行うことによって、入力信号（ｄ）に基づいて原画像（Ｏ）と同一又は類似した画像を生成できる。
【００４３】
図２及び数式１に示したように、ＲＦ信号（ｄ）から原画像（Ｏ）を取得する過程は、点広がり関数に対する逆関数を求めることである。
したがって、点広がり関数が２次元である場合に、２次元の点広がり関数に対する逆関数を直接的に演算すれば、計算の複雑さが増加せざるを得ない。
また、１次元の点広がり関数を用いる場合には演算量は減少するが、１次元という限界があるため、一定領域での解像度低下の問題が発生する。
例えば、一方向、例えば、垂直方向の場合には解像度が低下しないが、他の方向、例えば、側面方向の場合には解像度が低下する可能性がある。
【００４４】
しかし、画像生成部１０は、点広がり関数データベース１２から入力信号（ｄ）に適切な２次元の点広がり関数を選択した後、選択した点広がり関数を用いて２次元デコンボリューションを行うので、逆関数を取得する複雑な演算過程を避けることができる。
したがって、画像の復元において、演算の単純化と解像度の改善の効果を共に得ることができる。
【００４５】
図７は、本発明の実施形態による画像生成部の他の実施例のブロック図である。
図７に示すように、画像生成部１０は、本実施例では、第１点広がり関数選択部１１ａ、第２点広がり関数推定部１３、及びデコンボリューション部１４を含む。
第１点広がり関数選択部１１ａは、上述した点広がり関数選択部１１と同様に、別途の点広がり関数データベース１２から少なくとも一つの点広がり関数を選択し、選択した少なくとも一つの点広がり関数を第２点広がり関数推定部１３に伝達する。
【００４６】
ここで、点広がり関数データベース１２は、一実施例によれば、少なくとも一つの状況変数によって予め演算された１次元の点広がり関数で構築することができる。
例えば、図４ｃに示したように、音速及び深さという状況変数によって各々別途に取得された点広がり関数によって点広がり関数データベース１２を構築することができる。
【００４７】
第１点広がり関数選択部１１ａの一実施例によれば、第１点広がり関数選択部１１ａは、点広がり関数データベース１２から少なくとも一つの１次元の点広がり関数を選択することができる。この場合、第１点広がり関数選択部１１ａは、２次元の点広がり関数を迅速に推定できるように、軸方向の１次元の点広がり関数は類似すると仮定し、様々な状況変数、例えば、上述した超音波音速や目標部位（ｏｂ１）の深さによる側面方向の１次元の点広がり関数を選択することもできる。第１点広がり関数選択部１１ａで選択された１次元の点広がり関数は、第２点広がり関数推定部１３に伝達される。
【００４８】
第２点広がり関数推定部１３は、伝達された１次元の点広がり関数を用いて２次元の点広がり関数を推定する。
このとき、第２点広がり関数推定部１３に伝達される点広がり関数は、側面方向の１次元の点広がり関数であってもよい。
１次元の点広がり関数は、軸（ａｘｉａｌ）方向又は側面（ｌａｔｅｒａｌ）方向に区分することができる。ここで、軸方向の１次元の点広がり関数の場合には、深さ又は音速によってほぼ類似した形態を示す。
【００４９】
図８及び図９は、音速及び深さによる点広がり関数を説明するための図である。
図８及び図９の左側の図は、超音波画像装置で用いられる１次元の点広がり関数（１Ｄ  ＰＳＦ）の一例をグラフの形態で示すものである。
図８は、超音波の音速が１５４０ｍ／ｓである場合に、図９は、超音波の音速が１４５０ｍ／ｓである場合に、互いに異なる深さの三つの目標部位に超音波を照射して取得した超音波信号に対する１次元の点広がり関数を示すグラフである。
【００５０】
図８と図９の左側に示すように、超音波の音速は互いに異なるにもかかわらず、全体的な１次元の点広がり関数の形態は互いに非常に類似することがわかる。
また、目標部位の深さが互いに異なる場合（ｄｅｐｔｈ１〜ｄｅｐｔｈ３）にも、１次元点広がり関数はほぼ類似したパターンで形成されている。
したがって、軸方向の１次元の点広がり関数の形態は、音速と目標部位の深さが互いに異なるにもかかわらず、ほぼ類似することがわかる。
【００５１】
図８及び図９の右側には、深さによる２次元の点広がり関数を示す。
図８と図９の左側の１次元の点広がり関数とは異なり、２次元の点広がり関数は、深さによって非常に異なることがわかる。すなわち、軸方向の点広がり関数とは異なり、側面方向の点広がり関数の場合には、深さや音速によってその形状が非常に異なる。
【００５２】
軸方向の１次元の点広がり関数は、超音波の音速や深さなどによって大きく変化しないので、第２点広がり関数推定部１３は、軸方向の１次元の点広がり関数は一定であると仮定し、第１点広がり関数選択部１１ａが選択した側面方向の１次元の点広がり関数に基づいて、２次元の点広がり関数を演算などの方法を通じて推定する。
【００５３】
一実施例によれば、第２点広がり関数推定部１３は、第１点広がり関数選択部１１ａが選択した側面方向の１次元の点広がり関数だけでなく、入力信号（ｄ）の特徴や、又は予め入力された設定や各種変数、すなわち、第２点広がり関数推定変数などをさらに考慮して第２点広がり関数を推定することもできる。
例えば、第２点広がり関数推定部１３で用いる第２点広がり関数推定変数は、入力信号（ｄ）のビームフォーミングに用いられたビームフォーミング係数であってもよい。
【００５４】
仮に、超音波画像装置に適用された場合であれば、第２点広がり関数推定変数は、超音波の音速、音速の変化量、対象体内部の目標部位との距離、入力信号生成モジュールの配置形態、及び入力信号のチャネル別減衰率などのような変数であってもよい。
第２点広がり関数推定部１３は、このような変数を用いて、より正確な２次元の点広がり関数を推定することができる。
【００５５】
デコンボリューション部１４は、伝達を受けた推定された２次元の点広がり関数を用いてデコンボリューションを行って、画像信号（ｍ）を生成及び出力する。
したがって、画像生成部１０は、入力信号（ｄ）に応じて点広がり関数データベース１２から少なくとも一つの１次元の点広がり関数を選択し、選択した１次元の点広がり関数に基づいて２次元の点広がり関数を推定した後、推定された２次元の点広がり関数によってデコンボリューションを行うことによって、入力信号（ｄ）に基づいて原画像（Ｏ）と同一又は類似した画像を生成することができる。
【００５６】
図１０は、本発明の実施形態による画像生成部の更に他の実施例のブロック図である。
図１０に示すように、画像生成部１０は、第１点広がり関数推定部１１ｂ、第２点広がり関数推定部１３、及びデコンボリューション部１４を含む。
【００５７】
第１点広がり関数推定部１１ｂは、入力信号（ｄ）に基づいて適切な１次元の点広がり関数を推定する。
１次元の点広がり関数を推定する場合には、２次元の点広がり関数を推定する場合より入力画像（ｄ）の変数を相対的に少なく考慮する。したがって、２次元の点広がり関数を直接推定することより点広がり関数の決定に必要な演算量が減少する。
【００５８】
第１点広がり関数推定部１１ｂは、一実施例によれば、一方向の点広がり関数を推定してもよく、また、複数の方向の点広がり関数を推定してもよく、一方向に適用する複数の点広がり関数を推定したり、又は複数の方向のそれぞれに適用される複数の点広がり関数を推定するようにしてもよい。
【００５９】
第１点広がり関数推定部１１ｂの一実施例によれば、第１点広がり関数推定部１１ｂは、少なくとも一つのフィルターを用いて第１点広がり関数を推定することができる。
この場合、少なくとも一つのフィルターは、最小二乗フィルター（ｌｅａｓｔ  ｓｑｕａｒｅ  ｆｉｌｔｅｒ）のような最小二乗フィルター（ＬＳＦ）やケプストラムフィルターであってもよい。
一実施例によれば、第１点広がり関数推定部１１ｂは、側面方向の点広がり関数を推定するようにすることもできる。
【００６０】
第２点広がり関数推定部１３は、第１点広がり関数推定部１１ｂで推定された第１点広がり関数を用いて、第２点広がり関数を推定する。
第２点広がり関数推定部１３は、一実施例によれば、図８及び図９を通じて上述したように、軸方向の点広がり関数は音速及び深さによってほぼ一定であると仮定し、第１点広がり関数推定部１１ｂで推定した側面方向の点広がり関数に基づいて、２次元の点広がり関数を演算などの方法を通じて推定する。
【００６１】
第２点広がり関数推定部１３は、図１０に示すように、第１点広がり関数選択部１１ａが選択した側面方向の１次元の点広がり関数だけでなく、入力信号（ｄ）の特徴や、又は予め入力された設定や各種変数などをさらに考慮して第２点広がり関数を推定することも可能である。
上述と同様に、第２点広がり関数推定部１３は、ビームフォーミング係数を用いて第２点広がり関数を推定することでき、また、仮に超音波画像装置に用いられる場合であれば、超音波の音速、対象体内部の目標部位との距離などのような他の変数を用いるようにして第２点広がり関数を推定することができる。
【００６２】
デコンボリューション部１４は、伝達を受けた推定された２次元の点広がり関数を用いてデコンボリューションを行って、画像信号（ｍ）を生成及び出力する。
上述したところにしたがって、画像生成部１０は、まず、別途のデータベースを閲覧せずに、少なくとも一つの１次元の点広がり関数を入力信号（ｄ）に応じて推定し、次いで、推定された１次元の点広がり関数に基づいて２次元の点広がり関数を推定した後に、推定された２次元の点広がり関数を用いて入力信号（ｄ）に対するデコンボリューションを行うことによって、入力信号（ｄ）に基づいて原画像（Ｏ）と同一又は類似した画像を復元することができる。
【００６３】
以下、図１１〜図１８を参照して、本発明の実施形態による超音波画像装置の一実施例について説明する。
図１１は、本発明の実施形態による超音波画像装置の一実施例の斜視図であり、図１２は、超音波画像装置の一実施例の構成を示すブロック図である。
【００６４】
図１１及び図１２に示すように、超音波画像装置の一実施例によれば、超音波画像装置は、超音波を用いて対象体ｏｂの内部に対する超音波信号を収集する超音波探針部（Ｐ）と、超音波探針部（Ｐ）が収集した超音波信号を収集して超音波画像を生成する本体（Ｍ）とを含む。
【００６５】
超音波探針部（Ｐ）は、印加される電源によって超音波を生成し、生成した超音波を対象体ｏｂ内部の目標部位（ｏｂ１）に照射する超音波発生部（Ｐ１１）と、対象体ｏｂの目標部位（ｏｂ１）で反射されたエコー超音波を受信して電気的信号に変換する超音波受信部（Ｐ１２）とを含む。
超音波発生部（Ｐ１１）は、本体（Ｍ）などに設置された超音波発生制御部３１０の制御によって超音波発生部（Ｐ１１）に印加されるパルス信号又は交流電流に応じて振動して超音波を生成する。
【００６６】
超音波発生部（Ｐ１１）で発生した超音波は、対象体内部の目標部位（ｏｂ１）で反射される。
超音波受信部（Ｐ１２）は、エコー超音波を受信し、エコー超音波の周波数に応じて振動して、受信したエコー超音波を所定の電気的信号に変換させる。
その結果、超音波受信部（Ｐ１２）は、超音波信号ｘを出力する。仮に、超音波画像装置の一実施例によって、超音波画像装置が光音響画像装置（ｐｈｏｔｏａｃｏｕｓｔｉｃ  ｉｍａｇｉｎｇ  ａｐｐａｒａｔｕｓ）と結合されている異種画像装置である場合、超音波受信部（Ｐ１２）は、レーザーなどの照射に起因して目標部位（ｏｂ１）から発生した音波、例えば、超音波を受信するようにしてもよい。
【００６７】
このような超音波発生部（Ｐ１１）及び超音波受信部（Ｐ１２）の機能は、超音波探針部（Ｐ）の先端に配置された超音波トランスデューサ（Ｐ１０）によって行うことができる。
図１３は、超音波探針部の一実施例に係る平面図である。
超音波トランスデューサ（Ｐ１０）は、図１３に示すように、超音波探針部（Ｐ）の一端に設置される。
【００６８】
トランスデューサとは、ある形態のエネルギー、一例として、電気エネルギーを他の形態のエネルギー、一例として、波動エネルギーや光エネルギーなどに変換する装置のことをいう。
超音波トランスデューサ（Ｐ１０）は、波動エネルギーと電気的エネルギーを相互変換する。
具体的には、超音波トランスデューサ（Ｐ１０）は、入力される所定のパルス電流に応じて振動して超音波を生成し、また、外部から受信した超音波、例えば、エコー超音波に応じて振動して所定パルスの電気的信号（以下、超音波信号）を生成する。これによって、超音波トランスデューサ（Ｐ１０）は、上述した超音波発生部（Ｐ１１）及び超音波受信部（Ｐ１２）の機能を全て行うことができる。
【００６９】
より詳細には、超音波トランスデューサ（Ｐ１０）は、外部の電源供給装置や、又は内部の蓄電装置、例えば、バッテリーなどのような電源３１１から交流電流が供給され、印加される交流電流に応じて超音波トランスデューサ（Ｐ１０）の圧電振動子や薄膜などが振動することによって超音波を生成する。
逆に、超音波トランスデューサ（Ｐ１０）は、超音波の受信によって上述した圧電物質や薄膜が振動しながら振動周波数に対応する周波数の交流電流を生成して、超音波を超音波信号に変換する。変換された超音波信号は、図１２及び図１３に示すように、複数のチャネル（Ｃ１〜Ｃ１０）を通じて本体部（Ｍ）のビームフォーミング部１００に伝達する。
【００７０】
図１３に示すように、複数個の超音波トランスデューサ（Ｐ１０）が超音波探針部（Ｐ）の先端に設置される。例えば、超音波探針部（Ｐ）の先端に６４個又は１２８個の超音波トランスデューサ（Ｐ１０）が設置されていてもよい。このように、複数の超音波トランスデューサ（Ｐ１０）が超音波探針部（Ｐ）の一端に設置された場合には、伝達される超音波信号もまた超音波トランスデューサ（Ｐ１０）の個数に対応する複数のチャネル、例えば、６４個又は１２８個（Ｃ１〜Ｃ１０）のチャネルでビームフォーミング部１００に伝達される。
【００７１】
上述した超音波トランスデューサ（Ｐ１０）としては、磁性体の磁歪効果を用いる磁歪超音波トランスデューサ（Ｍａｇｎｅｔｏｓｔｒｉｃｔｉｖｅ  Ｕｌｔｒａｓｏｎｉｃ  Ｔｒａｎｓｄｕｃｅｒ）や、圧電物質の圧電効果を用いた圧電超音波トランスデューサ（Ｐｉｅｚｏｅｌｅｃｔｒｉｃ  Ｕｌｔｒａｓｏｎｉｃ  Ｔｒａｎｓｄｕｃｅｒ）、微細加工された数百又は数千個の薄膜の振動を用いて超音波を送受信する静電容量型微細加工超音波トランスデューサ（ｃＭＵＴ、Ｃａｐａｃｉｔｉｖｅ  Ｍｉｃｒｏｍａｃｈｉｎｅｄ  Ｕｌｔｒａｓｏｎｉｃ  Ｔｒａｎｓｄｕｃｅｒ）などを用いることができる。また、それ以外に、電気的信号によって超音波を生成したり、又は超音波によって電気的信号を生成することができる他の種類のトランスデューサもまた、上述した超音波トランスデューサ（Ｐ１０）の一例に該当する。
【００７２】
本体（Ｍ）は、一実施例によれば、ビームフォーミング部１００、画像生成部２００、点広がり関数データベース２２０、システム制御部３００、及び画像処理部３２０をさらに含むことができる。
【００７３】
本体（Ｍ）のビームフォーミング部１００は、超音波受信部（Ｐ１２）から伝達される超音波信号を受信し、受信した超音波信号をビームフォーミングする。
具体的には、ビームフォーミング部１００は、複数のチャネルデータに対して、特定空間の反射波の大きさを推定するためのビームフォーミングを行う。
【００７４】
図１４は、ビームフォーミング部の一実施例に係る構成図である。
図１４に示すように、ビームフォーミング部１００は、時差補正部１１０及び集束部１２０を含む。
図１４を参照すると、目標部位（ｏｂ１）で反射して戻ってくるエコー超音波は、超音波受信部（Ｐ１１）、一例として、超音波トランスデューサ（Ｐ１０）によって受信される。
【００７５】
ところで、超音波探針部（Ｐ）に設置されたそれぞれの超音波トランスデューサ（Ｔ１〜Ｔ８）と、目標部位（ｏｂ１）との間の距離は互いに異なり、音速はほぼ一定であるため、エコー超音波が同一の目標部位（ｏｂ１）で反射して戻ってくる場合であるとしても、それぞれの超音波トランスデューサ（Ｔ１〜Ｔ８）は、互いに異なる時間に同一の目標部位（ｏｂ１）のエコー超音波を受信するようになる。
言い換えると、同一の目標部位（ｏｂ１）のエコー超音波によってそれぞれの超音波トランスデューサ（Ｔ１〜Ｔ８）で出力する超音波信号間には所定の時間差が存在する。これによって、各超音波トランスデューサ（Ｔ１〜Ｔ８）によって発生した超音波信号間の時間差を先に補正する必要が発生する。
【００７６】
ビームフォーミング部１００の時差補正部１１０は、このような超音波信号間の時間差を補正する。
例えば、時差補正部１１０は、図１３及び図１４に示すように、特定のチャネルに入力される超音波信号の伝送を、予め定められたところにしたがって一定程度遅延させることで、各チャネルに入力される超音波信号が同じ時間に集束部１２０に伝達されるようにする。
【００７７】
ビームフォーミング部１００の集束部１２０は、時間差が補正された超音波信号を集束するようにする。
集束部１２０は、入力された各超音波信号毎に所定の加重値、すなわち、ビームフォーミング係数を付加して特定位置の信号を強調したり、他の位置の信号は相対的に減衰させて超音波信号を集束するようにする。
これによって、使用者の要求事項や使用者の便宜に基づく超音波画像を生成できるようになる。
【００７８】
この場合、集束部１２０は、超音波受信部（Ｐ１２）によって出力される超音波信号と関係なく、定められたビームフォーミング係数を用いて超音波信号を集束することができる（データ独立型ビームフォーミング方式）。
また、入力される超音波信号に基づいて最適のビームフォーミング係数を演算した後、演算されたビームフォーミング係数を用いて超音波信号を集束することもできる（適応型ビームフォーミング方式）。
【００７９】
このようなビームフォーミングプロセスは、以下の数式２のように表現することができる。
【数２】
2014023934.tif 000003
【００８０】
ここで、ｎは、目標部位（ｏｂ１）の位置、例えば、目標部位（ｏｂ１）の特定の深さを意味する値であり、ｍは、目標部位（ｏｂ１）の位置ｎにおいてｍ番目のチャネルを意味し、ｗｍは、ｍ番目のチャネルの超音波信号に付加されるビームフォーミング係数ｗを意味する。
一方、Δｍは、上述した時差補正部１１０で行われる超音波信号の伝送時間の遅延に用いられる時間遅延値である。上述した数式２によって、集束部１２０は、時差が補正された各チャネルの超音波信号を集束して、ビームフォーミングされた超音波信号（ｄ）を出力する。
ビームフォーミング部１００でビームフォーミングされた超音波信号（ｄ）は、図１２及び図１４に示したように、画像生成部２００に伝達される。
【００８１】
図１５は、本発明の実施形態による超音波画像装置の画像生成部の構成の一実施例を示すブロック図である。
図１５に示すように、画像生成部２００は、本実施例によれば、点広がり関数選択部２１０及びデコンボリューション部２４０を含む。
【００８２】
点広がり関数選択部２１０は、図５に示したような点広がり関数データベース１２から少なくとも一つの点広がり関数を選択する。
この場合、点広がり関数データベース２２０は、１次元の点広がり関数及び／又は２次元の点広がり関数を含んでいてもよい。
【００８３】
上述したように、超音波画像装置において、超音波画像の復元の際に点広がり関数に最も影響を及ぼす要素は、超音波の速度、すなわち、音速、及び超音波が発生又は反射される目標部位と収集手段との間の距離、例えば、人体内部の病変の深さである。
したがって、点広がり関数データベース２２０は、超音波画像復元の点広がり関数に重要な影響を及ぼす超音波の速度と目標部位の深さを状況変数として予め算出、測定された１次元又は２次元の点広がり関数を含み得る。
【００８４】
画像生成部２００の一実施例において、仮に、点広がり関数選択部２１０が１次元点広がり関数を選択する場合であれば、画像生成部２００は、１次元点広がり関数選択部２１０が選択した１次元点広がり関数に基づいて２次元点広がり関数を推定する第２点広がり関数推定部２３０をさらに含むことができる。
【００８５】
この場合、点広がり関数選択部２１０は、ビームフォーミングされた信号（ｄ）に基づいて対象体（ｏｂ）内の目標部位（ｏｂ１）の深さや超音波の音速によって、点広がり関数データベース２２０から適切な１次元の点広がり関数を選択する。
この場合、１次元の点広がり関数は、目標部位（ｏｂ１）の側面方向に対する点広がり関数であってもよい。
【００８６】
第２点広がり関数推定部２３０は、点広がり関数選択部２１０によって選択された１次元の点広がり関数に基づいて２次元の点広がり関数を推定する。
上述したように、軸方向の点広がり関数は、音速や目標部位（ｏｂ１）の深さが変化しても大きく相違しないので、第２点広がり関数推定部２３０は、軸方向の１次元の点広がり関数は一定であると仮定し、点広がり関数選択部２１０が選択した１次元の点広がり関数、一例として、側面方向の点広がり関数を用いて２次元の点広がり関数を推定して取得する。
【００８７】
図１６は、第２点広がり関数を説明するための図である。
第２点広がり関数推定部２３０は、図１６の左側に示すように、選択された１次元の点広がり関数と別途の係数、例えば、第２点広がり関数推定係数を組み合わせて、図１６の右側に示す第２点広がり関数を推定する。
【００８８】
第２点広がり関数推定係数は、図１６に示すように、ビームフォーミング係数であってもよく、より具体的に、このようなビームフォーミング係数は、超音波トランスデューサ（Ｐ１０）の配置構造、音速、深さ、信号の減衰率などの内の少なくとも一つを用いて演算されたものであってもよい。
また、第２点広がり関数推定係数は、ビームフォーミング係数の代わりに、上述した超音波トランスデューサｐ１０の配置構造、音速、深さ、信号の減衰率のうち少なくとも一つであってもよい。そのために、第２点広がり関数推定部２３０は、ビームフォーミングされた超音波信号（ｄ）を受信したり、またはビームフォーミングされた超音波信号（ｄ）に関連する情報をビームフォーミング部１００から受信することもできる。
【００８９】
第２点広がり関数推定部２３０は、上述したように推定された２次元点広がり関数をデコンボリューション部２４０に伝達して、ビームフォーミングされた信号（ｄ）に対する２次元デコンボリューションを行うことができるようにする。
画像生成部２００の他の実施例によれば、点広がり関数選択部２１０は、点広がり関数データベース２２０から入力されるビームフォーミングされた超音波信号（ｄ）に対応する適切な２次元の点広がり関数を選択することもできる。
この場合、上述した点広がり関数データベース２２０の少なくとも一つの２次元の点広がり関数は、軸方向の１次元の点広がり関数は類似すると仮定した状態で、様々な状況変数によって実測して取得した点広がり関数であってもよい。
【００９０】
点広がり関数選択部２１０は、２次元の点広がり関数を選択して決定した後、決定された点広がり関数をデコンボリューション部２４０に直接伝達して、ビームフォーミングされた信号（ｄ）に対する２次元デコンボリューションを行うことができるようにする。
【００９１】
画像生成部２００の更に他の実施例によれば、画像生成部２００は、第１点広がり関数推定部２１１及び第２点広がり関数推定部２３０を含むこともできる。
第１点広がり関数推定部２１１は、点広がり関数データベース２２０の閲覧なしに、ビームフォーミングされた超音波信号（ｄ）に基づいて１次元の点広がり関数を推定する。
この場合、第１点広がり関数推定部２１１は、最小二乗フィルターやケプストラムフィルターのような少なくとも一つのフィルターを用いて１次元の点広がり関数を推定することができる。
【００９２】
推定された１次元の点広がり関数は、第２点広がり関数推定部２３０に伝達され、第２点広がり関数推定部２３０は、推定された１次元の点広がり関数に基づいて２次元の点広がり関数を推定する。上述と同様に、第２点広がり関数推定部２３０は、軸方向の１次元点広がり関数は一定であると仮定し、第１点広がり関数推定部２１１によって推定された１次元の点広がり関数を用いて２次元の点広がり関数を推定することもできる。
第２点広がり関数推定部２３０は、デコンボリューション部２４０に、推定された２次元の点広がり関数を伝達する。
【００９３】
デコンボリューション部２４０は、点広がり関数選択部２１０又は第２点広がり関数推定部２３０から伝達された２次元の点広がり関数を用いて、ビームフォーミングされた信号（ｄ）に対する２次元デコンボリューションを行う。そして、原画像（Ｏ）と同一又は類似した画像を復元する。
復元された画像に対する画像信号（ｍ）は、画像処理部３２０、格納部３２１、又はディスプレイ部（ｄｐ）に伝達する。
画像処理部３２０は、画像生成部２００で生成された画像に基づいて、使用者が識別できるモードの超音波画像を生成したり、又は画像生成部２２０で生成された画像に所定の画像処理を行うことができる。
【００９４】
画像処理部３２０の一実施例によれば、画像処理部３２０は、画像生成部２００で出力される画像信号（ｍ）に基づいて様々なモードの超音波画像を生成することもできる。
このときに使用できる超音波画像モードとしては、Ａモード及びＢモードなどがあり得る。
Ａモードは、振幅（ａｍｐｌｉｔｕｄｅ）として超音波画像を表示する方法であって、目標部位（ｏｂ１）を超音波探針部（Ｐ）からの距離又は時間として表示し、反射の強度を振幅として表示するようにするモードである。
Ｂモードは、エコー超音波の大きさを明るさ（ｂｒｉｇｈｔｎｅｓｓ）として画面に表示するようにするモードである。超音波画像がＢモードで生成された場合には、使用者が超音波画像のみでも直観的に対象体（ｏｂ）内部の組織や構造を容易に把握できるので、広く使用される。図４ｃに示した画像がＢモード下での超音波画像である。
【００９５】
さらに、画像処理部３２０は、使用者の意図にしたがって、又は使用者の便宜のために画像生成部２００で生成された超音波画像を補正することができる。
例えば、画像処理部３２０は、使用者が超音波画像内の組織を明確に見ることができるように、明度や輝度、対照度（コントラスト）、色相を補正することもできる。仮に、画像生成部２００で複数の超音波画像が出力される場合には、出力される複数の超音波画像を用いて３次元立体超音波画像を生成することも可能である。
【００９６】
格納部３２１は、画像生成部２００で生成された画像信号（ｍ）、画像処理部３２０で補正された画像に対する画像信号などを一時的又は非一時的に格納する。
ディスプレイ部（ｄｐ）は、超音波画像を使用者に表示する。
実施例によっては、ディスプレイ部（ｄｐ）は、画像生成部２００でデコンボリューション処理された超音波画像を直接使用者に表示してもよく、画像生成部２００で生成された画像に対して画像処理部３２０が所定の画像処理を行った超音波画像を使用者に表示してもよい。この場合、ディスプレイ部（ｄｐ）に表示される超音波画像は、Ｂモードの超音波画像であってもよく、また、３次元立体超音波画像であってもよい。
【００９７】
一方、超音波画像装置の本体（Ｍ）は、超音波発生制御部３１０をさらに含むことができる。
超音波発生制御部３１０は、システム制御部３００などの制御命令を受信し、受信された制御命令に応じて所定のパルス信号を生成して超音波発生部（Ｐ１１）に伝達する。
超音波発生部（Ｐ１１）は、伝達されたパルス信号に応じて振動して超音波を生成する。また、超音波発生制御部３１０は、超音波発生部（Ｐ１１）と電気的に接続された電源３１１に対する別途の制御信号を生成し、制御信号を電源３１１に伝達することができる。電源３１１は、伝達された制御信号に応じて超音波発生部（Ｐ１１）に所定の交流電流を印加する。
【００９８】
システム制御部３００は、超音波探針部（Ｐ）、ビームフォーミング部１００、画像生成部２００、画像処理部３２０、格納部３２１、ディスプレイ部（ｄｐ）などの全般的な動作を制御する。
実施例によっては、システム制御部３００は、既に定められた設定によって、又は別途の入力部（ｉ）を通じて入力される使用者の指示又は命令に応じて上述した超音波探針部（Ｐ）などに対する所定の制御命令を生成して超音波画像装置を制御することができる。
【００９９】
入力部（ｉ）は、超音波画像装置の制御のために、使用者から所定の指示や命令が入力される。
入力部（ｉ）は、例えば、キーボード（ｋｅｙｂｏａｒｄ）、マウス（ｍｏｕｓｅ）、トラックボール（ｔｒａｃｋｂａｌｌ）、タッチスクリーン（ｔｏｕｃｈ  ｓｃｒｅｅｎ）、またはパドル（ｐａｄｄｌｅ）などの様々なユーザーインターフェースを含むことができる。また、入力部（ｉ）は、本体（Ｍ）と別途に接続されたワークステーション（ｗｏｒｋ  ｓｔａｔｉｏｎ）であってもよい。
【０１００】
以上、超音波画像装置の実施例を説明するにあたって、超音波トランスデューサ（Ｐ１０）などは超音波探針部（Ｐ）に設置され、ビームフォーミング部１００、画像生成部２００、超音波発生制御部３１０や画像処理部３２０などの構成要素は本体（Ｍ）に設置された一例を挙げて説明したが、実施例によっては、ビームフォーミング部１００、画像生成部２００、超音波発生制御部３１０や画像処理部３２０などの本体（Ｍ）に設置された一部の構成要素が超音波探針部（Ｐ）に設けられてもよい。
【０１０１】
以下、図１７〜図１９を参照して、入力画像に基づく原画像復元方法の様々な実施例について説明する。
図１７は、点広がり関数を用いた画像生成方法の一実施例を説明するためのフローチャートである。
【０１０２】
図１７に示すように、まず、入力信号（ｄ）が入力される。ここで、入力信号（ｄ）は、一実施例によれば、超音波画像装置などのビームフォーミングされた信号、一例として、ビームフォーミングされた超音波信号であってもよい（ステップｓ４１０）。
【０１０３】
次いで、入力された入力信号ｄを分析する（ステップｓ４１１）。
この場合、入力信号（ｄ）の分析のために、入力信号（ｄ）に別途に付加されるメタ情報を用いることも可能である。入力信号（ｄ）を分析するにあたって、入力信号（ｄ）を取得するために用いられた超音波の音速や、又はデータ収集手段と目標部位（ｏｂ１）との距離、一例として、目標部位（ｏｂ１）の深さなどの状況変数が分析の基準になることができる。
【０１０４】
入力信号（ｄ）の分析結果によって、適切な２次元の点広がり関数を選択する。
このときに選択される２次元の点広がり関数は、別途の点広がり関数データベース１２に格納された複数の点広がり関数のうち少なくとも一つであってもよい。
選択された２次元の点広がり関数は、点広がり関数データベース１２から呼び出される（ステップｓ４１２）。
【０１０５】
点広がり関数データベース１２に格納された点広がり関数は、図６に示したように、超音波の音速や、又はデータ収集手段と目標部位（ｏｂ１）との距離などの状況変数によって予め測定された点広がり関数であってもよい。
選択された２次元点の広がり関数を用いて、入力信号（ｄ）に対して２次元デコンボリューションを行い（ステップｓ４１３）、その結果、原画像と同一又はほぼ類似した画像信号（ｍ）を取得して出力する（ステップｓ４１４）。
【０１０６】
図１８は、点広がり関数に基づいて画像を生成する方法の他の実施例を説明するためのフローチャートである。
図１８を参照すると、まず、入力信号（ｄ）が入力され（ステップｓ４２０）、上述と同様に入力信号（ｄ）が分析される（ステップｓ４２１）。
同様に、入力信号（ｄ）はビームフォーミングされた信号であってもよい。また、入力信号（ｄ）の分析は、上述したような状況変数を基準に行われてもよい。状況変数の一例として、音速又は深さなどがあり得る。
【０１０７】
入力信号（ｄ）の分析結果によって、適切な１次元の点広がり関数を選択する。
このときに選択される１次元の点広がり関数は、点広がり関数データベース１２に格納された少なくとも一つの点広がり関数であってもよい。点広がり関数データベース１２に格納された点広がり関数は、一実施例によれば、超音波の音速や、又はデータ収集手段と目標部位（ｏｂ１）との距離などの状況変数によって予め測定された１次元の点広がり関数であってもよい。例えば、１次元の点広がり関数は、側面方向に対する点広がり関数であってもよい。
【０１０８】
選択された１次元の点広がり関数は、点広がり関数データベース１２から呼び出される（ステップｓ４２２）。
選択、呼び出された点広がり関数を用いて２次元の点広がり関数を推定する。
上述したように、軸方向の点広がり関数は、音速や距離などによって大きく変化しないので、側面方向の点広がり関数に基づいて２次元の点広がり関数を推定できる（ステップｓ４２３）。
【０１０９】
一方、実施例によっては、２次元の点広がり関数の推定のために予め定義したり、又は使用者によって選択される２次元の点広がり関数推定変数を用いることができる。
このような２次元の点広がり関数推定変数は、一例として、ビームフォーミング係数であってもよい。より具体的に、第２点広がり関数推定変数は、超音波の音速、音速の変化量、対象体内部の目標部位との距離、入力信号生成モジュールの配置形態、及び入力信号のチャネル別減衰率などであってもよい。
【０１１０】
推定された２次元の点広がり関数を用いて、入力信号（ｄ）に対するデコンボリューションを行い（ステップｓ４２４）、画像を復元する。
すなわち、復元された画像に対する画像信号を取得する（ステップｓ４２５）。
【０１１１】
図１９は、点広がり関数に基づいて画像を生成する方法の更に他の実施例を説明するためのフローチャートである。
図１９に示すように、まず、入力信号（ｄ）が入力される（ステップｓ４３０）。
上述したように、入力信号（ｄ）はビームフォーミングされた信号であってもよい。
【０１１２】
入力信号（ｄ）の分析結果によって、１次元の点広がり関数を推定する（ステップｓ４３１）。
このとき、入力信号（ｄ）に対して少なくとも一つのフィルターを用いて１次元の点広がり関数を推定することも可能である。
利用可能なフィルターとしては、最小二乗フィルターやケプストラムフィルターなどがあり得る。
【０１１３】
上述したように、１次元の点広がり関数が推定すると、推定した１次元の点広がり関数に基づいて２次元の点広がり関数を推定する（ステップｓ４３２）。
この場合、図１８を通じて説明した実施例と同様に、２次元の点広がり関数を推定するようにすることができる。
２次元点広がり関数を取得すると、取得した２次元の点広がり関数を用いて入力信号（ｄ）に対して２次元デコンボリューションを行う（ステップｓ４３３）。
その結果、画像は復元され、復元された画像に対する画像信号（ｍ）を取得できる（ステップｓ４３４）。
【０１１４】
次に、図２０〜図２２を参照して、超音波画像装置の制御方法の様々な実施例について説明する。
図２０は、超音波画像装置を制御する方法の一実施例を説明するためのフローチャートである。
【０１１５】
図２０に示すように、まず、超音波画像装置の超音波探針部（Ｐ）が対象体（ｏｂ）の目標部位（ｏｂ１）に超音波を照射する（ステップｓ４４０）。
目標部位（ｏｂ１）は、照射された超音波の全部又は一部を反射させる。反射して戻ってくる超音波、すなわち、エコー超音波を超音波探針部（Ｐ）が受信し、電気的信号に変換して超音波信号を出力する（ステップｓ４４１）。
【０１１６】
出力された超音波信号の時差を補正した後、集束して、超音波信号に対するビームフォーミングを行う（ステップｓ４４２）。
ビームフォーミングされた超音波信号を分析した後（ステップｓ４４３）、ビームフォーミングされた超音波信号（ｄ）に適切な点広がり関数を選択する。
【０１１７】
この場合、ビームフォーミングされた超音波信号（ｄ）の分析のために、一例として、超音波の音速や目標部位（ｏｂ１）の深さなどのような状況変数を用いることができる。
ビームフォーミングされた超音波信号（ｄ）の分析のために、ビームフォーミングされた超音波信号（ｄ）だけでなく、超音波信号（ｄ）に付加された各種情報を用いることができる。分析の結果、例えば、超音波信号の取得のために用いられた超音波の音速及び目標部位（ｏｂ１）の深さなどが決定され得る。
【０１１８】
次に、分析結果によって適切な点広がり関数を選択する。
選択される点広がり関数は、点広がり関数データベース２２０に格納されていてもよい。例えば、分析によって、超音波信号の取得のために用いられた超音波の音速及び目標部位（ｏｂ１）の深さが決定されると、図６に示したような点広がり関数データベースから、決定された音速及び深さに対応する点広がり関数を選択する。点広がり関数データベース２２０から点広がり関数を呼び出す（ステップｓ４４４）。
【０１１９】
呼び出された点広がり関数を用いて、ビームフォーミングされた超音波データ（ｄ）に対する２次元デコンボリューションを行う（ステップｓ４４５）。
その結果、超音波画像が復元及び取得される（取得ｓ４４６）。
必要に応じて、復元及び取得された超音波画像に対して画像の明度や対照度（コントラスト）などを調節する画像後処理をさらに行った後（ステップｓ４４７）、画像処理された超音波画像をモニタなどのようなディスプレイ部（ｄｐ）に出力する（ステップｓ４４８）。
【０１２０】
図２１は、超音波画像装置を制御する方法の他の実施例を説明するためのフローチャートである。
図２１に示すように、超音波探針部（Ｐ）が超音波を対象体（ｏｂ）の目標部位（ｏｂ１）に照射し（ステップｓ４５０）、目標部位（ｏｂ１）で反射されるエコー超音波を受信した後、これを超音波信号に変換する（ステップｓ４５１）。
次いで、超音波信号に対して時差補正及び集束などのビームフォーミングを行う（ステップｓ４５２）。
【０１２１】
ビームフォーミングされた超音波信号を分析し（ステップｓ４５３）、分析結果によって適切な１次元の点広がり関数を選択する。
１次元の点広がり関数は、点広がり関数データベース２２０に格納されている点広がり関数であってもよい。また、適切な１次元の点広がり関数を選択するために、目標部位（ｏｂ１）の位置、例えば、対象体（ｏｂ）内部に目標部位（ｏｂ１）が位置する深さ、一例として、人体内部の病変の深さに対する情報を用いることができる。また、超音波探針部（Ｐ）が生成した超音波又は受信したエコー超音波の音速に対する情報を用いることもできる。一方、選択される１次元の点広がり関数は側面方向の点広がり関数であってもよい（ステップｓ４５４）。
【０１２２】
１次元の点広がり関数に基づいて２次元の点広がり関数を推定する（ステップｓ４５５）。
上述したように、軸方向の点広がり関数は音速別にほぼ類似するので、実際の画像復元のための軸方向の点広がり関数は、予め定義された軸方向の点広がり関数と同一であると仮定し、選択された側面方向の点広がり関数を用いて２次元の点広がり関数を演算して推定する。
【０１２３】
２次元の点広がり関数を推定する一実施例によれば、予め定義したり、又は使用者によって選択される２次元の点広がり関数推定変数、一例として、ビームフォーミング係数を用いて２次元の点広がり関数を推定することができる。
超音波の音速、音速の変化量、対象体内部の目標部位との距離、入力信号生成モジュールの配置形態、及び入力信号のチャネル別減衰率などもまた第２点広がり関数推定変数として用いることができる。
【０１２４】
そして、２次元の点広がり関数を用いて、ビームフォーミングされた超音波信号（ｄ）に対する２次元デコンボリューションを行って（ステップｓ４５６）、超音波画像を復元する（ステップｓ４５７）。
復元された超音波画像は、所定の画像後処理を通じて補正され（ステップｓ４５８）、ディスプレイ部を通じて出力される（ステップｓ４５９）。
【０１２５】
図２２は、超音波画像装置を制御する方法の更に他の実施例を説明するためのフローチャートである。
図２２に示すように、超音波探針部（Ｐ）が超音波を対象体（ｏｂ）の目標部位（ｏｂ１）に照射し（ステップｓ４６０）、目標部位（ｏｂ１）で反射されるエコー超音波を受信した後、これを超音波信号に変換する（ステップｓ４６１）。
次いで、超音波信号に対して時差補正及び集束などのビームフォーミングを行う（ステップｓ４６２）。この過程は上述した実施例と同一である。
【０１２６】
ビームフォーミングされた超音波信号（ｄ）を分析して（ステップｓ４６３）、ビームフォーミングされた超音波信号（ｄ）に対する１次元の点広がり関数を推定する（ステップｓ４６４）。
ビームフォーミングされた超音波信号（ｄ）の分析及び１次元の点広がり関数の推定のために、少なくとも一つのフィルター、一例として、最小二乗フィルターやケプストラムフィルターを用いて１次元の点広がり関数を推定することができる。その他にも、１次元の点広がり関数の推定のための従来の様々な１次元の点広がり関数推定方法を用いることもできる。このとき、軸方向の１次元の点広がり関数の代わりに、ただ側面方向の１次元の点広がり関数のみを推定することもできる。
【０１２７】
推定された１次元の点広がり関数を用いて２次元の点広がり関数を推定する（ステップｓ４６５）。
上述した実施例と同様に、実際の画像復元のための軸方向の点広がり関数は、予め定義された軸方向の点広がり関数と同一であると仮定し、推定された側面方向の点広がり関数を用いて２次元の点広がり関数を演算して推定する。
この場合、２次元の点広がり関数推定のために、上述したような２次元の点広がり関数推定係数を用いるようにしてもよい。
【０１２８】
そして、２次元の点広がり関数を用いて、ビームフォーミングされた超音波信号（ｄ）に対する２次元デコンボリューションを行う（ステップｓ４６６）。
その結果、原画像（Ｏ）と同一又は非常に類似した超音波画像が復元される（ステップｓ４６７）。
復元された超音波画像は、所定の画像後処理を通じて補正することができる（ステップｓ４６８）。
復元された超音波画像又は後処理された超音波画像は、ディスプレイ部を通じて出力される（ステップｓ４６９）。
【０１２９】
尚、本発明は、上述の実施例に限られるものではない。本発明の技術的範囲から逸脱しない範囲内で多様に変更実施することが可能である。
【符号の説明】
【０１３０】
１０、２００    画像生成部
１１ａ    第１点広がり関数選択部
１１ｂ、２１１    第１点広がり関数推定部
１２、２２０    点広がり関数データベース
１３、２３０    第２点広がり関数推定部
１４、２４０    デコンボリューション部
１００    ビームフォーミング部
１１０    時差補正部
１２０    集束部
２１０    点広がり関数選択部
２１１    第１点広がり関数推定部
３００    システム制御部
３１０    超音波発生制御部
３１１    電源
３２０    画像処理部
３２１    格納部
Ｍ    本体
Ｐ    超音波探針部
Ｐ１０    超音波トランスデューサ
Ｐ１１    超音波発生部
Ｐ１２    超音波受信部
ｉ    入力部
ｄｐ    ディスプレイ部
【図１】
2014023934.tif 000004
【図２】
2014023934.tif 000005
【図３】
2014023934.tif 000006
【図４ａ】
2014023934.tif 000007
【図４ｂ】
2014023934.tif 000008
【図４ｃ】
2014023934.tif 000009
【図５】
2014023934.tif 000010
【図６】
2014023934.tif 000011
【図７】
2014023934.tif 000012
【図８】
2014023934.tif 000013
【図９】
2014023934.tif 000014
【図１０】
2014023934.tif 000015
【図１１】
2014023934.tif 000016
【図１２】
2014023934.tif 000017
【図１３】
2014023934.tif 000018
【図１４】
2014023934.tif 000019
【図１５】
2014023934.tif 000020
【図１６】
2014023934.tif 000021
【図１７】
2014023934.tif 000022
【図１８】
2014023934.tif 000023
【図１９】
2014023934.tif 000024
【図２０】
2014023934.tif 000025
【図２１】
2014023934.tif 000026
【図２２】
2014023934.tif 000027