JP2020018694A - 超音波診断装置及び超音波画像処理方法 - Google Patents

Abstract

【課題】断層画像に対して関心領域を自動的に設定する。【解決手段】断層画像５０には、乳腺像５６、大胸筋像５８及び境界像６０が含まれる。断層画像５０に対して複数の探索経路６８が設定される。個々の探索経路６８において、深いところから浅いところへ向かって境界探索が実行される。複数の境界点７４に基づいて、乳腺像５６を含むように関心領域８２が設定される。関心領域内の画像部分が画像解析の対象となる。【選択図】図４

Description

本発明は、超音波診断装置及び超音波画像処理方法に関し、特に、関心領域の生成に関する。

超音波診断装置は、生体への超音波の送受波により得られた受信信号に基づいて超音波画像を形成する装置である。典型的な超音波画像は断層画像であり、それは組織の断面を表した画像である。例えば、乳房の検査においては、超音波プローブが乳房表面に当接され、それによって表示された断層画像が観察され、また、その観察を通じて乳腺中の腫瘤の有無や腫瘤の態様等が診断される。

近時、コンピュータ支援診断（Computer-Aided Diagnosis：ＣＡＤ）機能を搭載した超音波診断装置や超音波画像処理装置が普及しつつある。そのような装置では、超音波画像の評価又は診断に際してＣＡＤ機能が利用される。例えば、乳腺診断においては、ＣＡＤ機能を利用して断層画像がリアルタイムで解析される。具体的には、断層画像中に含まれる低輝度の腫瘤像（あるいは低輝度の非腫瘤）が自動的に認識され、それがマーキングされる。ＣＡＤ機能として、腫瘤像ごとに悪性度が自動的に判定されるものもある。

例えば、特許文献１には、組織の輪郭に応じて関心領域を自動的に設定する超音波診断装置が開示されている。また、特許文献２には、断層画像上に設定された複数の探索経路上において深いところから浅いところへ境界探索を行う技術が開示されている。

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

乳房の超音波診断により得られた断層画像上においては、断層画像内の深いところに、具体的には乳腺像の下側に、大胸筋像が現れ、その中の低輝度部分が腫瘤像として誤認、誤検出され易い。通常、腫瘤は乳腺の中に生じるため乳腺像が専らの観察対象となる。このため、浅いところにある乳腺像を含み且つ深いところにある大胸筋像等が除外されるように関心領域を設定することが考えられるが、マニュアルでの関心領域の設定は煩雑である。プローブの送受波面の一部（例えば端部）が体表面から浮くことに起因して、断層画像内にシャドー（黒い領域）が現れることもある。そのようなシャドーが関心領域から除外されるように関心領域を定めることが望まれるが、マニュアルでの関心領域の設定はやはり煩雑である。

一方、乳房の超音波診断により得られた断層画像には、乳腺像と大胸筋像との間に横方向に流れる境界像が比較的に明瞭に現れる。その断層画像において画像解析を行いたい部分は、通常、境界像の上側に現れている乳腺像である。一方、境界像の下側の領域はおよそ一様性をもった比較的に低輝度の領域である。それらの性質又は特徴を利用して関心領域を自動的に生成することが望まれる。

なお、乳房以外の組織の超音波診断により得られた超音波画像の中にも、上記の性質又は特徴が認められる画像があり、そのような画像についても関心領域を自動的に生成することが望まれる。関心領域の設定に関して、特許文献１には、横方向に流れる境界像を利用することは開示されていない。また、特許文献２の技術では、複数の探索経路を設定する以前に関心領域を設定する必要がある。

本発明の目的は、超音波画像の内容に適した関心領域の自動的な生成を実現することにある。

本発明に係る超音波診断装置は、深さ方向に対して交差する方向に伸びた形態を有する境界像を含む超音波画像上において、前記境界像における複数の位置を横切るように複数の探索経路を設定し、前記複数の探索経路上において深いところから浅いところへ境界探索を行うことにより複数の境界点を検出する検出部と、前記複数の境界点に基づいて、前記境界像の浅い側に存在する注目組織像を包含する関心領域を生成する生成部と、を含むことを特徴とする。

本発明に係る超音波画像処理方法は、超音波画像上において、深さ方向に対して交差する方向に伸びた形態を有する境界像に基づいて、前記境界像の浅い側に存在する注目組織像を包含する関心領域を生成する生成工程と、前記超音波画像又は他の超音波画像における前記関心領域内の画像部分に対して行われた解析の結果を表す画像とともに前記関心領域を表す画像を表示する表示工程と、を含み、前記生成工程では、前記超音波画像における深いところから浅いところへ境界探索が実行され、その結果に基づいて前記関心領域が定められる、ことを特徴とする。

