▶ ユニヴェルシテ・ドゥ・ボルドーの特許一覧 ▶ アンスティチュ ポリテクニーク ドゥ ボルドーの特許一覧 ▶ サントル・オスピタリエ・ユニベルシテール・ドゥ・ボルドーの特許一覧 ▶ サントル ナシオナル ドゥ ラ ルシェルシェ シアンティフィクの特許一覧 ▶ アンスティチュ・ナシオナル・ドゥ・ルシェルシュ・アン・アンフォルマティク・エ・アン・オートマティクの特許一覧
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2024-02-06
(45)【発行日】2024-02-15
(54)【発明の名称】定位脳標的を判定する方法
(51)【国際特許分類】
A61N 1/36 20060101AFI20240207BHJP
A61N 5/10 20060101ALI20240207BHJP
A61N 7/02 20060101ALI20240207BHJP
【ＦＩ】
A61N1/36
A61N5/10 P
A61N7/02
【請求項の数】
14
(21)【出願番号】P 2021521382
(86)(22)【出願日】2019-11-18
(65)【公表番号】
(43)【公表日】2022-02-07
(86)【国際出願番号】 EP2019081571
(87)【国際公開番号】W WO2020099671
(87)【国際公開日】2020-05-22
【審査請求日】2022-08-15
(31)【優先権主張番号】1860588
(32)【優先日】2018-11-16
(33)【優先権主張国・地域又は機関】FR
(73)【特許権者】
【識別番号】514058706
【氏名又は名称】ユニヴェルシテ・ドゥ・ボルドー
(73)【特許権者】
【識別番号】513277212
【氏名又は名称】アンスティチュ ポリテクニーク ドゥ ボルドー
(73)【特許権者】
【識別番号】516324917
【氏名又は名称】サントル・オスピタリエ・ユニベルシテール・ドゥ・ボルドー
(73)【特許権者】
【識別番号】311016455
【氏名又は名称】サントル ナシオナル ドゥ ラ ルシェルシェ シアンティフィク
【氏名又は名称原語表記】ＣＥＮＴＲＥ ＮＡＴＩＯＮＡＬ ＤＥ ＬＡ ＲＥＣＨＥＲＣＨＥ ＳＣＩＥＮＴＩＦＩＱＵＥ
(73)【特許権者】
【識別番号】516308814
【氏名又は名称】アンスティチュ・ナシオナル・ドゥ・ルシェルシュ・アン・アンフォルマティク・エ・アン・オートマティク
(74)【代理人】
【識別番号】110001173
【氏名又は名称】弁理士法人川口國際特許事務所
(72)【発明者】
【氏名】エンゲルハート，ジュリアン
(72)【発明者】
【氏名】ゼムゼミ，ナジブ
(72)【発明者】
【氏名】キュニー，エマニュエル
【審査官】木村 立人
(56)【参考文献】
【文献】特表２０１２－５１０８７７（ＪＰ，Ａ）
【文献】KIM, Jinyoung et al.，“Automatic localization of the subthalamic nucleus on patient-specific clinical MRI by incorporating 7T MRI and machine learning: Application in deep brain stimulation”，bioRxiv [online]，Cold Spring Harbor Laboratory，2018年05月15日，[検索日 2023年10月05日]，Retrieved from <https://www.biorxiv.org/content/10.1101/322230v1.full.pdf>
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ａ６１Ｂ ５／０５５
Ａ６１Ｂ ６／００ ― ８／１５
Ａ６１Ｂ １７／００ ― ９０／９８
Ａ６１Ｎ １／００ ― ７／０２
(57)【特許請求の範囲】
【請求項1】
少なくとも１つの標的点（ＰＣ’）を含む定位脳標的を判定する方法であって、所与の症状についての前記標的における神経外科治療前にプロセッサで実施され、及び－ 少なくとも１つの標的点（ＰＣ）で前記症状について実行される前記治療後に測定される結果が閾値以上である患者の複数（ｎ）の臨床例を選択するステップであって、術後画像化は、前記患者の各々について実行されている、ステップ、
－ 数学的座標系、好ましくは正規直交デカルト座標系を選択するステップ、
－ それぞれの選択された臨床例について、前記選択された座標系における前記少なくとも１つの標的点（ＰＣ）の（ｐ）個の座標（ｘＰＣ、ｙＰＣ、ｚＰＣ、ｖｘＰＣ、ｖｙＰＣ、ｖｚＰＣ）の全部又は一部を判定するために前記術後画像化を処理するステップ、
－ 複数（ｍ）の脳マーカー点（ＰＲ）を選択するステップ、
－ それぞれの選択された臨床例について、前記脳マーカー点（ＰＲ）の（３×ｍ）個の座標（ｘＰＲ、ｙＰＲ、ｚＰＲ）の全部又は一部を判定するために前記術後画像化を処理するステップ、
－ 前記（ｎ）個の選択された臨床例の全部の前記標的点（ＰＣ）の前記判定された座標及び前記マーカー点（ＰＲ）の前記判定された座標を含む学習データベースを生成するステップ、
－ 前記学習データベース及び教師あり統計学習方法を使用することにより、前記マーカー点（ＰＲ）に従って少なくとも１つの標的点（ＰＣ）の前記座標を与える予測関数（Ｆ）を判定するステップ、
－ 前記症状について治療される新しい患者の術前画像化を、前記新しい患者の前記マーカー点（ＰＲ’）の座標（ｘＰＲ’、ｙＰＲ’、ｚＰＲ’）の全部又は一部を判定するために処理するステップ、
－ 前記新しい患者について判定される前記マーカー点（ＰＲ’）の前記座標に従い、前記新しい患者について少なくとも１つの標的点（ＰＣ’）の座標を得るために前記予測関数（Ｆ）を使用するステップ
を含む方法。
【請求項2】
－ 交差検証方法を使用して前記予測関数（Ｆ）を統合するステップであって、前記マーカー点（ＰＲ）に従って少なくとも１つの標的点（ＰＣ）の前記座標を与える統合予測関数（ＦＣ）をもたらすステップを更に含み、前記予測関数（Ｆ）を使用する前記ステップは、前記統合予測関数（ＦＣ）を使用するステップを含む、請求項１に記載の方法。
【請求項3】
前記マーカー点（ＰＲ、ＰＲ’）及び／又は前記少なくとも１つの標的点（ＰＣ）の前記座標を判定するために処理される前記画像化は、少なくとも１つのＭＲＩ画像、好ましくは複数のＭＲＩ画像である、請求項１又は２に記載の方法。
【請求項4】
前記教師あり統計学習方法は、再生カーネルヒルベルト空間におけるカーネルリッジ回帰方法の使用を含む、請求項１～３のいずれか一項に記載の方法。
【請求項5】
前記教師あり統計学習方法は、サポートベクトルマシン型の方法の使用を含む、請求項１～３のいずれか一項に記載の方法。
【請求項6】
前記教師あり統計学習方法は、ニューラルネットワーク型の方法の使用を含む、請求項１～３のいずれか一項に記載の方法。
【請求項7】
前記交差検証方法は、「リーブワンアウト交差検証」方法の使用を含む、請求項２～６のいずれか一項に記載の方法。
【請求項8】
前記交差検証方法は、「リーブｋアウト交差検証」方法の使用を含む、請求項２～６のいずれか一項に記載の方法。
【請求項9】
前記数学的座標系は、正規直交デカルト座標系であり、直線（ＣＡＣＰ）は、Ｏｙ軸を形成する前交連（ＣＡ）の上及び後縁並びに後交連（ＣＰ）の下及び前縁を通り、前記後交連（ＣＰ）の前記下及び前縁は、前記座標系の中心を形成し、及びＯｚ軸は、半球間平面において前記直線（ＣＡＣＰ）に垂直な直線である、請求項１～８のいずれか一項に記載の方法。
【請求項10】
前記マーカー点（ＰＲ）は、下記の１８個の点：－ 第３の軸平面上の乳頭視床束である第１のマーカー点（ＦＭＴ）、
－ 第１、第２及び第３の軸平面の各々の上の被殻の前点である第２、第３及び第４のマーカー点（ＰＡ１、ＰＡ２、ＰＡ３）、
－ 前記第１及び第２の軸平面上の前記被殻の内側点である第５及び第６のマーカー点（ＰＭ１、ＰＭ２）、
－ 前記第１及び第２の軸平面上の前記被殻の後点である第７及び第８のマーカー点（ＰＰ１、ＰＰ２）、
－ 前記第２の軸平面上の手綱交連である第９のマーカー点（ＣＨ）、
－ 前記第２の軸平面上の視床の前縁である第１０のマーカー点（ＢＡＴ）、
－ 前記第２の軸平面上の前記視床の後縁である第１１のマーカー点（ＢＰＴ）、
－ 前記前交連である第１２のマーカー点（ＣＡ）、
－ 中交連点における第３脳室の内側縁である第１３のマーカー点（Ａ）、
－ 前記第１３のマーカー点（Ａ）を通る矢状面上の前記視床の高さである第１４のマーカー点（Ｂ）、
－ 前記第１３及び第１４のマーカー点によって定義される線分［ＡＢ］の中点である第１５のマーカー点（Ｃ）、
－ 前記前交連の前記上及び後縁並びに前記後交連の前記下及び前縁を通る前記直線（ＣＡＣＰ）に平行であり、且つ前記第１５のマーカー点（Ｃ）を通る直線上の前記視床の前記前縁である第１６のマーカー点（Ｄ）、
－ 前記第５のマーカー点（ＰＭ１）を通る前頭面上の前記被殻の上縁である第１７のマーカー点（Ｐｃｕｌｍ）、及び
－ 前記第５のマーカー点（ＰＭ１）を通る前記前頭面上の前記被殻の外側縁である第１８のマーカー点（Ｐｌａｔ）
の中から選択される、請求項１～９のいずれか一項に記載の方法。
【請求項11】
前記学習データベースに機能データを追加する追加ステップを更に含み、前記機能データは、前記臨床例の前記標的点（ＰＣ）及び前記マーカー点（ＰＲ）における信頼度の少なくとも１つの指標を追加することができる、請求項１～１０のいずれか一項に記載の方法。
【請求項12】
深部脳刺激、ガンマナイフ又は集束超音波治療前にプロセッサで実施される、請求項１～１１のいずれか一項に記載の方法。
【請求項13】
請求項１～１２のいずれか一項に記載の方法の前記ステップの全部を実施するように構成されたプロセッサを含むデータ処理システム。
【請求項14】
コンピュータプログラムであって、前記プログラムがプロセッサによって実行されると、前記プロセッサに、請求項１～１２のいずれか一項に記載の方法の前記ステップの全部を実施させる命令を含むコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明の分野は、神経外科の分野であり、本発明は、より詳細には、神経外科的介入前の準備のための方法に関する。
【0002】
本発明は、より詳細には、定位脳標的を正確に判定する方法に関する。
【0003】
本発明は、特に、深部脳刺激のための電極の埋め込み前に準備ステップとして刺激される脳標的を判定する方法に関し得る。
【背景技術】
【0004】
脳標的の判定は、任意の神経外科治療手順、特に定位固定において重要である。定位固定は、脳内に位置する点に基づいて、神経外科医が介入する必要がある脳の部位又は病巣の体積及び位置を正確に判定するために、医学画像化によって支援される頭蓋内構造体を三次元的に同定する方法を使用する神経外科技法である。
【0005】
１つの介入は、深部脳刺激（ＤＢＳ）であり得る。深部脳刺激は、神経学的病状（例えば、パーキンソン病、本態性振戦病、ジストニー、強迫性障害、ツレット症候群、不応性てんかん又は更に治療の効かない重度のうつ病）を有する患者を軽減する非常に効果的な外科技法である。
【0006】
深部脳刺激は、脳の奥深くにある過剰な構造体を刺激することを含む。刺激は、ニューロンを抑制又は活性化し得、目的は、ニューラルネットワークの機能を修正することである。刺激は、症状を軽減又は更に解消して、通常の治療に対して抵抗性がある患者の生活の質を向上させ得る。
【0007】
深部脳刺激は、一般的に、核、一般的に、通常、灰白質（神経細胞体）における核の一部に対応する脳構造体に電極を埋め込むことによって実行される。残りの説明の全体にわたり、用語「核」は、核又は核の一部を意味し得る。
【0008】
電極は、このように標的核に送出される低強度の電流を受け取る。
【0009】
深部脳刺激で標的にされる脳構造体（パーキンソン病を治療するための視床下核（ＳＴＮ）、本態性振戦を治療するための視床の中間腹側核（ＶＩＭ）、ジストニー及びある種のパーキンソン病を治療するための内部淡蒼球（ＧＰｉ））は、治療される病状によって異なる。これらの核の寸法は、内部淡蒼球を除いて数ミリメートルを超えない。
【0010】
この技法によって提示される課題の中には、刺激される核を同定すること（即ち定義する及び／又は突き止めること）、同定核に印を付けること（即ち標的にすること）及び刺激電極を位置決めすることが含まれる。
【0011】
第１の課題は、刺激される核を同定することに関する。例えば、ＶＩＭに関して、人間の視床に対して多くの用語があり、ＶＩＭの同定に関して不一致をもたらす。解剖学的定義に関する問題は、ＶＩＭよりもＳＴＮに対して重要でない。しかし、ＳＴＮのいずれの下部分を刺激すべきかに関して又は更にＳＴＮの上の線維（「不確帯」）が実際に最良の標的でないかどうかに関して不一致が残る。
【0012】
この同定問題を解決した場合でも、同定核を突き止める、従って刺激電極を位置決めする他の課題が依然としてある。例えば、磁気共鳴画像化（ＭＲＩ）を用いて、ジストニーを治療するためのＧＰｉ核を突き止めることができるが、現在の定位診療で使用可能なＭＲＩ（典型的には１．５又は３テスラの磁場を有するＭＲＩ）を用いて、ＳＴＮ核及びＶＩＭ核を突き止めるのは、非常に困難であるか又は更に不可能である。従って、特定の核の突き止めは、ＭＲＩの進歩によって確実に向上しており、例えば、ＭＲＩシーケンス（３Ｄ Ｔ１、Ｔ２、ＳＷＩ（「感受性重み付け画像化」）、ＦＧＡＴＩＲ（「高速灰白質取得反転回復」）、ＦＬＡＩＲ（「流体減衰反転回復」）のシーケンス）の組み合わせは、ＳＴＮ核を視覚化することができるが、ＭＲＩによって放射線的に視覚化されたＳＴＮ核は、（後述の術中電気生理学的記録技法を用いて測定され得る）刺激が有効であるＳＴＮ核に正確に対応せず、５ｍｍを超える誤差が個々のレベルで見られることが実証されている。
【0013】
これらのあまり見えない又は突き止め不可能な領域に対して、定位地図に基づく「間接」突き止め技法と呼ばれるものを実施することが既知の診療である。これらの間接突き止め技法の原理は、ＭＲＩによって容易に同定可能な解剖学的マーカー点を含むデカルト基準フレーム及び様々な脳領域（例えば、様々な脳核）の位置を与える地図を有することであり、解剖学的マーカー点の全部をＭＲＩによって同定することができるとは限らない。原理は、この地図を用いて、マーカー点の全部又は一部に基づいて任意の個人の脳における任意の領域の位置を突き止めることである。マーカー点と地図における点との間の比例の関数のセットを使用し、従って脳のサイズ又は脳の解剖学的特異性にかかわらず、任意の脳を再調整することができる。
【0014】
定位地図は、ＤＢＳへの用途をはるかに越えて、脳を調べる場合により広く使用されている。１つの既知の例は、タライラッハ地図又は座標である。この基準フレームの中心（即ち座標ｘ＝０、ｙ＝０、ｚ＝０を有する点）は、正中矢状面（即ち半球間平面とも呼ばれる脳の内面に平行な平面）における前交連の上及び後縁であり、「ＣＡ」と呼ばれる（「ＡＣ」とも呼ばれる）この点は、ＭＲＩで非常に見やすい。別の基準点、（更に正中矢状面における）後交連の下及び前縁を使用し、「ＣＰ」と呼ばれる（「ＰＣ」とも呼ばれる）この点もＭＲＩで非常に見やすい。次に、座標系の３軸を下記のように定義する。Ｏｙ軸は、ＣＡ及びＣＰを通り、頭蓋骨の前の方に方向付けられる。Ｏｚ軸は、Ｏｙに垂直にＣＡを通る矢状面の軸である。Ｏｚ軸は、頭蓋骨の頂点の方に方向付けられる。最後に、Ｏｘ軸は、ＣＡを通るＯｙ及びＯｚに直交する軸である。Ｏｘ軸は、左から右に方向付けられる。別の例は、定位基準フレームの中心として、中交連点と呼ばれる線分［ＣＡＣＰ］の中点を取ることにより、タライラッハ系を単純化することができるシャルテンブラント地図である。点ＣＡ及びＣＰは、前交連及び後交連の解剖学的中心をそれぞれ表す。
【0015】
刺激される核が、治療される病状に従って同定されると、ＭＲＩによって同定可能なマーカー点を用いて、この核は、使用される地図で突き止められる。深部刺激において、刺激電極は、このようにして間接的に同定される核に位置決めされるべきである。
【0016】
既知の間接突き止め技法に関する問題は、事実と異なる、あたかも２つの点間の距離の変動が全ての観察解剖学的変動に相当するかのように、間接突き止め技法が脳の比例の想定に基づくことである。従って、これは、深部脳刺激で非常に悪影響をもたらすことがある標的化不正確度を必ず引き起こす。
【0017】
これらの標的化不正確度を改善するため、電極が標的に正しく位置決めされることを保証する目的で、電気生理学的記録及び術中臨床試験なしで現在の外科手術を行うことはできない。術中電気生理学的記録の原理は、部位でニューロンを刺激し、術中に刺激の効果を記録するために、解剖学的標的を含むと推定される部位に平行微小電極を埋め込むことである。この技法は、電極の最終的な位置を最適化し、従って標的化を改良することができる。この技法は、標的化精度を向上させるが、非常に長く（７～１０時間のオーダーの動作時間）、感染リスク及び使用される記録微小電極の数に比例する出血のリスクの増加を伴う。更に、電極の最終的な埋め込みを正確に実現することは、必ずしも容易であるとは限らない。最後に、刺激の臨床効果を試験し、従って電極の位置を最適化することができるように、外科手術は、患者の参加を必要とする。従って、患者にとって非常に不快な局所麻酔下で外科手術を行う必要がある。
【0018】
上述の技法は、ＤＢＳの場合、刺激される核を定義し、次に前記核を標的にし、次に刺激電極を位置決めすることを含む手法に基づく。得られた結果の解析により、刺激電極の位置及び／又は使用される基準フレーム（一般的にはシャルテンブラント又はタライラッハ基準フレーム）における解剖学的標的の突き止めを移動させることになる。
【0019】
別の技法は、深部脳刺激が成功している患者を同定し、これらの患者における電極の位置を同定し、ＭＲＩで容易に同定可能な基準点に電極の位置を相関させ、次に基準点を用いて、治療される新しい患者にＭＲＩを適用することを含む間接手法からなる。この手法により、刺激される核を同定し、標的を突き止め（標的にし）、標的における刺激電極を位置決めしている不正確要因の低下を回避することができる。この場合、効果が挙がっている臨床例に基づいて刺激標的を判定することを考えると、「臨床標的」又は「機能標的」を参照する。
【0020】
このような手法は、Ｆｒａｎｃｏｉｓ Ｃａｉｒｅ“Ｉｎｔｒａｏｐｅｒａｔｉｖｅ ｉｍａｇｉｎｇ ｏｆ ｄｅｅｐ ｂｒａｉｎ ｓｔｉｍｕｌａｔｉｏｎ ｅｌｅｃｔｒｏｄｅｓ ａｎｄ ｐｒｏｐｏｓａｌ ｏｆ ａ ｎｅｗ ｍｅｔｈｏｄ ｆｏｒ ｉｎｄｉｒｅｃｔ ｓｔｅｒｅｏｔａｃｔｉｃ ｌｏｃａｔｉｏｎ ｏｆ ｔｈｅ ｓｕｂｔｈａｌａｍｉｃ ｔａｒｇｅｔ”［元の名称“Ｉｍａｇｅｒｉｅ ｐｅｒ－ｏｐｅｒａｔｏｉｒｅ ｄｅｓ ｅｌｅｃｔｒｏｄｅｓ ｄｅ ｓｔｉｍｕｌａｔｉｏｎ ｃｅｒｅｂｒａｌｅ ｐｒｏｆｏｎｄｅ ｅｔ ｐｒｏｐｏｓｉｔｉｏｎ ｄ’ｕｎｅ ｎｏｕｖｅｌｌｅ ｍｏｄａｌｉｔｅ ｄｅ ｒｅｐｅｒａｇｅ ｓｔｅｒｅｏｔａｘｉｑｕｅ ｉｎｄｉｒｅｃｔ ｄｅ ｌａ ｃｉｂｌｅ ｓｕｂｔｈａｌａｍｉｑｕｅ”の翻訳］による論文文献に記載されている。この論文で提案された解決策は、
－ 視床下深部刺激に対する非常に良好な反応を有する患者を選択するステップ、
－ 前記患者の術後ＭＲＩを検索するステップ、
－ 標準ＣＡＣＰ定位空間を各術後ＭＲＩについて構成するステップであって、このために、点ＣＡ及びＣＰ及び半球間平面は、同定される、ステップ、
－ 活性シミュレーション接点の座標を術後ＭＲＩで判定するステップ、
－ 様々な解剖学的マーカー点を術後ＭＲＩで判定するステップ、
－ 線形回帰モデルを用いて、活性接点の各座標と、各マーカー点の対応する座標との間の相関関係を計算するステップ、
－ 相関関係が正しい場合、活性接点のｘ、ｙ及びｚ座標と最良の相関関係を与えるマーカー点に基づいて回帰線の方程式を計算するステップ、
－ この方程式を用いて、理論標的の座標を計算するステップ、
－ この理論標的の座標を活性接点の実際の座標と比較するステップ
を含む。
－ 視床下深部刺激に対する非常に良好な反応を有する患者を選択するステップ、
－ 前記患者の術後ＭＲＩを検索するステップ、
－ 標準ＣＡＣＰ定位空間を各術後ＭＲＩについて構成するステップであって、このために、点ＣＡ及びＣＰ及び半球間平面は、同定される、ステップ、
－ 活性シミュレーション接点の座標を術後ＭＲＩで判定するステップ、
－ 様々な解剖学的マーカー点を術後ＭＲＩで判定するステップ、
－ 線形回帰モデルを用いて、活性接点の各座標と、各マーカー点の対応する座標との間の相関関係を計算するステップ、
－ 相関関係が正しい場合、活性接点のｘ、ｙ及びｚ座標と最良の相関関係を与えるマーカー点に基づいて回帰線の方程式を計算するステップ、
－ この方程式を用いて、理論標的の座標を計算するステップ、
－ この理論標的の座標を活性接点の実際の座標と比較するステップ
を含む。
【0021】
この方法の欠点は、比例システム（ここでは線形回帰方法）を保ちながら、３軸のみに対する脳の変動（又は相似）を想定することである。これは、２つの点（ＣＡ及びＣＰ）間の距離の変動が全ての観察解剖学的変動に相当することがあると想定することよりも優れているが、依然として不十分である。従って、この方法における平均誤差は、２．５±０．６ｍｍであり、標的にされる構造体のサイズに対して非常に大きい。
【先行技術文献】
【非特許文献】
【0022】
【文献】Ｆｒａｎｃｏｉｓ Ｃａｉｒｅ“Ｉｎｔｒａｏｐｅｒａｔｉｖｅ ｉｍａｇｉｎｇ ｏｆ ｄｅｅｐ ｂｒａｉｎ ｓｔｉｍｕｌａｔｉｏｎ ｅｌｅｃｔｒｏｄｅｓ ａｎｄ ｐｒｏｐｏｓａｌ ｏｆ ａ ｎｅｗ ｍｅｔｈｏｄ ｆｏｒ ｉｎｄｉｒｅｃｔ ｓｔｅｒｅｏｔａｃｔｉｃ ｌｏｃａｔｉｏｎ ｏｆ ｔｈｅ ｓｕｂｔｈａｌａｍｉｃ ｔａｒｇｅｔ”
【発明の概要】
【発明が解決しようとする課題】
【0023】
本発明は、先行技術の上述の欠点を克服することを目的とする。
【0024】
より詳細には、本発明は、より高い標的化精度を可能にし、且つ使用が容易である、定位脳標的を判定する方法を提供することを目的とする。本発明は、神経外科治療の効果を改善することができる、神経外科治療前に準備ステップとして使用される定位脳標的を判定する方法を提供することを目的とする。
【課題を解決するための手段】
【0025】
この目的を達成することができる本発明の１つの主題は、少なくとも１つの標的点を含む定位脳標的を判定する方法であって、
－ 少なくとも１つの標的点における治療後に測定される結果が閾値以上である患者を選択するステップであって、術後画像化は、前記患者の各々について実行されている、ステップ、
－ それぞれの選択された患者について、少なくとも１つの標的点の座標の全部又は一部を判定するために前記術後画像化を処理するステップ、
－ 脳マーカー点を選択するステップ、
－ それぞれの選択された患者について、選択マーカー点の座標の全部又は一部を判定するために前記術後画像化を処理するステップ、
－ 選択された患者の全部について判定される標的点の座標及びマーカー点の座標を含む学習データベースを生成するステップ、
－ 学習データベース及び教師あり統計学習方法を使用することにより、マーカー点の座標に従って少なくとも１つの標的点の座標を与える予測関数を判定するステップ、
－ 新しい患者の術前画像化を、前記新しい患者のマーカー点の座標の全部又は一部を判定するために処理するステップ、
－ 前記新しい患者について判定されるマーカー点の座標に従い、前記新しい患者について少なくとも１つの標的点の座標を得るために予測関数を使用するステップ
を含む方法である。
－ 少なくとも１つの標的点における治療後に測定される結果が閾値以上である患者を選択するステップであって、術後画像化は、前記患者の各々について実行されている、ステップ、
－ それぞれの選択された患者について、少なくとも１つの標的点の座標の全部又は一部を判定するために前記術後画像化を処理するステップ、
－ 脳マーカー点を選択するステップ、
－ それぞれの選択された患者について、選択マーカー点の座標の全部又は一部を判定するために前記術後画像化を処理するステップ、
－ 選択された患者の全部について判定される標的点の座標及びマーカー点の座標を含む学習データベースを生成するステップ、
－ 学習データベース及び教師あり統計学習方法を使用することにより、マーカー点の座標に従って少なくとも１つの標的点の座標を与える予測関数を判定するステップ、
－ 新しい患者の術前画像化を、前記新しい患者のマーカー点の座標の全部又は一部を判定するために処理するステップ、
－ 前記新しい患者について判定されるマーカー点の座標に従い、前記新しい患者について少なくとも１つの標的点の座標を得るために予測関数を使用するステップ
を含む方法である。
【0026】
好ましくは、本発明は、少なくとも１つの標的点を含む定位脳標的を判定する方法であって、所与の症状についての前記標的における神経外科治療前に実施され、及び
－ 少なくとも１つの標的点で前記症状について実行される治療後に測定される結果が閾値以上である患者の複数の臨床例を選択するステップであって、術後画像化は、前記患者の各々について実行されている、ステップ、
－ 数学的座標系、好ましくは正規直交デカルト座標系を選択するステップ、
－ それぞれの選択された臨床例について、選択された座標系における少なくとも１つの標的点の座標の全部又は一部を判定するために前記術後画像化を処理するステップ、
－ 複数の脳マーカー点を選択するステップ、
－ それぞれの選択された臨床例について、脳マーカー点の座標の全部又は一部を判定するために前記術後画像化を処理するステップ、
－ 選択された臨床例の全部の標的点の判定された座標及びマーカー点の判定された座標を含む学習データベースを生成するステップ、
－ 学習データベース及び教師あり統計学習方法を使用することにより、マーカー点に従って少なくとも１つの標的点の座標を与える予測関数を判定するステップ、
－ 前記症状について治療される新しい患者の術前画像化を、前記新しい患者のマーカー点の座標の全部又は一部を判定するために処理するステップ、
－ 前記新しい患者について判定されるマーカー点の座標に従い、前記新しい患者について少なくとも１つの標的点の座標を得るために予測関数を使用するステップ
を含む方法に関する。
－ 少なくとも１つの標的点で前記症状について実行される治療後に測定される結果が閾値以上である患者の複数の臨床例を選択するステップであって、術後画像化は、前記患者の各々について実行されている、ステップ、
－ 数学的座標系、好ましくは正規直交デカルト座標系を選択するステップ、
－ それぞれの選択された臨床例について、選択された座標系における少なくとも１つの標的点の座標の全部又は一部を判定するために前記術後画像化を処理するステップ、
－ 複数の脳マーカー点を選択するステップ、
－ それぞれの選択された臨床例について、脳マーカー点の座標の全部又は一部を判定するために前記術後画像化を処理するステップ、
－ 選択された臨床例の全部の標的点の判定された座標及びマーカー点の判定された座標を含む学習データベースを生成するステップ、
－ 学習データベース及び教師あり統計学習方法を使用することにより、マーカー点に従って少なくとも１つの標的点の座標を与える予測関数を判定するステップ、
－ 前記症状について治療される新しい患者の術前画像化を、前記新しい患者のマーカー点の座標の全部又は一部を判定するために処理するステップ、
－ 前記新しい患者について判定されるマーカー点の座標に従い、前記新しい患者について少なくとも１つの標的点の座標を得るために予測関数を使用するステップ
を含む方法に関する。
【0027】
本発明によれば、マーカー点は、所与の数学的座標系で突き止め可能な脳における解剖学的構造体の特徴点であり、特に従来の画像化を用いて見える解剖学的点である。
【0028】
本発明によれば、標的点は、神経外科治療が臨床例に適用されているか、又は新しい患者に適用される必要がある点である。
【0029】
本発明による方法は、治療の効果が判定されている患者の臨床例からのデータに基づく。従って、学習を行う標的点は、臨床的に検証される標的である。これにより、本発明による方法に従い得る神経外科治療を改善することができる。
【0030】
これを行うために、本発明による方法は、これらの患者に利用できる術後画像化に基づいており、治療が適用された点（標的点）及び特に従来の画像化を用いて見える脳における解剖学的構造体の特徴点の中から選択されるマーカー点から推定するために、この画像化を処理する。これらの標的点及びこれらのマーカー点の座標は、予測関数を判定することができる学習データベースに統合される。予測関数は、マーカー点に応じて少なくとも１つの標的点を与えることができる。教師あり統計学習方法を使用して予測関数を学習データベースから構成する。
【0031】
本発明による方法は、従来使用されるＭＲＩ（例えば、１．５テスラＭＲＩ）から抽出される画像化に基づき得る。よりはるかに高価であり（特にハードウェアの点で）、１．５テスラＭＲＩよりも高い放射線技能を必要とするより強力なＭＲＩ（例えば、３テスラ又は７テスラＭＲＩ）を使用する必要はない。マーカー点の座標を判定するためにこれらの画像を処理する。
【0032】
本発明による方法は、１軸、２軸又は３軸ではなく、多次元方法で脳変動を制御しようとする。教師あり統計学習技法を選択することにより、方法は、この脳変動を克服することができ、従ってより高い標的化精度を可能にする。
【0033】
更に、刺激される脳構造体（ＳＴＮ、ＶＩＭなど）をまず同定し、次に突き止め、その後、このようにして同定され且つ突き止められた脳構造体に刺激電極を埋め込もうとする他の方法と異なり、標的点は、臨床的に検証される標的として判定される。従って、これにより、不正確要因の低下を回避し、従って標的化精度を向上させることができる。
【0034】
この意味において、本発明による方法は、マーカー点に応じて脳構造体（ＳＴＮ、ＶＩＭなど）を視覚化するために解剖学的地図を構成しようとしない。むしろ、本発明による方法は、治療を適用する場所を直接知るために標的点を求め、患者に有効である標的点を使用することによってそれを行う。実際に、本発明者らは、解剖学的構造体（ＳＴＮ、ＶＩＭなど）の中心を突き止めても、上述のように、この中心が、刺激される正確な標的であることを保証しないことを特定している。
【0035】
更に、本発明による方法は、１つの点だけでなく、標的（複数の点）を判定し、従って三次元で標的を構成することができる。
【0036】
標的が新しい患者について判定されると、ＭＲＩにおける標的の位置に印（例えば、ばつ印）を付けることができる。ばつ印を付けるために、最大値（白色に対応する）を５つの画素に割り当てることにより、数学的座標系の３方向の各々で５つの画素の値を変更することができる。これにより、白色のばつ印を形成することができ、この白色のばつ印の中心により、専門家は、標的を視覚化することができる。
【0037】
最後に、本発明によって定位標的を判定する方法は、決定支援として、特に判定される少なくとも１つの正確な脳標的を必要とする任意の神経外科治療に備えて使用され得る。
【0038】
特に深部脳刺激の場合、本発明の方法は、刺激される脳構造体に電極を正確に位置決めすることを保証するために一般的に使用される電気生理学を不要にすることができる。従って、本発明の方法は、患者に対して侵襲的でない。
【0039】
予測関数（「メタモデル」とも呼ばれる）のパラメータを選択及び最適化する適切な技法を使用することにより、方法は、標的化精度を更に向上させることができる。
【0040】
従って、１つの有利な実施形態によれば、方法は、
－ 交差検証方法を使用して予測関数を統合するステップであって、マーカー点に従って少なくとも１つの標的点の座標を与える統合予測関数をもたらすステップ
を更に含み得、予測関数を使用するステップは、統合予測関数を使用するステップを含む。
－ 交差検証方法を使用して予測関数を統合するステップであって、マーカー点に従って少なくとも１つの標的点の座標を与える統合予測関数をもたらすステップ
を更に含み得、予測関数を使用するステップは、統合予測関数を使用するステップを含む。
【0041】
一実施形態によれば、マーカー点及び少なくとも１つの標的点の座標を判定するために処理される画像化は、少なくとも１つのＭＲＩ画像、好ましくは複数のＭＲＩ画像である。
【0042】
従って、方法は、選択された患者における術後ＭＲＩを実行するステップを含み得、マーカー点の判定は、ＭＲＩによって得られる画像を処理するステップを含み得、且つ／又は少なくとも１つの標的点の判定は、ＭＲＩによって得られる画像を処理するステップを含み得る。
【0043】
代わりに又は更に、方法は、選択された患者における術後コンピュータ断層撮影（ＣＴ）スキャンを実行するステップを含み得、少なくとも１つの標的点の判定は、ＣＴによって得られる画像を処理するステップを含み得、且つ／又はマーカー点の判定は、ＣＴによって得られる画像を処理するステップを含み得る。
【0044】
一実施形態によれば、教師あり統計学習方法は、再生カーネルヒルベルト空間におけるカーネルリッジ回帰方法の使用を含む。
【0045】
別の実施形態によれば、教師あり統計学習方法は、サポートベクトルマシン型の方法の使用を含む。
【0046】
別の実施形態によれば、教師あり統計学習方法は、ニューラルネットワーク型の方法の使用を含む。
【0047】
一実施形態によれば、交差検証方法は、「リーブワンアウト交差検証」方法の使用を含む。
【0048】
別の実施形態によれば、交差検証方法は、「リーブｋアウト交差検証」方法の使用を含む。
【0049】
１つの特定の実施形態によれば、選択数学的座標系は、正規直交デカルト座標系であり、直線は、Ｏｙ軸を形成する前交連の上及び後縁並びに後交連の下及び前縁を通り、後交連の下及び前縁は、座標系の中心を形成し、及びＯｚ軸は、半球間平面においてＯｙ軸に垂直な直線である。
【0050】
１つの特定の実施形態によれば、マーカー点は、好ましくは、事前に選択される数学的座標系に対して定義される下記の１８個の点：
－ 第３の軸平面上の乳頭視床束である第１のマーカー点、
－ 第１、第２及び第３の軸平面の各々の上の被殻の前点である第２、第３及び第４のマーカー点、
－ 第１及び第２の軸平面上の被殻の内側点である第５及び第６のマーカー点、
－ 第１及び第２の軸平面上の被殻の後点である第７及び第８のマーカー点、
－ 第２の軸平面上の手綱交連である第９のマーカー点、
－ 第２の軸平面上の視床の前縁である第１０のマーカー点、
－ 第２の軸平面上の視床の後縁である第１１のマーカー点、
－ 前交連である第１２のマーカー点、
－ 中交連点における第３脳室の内側縁である第１３のマーカー点、
－ 第１３のマーカー点を通る矢状面上の視床の高さである第１４のマーカー点、
－ 第１３及び第１４のマーカー点によって定義される線分の中点である第１５のマーカー点、
－ 前交連の上及び後縁並びに後交連の下及び前縁を通る直線に平行であり、且つ第１５のマーカー点を通る直線上の視床の前縁である第１６のマーカー点、
－ 第５のマーカー点を通る前頭面上の被殻の上縁である第１７のマーカー点、及び
－ 第５のマーカー点を通る前頭面上の被殻の外側縁である第１８のマーカー点
の中から選択される。
－ 第３の軸平面上の乳頭視床束である第１のマーカー点、
－ 第１、第２及び第３の軸平面の各々の上の被殻の前点である第２、第３及び第４のマーカー点、
－ 第１及び第２の軸平面上の被殻の内側点である第５及び第６のマーカー点、
－ 第１及び第２の軸平面上の被殻の後点である第７及び第８のマーカー点、
－ 第２の軸平面上の手綱交連である第９のマーカー点、
－ 第２の軸平面上の視床の前縁である第１０のマーカー点、
－ 第２の軸平面上の視床の後縁である第１１のマーカー点、
－ 前交連である第１２のマーカー点、
－ 中交連点における第３脳室の内側縁である第１３のマーカー点、
－ 第１３のマーカー点を通る矢状面上の視床の高さである第１４のマーカー点、
－ 第１３及び第１４のマーカー点によって定義される線分の中点である第１５のマーカー点、
－ 前交連の上及び後縁並びに後交連の下及び前縁を通る直線に平行であり、且つ第１５のマーカー点を通る直線上の視床の前縁である第１６のマーカー点、
－ 第５のマーカー点を通る前頭面上の被殻の上縁である第１７のマーカー点、及び
－ 第５のマーカー点を通る前頭面上の被殻の外側縁である第１８のマーカー点
の中から選択される。
【0051】
一実施形態によれば、これらの１８個のマーカー点の全部が使用される。
【0052】
他のマーカー点が使用され得る。
【0053】
他の数学的座標系、特に既知の定位地図でよく使用される座標系を使用し得る。
【0054】
一実施形態によれば、方法は、学習データベースに機能データを追加する追加ステップを含み、前記機能データは、臨床例の標的点及びマーカー点における信頼度の少なくとも１つの指標を追加することができる。
【0055】
一実施形態によれば、方法は、深部脳刺激前に実施される。別の実施形態によれば、方法は、ガンマナイフ前に実施される。別の実施形態によれば、方法は、集束超音波治療前に実施される。
【0056】
マーカー点又は標的点は、計画的治療に応じて、従って治療される症状に応じて適合される。
【0057】
本発明の第２の主題は、方法のステップの全部又は一部を実施するように構成されたプロセッサを含むデータ処理システムである。データ処理システムは、コンピュータ、タブレット、スマートフォンなどであり得る。
【0058】
本発明の第３の主題は、コンピュータプログラムであって、プログラムがプロセッサによって実行されると、このプロセッサに、方法のステップの全部又は一部を実施させる命令を含むコンピュータプログラムである。
【0059】
本発明の他の特徴及び利点は、添付図面を参照して、非限定的な例として与えられる下記の説明によって明瞭になるであろう。
【図面の簡単な説明】
【0060】
【図1】２つの異なる点及び２つの異なる方向における２つの刺激電極の埋め込みを例示する合成画像である。
【図2】標的に埋め込まれ、複数の活性接点を含む電極を例示する合成画像である。
【図3A】１８個のマーカー点を例示する様々な断面に沿って見られるＭＲＩ画像である。
【図3B】１８個のマーカー点を例示する様々な断面に沿って見られるＭＲＩ画像である。
【図3C】１８個のマーカー点を例示する様々な断面に沿って見られるＭＲＩ画像である。
【図3D】１８個のマーカー点を例示する様々な断面に沿って見られるＭＲＩ画像である。
【図3E】１８個のマーカー点を例示する様々な断面に沿って見られるＭＲＩ画像である。
【図3F】１８個のマーカー点を例示する様々な断面に沿って見られるＭＲＩ画像である。
【図4A】数学的座標系ＣＡ－ＣＰの軸を例示する３つのＭＲＩ画像である。
【図4B】数学的座標系ＣＡ－ＣＰの軸を例示する３つのＭＲＩ画像である。
【図4C】数学的座標系ＣＡ－ＣＰの軸を例示する３つのＭＲＩ画像である。
【図4D】三次元における数学的座標系ＣＡ－ＣＰを示す。
【図5】本発明による例示的な方法及び変型例を示すフローチャートである。
【発明を実施するための形態】
【0061】
本発明は、神経外科治療のための標的点を判定することを含み、これらの標的点は、一般的に、従来の画像化によって見えない深部脳内構造体に位置決めされる。
【0062】
このために、本発明は、従来の画像化によって見えるマーカー点及び予測関数を使用することを含む。治療の効果が観察されている患者の臨床例によって生じるマーカー点及び標的点を含む学習データベースを用いて、前記予測関数を構成する。更に後述のように、少なくとも１つの測定結果によって効果を判定する。教師あり統計学習技法を使用して前記予測関数（又はメタモデル）を構成し、従ってマーカー点に応じて少なくとも１つの標的点の座標を与えることができる。
【0063】
下記の実施形態の説明は、本発明による方法を深部脳刺激による治療に対する準備として実施する場合において一般的に例示され、本発明、特に前記実施形態は、判定される少なくとも１つの正確な脳標的を必要とする任意の神経外科治療に対する準備としてより一般的に実施され得ることが理解される。
【0064】
深部脳刺激への特定の用途において、標的点を刺激点と呼ぶこともできる。
【0065】
臨床例（患者）の選択
少なくとも１つの刺激点における電極刺激が同じ病状を治療するのに効果的である患者を選択する。少なくとも１つの測定結果によって効率を判定する。
少なくとも１つの刺激点における電極刺激が同じ病状を治療するのに効果的である患者を選択する。少なくとも１つの測定結果によって効率を判定する。
【0066】
病状は、パーキンソン病、本態性振戦病、ジストニー、強迫性障害、ツレット症候群、不応性てんかん又は更に治療の効かない重度のうつ病であり得る。
【0067】
刺激に対する標的構造体は、パーキンソン病を治療するための視床下核（ＳＴＮ）、本態性振戦を治療するための視床の中間腹側核（ＶＩＭ）、ジストニー及びある種のパーキンソン病を治療するための内部淡蒼球（ＧＰｉ）であり得る。
【0068】
神経外科治療後に測定されるこの結果は、術後臨床評価、解剖学的又は放射線局在による評価又は電気生理学的評価であり得る。
【0069】
治療される病状に特有のスケール及び／又は生活の質のスケール（例えば、治療後３カ月）を用いて術後臨床評価を行う。これにより、第１の基準、従って臨床例（患者）の選択のための第１の閾値を判定することができる。
【0070】
解剖学的又は放射線術前標的化及び前記術前標的化と術後治療点（例えば、深部脳刺激（ＤＢＳ）の場合の電極の位置）との間の距離により、解剖学的又は放射線局在による評価を定義する。これにより、第２の基準、従って臨床例（患者）の選択のための第２の閾値を判定することができる。
【0071】
標的の術中同定及び前記術中同定と術後治療点（例えば、ＤＢＳの場合の電極の位置）との間の距離により、電気生理学的評価を定義する。これにより、第３の基準、従って臨床例（患者）の選択のための第３の閾値を判定することができる。
【0072】
結合閾値を定義するために様々な評価基準が結合され得る。
【0073】
選択された患者は、少なくとも１つの病状についてマーカー点／標的点の対を含む統計学習データベースを構築することができる。
【0074】
更に、方法は、患者の学習データベースに機能データを追加するステップを含み得る。これらのデータは、術前に病状の兆候を評価するパラメータ及び術後に病状の兆候の改善の程度を評価するパラメータを含み得る。一例として、パーキンソン病に対するＵＰＤＲＳスケール（パーキンソン病統一スケール）又は本態性振戦に対するＦａｈｎ－Ｔｏｌｏｓａ－Ｍａｒｉｎスケールを使用し得る。
【0075】
更に、深部脳刺激の場合、これらのデータは、刺激の電気的パラメータ（振幅、周波数及び持続時間）を含み得る。
【0076】
従って、これらの機能データの全部は、標的点及びマーカー点に加えて、学習データベース及びメタモデルに入力され得る。これらの機能データの全部又は一部は、患者データにおける信頼度の少なくとも１つの指標として費用関数に使用され得る。この場合、これらのデータは、最小化される費用関数を重み付けする重み関数としてメタモデルに数学的に変換される。
【0077】
マーカー点は、使用される数学的座標系で突き止め可能な脳における解剖学的構造体の特徴点であり、従来の画像化（特にＭＲＩ）を用いて見える解剖学的点である。
【0078】
従って、既存の術後ＭＲＩによって生じる画像化を使用することにより、マーカー点を得ることができる。
【0079】
記載の特定の実施形態において、標的点（又は刺激点）は、電極の埋め込みの点である。刺激点は、使用される数学的座標系における座標によって特徴付けられる。刺激点は、図１及び図２に例示のように、この点（ＰＣ１、ＰＣ２）に埋め込まれた電極（Ｅ１、Ｅ２）の方向によって特徴付けられ得る。
【0080】
例えば、ＭＲＩ又はコンピュータ断層撮影法（「ＣＴ」）による術後画像化を実行及び使用することにより、且つ末端アーチファクト及び様々な接触によって生成されたアーチファクトに基づいて基準フレームで電極を再構成することにより各刺激点を判定する。
【0081】
本明細書において、標的点又は標的点のセットは、深部脳刺激の特定の用途で「標的」又は「刺激標的」を構成する。標的は、各標的点の座標によって特徴付けられる。深部脳刺激の特定の用途において、標的は、刺激点における電極の方向（電極の方向ベクトルの座標）によって更に特徴付けられ得る。
【0082】
従って、選択された患者について判定されたこれらのマーカー点及び標的点は、他のデータで拡充可能な統計学習データベースを構成する。
【0083】
本発明による方法は、マーカー点と１つ又は複数の標的点との間の予測関数を確定するように、前記学習データベースを用いた、教師あり統計学習に基づく回帰方法の使用を含む。
【0084】
【0085】
学習データベースを構成及び／又は初期化するステップ
使用される数学的座標系
「ＣＡ」は、前交連の上及び後縁を示し、「ＣＰ」は、後交連の下及び前縁を示す。
使用される数学的座標系
「ＣＡ」は、前交連の上及び後縁を示し、「ＣＰ」は、後交連の下及び前縁を示す。
【0086】
【0087】
例示の座標系は、点ＣＡ及びＣＰ（記号「ＣＡＣＰ」）を通る直線であるＯｙ軸並びに点ＣＰを原点として有する正規直交デカルト座標系である。Ｏｚ軸は、半球間平面において線（ＣＡＣＰ）に垂直な直線であり、Ｏｘ軸は、推定される。このような座標系は、本明細書ではＣＡ－ＣＰによって記号で表す。
【0088】
使用可能なマーカー点
１つの例示的な実施形態において、使用可能なマーカー点は、特に、図３Ａ～図３Ｆに示す点に対応する１８個（ｍ＝１８）に番号を付ける。１８個の使用可能なマーカー点は、下記の通りである。
－ 第１のマーカー点（ＦＭＴ）：第３の軸平面上の乳頭視床束、
－ 第２、第３及び第４のマーカー点（ＰＡ１、ＰＡ２、ＰＡ３）：第１、第２及び第３の軸平面の各々の上の被殻の前点、
－ 第５及び第６のマーカー点（ＰＭ１、ＰＭ２）：第１及び第２の軸平面上の被殻の内側点、
－ 第７及び第８のマーカー点（ＰＰ１、ＰＰ２）：第１及び第２の軸平面上の被殻の後点、
－ 第９のマーカー点（ＣＨ）：第２の軸平面上の手綱交連、
－ 第１０のマーカー点（ＢＡＴ）：第２の軸平面上の視床の前縁、
－ 第１１のマーカー点（ＢＰＴ）：第２の軸平面上の視床の後縁、
－ 第１２のマーカー点（ＣＡ）：前交連、
－ 第１３のマーカー点（Ａ）：「ＭＣＰ」点と呼ばれる中交連点における第３脳室の内側縁、
－ 第１４のマーカー点（Ｂ）：第１３のマーカー点（Ａ）を通る矢状面上の視床の高さ、
－ 第１５のマーカー点（Ｃ）：第１３及び第１４のマーカー点によって定義される線分［ＡＢ］の中点、
－ 第１６のマーカー点（Ｄ）：（ＣＡＣＰ）に平行であり、第１５のマーカー点（Ｃ）を通る直線上の視床の前縁、
－ 第１７のマーカー点（Ｐｃｕｌｍ）：第５のマーカー点（ＰＭ１）を通る前頭面上の被殻の上縁、及び
－ 第１８のマーカー点（Ｐｌａｔ）：第５のマーカー点（ＰＭ１）を通る前頭面上の被殻の外側縁。
１つの例示的な実施形態において、使用可能なマーカー点は、特に、図３Ａ～図３Ｆに示す点に対応する１８個（ｍ＝１８）に番号を付ける。１８個の使用可能なマーカー点は、下記の通りである。
－ 第１のマーカー点（ＦＭＴ）：第３の軸平面上の乳頭視床束、
－ 第２、第３及び第４のマーカー点（ＰＡ１、ＰＡ２、ＰＡ３）：第１、第２及び第３の軸平面の各々の上の被殻の前点、
－ 第５及び第６のマーカー点（ＰＭ１、ＰＭ２）：第１及び第２の軸平面上の被殻の内側点、
－ 第７及び第８のマーカー点（ＰＰ１、ＰＰ２）：第１及び第２の軸平面上の被殻の後点、
－ 第９のマーカー点（ＣＨ）：第２の軸平面上の手綱交連、
－ 第１０のマーカー点（ＢＡＴ）：第２の軸平面上の視床の前縁、
－ 第１１のマーカー点（ＢＰＴ）：第２の軸平面上の視床の後縁、
－ 第１２のマーカー点（ＣＡ）：前交連、
－ 第１３のマーカー点（Ａ）：「ＭＣＰ」点と呼ばれる中交連点における第３脳室の内側縁、
－ 第１４のマーカー点（Ｂ）：第１３のマーカー点（Ａ）を通る矢状面上の視床の高さ、
－ 第１５のマーカー点（Ｃ）：第１３及び第１４のマーカー点によって定義される線分［ＡＢ］の中点、
－ 第１６のマーカー点（Ｄ）：（ＣＡＣＰ）に平行であり、第１５のマーカー点（Ｃ）を通る直線上の視床の前縁、
－ 第１７のマーカー点（Ｐｃｕｌｍ）：第５のマーカー点（ＰＭ１）を通る前頭面上の被殻の上縁、及び
－ 第１８のマーカー点（Ｐｌａｔ）：第５のマーカー点（ＰＭ１）を通る前頭面上の被殻の外側縁。
【0089】
軸方向基準面（又は横方向平面）は、ＣＡ－ＣＰ座標系の中心及びＯｘ軸及びＯｙ軸を通る平面である。第１、第２及び第３の軸平面を軸方向基準面に対して定義する。
【0090】
第１の軸平面は、ｚ（Ｃ）－５ｍｍの座標（第１５のマーカー点のｚ座標から５ｍｍ下に位置する点並びにゼロに等しいｘ及びｙ座標）を通る軸方向基準面に平行な平面に対応する。
【0091】
第２の軸平面は、ｚ（Ｃ）座標（第１５のマーカー点のｚ座標並びにゼロであるｘ及びｙ座標）を通る軸方向基準面に平行な平面に対応する。
【0092】
第３の軸平面は、ｚ（Ｃ）－１０ｍｍの座標（第１５のマーカー点のｚ座標から１０ｍｍ下並びにゼロであるｘ及びｙ座標）を通る軸方向基準面に平行な平面に対応する。
【0093】
矢状面によって理解されるべきことは、Ｏｙ軸及びＯｚ軸によって形成される平面と一致するか又はこの平面に平行な任意の平面である。更に、実施形態で選択される矢状面は、第１３のマーカー点（Ａ）のｘ座標を通る。
【0094】
前頭面によって理解されるべきことは、Ｏｘ軸及びＯｚ軸によって形成される平面と一致するか又はこの平面に平行な任意の平面である。更に、実施形態で選択される前頭面は、第５のマーカー点（ＰＭ１）のｙ座標を通る。
【0095】
図３Ａ～図３Ｆは、様々な断面に沿って見られるＭＲＩ画像を示す。図３Ａは、軸方向断面（軸方向基準面）である。図３Ｂは、上述の矢状面に沿った矢状断面である。図３Ｃは、軸方向断面（第１の軸平面）である。図３Ｄは、軸方向断面（第２の軸平面）である。図３Ｅは、軸方向断面（第３の軸平面）である。図３Ｆは、上述の前頭面に沿った冠状断面である。
【0096】
計画的神経外科的介入により、他のマーカー点を使用し得る。従って、標的点は、治療される病状に左右されるため、一般的に、前記標的点と相関可能であり、従来の画像化（例えば、ＭＲＩ）で容易に同定可能であるマーカー点を選択し得る。
【0097】
特定の実施形態によれば、マーカー点ＰＲは、行列Ｘ（ここで、
であり、ｍは、マーカー点の数を表し、ｍ×３は、これらの全マーカー点の全座標を表す）で表され得る。
【数1】
【0098】
行列Ｘの各行は、所与のマーカー点を表す。各列は、３つの座標の１つを表す。
【0099】
マーカー点の座標を連結することにより、行列Ｘは、ベクトル
に再形成され得る。残りの説明の全体にわたり、複数（３×ｍ）のスカラー要素を含むこのようなベクトルＸを考える。
【数2】
【0100】
１８個のマーカー点の各々は、使用される数学的座標系における３つの座標ｘＰＲ、ｙＰＲ、ｚＰＲによって特徴付けられる。これは、合計５４個の座標を形成する。
【0101】
一般的に、半球（右半球（本明細書では「ＨＤ」で示す）及び左半球（本明細書では「ＨＧ」で示す））によって判定を行う。次に、ＨＤで１８個の点又は５４個の座標があり、ＨＧで１８個の点又は５４個の座標がある。この場合、ｍ＝３６及び合計１０８個の座標が得られる。
【0102】
従って、左側の標的を判定することができる、左半球に対する１つの予測関数ＦＧを構成することができ、且つ右側の標的を判定することができる、右半球に対する別の予測関数ＦＤを構成することができるか、又は両方の標的（右側の標的及び左側の標的）を判定することができる単一の関数Ｆを構成することができる。
【0103】
標的点（又は刺激点）
標的（１つ又は複数の標的点ＰＣをグループ化する）は、ベクトルＹ（ここで、
である）で表される。ベクトルＹは、複数（ｐ）のスカラー要素を含む。
標的（１つ又は複数の標的点ＰＣをグループ化する）は、ベクトルＹ（ここで、
【数3】
【0104】
標的における標的点の数は、可変である。値ｐは、点毎の座標の数が掛けられた標的点の数に左右される。座標は、標的点のデカルト座標ｘＰＣ、ｙＰＣ、ｚＰＣであり得る。
【0105】
深部脳刺激の場合、座標は、この点における刺激電極の方向を特徴付ける座標を更に含み得る。刺激点及びこの点における電極の方向を知ることにより、同じ電極上の他の刺激点を見付けることができる。刺激点は、「刺激ゾーン」とも呼ばれる。各電極Ｅ１、Ｅ２は、例えば、図２に示すように、例えば４つのゾーンを有し得る。点０は、標的点ＰＣ１、ＰＣ２を表し得る。従って、半球毎に１つ又は２つのゾーンを介して患者を刺激し得る。図２の例において、ゾーン０及び１は、標的内に全部あり、ゾーン２は、標的内に部分的にあり、ゾーン３は、標的内にない。例えば、各ゾーンの長さは、１．５ｍｍであり、２つのゾーンは、０．５ｍｍ間隔で配置されている。
【0106】
標的点に対する座標のみを求める場合、使用される数学的座標系で標的点の３つの座標ｘＰＣ、ｙＰＣ、ｚＰＣに対応するｐ＝３である。
【0107】
深部脳刺激の場合、刺激点に対する電極の方向を更に求める場合、標的ｘＰＣ、ｙＰＣ、ｚＰＣ及び使用される数学的座標系の３つの軸に沿った電極の方向ベクトルの座標ｖｘＰＣ、ｖｙＰＣ、ｖｚＰＣに対応するｐ＝６である。
【0108】
計画的神経外科的介入により、他の座標を追加し得る。例えば、ガンマナイフ又は集束超音波治療に対して標的点の座標を適合させ得る。
【0109】
一方では右半球に対して、他方では左半球に対して標的点の座標を確定することを望む場合、ｐ＝６（デカルト座標のみ）又はｐ＝１２（ベクトル座標を含む）である。
【0110】
従って、標的において１つ又は複数の標的点（又はゾーン）があり得るため、値ｐは、各点に対して使用される座標の数が掛けられた点の数に等しい。
【0111】
学習データベース
患者ｉ（ｉ＝１．．．ｎ）に対して、マーカー点の（３×ｍ値ｘｉの）ベクトルＸｉ及び標的点の（ｐ値ｙｉの）ベクトルＹｉがある。ｎ人の患者から抽出されたベクトル（Ｘｉ、Ｙｉ）ｉ＝１．．．ｎは、全学習データベースを定義する。
患者ｉ（ｉ＝１．．．ｎ）に対して、マーカー点の（３×ｍ値ｘｉの）ベクトルＸｉ及び標的点の（ｐ値ｙｉの）ベクトルＹｉがある。ｎ人の患者から抽出されたベクトル（Ｘｉ、Ｙｉ）ｉ＝１．．．ｎは、全学習データベースを定義する。
【0112】
学習ステップ（関数Ｙ＝Ｆ（Ｘ）を構成する）
教師あり統計学習方法を使用して、本発明による方法は、Ｙ≒Ｆ（Ｘ）であるように予測関数Ｆを構成する。
教師あり統計学習方法を使用して、本発明による方法は、Ｙ≒Ｆ（Ｘ）であるように予測関数Ｆを構成する。
【0113】
予測関数は、本明細書では単に「関数」と呼ばれ得る。
【0114】
関数Ｆは、Ｙを生成する引数Ｘ又はＸの一部を取る。
【0115】
標的を判定する必要がある新しい患者について、マーカー点のベクトルＸ’をＭＲＩ（例えば、手術前に実行される術前ＭＲＩ）から判定する必要がある。
【0116】
次に、下記のような関数Ｆを用いて標的の位置Ｙ’を判定する。
【数4】
【0117】
【数5】
【0118】
１つの特定の実施形態によれば、教師あり統計学習方法は、再生カーネルヒルベルト空間におけるカーネルリッジ回帰方法の使用を含む。
【0119】
下記の記載は、この特定の実施形態を説明する。
【0120】
カーネルリッジ回帰は、比較的単純な方法である。カーネルリッジ回帰は、一般的に、全学習データベースを使用し、「小さい」学習データベース（患者（ｎ）の数が２０未満である）又は「中程度の」データベース（患者（ｎ）の数が２０～１００である）に適している。
【0121】
正規化項
（ここで、λは、正規化係数である）を有する最小二乗法を用いて、関数Ｆをヒルベルト空間
で求める。正規化は、ノルム
であり、従って関数Ｆが退化しないことを保証する。
【数6】
【数7】
【数8】
【0122】
従って、選択ヒルベルト空間
は、下記のようなガウスカーネルに基づく再生カーネルヒルベルト空間である。
【数9】
【数10】
【0123】
カーネルは、下記の式（１）によって与えられる特性により、空間
で「再生する」と言われる。
ここで、
は、
におけるスカラー積である。
【数11】
【数12】
【数13】
【数14】
【0124】
下記の正規化最小二乗問題をヒルベルト空間
で考える。学習データベース
（ここで、ｎは、データベースにおける要素の数である）を考えると、目的は、下記の式（２）を解くことである。
ここで、λ＞０は、所与の正規化係数である。
【数15】
【数16】
【数17】
【0125】
リースの表現定理は、式（２）に対する解を下記のように記述し得ることを示す。
ここで、
は、新しい未知数である。
【数18】
【数19】
【0126】
Ｋｉｊ＝Ｋ（ｘｉ，ｘｊ）でグラム行列Ｋ＝（Ｋｉｊ）ｉ，ｊ＝１．．．ｎを定義することにより、解かれる式（２）は、下記の式（３）を解くことと同等である。
解は、式（４）によって簡単に与えられる。
α＝（Ｋ＋λｎＩ）－１ｙ （４）
Ｉは、
における恒等行列である。
【数20】
α＝（Ｋ＋λｎＩ）－１ｙ （４）
Ｉは、
【数21】
【0127】
標的（１つ又は複数の標的点をグループ化する）を、スカラー特性のみによって定義するのでなく、各標的点の座標、及び／又は各点で位置決めされる刺激電極の方向（３つの軸におけるベクトルの座標によって特徴付けられる）、及び／又は計画的神経外科的介入に応じた他の座標を含むサイズｐのＹベクトルによって定義するため、これは、下記の式（５）を解く問題である。
Ａ＝（Ｋ＋λｎＩ）－１Ｙ （５）
Ｉは、
における恒等行列である。
である。
Ａ＝（Ｋ＋λｎＩ）－１Ｙ （５）
Ｉは、
【数22】
【数23】
【0128】
Ａの各列は、標的（１とｐとの間）の所与の座標に対する係数αに対応する。Ｙの各列は、標的（１とｐとの間）の所与の座標に対する値に対応する。
【0129】
上述のように、これを知ると、Ｋは、下記のようにグラム行列によって定義されるガウスカーネルである。
【数24】
【0130】
再構成ステップ（新しい患者）
学習ステップが終了する、即ち関数Ｆが構成されると、標的を再構成する、即ち下記の式（７）を用いて、例えば患者のＭＲＩから判定されるマーカー点
のセットに基づいて、１つ又は複数の標的点の座標のベクトル
を判定することができる。
ここで、
である。
学習ステップが終了する、即ち関数Ｆが構成されると、標的を再構成する、即ち下記の式（７）を用いて、例えば患者のＭＲＩから判定されるマーカー点
【数25】
【数26】
【数27】
【数28】
【0131】
データベースに機能データを追加する場合、下記の表記法を導入する。
【0132】
少なくとも１つのベクトルＸｉは、１８個のマーカー点などの術前解剖学的データの全部及び病状の兆候に対する評価スコアなどの機能データを考慮する。従って、ベクトルＸｉのサイズは、３×ｍ＋Ｍ（ここで、Ｍは、機能データ、例えば患者の術前兆候（又は「術前スコア」）の評価によって生じるパラメータの数である）である。
【0133】
データベースにおける各患者ｉに対して、ベクトルｐｏｓｔＯｐＶｅｃｔｏｒｉを定義する。このベクトルは、例えば、術後兆候（又は「術後スコア」）の評価によって生じるパラメータ及び／又は刺激の電気的パラメータ（周波数、振幅及び持続時間）を含み得る。
【0134】
機能データを考慮する場合に最小化される関数は、下記の通りである。
ここで、Ｗは、重み関数であり、患者データにおける信頼度測定する。この重み関数は、複雑又は単純であるように選択され得る。重み関数の１つの例は、刺激振幅の逆数又は１秒毎に注入される電流の量の逆数であり、刺激パラメータの積の逆数に対応する。
【数29】
【0135】
刺激パラメータの積の逆数の場合、下記の通りである。
【数30】
【0136】
重み関数の別の例は、スコアが低いほど、電極がより良好に配置されているという事実を反映する患者の術後スコアの和の逆数である。
【0137】
重み関数Ｗは、術後スコア及び電気的パラメータを組み合わせた別の表現を取り得る。
【0138】
選択ステップ
式（５）を解くことによって得られる関数は、ＭＲＩから抽出される１８個のマーカー点の５４個の座標（又は脳の両方の半球を考慮する場合、３６個の点の１０８個の座標）を学習及び再構成ステップで使用するという意味において、マーカー点又はこれらの点の座標の任意の選択を含まない。しかし、幾つかの座標は、標的と相関性がない。これは、標的を判定する際の精度を損なう。
式（５）を解くことによって得られる関数は、ＭＲＩから抽出される１８個のマーカー点の５４個の座標（又は脳の両方の半球を考慮する場合、３６個の点の１０８個の座標）を学習及び再構成ステップで使用するという意味において、マーカー点又はこれらの点の座標の任意の選択を含まない。しかし、幾つかの座標は、標的と相関性がない。これは、標的を判定する際の精度を損なう。
【0139】
好ましくは、方法は、標的と相関性があるマーカー点の座標のみを使用するために選択ステップを含む。
【0140】
選択ステップは、射影演算子
（ここで、
は、
における全射影機のセットである）の使用を含み得る。
【数31】
【数32】
【数33】
【0141】
下記の式（１０）を解くことにより、予測関数Ｆの構成と組み合わされた最適マーカー点の座標の選択を行うことができる。
【数34】
【0142】
選択は、射影演算子πの最適化に対応する。実際に、これは、考慮される各特性（マーカー点の座標）に対する最小化される関数の感度を計算することによって行われる。
【0143】
これにより、標的を判定する際の精度を向上させることができる。
【0144】
パラメータ最適化ステップ
方法は、可能な限りの最も正確な関数を得るために、従って最も正確な（従って最も効果的な）標的を得るために、最適化される必要がある比較的少ない数のパラメータ（正規化係数λ、ガウスカーネルの幅σ及び射影演算子π）に依存する関数Ｆの構成を含む。
方法は、可能な限りの最も正確な関数を得るために、従って最も正確な（従って最も効果的な）標的を得るために、最適化される必要がある比較的少ない数のパラメータ（正規化係数λ、ガウスカーネルの幅σ及び射影演算子π）に依存する関数Ｆの構成を含む。
【0145】
これらのパラメータの最適化を反復的に解く。
【0146】
「リーブワンアウト」型の「交差検証」手順を用いて、パラメータλ、σ及びπを最適化する。
【0147】
「リーブワンアウト」型の交差検証の原理は、１人（患者）を除いて学習データベースにおける例（ここでは患者）の全部を用いて予測関数を構成し、学習データベースから除去された例に対する構成関数の性能を評価することである。全学習データベースの患者の各々に対してこの手順を再生することにより、手順を再生する。
【0148】
一実施形態によれば、「リーブワンアウト交差検証誤差」（「ＬＯＯＥ」とも呼ばれる）は、λ、σ及び演算子πの設定値に対して使用され、下記の式（１１）によって下記のように得られる。
【数35】
【0149】
下記の通りであるため、行列Ｋは、射影演算子πを統合する。
Ｋｉ，ｊ＝ｋ（π（ｘｉ），π（ｘｊ））の場合、Ｋ＝（Ｋｉｊ）ｉ，ｊ＝１．．．ｎ
Ｋｉ，ｊ＝ｋ（π（ｘｉ），π（ｘｊ））の場合、Ｋ＝（Ｋｉｊ）ｉ，ｊ＝１．．．ｎ
【0150】
第１の変型例によれば、進化戦略による大域的最適化のためのアルゴリズムを用いて、２つのパラメータλ及びσ及び演算子πを直接最適化する。このアルゴリズムは、ＣＭＡ－ＥＳ（「共分散行列適応－進化戦略」の略語）方法に基づいており、下記の式（１２）を解くことができる。
【数36】
【0151】
この第１の変型例によれば、方法は、直接最適化に対応する下記のステップを含む。
－ ステップ１：学習データベースを構成／初期化する
○データをロードする
○入力データを定義する（マーカー点ＰＲ）
○出力データを定義する（標的点ＰＣ）
○予備関数Ｆ、ＦＰを構成する
－ ステップ２：パラメータを最適化する
○式（１２）を解く
－ ステップ３：最適パラメータλ、σ及びπ（πは、マーカー点ＰＲの最適座標を統合する）を用いて、統合関数ＦＣを構成する
－ ステップ４：新しい患者について再構成する
○術前ＭＲＩ（例えば、１．５テスラＭＲＩ）を実行する
○ＭＲＩで見える１８個のマーカー点ＰＲ’を判定する
○統合関数ＦＣ及びマーカー点ＰＲ’を用いて、標的点ＰＣ’（標的）を判定する
－ ステップ１：学習データベースを構成／初期化する
○データをロードする
○入力データを定義する（マーカー点ＰＲ）
【数37】
【数38】
－ ステップ２：パラメータを最適化する
○式（１２）を解く
－ ステップ３：最適パラメータλ、σ及びπ（πは、マーカー点ＰＲの最適座標を統合する）を用いて、統合関数ＦＣを構成する
－ ステップ４：新しい患者について再構成する
○術前ＭＲＩ（例えば、１．５テスラＭＲＩ）を実行する
○ＭＲＩで見える１８個のマーカー点ＰＲ’を判定する
○統合関数ＦＣ及びマーカー点ＰＲ’を用いて、標的点ＰＣ’（標的）を判定する
【0152】
この第１の変型例によれば、選択ステップ及び最適化ステップを同時に実行し、射影演算子πを最適化することによって選択を行う。
【0153】
この第１の変型例の実施形態の問題は、１８×３＝５４又は標的の再構成のために両方の半球を結合する場合、１８×３×２＝１０８である射影演算子πのサイズである。この場合、最適化されるパラメータの数は、演算子πに対する可能性の数が、２つのパラメータλ及びσに加えて演算子のサイズに応じて２５４又は２１０８であるという意味で大き過ぎる。
【0154】
特に図５を参照して説明される第２の変型例によれば、２つのパラメータλ及びσ及び演算子πを間接的に最適化する。両方のパラメータλ及びσ及び射影演算子πを最適化することができる反復方法を使用する。
【0155】
まず、下記の式（１３）を解く。
【数39】
【0156】
次に、マーカー点の座標（ｋ）の１つの除去に対するＬＯＯＥ関数の感度を試験する。ｋ番目の座標を除去した場合、下記のように２つの可能な結果が得られて処理される。
－ ＬＯＯＥが増加する場合、標的点及びｋ番目の座標は、相関性があり、この座標を保持する必要があることを意味する。
－ ＬＯＯＥが一定であるか又は減少する場合、標的点及びｋ番目の座標は、相関性がなく、このｋ番目の座標を潜在的に除去できることを意味する。
－ ＬＯＯＥが増加する場合、標的点及びｋ番目の座標は、相関性があり、この座標を保持する必要があることを意味する。
－ ＬＯＯＥが一定であるか又は減少する場合、標的点及びｋ番目の座標は、相関性がなく、このｋ番目の座標を潜在的に除去できることを意味する。
【0157】
潜在的に除去可能な座標の中で、標的と最も相関性が少ない座標を除去する。
【0158】
座標を除去するごとに射影演算子πを（より少ない座標で）更新し、ＬＯＯＥ関数を用いてパラメータσ及びλを最適化する。
【0159】
座標を除去する必要がなくなるまでこの手順を繰り返す。この場合、ＬＯＯＥの意義の範囲内で最適解を保持するために、残りの座標の全部が必要である。
【0160】
この第２の変型例によれば、方法は、間接最適化に対応する下記のステップを含む。
－ ステップ１：学習データベースを構成／初期化する
○データをロードする
○入力データを定義する（マーカー点ＰＲ）
○出力データを定義する（標的点ＰＣ）
○射影演算子πを初期化する
○πが恒等関数に等しいと仮定して、予備関数Ｆ、ＦＰを構成する
－ ステップ２：パラメータを最適化する
○式（１３）でσ及びλの第１の最適化を実行する
○演算子πに存在する座標が（ＬＯＯＥの意義の範囲内で）最適でない間、ループ操作を実行する
●マーカー点の各座標ｋ（ｋ＝１～３×ｍ）に対して、
■座標ｋなしでＬＯＯＥの計算を実行する
■ＬＯＯＥが一定であるか又は減少する場合、学習データベースの座標ｋは、潜在的に除去可能であることに留意されたい
■他の場合、座標ｋを保持し、次の座標ｋ＋１を考える
●潜在的に除去可能な座標の中から座標ｋｍｉｎを除去する。座標ｋｍｉｎは、ＬＯＯＥの意義の範囲内で標的と最も相関性が少ない座標であり、換言すれば、ＬＯＯＥを考慮する場合、最もＬＯＯＥを低下させる
●座標ｋｍｉｎなしで演算子πを更新し、式（１３）でσ及びλの最適化を再実行する
○ｗｈｉｌｅループの終了
－ ステップ３：最適パラメータλ、σ及び演算子π（マーカー点ＰＲの最適座標を統合する）を用いて、統合関数ＦＣを構成する
－ ステップ４：新しい患者について再構成する
○術前ＭＲＩ（例えば、１．５テスラＭＲＩ）を実行する
○ＭＲＩで見える１８個のマーカー点ＰＲ’を判定する
○統合関数ＦＣ及びマーカー点ＰＲ’を用いて、標的点ＰＣ’（標的）を判定する
－ ステップ１：学習データベースを構成／初期化する
○データをロードする
○入力データを定義する（マーカー点ＰＲ）
【数40】
【数41】
○πが恒等関数に等しいと仮定して、予備関数Ｆ、ＦＰを構成する
－ ステップ２：パラメータを最適化する
○式（１３）でσ及びλの第１の最適化を実行する
○演算子πに存在する座標が（ＬＯＯＥの意義の範囲内で）最適でない間、ループ操作を実行する
●マーカー点の各座標ｋ（ｋ＝１～３×ｍ）に対して、
■座標ｋなしでＬＯＯＥの計算を実行する
■ＬＯＯＥが一定であるか又は減少する場合、学習データベースの座標ｋは、潜在的に除去可能であることに留意されたい
■他の場合、座標ｋを保持し、次の座標ｋ＋１を考える
●潜在的に除去可能な座標の中から座標ｋｍｉｎを除去する。座標ｋｍｉｎは、ＬＯＯＥの意義の範囲内で標的と最も相関性が少ない座標であり、換言すれば、ＬＯＯＥを考慮する場合、最もＬＯＯＥを低下させる
●座標ｋｍｉｎなしで演算子πを更新し、式（１３）でσ及びλの最適化を再実行する
○ｗｈｉｌｅループの終了
－ ステップ３：最適パラメータλ、σ及び演算子π（マーカー点ＰＲの最適座標を統合する）を用いて、統合関数ＦＣを構成する
－ ステップ４：新しい患者について再構成する
○術前ＭＲＩ（例えば、１．５テスラＭＲＩ）を実行する
○ＭＲＩで見える１８個のマーカー点ＰＲ’を判定する
○統合関数ＦＣ及びマーカー点ＰＲ’を用いて、標的点ＰＣ’（標的）を判定する
【0161】
図５に例示のように、下記のような他のステップを方法に追加し得る。
－ 本発明によって判定される標的で実行される刺激（より広義には神経外科治療）の効果を評価するステップ５。このステップを治療後の数カ月に実行する。例えば、パーキンソン病に対するＵＰＤＲＳ－３スケール（ＳＴＮ及びＧＰｉ標的）及び本態性振戦に対するＦａｈｎ－Ｔｏｌｏｓａ－Ｍａｒｉｎスケール（ＶＩＭ標的）を使用することができる。
－ 学習データベースに入力するステップ０。特に（本発明による方法又は別の方法によって標的が判定された）術後評価に従って臨床転帰が最適であると考えられている患者の新しい症例を用いて、学習データベースに入力し得る。
－ 本発明によって判定される標的で実行される刺激（より広義には神経外科治療）の効果を評価するステップ５。このステップを治療後の数カ月に実行する。例えば、パーキンソン病に対するＵＰＤＲＳ－３スケール（ＳＴＮ及びＧＰｉ標的）及び本態性振戦に対するＦａｈｎ－Ｔｏｌｏｓａ－Ｍａｒｉｎスケール（ＶＩＭ標的）を使用することができる。
－ 学習データベースに入力するステップ０。特に（本発明による方法又は別の方法によって標的が判定された）術後評価に従って臨床転帰が最適であると考えられている患者の新しい症例を用いて、学習データベースに入力し得る。
【0162】
ステップ５及び／又はステップ０により、治療が効果的であると評価されている新しい患者からのデータで学習データベースを拡充することができる。
【0163】
他の教師あり統計学習方法、例えばカーネルリッジ回帰と異なり、全学習データベースを使用せず、サブセット（サポートベクトル）のみを使用するサポートベクトルマシン（ＳＶＭ）型の方法を使用し得る。学習データベースが「大きい」場合、即ちデータベースにおける患者の数（ｎ）が１００以上である場合、このような方法を使用することができる。この場合、より安価な方法を有することができる。
【0164】
学習データベースが大きい場合、ニューラルネットワーク型の方法を使用することもできる。
【0165】
有利には、１００人の患者の臨床サイズに達するまで新しい患者からのデータで学習データベースを拡充する場合、サポートベクトルマシン型又はニューラルネットワーク型の方法は、関数Ｆ（Ｘ）＝Ｙを構成する教師あり統計学習方法として使用され得る。
【0166】
有利には、学習データベースで複数のクラスを定義するために、統計的分類（教師ありであろうと教師なしであろうと）のステップを導入することもでき、各クラスは、別のクラスのデータと異なる患者データ（例えば、症状、年齢、他の病気などに応じて）のサブセットを含む。次に、各クラスに対して、例えば図５を参照して第２の変型例に記載のように、メタモデル、即ち予測関数Ｆを構成する。
【0167】
統計的分類（教師ありであろうと教師なしであろうと）は、データベースにおける様々な患者のマーカー点の座標に従って実施され得る。データベースの場合、クラスは、同様の座標で患者をグループ化し得、その結果、標的は、同様であり得る。
【0168】
所与の新しい患者について標的を構成するために、実行される必要がある全てのことは、新しい患者がいずれのクラスに属するかを同定し、このクラスに対応する関数Ｆを使用して標的を予測することである。これにより、標的化精度を向上させることができる。
【0169】
下記の方法の１つを用いて統計的分類を実行し得る。
－ ロジスティック回帰
－ ＳＶＭ
－ 階層分類方法
－ ニューラルネットワーク
－ ランダムフォレスト
－ ロジスティック回帰
－ ＳＶＭ
－ 階層分類方法
－ ニューラルネットワーク
－ ランダムフォレスト
【0170】
任意の他の適切な統計的分類方法が使用され得る。
【0171】
提示の様々な実施形態が互いに組み合わされ得る。
【0172】
深部脳刺激以外の方法の任意の用途のために、一般的に、判定される少なくとも１つの正確な脳標的を必要とする任意の神経外科治療に備えて、下記のステップ（学習データベースを構成及び／又は初期化するステップ、学習ステップ、選択ステップ、パラメータ最適化ステップ、再構成ステップ、実施形態に記載の数学的座標系及び／又はマーカー点を定義するステップ）の全部又は一部を使用し得る。
【0173】
更に、本発明は、上述の実施形態に限定されず、むしろ特許請求の範囲の範囲内に入る任意の実施形態に及ぶ。
【0174】
本発明による方法、システム及びコンピュータプログラムは、深部脳刺激のための準備以外の用途を見出し得る。本発明は、電離ビームを当てる領域を判定するために、「ガンマナイフ」型の治療のための準備ステップで使用され得る。更に、本発明は、超音波ビームを当てる領域を判定するために、高強度集束超音波（ＨＩＦＵ）治療のための準備ステップで使用され得る。
【0175】
一般的に、本発明は、判定される少なくとも１つの脳標的を必要とする任意の神経外科治療に備える用途を見出し得る。
【図1】
【図2】
【図3A】
【図3B】
【図3C】
【図3D】
【図3E】
【図3F】
【図4A】
【図4B】
【図4C】
【図4D】
【図5】