WO2006068271A1

WO2006068271A1 - 心臓弁データ計測方法および装置

Info

Publication number: WO2006068271A1
Application number: PCT/JP2005/023797
Authority: WO
Inventors: Nozomi Watanabe; Yasuo Ogasawara; Masashi Sakurai
Original assignee: Seikotec Co Ltd; YD Ltd
Current assignee: Seikotec Co Ltd; YD Ltd
Priority date: 2004-12-24
Filing date: 2005-12-26
Publication date: 2006-06-29
Anticipated expiration: 2007-06-24
Also published as: US20080085043A1; KR20070110837A; EP1832233A1; CN101111194A; JPWO2006068271A1; JP4452281B2; EP1832233A4

Abstract

　本発明が解決しようとする課題は、心臓における僧帽弁のtenting volume、tenting area、tenting height、弁輪の面積、周囲長、高さ（最後部－最低部の差）など、臨床で必要とされる心臓弁に関する情報を得ることである。　心臓弁について臨床で必要なデータを計測可能な３次元心臓弁画像を得るため、心エコー図検査装置でスキャンした２次元の心エコー断層画像から３次元の心エコー図を形成し、前記３次元の心エコー図から３次元心臓弁画像をコンピュータ処理により自動抽出する方法において、心臓および弁輪の物理的形状を考慮したフィッティングモデルにおける前記弁輪のモデルのフィッティング評価関数（ポテンシャルエネルギー）をレプリカ交換法・拡張徐冷法で最適化することを特徴とする臨床で必要なデータを計測可能な３次元心臓弁画像自動抽出方法である。

Description

明細書

心臓弁データ計測方法および装置

技術分野

[0001] 本発明は、臨床に用いる心臓弁のデータを計測する方法及び装置に関する。更に詳細には、心臓弁の各種データを計測可能な鮮明な 3次元心臓弁画像を自動で抽出する方法および装置に関する。

背景技術