本発明に係るプログラムは、超音波画像処理装置において実行されるプログラムであって、深さ方向に対して交差する方向に伸びた形態を有する境界像を含む超音波画像上において、前記境界像における複数の位置を横切るように複数の探索経路を設定し、前記複数の探索経路上において深いところから浅いところへ境界探索を行うことにより複数の境界点を検出する機能と、前記複数の境界点に基づいて、前記境界像の浅い側に存在する注目組織像を包含する関心領域を生成する機能と、を含むことを特徴とする。

本発明によれば、超音波画像の内容に適した関心領域の自動的な生成を実現できる。

実施形態に係る超音波診断装置を示すブロック図である。 ＲＯＩ設定部の構成例を示すブロック図である。 比較例を示す図である。 実施形態に係るＲＯＩ生成方法を示す図である。 境界検出方法を示す図である。 ＲＯＩの生成及び適用を示す図である。 平滑化の前後を示す図である。 時間軸上の平滑化を説明するための図である。 表示例を示す図である。 超音波画像処理方法の第１例を示す図である。 超音波画像処理方法の第２例を示す図である。 シャドー及び肋骨像を含む断層画像を示す図である。 超音波画像処理方法の第１例の処理結果を示す図である。 超音波画像処理方法の第２例の処理結果を示す図である。 境界判定方法を説明するための図である。

以下、本発明の実施形態を図面に基づいて説明する。

（１）実施形態の概要
実施形態に係る超音波診断装置は、深さ方向に対して交差する方向に伸びた形態を有する境界像を含む超音波画像を処理する機能を有し、その処理のために、検出部、及び、生成部を有する。検出部は、境界像における複数の位置を横切るように複数の探索経路を設定し、複数の探索経路上において深いところから浅いところへ境界探索を行うことにより複数の境界点を検出するものである。生成部は、複数の境界点に基づいて、境界像の浅い側に存在する注目組織像を包含する関心領域を生成するものである。

上記構成は、超音波画像において比較的に明瞭に境界像が現れていること、及び、境界像の深い側に境界像とは異なる領域が広がっていることを前提とし、それらの性質又は特徴を利用して、境界像の手前側に存在する注目組織像を包含する関心領域を設定するものである。境界像を含むように関心領域が設定されてもよいし、境界像を除外するように関心領域が設定されてもよい。超音波画像の一部にシャドーが生じた場合、シャドー内においては境界像も部分的に消失する。その場合、境界探索がプローブ近傍まで進行することになるので、シャドーを避けて関心領域が設定される。

実施形態において、注目組織像は乳腺像であり、境界像は乳腺像とそれよりも深いところにある大胸筋像との間の境界像である。実施形態において、関心領域は、超音波画像又は他の超音波画像において解析対象となる部分を画定するものである。他の目的で関心領域が利用されてもよい。

実施形態において、複数の境界点からなる境界点列に基づいて関心領域の下辺が定められ、下辺を空間的に平滑化する平滑化部が設けられる。この構成によれば、下辺が複雑な形態を有する場合にそれを滑らかにできる。また、複数の境界点の検出において誤認があってもそれによる影響を緩和できる。更に、関心領域を表示する場合、その見栄えを良好にできる。境界点列が平滑化され、その結果として、下辺が平滑化されてもよいし、下辺それ自体が平滑化されてもよい。

実施形態において、平滑化部は、空間的に平滑化された下辺を更に時間的に平滑化するものである。この構成によれば時間軸方向に関心領域の下辺の形態が激しく変化するような場合においてその変化を抑えることができる。結果として関心領域の設定精度を向上できる。

実施形態において、検出部は、各探索経路上において輝度条件が充足されたことに基づいて各境界点を検出する。境界検出は、例えば、下方から上方へエコーデータ単位で実行される。その場合において、１つのエコーデータが参照されてもよいし、一次元窓内の複数のエコーデータが参照されてもよいし、二次元窓内の複数のエコーデータが参照されてもよい。

実施形態において、検出部は、各探索経路上において輝度条件が連続して充足されたことに基づいて各境界点を検出する。この構成は境界像がある程度の厚みを有することを前提とするものであり、この構成によれば、境界誤認を防止又は軽減できる。例えば、探索経路上において輝度条件がｎ回連続して充足されたことに基づいて境界点が検出されてもよい。その場合、ｎが診断レンジ（診断深さ範囲）に基づいて適応的に設定されてもよい（但しｎは１以上の整数である）。診断レンジに応じて境界像の厚みが変化するので、そのような変化に合わせてｎの値を適応的に設定するものである。

