(57)【要約】【課題】生体内の結石の位置を超音波断層像に基づいて検出する際に、超音波断層像におけるノイズが検出精度を低下させる。【解決手段】超音波断層像の画素データの原データ分布４０からレイリー分布に従う成分４２を減じて残差データ分布４４を求める。原データ分布４０の強度に対する残差データ分布４４の強度が所定の割合以上となる抽出対象データ範囲を求め、超音波断層像にて当該抽出対象データ範囲に属する画素からなる高反射領域を結石の候補領域として抽出する。【選択図】図３
2010207427.tif 000002
【特許請求の範囲】
【請求項１】
送出した超音波に対する反射波の強度に応じた画素データを有する生体内の超音波断層像に基づいて、前記生体内の結石を検出する生体内結石検出装置であって、
前記画素データの原データ分布からレイリー分布に従う成分を減じて、前記原データ分布のピーク位置での強度を抑圧し、これにより残る残差データ分布を求める残差分布抽出手段と、
前記ピーク位置よりも前記反射波の強度が大きなデータ範囲であって、前記原データ分布の強度に対する前記残差データ分布の強度が所定の割合以上となる抽出対象データ範囲を求め、当該抽出対象データ範囲に属する画素からなる高反射領域を前記結石の候補領域として抽出する候補領域抽出手段と、
を有することを特徴とする生体内結石検出装置。
【請求項２】
請求項１に記載の生体内結石検出装置において、
前記候補領域抽出手段は、前記残差データ分布と前記原データ分布から減じた前記レイリー分布とが等しい強度となる前記画素データを閾値とし、当該閾値以上を前記抽出対象データ範囲とすること、を特徴とする生体内結石検出装置。
【請求項３】
請求項１又は請求項２に記載の生体内結石検出装置において、
前記残差分布抽出手段は、前記原データ分布と前記レイリー分布とがそれぞれのピーク点と前記ピーク位置を超えるデータ範囲に存在する変曲点とで互いに一致するように前記レイリー分布を定めること、を特徴とする生体内結石検出装置。
【請求項４】
請求項１から請求項３のいずれか１つに記載の生体内結石検出装置において、
前記超音波の送受波位置から見た前記高反射領域の背後領域からの前記反射波の強度を、その隣接領域からの前記反射波の強度と比較する背後領域反射強度判定手段を有し、
前記候補領域抽出手段は、前記隣接領域よりも前記反射波の強度が小さい前記背後領域を伴う前記高反射領域を前記結石の候補領域として抽出すること、
を特徴とする生体内結石検出装置。
【請求項５】
請求項４に記載の生体内結石検出装置において、
前記背後領域反射強度判定手段は、前記高反射領域を通過する超音波ビームの走査線における前記背後領域に属する注目区間での前記画素データの平均値である背後データ強度と、前記高反射領域の両脇を通過する２本の走査線についての前記注目区間に対応する区間での前記画素データの平均値である隣接データ強度との比を求め、
前記候補領域抽出手段は、前記隣接データ強度に対する前記背後データ強度の比が所定の閾値未満となる前記高反射領域を前記候補領域として抽出すること、
を特徴とする生体内結石検出装置。
【請求項６】
請求項１から請求項５のいずれか１つに記載の生体内結石検出装置において、
検出対象とする前記結石であって予め輪郭が特定されている特定結石について当該輪郭を表すテンプレートを記憶するテンプレート記憶手段を有し、
前記候補領域抽出手段は、前記テンプレートと前記高反射領域の輪郭とが適合する前記テンプレートの配置を探索し、当該配置の前記テンプレートに対応する前記高反射領域を前記特定結石の候補領域として抽出するテンプレート照合手段を有し、
前記テンプレート記憶手段は、前記テンプレートを表す情報として、前記特定結石の輪郭に沿って所定の順序で配置される複数の代表点それぞれについて、当該代表点での前記結石の輪郭の形状に関する形状特徴情報と、前記結石の輪郭に対して所定の位置に設定される基準点と当該代表点との相対的な位置関係を表す特徴ベクトル情報とを記憶し、
前記テンプレート照合手段は、前記各代表点毎に、前記高反射領域の輪郭上にて前記形状特徴情報に適合する適合点を探索し、当該適合点が当該代表点であるとした場合に前記基準点に対応する点となる仮想基準点を前記特徴ベクトル情報に基づいて求め、前記複数の代表点に対応して複数得られる前記仮想基準点の分布から前記超音波断層像における前記基準点の位置を推定し、当該推定された基準点に基づいて前記超音波断層像での前記テンプレートの配置を決定すること、
を特徴とする生体内結石検出装置。
【請求項７】
請求項６に記載の生体内結石検出装置において、
前記複数の代表点は、前記結石の輪郭形状に応じた凸多角形の頂点である凸頂点を含み、
前記形状特徴情報は、前記各凸頂点に接続される前記凸多角形の辺の向きを含むこと、
を特徴とする生体内結石検出装置。
【請求項８】
請求項７に記載の生体内結石検出装置において、
前記凸頂点は、前記結石の輪郭上に設定される探査の始点、及び当該始点から前記結石の輪郭を所定の凸頂点探査手順で所定の回転方向に一周順次探査して得られる複数の点であり、
前記凸頂点探査手順は、前記回転方向を示すベクトルをｈＲとし、先行して得られた前記凸頂点を基準凸頂点とし、当該基準凸頂点から次の前記凸頂点の候補点へ向かうベクトルｈＣ１と前記基準凸頂点から次の探査点へ向かうベクトルｈＳ１との外積（ｈＣ１×ｈＳ１）を前記ベクトルｈＲのｋ１倍（ここでｋ１は実数である）と置いたときに、前記ｋ１がε１＞０なる所定の閾値ε１以下である場合には、前記候補点を前記探査点で置き換えて探査を継続し、前記ｋ１が閾値ε１より大きい場合には、現在の前記候補点を前記基準凸頂点の次の前記凸頂点と定めるものであること、
を特徴とする生体内結石検出装置。
【請求項９】
請求項７又は請求項８に記載の生体内結石検出装置において、
前記複数の代表点はさらに、前記凸頂点で前記結石の輪郭を区切って得られる各輪郭区間における凹頂点を含み、
前記形状特徴情報は、前記各輪郭区間における前記凸頂点及び前記凹頂点を順次接続する辺の向きを含むこと、
を特徴とする生体内結石検出装置。
【請求項１０】
請求項９に記載の生体内結石検出装置において、
前記凹頂点は、前記結石の輪郭を周回する所定の回転方向を探査順序の方向として、前記輪郭区間をその一方端から所定の凹頂点探査手順で順次探査して得られる点であり、
前記凹頂点探査手順は、当該探査の始点となる前記凸頂点又は先行して得られた前記凹頂点を基準頂点とし、当該基準頂点から次の前記凹頂点の候補点へ向かうベクトルｈＣ２と前記基準頂点から次の探査点へ向かうベクトルｈＳ２との外積（ｈＣ２×ｈＳ２）を前記ベクトルｈＲのｋ２倍（ここでｋ２は実数である）と置いたときに、前記ｋ２がε２＜０なる所定の閾値ε２以上である場合には、前記候補点を前記探査点で置き換えて探査を継続し、前記ｋ２が閾値ε２未満である場合には、現在の前記候補点を前記凹頂点と定めるものであること、
を特徴とする生体内結石検出装置。
【請求項１１】
請求項６から請求項１０のいずれか１つに記載の生体内結石検出装置において、
前記テンプレートの配置に基づいて把握される、過去の時系列での前記結石の位置からなる結石軌跡情報を記憶する結石軌跡記憶手段と、
前記結石の位置に関する前記生体の呼吸に応じた周期的な時間変化を表す運動モデルを記憶する運動モデル記憶手段と、
現時刻における前記結石軌跡情報を前記運動モデルと照合し、前記運動モデルの周期における現時刻の位相を検出し、所定時間後の前記結石の予測位置を求める結石位置予測手段と、
を有することを特徴とする生体内結石検出装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、生体を超音波ビームで走査して得られる超音波断層像に基づいて、生体内に存在する腎結石等の結石の位置を検出する装置に関する。
【背景技術】
【０００２】
超音波は医療分野にて診断や治療に用いられている。例えば、超音波診断装置は、超音波ビームで生体内を走査し、生体内の音響インピーダンスの差に応じて生成される断層画像を診断の用に供する。また、超音波のエネルギーを患部等に集中させて破壊・壊死させる治療装置も存在する。例えば、近年では、高密度焦点式超音波（High intensity focused ultrasound：ＨＩＦＵ）法を用いた癌治療が行われており、この治療では、患部に集中させた超音波により患部だけを高温度に加熱することができ、癌組織を熱凝固、壊死させる。
【０００３】
また、ＨＩＦＵ法を用いた腎臓結石の非侵襲な破砕の研究も行われている。ＨＩＦＵ法で破砕された結石は現存の体外衝撃波結石破砕術（ＥＳＷＬ）よりも細かくなり、破砕片を無痛で尿道から排出させることが可能となる。
【０００４】
ＥＳＷＬ等においては、衝撃波を結石に集束させるためなどのため、体内における結石の位置を検出することが必要である。この結石の位置検出はＸ線を用いて行われている。また、Ｘ線とは異なり放射線被曝がないという利点から、超音波断層画像を用いた結石位置の検出も提案されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００５−０００６３３号公報
【特許文献２】特開平０６−１８９９７３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
ＥＳＷＬやＨＩＦＵによって体外から結石を破砕する治療においてはその安全性及び精度を確保するために結石の位置をリアルタイムで精度良く自動検出できることが重要な課題となる。特に、ＨＩＦＵは超音波の集束度合いが高いので、治療を効果的に行う観点から、高い位置検出精度が要求される。また、ＨＩＦＵは焦点領域に高いエネルギーを集中させるため、他の組織への誤照射を避ける必要があり、この点でも高い位置検出精度が求められる。
【０００７】
超音波断層像を用いて、結石の位置を検出する際には、結石以外の領域に生じるスペックルパターンなどのノイズの中から結石に起因する領域を抽出する必要があり、ノイズの影響の除去が精度・信頼性の向上のための課題となる。
【０００８】
特にＨＩＦＵによる治療中は集束超音波がその焦点に位置する結石の周りにマイクロ・バブルを発生させ、これは超音波画像上にて高輝度のノイズ領域を形成する。この高輝度のノイズ領域は、結石に起因する高輝度領域と一体的に超音波断層像から抽出され得る結果、超音波断層像上での結石の追跡を困難にする要因になる。
【０００９】
また、結石が存在する腹部は呼吸等によって運動しており、それに伴い結石も運動し得る。よって、この運動に対する補償も検出位置の精度向上に関する課題となる。
【００１０】
本発明は上記課題を解決するためになされたものであり、超音波断層像を用いて、生体内における結石を好適なリアルタイム性及び位置精度で自動検出することができる生体内結石検出装置を提供することを目的とする。
【課題を解決するための手段】
【００１１】
本発明に係る生体内結石検出装置は、送出した超音波に対する反射波の強度に応じた画素データを有する生体内の超音波断層像に基づいて、前記生体内の結石を検出するものであって、前記画素データの原データ分布からレイリー分布に従う成分を減じて、前記原データ分布のピーク位置での強度を抑圧し、これにより残る残差データ分布を求める残差分布抽出手段と、前記ピーク位置よりも前記反射波の強度が大きなデータ範囲であって、前記原データ分布の強度に対する前記残差データ分布の強度が所定の割合以上となる抽出対象データ範囲を求め、当該抽出対象データ範囲に属する画素からなる高反射領域を前記結石の候補領域として抽出する候補領域抽出手段と、を有する。この生体内結石検出装置における前記候補領域抽出手段は、前記残差データ分布と前記原データ分布から減じた前記レイリー分布とが等しい強度となる前記画素データを閾値とし、当該閾値以上を前記抽出対象データ範囲とする構成とすることができる。また、上記生体内結石検出装置における前記残差分布抽出手段は、前記原データ分布と前記レイリー分布とがそれぞれのピーク点と前記ピーク位置を超えるデータ範囲に存在する変曲点とで互いに一致するように前記レイリー分布を定める構成とすることができる。
【００１２】
別の本発明に係る生体内結石検出装置は、さらに、前記超音波の送受波位置から見た前記高反射領域の背後領域からの前記反射波の強度を、その隣接領域からの前記反射波の強度と比較する背後領域反射強度判定手段を有し、前記候補領域抽出手段が、前記隣接領域よりも前記反射波の強度が小さい前記背後領域を伴う前記高反射領域を前記結石の候補領域として抽出するものである。この生体内結石検出装置において、前記背後領域反射強度判定手段は、前記高反射領域を通過する超音波ビームの走査線における前記背後領域に属する注目区間での前記画素データの平均値である背後データ強度と、前記高反射領域の両脇を通過する２本の走査線についての前記注目区間に対応する区間での前記画素データの平均値である隣接データ強度との比を求め、前記候補領域抽出手段は、前記隣接データ強度に対する前記背後データ強度の比が所定の閾値未満となる前記高反射領域を前記候補領域として抽出する構成とすることができる。
【００１３】
さらに別の本発明に係る生体内結石検出装置は、検出対象とする前記結石であって予め輪郭が特定されている特定結石について当該輪郭を表すテンプレートを記憶するテンプレート記憶手段を有し、前記候補領域抽出手段は、前記テンプレートと前記高反射領域の輪郭とが適合する前記テンプレートの配置を探索し、当該配置の前記テンプレートに対応する前記高反射領域を前記特定結石の候補領域として抽出するテンプレート照合手段を有し、前記テンプレート記憶手段は、前記テンプレートを表す情報として、前記特定結石の輪郭に沿って所定の順序で配置される複数の代表点それぞれについて、当該代表点での前記結石の輪郭の形状に関する形状特徴情報と、前記結石の輪郭に対して所定の位置に設定される基準点と当該代表点との相対的な位置関係を表す特徴ベクトル情報とを記憶し、前記テンプレート照合手段は、前記各代表点毎に、前記高反射領域の輪郭上にて前記形状特徴情報に適合する適合点を探索し、当該適合点が当該代表点であるとした場合に前記基準点に対応する点となる仮想基準点を前記特徴ベクトル情報に基づいて求め、前記複数の代表点に対応して複数得られる前記仮想基準点の分布から前記超音波断層像における前記基準点の位置を推定し、当該推定された基準点に基づいて前記超音波断層像での前記テンプレートの配置を決定するものである。
【００１４】
本発明の好適な態様は、前記複数の代表点が、前記結石の輪郭形状に応じた凸多角形の頂点である凸頂点を含み、前記形状特徴情報が、前記各凸頂点に接続される前記凸多角形の辺の向きを含む生体内結石検出装置である。この態様において、前記凸頂点は、前記結石の輪郭上に設定される探査の始点、及び当該始点から前記結石の輪郭を所定の凸頂点探査手順で所定の回転方向に一周順次探査して得られる複数の点であり、前記凸頂点探査手順は、前記回転方向を示すベクトルをｈＲとし、先行して得られた前記凸頂点を基準凸頂点とし、当該基準凸頂点から次の前記凸頂点の候補点へ向かうベクトルｈＣ１と前記基準凸頂点から次の探査点へ向かうベクトルｈＳ１との外積（ｈＣ１×ｈＳ１）を前記ベクトルｈＲのｋ１倍（ここでｋ１は実数である）と置いたときに、前記ｋ１がε１＞０なる所定の閾値ε１以下である場合には、前記候補点を前記探査点で置き換えて探査を継続し、前記ｋ１が閾値ε１より大きい場合には、現在の前記候補点を前記基準凸頂点の次の前記凸頂点と定めるものとすることができる。
【００１５】
本発明の他の好適な態様は、前記複数の代表点がさらに、前記凸頂点で前記結石の輪郭を区切って得られる各輪郭区間における凹頂点を含み、前記形状特徴情報が、前記各輪郭区間における前記凸頂点及び前記凹頂点を順次接続する辺の向きを含む生体内結石検出装置である。この態様において、前記凹頂点は、前記結石の輪郭を周回する所定の回転方向を探査順序の方向として、前記輪郭区間をその一方端から所定の凹頂点探査手順で順次探査して得られる点であり、前記凹頂点探査手順は、当該探査の始点となる前記凸頂点又は先行して得られた前記凹頂点を基準頂点とし、当該基準頂点から次の前記凹頂点の候補点へ向かうベクトルｈＣ２と前記基準頂点から次の探査点へ向かうベクトルｈＳ２との外積（ｈＣ２×ｈＳ２）を前記ベクトルｈＲのｋ２倍（ここでｋ２は実数である）と置いたときに、前記ｋ２がε２＜０なる所定の閾値ε２以上である場合には、前記候補点を前記探査点で置き換えて探査を継続し、前記ｋ２が閾値ε２未満である場合には、現在の前記候補点を前記凹頂点と定めるものとすることができる。
【００１６】
他の本発明に係る生体内結石検出装置は、さらに、前記テンプレートの配置に基づいて把握される、過去の時系列での前記結石の位置からなる結石軌跡情報を記憶する結石軌跡記憶手段と、前記結石の位置に関する前記生体の呼吸に応じた周期的な時間変化を表す運動モデルを記憶する運動モデル記憶手段と、現時刻における前記結石軌跡情報を前記運動モデルと照合し、前記運動モデルの周期における現時刻の位相を検出し、所定時間後の前記結石の予測位置を求める結石位置予測手段と、を有するものである。
【発明の効果】
【００１７】
本発明によれば、原データ分布に占める残差データ分布の割合が所定値以上となるデータ範囲の画素を抽出することで、生体内の軟組織が生じやすいスペックルノイズに対応する画素を好適に除去し、結石が存在し得る高反射領域を好適に抽出することができる。
【００１８】
また、結石は超音波を透過しにくく、その背後領域にいわば影となる反射波強度が低い領域を生じ得る。本発明によれば、当該影の有無の判断に基づいて、高反射領域のうち結石に対応する部分をより好適に抽出することができる。
【００１９】
本発明によれば、輪郭が特定されている特定結石について、高反射領域のうち当該特定結石の輪郭に適合する部分を抽出することによって、例えば、結石の周囲の一部にマイクロ・バブルが生じている場合などのように、結石による高反射領域の近傍に高反射のノイズ領域が付随している場合でも、結石の位置を精度良く検出することができる。
【００２０】
本発明によれば、呼吸に応じて運動する結石の位置を好適に検出することができる。
【図面の簡単な説明】
【００２１】
【図１】本発明の実施形態である超音波診断・治療システムの概略のブロック図である。
【図２】原データ分布、及びこれに対する一般化レイリー曲線を示す模式図である。
【図３】原データ分布、一般化レイリー曲線及び残差データ分布を示す模式図である。
【図４】高反射領域を近似する楕円形のモデルを示す模式図である。
【図５】楕円モデルのパラメータを説明する模式図である。
【図６】ターゲットの輪郭に配置された代表点の一つ、基準点、及び特徴ベクトルを示す模式図である。
【図７】輪郭、及び近似的な凸包モデルを示す模式図である。
【図８】凸頂点探査手順を説明するための模式図である。
【図９】図７に示した輪郭及び凸包モデルに対して得られた簡易輪郭モデルの例を示す模式図である。
【図１０】凹頂点を有する輪郭区間ｐｋｐｋ＋１の模式図である。
【発明を実施するための形態】
【００２２】
以下、本発明の実施の形態（以下実施形態という）について、図面に基づいて説明する。
【００２３】
図１は実施形態である超音波診断・治療システム２の概略のブロック図である。本システムは、超音波断層像に基づいて患者の腎臓結石を検出し、検出した結石をＨＩＦＵ法により破砕する。本システム２は、エンドエフェクタ４、移動機構６、本体８を有する。エンドエフェクタ４はプローブ１０と圧電トランスデューサ１２とを備え、本システム２による診断・治療時に患者に当接される。移動機構６は、本体８からの制御に応じて、エンドエフェクタ４の位置・姿勢を変える。
【００２４】
プローブ１０は、超音波パルスの送出及びエコーの受波を行う超音波探触子である。このプローブ１０は振動子アレイを有している。振動子アレイは患者の体内へ超音波ビームを送出する。超音波ビームは電子的にアレイ方向に走査される。超音波ビームの送出はパルス状に行われ、振動子アレイは超音波ビームの送出後、体内からのエコーを受波する。体内での結石の三次元的な移動に対応できるように、２つの振動子アレイが互いの走査面を基本的に直交させてエンドエフェクタ４に配置される。２つの振動子アレイは、交互に送受波を行う。
【００２５】
圧電トランスデューサ１２は、ＨＩＦＵ治療のための集束超音波を発生する。プローブ１０及び圧電トランスデューサ１２は、プローブ１０の２つの振動子アレイの走査面の交線上に集束超音波の焦点位置が存在するように、互いの位置関係が設定されている。
【００２６】
プローブ１０及び圧電トランスデューサ１２は、本体８からの駆動信号に応じて超音波を送出する。またプローブ１０の受信信号は本体８へ送られる。
【００２７】
本体８は、送受波回路２０、ＤＳＣ（Digital Scan Converter）２２、パルス発生回路２４、演算処理装置２６、記憶装置２８を含んで構成されている。演算処理装置２６は記憶装置２８に格納されたプログラム等に基づいて動作し、例えば、結石検出部３０や制御部３２として機能する。
【００２８】
送受波回路２０は、制御部３２による制御に従って、プローブ１０から体内への超音波の送信、及びプローブ１０が受信したエコー信号の処理を行う。送信時には、送受波回路２０は、振動子アレイの各チャネルごとに、対応するチャネルの振動子を励振駆動するための送信パルスを生成しプローブ１０へ出力する。送受波回路２０は、チャネル毎に送信パルスに与える遅延量を調整して、プローブ１０から送信される超音波が所望の方向に送信ビームを形成するように各振動子の励振タイミングを制御する。
【００２９】
一方、受信時には、送受波回路２０は、プローブ１０から振動子アレイの各チャネルごとの受信信号を入力される。送受波回路２０は、チャネルごとの受信信号を増幅した後、チャネル間の受信信号の位相差を調整して互いに加算する整相加算処理を行って受信ビームを形成する。また、送受波回路２０は、受信信号をアナログ信号からデジタル信号に変換し、受信信号を超音波ビームの方向に沿ったエコーデータ列として出力する。
【００３０】
振動子アレイの電子走査によって超音波ビームが電子的に走査されることによって一つの走査面が形成され、この走査面から１フレームのエコーデータが取得される。ＤＳＣ２２は、送受波回路２０から出力される超音波ビームに沿ったエコーデータ列を、内蔵のフレームメモリに一旦格納する。ＤＳＣ２２は、エコーデータが定義される空間座標系を、超音波ビームの走査方向（ｉ座標とする。）とビーム方向の深さ（ｊ座標とする。）とで指定される送受波座標系（円座標系）から、断層画像表示を行う表示装置（図示せず。）の走査方式に適合した二次元直交座標系（水平方向をｕ座標、垂直方向をｖ座標とする。）に変換する。
【００３１】
パルス発生回路２４は、制御部３２の制御に従い、圧電トランスデューサ１２への駆動信号を生成し出力する。
【００３２】
結石検出部３０は、本発明に係る生体内結石検出装置に相当する。結石検出部３０は、ＤＳＣ２２からエコー強度に応じた画素データを有する生体内の超音波断層像を読み込み、当該超音波断層像に基づいて演算処理を行い生体内の結石を検出する。なお、結石検出部３０での演算に用いる超音波断層像のデータを、座標変換前のものと変換後のものとのいずれにするかは任意に設計することができる。結石検出部３０の演算処理については後述する。
【００３３】
制御部３２は、本システム２の各部の動作を制御する。例えば、制御部３２は、送受波回路２０やパルス発生回路２４を制御する。また制御部３２は、結石検出部３０から得られる結石の位置に基づいて移動機構６を制御することができ、例えば、治療時には、結石の運動に応じてエンドエフェクタ４を移動させ、圧電トランスデューサ１２から発せられる集束超音波の焦点を結石に追従させる。
【００３４】
次に結石検出部３０の機能・処理等について説明する。結石検出部３０は、（１）超音波断層像の画素データの分布に基づいて所定の高反射の領域を結石の候補領域として抽出する処理（以下、高反射領域抽出処理）、（２）結石の背後ではエコーが弱くなるという特徴に基づいて結石の候補領域を特定する処理（以下、影判定処理）、（３）予め得られた結石の輪郭に適合する候補領域を探索する処理（以下、テンプレート照合処理）、及び、（４）呼吸に応じて動く結石の位置を予測する処理（以下、位置予測処理）を備えている。以下、各処理について順番に説明する。
【００３５】
（１）高反射領域抽出処理
腎臓結石は、その周囲の生体の軟組職に比べて音響インピーダンスが高く、超音波に対して高い反射を示す。なお、一般に、超音波断層像の画素データはエコー強度が強いほど大きな値で定義され、表示装置上にて高い輝度で表現される。ちなみに、結石の周囲には超音波断層像にて背景となるものとして腎臓、筋肉、血液、水分が存在し、これらと比較した結石の反射率は、7.09〜12.75倍であるとの研究結果が報告されている。
【００３６】
結石検出部３０は、超音波断層像中の結石を背景から分離するために、画素データのヒストグラム上で当該背景の部分に起因する分布をモデリングして、結石に起因し得る高反射の領域を抽出するための閾値を決定する。
【００３７】
以下、この処理を説明する。結石検出部３０は、関心領域（Region of Interest：ＲＯＩ）の画素データを抽出する。なお、ＲＯＩの設定は、例えば、表示装置に表示される超音波断層画像に基づいて、医師等の操作者が本システム２のキーボード、マウス、タッチスクリーン等の操作手段（図示せず。）を用いて行うことができる。
【００３８】
結石検出部３０は、抽出したＲＯＩの画素データのヒストグラムを生成し、さらに、当該画素データの分布の概形を抽出するために、ヒストグラムを平滑化する。例えば、重み付き移動平均等で平滑化する。
【００３９】
超音波断層画像においてスペックルパターンを生じる軟組織の画素データは、レイリー分布（Rayleigh probability density function）に好適に従う。そこで、結石検出部３０は、平滑化したヒストグラムを画素データの原データ分布として、当該原データ分布からレイリー分布に従う成分を減じて、背景以外の画素データに対応する残差データ分布を求める。ここで、一般に、結石以外の背景がＲＯＩ中にて主要な比率を占めるので、原データ分布のピークは基本的に背景に起因する。よって、結石検出部３０は、原データ分布のピーク位置での強度を抑圧するようにレイリー分布と相似の関数を定める。具体的には、レイリー分布と相似の関数として下記（１）式で表現される一般化レイリー曲線の関数を仮定し、そのパラメータσ，Ｓ，ｘ０を定める。なお、変数ｘが画素データを表す。
【００４０】
【数１】
2010207427.tif 000003
【００４１】
パラメータは、原データ分布のピーク点ＰＨＰ（≡（Ｇｐｋ，Ｎｐｋ））及び変曲点ＰＨＩ（≡（Ｇｉｎｆ，Ｎｉｎｆ））と一般化レイリー曲線のピーク点ＰＲＰ及び変曲点ＰＲＩとが互いに一致するように定められる。この条件は下記（２）式及び（３）式で表される。
【００４２】
【数２】
2010207427.tif 000004
【００４３】
【数３】
2010207427.tif 000005
【００４４】
（２）式，（３）式それぞれの左辺は一般化レイリー曲線のピーク点ＰＲＰ及び変曲点ＰＲＩの座標である。ピーク点ＰＲＰは（１）式の１階微分が０となることに基づいて、また変曲点ＰＲＩは２階微分が０となることに基づいてそれぞれ定められる。
【００４５】
原データ分布のピーク点ＰＨＰは原データ分布が最大となる点に定めることができる。原データ分布を表す曲線は必ずしも解析的な関数ではないので、その変曲点ＰＨＩを一般化レイリー曲線におけると同様に２階微分が０となるという条件から決めることは困難である。そこで本実施形態では、一般化レイリー曲線におけるピーク点ＰＲＰでの度数（画素数）に対する変曲点ＰＲＩでの度数の比が、（２）式，（３）式より、一定値ｋ（ｋ≒0.637）となるという特徴を利用して変曲点ＰＨＩを求める。すなわち、原データ分布上にてピーク点ＰＨＰでの度数Ｎｐｋに対する変曲点ＰＨＩでの度数Ｎｉｎｆの比（Ｎｉｎｆ/Ｎｐｋ）がｋに最も近い値となる点を変曲点ＰＨＩと定義する。なお、一般化レイリー曲線の変曲点ＰＲＩはピーク点ＰＲＰより大きな画素データ範囲に存在することに対応して、原データ分布の変曲点ＰＨＩもＧｐｋを超える画素データ範囲にて探索される。
【００４６】
（２）式，（３）式の左辺に現れるＰＲＰ，ＰＲＩの座標を用いて一般化レイリー曲線のピーク点ＰＲＰと変曲点ＰＲＩとを結ぶ直線の傾きａＲを計算すると、パラメータＳ，ｘ０を消去できａＲはσだけに依存する。よって、例えば、原データ分布のピーク点（Ｇｐｋ，Ｎｐｋ）と変曲点（Ｇｉｎｆ，Ｎｉｎｆ）とを結ぶ直線の傾きａＨがａＲに一致するという条件からまずσを決定し、そのσを（２）式，（３）式に代入して、Ｓ及びｘ０を決定することができる。
【００４７】
図２は原データ分布４０と一般化レイリー曲線４２とを示す模式図であり、横軸が画素データの値、縦軸が画素数である。軟組織からの弱いエコーは一般化レイリー曲線４２に好適に従うので、結石が存在しないＲＯＩから得られた原データ分布４０は一般化レイリー曲線４２によく一致し得る。これに対し、ＲＯＩに結石が存在する場合には、画素データが小さい領域（例えば、変曲点ＰＨＩより左側の領域）では原データ分布４０と一般化レイリー曲線４２とは結石が存在しない場合と同様によく一致し得るが、画素データが大きい領域では、結石に起因する大きなデータ値の画素の存在により、原データ分布４０は一般化レイリー曲線４２から乖離し得る。
【００４８】
結石検出部３０は上述のように、原データ分布に適合する一般化レイリー曲線を決定すると、原データ分布から一般化レイリー曲線が表す強度（度数）を減じて、残差データ分布を求める。一般化レイリー曲線はスペックルパターンなどに対応し、その成分を減算することで、結石検出におけるノイズ成分が除去される。すなわち、残差データ分布を形成する画素は、結石に対応するものである可能性が高い。
【００４９】
図３は原データ分布４０、一般化レイリー曲線４２、及び残差データ分布４４を示す模式図であり、図２と同様、横軸が画素データの値、縦軸が画素数である。結石検出部３０は、原データ分布のピーク位置Ｇｐｋよりも値が大きな範囲にて、原データ分布の強度に対する残差データ分布の強度が所定の割合以上となる画素データを閾値として定め、当該閾値以上を抽出対象データ範囲とし、当該抽出対象データ範囲に属する画素からなる高反射領域を結石の候補領域として抽出する。
【００５０】
本実施形態では、結石検出部３０は、残差データ分布と一般化レイリー曲線とが等しい強度となる画素データを閾値Ｇｔｈ（図３参照）とする。この閾値Ｇｔｈは一般化レイリー曲線４２と残差データ分布４４とが交差するデータ値であり、閾値Ｇｔｈを超えると、残差データ分布４４が一般化レイリー曲線４２を上回る。このことから、閾値Ｇｔｈ以上の高反射領域にはノイズに起因する部分は少なく結石の像が占める割合が高いことが期待できる。そこで、結石検出部３０は、閾値Ｇｔｈ以上の画素からなる高反射領域を、結石の候補領域として抽出する。
【００５１】
（２）影判定処理
一方、マイクロ・バブルのように、結石以外にも高反射を生じる場合がある。そこで、結石検出部３０は、上述の高反射領域抽出処理で抽出される高反射領域に基づく候補領域から結石の像以外の部分を除外し、候補領域を絞り込むフィルタリング処理の一つとして、影判定処理を行う。
【００５２】
結石は、その周囲の組織等に比べて高い音響インピーダンスを有するので、結石に向けて送出された超音波のほとんどはその表面で反射され、透過する成分はわずかとなる。そのため、結石の背後にはエコー強度の低い低反射領域が生じる。この現象をここでは「音響的影効果」と称し、低反射強度である背後領域を「影」領域と称する。
【００５３】
本発明における影判定処理の原理は、プローブ１０から見た高反射領域の背後領域からのエコー強度を、その隣接領域からのエコー強度と比較するというものである。結石検出部３０は背後領域反射強度判定手段としてこの影判定処理を行い、候補領域抽出手段として、隣接領域よりも低反射強度である背後領域を伴う高反射領域を結石の候補領域として抽出する処理を行う。
【００５４】
結石検出部３０は、背後領域反射強度判定手段として、例えば、（４）式に示すように、高反射領域を通過する超音波ビームの走査線（ｉ＝ｉＣ）における背後領域に属する注目区間（ｊｓｔ≦ｊ≦ｊｅｄ）での画素データＩ(ｉＣ,ｊ)の平均値である背後データ強度と、高反射領域の両脇を通過する２本の走査線（ｉ＝ｉＲ及びｉＬ）についての注目区間に対応する区間での画素データＩ(ｉＲ,ｊ)及びＩ(ｉＬ,ｊ)の平均値である隣接データ強度との比で定義する影係数ｆｓｈａｄｏｗを求める。そして、結石検出部３０は、影係数ｆｓｈａｄｏｗが所定の閾値未満となる高反射領域を結石の候補領域として抽出する。
【００５５】
【数４】
2010207427.tif 000006
【００５６】
結石検出部３０は、例えば、図４に示すように高反射領域５０を近似する楕円形のモデル５２を求め、その楕円の中心（ｕＣ，ｖＣ）を通過するビーム５４の座標をｉ＝ｉＣとし、当該楕円に接する２本のビームの座標５６，５８をｉ＝ｉＲ，ＩＬとする。また、楕円中心（ｕＣ，ｖＣ）から長軸半径(ｌ/２)深い深さをｊｓｔと設定することができる。なお、ｊｅｄは（ｊｅｄ−ｊｓｔ）が所定の値となることに基づいて定めることができ、例えば、（ｊｅｄ−ｊｓｔ）は楕円の長軸の長さや、ｉＲとｉＬとの差など、高反射領域の大きさに応じた値に設定することができる。
【００５７】
図５は、楕円モデルのパラメータを説明する模式図である。楕円モデルは、長軸ｌ、短軸ｗ、中心座標ｕＣ，ｖＣ、及び傾斜角θからなる５つのパラメータで表すことができる。これらのパラメータは、近似対象とする高反射領域を構成する各画素（ｕｉ，ｖｉ）について（５）式で定義されるモーメントｍｐｑに基づき、（６）〜（８）式により求められる。
【００５８】
【数５】
2010207427.tif 000007
【００５９】
【数６】
2010207427.tif 000008
【００６０】
【数７】
2010207427.tif 000009
【００６１】
【数８】
2010207427.tif 000010
【００６２】
なお、結石検出部３０は、上述の楕円モデルそのものに基づいて別途のフィルタリングを行って、これを上述の音響的影効果に基づくフィルタリングと併用して結石の候補領域を絞り込むことができる。具体的には、０次モーメントｍ００で与えられる楕円の面積Ｓや、（１−ｗ/ｌ）で定義される扁平率Ｃを閾値判定して、候補領域をフィルタリングすることができる。
【００６３】
例えば、ＨＩＦＵ法などによる破砕が行われるのは４ｍｍ程度より大きな腎臓結石であることに対応して、面積Ｓの閾値を定めることができる。例えば、結石検出部３０は、Ｓが6.45ｍｍ２を超える領域を抽出する。また、結石の像は一般に円形に近く、扁平率はそれほど大きくならないことから、結石検出部３０は、Ｃが0.6未満の候補領域を抽出する。また、結石検出部３０は上述の影係数ｆｓｈａｄｏｗが0.7未満の影領域を伴う候補領域を抽出するように設定することができる。
【００６４】
（３）テンプレート照合処理
既に述べたように、超音波治療の際に集束超音波によって発生するマイクロ・バブルは超音波画像上に高輝度領域を形成する。集束超音波の焦点は破砕対象の結石、すなわち、ターゲットの位置に設定しているため、焦点の周りに発生した高輝度のノイズは実際の結石の像と一体の高反射領域として捉えられる。その結果、例えば、上述のように楕円モデルなどを用いて高反射領域の重心点を結石の位置として特定する際に当該位置の精度が低下し、ターゲットの追跡が困難になり得る。
【００６５】
テンプレート照合処理は、この問題への解決策を与えるものである。集束超音波の焦点領域は１ｍｍ程度と小さく、ノイズもその近傍で生成されるため、ターゲット物体の像は基本的にノイズによって部分的にしか影響を受けない。そこで、ターゲット物体においてノイズによって影響されていない部分の形状を認識することにより、実際のターゲット像の重心点を検出する。
【００６６】
ＨＩＦＵによる治療対象とする結石は、集束超音波をしない状態での超音波断層像にて特定できる。結石検出部３０は、当該超音波断層像を二値化して、特定結石（ターゲット）の高反射領域の輪郭を、集束超音波によるノイズの影響を受けていない状態で抽出し、得られたターゲット像の形状を、後述する照合処理に用いるテンプレートとする。得られたテンプレートを表す情報（ターゲットモデルと称する）は記憶装置２８に記録される。
【００６７】
集束超音波の照射による治療が開始されると、現在の超音波断層像の二値化画像においてターゲットに対応する高反射領域の位置が、ターゲットの過去の運動軌跡から推定される。結石検出部３０は、当該位置の高反射領域がターゲットの像を包含するものであるとして、その輪郭に好適に適合するテンプレートの配置を探索する。そして、高反射領域のうち、探索により決定された配置のテンプレートに対応する部分をターゲットの候補領域として抽出する。
【００６８】
記憶装置２８に格納されるターゲットモデルは、ターゲットの輪郭に沿って所定の順序で配置される複数の代表点それぞれについて、当該代表点でのターゲットの輪郭の形状に関する形状特徴情報と、ターゲットの輪郭に対して所定の位置に設定される基準点と当該代表点との相対的な位置関係を表す特徴ベクトル情報とを含んでいる。
【００６９】
このターゲットモデルに基づいて、結石検出部３０は、テンプレート照合処理として、各代表点毎に、高反射領域の輪郭上にて形状特徴情報に適合する適合点を探索し、当該適合点が代表点であるとした場合に基準点に対応する点となる仮想基準点を特徴ベクトル情報に基づいて求める。そして、複数の代表点に対応して複数得られる仮想基準点の分布から現在の超音波断層像における基準点の位置を推定し、当該推定された基準点に基づいて超音波断層像でのテンプレートの配置を決定する。
【００７０】
図６は、ターゲットの輪郭６０に配置された代表点６２の一つ、基準点６４、及び特徴ベクトル６６を示す模式図である。
【００７１】
輪郭６０は、二値化された超音波断層像における高反射領域の境界であり、その境界上に位置する画素の連鎖からなる。代表点６２は、当該輪郭６０上の各画素とすることもできるが、一般には当該画素の数は大きく、処理負荷が高くなる。そこで、輪郭６０の形状を凸包モデル（convex hull model）で近似し、その凸多角形の頂点である複数の凸頂点を代表点６２とすることで、代表点６２の数を減らすことができる。また、この場合の各代表点での輪郭６０の形状の特徴は、各凸頂点に接続される凸多角形の辺の向きを含む形状特徴情報で表すことができる。
【００７２】
ここで、凸包モデルを厳密に適用し、内側に輪郭６０が完全に包含されるように凸多角形を定めると、凸頂点の数は一般に多くなる。そこで、本実施形態では、近似的に凸包モデルを適用することで凸頂点の数を減らして、演算処理装置２６の処理負荷を軽減し、治療中のリアルタイムでのターゲットの追跡を可能とする。代表点の数は、処理負荷の軽減の観点からは少なくされるが、テンプレート照合で決定される位置精度の向上の観点からは或る程度の数の確保が必要であり、これらのトレードオフに基づいて適切な数の凸頂点の数が設定される。
【００７３】
結石検出部３０は、凸頂点として、ターゲットの輪郭６０上に設定される探査の始点、及び当該始点から輪郭６０を所定の凸頂点探査手順で所定の回転方向に一周順次探査して得られる複数の点を抽出する。
【００７４】
図７は輪郭６０と、その近似的な凸包モデル６８とを示す模式図である。図７に示す凸包モデル６８は、始点ｐ１から反時計回りに凸頂点探査を行って１０個の凸頂点ｐｉ（ｉ＝１〜１０）が得られた例を示している。図８は、凸頂点探査手順を説明するための模式図である。超音波断層像上での水平方向の座標軸ｕは右向きが正の方向、垂直方向の座標軸ｖは下向きが正の方向に設定される。ここで、超音波断層像に垂直な軸ξを導入して三次元直交座標（ｕ，ｖ，ξ）を定義する。ξ軸は超音波断層像を臨む視点から遠ざかる方向を正に設定される。輪郭６０上を反時計回りに周回して探査を行う場合、その周回方向（回転方向）を示すベクトルｈＲはξ軸の負の方向に向かう。結石検出部３０は、先行して得られた凸頂点ｐｋ（ｋ＝１〜９）を基準凸頂点とし、当該基準凸頂点から次の凸頂点ｐｋ＋１の候補点ｐＣへ向かうベクトルｈＣ１と基準凸頂点から次の探査点ｐＳへ向かうベクトルｈＳ１との外積（ｈＣ１×ｈＳ１）をベクトルｈＲのｋ１倍（ここでｋ１は実数である）と置いたときに、ｋ１がε１＞０なる所定の閾値ε１以下である場合には、候補点ｐＣを現在の探査点ｐＳで置き換えて探査を継続し、ｋ１が閾値ε１より大きい場合には、現在の候補点ｐＣを基準凸頂点の次の凸頂点ｐｋ＋１と定める。
【００７５】
具体的には、ベクトルｈＣ１，ｈＳ１のξ成分は０であることから、外積（ｈＣ１×ｈＳ１）の成分のうち非零となり得るのはξ成分のみであり、これをＣＰ１と表す。ｐｋ，ｐＣ，ｐＳのｕ−ｖ平面上の座標をそれぞれ（ｘ０，ｙ０），（ｘ１，ｙ１），（ｘ２，ｙ２）と表すと、ＣＰ１は次式で与えられる。
ＣＰ１＝(ｘ１−ｘ０）(ｙ２−ｙ０）−(ｘ２−ｘ０）(ｙ１−ｙ０）      （９）
【００７６】
便宜上、ベクトルｖＲの大きさを１とすれば、ｖＲ＝（０，０，−１）であり、ｋ１＝−ＣＰ１となる。ｐＳがベクトルｈＣ１と同じ方向にあればＣＰ１＝０（又はｋ１＝０）であり、ｐＳがベクトルｈＣ１の向きにより右側に位置すればＣＰ１＞０（又はｋ１＜０）、左側に位置すればＣＰ１＜０（又はｋ１＞０）である。例えば、ε１が０であるとするならば、上述の探査手順は、輪郭６０を反時計回りに移動したときに左に曲がる点、すなわち凸頂点を探査することになり、この場合、凸包モデルを好適な精度で求めることができる。本実施形態では、ε１＞０なる閾値を設定することで、ε１の大きさに応じた或る程度の左曲がりは許容して凸頂点として検出しないことで凸頂点間の距離を広げ、凸頂点の数を減らした近似的な凸包モデル６８を得ることができる。なお、凸頂点の数は閾値ε１によって調節することができる。
【００７７】
結石検出部３０は、上述の凸頂点探査手順により、輪郭６０における凸頂点ｐｋを求める。また、結石検出部３０は各辺ｐｋｐｋ＋１の向きφを求め、各頂点ｐｋのインデックスｋにｐｋの座標及び、形状特徴情報として辺ｐｋｐｋ＋１の向きを対応付けて、記憶装置２８に格納する。なお、辺の向きφは、例えば、点ｐｋから点ｐｋ＋１に向かうベクトルがｕ軸の正の方向に対してなす角度で表すことができる。例えば、当該角度φ［rad］は−π＜φ≦πなる範囲で定義することができる。
【００７８】
また、結石検出部３０は、各凸頂点ｐｋについて、特徴ベクトル６６を表す情報として、ｐｋから基準点６４に向かうベクトルの長さｒ及び向きαを求め、記憶装置２８に格納する。例えば、基準点６４を楕円モデル５２の中心（ｕＣ，ｖＣ）に設定して特徴ベクトルを定義することができる。向きαはφと同様にｕ軸の正の方向に対してなす角度で表すことができる。
【００７９】
以上、凸包モデルに基づくターゲットモデルを説明した。上述の凸包モデルはテンプレートの一例であり、結石検出部３０はより好適に輪郭６０の形状特徴を表すテンプレートを採用することができる。例えば、上述の凸頂点に加えて、凸頂点で輪郭６０を区切って得られる各輪郭区間における凹頂点を代表点６２に含めたモデル（簡易輪郭モデルと称する）を用いれば、より精度の高いテンプレート照合が可能となる。
【００８０】
以下、簡易輪郭モデルについて説明する。図９は図７に示した輪郭６０及び凸包モデル６８に対して得られた簡易輪郭モデル８０の例を示す模式図である。
【００８１】
結石検出部３０は、輪郭６０を所定の回転方向に辿り、各輪郭区間ｐｋｐｋ＋１にて凹頂点を探査する。図１０は凹頂点を有する輪郭区間ｐｋｐｋ＋１の模式図である。ここでは、輪郭６０を反時計回りに周回する方向に凹頂点探査を行う例を説明する。この場合、凸頂点ｐｋとｐｋ＋１とに挟まれた輪郭区間において、ｐｋからｐｋ＋１へ向けて探査が行われる。当該輪郭区間にて設定された凹頂点をｐｋ側からｐｋ＋１側へ向けて順に付すインデックスｊ（ｊ＝１〜ｊｍ）を用いてｑｊと表す。また、当該輪郭区間にて、ｐｋからｑ１へ向かうベクトルをｇ０，ｑｊからｑｊ＋１（１≦ｊ≦ｊｍ−１）に向かうベクトルをｇｊ，ｑｊｍからｐｋ＋１に向かうベクトルをｇｊｍと表す。
【００８２】
結石検出部３０は、探査の始点となる凸頂点ｐｋ又は先行して得られた凹頂点ｑｊを基準頂点とし、基準頂点から次の凹頂点の候補点ｑＣへ向かうベクトルｈＣ２と基準頂点から次の探査点ｑＳへ向かうベクトルｈＳ２との外積（ｈＣ２×ｈＳ２）を、探査を行う周回方向を示すベクトルｈＲのｋ２倍（ここでｋ２は実数である）と置いたときに、ｋ２がε２＜０なる所定の閾値ε２以上である場合には、候補点ｑＣを現在の探査点ｑＳで置き換えて探査を継続し、ｋ２が閾値ε２未満である場合には、現在の候補点ｑＣを凹頂点と定める。
【００８３】
この凹頂点探査手順は、上述の凸頂点探査手順と比較すれば、当該凸頂点探査手順と対称なものであることが理解される。すなわち、凸頂点探査手順が基本的に輪郭６０を反時計回りに移動したときに左に曲がる点を探査することにより凸頂点を見つけるものであったのに対し、上述の凹頂点探査手順は基本的に右に曲がる点を探査するものであり、これにより凹頂点を見つけることができる。
【００８４】
凸頂点と凹頂点とを合わせた代表点の全体について、輪郭６０を例えば反時計回りに周回する方向に順にインデックスｋを付し、各代表点について、そのインデックスｋと当該代表点の座標及び、形状特徴情報を対応付けて記憶装置２８に格納する。例えば、簡易輪郭モデルにおける形状特徴情報は、凸包モデルと同様、代表点からなる多角形の各辺の向きφを含むことができる。
【００８５】
また、当該φと併せて、各代表点が属する輪郭区間の両端を結ぶ直線の向き、すなわち凸包モデルにおける辺の向きを表すφｃｖｘを形状特徴情報とすることができる。この場合、例えば、頂点ｐｋとｐｋ＋１とに挟まれる隣接区間内の各代表点に対しては、当該代表点を始点とするベクトルｇｊ（０≦ｊ≦ｊｍ−１）の角度φと、ｐｋからｐｋ＋１へ向かうベクトルの角度φｃｖｘとが形状特徴情報として対応付けられる。
【００８６】
さらに、結石検出部３０は、各代表点について、特徴ベクトル６６を表す情報として、凸包モデルの場合と同様、各代表点６２から基準点６４に向かうベクトルの長さｒ及び向きαを求め、記憶装置２８に格納する。
【００８７】
結石検出部３０は、集束超音波を照射している状態において得られる現在の超音波断層像の高反射領域と、記憶装置２８に格納されたターゲットモデルとを照合し、当該高反射領域中におけるテンプレートの配置を決定する。
【００８８】
このテンプレート照合処理では、結石検出部３０は、集束超音波を照射している状態である現在の高反射領域の輪郭の代表点及び形状特徴情報を、ターゲットモデルを決定したのと同様の方法で求める。そして、現在の各代表点に対応付けられた形状特徴情報とターゲットモデルの形状特徴情報とを照合することにより、ターゲットモデルの代表点それぞれについて、現在の高反射領域の代表点のうちで対応するものを適合点として求める。なお、この照合に際しては、ターゲットの回転や移動を考慮して、ターゲットモデルと現在の高反射領域との角度調整や倍率の調整を行うことが好適である。
【００８９】
得られた各適合点には、ターゲットモデルの特徴ベクトル６６が対応付けられ、適合点を始点とする特徴ベクトルの終点は仮想的な基準点となる。現在の高反射領域の輪郭のうち、マイクロ・バブル等のノイズの影響を受けていない部分の適合点からの特徴ベクトルの終点は一点に集中することが期待され、一方、ノイズの影響を受けて変形した輪郭部分にて得られた適合点については、当該適合点からの特徴ベクトルの終点は基本的にはばらつく。そこで、結石検出部３０は、既に述べたように、仮想基準点の分布から現在の超音波断層像における基準点の位置を推定し、当該推定された基準点に基づいて超音波断層像でのテンプレートの配置を決定する。例えば、終点にピークを有し、終点からの距離が大きくなるにつれて値が小さくなる二次元の分布関数の値を各終点の近傍画素にて定義し、各終点に起因する分布関数値の合計が最も高くなる画素を基準点の位置とすることができる。
【００９０】
（４）位置予測処理
本システム２は、集束超音波の照射により結石の破砕を行う。照射する期間内にて結石は移動し得るので、制御部３２は、上述のように検知された結石の位置に基づき、移動機構６を駆動してエンドエフェクタ４を移動させ、集束超音波の焦点位置をターゲットの運動に対して補償する。
【００９１】
その際、超音波断層像の取得、結石位置の決定、移動機構６の制御という一連の過程により制御遅延が生じて追随精度が低下し得る。そこで、本システム２では、フィードバック制御に加えてフィードフォワード制御を利用し、現在までに得られた結石の位置運動の周期的な成分に対してはフィードフォワード制御を適用する。ちなみに、フィードバック制御は、主に、運動の非周期的な成分に対して適用される。
【００９２】
結石の運動に影響を与える大きな要因の一つは患者の呼吸である。呼吸によるターゲットの運動は周期的な運動成分と非周期的な運動成分とから構成される。制御部３２は、結石検出部３０により得られたテンプレートの配置に基づいて把握される結石の位置を記憶装置２８に蓄積する。制御部３２は、記憶装置２８に蓄積された結石位置から、呼吸に応じた周期的な時間変化を表す運動モデルを抽出し、演算処理装置２６の内部メモリや記憶装置２８に記憶する。また、結石検出部３０は当該運動モデルは、過去の時系列での結石の位置からなる結石軌跡情報を記憶装置２８に記録する。
【００９３】
制御部３２は、現時刻における結石軌跡情報を運動モデルと照合し、運動モデルの周期における現時刻の位相を検出し、移動機構６の移動完了までの遅延時間に応じた所定時間後の結石の予測位置を求める。
【００９４】
具体的に、まず、制御部３２は、記憶装置２８に記録したターゲットの運動履歴をフーリエ解析し、例えば、最も大きな振幅を有する周波数を求め、その運動周期と振幅とで表される正弦波振動を運動モデルとする。次に、現在の患部の位置を周期的な運動モデルと同期させるため、結石軌跡情報として記憶装置２８に記録されている直近の運動履歴と、運動モデルとの相関値を計算し、当該相関値が最大となる位相を現在の運動の位相とする。そして、超音波断層像の取得から移動機構６の移動完了までの想定される遅延時間に対応したフィードフォワード項を運動モデルから生成し、フィードバック項と合成して移動機構６に対する操作量を決定する。
【符号の説明】
【００９５】
２  超音波診断・治療システム、４  エンドエフェクタ、６  移動機構、８  本体、１０  プローブ、１２  圧電トランスデューサ、２０  送受波回路、２２  ＤＳＣ、２４  パルス発生回路、２６  演算処理装置、２８  記憶装置、３０  結石検出部、３２  制御部、４０  原データ分布、４２  一般化レイリー曲線、４４  残差データ分布、５０  高反射領域、５２  楕円モデル、６２  代表点、６４  基準点、６６  特徴ベクトル、６８  凸包モデル、８０  簡易輪郭モデル。
【図１】
2010207427.tif 000011
【図２】
2010207427.tif 000012
【図３】
2010207427.tif 000013
【図４】
2010207427.tif 000014
【図５】
2010207427.tif 000015
【図６】
2010207427.tif 000016
【図７】
2010207427.tif 000017
【図８】
2010207427.tif 000018
【図９】
2010207427.tif 000019
【図１０】
2010207427.tif 000020