[0002] 僧帽弁逆流 (僧帽弁閉鎖不全）は弁膜疾患の中でも頻度が高ぐ重症の逆流例では左心不全を引き起こす重大な疾患である。重症僧帽弁逆流の治療は基本的に外科治療であり、従来は人工弁を用いた僧帽弁置換術が行われてきた。しかし、人工弁置換後には心機能の低下や抗凝固療法に伴う合併症などの問題が多ぐ最近自己弁を温存する僧帽弁形成術が広く行われるようになつてきた。

[0003] 僧帽弁形成術は僧帽弁輪 ·弁尖 ·腱索などのうち逆流の原因となっている部分を選択的に外科的形成する方法であり、その成功のためにも心エコー図検査による術前の正確な病因 ·病変診断が不可欠である。

[0004] しかし、現在広く用いられている心エコー図検查は 2次元画像による診断を行うものであり、 3次元的に複雑な構造をしている僧帽弁と僧帽弁周囲の解剖学的 ·位置的関係を明らかにすることは困難であった。つまり、鞍馬型 (サドル型）の曲線を描く弁輪や、精巧なカーブを有する弁尖 '弁葉、さらに腱索から乳頭筋'左心室に至る弁下部支持組織でなされる僧帽弁の機能的複合体 (僧帽弁機構)の立体的な構造を把握するためには、 2次元画像では不十分であり、 3次元的な画像診断が必要とされてきた。

[0005] 最近開発された 3次元心エコー図装置を用いることにより、非侵襲的かつ簡便に心臓全体をリアルタイムにスキャンし、画像を取り込むことが可能となった。 3次元心ェコ一画像では、外科医が心臓を見るのと同じように心筋 ·弁膜などの構造を観察することができ、これまでの 2次元画像による診断に比べて詳細な術前診断ができることが期待されている。 [0006] しかし、この 3次元画像を用いた 3次元解析 ·計測は未だ困難であり、実際の詳細な形態や位置関係を定量ィ匕することができないため、 3次元心エコー図は臨床的実用化されてレヽなレ、のが現状である。

非特許文献 1 :山田博三、「大局解抽出と局所解追跡の協調による心臓動画像認識」、 [online]、電子技術総合研究所彙報 62卷 7号、 [平成 17年 12月 22日検索]、インタ¹ ~ネット < URL:http://www.etl. go.jp/jp/results/bulletin/pdf/62_7/yamaaa72.pdf >

発明の開示

発明が解決しょうとする課題

[0007] CTや MRIは装置が大きぐ高価であるために、どこの病院でも利用できるというものではない。これに対して超音波を利用した心エコー診断装置は小型で、手軽に利用でき、広く普及している。このため心エコー（心臓超音波)検査は、心臓病や高血圧をはじめとする循環器疾患の診断、治療には欠かせない存在となっている。ひと昔前までは心臓カテーテル検查でしかわからなかった情報が、患者の痛みを伴うことなぐ瞬時に把握できるようになった。また装置の軽量化に伴い、レ、までは聴診器のように携帯し、出向先でも検査することが可能となった。

[0008] 心エコー装置から得られる心エコー図には、 Mモード心エコー図、断層心エコー図、ドッブラ心エコー図がある。 Mモード心エコー図では、心臓構造の経時的な動きをグラフィックに記録することができ、弁膜、心室壁、大動脈などの運動がそれぞれ特徴的なパターンとして描かれる。断層心エコー図では、超音波ビームを高速走査することによって心臓の 2次元断層像 (Bモード）が得られる。

[0009] 超音波ビームの走査方法として高速機械走査法と電子走査法がある。断層法には心臓の形態、または心臓の動きを容易に観察でき、異常所見の有無、部位、程度の診断に役立つ。ドッブラ心エコー図にはパルスドッブラ法、連続波ドッブラ法、ドッブラ断層法、 2次元血流映像法、カラードップラ法などがあり、狭窄流や弁逆流などの心腔内の異常血流の検査による定性診断とともに、血流計測や圧推定などの定量診断や心機能評価にも応用されている。

[0010] し力、し、心エコー図（echocardiogram)を用いれば、すべてが解決するというものではない。これまでの閾値によって判定するウィンドウ法、輝度変化の大きいところを抽出するエッジ抽出法などで行われてレ、る心臓弁自動抽出装置は誤認識が多レ、。ェコ一画像は CTや MRIなどの境界がはっきりした画像とは違レ、、境界が不明瞭である。ウィンドウ法やエッジ抽出法は、 CTや MRIなどのような境界がはっきりしたものには適用できるが、心エコー画像のような境界が不明瞭なものには適用できない。

[0011] 前記のウィンドウ法やエッジ抽出法の他に、心臓弁を曲線のフィッティングによりモデル化し、適当な最適化手法を用いて境界のはっきりしない画像から輪郭像を得る方法も提案されている。しかし、ニュートン法、最急勾配降下法では複雑な図形には適用できず、また自由度を大きくしても GA (遺伝的アルゴリズム）、 SA (徐冷法）などを用いても局所最小値に簡単に捕われてしまって、最適解を見つけることが難しいとレ、う問題点がある。

[0012] 以上の点を踏まえ、本発明が解決しょうとする課題は、心臓における僧帽弁の tenti ng volume, tenting area, tenting height,弁輪の面積、周囲長、高さ（最後部一最低部の差)など、臨床で必要とされる心臓弁に関する情報を得ることである。

[0013] 心エコー図検査装置で取り込んだ心エコー図をもとに、心臓弁を自動抽出して、鮮明な 3次元画像を得て、この画像から必要な量を計測する。本発明で用いる方法及び装置では、心エコー図における鮮明な 3次元画像の心臓弁輪を自動抽出することのほかに、エコー画像に写っていない組織の境界を再現することも可能となる。すなわち、これまで熟練した医師の目でしか認識できなかった心臓弁輪抽出と機能を、自動化する方法及び装置である。

課題を解決するための手段

[0014] 上記の発明が解決しょうとする課題を解決するために、請求項 1に記載された発明は、心臓弁について臨床で必要なデータを計測可能な 3次元心臓弁画像を得るため、心エコー図検查装置でスキャンした 2次元の心エコー断層画像から 3次元の心ェコ一図を形成し、前記 3次元の心エコー図から 3次元心臓弁画像をコンピュータ処理により自動抽出する方法において、心臓および弁輪の物理的形状を考慮したフイツティングモデルにおける前記弁輪のモデルのフィッティング評価関数 (ポテンシャルエネルギ一）をレプリカ交換法'拡張徐冷法で最適化することを特徴とする臨床で必要なデータを計測可能な 3次元心臓弁画像自動抽出方法である。

[0015] 請求項 2に記載された発明は、心臓弁について臨床で必要なデータを計測可能な 3次元心臓弁画像を得るため、心エコー図検査装置でスキャンした 2次元の心エコー断層画像から 3次元の心エコー図を形成し、前記 3次元の心エコー図から 3次元心臓弁画像をコンピュータ処理により自動抽出する装置において、心臓および弁輪の物理的形状を考慮したフィッティングモデルにおける前記弁輪のモデルのフイツティング評価関数 (ポテンシャルエネルギー）をレプリカ交換法 ·拡張徐冷法で最適化する手段を備えたことを特徴とする臨床で必要なデータを計測可能な 3次元心臓弁画像自動抽出装置である。

[0016] 本発明では、弁輪の抽出とそのフィッティングのために、具体的には以下のような手段をとること力 Sできる。弁輪抽出処理の流れは、以下の 2つの工程からなる。まず心臓の物理的形状を考慮したフィッティングモデルを用意して、心筋の輝度の高い部分をフィットさせる。次に、そのフィットした形状の上で、弁輪らしき場所を探す。

[0017] 弁輪のモデルは、弾性体でできた円筒状かつ網状の構造を用いることができる。たとえば、 1周 40個、高さ 40個、計 1600個の制御点を取り、各点間は適当な力のパネでつなぐ。このとき、制御点はなるべく輝度の高い位置に指定する。この円筒状の弁輪のモデルのフィッティング評価関数 (ポテンシャルエネルギー）をレプリカ交換法 '拡張徐冷法で最適化する。

[0018] レプリカ交換法 (RE)はタンパク質などの立体分子構造を解明する手法として広く用いられている。この手法では、等価な相互作用のない系（レプリカ）複数個で構成される全体系を考え、各レプリカ（コピー）には異なる温度（エネルギー）を割り当て、最初はすべてのレプリカに同じ分子を配置する。各レプリカ系で独立にメトロポリスシミユレーシヨンを行い、定期的に隣接するレプリカ間で分子配置を交換する。

[0019] さらに徐冷法（SA)を用レ、、高温（高工ネルギー）から低温 (低エネルギー）に徐冷して最適解を見つける。この方法では、構造の、エネルギーのポテンシャル面が最小 (または極小）になる点 (最適解）を見つけ、安定構造として最終的な分子構造を決定している。

[0020] レプリカ間分子配置交換には、まったくランダム交換するモンテカルロ法や、遺伝子組み替えのように距離的に近い分子間（隣接分子間)で交換を行う遺伝子的アルゴリズム (GA)などが用いられる。分子構造をモデル化するに際しては、分子は点と考え、隣接する分子間はクーロン力、パネ相互作用などが働き、これらの力の和が分子内ポテンシャルエネルギーとして表現されている。

[0021] 弁輪探索は以下のルールに従って行う。

•なるべく輝度の高いところ

•下から上の方を見て 2次微分が正の場所 (窪んでレ、る所）

•隣の弁輪の制御点とあまり離れないところ

このルールにも適当な評価関数を定義して同様に最適化を行い、ポテンシャルエネルギ一が最小となる構造が、最終的に弁輪らしき場所として抽出される。弁輪の自動抽出力まくいかない場合に備え、マニュアルにて修正することのできる道も用意することちできる。

発明の効果

[0022] 本発明の装置の特徴は、以下のとおりである。（1)従来の 2次元心エコー図では不可能であった、僧帽弁複合体の 3次元表示および定量解析ができること。（2)従来の 2 次元画像からの再構築には多大な労力と時間を要していたが、本発明の装置を用いた場合、現時点で僧帽弁 3次元解析にいたるすべてのプロセス（エコー画像の収集- 画像のトレース · 3次元画像再構築 · 3次元データ定量解析）に要する時間を計約 15 分に短縮できる。

[0023] 3次元心エコー図を用いた 3次元的定量解析は未だ実現されておらず、本発明を用いた研究成果は世界的にも初めてのものである。とくに現在、僧帽弁尖'弁葉に異常がないにもかかわらず乳頭筋や左心室の機能不全によって引き起こされる「機能的僧帽弁逆流」のメカニズム解明や治療法の開発が世界的に注目されており、これまで 2次元心エコー図画像での解析に頼っていたこれらのテーマに対する研究は、 3 次元解析が可能となることで大きく飛躍することが予想される。本発明の装置は、僧帽弁逆流の術前診断'外科治療にまで密接した、臨床的に大いに有意義なものである。

図面の簡単な説明 [0024] [図 1]心臓収縮時の心臓画像を心エコー図診断装置を用いて等間隔で 18枚の画像を撮影し、僧帽弁のリーフレットと弁輪を 3-Dイメージに再構築する手法を説明するための図である。

[図 2]2-D画像から 3-Dイメージ画像を作成するときの説明図である。

[図 3]本発明の装置により 3次元心臓弁画像を再構築した、正常な人の僧帽弁のリーフレットと弁輪の、（A)外観を示す斜視図、（B) LV方向から見た僧帽弁のリーフレツト（上面図）および側面図、（C)前図を修正した僧帽弁のリーフレットの上面図と側面図である。

[図 4]本発明の装置により 3次元心臓弁画像を再構築した、虚血性 MRの人の僧帽弁のリーフレットと弁輪の、（A)外観を示す斜視図、（B) LV方向から見た僧帽弁のリーフレット（上面図）および側面図、（C)前図を修正した僧帽弁のリーフレット上面図と側面図である。