実施形態においては、超音波画像とともに関心領域を示す関心領域画像を表示する表示器を含む。この構成によれば、関心領域に注目組織像が含まれていることを確認することが可能となる。超音波画像に対して関心領域画像が重畳表示されてもよいし、それらの画像が並列表示されてもよい。

実施形態に係る超音波画像処理方法は、生成工程及び表示工程を含むものである。生成工程では、超音波画像上において、深さ方向に対して交差する方向に伸びた形態を有する境界像に基づいて、境界像の浅い側に存在する注目組織像を包含する関心領域が生成される。具体的には、生成工程では、超音波画像上における深いところから浅いところへの境界探索が実行され、その結果に基づいて関心領域が定められる。表示工程では、超音波画像又は他の超音波画像における関心領域内の画像部分に対して行われた解析の結果を表す画像とともに関心領域を表す画像が表示される。

実施形態に係るプログラムは、超音波画像処理装置において実行されるプログラムである。ここで、超音波画像処理装置は、超音波診断装置、情報処理装置等を含む概念である。プログラムは、可搬型記憶媒体を介して、又は、ネットワークを介して超音波画像処理装置にインストールされる。

（２）実施形態の詳細
図１には、実施形態に係る超音波診断装置の構成がブロック図として示されている。超音波診断装置は、病院等の医療機関に設置され、生体（被検者）に対する超音波の送受波により得られた受信データに基づいて超音波画像を形成する医療用の装置である。実施形態に係る超音波診断装置は、後に詳述するように、関心領域（ＲＯＩ）を自動的に設定する機能、及び、超音波画像を自動的に解析する機能（ＣＡＤ機能）を備えている。実施形態において、超音波診断対象となる組織は乳房であり、より詳しくは乳腺である。もちろん、他の組織が超音波診断対象となってもよい。

超音波プローブ１０は、超音波を送受波する手段又は送受波器として機能するものである。乳腺の超音波診断に際しては、超音波プローブの１０の送受波面（音響レンズ表面）が被検者の乳房表面に当接され、その状態で超音波が送受波される。超音波プローブ１０は、一次元配列された複数の振動素子からなる振動素子アレイを備えている。振動素子アレイによって超音波ビーム１１が形成され、超音波ビーム１１の電子的な走査により走査面が形成される。走査面は観察面であり、すなわち二次元データ取込領域である。超音波ビームの電子走査方式として、電子セクタ走査方式、電子リニア走査方式等が知られている。超音波ビームのコンベックス走査が行われてもよい。超音波プローブ内に２Ｄ振動素子アレイを設け、生体内からボリュームデータが取得されてもよい。

送信部１２は、送信時において、複数の振動素子に対して複数の送信信号を並列的に供給する送信ビームフォーマーであり、それは電子回路として構成される。受信部１４は、受信時において、複数の振動素子から並列的に出力される複数の受信信号を整相加算（遅延加算）する受信ビームフォーマーであり、それは電子回路として構成される。受信部１４は、複数のＡ／Ｄ変換器、検波回路等を備えている。受信部１４での複数の受信信号の整相加算によりビームデータが生成される。ちなみに、１回の電子走査当たり、電子走査方向に並ぶ複数のビームデータが生成され、それらが受信フレームデータを構成する。個々のビームデータは深さ方向に並ぶ複数のエコーデータにより構成される。

ビームデータ処理部１６は、受信部１４から出力された各ビームデータを処理する電子回路である。その処理には、対数変換、相関処理、等が含まれる。処理後の各ビームデータが断層画像形成部１８へ送られている。

断層画像形成部１８は受信フレームデータに基づいて断層画像（Ｂモード断層画像）を形成する電子回路である。それはＤＳＣ（Digital Scan Converter）を有している。ＤＳＣは、座標変換機能、補間機能、フレームレート変換機能等を有し、ビーム走査方向に並ぶ複数のビームデータに基づいて断層画像を形成する。断層画像のデータが表示処理部２０及びＲＯＩ設定部２２へ送られている。

なお、表示処理部２０、ＲＯＩ設定部２２及び画像解析部２４によって画像処理モジュール２６が構成されている。画像処理モジュール２６は、プログラムに従って動作する１又は複数のプロセッサにより構成され得る。画像処理モジュール２６がＰＣ等の情報処理装置上に構築されてもよい。その場合、断層画像データが超音波診断装置から情報処理装置へ転送される。情報処理装置は超音波画像処理装置として機能する。

