(57)【特許請求の範囲】
【請求項１】
診断領域に対して超音波を送受するプローブと、
プローブを制御することにより前記診断領域から受信信号を得る送受信部と、
受信信号に基づいて形成される前記診断領域の画像データ内において対象物を特定する対象物特定部と、
前記画像データ内において対象物を複数の対象部位に区分けする区分処理部と、
前記画像データ内において複数の対象部位の各々に対して各対象部位に応じた処理回数だけ段階的な画像処理を実行する画像処理部と、
を有する、
ことを特徴とする超音波診断装置。
【請求項２】
請求項１に記載の超音波診断装置において、
前記画像データは、受信信号から得られる複数のサンプルデータを割り当てられた複数の画素により構成され、
前記各対象部位に応じた処理回数は、その対象部位を構成する画素の大きさに基づいて決定される、
ことを特徴とする超音波診断装置。
【請求項３】
請求項２に記載の超音波診断装置において、
前記各対象部位内の代表位置における画素の大きさに基づいて、その対象部位に応じた前記処理回数を設定する処理回数設定部をさらに有する、
ことを特徴とする超音波診断装置。
【請求項４】
請求項３に記載の超音波診断装置において、
前記処理回数設定部は、前記画像データ内の基準位置における画素の体積と前記各対象部位内の代表位置における画素の体積とから算出される体積比に基づいて、その対象部位に応じた前記処理回数を設定する、
ことを特徴とする超音波診断装置。
【請求項５】
請求項１から４のいずれか１項に記載の超音波診断装置において、
前記画像処理部は、複数の対象部位の各々に対して各対象部位に応じた前記処理回数だけ段階的な収縮処理を実行する、
ことを特徴とする超音波診断装置。
【請求項６】
請求項５に記載の超音波診断装置において、
前記画像処理部は、収縮処理された複数の対象部位の各々に対して各対象部位に応じた前記処理回数だけ段階的な膨張処理を実行する、
ことを特徴とする超音波診断装置。
【請求項７】
超音波の画像データ内において対象物を特定する対象物特定部と、
前記画像データ内において対象物を複数の対象部位に区分けする区分処理部と、
前記画像データ内において複数の対象部位の各々に対して各対象部位に応じた処理回数だけ段階的な画像処理を実行する画像処理部と、
を有する、
ことを特徴とする超音波画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波診断装置に関し、特に画像処理の適正化に関する。
【背景技術】
【０００２】
超音波画像内において組織等に対して診断に適した画像処理を施す技術が従来からよく知られており、その画像処理は、組織等の種類や診断内容に応じて多種多様である。
【０００３】
例えば、特許文献１には、超音波画像内において比較的大きなノイズの塊を抽出して除去するためにラベリング処理を利用する旨の発明が記載されている。また、特許文献２には、複数のサンプル点から得られる体積を累積加算することにより対象組織等の体積を算出する旨の発明が記載されている。そして、特許文献３には、画像データに対して収縮処理や膨張処理を行う技術が記載されている。
【０００４】
こうした状況において、本願の発明者は、超音波の画像データに対する画像処理について研究開発を重ねてきた。例えば、超音波を送受して得られるエコーデータは、超音波ビームの走査態様などに応じて、必ずしも均一な間隔で得られるとは限らない。コンベックス走査においては、深い位置におけるエコーデータほど隣接するエコーデータとの間の距離が広がる。そのため、深い位置におけるエコーデータほど比較的広い画素に対応付けられ、座標変換等の処理前においては、画像データ内において画素の大きさなどが均一ではなくなる。例えば、こうした不均一な画像データ内においても、画像処理の影響が全体的に比較的均一に及ぶことが望ましい。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００４−２６７５８４号公報
【特許文献２】特開２００７−５４３３４号公報
【特許文献３】特開２００８−１８０１６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明は、上述した研究開発の過程において成されたものであり、その目的は、超音波の画像データに対する画像処理を適正化することにある。
【課題を解決するための手段】
【０００７】
上記目的にかなう好適な超音波診断装置は、診断領域に対して超音波を送受するプローブと、プローブを制御することにより前記診断領域から受信信号を得る送受信部と、受信信号に基づいて形成される前記診断領域の画像データ内において対象物を特定する対象物特定部と、前記画像データ内において対象物を複数の対象部位に区分する区分処理部と、前記画像データ内において複数の対象部位の各々に対して各対象部位に応じた処理回数だけ段階的な画像処理を実行する画像処理部と、を有することを特徴とする。
【０００８】
上記構成では、各対象部位に応じた処理回数だけ段階的な画像処理が実行される。そのため、例えば、画像データ内において画像処理の影響が全体的に比較的均一に及ぶように各対象部位に関する処理回数を調整することなどが可能になる。
【０００９】
望ましい具体例において、前記画像データは、受信信号から得られる複数のサンプルデータを割り当てられた複数の画素により構成され、前記各対象部位に応じた処理回数は、その対象部位を構成する画素の大きさに基づいて決定される、ことを特徴とする。
【００１０】
望ましい具体例において、前記各対象部位内の代表位置における画素の大きさに基づいて、その対象部位に応じた前記処理回数を設定する処理回数設定部をさらに有する、ことを特徴とする。
【００１１】
望ましい具体例において、前記処理回数設定部は、前記画像データ内の基準位置における画素の体積と前記各対象部位内の代表位置における画素の体積とから算出される体積比に基づいて、その対象部位に応じた前記処理回数を設定する、ことを特徴とする。
【００１２】
望ましい具体例において、前記画像処理部は、複数の対象部位の各々に対して各対象部位に応じた前記処理回数だけ段階的な収縮処理を実行する、ことを特徴とする。
【００１３】
望ましい具体例において、前記画像処理部は、収縮処理された複数の対象部位の各々に対して各対象部位に応じた前記処理回数だけ段階的な膨張処理を実行する、ことを特徴とする。
【００１４】
また、上記目的にかなう好適な超音波画像処理装置は、超音波の画像データ内において対象物を特定する対象物特定部と、前記画像データ内において対象物を複数の対象部位に区分けする区分処理部と、前記画像データ内において複数の対象部位の各々に対して各対象部位に応じた処理回数だけ段階的な画像処理を実行する画像処理部と、を有することを特徴とする。
【００１５】
例えば、上述した対象物特定部と区分処理部と画像処理部の機能を実現するためのプログラムを利用することにより、コンピュータを前記超音波画像処理装置として機能させてもよい。なお、そのプログラムは、例えばディスクやメモリなどの記憶媒体に記憶され、その記憶媒体を介してコンピュータに提供されてもよい。
【発明の効果】
【００１６】
本発明により、超音波の画像データに対する画像処理が適正化される。例えば、本発明の好適な態様によれば、画像データ内において画像処理の影響が全体的に比較的均一に及ぶように画像処理を適正化することが可能になる。
【図面の簡単な説明】
【００１７】
【図１】本発明の実施において好適な超音波診断装置の全体構成を示す図である。
【図２】図１の超音波診断装置による各種の処理を説明するための図である。
【図３】三次元データ空間内におけるフィルタの走査を説明するための図である。
【図４】各部位ごとの処理回数の制御を説明するためのフローチャートである。
【発明を実施するための形態】
【００１８】
図１は、本発明の実施において好適な超音波診断装置の全体構成を示す図である。プローブ１０は、診断対象を含む領域に対して超音波を送受波する超音波プローブである。プローブ１０は、超音波を送受する複数の振動素子を備えており、複数の振動素子が送受信部１２によって送信制御されて送信ビームが形成される。また、複数の振動素子が診断対象を含む領域内から得られる超音波を受波し、これにより得られた信号が送受信部１２へ出力され、送受信部１２が受信ビームを形成して受信ビームに沿ってエコーデータが収集される。
【００１９】
プローブ１０は、超音波ビーム（送信ビームと受信ビーム）を三次元空間内において走査して立体的にエコーデータを収集する三次元プローブが好適である。例えば、一次元的に配列された複数の振動素子（１Ｄアレイ振動子）によって電子的に形成される走査面を機械的に動かすことにより超音波ビームが立体的に走査される。また、二次元的に配列された複数の振動素子（２Ｄアレイ振動子）を電子的に制御して超音波ビームを立体的に走査してもよい。もちろん、超音波ビームを断層面内において走査する二次元の超音波プローブが利用されてもよい。
【００２０】
三次元空間内において超音波ビームが走査されてエコーデータが収集されると、その三次元空間に対応した三次元データ空間を構成する複数のボクセルについてのエコーデータ（ボクセルデータ）が図示しないメモリなどに記憶される。そして、三次元データ空間を構成する複数のボクセルに対して、二値化処理部２０以降の各部において各種の処理が実行される。そこで、それら各種の処理について説明する。なお、図１に示した部分（構成）については、以下の説明において図１の符号を利用する。
【００２１】
図２は、図１の超音波診断装置による各種の処理を説明するための図である。図２においては各画像データが二次元的に描かれているものの、各種の処理は三次元データ空間内において三次元的に実行される。
【００２２】
図２（１）は、二値化処理部２０（図１）における二値化処理を示している。二値化処理部２０は三次元データ空間を構成する複数のボクセルに対して二値化処理を施すことにより、図２（１）に示す二値化処理後の画像データを形成する。本実施形態における診断対象としては、生体内の卵胞が好適である。二値化処理部２０は、所定の閾値と各ボクセルのボクセル値（エコーデータの大きさ）を比較することにより、卵胞に対応するボクセルとそれ以外のボクセルとを識別する。そして、例えば、卵胞に対応するボクセルのボクセル値を「１」とし、それ以外のボクセルのボクセル値を「０」とする。
【００２３】
なお、図２（１）の二値化画像内には、参考として、走査される超音波ビームの位置が破線で示されている。例えばコンベックス走査では、図２（１）に示すように超音波ビームが扇状に走査され、超音波ビームの深さ方向に沿って、浅い部分に比べて深い部分のビーム密度が低くなる。つまり、浅い部分に比べて、深い部分ではビーム間隔が広くなり、エコーデータの間隔が広くなる。そのため、図２（１）の二値化画像では、深い部分におけるエコーデータほど、比較的大きい画素（微小空間）が割り当てられ、画像内において画素の大きさが不均一となる。本実施形態においては、このように画素の大きさが不均一な場合においても、画像処理の影響が全体的に比較的均一に及ぶように制御される。
【００２４】
図２（２）は、ラベリング処理部３０（図１）におけるラベリング処理を示している。ラベリング処理部３０は、三次元データ空間を構成する二値化処理後のボクセルデータ内において、つまり図２（１）に示す二値化処理後の画像データ内において、ラベリング処理を施し、複数の部位に対して互いに異なるラベルを割り当て、卵胞を複数の部位に区分けする。ラベリング処理としては公知の手法を用いることができ、例えば三次元データ空間内において同じボクセル値を持つ複数のボクセルの塊が検出され、各塊ごとにラベル番号が付与される。例えば、図２（１）に示すように、ボクセル値「１」の塊である複数の部位に対して、それぞれラベル１〜３が割り当てられる。なお、ボクセル値「０」の塊である背景部分に対しては、例えばラベル０が割り当てられる。
【００２５】
ラベリング処理が施されると、収縮膨張処理部５０（図１）は、三次元データ空間を構成するラベリング処理後のボクセルデータ内において、つまり図２（２）に示すラベリング処理後の画像データ内において、複数の部位に対して収縮膨張処理（収縮処理と膨張処理）を施す。その収縮膨張処理においては、複数の部位の各々に対して各部位に応じた処理回数だけ段階的な収縮処理が実行され、さらに収縮処理が終了してから、複数の部位の各々に対して収縮処理と同じ回数だけ膨張処理が実行される。これにより、図２（３）に示すように、卵胞の表面などに存在していたノイズが除去され、さらに図２（２）においてラベル３が付されたノイズも除去される。本実施形態においては、収縮膨張処理において各部位に応じた処理回数を利用しているため、比較的浅い位置にあるラベル１の部位と比較的深い位置にあるラベル２の部位が、共に、図２（２）のものと同程度の形状かつ同程度の大きさに復元される。
【００２６】
これに対し、図２（Ｃ）に示す比較例は、図２（２）に示す画像データ内の全域において同一の処理回数を利用して収縮膨張処理を施したものである。先に説明したように、深い部分ほど、比較的大きい画素が割り当てられているため、画素を単位とする収縮膨張処理を浅い部分と深い部分で同じ処理回数だけ実行すると、深い部分ほど処理の影響が強くなる。そのため、図２（Ｃ）に示す比較例では、比較的深い位置にあるラベル２の部位に対して、収縮膨張処理の影響が必要以上に強く及んでしまい、元の形状や大きさが復元できていない。
【００２７】
本実施形態では、各部位に応じた処理回数を利用している。その処理回数は、図１の処理回数設定部４０において設定される。処理回数設定部４０は、画像データ内の基準位置における画素の体積Ｖｓと各部位内の代表位置における画素の体積Ｖｇに基づいて、その部位の処理回数を決定する。
【００２８】
画像データを構成する複数の画素（ボクセル）の体積は、体積値テーブル４２に記憶されている。例えば、超音波ビーム方向やビーム間隔やエコーデータのサンプリング間隔などの情報に基づいて、三次元の画像データを構成する複数の画素の体積が予め算出され、各画素の位置（座標）とその画素の体積値が対応付けられて、体積値テーブル４２に格納される。
【００２９】
画像データ内の基準位置は、例えば画像データの中心位置とされる。画像データ内に関心領域が設定され、その関心領域の中心位置を基準位置としてもよい。基準位置をユーザが決定するようにしてもよい。
【００３０】
また、各部位内の代表位置は、代表位置決定部３２により決定される。代表位置決定部３２は、例えば、各部位ごとにその部位の重心位置の座標を算出してその部位の代表位置とする。もちろん、重心位置以外の座標を代表位置としてもよいし、代表位置をユーザが決定するようにしてもよい。
【００３１】
処理回数設定部４０は、各部位の処理回数を決定するにあたり、代表位置決定部３２において決定される各部位の代表位置における画素の体積Ｖｇと、基準位置における画素の体積Ｖｓを体積値テーブル４２から読み出し、これらの体積の比Ｖｓ／Ｖｇを算出する。そして、基準位置における基準処理回数Ｎｓから、各部位の処理回数をＮｓ×Ｖｓ／Ｖｇとし、非整数値を四捨五入または切り捨てするなどして、各部位の処理回数を整数値とする。なお、基準処理回数Ｎｓは、予め基準として装置に記憶された値を利用してもよいしユーザが適宜に修正または設定するようにしてもよい。こうして、処理回数設定部４０は複数の部位の各々について、各部位ごとに処理回数を算出し、各部位のラベル値（図２参照）とその部位の処理回数とを対応付ける。
【００３２】
そして、収縮膨張処理部５０は、各部位ごとに設定された処理回数に応じて、収縮処理と膨張処理を実行する。その収縮処理と膨張処理においては、処理用のフィルタが利用され、そのフィルタが三次元データ空間内の全域に亘って走査される。
【００３３】
図３は、三次元データ空間１００内におけるフィルタ１２０の走査を説明するための図である。図３において、三次元データ空間１００は、ｒθφ座標系で示されている。ｒ方向は超音波ビームの深さ方向であり、その超音波ビームがθ方向とφ方向に二次元的に走査され、三次元データ空間１００を構成する複数のボクセル（画素）のエコーデータが三次元的に収集される。
【００３４】
フィルタ１２０は、ｒ方向とθ方向とφ方向のそれぞれの長さが３個のボクセルに相当し、合計２７個のボクセルに相当する立体的な構成となっている。フィルタ１２０の中心に位置するボクセルが注目ボクセルであり、それを取り囲む２６個のボクセルが周辺ボクセルである。そして、三次元データ空間１００内の全てのボクセルが注目ボクセルとなるように、フィルタ１２０がｒ方向とθ方向とφ方向に移動されて三次元データ空間１００内の全域に亘って走査される。
【００３５】
なお、走査の過程において、フィルタ１２０は、常に２７個のボクセルを対象とする。したがって、ｒ方向に沿って深くなるにつれてボクセルの大きさが増大すると、フィルタ１２０の大きさもそれに応じて大きくなる。
【００３６】
収縮処理においては、各走査位置でフィルタ１２０内の２６個の周辺ボクセルの中に１つでもボクセル値「０」のボクセルがあれば、フィルタ１２０の中心に位置する注目ボクセルのボクセル値が「０」とされる。例えば、注目ボクセルが「１（卵胞）」であり、周辺ボクセルの中に１つでもボクセル値「０（背景）」のボクセルがあれば、その注目ボクセルが「０（背景）」に変換される。
【００３７】
そして、フィルタ１２０が三次元データ空間１００内の全域に亘って一通り走査され、各走査位置においてフィルタ処理が施されることにより、１段階の収縮処理が終了する。なお、注目ボクセルに関するボクセル値の変換は、フィルタ１２０が三次元データ空間１００内の全域に亘って一通り走査された後に実行される。つまり、フィルタ１２０が走査されている途中においては、ボクセル値の変換は行われず、どの走査位置においても変換前のボクセル値を対象としてフィルタ処理が実行される。
【００３８】
１段階の収縮処理が終了し、その結果に基づいてボクセル値の変換が行われると、その変換後のボクセル値で構成される三次元データ空間１００に対して、２段階目の収縮処理が実行される。２段階目の収縮処理においても、１段階目と同じフィルタ処理が実行される。つまり、各走査位置でフィルタ１２０内の２６個の周辺ボクセルの中に１つでもボクセル値「０」のボクセルがあれば、フィルタ１２０の中心に位置する注目ボクセルのボクセル値が「０」とされ、フィルタ１２０が三次元データ空間１００内の全域に亘って一通り走査された後にボクセル値が変換される。
【００３９】
こうして、収縮処理が繰り返し実行される際に、本実施形態では、三次元データ空間１００の全域に亘る各段階の走査において、各部位（各ラベル）ごとに、収縮処理を行うか否かが適宜制御され、これにより各部位に対する処理回数が制御される。
【００４０】
一方、膨張処理においては、各走査位置でフィルタ１２０内の２６個の周辺ボクセルの中に１つでもボクセル値「１」のボクセルがあれば、フィルタ１２０の中心に位置する注目ボクセルのボクセル値が「１」とされる。例えば、注目ボクセルが「０（背景）」であり、周辺ボクセルの中に１つでもボクセル値「１（卵胞）」のボクセルがあれば、その注目ボクセルが「１（卵胞）」に変換される。
【００４１】
そして、膨張処理においても、フィルタ１２０が三次元データ空間１００内の全域に亘って一通り走査され、各走査位置においてフィルタ処理が施され、フィルタ１２０が三次元データ空間１００内の全域に亘って一通り走査された後にボクセル値が変換される。また、膨張処理においても、三次元データ空間１００の全域に亘る各段階の走査において、各部位（各ラベル）ごとに、膨張処理を行うか否かが適宜制御され、これにより各部位に対する処理回数が制御される。
【００４２】
ちなみに、二次元の画像データの場合には、図３に示すフィルタ１２０に代えて、縦と横の長さがそれぞれ３個のボクセルに相当する合計９個のボクセルに対応した二次元フィルタを利用し、中心に位置するボクセルを注目ボクセルとし、それを取り囲む８個のボクセルを周辺ボクセルとすればよい。
【００４３】
図４は、各部位ごとの処理回数の制御を説明するためのフローチャートである。図４に示すフローチャートは、収縮処理と膨張処理の両方に適用される。つまり、図４のフローチャートに従って収縮処理が実行された後、収縮処理後のボクセルを対象として、図４のフローチャートに従って膨張処理が実行される。以下に、図４のフローチャートに沿って処理回数の制御を説明する。
【００４４】
まず、現在の処理回数を示す現回数Ｎｃが１に初期化され（Ｓ４０１）、さらに、三次元データ空間１００（図３）内においてフィルタ１２０（図３）が走査の初期位置に配置され、三次元データ空間１００内の先頭ボクセルが注目ボクセルとされる（Ｓ４０２）。
【００４５】
次に、現回数Ｎｃと最小処理回数Ｎminが比較される（Ｓ４０３）。ここで、最小処理回数Ｎminは、処理回数設定部４０（図１）において設定された複数の部位の各々についての処理回数の中で最小のものである。
【００４６】
そして、現回数Ｎｃが最小処理回数Ｎmin以下の場合には、収縮または膨張のフィルタ処理が実行される（Ｓ４０４）。なお、収縮処理の場合には、本フローチャートが終了するまでＳ４０４において常に収縮処理が実行され、膨張処理の場合には、本フローチャートが終了するまでＳ４０４において常に膨張処理が実行される。
【００４７】
次に、現在のフィルタ位置における注目ボクセルが三次元データ空間１００（図３）内の最終ボクセルか否かが確認され（Ｓ４０５）、最終ボクセルでなければ、フィルタ１２０（図３）が１ボクセルだけ移動され、注目ボクセルが１つ移動される（Ｓ４０６）。
【００４８】
こうして、フィルタ１２０が三次元データ空間１００の全域に亘って走査されるまで、Ｓ４０４からＳ４０６までの処理が繰り返し実行され、注目ボクセルが三次元データ空間１００内の最終ボクセルに達すると（Ｓ４０５）、現回数Ｎｃと最大処理回数Ｎmaxが比較される（Ｓ４１３）。ここで、最大処理回数Ｎmaxは、処理回数設定部４０（図１）において設定された複数の部位の各々についての処理回数の中で最大のものである。
【００４９】
以上により、三次元データ空間１００の全域に亘ってフィルタ１２０が走査されると、１段階目の処理が終了し、現回数Ｎｃに１が加えられ（Ｓ４１４）、Ｓ４０２に戻って次の段階の処理が実行される。
【００５０】
現回数Ｎｃが最小処理回数Ｎmin以下の場合には、Ｓ４０３からＳ４０４へ進む処理が実行される。つまり、現回数Ｎｃが最小処理回数Ｎmin以下の場合には、三次元データ空間１００の全域がフィルタ処理（収縮処理または膨張処理）の対象とされ、各部位ごとの制御は行わずに処理を簡略化して演算量などを減らしている。
【００５１】
三次元データ空間１００の全域に亘るフィルタ１２０の走査が繰り返され、現回数Ｎｃが最小処理回数Ｎminより大きくなると、Ｓ４０３からＳ４０７へ進み、各部位ごとの処理回数の制御が行われる。
【００５２】
まず、注目ボクセルに対応した処理回数Ｎｔが確認される（Ｓ４０７）。例えば、注目ボクセルのラベルが参照され、そのラベルの部位について、処理回数設定部４０において設定された処理回数Ｎｔが確認される。そして、現回数Ｎｃと処理回数Ｎｔが比較される（Ｓ４０８）。
【００５３】
現回数Ｎｃが処理回数Ｎｔ以下の場合には、収縮または膨張のフィルタ処理が実行される（Ｓ４１０）。なお、収縮処理の場合には、本フローチャートが終了するまでＳ４１０において常に収縮処理が実行され、膨張処理の場合には、本フローチャートが終了するまでＳ４１０において常に膨張処理が実行される。
【００５４】
一方、現回数Ｎｃが処理回数Ｎｔよりも大きい場合には、注目ボクセルが属する部位に対するフィルタ処理は既に終了していると判断され、フィルタ処理は実行せずに注目ボクセルのボクセル値がそのまま維持される（Ｓ４０９）。
【００５５】
次に、現在のフィルタ位置における注目ボクセルが三次元データ空間１００（図３）内の最終ボクセルか否かが確認され（Ｓ４１１）、最終ボクセルでなければ、フィルタ１２０（図３）が１ボクセルだけ移動され、注目ボクセルが１つ移動される（Ｓ４１２）。
【００５６】
こうして、フィルタ１２０が三次元データ空間１００の全域に亘って走査されるまで、Ｓ４０７からＳ４１２までの処理が繰り返し実行され、注目ボクセルが三次元データ空間１００内の最終ボクセルに達すると（Ｓ４１１）、現回数Ｎｃと最大処理回数Ｎmaxが比較される（Ｓ４１３）。
【００５７】
以上の処理は、現回数Ｎｃが最大処理回数Ｎmaxを超えるまで繰り返し実行される。最大処理回数Ｎmaxは、複数の部位の処理回数の中で最大のものであるため、現回数Ｎｃが最大処理回数Ｎmaxを超えると、全ての部位に対する処理が終了したと判断され、本フローチャートが終了する。最大処理回数Ｎmaxを利用することにより、比較的容易に処理の終了を判定できる。
【００５８】
図１に戻り、体積演算部６０は、収縮膨張処理により復元された卵胞について、各部位ごとに体積を算出する。体積演算部６０は、体積値テーブル４２に格納された各部位に属するボクセル（画素）の体積を参照し、各部位ごとに体積を算出する。もちろん、複数の部位に関する体積の合計値（卵胞全体の体積）や、複数の部位に関する体積についての統計的な数値などが算出されてもよい。
【００５９】
体積演算部６０により算出された体積値等は、表示部８０に表示される。なお、二値化処理前の複数のボクセルのボクセル値に基づいて、表示画像形成部７０が卵胞を立体的に映し出した三次元の超音波画像を形成し、その超音波画像が表示部８０に表示されてもよい。三次元の超音波画像としては、例えば、ボリュームレンダリングによる画像などが好適である。
【００６０】
以上、本発明の好適な実施形態である超音波診断装置について説明したが、例えば、図１に示した二値化処理部２０とラベリング処理部３０と代表位置決定部３２と処理回数設定部４０と体積値テーブル４２と収縮膨張処理部５０と体積演算部６０のうちの少なくとも一つをコンピュータにより実現し、そのコンピュータを超音波画像処理装置として機能させてもよい。
【００６１】
なお、上述した実施形態は、あらゆる点で単なる例示にすぎず、本発明の範囲を限定するものではない。本発明は、その本質を逸脱しない範囲で各種の変形形態を包含する。
【符号の説明】
【００６２】
１０  プローブ、２０  二値化処理部、３０  ラベリング処理部、３２  代表位置決定部、４０  処理回数設定部、５０  収縮膨張処理部。
【図１】
0005467982.tif 000002
【図３】
0005467982.tif 000003
【図４】
0005467982.tif 000004
【図２】
0005467982.tif 000005