[図 5]僧帽弁のリーフレット上に、 12人の局所貧血 MRを持った患者の、最大テンティングサイトの人数分布を示した図である。

[図 6]制御点と弾性パネの説明図である。

[図 7]評価関数と積分領域の説明図である。

[図 8]ポテンシャル面と局所最小値の説明図である。

発明を実施するための最良の形態

[0025] 本発明を実施するための最良の形態を図を用いて説明する。以下の説明では次のような省略リストを用いる。

MR =僧帽弁逆流

3-D = 3次元の

2-D = 2次元の

LV =左心室

LA =左心房

ROA =逆流開口部面積

EF =駆出分画

PISA =近位部等流速面積 EDV =拡張末期容量 (容積）

ESV =収縮末期容量 (容積）

[0026] 2次元の心エコー図検査では、以下のような検査が行える。すべての被験者の標準的な 2-D心エコー図検査を行うことにより、拡張末期容量 (EDV)と収縮末期容量 (ES V)がモディファイドシンプソン方法 (左心室全体を円筒形の積み重ねとして近似する )によって測れる。その結果駆出分画（％)が方程式 100 X (EDV - ESV)/EDVによつて計算できる。 MRはカラードップラ心エコー図検查によって評価し、 MRの度合いが ROAを使っている PISA方法によって数量化することができる。しかし、僧帽弁逆流診断や僧帽弁手術を行うには、正確な弁輪の位置を見分け、立体的な弁膜画像が必要になる。本発明の装置では僧帽弁を以下の手順の実行により抽出し、鮮明な画像を再現する。

[0027] 立体画像イメージ (volumetic image)を得るために、リアルタイム 3-D心エコー図システムを利用して、被験者の心尖部ビューで胸腔を通しての容積測定のイメージ (フルのボリュームモード）を得る。容積測定のフレームレートは 12から 16センチメートノレの深度で、毎秒 16から 22フレームを撮影する。（このフレーム数は深度に依存する）フルの立体画像イメージを獲得する前に、プローブが 2D画像の様式で僧帽弁の中心部の頂点に位置するように調整する。すべての立体画像イメージはディジタル方式でコンパクトディスク上に記録し、オフライン分析のためにパーソナルコンピュータへ転送する。

[0028] 図 1は、 3次元心エコー図診断装置を用い、 18枚の等間隔の放射状平面画像を自動撮影 (スキャン）し、その平面画像をもとに 3次元（3-D)の立体画像イメージを形成する過程を示している。心臓収縮時に、スキャンされた各平面画像に対して、手作業で僧帽弁の弁輪とリーフレットにマークを付ける。これらのデータから、僧帽弁のリーフレットと弁輪の 3-Dイメージが再構築される。具体的には、図 2に示す。

[0029] 図 2はスキャンされた各平面画像を順番に並べ、チェック点（画像間で対応する点）を線でつなぎ、スムージングをかけ、レンダリングを行うと滑らかな対応物の 3次元立体画像が得られる（図の例は 18フレームの場合)。し力背景技術でも述べたように、エコー画像は MRIや CTのような鮮明な輪郭が得られない。とくに弁輪は複雑であり、入り組んでいるために、本発明の装置ではフィッティングモデルを使って弁輪の抽出処理を行う。そのために、心臓の物理的形状を考慮して、心筋の輝度の高い部分をフィットさせる。さらに、フィットした形状の上で、弁輪らしき場所を探す。この弁輪探

[0030] この例では、フィッティングモデルには、弾性体で出来た円筒状かつ網状の構造を用いる。さらに 1周 40個、高さ 40個の、計 1600個の制御点を設け、各制御点を適当なパネでつなぐ。このとき、制御点はなるべく輝度の高い所が選ばれるようにしてある。このような構造物（レプリカ）を複数用意し、それぞれのレプリカには異なる輝度を割り当てる。この輝度をポテンシャルエネルギーとし、フィッティング評価関数を用いてポテンシャルエネルギーが最も安定した（最小値となる）構造を決定する。この際に用いる方法が、レプリカ交換法と呼ばれる拡張徐冷法である。すなわち、輝度の高いところから始め、レプリカ間で制御点の交換を行レ、、そのつどポテンシャルエネルギーを求める。このようなシミュレーションを行うことで、ポテンシャルエネルギーが最小（安定構造)となるものが、構造の形状 (弁輪)として抽出される (最適化)。