ＲＯＩ設定部２２は、乳腺像を含む断層画像に基づいてＲＯＩ（つまり関心領域）を自動的に生成するものである。実施形態においては、フレーム単位でＲＯＩが生成されている。もちろん、所定数のフレームを単位としてＲＯＩが生成されてもよい。ＲＯＩ設定部２２は、生成されたＲＯＩを表す座標情報を画像解析部２４に送っている。また、ＲＯＩ設定部２２は、枠としてのＲＯＩ及びその内部画像を含むＲＯＩ画像を生成し、そのデータを表示処理部に送っている。ＲＯＩ設定部２２の詳細については後に詳述する。

画像解析部２４は、画像解析手段として機能し、断層画像の内で関心領域に含まれる画像部分に対して画像解析を実行するものである。すなわち、画像解析部２４はＣＡＤ機能を発揮するものである。画像解析部２４は、フレーム単位で画像解析を行う。もちろん、所定数のフレームを単位として画像解析が実行されてもよい。画像解析部２４は、ＣＮＮ（Convolutional Neural Network）等の機械学習型解析器によって構成され得る。画像解析部２４は、低輝度の腫瘤や低輝度の非腫瘤等を認識、抽出又は弁別する機能を有している。画像解析部２４が腫瘤の悪性度を評価する機能を備えていてもよい。画像解析部２４は、実施形態において、断層画像を解析して腫瘤等を特定し、それを指し示すマーカーを生成する。マーカーが付与された断層画像（つまり解析結果画像）のデータが表示処理部２０へ送られている。

表示処理部２０は、グラフィック画像生成機能、カラー演算機能、画像合成機能等を有する。具体的には、表示処理部２０は、解析結果画像とＲＯＩ画像とを並べて合成し、これにより表示画像を生成する機能を有している。解析結果画像は、上記のように画像解析において見付かった部位に対してマーキングが施されている画像である。ＲＯＩ画像はＲＯＩの形態及び中身を表す画像である。解析結果画像又は断層画像に対して、ＲＯＩ画像が重畳されてもよい。表示器２８には断層画像等の表示画像が表示される。それはＬＣＤ、有機ＥＬ表示デバイス等によって構成される。

画像解析部２４において、弾性情報画像（エラスト画像）２５が解析されてもよい。弾性情報画像は、組織変位から計算される組織弾性情報を表した画像である。その場合においても、弾性情報画像上にＲＯＩが設定され、ＲＯＩ内部が解析対象となる。他の超音波画像が解析対象となってもよい。なお、実施形態においては、座標変換後の断層画像が解析対象とされているが、座標変換前のフレームデータが解析対象とされてもよい。そのフレームデータも組織断面を表すデータであるから、広義の断層画像に含まれ得るものである。

制御部３４は、図１に示されている各構成要素の動作を制御するものである。制御部３４は、実施形態において、ＣＰＵ及びプログラムによって構成される。制御部３４が上記の画像処理モジュール２６として機能してもよい。操作パネル３２は入力デバイスであり、それは複数のスイッチ、複数のボタン、トラックボール、キーボード等を有する。なお、図１において、断層画像形成部１８以外の超音波画像形成部については図示省略されている。例えば、弾性情報（エラスト）画像形成部、血流画像形成部、その他が設けられてもよい。

図２には、ＲＯＩ設定部２２の構成例が図示されている。ＲＯＩ設定部２２は、図示の例において、前処理部３６、境界検出部３８及びＲＯＩ生成部４２を有している。前処理部３６は、断層画像に対して必要な前処理を施すものである。前処理として、平滑化処理、最小値抽出処理、最大値抽出処理、メディアン（中央値抽出）処理、エッジ強調処理等があげられる。後に行う境界検出のために断層画像の外側に０を並べるゼロパディングが実行されてもよい。

境界検出部３８は、後に説明するように、断層画像に対して複数の探索経路を設定し、各探索経路上において境界探索を実行するものである。境界探索の開始点は、実施形態において、各探索経路上において最も深い地点とされており、その開始点から浅い方へ境界探索が順番に進められる。乳房の断層画像においては、通常、乳腺像と大胸筋像との間に境界像が明瞭に現れ、且つ、境界像の奥側（深い側）が概ね一様性をもった低輝度の領域となる。それらの性質又は特徴を前提として、深いところから浅いところへ、境界探索が順次実行される。ちなみに、実施形態において観察対象は乳腺像であり、それは境界像の手前側つまり浅い側に存在する。境界像を基準として、乳腺像をカバーするようにＲＯＩが自動的に設定される。これについては後に更に詳しく説明する。