[0031] なお、弁輪自動抽出探索を行う際の主なルールは以下のとおりである。すなわち、 •探索はなるべく輝度の高いところ

•隣の弁輪の制御点とあまり離れていない所

このルールにも適当な評価関数を定義して、心臓の形状探索と同様に最適化を行うと、最終的に弁輪らしき場所が抽出される。なお、弁輪の自動抽出がうまく行かなかつたり、あるいは曖昧な点があるような場合には、マニュアル修正のルートも設けてお <。

[0032] 以上のようにして抽出された弁輪が図 3である。なお以下では、抽出された僧帽弁膜を「リーフレット（leaflet :弁葉）」と呼び、リーフレットの付け根を「弁輪」と呼ぶ。また、リーフレットが天幕を張ったように膨らみをもつこと、またはその膨らみを「テンティング (tenting)」と呼ぶことにする。肺で綺麗になった血液は左心房（LA)に流れ込み、 L A力僧帽弁を通して左心室 (LV)に送られ、さらに LVから大動脈を通して体全体に送出される。このため、 LAの圧力よりも LVの圧力の方が高くなる。 [0033] 弁に物理的あるいは機能的な障害が生じると、狭窄症や虚血症を引き起こすことになる。例えば、僧帽弁の開きが悪ぐ LVに十分に血液が送られないと、狭窄症が起きる。逆に僧帽弁が緩み、僧帽弁閉鎖不全になると、 LVから LAの血液の逆流が起き、十分な動脈血が体に供給されないために、虚血が起きる。最近は、僧帽弁の弛緩などによる僧帽弁閉鎖不全は、人工弁を使わずに弁形成手術で治療することが多くなっている。そのためには僧帽弁の確力、な形状を得ることが重要であり、本発明の装置の 3次元心臓弁画像取得が有効となる。

[0034] なお、図に示した記号は以下の意味を表している。

A 冃 ij尖部 (anterior)

P 後尖部（posterior)

しし父連 (antero- lateral commissure)

し M 後父連、 postero-medial commissure)

LV 左心室（left ventricle)

LA 左、房 (left atnum)

an皿 lar height 弁輪の高さ（湾曲度合い）

tenting length テンティングの長さ（テンティング長）

[0035] 図 3は、本発明の装置で行われる 3次元心臓弁画像抽出による、正常な人の僧帽弁のリーフレットの立体画像とその形状を示している。図の（A)は、リーフレットを異なる方向から見た 3-Dイメージである。ここでは、僧帽弁の弁輪 (付け根部）が「サドノレ形」をしている。僧帽弁のリーフレットが LVの中にわずかに食い込んでいる力ほとんど平らに見える。

[0036] 図の（B)は、実際の 3-Dテンティングイメージである。弁輪は 3-Dで測定するための、およその外観である。左図は、 LV方向から見た僧帽弁のリーフレット形状であり、テンティングの度合いを等高線として表している。右図は水平方向から見たリーフレツトで、弁輪とリーフレットのテンティング度合いを正確に測ることができる。また、僧帽弁の弁輪の円周やエリアが、これらの 3-Dデータから測れる。図の弁輪の高さが、弁輪の湾曲度合レ、を示してレ、る。図の黒レ、点は接合ラインである（弁の交連部すなわち LV収縮時の前尖と後尖のかみ合う位置)。正常な場合、 LVの収縮期には僧帽弁の腱索に支えられて、前尖と後尖はきちんと合わさり、左心室から左心房への血流は遮断される。この接合ラインが黒点である。

[0037] 図の（C)は、 3-Dテンティングイメージを修正したものである。図の太線で示した曲線は僧帽弁の弁輪を表し、輪状の表面からリーフレットまでの距離を一定に保ちながら、平面上で滑らかに描いたものである。左図は LV方向から見た僧帽弁のリーフレツト図であり、テンティングの度合いを等高線で表している。右図は水平方向から見た図で、僧帽弁の弁輪のテンティング度合いを定量的に測定することができる。最大テンティング長、平均テンティング長およびテンティングの容積も、これらの 3-Dデータ力測定することができる。なお、黒い点は接合ラインを示す。