境界平滑化部４０は、複数の境界点からなる境界点列又はそれに基づくＲＯＩ下辺に対して平滑化を行うものである。この平滑化により境界点列又は下辺が滑らかとなるので、ＲＯＩを表示する場合に見栄えがよくなる。また、この平滑化によりＲＯＩの形態が不自然なものとなることを防止又は軽減できる。その平滑化は空間的な平滑化である。実施形態に係る境界平滑化部４０は、更に時間軸方向に境界点列又は下辺を平滑化する機能も備えている。時間的な平滑化により、ＲＯＩの形態がフレーム単位で激しく変化することを抑圧することが可能となり、目障り感を緩和することができる。もっとも、超音波プローブを動かしている最中においては、時間的な平滑化をオフにしてもよい。

ＲＯＩ生成部４２は、ＲＯＩ生成手段として機能し、それは境界点列に基づいてＲＯＩを定めるものである。具体的には、境界点列に基づいてＲＯＩの下辺が定められる。ＲＯＩの２つの側辺（左辺及び右辺）を、断層画像の２つの側辺に揃えてもよい。平滑化によって生じた線で側辺を構成してもよい。ＲＯＩの上辺は断層画像の上辺に一致している。但し、断層画像の上辺から離れた位置にＲＯＩの上辺を設定するようにしてもよい。いずれにしても、画像解析対象となる組織ができるだけカバーされるようにＲＯＩが生成される。矩形を基本形としてＲＯＩが生成されているが、他の形態をもったＲＯＩが生成されてもよい。画像解析部２４は、生成されたＲＯＩに基づき、断層画像等から部分画像を切り出し、それに対して解析を実行する。あるいは、画像解析部２４は、生成されたＲＯＩに基づき、断層画像等において解析対象範囲を画定し、その範囲内において解析を実行する。

図３には、比較例が示されている。この比較例ではＲＯＩが設定されておらず、断層画像５０の全体が解析対象とされている。断層画像５０には、表層像５２、脂肪層像５４、乳腺像５６、及び、大胸筋像５８が含まれる。乳腺像５６と大胸筋像５８の間には境界像６０が現れている。境界像６０は概ね明瞭に現れる像である。このような断層画像５０に対して画像解析を実行すると、図示されるようなマーカー６４，６６が自動的に表示される。マーカー６４は、低輝度の腫瘤６２を指し示している。一方、マーカー６６は、大胸筋内における低輝度部分を指し示しており、つまり誤認が生じている。大胸筋像５８には肋骨像が現れることもあり、また、アーチファクトが現れることもある。そのような要素によって誤認が引き起こされる。ユーザー（医師、検査技師等）が観察したいものは乳腺像５６に含まれる腫瘤等であり、それは境界像６０の上側に現れるものである。一方、境界像６０の下側には比較的低輝度で概ね一様性をもった領域が広がっている。乳房の断層画像には以上のような幾つかの特有性が認められる。実施形態は、そのような特有性を利用して関心領域を自動的にしかも適切に定めるものである。以下に詳述する。

図４には、実施形態に係る境界像検出方法が示されている。図４は図３に示した断層画像に相当する模式図である。xは水平方向（横方向）を示しており、それは実施形態において電子走査方向である。ｙは垂直方向（縦方向）を示しており、それは実施形態において深さ方向である。

断層画像において、境界像６０は、それ全体として、深さ方向であるｙ方向に交差する方向に伸びている形態を有する。境界像６０が流れている方向は概ねｘ方向であるが、ｘ方向に対して傾斜している場合も多い。

実施形態においては、断層画像に対して複数の探索経路６８が設定される。複数の探索経路６８は、横方向に流れる境界像６０における複数の位置を横断している。具体的には、各探索経路６８は、断層画像を構成する個々の垂直画素列に相当する。もっとも、水平方向に離間しつつ複数の探索経路６８を設定してもよい。探索経路に沿った探索と境界点からの横方向の探索とを併用してもよい。例えば、１番から最終番までの探索経路６８において順番に境界の探索が実行される。複数の探索経路をｙ方向に対して傾斜させてもよい。コンベックス走査においては、個々のビーム方向（深さ）に平行に放射状に探索経路を設定してもよいし、画像における垂直方向に平行に探索経路を設定してもよい。

各探索経路６８上においては、図示の例では、最も深い座標が探索の開始点７０とされる。符号７２で示されているように、開始点７０から浅い側へ、つまり超音波プローブに近付く方向へ、境界探索が繰り返し実行される。具体的には、各エコーデータを注目データとして、注目データごとに、輝度条件が満たされるか否かが判定される。輝度条件が満たされた場合にその注目データの座標が境界点７４として認定される。輝度条件は、後述するように、例えば、注目データを中心とした窓内における輝度の最小値が所定の閾値以上であるという条件である。各探索経路６８における境界点７４は、境界像６０上の点であると同時に、関心領域に属する部分７６と非関心領域に属する部分７８とを分かつ点である。開始点を最深点よりも浅い点としてもよい。