[0038] 図 4は、僧帽弁の弁輪の立体画像で、虚血性の僧帽弁逆流（MR)患者の僧帽弁を表すリーフレットである。図の（A)は、リーフレットを異なる方向から見た 3-Dイメージである。僧帽弁の弁輪が外見上も明らかにテンティングによって、平滑化している。また、山の形をしたリーフレットが凸状になり、全体的に LVへと食い込んでいる。

[0039] 図の（B)は、実際の 3-Dテンティングイメージを表している。僧帽弁のリーフレット全体が明らかに LV方向に膨らみ、弁輪の高さも正常と比較してみても低いことがわかる。また僧帽弁の弁輪も広がっている。なお、黒い点は僧帽弁の接合ラインを示している。

[0040] 図の（C)は、 3-Dテンティングイメージを修正したものである。左図からもわかるように、僧帽弁のリーフレットが、僧帽弁の弁輪で見た場合、 A-Pに対してほとんど対称的になっている。右図力らもわかるように、最大テンティング長は正常なものよりも長くなっている。黒い点は接合ラインを示している。カラー表示した場合、緑のマーク (右図の薄く印刷された部分）がリーフレットの最大テンティングサイトを示している。この患者の場合、最大テンティングサイトはリーフレット前尖部 Aの中央に位置している（左図の等高線が最も高い位置を示している場所 <山の頂上に相当した位置 >)。

[0041] 図 5は、局所貧血 MRを持っている 12人の患者について最大テンティングサイトの場所を調べた結果を、リーフレット上に人数分布として表した図である。図の英字 'A' は前尖部（anterior)、 'P'は後尖部（posterior)を表し、 'L'は側部（lateral)、 'C，は中央部（central)、 'M'は中間部（medial)を表す。また、括弧内の数字は人数である。図が示すように、最大のテンティングサイトは、 12人すべてにおいてリーフレット前方部に位置していた。内訳は AMが 3人、 ACが 5人、 ALが 4人であった。

[0042] 虚血性の MRを持っている 12人の患者の内訳は、単一の脈管疾患を持つ 3人の患者、二つの脈管疾患を持つ 6人の患者、三つの脈管疾患を持つ 3人の患者である。

LV機能障害は広範囲（EF33.9±9.1%、幅 18%〜47%)でひどかった。 R〇Aは 0.29 ± 0 .15cm²であって、 0.15〜0.62cm²に及んだ。 10件の実験照查基準と比較して、虚血性の MRを持っている患者が年齢、性別、あるいは体表面エリアにおける相違は見いだせな力、つた。しかし LVの容積力正常な人と比較して虚血性の MRの患者では際立つて増カロしていた。

[0043] 以上のように、本発明の 3次元心臓弁画像取得方法を利用して開発された、リアルタイム 3-D心エコー図検查用イメージ作成ソフトウェアシステムは、 1)僧帽弁のリーフレットと弁輪の 3-D幾何学的な奇形、 2)僧帽弁のリーフレットの最大テンティングサイト、および 3)虚血性の MRを持っている患者での僧帽弁テンティングと弁輪の奇形の量的な測定をすることができた。

[0044] 本発明で用いるフィッティングモデルについて、更に詳細な例を説明する。エコー測定装置の特性により、画像にノイズや影が発生し、その画像の情報だけから正確な組織の像を得ることは難しレ、。医師は実際の組織の理想的なイメージを知っており、その理想的イメージといろいろな角度や時間のエコー画像を頭の中で組み合わせ、不鮮明なエコー画像を補完して組織の境界線を引くという作業を行っている。物理モデリングはこの医師の頭の中で行っている画像補完作業をコンピューターに肩代わりさせるものである。

[0045] コンピューター上でのモデルの構築は、図 6および図 7に示すような制御点間を繋ぐパネと制御点によって張られる境界評価関数を用いる。制御点間のパネは組織の物理的な構造を維持する。一方、境界評価関数は画像から組織の境界情報を拾う。

[0046] モデルの評価には、パネの弾性エネルギーと境界評価による評価エネルギーを足し合わせたポテンシャル関数を用いる。 i番目の制御点の位置を rとして、制御点のセット r，r，· · · Γを r^Nのように書く。このとき、パネの弾性エネルギー関数 S(r^N)を以下のよ

1 2 N

うに定義する。 [0047] [数 1]

[0048] ここではパネの弾性強度で、制御点間の組織強度などから経験的に決定される。