最終的に、境界像６０に沿って並ぶ複数の境界点７４からなる境界点列８０が構成される。境界点列８０は、ＲＯＩの下辺に相当する。下辺はＲＯＩ８０と非ＲＯＩとを分ける境界線である。画像の局所的な状態によっては、境界点を検出できないことも生じ得る。他の点を境界点であるとして誤認してしまうことも生じ得る。誤認を防止するために、輝度条件に加えて他の条件を組み合わせるのが望ましい。

図５には、境界判定方法の一例が示されている。符号５０は断層画像を示しており、その外側にはゼロ値で満たされた領域８８が存在している。ここでは、断層画像５０及び領域８８からなる画像が処理対象画像９０とされている。

図示の例では、探索経路６８上のデータ（注目データ）ａごとに、窓８６が設定される。窓８６は３×３のサイズを有する（ａ〜ｉを参照）。例えば、窓８６内における輝度の最小値が閾値未満である場合に非境界が判定され、窓８６内における輝度の最小値が閾値以上である場合に境界が判定される。輝度条件を満たした時点での注目データの座標が境界点とされる。その座標を基準にそれよりも浅い地点又はそれよりも深い地点を境界点と定めるようにしてもよい。大胸筋像の内部においては、輝度の最小値がゼロに近くなり、一方、境界像においては輝度の最小値が大きくなる。そのような違いを利用して境界の判定を行うものである。輝度の最小値に代えて、輝度の平均値等を利用してもよい。また、二次元の窓に代えて一次元の窓を利用してもよい。あるいは、注目データのみ参照するようにしてもよい。あるいは、窓内の複数のデータに基づいてヒストグラム等を作成し、それに基づいて境界判定を行ってもよい。

図６において、元画像９２は、断層画像であり、あるいは、それに対してゼロパディングを施した画像である。元画像９２に対して上記のように生成されたＲＯＩ９４が適用される。ＲＯＩ９４の内部が画像解析の対象となる。ＲＯＩ９４の外部はマスク領域であり、図６にはそれがマスク画像９６と表現されている。符号９２Ａは画像解析結果を示す画像である。その画像９２Ａにおいては乳腺像に含まれる腫瘤像がマーカー９８によってマーキングされている。大胸筋像中には低輝度部分が含まれるものの、そこに対してマーカーは与えられていない。つまり、図３で説明した誤認が防止されている。

図７には、平滑化前後のＲＯＩ画像１００，１０４が示されている。ＲＯＩ画像１００は下辺１０２を有し、それはかなり複雑な形態を有している。検出点列の形態が下辺１０２の形態にそのまま反映されている。下辺１０２に対して平滑化を施したものがＲＯＩ画像１０４の下辺１０６である。下辺１０６はそれ全体として滑らかな形態を有している。なお、画像１０４の下辺１０６は、平滑化の作用を模式的にあるいは誇張して示したものであり、実際の平滑化処理により生成されたものではない。

図８には、平滑化方法が示されている。ｘ方向は水平方向であり、その軸上には複数の探索経路上で検出された複数の境界点のｙ座標が示されている。それらの内で、注目ｘ座標（ｙ座標はym）（符号１０８を参照）を中心とした一定の区間１１０内に含まれるｙ座標（ym-k〜ym〜ym+k）が特定され、それらの平均値y’mが演算され（符号１１２を参照）、その平均値y’mが注目ｘ座標に与えられる。この処理が区間１１０を移動させながら繰り返し実行される（符号１１４を参照）。単純平均に代えて重み付け平均等を利用してもよい。更に、個々のｘ座標において、時間軸方向にｙ座標を平滑化して平均値y”m（符号１１６を参照）を演算し、それを個々のｘ座標に与えてもよい。空間的平滑化によればフレーム内においてＲＯＩの見栄えを良好にでき、時間的平滑化によればフレーム間にわたって見栄えを良好にできる。

図９には、表示例が示されている。表示画像１２０は、左右方向に並んだ解析結果画像１２２及びＲＯＩ画像１２４により構成されている。解析結果画像１２２は断層画像及びマーカーにより構成されるものである。ＲＯＩ画像１２４は断層画像から切り出された画像部分により構成される参照画像である。ＲＯＩ画像１２４を通じて画像解析対象となった部分を確認することが可能となる。解析結果画像１２２上にＲＯＩを表す図形を重畳表示してもよい。

図１０には、実施形態に係る超音波画像処理方法の第１例がフローチャートとして示されている。このフローチャートは、図１に示した超音波診断装置に含まれる画像処理モジュールの動作を示すものである。

Ｓ１０においては、初期設定として、カウンタとしてのｊに１が与えられ、カウンタとしてのｉに最大値Dmaxが与えられる。ｊはライン番号（垂直画素列番号）を示すものであり、ｉは深さ番号を示すものである。Ｓ１２は、境界探索経路としての注目ラインがｊ番目のラインとされる。Ｓ１４では、注目データがｉ番目のデータとされる。Ｓ１６では上記の輝度条件が充足されたか否かが判断される。輝度条件が充足されなかった場合、処理がＳ２０へ移行し、ｉが１つインクリメントされる。

一方、Ｓ１６において、輝度条件が充足された場合、Ｓ２２において、注目データの座標が境界点とされ、境界点よりも手前側のライン部分が関心領域の一部とされ、境界点よりも奥側のライン部分が非関心領域の一部とされる。境界点それ自体はいずれの領域とされてもよい。Ｓ２４ではｊが１つインクリメントされる。Ｓ２６では、ｊがNmaxを超えたか否かが判断される。Nmaxは最大ライン番号である。ｊがNmaxを超えていなければ、Ｓ１２からの各工程が繰り返し実行される。ｊがNmaxを超えていれば、本処理が終了する。

図１１には、実施形態に係る超音波画像処理方法の第２例がフローチャートとして示されている。図１１に示した工程と同様の工程には同一の工程番号を付し、その説明を省略する。

この第２例では、Ｓ１６の後にＳ１８が設けられている。Ｓ１８では、輝度条件の充足がｎ回連続したか否かが判断される。つまり、連続条件の充足の有無が判断される。連続条件が満たされない場合には処理がＳ２０へ移行し、連続条件が満たされた場合には処理がＳ２２へ移行する。大胸筋像の内部においては、横筋状の像が生じるものの、その垂直方向の厚みは薄く、連続条件は満たされ難くなる。一方、境界像は深さ方向にある程度の厚みをもっており、その内部においては高輝度画素が多く存在する。よって、乳腺像の中へ境界探索が進行した場合、連続条件が満たされやすくなる。第２例によれば、より正確に境界像を検出することが可能となる。上記のｎは１以上の整数である。ｎを診断レンジ（診断深さ範囲）に応じて適応的に可変設定してもよい。

図１２〜図１４を用いて上記第２例の効果を説明する。図１２には乳房の断層画像１３０が示されている。断層画像１３０には境界像１３６及び大胸筋像１３４が含まれ、大胸筋像１３４の中には肋骨像１３８が認められる。また、断層画像１３０の左端部分にはシャドー１３２が生じている。例えば、超音波プローブの送受波面の端部が体表から離れ、又は、それらの間での密着性が低下した場合、このようなシャドーが生じ易くなる。

図１３には、上記第１例の処理結果が示されている。ここでは非ＲＯＩ（マスク領域）１４０が破線で示されている。残りの部分がＲＯＩに相当する。非ＲＯＩ１４０にはシャドーの大半が取り込まれており、シャドーが解析対象となることが回避されている。しかし、境界像の端部が低輝度ながらシャドーに入り込んでおり、そこまでが境界として検出された結果、シャドーの一部が非ＲＯＩ１４０から外れて、ＲＯＩに入り込んでいる。また、ＲＯＩの下辺に相当する、非ＲＯＩ１４０の上辺１４０ａが、肋骨像１３８に沿って設定されており、誤検出が生じている。

図１４には、上記第２例の処理結果が示されている。第２例では、連続条件が付加されるため、非ＲＯＩ１４２の形態が妥当性なものとなっている。すなわち、境界像１３６の内で端部を除く部分に沿って上辺１４２ａが設定されており、また、上辺１４２ａは肋骨像１３８には影響されていない。このように第２例によればＲＯＩをより正確に設定することが可能である。

図１５には、断層画像に基づく境界判定方法が整理されている。ここでは、輝度条件１４４、空間（連続）条件１５８、及び、時間（連続）条件１６０が示されている。各条件１４４，１５８，１６０はそれぞれ単独で又は組み合わせて利用することが可能である。

輝度条件１４４の適用に際しては、窓内の最小値１４５、平均値１４６、他の代表値１４８等を参照し得る。あるいは、窓内の輝度に基づくヒストグラム１５０や他の特徴量１５２が参照されてもよい。１又は複数の参照値が１又は複数の閾値１５６と比較され、これによって輝度条件の充足の有無が判定される（符号１５４を参照）。空間条件１５８は、空間的な連続性を評価するための条件である。例えば、それまでに検出された１又は複数の検出点に対して、垂直方向に大きく離れていない場合に空間条件１５８の充足が判定される。時間条件１６０は、時間的な連続性を評価するための条件である。例えば、時間軸上に並ぶ過去に検出された１又は複数の検出点に対して、垂直方向に大きく離れていない場合に時間条件の充足が判定される。１又は複数の条件１４４，１５８，１６０が満たされた場合に境界が判定される（符号１６２を参照）。