繋がっていない制御点間の場合は 0とする。 σは制御点の排除半径で、制御点同土が重ならないように選ぶ。このとき、パネの自然長は

[0049] [数 2]

[0050] となる。理想的な物理形状の時にエネルギーが最低になるようにこれらのパラメータ一を設定する。評価エネルギー関数

[0051] [数 3]

E{ r^N)

[0052] は、画像のベクトル r地点の輝度を返す関数

[0053] [数 4]

M{r)

[0054] を用いて以下のように定義する。

[0055] [数 5]

[0056] ここで

[0057] [数 6] [0058] は、

[0059] [数 7] = ( 一一 I

[0060] で定義される長さ 1のベクトル、 cは結合の定数、

ij

[0061] [数 8]

[0062] は制御点 i， j間の評価関数 (後述）、線積分は rと rの間の最短ルートで定義される

[0063] 各制御点間の評価関数は、物理的特性を反映した関数を経験的に選ぶ。例えば、輝度値が大きく変化していることを境界として認識する（境界付近でエネルギーが下がる）ような関数は以下のように書ける。

[0064] [数 9]

[0065] ここで、

[0066] [数 10]

[0067] は、

[0068] [数 11] [0069] に垂直な長さ 1のベクトル、

[0070] [数 12]

[0071] は関数 Mの勾配である。

[0072] 同様に、以下のような関数を評価関数として用いることが出来る。

[0073] [数 13] 空洞： ,_void(A ,p,?)=-M( »)

組織： f^(M,p,f)=M(p)

評価関数なし： )-0

[0074] これらの関数を組み合わせて、最終的に以下のようなエネルギー関数を定義する。

[0075] [数 14]

F(r^N)=W_s-S(r^N)+W_E'E(r^N)

[0076] W ,Wは弾性エネルギーと評価エネルギーの重みで、構造と境界評価のどちらを重視したレ、かを調節できる。この関数 Fが最低になるような制御点のセット

[0077] [数 15]

[0078] を探すことで、組織の境界を抽出することが出来る。この関数 Fの特徴は、関数 S による物理形状を維持したまま、関数 Eによって境界を探すことができることである。エコー画像に多少のノイズや影があっても物理形状による補完が行われるため、画像だけから形状が推測できない場合でも、もっともらしい組織の境界を抽出することが出来る。

[0079] 拡張徐冷法による最適化処理について説明する。組織の境界を正確に求めるためには多数の制御点が必要であるが、評価関数は非線形のため、制御点の数が多くなると評価関数の最小点を見つけることは難しい。一般的な最適化手法である二ユートン法、最急勾配降下法、 GA (遺伝的アルゴリズム）、 SA (徐冷法）などを用いても局所最小値に簡単に捕われてしまって、最適解を見つけることが出来ない。

[0080] 拡張徐冷法は、スピングラスやタンパク質のフォールデイング問題など、物理'化学の分野で近年注目されている強力な最適化手法で、多自由度の複雑な最適化問題を効率良く解くことの出来る計算手法である。我々は、この拡張徐冷法の一種であるレプリカ交換モンテカルロ法を用いて評価関数を最適化した。

[0081] まず、レプリカ交換モンテカルロ法の基礎になる、モンテカルロ法による徐冷法を説明する。モンテカルロ法は確率的アルゴリズムを使って計算機シミュレーションを行う手法である。制御点のセット

[0082] [数 16]

' \, ^r2, " ' ^rN

[0083] について、最初のステップの制御点のセットを

[0084] [数 17]

N

' 0

[0085] とする。次に、制御点の中からランダムに 1つ選んで、ランダムな方向'大きさ

[0086] [数 18]

[0087] だけずらす。その制御点のセットを

[0088] [数 19]

, ,Ν

r 0

[0089] と書くと、最初のステップとずらした時の評価エネルギーはそれぞれ [0090] [数 20]

E(r!) 、 E(r'S)

[0091] と書ける _c のずらした制御点のセットは以下の確率で次のステップに採用される _c [0092] [数 21]

[0093] の時、採用 _c

[0094] [数 22]

E(^)<E(r^)

[0095] の時

[0096] [数 23]

[0097] の確率で採用。

[0098] 採用された場合は、

[0099] [数 24]

[0100] として次のステップに進む。棄却された場合は

[0101] [数 25]

[0102] として次のステップに進む。

[0103] ここで βはシステムの最適化具合を決定するパラメーターで、統計力学ではボルツマン定数を k、温度を Tとしたとき = l/k Tとして知られるパラメーターである。 β

Β Β

が十分大きい（温度が低い）とき、シミュレーションが進んでいくと評価関数は小さい方向のみに進んでいく。一方、が小さい（温度が高い）とき、評価関数は大きくなる方向にも動けるようになるため、大きな変動を示すようになる。

[0104] 評価関数は 3N+1次元空間の面（ポテンシャル面）として表現でき、シミュレーションはその表面の極小点を飛び移りながら進んでいく。

[0105] このモンテカルロ法による最適化探索を行うには、まず j3を小さくして制御点のセットを十分ランダムな状態にして力 ^混ぜた後、 βを徐々に大きくして評価関数の値を収束させていく。十分 j3を大きくしたところでしばらくシミュレーションを走らせて、評価関数が最小になる制御点のセット

[0106] [数 26]

N

0

[0107] を探す。制御点の数が少ない場合や、評価関数が複雑でない場合はこの方法で最適化点を見つけることが出来る。しかし、一般的にはポテンシャル面は複雑な形状をしており、温度を単純に下げただけではすぐに局所最小値にトラップされてしまい、レ、くら時間をかけても十分な最適化を行うことが出来ない。（図 8参照）

[0108] レプリカ交換モンテカルロ法は、上で説明したモンテカルロシミュレーションを異なる温度で並列に走らせることで局所最小値へのトラップを防ぎ、効率良く最適化を行う手法である。 M個の制御点のセット（レプリカ）を用意し、 m番目の制御点のセット

[0109] [数 27]

-JV

' M

[0110] をパラメーター β でシミュレーションを進める。ここで、各レプリカの温度パラメータ

m

一は β く β のように温度の高い順に並んでいるものとする。適当なステップごとに m m+1

、以下の方法でレプリカ間の制御点の位置を入れ替える。

[0111] [数 28] [0112] を計算し、

Δ < 0の時、入れ替え

Δ >0の時、 exp (— Δ )の確率で入れ替え

[0113] この温度交換の方法により、温度の低い制御点のセットが局所最小値にトラップされていても、適度に温度の高いセットと入れ替えられ、局所最小値を脱出することが出来る。ノメーターを適切にセットすれば、時間をかければかけるほど最適化が進む。

[0114] 一般的にレプリカの数を増やし、隣り合う /3 の差を狭くすると最適化解を見つけや m

すくなるが、レプリカの数を増やすと計算コストがかさむため、評価関数の分散を調べながら計算効率が最高になるように調節する必要がある。計算効率を上げるためには、レプリカ間で十分な頻度で交換が起こり、 1つのレプリカが温度空間でランダムゥオーク出来るような状況が必要である。そのためには、隣り合うレプリカの評価関数の分散が同じ面積で重なり合うように、レプリカの数やの値を調節する必要がある。

m

[0115] この方法はパラメーターの調節や、計算アルゴリズムの実装に非常に手間がかかるが、時間さえかければかなりの確率で最適化点を見つけることができる。そのため、一旦パラメーターを調節すれば、他の最適化アルゴリズムのように初期値を注意深く選んで何度も試行する必要がなレ、。そのため、人間がほとんど手をカ卩えずに自動的に精度の高い境界線を抽出することが出来る。

Claims

請求の範囲

[1] 心臓弁について臨床で必要なデータを計測可能な 3次元心臓弁画像を得るため、心エコー図検查装置でスキャンした 2次元の心エコー断層画像から 3次元の心ェコ一図を形成し、前記 3次元の心エコー図から 3次元心臓弁画像をコンピュータ処理により自動抽出する方法において、

心臓および弁輪の物理的形状を考慮したフィッティングモデルにおける前記弁輪のモデルのフィッティング評価関数 (ポテンシャルエネルギー）をレプリカ交換法 '拡張徐冷法で最適化することを特徴とする臨床で必要なデータを計測可能な 3次元心臓弁画像自動抽出方法。

[2] 心臓弁について臨床で必要なデータを計測可能な 3次元心臓弁画像を得るため、心エコー図検査装置でスキャンした 2次元の心エコー断層画像から 3次元の心ェコ一図を形成し、前記 3次元の心エコー図から 3次元心臓弁画像をコンピュータ処理により自動抽出する装置において、

心臓および弁輪の物理的形状を考慮したフィッティングモデルにおける前記弁輪のモデルのフィッティング評価関数 (ポテンシャルエネルギー）をレプリカ交換法 '拡張徐冷法で最適化する手段を備えたことを特徴とする臨床で必要なデータを計測可能な 3次元心臓弁画像自動抽出装置。