上記実施形態においては、乳房の断層画像が処理対象となっていたが、他の断層画像又は他の超音波画像が処理対象とされてもよい。例えば、腹部臓器の断層画像の中には、横方向に流れる境界像及びその下側の低輝度領域を含むものがある。そのような画像に対して上記処理が適用されてもよい。上記実施形態においては、関心領域によって画像解析対象が画定されていたが、関心領域が他の目的で利用されてもよい。例えば、エラスト画像を作成する範囲が上記関心領域として自動的に設定されてもよい。

１０ 超音波プローブ、１８ 断層画像形成部、２０ 表示処理部、２２ ＲＯＩ設定部、２４ 画像解析部、３６ 前処理部、３８ 境界検出部、４０ 境界平滑化部、４２ ＲＯＩ生成部。

Claims (11)

1. 深さ方向に対して交差する方向に伸びた形態を有する境界像を含む超音波画像上において、前記境界像における複数の位置を横切るように複数の探索経路を設定し、前記複数の探索経路上において深いところから浅いところへ境界探索を行うことにより複数の境界点を検出する検出部と、
   前記複数の境界点に基づいて、前記境界像の浅い側に存在する注目組織像を包含する関心領域を生成する生成部と、
   を含むことを特徴とする超音波診断装置。
2. 請求項１記載の超音波診断装置において、
   前記注目組織像は乳腺像であり、
   前記境界像は前記乳腺像とそれよりも深いところにある大胸筋像との間の境界像である、
   ことを特徴とする超音波診断装置。
3. 請求項１記載の超音波診断装置において、
   前記関心領域は、前記超音波画像又は他の超音波画像において解析対象となる部分を画定するものである、
   ことを特徴とする超音波診断装置。
4. 請求項１記載の超音波診断装置において、
   前記複数の境界点からなる境界点列に基づいて前記関心領域の下辺が定められ、
   前記下辺を空間的に平滑化する平滑化部が設けられた、
   ことを特徴とする超音波診断装置。
5. 請求項４記載の超音波診断装置において、
   前記平滑化部は、前記空間的に平滑化された下辺を更に時間的に平滑化する、
   ことを特徴とする超音波診断装置。
6. 請求項１記載の超音波診断装置において、
   前記検出部は、前記各探索経路上において輝度条件が充足されたことに基づいて前記各境界点を検出する、
   ことを特徴とする超音波診断装置。
7. 請求項６記載の超音波診断装置において、
   前記検出部は、前記各探索経路上において前記輝度条件が連続して充足されたことに基づいて前記各境界点を検出する、
   ことを特徴とする超音波診断装置。
8. 請求項７記載の超音波診断装置において、
   前記検出部は、前記各探索経路上において前記輝度条件がｎ回連続して充足されたことに基づいて前記各境界点を検出し、
   前記ｎは診断レンジに基づいて適応的に設定される、
   ことを特徴とする超音波診断装置。
9. 請求項１記載の超音波診断装置において、
   前記超音波画像とともに前記関心領域を示す関心領域画像を表示する表示器を含む、
   ことを特徴とする超音波診断装置。
10. 超音波画像上において、深さ方向に対して交差する方向に伸びた形態を有する境界像に基づいて、前記境界像の浅い側に存在する注目組織像を包含する関心領域を生成する生成工程と、
    前記超音波画像又は他の超音波画像における前記関心領域内の画像部分に対して行われた解析の結果を表す画像とともに前記関心領域を表す画像を表示する表示工程と、
    を含み、
    前記生成工程では、前記超音波画像における深いところから浅いところへ境界探索が実行され、その結果に基づいて前記関心領域が定められる、
    ことを特徴とする超音波画像処理方法。
11. 超音波画像処理装置において実行されるプログラムであって、
    深さ方向に対して交差する方向に伸びた形態を有する境界像を含む超音波画像上において、前記境界像における複数の位置を横切るように複数の探索経路を設定し、前記複数の探索経路上において深いところから浅いところへ境界探索を行うことにより複数の境界点を検出する機能と、
    前記複数の境界点に基づいて、前記境界像の浅い側に存在する注目組織像を包含する関心領域を生成する機能と、
    を含むことを特徴とするプログラム。