(57)【要約】      （修正有）【課題】特に視床下核（ＳＴＮ）等の大脳基底核の領域であるがこれらに限定されない脳領域分析の体積分割方法を提供する。【解決手段】方法は、脳深部刺激処置で関心領域の予測の例として、大脳基底核の皮質下領域内を視覚化し局所化するように機能する。統計学的形状モデルが、高磁場、例えば、７ＴのＭＲ撮像から得られた高品質トレーニングセットでの様々なモードのＳＴＮ又は対応する関心領域及びその予測子に適用される。部分最小二乗回帰（ＰＬＳＲ）法が適用されて、予測すべき領域、例えば、ＳＴＮと、その予測子との空間関係が誘導される。本発明を検証する予測正確性は、手動でセグメント化されたＳＴＮとその予測されたものとの形状類似性並びに位置、サイズ、及び向きの誤差を測定することによって評価される。【選択図】図２
2020028721.tif 000002
【特許請求の範囲】
【請求項１】
電子デバイスを作動させる脳画像パイプライン方法であって、
予測子領域および関心領域を有する患者の脳画像と当該脳画像の特性を示す患者情報とを受信するステップと、
前記患者の脳画像とは異なる複数の脳画像と前記患者の脳画像とは異なる当該複数の脳画像の各々の特性を示すデータベース画像情報とを有するとともに、複数の異なるタイプの画像を含むデータベース（２１４）にアクセスするステップと、ここで、前記データベース画像情報は、撮像モダリティの１つまたは複数を含み、
前記データベース（２１４）から前記患者の脳画像の前記患者情報に一致するデータベース画像情報を有する複数の脳画像の患者固有のトレーニングセットを取得するステップであって、前記トレーニングセットの各脳画像は前記予測子領域とは解剖学的に異なる関心領域に関連付けられた予測子領域を有し、前記トレーニングセットの各脳画像の前記予測子領域は、前記患者の脳画像の前記予測子領域に対応しており、前記トレーニングセットの各脳画像の前記関心領域は、前記患者の脳画像の前記関心領域に対応している、前記複数の脳画像の患者固有のトレーニングセットを取得するステップと、
前記複数の脳画像の患者固有のトレーニングセットを処理して、前記複数の脳画像の前記トレーニングセット内の前記予測子領域の形状、位置、サイズ、または向きに対する前記複数の脳画像の前記トレーニングセット内の前記関心領域の形状、位置、サイズ、または向きの関係に基づいて、予測形状、位置、サイズ、または向きのうちの１つまたは複数を有する予測された関心領域を表す予測子情報を抽出するステップと、
前記予測子情報を用いて前記患者の脳画像を処理して、前記患者の脳画像の前記予測子領域に対して前記患者の脳画像内に前記予測形状、位置、サイズ、または向きのうちの１つまたは複数を有する前記予測された関心領域を組み込むことにより、患者固有アトラスを生成するステップとを含む方法。
【請求項２】
移植挿入後の前記患者の脳を示す移植後画像を受信するステップと、
前記移植後画像を前記患者固有アトラスとマージして、移植後複合画像を作成するステップとをさらに含む、請求項１に記載の方法。
【請求項３】
移植後画像を受信することは、電極を前記患者の視床下核に挿入した後、前記患者の脳を示す画像を受信することを含み、
方法は、前記電極のプログラミングの支援として、前記電極を含む移植後画像を表示するステップをさらに含む、請求項１または２に記載の方法。
【請求項４】
前記患者固有アトラスに、標的および標的に到達するために辿る外科処置の経路のうちの一方または両方を組み込むステップと、
組み込まれた標的および経路のうちの一方または両方を使用して、前記患者固有アトラスを表示するステップとをさらに含む、請求項１乃至３のいずれか一項に記載の方法。
【請求項５】
前記患者の脳画像を処理して前記患者固有のアトラスを生成するステップは、視床下核を含む患者固有のアトラスを生成することを含み、
前記患者固有アトラスに標的および経路のうちの一方または両方を組み込むステップは、視床下核を含む前記患者固有のアトラスに標的および経路の一方または両方を組み込むことを含む、請求項４に記載の方法。
【請求項６】
脳画像の前記患者固有のトレーニングセットを処理するステップは、脳画像の前記患者固有のトレーニングセットを処理して、予測形状および位置を有する予測された関心領域を表す予測子情報を抽出することを含み、
前記患者の脳画像を処理することは、前記予測子情報を使用して前記患者の脳画像を処理して、予測形状および位置を有する前記予測された関心領域を前記患者の脳画像に組み込むことを含む、請求項１乃至５のいずれか一項に記載の方法。
【請求項７】
電子デバイスを動作させる脳画像パイプライン方法であって、
患者の脳画像を受信するステップと、
磁気共鳴画像を含む複数の脳画像を含むデータベースから１つ又は複数の追加の脳画像を受信するステップであって、各追加の脳画像は前記患者の脳画像と異なり、且つ前記患者の脳の解剖学的構造の表現の正確性を最大化するような特徴を有する、当該追加の脳画像を受信するステップと、
前記患者の脳画像及び各追加の脳画像を処理するステップであって、前記患者の脳画像の少なくとも一部分及び各追加の脳画像の少なくとも一部分をマージして患者固有アトラスを生成することを含む、前記脳画像を処理するステップと
を含む方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示の実施形態は、概して医療撮像システムに関し、特に、脳深部刺激（ＤＢＳ：Deep brain stimulation）を含む脳手術に含まれるが、これに限定されない医療処置をサポートするための患者固有のマルチモーダル脳高品質アトラスの生成及び使用に関する。
【背景技術】
【０００２】
脳神経外科処置をガイドするための医療画像の使用は、どこでも見られる。脳は複雑な器官であり、患者の生存及び健康にとって即座に極めて重要であり、撮像技術によらずに把握することが難しい構造及び血管に満ちている。三次元撮像技法、特に利用できるようになった様々なコントラストウェイトを使用した磁気共鳴撮像技法は、脳の様々な内部特徴の区別をより容易にし、更に繊細で安全且つより効率的な処置への道を開く。しかし、現在の一連の撮像技術、特に一般的な臨床設定で利用可能な撮像技術であっても、最良の結果を保証するために医師に必要な情報を描写するには十分ではない。
【０００３】
視床下核（ＳＴＮ：subthalamic nucleus）は、視床の腹側にある脳内の小さな楕円形（レンズ形）の核である。解剖学に関して、視床下核は視床腹部の主要部分であり、大脳基底核系の部分である。多くの場合、医師が利用可能な現在の臨床撮像方法の感度制限等の問題により、パーキンソン病を治療する脳深部刺激の標的として使用されるＳＴＮ等の標的部位を示す画像を得ることは難しい。ＳＴＮは、本明細書において本発明の例として使用されることが多いが、本発明は、ＤＢＳに必須の構造もあれば、他の脳処置用の構造もある他の脳構造にも同様に適用される。
【０００４】
患者の画像の分解能を高めるために、保存された画像又は図を使用して作成されるアトラスとして知られているモデルが照合され、患者の画像にフィッティングされる。しかし、これらのアトラスはそれら自体、制限されることが多く、患者固有ではなく、且つ／又は特定の患者に向けて適合されず、一般に、患者の脳及び特定の処置への有用なマッチングを保証するのに十分な範囲の画像タイプ及び脳の撮像を提供しない。死後サンプリングから導出されるアトラスは、場合によっては高分解能であるが、通常、一人の個人から導出され、それにより、現在評価されている患者の脳の解剖学的構造の正確な表現を反映しない。複数の脳の平均から導出されるアトラスは、特定の患者の解剖学的特徴を描写せず、したがって、正確性がより低い。
【０００５】
大脳基底核の皮質下領域内のＳＴＮは、特にパーキンソン病（ＰＤ：Parkinson’s diseaes）の症状を軽減するための脳深部刺激（ＤＢＳ）手術の極めて重要な標的構造である。それにより、臨床ＭＲＩプロトコルでは見つけにくいそのような小さく複雑な構造の体積分割は、ＰＤの確実なＤＢＳ標的のための必須のプロセスである。ＳＴＮの直接的な視覚化及び位置特定は、高度な高磁場７ＴのＭＲ撮像を用いて促進されるが、そのような高磁場は臨床的に常に利用可能であるわけではない。ＳＴＮは、良好に位置特定する必要があり、臨床設定では容易に利用可能ではない関心領域の例である。
【０００６】
述べたように、脳深部刺激（ＤＢＳ）は、パーキンソン病（ＰＤ）、本態性振戦、及び精神疾患等の神経変性疾患を治療するために広く使用されている脳神経外科的介入である。特に、ＳＴＮ構造のＤＢＳが進行したＰＤに重要な臨床的有効性を有することが報告されている。ＳＴＮ領域内であっても、わずかな配置ずれが深刻な副作用を生じさせるおそれがあるため、ＳＴＮ内への電極の正確な位置決めは、ＤＢＳ処置の成功にとって極めて重要である。
【０００７】
ＤＢＳ処置でのＳＴＮ位置特定に、様々な標的手法が報告されている。それらの方法の幾つかは、アトラスに基づいて交連間距離及び脳室等の解剖学的情報を参照し、ＳＴＮを位置特定する。しかし、そのような間接的な汎用アトラスベースの標的手順では、ＳＴＮの現在使用されているアトラスの適用性を評価するために、個々の被験者にわたるＳＴＮの位置及びサイズのばらつきが、大きな母集団の状況で更に分析される必要がある。
【０００８】
最近の研究では、ＰＤ患者にわたる交連中点に基づいてＳＴＮの位置及びサイズのばらつきを査定することにより、この問題に対処している。より確実な標的手法は、個々の被験者のＭＲＩでのＳＴＮの直接視覚化及び位置特定に基づく。
【０００９】
高磁場強度、例えば７テスラ（Ｔ）のＭＲ撮像技法の近年の進歩により、優れたコントラスト及び高分解能のため、小さく複雑な解剖学的構造の直接識別が可能である。さらに、特に磁化率強調画像（ＳＷＩ：susceptibility-weighted image）及びＴ２強調（Ｔ２Ｗ：T2-weighted）画像等のマルチコントラストＭＲ画像の組合せで７Ｔ-ＭＲＩの恩恵を利用して、被験者固有の３Ｄ構造並びに大脳基底核及び視床内のそれらの構造の結び付きがモデリングされている。不都合なことに、そのような高品質視覚化は、標準の臨床的な１．５Ｔ（又は３Ｔ）のＭＲＩプロトコルを用いて常に可能であるわけではない。位置特定に加えて、ＤＢＳ標的構造の位置、寸法、及び向き等の空間情報を三次元で提供する正確な分割も、最適な電極配置にとって必須である。これは部分的に、ＤＢＳ電極をＳＴＮ（ＰＤの場合）内又は他の疾患の他の部位に不適切に位置決めすることから生じる治療有効性の低下又は副作用の潜在性による。
【００１０】
手動での分割は、時間がかかるとともに、部分的に、臨床の、多くの場合は低磁場のＭＲＩプロトコルでの明確な視覚化の欠如に起因して、解剖学的主観性によって主に行われる。ＳＴＮ構造の解剖学的（又は半解剖学的）分割は、ＳＴＮ構造が、隣接構造と比較的不明確な境界を有する小さく複雑な形状であるため、今もなお問題であるが、上述したように、７Ｔ-ＭＲＩでの優れたコントラスト及び高分解能により、ＳＴＮ構造の位置及び形状を直接視覚化し識別することが可能である。他の脳構造も、臨床設定では位置特定が難しいという同様の特徴を有し、高品質アトラスにより促進することができる。
【００１１】
多くの自動分割技法の中でも、動的形状モデル（active shape model）等の統計学的形状モデルベースに分割は、医療撮像の様々な用途での有効性を示している。これらの手法は、被験者にわたる標的構造のばらつきを統計学的にモデリングし、形状モデルと被験者固有情報とを結合する基準を最小化することにより、最良フィッティングを探す。これらの分割の正確性は、実際の入力（すなわち、被験者）データの初期化及び品質に依存する。さらに、母集団にわたる脳構造の形態学的ばらつきが、統計学的形状モデルを使用して分析された。
【００１２】
近年の研究により、構造の異なる形状間の相関を考慮した、統計学的形状モデルを使用する回帰ベースの形状予測手法が提案されている。皮質下脳構造は、正準相関分析と部分最小二乗回帰（ＰＬＳＲ：partial least squares regression）とを結合することによって予測される。形状予測の幾つかの回帰モデルは、ランドマーク点での不確実性を考慮して比較され、関連する予測子を組み込んで、回帰手法の代わりに疎な観測に基づいて大腿骨及び脛骨予測を更に改善することにより、形状予測の様々な不確実性を有する統計学的形状モデルに基づく関節入力−出力分布モデルが構築された。さらに、予測された標的形状の信頼領域の推定が調べられている。これらの形状予測方法では、関心領域内の限られた情報を有するデータであっても、標的構造を推定することが可能である。
【発明の概要】
【００１３】
本開示は、近傍構造との空間関係を明示的に利用して、本明細書では例として使用されるＳＴＮの位置及び形状を予測するように、上記研究を拡張する。ＳＴＮと解剖学的近接性を有し、７Ｔ-ＭＲＩ又は従来の１．５Ｔ／３ＴのＭＲ撮像でも標準標的プロセスを使用して分割可能な皮質下構造は、ＳＴＮの予測子として見なされる。本発明は、統計学的形状モデルを使用して、７Ｔ-ＭＲデータから得られるトレーニングセットにわたるＳＴＮの形状及び姿勢及びそれらの予測子のばらつきをモデリングする。次に、ＳＴＮとその予測子との空間依存性が、ＰＬＳＲ法によって抽出されるが、そのような関係を利用する他の方法も同様に利用することができる。特定の被験者からの実際の７Ｔ又は１．５Ｔ-ＭＲデータでのＳＴＮの事前分割された予測子の形状及び姿勢を所与として、学習された空間依存性を使用してその被験者のＳＴＮの完全な形状及び姿勢を推定する。ここでも、これは本明細書では主にＳＴＮに関して例示されるが、この手法は他の関心部位に対しても汎用である。
【００１４】
ＳＴＮ等の大脳基底核の構成要素及び他の脳構造の形状、位置、サイズ、及び向きの高分解能自動推定及び予測の方法が開示される。本方法は、好ましくは高分解能且つ高磁場の磁気共鳴（単一コントラスト又はマルチコントラスト）画像のトレーニングセットと、個々の患者のデータセット（臨床磁気共鳴）と、場合によっては予測分析のためのＣＴ画像との受信に基づく。
【００１５】
本方法はまた、トレーニングセット及び患者のセットからの情報を結合して、患者のセット内の構成要素を予測すること、大脳基底核構造又は他の関心領域、変動モード、及びその予測子をトレーニングセットから取得することを含む。さらに、本方法は、統計学的形状モデルを脳関心領域（例えば、大脳基底核構造）、変動モード、及び視床下核予測子に適用すること、部分最小二乗回帰法又は他の統計学的分析ツールを変動モード及びその予測子に適用することを含む。
【００１６】
加えて、本方法は、視床下核の形状及び位置と予測子の形状及び位置との関係を抽出すること、そのような関係を利用して、患者のデータでの関連する形状（例えば、視床下核、ＳＴＮ）及びその位置を予測すること、予測された形状及び位置を患者の臨床画像に組み込むこと、電子媒体を介して、臨床設定に組み込まれた予測を有する患者のデータを転送することを含む。
【００１７】
本方法の変形形態は、予測子が大脳基底核、黒核、赤核、淡蒼球内節、及び視床、又はそれらの組合せに属する場合を含む。
電子デバイスが、データセット及び患者のデータの両方に７テスラ、３テスラ、又は１．５テスラ磁石のうちの少なくとも１つを含む、請求項１に記載の方法。
【００１８】
本発明の関連する組では、少なくとも１つの追加の脳画像は、同じ患者の脳の画像を含む。代替的には、脳画像セット内の少なくとも１つの画像は、Ｔ１強調磁気共鳴画像を含む。幾つかの実施形態によれば、脳画像セット内の少なくとも１つの画像は、Ｔ２強調磁気共鳴画像を含む。幾つかの追加の実施形態によれば、脳画像セット内の少なくとも１つの画像は、磁化率強調磁気共鳴画像を含む。別の実施形態では、脳画像セット内の少なくとも１つの画像は、拡散強調磁気共鳴画像を含む。幾つかの実施形態では、脳画像セット内の少なくとも１つの画像は、機能的磁気共鳴画像を含む。脳画像セット内の少なくとも１つの画像は、別の実施形態では、コンピュータ断層画像を含む。追加の実施形態下では、脳画像セット内の少なくとも１つの画像は、ポジトロン放出断層画像を含む。別の実施形態では、脳画像セット内の少なくとも１つの画像は、３テスラ以上の磁場強度で撮影された磁気共鳴画像を含む。別の実施形態によれば、脳画像セット内の少なくとも１つの画像は、示された脳内の動脈の血管造影画像を含む。脳画像セット内の少なくとも１つの画像は、別の実施形態によれば、示された脳内の血管の静脈造影画像を含む。
【００１９】
別の実施形態は、電子デバイスにより、脳画像のデータベースを保持すること、そのデータベースから少なくとも１つの追加の脳画像を検索することを含む。関連する実施形態では、データベースは、複数のモダリティ及びコントラストを含め、少なくとも２つの画像タイプに属する画像を含み、データベースから脳画像を検索することは、少なくとも２つの画像タイプのいずれが、意図される臨床目的で最良の患者固有アトラスを生成するかを判断すること、データベースからその画像タイプに属する画像を検索することを更に含む。別の実施形態では、データベースは、少なくとも２つの画像タイプに属する画像を含み、データベースから脳画像を検索することは、少なくとも２つの画像タイプのいずれが、患者脳の特定の部分の最良の患者固有アトラスを生成するかを判断すること、データベースからその画像タイプに属する画像を検索することを更に含む。別の実施形態では、データベースから脳画像を検索することは、患者脳画像に関する患者情報を保持すること、データベースに保持される各画像に関する患者情報を保持すること、患者脳画像に関する患者情報に一致する患者情報を有する画像をデータベースから検索することを更に含む。更に別の実施形態では、データベースから脳画像を検索することは、患者画像に一致する画像をデータベースにおいて見つけること、一致する画像を検索することを更に含む。別の実施形態下では、データベースから脳画像を検索することは、患者画像に一致する２つ以上画像をデータベースにおいて見つけること、一致する画像を検索することを更に含む。別の実施形態下では、データベースから脳画像を検索することは、少なくとも１つの電子デバイスにより、データベース内の特定の画像を指定する命令を受信すること、指定された画像を検索することを更に含む。
【００２０】
代替的には、患者固有アトラスを表示することは、患者固有アトラスにおいて関心領域を識別すること、表示されたアトラス内のその関心領域を強調することを更に含む。別の実施形態は、少なくとも２つの関心領域を識別すること、少なくとも２つの関心領域の画像を結合して、患者固有アトラスを作成することを含む。別の実施形態下では、患者固有アトラスが使用されて、外科処置をガイドする。別の実施形態は、画像を使用して、移植処置をガイドすることを含む。追加の実施形態は、移植挿入後の患者の脳を示す少なくとも１つの移植後画像を受信すること、移植後画像を複合画像とマージして、移植後患者固有複合画像を形成すること、移植後複合画像を表示することを含む。更に他の実施形態は、患者固有アトラスを使用して、患者の脳との非侵襲的エネルギー交換を含む処置をガイドすることを含む。
【００２１】
通常、高分解能脳画像パイプラインのシステムであって、システムは、電子デバイスと、電子デバイスで実行され、患者の脳の画像と、少なくとも１つの追加の脳画像とを含む脳画像セットを受信する脳画像検索構成要素と、電子デバイスで実行され、患者の脳画像を少なくとも１つの追加の脳画像とマージして、患者固有アトラスを生成するアトラス生成構成要素と、電子デバイスで実行され、複合画像を表示する表示構成要素とを備える。
【００２２】
本システム及び本方法の他の態様、実施形態、及び特徴は、添付図と併せて検討される場合、以下の詳細な説明から明らかになろう。添付図は概略を目的とし、一定の縮尺で描かれることが意図されない。図中、様々な図に示される同一又は略同一の各構成要素は、単一の番号又は表記で表される。明確にするために、あらゆる構成要素があらゆる図に記されているわけではない。また、当業者が本システム及び本方法を理解できるようにするために示される必要がない場合、システム及び方法の各実施形態のあらゆる構成要素が示されているわけではない。
【００２３】
高分解能脳画像パイプラインの方法が開示され、本方法は、電子デバイスにより、患者の脳の画像と、少なくとも１つの追加の脳画像とを含む脳画像セットを受信すること、電子デバイスにより、患者の脳画像を少なくとも１つの追加の脳画像とマージして、患者固有アトラスを生成すること、電子デバイスにより患者固有アトラスを表示することを含む。
【００２４】
関連する組の実施形態では、少なくとも１つの追加の脳画像は、同じ患者の脳の画像を含む。他の実施形態では、脳画像セット内の少なくとも１つの画像は、Ｔ１強調磁気共鳴画像を含む。幾つかの実施形態によれば、脳画像セット内の少なくとも１つの画像は、Ｔ２強調磁気共鳴画像を含む。幾つかの追加の実施形態によれば、脳画像セット内の少なくとも１つの画像は、磁化率強調磁気共鳴画像を含む。別の実施形態では、脳画像セット内の少なくとも１つの画像は、拡散強調磁気共鳴画像を含む。幾つかの実施形態では、脳画像セット内の少なくとも１つの画像は、機能的磁気共鳴画像を含む。脳画像セット内の少なくとも１つの画像は、別の実施形態では、コンピュータ断層画像を含む。追加の実施形態下では、脳画像セット内の少なくとも１つの画像は、ポジトロン放出断層画像を含む。別の実施形態では、脳画像セット内の少なくとも１つの画像は、３テスラ以上の磁場強度で撮影された磁気共鳴画像を含む。別の実施形態によれば、脳画像セット内の少なくとも１つの画像は、示された脳内の動脈の血管造影画像を含む。脳画像セット内の少なくとも１つの画像は、別の実施形態によれば、示された脳内の血管の静脈造影画像を含む。
【００２５】
別の実施形態は、電子デバイスにより、脳画像のデータベースを保持すること、そのデータベースから少なくとも１つの追加の脳画像を検索することを含む。関連する実施形態では、データベースは、複数のモダリティ及びコントラストを含め、少なくとも２つの画像タイプに属する画像を含み、データベースから脳画像を検索することは、少なくとも２つの画像タイプのいずれが、意図される臨床目的で最良の患者固有アトラスを生成するかを判断すること、データベースからその画像タイプに属する画像を検索することを更に含む。別の実施形態では、データベースは、少なくとも２つの画像タイプに属する画像を含み、データベースから脳画像を検索することは、少なくとも２つの画像タイプのいずれが、患者脳の特定の部分の最良の患者固有アトラスを生成するかを判断すること、データベースからその画像タイプに属する画像を検索することを更に含む。別の実施形態では、データベースから脳画像を検索することは、患者脳画像に関する患者情報を保持すること、データベースに保持される各画像に関する患者情報を保持すること、患者脳画像に関する患者情報に一致する患者情報を有する画像をデータベースから検索することを更に含む。更に別の実施形態では、データベースから脳画像を検索することは、患者画像に一致する画像をデータベースにおいて見つけること、一致する画像を検索することを更に含む。別の実施形態下では、データベースから脳画像を検索することは、患者画像に一致する２つ以上画像をデータベースにおいて見つけること、一致する画像を検索することを更に含む。別の実施形態下では、データベースから脳画像を検索することは、少なくとも１つの電子デバイスにより、データベース内の特定の画像を指定する命令を受信すること、指定された画像を検索することを更に含む。
【００２６】
追加の実施形態によれば、患者固有アトラスを表示することは、患者固有アトラスにおいて関心領域を識別すること、表示されたアトラス内のその関心領域を強調することを更に含む。別の実施形態は、少なくとも２つの関心領域を識別すること、少なくとも２つの関心領域の画像を結合して、患者固有アトラスを作成することを含む。別の実施形態下では、患者固有アトラスが使用されて、外科処置をガイドする。別の実施形態は、画像を使用して、移植処置をガイドすることを含む。追加の実施形態は、移植挿入後の患者の脳を示す少なくとも１つの移植後画像を受信すること、移植後画像を複合画像とマージして、移植後患者固有複合画像を形成すること、移植後複合画像を表示することを含む。更に他の実施形態は、患者固有アトラスを使用して、患者の脳との非侵襲的エネルギー交換を含む処置をガイドすることを含む。
【００２７】
高分解能脳画像パイプラインのシステムも特許請求され、本システムは、電子デバイスと、電子デバイスで実行され、患者の脳の画像と、少なくとも１つの追加の脳画像とを含む脳画像セットを受信する脳画像検索構成要素と、電子デバイスで実行され、患者の脳画像を少なくとも１つの追加の脳画像とマージして、患者固有アトラスを生成するアトラス生成構成要素と、電子デバイスで実行され、複合画像を表示する表示構成要素とを備える。
【００２８】
本システム及び本方法の他の態様、実施形態、及び特徴は、添付図と併せて検討される場合、以下の詳細な説明から明らかになろう。添付図は概略を目的とし、一定の縮尺で描かれることが意図されない。図中、様々な図に示される同一又は略同一の各構成要素は、単一の番号又は表記で表される。明確にするために、あらゆる構成要素があらゆる図に記されているわけではない。また、当業者が本システム及び本方法を理解できるようにするために示される必要がない場合、システム及び方法の各実施形態のあらゆる構成要素が示されているわけではない。
【００２９】
先の概要並びに開示されるシステム及び方法の以下の詳細な説明は、添付図面と併せて読んだ場合によりよく理解されよう。システム及び方法を例示するために、現在好ましい実施形態が図面に示される。しかし、システムも方法も、示される厳密な構成及び手段に限定されないことを理解されたい。
【図面の簡単な説明】
【００３０】
【図１】高分解能脳画像パイプライン方法を示すフローチャートである。
【図２】ＳＴＮ（又は他の脳構造）予測の概略全体図であり、（ａ）高磁場（例えば、７Ｔ）テストデータでの予測の回帰プロセスであり、（ｂ）臨床（例えば、１．５Ｔ）テストデータでの予測の回帰プロセスであり、７Ｔ及び１．５Ｔでの参照データが同じ被験者からスキャンされたことに留意し、（ｃ）予測の枠組みである。
【図３】本明細書に記載されるような電子デバイスの例を示す概略図である。
【図４】本開示におけるネットワークベースのプラットフォームの概略図である。
【図５】開示されるシステムの一実施形態を示すブロック図である。
【図６】標的とされる解剖学的構造が強調された、標的取得のための患者の脳の可能な描写を示す。
【図７】標的とされる解剖学的構造が強調され、移植の経路が示された、標的取得のための患者の脳の別の可能な描写を示す。
【図８】影響を受けた解剖学的構造及び電極それ自体の位置を強調した、ＤＢＳ電極の移植後の患者の脳の可能な描写を示す。
【図９】Ｔ２Ｗ-ＭＲＩ（一人の被験者からの）上の手動分割された構造の３Ｄ形状を示す。（ａ）冠状方向でのＳＮ及びＳＴＮである。（ｂ）冠状方向でのＴｈａ及びＳＴＮである。（ｃ）軸方向でのＧＰｉ及びＳＴＮである。（ｄ）軸方向でのＲＮ及びＳＴＮである。赤はＳＴＮ構造を表し、青、金、緑、及びグレーはＳＮ、Ｔｈａ、ＧＰｉ、及びＲＮ構造をそれぞれ表す。矢印は前方を向く。
【図１０】座標空間でのトレーニングデータ１（ａ）〜５（ｅ）からのＲＮ、ＳＮ、ＳＴＮ、ＧＰｉ、及びＴｈａの対応する３Ｄ表面及びそのランドマーク点を示す。グレー、青、赤、緑、及び黄はＲＮ、ＳＮ、ＳＴＮ、ＧＰｉ、及びＴｈａ構造をそれぞれ表す。各構造の表面に沿って分散する小さな黒点は、表面ランドマーク点である。
【図１１】Ｔ２Ｗ-ＭＲ撮像での１．５Ｔ及び７Ｔの視覚的比較を示す。（ａ）及び（ｂ）は、軸方向での１．５Ｔ及び７ＴのＭＲデータでのＧＰｅ、ＧＰｉ、及びＴｈａ領域をそれぞれ示す。ＳＮ、ＳＴＮ、及びＲＮ領域は、（ｃ）では軸方向平面での１．５Ｔ-ＭＲ撮像の場合が示され、（ｄ）では軸方向平面での７Ｔ-ＭＲ撮像の場合が示される。
【図１２】７Ｔ-Ｔ２Ｗ-ＭＲトレーニングデータセットでの予測子としてＳＮを使用して予測されたＳＴＮの３Ｄ構造及び慣性楕円体を手動分割されたものを重ねた状態で示す。予測されたＳＴＮの３Ｄ形状及び対応する慣性楕円体（セントロイドを有する）は、軸方向平面（（ａ）及び（ｂ））並びに冠状平面（（ｃ）及び（ｄ））においてそれぞれ示される。青は手動分割の３Ｄ形状及び慣性楕円体（セントロイドを有する、白色矢印参照）を表す。赤は予測されたＳＴＮの３Ｄ形状及び慣性楕円体（セントロイドを有する、白色矢印参照）を表す。緑の矢印は前方を向く。この図は、本開示の実施形態の１つの性能の説明に役立つ。
【図１３】ＲＮを使用して予測されたＳＴＮ、トレーニングデータセットにわたるＳＴＮの平均形状、及び１．５Ｔ-Ｔ２Ｗ-ＭＲＩでの手動分割ＳＴＮを示す。予測されたＳＴＮ、ＳＴＮの平均、及び手動分割ＳＴＮの２Ｄ輪郭（青点としてのセントロイドを有する）は、左構造（ａ）及び右構造（ｂ）の軸方向平面（上）、冠状平面（中央）、及び矢状平面（下）においてそれぞれ示される。予測されたＳＴＮ、ＳＴＮの平均、及び手動分割ＳＴＮの３Ｄ表面（ｃ）及び慣性楕円体（セントロイドを有する、白色矢印参照）（ｄ）は、軸方向平面においても示される。青、赤、及び白は、手動分割、予測されたＳＴＮ、及びＳＴＮの平均をそれぞれ表す。緑の矢印は前方を向く。この図は、本開示の実施形態の１つの性能の説明に役立つ。
【図１４】ＲＮを使用して予測されたＳＴＮ、トレーニングデータセットにわたるＳＴＮの平均形状、及び７Ｔ-Ｔ２Ｗ-ＭＲＩでの手動分割されたＳＴＮを示す。この図は、本開示の実施形態の１つの性能の説明に役立つ。
【発明を実施するための形態】
【００３１】
高磁場ＭＲＩ技術の近年の進歩により、ＰＤ等の運動疾患のＤＢＳ標的に重要なＳＴＮ及びＧＰｉ等の皮質下構造の生体内直接視覚化及び位置特定が可能である。三次元でのこれらの領域の視覚化は、外科ＤＢＳ標的及び手術後プログラミングをより確実且つ効率的に提供する。しかし、例えば、ＳＴＮ領域を自動的に線引きすることは、ＳＴＮ領域がＳＮ等の隣接構造間で不明確な分離を有する非常に小さく複雑であるため、難問のままである。特に、標準の低磁場臨床プロトコル下で取得される画像を使用する場合又はデータが領域内に限られた視野を有する場合、そのような直接的な方法でＳＴＮを分割することは可能ではない。
【００３２】
脳処置に関連するＳＴＮ又は他の脳構造の（高分解能）自動形状予測の方法が開示され、ＳＴＮは、ここで及び本明細書全体で、説明のための例として使用される。本方法は、隣接構造へのＳＴＮの空間依存性を予測子として利用し、予測子の幾つかは、標準の臨床処置を用いて容易に位置特定される。
【００３３】
第１に、統計学的形状モデルは、高磁場、例えば、７Ｔ-ＭＲ撮像又は他の高磁場データから得られる１組の高品質トレーニングセットでのＳＴＮ及びその予測子の様々なモードについて開示されている。
【００３４】
第２に、ＳＴＮとその予測子との空間関係を誘導するために、部分最小二乗回帰（ＰＬＳＲ）方法が開示される。そのような関係を誘導するために他の方法も使用可能であり、ＳＴＮはここでは、脳外科、特に脳深部刺激の重要な部位の例として使用されている。
【００３５】
トレーニング時での予測正確性は、手動分割構造（例えば、ＳＴＮ）と予測された構造との形状類似性並びに位置、サイズ、及び向きの誤差を測定することによって評価することができる。これは、所与の関心構造の予測子の定義に使用することができる。
【００３６】
定義。この説明及び添付の特許請求の範囲で使用されるように、以下の用語は、文脈により別段のことが求められる場合を除き、示される意味を有するものとする。
「電子デバイス」は、パーソナルコンピュータ、ラップトップ、タブレット、スマートフォン、及び本明細書に記載される用途をサポート可能な任意の他の電子デバイスを含むものとして定義される。
【００３７】
デバイス又は構成要素は、製品又は手段及びデバイスを１つのマシンとして一緒に動作し得るように関連する場合、電子デバイスに「結合」される。特に、１つの電子機器は、電子デバイスに組み込まれる（例えば、スマートフォンの内蔵カメラ）場合、機器とデバイスとの間で信号を伝搬可能なワイヤによってデバイスに取り付けられる（例えば、コンピュータのポートの１つに差し込まれたワイヤによってパーソナルコンピュータに接続されたマウス）場合、信号を伝搬するワイヤの能力を置換するワイヤレス技術によってデバイスに繋がれる（例えば、モバイル電話のワイヤレスブルートゥース（BLUETOOTH）（登録商標）ヘッドセット）場合、又は複数のマシン間のワイヤレス接続及び有線接続からなる幾つかのネットワークでの共有メンバシップによって電子デバイスに関連する（例えば、コンピュータが及びプリンタをインターネットに接続することができる限り、コンピュータがどこにあるかに関係なく、オフィスに属するコンピュータに文書をプリントするそのオフィス内のプリンタ）場合、電子デバイスに結合される。
【００３８】
「データ入力デバイス」は、データをデバイスに入力するために使用し得るその電子デバイスに結合された全ての機器の総称である。この定義は、限定せずに、キーボード、コンピュータマウス、タッチスクリーン、デジタルカメラ、デジタルビデオカメラ、ワイヤレスアンテナ、全地球測位システムデバイス、オーディオ入力及び出力デバイス、ジャイロスコープ向きセンサ、近接センサ、コンパス、スキャナ、指紋又は網膜スキャナ等の専用読み取りデバイス、及び電磁放射、電磁場、重力、電磁力、温度、振動、又は圧力を検知可能な任意のハードウェアデバイスを含む。
【００３９】
電子デバイスの「手動データ入力デバイス」は、手動操作を使用してユーザがデータを電子デバイスに入力できるようにする、電子デバイスに結合される全てのデータ入力デバイスの集合である。手動入力デバイスは、限定せずに、キーボード、キーパッド、タッチスクリーン、トラックパッド、コンピュータマウス、ボタン、及び他の同様の構成要素を含む。
【００４０】
電子デバイスの「ディスプレイ」は、電子デバイスに結合されるデバイスであり、それにより、電子デバイスは、画像及び分析からの結果を表示することができる。ディスプレイは、限定せずに、画面、テレビデバイス、プリント資料、及びプロジェクタを含む。データを電子デバイスのメモリに「保持」することは、問題となっているアルゴリズムによって要求される検索に好都合な形態でそのデータをそのメモリに記憶すること、及び必要に応じてデータを検索、更新、又は削除することを意味する。
【００４１】
開示される方法は、高磁場（例えば、７Ｔ）ＭＲトレーニングデータから得られるかなり詳細な情報を使用して、ＳＴＮとその予測子−黒質（ＳＮ：substantia nigra）、赤核（ＲＮ：red nucleus）、淡蒼球内節（ＧＰｉ：internal globus pallidus）、及び視床（Ｔｈａ：Thalamus）との学習される空間情報を利用することにより、予測問題、特に、ＳＴＮについての例に対処する。ＳＴＮに高度に相関し、個々の被験者データ（７Ｔ又は臨床１．５Ｔ磁石からの）で容易に線引きすることができる近傍構造は、予測子として利用され、次に、それらの間の空間依存性が使用されて、ＳＴＮの完全形状が予測される。したがって、ＳＮ、ＲＮ、ＧＰｉ、及びＴｈａ構造は、ＳＴＮに隣接し、７Ｔ-ＭＲ撮像で完全に可視であり、７Ｔ-ＭＲデータにおいてＳＴＮを予測する予測子として利用される。しかし、１．５Ｔ-ＭＲ（臨床）撮像では、ＳＮ、ＧＰｉ、及びＴｈａを識別することが難しいため、それでもなお可視であるＲＮが、ＳＴＮの予測子として選ばれる。
【００４２】
提案される開示では、ＳＴＮ予測問題は、統計学的形状モデルと、例としてＰＬＳＲ法を使用して対処される。トレーニングセットにわたる対応する表面点座標ベクトルは、ＳＴＮ及びその予測子を表す。ＳＴＮ及びその予測子の表面点座標の一般化プロクラステス分析により、姿勢パラメータが与えられ、トレーニングセットにわたり各構造を局所的に位置合わせする。次に、カーネル主成分分析が、トレーニングセットにわたるＳＴＮ及びその予測子の位置合わせされた表面点座標に対して実行されて、各形状パラメータをより低次元の空間で表す。
【００４３】
次に、ＰＬＳＲにより、トレーニングセットにわたるＳＴＮ及びその予測子の形状パラメータと姿勢との関係が特定される。実際のデータセット（７Ｔ又は１．５ＴのＭＲＩ）でのＳＴＮの予測子の形状パラメータ及び姿勢を所与として、ＳＴＮの形状パラメータ及び姿勢が、学習された関係を使用して予測される。最後に、ＳＴＮの完全な形状が推定され、予測形状パラメータが元の空間にマッピングされ、予測された姿勢を用いて変換する。ここで、ＳＴＮは、実証例として本明細書の残りの部分で使用される。
【００４４】
図１は、本開示と組み合わせて使用することができる高分解能脳画像パイプライン技法１００を示す。したがって、方法１００は、電子デバイスにより、患者の脳の画像と、少なくとも１つの追加の脳画像（一実施形態では、同じ被験者の、別の実施形態ではデータベースからの）とを含む脳画像セットを受信すること（１０１）を含む。さらに、方法１００は、電子デバイスにより、患者の脳画像を少なくとも１つの追加の脳画像とマージして、患者固有アトラスを生成すること（１０２）を含む。方法１００は、電子デバイスにより患者固有アトラスを表示すること（１０３）も含む。
【００４５】
方法１００は、電子デバイスにより、患者の脳の画像と、少なくとも１つの追加の脳画像とを含む脳画像セットを受信すること（１０１）を含む。幾つかの実施形態では、脳画像検索構成要素４０２がこのステップを実行する。幾つかの実施形態では、少なくとも１つの追加の脳画像は、同じ患者の脳の画像を含む。幾つかの実施形態によれば、脳画像セット内の少なくとも１つの画像は、Ｔ１強調磁気共鳴画像を含む。幾つかの実施形態では、Ｔ１強調ＭＲ画像は１×１×１ｍｍ３又はそれよりも良好な分解能を有し、ここで、良好は、より高い分解能又は改善された信号対雑音比を示す。幾つかの実施形態では、脳画像セット内の少なくとも１つの画像は、Ｔ２強調磁気共鳴画像を含む。幾つかの実施形態でのＴ２強調画像は、０．４×０．４×１ｍｍ３又はそれよりも良好な分解能を有する。幾つかの追加の実施形態によれば、脳画像セット内の少なくとも１つの画像は、磁化率強調（ＳＷＩ）磁気共鳴画像を含む。幾つかの実施形態でのＳＷＩは、０．４×０．４×０．８ｍｍ３又はそれよりも良好な分解能を有する。幾つかの実施形態では、脳画像セット内の少なくとも１つの画像は、拡散強調磁気共鳴画像（ＤＷＩ：diffusion-weighted magnetic resonance image）を含む。幾つかの実施形態でのＤＷＩ画像は、１．５×１．５×１．５ｍｍ３又はそれよりも良好な分解能を有する。ＤＷＩ画像は、トラクトグラフィを使用して強化して、関心のある神経経路を示すこともできる。幾つかの実施形態では、そのようなトラクトグラフィは、全データセットで行われる。幾つかの実施形態では、トラクトグラフィは関心のある領域間で行われる。幾つかの実施形態では、トラクトグラフィは、ＤＢＳ電極の位置から脳全体又は特定の関心領域まで行われる。幾つかの実施形態では、ＤＷＩ画像は、多数の拡散勾配軌跡を用いて撮影され、より正確なトラクトグラフィを作成する。幾つかの実施形態では、ＤＷＩ画像は、５５の別個の勾配軌跡を用いて撮影される。幾つかの実施形態では、ＤＷＩ画像は、５５を超える別個の勾配軌跡を用いて撮影される。幾つかの実施形態では、脳画像セット内の少なくとも１つの画像は、機能的磁気共鳴画像（ｆＭＲＩ）を含む。幾つかの実施形態によるｆＭＲＩは、１．５×１．５×１．５ｍｍ３又はそれよりも良好な分解能を有する。幾つかの実施形態では、上述した分解能は、ＭＲ画像を取得する際により高い強度の磁場を使用することによって達成される。幾つかの実施形態では、脳画像セット内の少なくとも１つの画像は、３Ｔ以上の磁場強度を用いて撮影された磁気共鳴画像を含む。幾つかの実施形態では、磁場は約７テスラであり、６．５テスラ〜７．５テスラの範囲を有する。幾つかの実施形態では、磁場は７テスラよりも高く、例えば、磁場は１１．７テスラであり得る。磁場は１１．７テスラより高くてもよい。
【００４６】
（ここでは実証例としてＳＴＮが使用される）予測枠組みの概観が図２に提示される。本開示は、高磁場（７Ｔ）ＭＲＩ（マグネックス・サイエンティフィック（Magnex Scientific）社、英国）を使用して得られたＭＲＩのトレーニングデータセットを使用するＭＲＩデータ取得を提案する。データセット毎に、Ｔ２Ｗ画像及びＳＷＩが利用されて、トレーニングセット（ＳＴＮ、ＳＮ、ＲＮ、及び全ての関連する構造を含む）にわたり、ＳＴＮと、大脳基底核及び視床領域内の予測子としての近傍構造との手動分割を取得した。さらに、別の組の被験者も１．５Ｔ又は臨床ＭＲＩ（Ｔ１Ｗ及びＴ２Ｗ画像）及び（検証のために同様に）７Ｔ-ＭＲＩの両方でスキャンした。７Ｔ-ＭＲＩ、冠状Ｔ２Ｗ画像、及びＳＷＩも利用して（検証のために）、ＳＴＮ及びその近傍構造を手動分割した。次に、７Ｔ-ＭＲ画像を同じ被験者のＴ１Ｗ１．５Ｔ又は臨床ＭＲデータに位置合わせして、ＳＴＮを予測し、検証した。
【００４７】
各７Ｔ-ＭＲＩの取得プロトコルは、ｉ）パラメータとして、視野（ＦＯＶ）＝２０５×２０５×３６ｍｍ３、分解能０．４×０．４×２．０ｍｍ３、繰り返し時間／エコー時間（ＴＲ／ＴＥ）＝５０００／５７ｍｓｅｃ、フリップ角＝１２０°、帯域幅＝２２０Ｈｚ／ピクセル、位相符号化方向に沿った加速係数＝２（一般化自動較正型部分的並列取得（ＧＲＡＰＰＡ：generalized auto-calibrating partially parallel acquisition））を使用した２Ｄターボスピンエコーシーケンスを用いてのＴ２Ｗ画像取得を含む。合計取得時間は、１つにつき平均約７分であった。このプロトコルを２回繰り返し、中脳の軸方向画像及び冠状画像の両方を得た。また、各７Ｔ-ＭＲＩの取得プロトコルは、ｉｉ）パラメータとして、ＦＯＶ＝１８０×１８０×６０ｍｍ３、分解能０．４×０．４×１．０ｍｍ３、ＴＲ／ＴＥ＝２８／２０ｍｓｅｃ、フリップ角１５°、帯域幅＝１２０Ｈｚ／ピクセル、位相符号化方向に沿った加速係数＝２（ＧＲＡＰＰＡ）を使用した３Ｄフロー補償勾配エコーシーケンスを用いてのＳＷＩ取得を含む。合計取得時間約７分で１つの平均を使用した。このプロトコルも２回繰り返して、中脳の軸方向画像及び冠状画像の両方を得た。
【００４８】
本開示の他の実施形態と組み合わせて、他のプロトコルも同様に使用することができる。
１．５Ｔ-ＭＲＩが、以下の標準臨床プロトコルを用いて取得された。まず、ｉ）画像パラメータとして、ＦＯＶ＝１９２×２５６×１７６ｍｍ３、分解能０．９８×０．９８×１ｍｍ３、ＴＲ＝１６５０ｍｓ、ＴＥ＝３．０２ｍｓ、公称フリップ角１５°、帯域幅１７９Ｈｚ／ピクセルを用い、シーメンス（Siemens）社のＴ１Ｗ-ＡＸ-Ｔ１-ＭＰＲ臨床プロトコルを使用してＴ１Ｗ画像を取得した。また、ｉｉ）画像パラメータとして、ＦＯＶ：１７２×２３０×１９２ｍｍ３、分解能０．７２×０．７２×２．０ｍｍ３、ＴＲ／ＴＥは２５００／２４９ｍｓｅｃ、フリップ角１２０度、帯域幅５３９Ｈｚ／ピクセル、及び２平均を用いて、市販のシーメンス（Siemens）社のＴ２Ｗ-ＡＸ-Ｔ２-３Ｄ-ＳＰＣ（スピンエコー）臨床プロトコルを使用してＴ２Ｗ画像を取得した。
【００４９】
高分解能脳画像パイプライン技法。本開示のシステム及び方法は、開示される用途をサポートする電子デバイスに関するとともに、用途全般の性質に関する以下の観測に鑑みてよりよく理解することができる。例示的な電子デバイスが図３に示される。プロセッサ２００は、専用又は汎用プロセッサデバイスであり得る。関連分野の当業者によって理解されるように、プロセッサデバイス２００は、マルチコア／マルチプロセッサシステム内の単一のプロセッサであってもよく、そのようなシステムは単独で、又はクラスタ若しくはサーバファーム内で動作する計算デバイスのクラスタで動作し得る。プロセッサ２００は、通信基盤２０１、例えば、バス、メッセージキュー、ネットワーク、又はマルチコアメッセージパッシング方式に接続される。
【００５０】
電子デバイスは、ランダムアクセスメモリ（ＲＡＭ）等のメインメモリ２０２も含み、セカンダリメモリ２０３を含むこともできる。セカンダリメモリ２０３は、例えば、ハードディスクドライブ２０４、リムーバブルストレージユニット２０６に接続されたリムーバブルストレージドライブ若しくはインターフェース２０５、又は他の同様の手段を含み得る。関連分野の当業者によって理解されるように、リムーバブルストレージユニット２０６は、コンピュータソフトウェア及び／又はデータを記憶したコンピュータ使用可能記憶媒体を含む。セカンダリメモリ２０３を生成する追加の手段の例としては、プログラムカートリッジ及びカートリッジインターフェース（ビデオゲームデバイスで見られる等）、リムーバブルメモリチップ（ＥＰＲＯＭ又はＰＲＯＭ等）及び関連付けられたソケット、並びにソフトウェア及びデータをリムーバブルストレージユニット２０６からコンピュータシステムに転送できるようにする他のリムーバブルストレージユニット２０６及びインターフェース２０５が挙げられ得る。
【００５１】
電子デバイスは、通信インターフェース２０７を含むこともできる。通信インターフェース２０７は、ソフトウェア及びデータを電子デバイスと外部デバイスとの間で転送できるようにする。通信インターフェース２０７は、モデム、ネットワークインターフェース（イーサネット（Ethernet）（登録商標）カード等）、通信ポート、ＰＣＭＣＩＡスロット及びカード、又は電子デバイスを外部デバイスに結合する他の手段を含み得る。通信インターフェース２０７を介して転送されるソフトウェア及びデータは、信号の形態であり得、信号は、通信インターフェース２０７によって受信可能な電子信号、電磁信号、光学信号、又は他の信号であり得る。これらの信号は、ワイヤ又はケーブル、光ファイバ、電話回線、セルラ電話リンク、及び無線周波数リンク又は他の通信チャネルを介して通信インターフェース２０７に提供し得る。本明細書で考察されるシステム実施形態の通信インターフェースは、有線又は無線２１３に関係なく、データ入力デバイス２０８、デバイスのディスプレイ２１０、及びネットワーク接続との電子デバイスの結合に役立つ。これらの手段のそれぞれが、デバイスそれ自体内に埋め込まれてもよく、ポートを介して取り付けられてもよく、又はブルートゥース（Bluetooth）（登録商標）等のワイヤレス技術を使用して繋がれてもよいことに留意されたい。
【００５２】
コンピュータプログラム（コンピュータ制御論理とも呼ばれる）は、メインメモリ２０２及び／又はセカンダリメモリ２０３に記憶される。コンピュータプログラムは、通信インターフェース２０７を介して受信することもできる。そのようなコンピュータプログラムは、実行されると、プロセッサデバイス２００が後述するシステム実施形態を実施できるようにする。したがって、そのようなコンピュータプログラムは、システムのコントローラを表す。実施形態がソフトウェアを使用して実施される場合、ソフトウェアは、コンピュータプログラム製品に記憶し、リムーバブルストレージドライブ又はインターフェース２０５、ハードディスクドライブ２０４、又は通信インターフェース２０７を使用して電子デバイスにロードし得る。
【００５３】
電子デバイスは、デバイスがアクセス可能なデータベース２１２にデータを記憶することもできる。データベース２１２は、データの任意の構造化された集まりである。本明細書で使用される場合、データベースは「ＮｏＳＱＬ」データストアを含むことができ、これは、既知の組のキー（例えば、アレイインデックス）を使用しての高速検索のために、アレイ等の少数のキー値構造にデータを記憶する。別の可能性は関係データベースであり、これは、データの有用なカテゴリを表すフィールドに記憶データを分割することができる。その結果、記憶されたデータ記録は、データベース２１２内のその既知のデータのカテゴリ内を検索することにより、その記録に記憶されたデータの任意の既知の部分を使用して素早く検索することができ、構造化クエリ言語等の言語を使用してより複雑なクエリによってアクセスすることができ、構造化クエリ言語等の言語は、パラメータとして渡される制限値と、検索されているデータの関係とに基づいてデータを検索する。幾つかのデータベースを検索するために、画像及びメタデータ照合クエリ等のより特殊なクエリの使用も可能である。データベースは任意のデジタルメモリ内に作成することができる。
【００５４】
関係分野の当業者は、任意のデバイスがプロセッサ２００、通信基盤２０１、少なくともメインメモリ２０２、及び通常は通信インターフェース２０７の機能を実行する機能を必ず備えるはずであるが、全てのデバイスが必ずしもこれらの機能を別個に収容するわけではないことにも気付くであろう。例えば、上で定義されたような幾つかの形態の電子デバイスでは、プロセッサ２００及びメモリ２０２は、ニューラルネットと同様に、同じハードウェアデバイスを通して分散することもでき、したがって、通信基盤２０１は、その特定のハードウェアデバイスの構成のプロパティであることができる。多くのデバイスは上述したようなタスクの物理的分割を実施するが、当業者は、物理的構成要素が融合される場合であっても規定通りにタスクの概念的分離を理解しよう。
【００５５】
システムは、スタンドアロン電子デバイス、ネットワーク内で一緒に機能する１組の電子デバイス、又はウェブアプリケーションを含め、幾つかの方法で展開し得る。ウェブアプリケーションをインターネット等のネットワークにわたって機能するように設計された特定の種類のコンピュータプログラムシステムとして当業者は認識しよう。ウェブアプリケーションプラットフォームの概略図は図４に提供される。ウェブアプリケーションプラットフォームは通常、少なくとも１つのクライアントデバイス３００を含み、クライアントデバイス３００は上述したような電子デバイスである。クライアントデバイス３００は、何らかの形態のネットワーク接続を介してインターネット等のネットワーク３０１に接続する。ネットワーク３０１は、電子デバイス３００、３０２に一緒にリンクする任意の構成であり得、限定せずに、電話、ケーブル、及び光ファイバネットワークを含むローカル及び国際有線ネットワーク、セルラ通信及びデータネットワークを含む、電磁放射の信号を使用して情報を交換するワイヤレスネットワーク、並びにそれらの有線及びワイヤレスネットワークの任意の組合せを含む。ネットワーク３０１には、少なくとも１つのサーバデバイス３０２も接続され、サーバデバイス３０２も上述したような電子デバイスである。当然、ウェブアプリケーションが、幾つかのサーバデバイス３０２及び膨大で常時変化する集団のクライアントデバイス３００で実行可能であり、通常、実行されることを当業者は認識し得る。クライアントデバイス３００及びサーバデバイス３０２の両方のコンピュータプログラムは、ウェブアプリケーション３０４の所要機能を実行するように両デバイスを構成する。ウェブアプリケーション３０４は、大量の処理タスクが、ウェブアプリケーションプログラムによってそれらのタスクを実行するように構成されたサーバデバイス３０２により、又は代替的にはクライアントデバイス３００によって達成されるように設計することができる。しかし、ウェブアプリケーションは本質的に、各デバイスにいくらかのプログラミングを含まなくてはならない。
【００５６】
幾つかの態様では、本明細書に記載されるシステム及び方法では、患者固有のデジタル脳アトラスの作成が可能である。アトラスは、その組成に含まれる画像の選択により、現在の臨床状況に基づいて、患者の特定のニーズに合わせ得る。したがって、患者の大脳動脈及び大脳静脈の詳細な画像を含み、ＤＢＳ電極の安全な移植を保証することができる。手術の特定の解剖学的特徴の識別及び標的をより容易にする他の画像が、含まれることもできる。様々な他の脳からの全ての種類の画像のデータベースにより、構造的及び機能的に関連する、他の患者に属する脳のポジトロン放出断層画像（ＰＥＴ）スキャン、機能的ＭＲＩ（ｆＭＲＩ）スキャン、単光子放出コンピュータ断層画像（ＳＰＥＣＴ：single-photon emission computed tomography）スキャン等であるが、これらに限定されない画像を用いて患者固有アトラスを強化し、それにより、患者データを増補することが可能である。データベース画像は、患者脳の特定の画像の非可用性を補償することもできる。
【００５７】
図５は、開示されるシステム４００の幾つかの実施形態を示す。第１の要素は少なくとも１つの電子デバイス４０１である。電子デバイス４０１は、図２及び図３を参照して上で開示した電子デバイスであり得る。電子デバイス４０１は、図３及び図４を参照して上で開示されたように一緒に機能する１組の電子デバイスであり得る。電子デバイスは、アプリケーション３０４のタスクを実行するようにプログラム又は構成される。アプリケーション３０４は、図３及び図４を参照して上述した１つ又は複数のコンピュータプログラムで構成され得る。アプリケーション３０４が実行するタスクは、脳画像検索構成要素４０２、アトラス生成構成要素４０３、及び表示構成要素４０４に概念的に分割し得る。それらの３つの構成要素へのタスクの編成は、実行されるタスクの分類を単に反映したものであり、システム４００の特定の実装形態のアーキテクチャを指示するものではない。
【００５８】
脳画像検索構成要素４０２は、システムによって使用される二次元及び／又は三次元脳画像を取得する。幾つかの実施形態では、脳画像検索構成要素４０２は、画像取得システム４０５と通信する。画像取得システム４０５は、患者の脳をスキャンし、スキャンの結果を捕捉し、それらの結果をデジタル記録に変換し、そのデジタル記録を電子デバイス４０１等の他のデバイスに提供することにより、画像データを取得するマシンである。画像取得システム４０５は、脳の少なくとも一部の画像をレンダリング可能な任意のシステムであり得る。幾つかの実施形態では、画像取得システム４０５は磁気共鳴（ＭＲ）スキャナを含む。幾つかの実施形態では、画像取得システム４０５はコンピュータ断層（ＣＴ）スキャナを含む。幾つかの実施形態では、画像取得システム４０５はポジトロン放出断層（ＰＥＴ）スキャナを含む。他の実施形態では、画像取得システム４０５は超音波スキャナを含む。画像取得システム４０５は高周波超音波（ＨＩＦＵ：high-frequency ultrasound）スキャナを含むこともできる。画像取得システム４０５は光学コヒーレンス断層（ＯＣＴ：optical coherence tomography）スキャナを含み得る。幾つかの実施形態では、画像取得システム４０５は、血管内超音波システム（ＩＶＵＳ：intravascular ultrasound system）を含む。他の実施形態によれば、画像取得システム４０５は、単光子放出コンピュータ断層画像（ＳＰＥＣＴ）システムを含む。画像取得システム４０５の他の実施形態は、プラナーガンマシンチグラフィシステムを含む。
【００５９】
脳画像検索構成要素４０２は、幾つかの実施形態では、電子デバイス４０１がアクセス可能なメモリから脳画像データを検索する。幾つかの実施形態では、脳画像検索構成要素４０２は、電子デバイス４０１がアクセス可能なデータベース２１３と通信する。データベース２１３は、上記の画像取得システム４０５に関連して上述した任意の種類の撮像技術及びプロトコルにより生成される画像を含み得る。データベース２１３は、画像取得システム４０５に関連して上述した画像の任意の組合せを含み得る。データベース２１３は、性別、年齢、及び健康記録等のメタデータを含むこともできる。
【００６０】
システム４００はアトラス生成構成要素４０３も含み、この構成要素は、電子デバイス４０１で実行され、検索された画像全体又は検索された画像の部分を一緒に結合して、図１を参照してより詳細に後述する患者固有アトラスを作成する。
【００６１】
システム４００は表示構成要素４０４も含み、この構成要素は、電子デバイス４０１で実行され、患者固有アトラスを表示する。表示構成要素４０４は、幾つかの実施形態では、電子デバイスに結合されたディスプレイ２１０と通信する。幾つかの実施形態では、ディスプレイ２１０はモニタである。幾つかの実施形態では、表示構成要素４０４は、ゼウス（ZEUS）（登録商標）外科システム等のコンピュータ支援外科ユーザインターフェースシステムと通信し、外科医がロボットアームを使用して外科処置を実行し、これもまた外科医の制御下にあるディスプレイ２１０で手術の進行を見られるようにする。幾つかの実施形態では、表示構成要素４０４は、外科計画に使用される臨床器具と通信する。幾つかの実施形態では、表示構成要素４０４は、脳深部刺激での術後プログラミングに使用されるシステムと通信する。
【００６２】
脳画像検索構成要素４０２は、幾つかの実施形態では、ＭＲ画像に加えて他の形態の画像を受信する。一実施形態によれば、脳画像セット内の少なくとも１つの画像はコンピュータ断層（ＣＴ）画像を含む。他の実施形態では、脳画像セット内の少なくとも１つの画像は、ポジトロン放出断層（ＰＥＴ）画像を含む。更に他の実施形態によれば、脳画像セット内の少なくとも１つの画像は、示された脳内の動脈の画像を含む。幾つかの実施形態では、動脈の画像は動脈の飛行時間磁気共鳴血管造影を含む。脳画像セット内の少なくとも１つの画像は、一実施形態によれば、示された脳内の静脈を示す。一実施形態では、静脈の画像は静脈のＴ２−ｓｔａｒ勾配エコー画像を含む。別の実施形態では、脳画像セット内の少なくとも１つの画像は、超音波画像である。幾つかの実施形態では、脳画像セット内の少なくとも１つの画像は、ＨＩＦＵ画像である。他の実施形態では、脳画像セット内の少なくとも１つの画像はＯＣＴ画像である。他の実施形態では、脳画像セット内の少なくとも１つの画像は、ＩＶＵＳ画像である。他の実施形態によれば、脳画像セット内の少なくとも１つの画像は、ＳＰＥＣＴ画像である。少なくとも１つの画像は、幾つかの実施形態では、プラナーガンマシンチグラフィ画像である。
【００６３】
幾つかの実施形態は、電子デバイスにより、脳画像のデータベースを保持すること、そのデータベースから少なくとも１つの追加の脳画像を検索することを含む。幾つかの実施形態では、データベース内の画像は上述した任意の画像を含む。幾つかの実施形態では、データベース内の画像は、画像が表す脳に従って編成される。幾つかの実施形態では、示される各脳に対応する画像は、より詳細に後述する方法に従って結合されて、三次元アトラスを形成する。幾つかの実施形態では、データベースは、少なくとも２つの画像タイプに属する画像を含み、脳画像をデータベースから検索することは、画像タイプのいずれが、意図される臨床目的で最良の患者固有アトラスを生成するかを判断すること、データベースからその画像タイプに属する画像を検索することを更に含む。「画像タイプ」は、本明細書で使用される場合、特定の撮像技法によって生成される画像のカテゴリである。例えば、１つの画像タイプはｆＭＲＩ画像であり、別の画像タイプはＴ１強調画像である。幾つかの実施形態では、各画像タイプの部分が結合される。幾つかの実施形態では、画像タイプ全体が結合される。幾つかの実施形態では、患者の脳の画像で明確には見えない患者の脳の特定の特徴又は領域は、より詳細に後述するように、欠けた特徴又は領域を強調する可能性が高い特定の画像タイプに属するデータベースからの画像を含むように患者固有アトラスが作成される場合に最も可視である。このステップを参照して上述した全ての画像タイプは、限定なしで、高強度又は低強度磁場を用いて撮影されるＴ１、Ｔ２、ＳＷＩ、ＦＭＲＩ、及びＤＷＩ画像を含め、データベースに含まれ得る。患者の脳に関して以下により詳細に記載するように、データベースに示される脳の移植後画像を含むこともできる。
【００６４】
幾つかの実施形態では、脳画像検索構成要素４０２は、患者脳画像に関する患者情報を保持し、データベースに保持される各画像に関する患者情報を保持し、患者脳画像に関する患者情報に一致する患者情報を有する画像をデータベースから検索することにより、データベースから画像を検索する。幾つかの実施形態では、患者脳画像及びデータベース画像に関する患者情報は、人口統計学的情報である。人口統計学的情報は、脳が画像にある人物の年齢であり得る。人口統計学的情報は、脳が画像にある人物の性別であり得る。幾つかの実施形態では、患者情報は、示された脳に関する測定データを含む。幾つかの実施形態では、測定データは示された脳のサイズを記述する。他の実施形態では、測定データは、示された脳の特徴の寸法を記述する。例えば、患者の脳に関して記録された患者データは、脳室の寸法を記述し得る。幾つかの実施形態では、患者情報は病歴データである。病歴データは、特定の疾患の診断を含み得る。診断された疾患は、撮像された脳の形態に影響する疾患であり得る。診断は、疾患の進行度を記述することもできる。診断は、疾患に対応する撮像された脳内の特定の特徴を記述し得る。病歴データは、治療履歴を記述することもできる。治療履歴は、治療の過程で試みられた処置を含み得る。治療履歴は、試みられた処置の結果を含むこともできる。幾つかの実施形態では、患者情報は、撮影された画像に関するデータである。データは、撮影された画像タイプを含み得、例えば、データは、所与の人物の脳が幾つかのＭＲタイプ及びＰＥＴスキャンを使用して撮像されたことを示し得る。データは、撮影された画像の分解能レベルを示し得る。データは、画像の取得に使用される取得パラメータを記述することもできる。
【００６５】
幾つかの実施形態では、照合プロセスは、患者情報の幾つかの項目を連続して探すことを含む。例えば、脳画像検索構成要素４０２はまず、患者と同じ年齢範囲内の人物に対応するデータベース内の全ての画像を探し、次に、患者と同じ疾患を有する人物に対応する全ての画像の結果セット内を探し得る。そして、その結果セットは、患者の脳の利用可能な画像を補足する可能性が高い画像タイプを更に探し得る。幾つかの実施形態では、探すことは、患者の画像に使用されたタイプと同じタイプのデータベース内の画像を使用することを含み、例えば、患者データが１．５テスラ（又は３テスラ）ＭＲＩを含む場合、探すことは、データベース内の一致する１．５テスラ（又は３テスラ）ＭＲＩ画像を探すことに基づき得る。
【００６６】
データベース内の画像は、線形法及び非線形法を使用して座標系に位置合わせし得る。座標系は患者画像の座標系であり得る。座標系はデータベース画像の座標系であり得る。座標系は外科処置によって指定され得る。位置合わせは線形であり得、位置合わせプロセスが、線形変換系を使用して記述し得るマッピングに画像を位置合わせすることを意味する。線形変換は並進移動を含み得る。線形並進移動は全体画像のスケーリングを含み得る。線形並進移動は画像の回転を含み得る。線形変換は、アフィンレジストレーションを含み得る。位置合わせは非剛体であり得、一次方程式系に関して記述することができず、元画像の比率を必ずしも保存しないことを意味する。非剛体レジストレーションプロセスは、画像の異なる部分に異なる線形変換を受けさせ得、領域間を補間して、画像全体の完全性を維持する。非剛体レジストレーションプロセスは、画像を弾性変形させ得る。非剛体変換は、位置合わせされている画像内の特徴又はランドマークを他の画像の同様の特徴と並べるために使用され得る。ランドマークはユーザ提供であり得る。ランドマークは解剖学的であり得る。ランドマークは自動的に発見され得、例えば、ランドマークは、自動大脳溝検出手順によって発見され得る。位置合わせは、線形レジストレーションと非剛体レジストレーションの組合せであり得、例えば、線形レジストレーションは、まず、画像を適切な向きで座標系内に配置するために使用され得、非剛体レジストレーションが続き、画像の特定の特徴を標準位置に位置合わせし得る。線形レジストレーションは、将来の非線形レジストレーションのためにデータベースのサブセットを選択するために使用することもできる。レジストレーションはユーザガイドされ得る。幾つかの実施形態では、ユーザは、座標系内で画像を回転させるように電子デバイスに指示する。幾つかの実施形態では、ユーザは、座標系内で画像のサイズを変更するように電子デバイスに指示する。幾つかの実施形態では、ユーザは、撮像された脳の少なくとも１つの特徴を撮像された別の脳の特徴と位置合わせするように電子デバイスに指示する。位置合わせは、自動プロセス及び手動プロセスの組合せを含み得る。例えば、ある画像の位置合わせは、画像の特定の特徴を別の画像の特定の特徴に並べるように自動的に変更し得、それに続けて、ユーザは、更なる線形又は非剛体レジストレーション変更を行い、ユーザが満足するまで画像を並べるように電子デバイスに指示し得る。位置合わせは、平均二乗誤差又は相互情報等のコストを最小化するように行うことができる。位置合わせは、既存のアルゴリズム又は適合された技法を用いて実行することができる。
【００６７】
幾つかの実施形態では、脳画像検索構成要素４０２は、患者画像に一致するデータベース内の画像を見つけ、一致する画像を検索することにより、データベースから画像を検索する。幾つかの実施形態では、照合は、データベース内の同じモダリティ及びコントラストの全ての画像の中の所与のコストを最小化する、データベース内の画像によって決定される。幾つかの実施形態では、コストは関心領域のみで計算される。幾つかの実施形態では、コストは、異なるモダリティ及びコントラストについて独立して計算され、それぞれ１つに、最小化画像が選択される。幾つかの実施形態では、画像選択は線形レジストレーションに基づく。幾つかの実施形態では、画像選択は非線形レジストレーションに基づく。幾つかの実施形態では、画像選択は、線形レジストレーション及び非線形レジストレーションの組合せに基づく。幾つかの実施形態では、脳画像検索構成要素４０２は、患者の脳の類似特徴に一致する特徴を有する脳画像を見つける。例えば、脳画像検索構成要素４０２は、患者の脳内の脳室と略同じサイズ及び形状の脳室を有する画像を見つけ得る。脳画像検索構成要素４０２は、患者の脳内の類似特徴と略同じ幾つかの特徴を有する画像を見つけ得る。脳画像検索構成要素４０２は、患者情報が患者の画像の患者情報に一致する画像の先の検索の結果セット内でこの画像照合検索を実行することもできる。脳画像検索構成要素は、２つ以上の一致結果を生成することができる。幾つかの実施形態では、データベース内の同じ被験者に対応する全てのデータが、検索が照合にそのサブセットを使用した場合であっても、患者固有アトラスに使用される。幾つかの実施形態では、患者からの異なる画像がデータベースの検索に使用され得る。例えば、患者の脳のＴ１強調画像が、最良一致Ｔ１強調画像の検索に使用され得、一方、患者の脳のＴ２強調画像が、最良一致Ｔ２の検索に使用され得る。
【００６８】
他の実施形態では、データベースから脳画像を検索することは、少なくとも１つの電子デバイスにより、データベース内の特定の画像を指定する命令を受信すること、指定された画像を検索することを更に含む。幾つかの実施形態では、ユーザはユーザインターフェースを通して命令を入力する。ユーザインターフェースは、ユーザが、画面上でデータベースからの画像を閲覧し、患者画像と結合する画像を選択できるようにし得る。ユーザインターフェースは、ユーザが、データベースから幾つかの画像を組み合わせて閲覧し、患者画像と結合する組合せを選択できるようにし得る。ユーザは、ユーザインターフェースを使用して、患者情報照合結果セット内の画像を閲覧し得る。ユーザは、ユーザインターフェースを使用して、画像照合結果セット内の画像を閲覧し得る。ユーザは、任意の形態の手動データ入力デバイスを使用してこれらの指示を入力し得る。幾つかの実施形態では、ユーザインターフェースはカスタマイズされる。幾つかの実施形態では、ユーザインターフェースは既存の外科プラットフォームの１つである。
【００６９】
加えて、方法１００は、電子デバイスにより、患者の脳画像を少なくとも１つの追加の脳画像とマージして、患者固有アトラスを生成すること（１０２）を含む。幾つかの実施形態では、患者画像及び少なくとも１つの追加の脳画像が両方とも同じ座標系に位置合わせされる場合、アトラス生成構成要素４０３は、同じ座標系で一方の画像を他方の画像に重ねることによって２つの画像をマージし得る。マージするプロセスは、幾つかの実施形態では、別の画像に一致するように、一方の画像の位置合わせを変更することを含む。位置合わせの変更は自動であり得る。幾つかの実施形態では、自動位置合わせ変更は、２つの画像間の向きの差を検出することと、一方の画像を線形変換させて、画像を回転させることとを含む。幾つかの実施形態では、自動位置合わせ変更は、２つの画像間のスケールの差を検出することと、一方の画像を線形変換させて、画像をスケーリングすることを含む。幾つかの実施形態では、自動位置合わせ変更は、２つの画像間のサイズ差を検出することと、一方の画像を線形変換させて、画像をスケーリングすることとを含む。幾つかの実施形態では、自動位置合わせ変更は、２つの画像間の分解能差を検出することと、一方の画像を線形変換させて、画像をスケーリングすることとを含む。幾つかの実施形態では、自動位置合わせ変更は、２つの画像間の位置差を検出することと、一方の画像を線形変換させて、画像を移動させることとを含む。他の実施形態では、自動位置合わせ変更は、２つの画像間の向きの差を検出することと、一方の画像を非剛体変換させて、他方の画像に合わせることとを含む。非剛体変換は、両画像に存在する少なくとも１つの解剖学的特徴の検出及び一致によりガイドし得る。
【００７０】
位置合わせ変更は手動であり得る。幾つかの実施形態では、ユーザは、座標系内で画像を回転させるように電子デバイスに指示する。幾つかの実施形態では、ユーザは、座標系内で画像のサイズを変更するように電子デバイスに指示する。幾つかの実施形態では、ユーザは、撮像された脳の少なくとも１つの特徴を撮像された別の脳の特徴と位置合わせするように電子デバイスに指示する。位置合わせ変更は、自動プロセス及び手動プロセスの組合せを含み得る。例えば、一方の画像の位置合わせは、画像の特定の特徴が別の画像の特定の特徴と並ぶように自動的に変更し得、それに続き、ユーザは、更なる線形又は非剛体レジストレーション変更を行い、ユーザが満足するまで画像を並べるように、電子デバイスに指示し得る。ユーザは、任意の形態の手動データ入力デバイスを使用してこれらの指示を入力し得る。
【００７１】
幾つかの実施形態では、アトラス生成構成要素４０３は、示された脳の特定の特徴を識別する。幾つかの実施形態では、アトラス生成構成要素は、識別された特徴又は領域に対応する１組のボクセルを記録し、その情報を患者固有アトラス又は患者固有臨床画像に組み込む。幾つかの実施形態では、アトラス生成構成要素４０３は、脳画像セット内の少なくとも１つを選択し、少なくとも１つの選択された画像内の関心領域を識別し、関心領域を除く全ての画像データを少なくとも１つの選択された画像から除去することにより、患者固有データを増補する。例えば、アトラス生成構成要素４０３は、データベースからの画像を患者画像に位置合わせした後、ＳＴＮ又は関連する形状以外の全てをデータベース画像からなくし、そうして変更された画像のみを患者画像と結合して、患者固有アトラスを生成し得る。図６は、脳の特定の構造を強調することによって作成された画像を示す。表示された画像内の脳の大半は、関心領域の背後及び下の２つの断面平面図５０１に制限される。三次元でレンダリングされるのは関心部位であり、関心部位はＳＴＮ５０２、尾状核５０３、及び被殻５０４を含む。脳の三次元構造の残りはこの画像では見えないため、関心部位は、画像の視角から明確に見える。
【００７２】
再び図１を参照すると、方法１００は、電子デバイスにより、患者固有アトラスを表示すること（１０３）も含む。幾つかの実施形態では、表示構成要素４０４は、示された脳の特定の角度からのビューをシミュレートするように、電子デバイス４０１が利用可能な二次元ディスプレイに患者固有アトラスを投射する。幾つかの実施形態によれば、表示構成要素４０４は、患者固有アトラス内の関心領域を識別し、図６に示されるように、表示されたアトラスにおいてその関心領域を強調する。幾つかの実施形態では、投射は、示された脳内の特定の構造のビューを提供し、例えば、図６によって表される二次元投射は、関心形状の例として、ＳＴＮ５０２の特定の角度からのビューを提供する。幾つかの実施形態では、投射は、示された脳内の２つ以上の構造のビューを提供する。幾つかの実施形態では、示された各解剖学的特徴は、異なる色で表示されて、示された他の解剖学的特徴（図示せず）から区別される。例えば、ＳＴＮ５０２又は特定の関心形状は、黄色で表示され得、一方、尾状核５０３はターコイズ又は施術者に有用な任意の色で表示され得る。幾つかの実施形態では、ユーザは、投射によって示される視角を変更するように、表示構成要素４０４に指示する命令を入力することができる。幾つかの実施形態では、表示構成要素４０４は、特定の解剖学的特徴を表示するように指示する命令を受信する。幾つかの実施形態では、表示構成要素４０４は、特定の解剖学的特徴を特定の色で描写するように指示する命令を受信する。幾つかの実施形態では、表示構成要素４０４は、外科標的を表示する。外科標的は、周囲と対比するために別個の色で表示され得る。幾つかの実施形態では、表示構成要素４０４は、外科標的への経路を表示する。幾つかの実施形態では、経路は、血管の破損を避けるように計算される。幾つかの実施形態では、経路は、神経構造を避けるように計算される。
【００７３】
方法の幾つかの実施形態は、医療処置の標的を見つけることを含む。例として、図６のＳＴＮ５０２の描写は、医療処置でＳＴＮ５０２を標的とするために使用される１つのビューであり得る。幾つかの実施形態では、処置は外科処置である。幾つかの実施形態では、表示構成要素４０４は、処置前に、外科処置の標的を表示する。幾つかの実施形態では、表示構成要素４０４は、処置前に、標的に到達するために辿る外科処置の経路を表示する。例えば、図７の画像は、ＤＢＳ電極挿入の解剖学的標的６０１と、辿るべき電極挿入の挿入経路６０２とを示す。表示構成要素４０４は、幾つかの実施形態では、手術中の標的の更新画像を定期的に表示する。幾つかの実施形態では、表示構成要素４０４は、手術中の経路の更新画像を定期的に表示する。幾つかの実施形態では、表示構成要素４０４は、手術中の標的の更新画像を連続表示する。幾つかの実施形態では、表示構成要素４０４は、手術中の経路の更新画像を連続表示する。幾つかの実施形態では、表示構成要素４０４は、内視鏡処置で使用される内視鏡画像を補足する画像を表示する。幾つかの実施形態では、表示構成要素４０４は、腹腔鏡手術を支援する画像を表示する。幾つかの実施形態では、表示構成要素４０４は、血管内画像誘導処置を支援する画像を生成する。表示構成要素４０４は、幾つかの実施形態では、ロボット手術を誘導する画像を提供する。幾つかの実施形態では、ロボット手術は、外科医が、外科ツールを直接制御する電子デバイスに結合されたデータ入力デバイスを操作することにより、幾つかの外科ツールを間接的に制御する処置である。幾つかの実施形態では、ロボット手術は、手術部位及び周囲の解剖学的構造で何が行われているかを外科医に示す表示を組み込む。
【００７４】
処置は、物質を患者の脳内に移植することを含み得る。幾つかの実施形態では、患者の脳内に移植される物質は生体組織である。幾つかの実施形態では、患者の脳内に移植される物質は１つ又は複数の生体細胞を含む。幾つかの実施形態では、患者の脳内に移植される物質は遺伝子材料である。遺伝子材料は、患者の脳内の細胞に直接配置し得る。遺伝子材料は、仲介エージェントにより患者の脳に送達され得る。仲介エージェントはマシンであり得る。仲介エージェントは、「ナノボット」等の小型ロボットデバイスであり得る。仲介エージェントは自律的に機能し得る。仲介エージェントは生物であり得る。仲介エージェントは原核単細胞生物であり得る。仲介エージェントはウィルスであり得る。仲介エージェントは真核単細胞生物であり得る。仲介エージェントは多細胞生物であり得る。仲介エージェントは化合物であり得る。仲介エージェントはタンパク質であり得る。仲介エージェントは、１組のタンパク質である。幾つかの実施形態では、移植処置は、幹細胞を挿入する処置である。幾つかの実施形態では、処置は、移植物を患者の脳内に配置することを含み得る。幾つかの実施形態では、表示構成要素４０４は、移植前に、患者の脳内の移植標的を示す画像を生成する。幾つかの実施形態では、表示構成要素４０４は、移植前に患者の脳内の移植経路を示す画像を生成する。幾つかの実施形態では、表示構成要素４０４は、移植中、移植標的の更新画像を定期的に表示する。幾つかの実施形態では、表示構成要素４０４は、移植中、移植経路の更新画像を定期的に表示する。幾つかの実施形態では、表示構成要素４０４は、移植中、移植標的の更新画像を連続表示する。幾つかの実施形態では、表示構成要素４０４は、移植中、移植経理の更新画像を連続表示する。幾つかの実施形態では、移植処置は、ＤＢＳ電極及びリードを挿入する処置である。
【００７５】
幾つかの実施形態では、方法は、移植の挿入後の患者の脳を示す少なくとも１つの移植後画像を受信することと、移植後画像を患者固有構築アトラスとマージして、移植後複合画像を形成することと、移植後複合画像を表示することとを更に含む。移植は、脳内に完全又は部分的に挿入される任意の異物であり得る。移植物は、限定せずに、ＤＢＳ電極、ＤＢＳリード、マイクロチップ、ニューロンと相互作用する電極、無線周波数送信器、無線周波数受信器、感覚代行用の神経移植、及び脳幹人工内耳を含み得る。移植後画像はＣＴ画像であり得る。移植後画像は、図４を参照して上述したような任意の他の医療画像であり得る。幾つかの実施形態では、アトラス生成構成要素４０３は、ステップ１０２を参照して上述した技法を使用して、移植後画像を患者固有アトラスと結合する。幾つかの実施形態では、アトラス生成構成要素４０３は、移植後画像を任意の移植前患者データと位置合わせすることにより、結合画像を作成する。幾つかの実施形態では、アトラス生成構成要素４０３は、移植後画像を対応する患者固有アトラスと位置合わせすることにより、結合画像を作成し、そのような結合を行う。幾つかの実施形態では、表示される電極位置は、移植の再位置決めする作業をガイドするために使用される。幾つかの実施形態では、リアルタイムＣＴ撮像は患者固有アトラスと結合されて、移植物の初期挿入及び再位置決めをガイドする。幾つかの実施形態では、表示される電極位置は、ＤＢＳの最適恩恵に最も関連する電極の１つ又は複数の接点を決定するような電極のプログラムの支援に使用される。図８は、一実施形態により生成することができる結合画像を示し、ＤＢＳ電極７０１及びリード７０２の捕捉位置が、電極７０１が移植された患者のＳＴＮの画像（又はアトラスからの対応する画像）７０３に重ねられる。
【００７６】
幾つかの実施形態では、表示構成要素４０４は、移植のない状態の患者の脳と共に、エネルギーを交換する処置の標的を表示する。幾つかの実施形態では、表示構成要素４０４は放射線治療の標的を表示し、例えば、表示構成要素４０４は、放射線治療の標的の主要を表示し得る。幾つかの実施形態では、表示構成要素４０４は放射線手術の標的を表示する。放射線手術は、少数の別個の放射線量の放射線を定位的に患者の脳内の特定の標的に集束させることを含み得る。幾つかの実施形態では、表示構成要素４０４は、集束超音波手術の標的を表示する。集束超音波手術は、幾つかの実施形態では、高強度集束超音波ビームを使用して、標的組織を破壊する。
【００７７】
開示されるシステム及び方法の実施形態は、神経外科医が、はるかに洗練され正確な手術及びＤＢＳ移植処置に従事できるようにする。特に、血管及び重要な神経構造を回避する能力並びにＤＢＳ移植のより正確な初期標的及び経路を選択する能力は、ＤＢＳ移植に必要な時間及び合併症率を低減する潜在性を有する。開示されるシステム及び方法は、ロボット手術等の自動処置に従事する能力を強化する。そのような状況では、患者固有アトラスは、ロボットシステムに標的及び全侵入経路を指示し得る。遠隔手術もより高品質でより柔軟な撮像技法の使用から恩恵を受ける。白質等の撮像が難しい組織の描写を改善する能力により、神経外科医は従来は届かなかった組織を標的とすることができる。
【００７８】
［ＳＴＮ位置予測方法］
事前処理。大脳基底核及び視床領域内の皮質下構造−ＳＮ、赤核（ＲＮ）、淡蒼球内節（ＧＰｉ）、及び視床（Ｔｈａ）は、ＳＴＮへの隣接性及び分割するための可視性を考慮して、７Ｔ-ＭＲＩでのＳＴＮの潜在的な予測子として選ばれる。１．５Ｔ（又は臨床）ＭＲＩデータでは、唯一検出可能なＲＮ構造が、ＳＴＮの予測子として使用される（その他の構造は、典型的な臨床撮像で識別することは難しいことがあり、識別可能な場合には使用することができる）。
【００７９】
高磁場（７Ｔ）ＭＲトレーニングデータセット上のＳＮ、ＳＴＮ、ＲＮ、ＧＰｉ、及びＴｈａの左側及び右側は、被験者固有分析パイプラインに従ってコンピュータプログラム又は解剖学専門家によって手動分割される。７Ｔ-ＭＲＩデータセットの軸方向及び冠状Ｔ２Ｗ内のＳＮ、ＳＴＮ、ＲＮ、ＧＰｉ、及びＴｈａの３Ｄ分割は、図９のアミラ（Amira）ソフトウェア（マーキュリーコンピュータシステム（Mercury Computer System）社）を用いて３Ｄ視覚化され、このステップにおいて、市販又は組織内のＧＵＩ及び視覚化ツールが使用可能である。ＧＰｉ、ＲＮ、及びＴｈａはここでは、高分解能ＳＷＩ及びＴ２Ｗの軸方向で分割されたため、ＳＮ及びＳＴＮは、高分解能ＳＷＩ及びＴ２Ｗの冠状方向において分割されたが、各データセット上のＳＮ及びＳＴＮはまず、ＦＳＬ-ＦＬＩＲＴ（ＦＭＲＩＢの線形画像レジストレーションツール）を使用して軸方向ＳＷＩに位置合わせされ、全ての構造を共通座標空間に変換する。この位置合わせに他の市販又は組織内のソフトウェアパッケージが使用可能であり、位置合わせに他の方向が使用可能である。
【００８０】
臨床（１．５Ｔ）ＭＲデータセットでの予測テスト及び検証のために、ＲＮ及びＳＴＮは手動で、対応する７Ｔ-ＭＲデータセットで取得され、ＦＳＬ-ＦＬＩＲＴ（ここでも、他の位置合わせ技法が使用可能である）を使用して軸方向Ｔ２Ｗデータに位置合わせされた冠状ＳＷＩ及びＴ２Ｗの組合せ上で分割された。次に、位置合わせされたＲＮ及びＳＴＮ構造は、７Ｔでの軸方向Ｔ２Ｗを１．５Ｔでの軸方向Ｔ１Ｗに位置合わせすることにより得られる変換行列を使用して、１．５ＴでのＴ１Ｗデータに位置合わせされ、７Ｔ-ＭＲＩ及び１．５Ｔ-ＭＲＩでは同じ位置合わせプロトコルを用いた。７Ｔ-ＭＲＩ撮像が利用可能ではない実際の臨床の場合、ＲＮ構造は、１．５Ｔ-ＭＲテストデータで取得した軸方向Ｔ２Ｗ上で自動的又は手動で分割し、次に、一対の対応する７Ｔ-ＭＲ画像を有する１．５Ｔ参照データに位置合わせすることができる（図２（ｂ）参照）。
【００８１】
テストデータでの予測子と７Ｔ-ＭＲトレーニングデータでの予測子との間での姿勢の変動（例えば、スキャン中の患者の異なる位置及び向き並びに異なるＭＲモダリティ又は磁場間での分解能不一致）が、不正確な予測結果を生じさせるおそれがあるため、高磁場（７Ｔ）で得られた各トレーニングセットのＳＴＮ及び予測子が、テストデータ（７Ｔ又は１．５Ｔ）と同じ座標空間に正規化される必要があることに留意されたい。そのために、相互情報（ＭＩ）類似性尺度を用いるアフィン線形レジストレーションを使用して、７Ｔ-ＭＲトレーニングセットの座標空間及び位置合わせされた軸方向Ｔ２Ｗ画像を７Ｔ-ＭＲテストデータの軸方向Ｔ２Ｗ画像に正規化した。次に、７Ｔ-ＭＲトレーニングセット上のＳＴＮ及び予測子は、７Ｔ-ＭＲトレーニングセットとテストセットとの間で得られた変換行列（Ｍ＿（ｉ，７Ｔ））を使用して、７Ｔ-ＭＲテストデータに位置合わせされた（図２（ａ）参照）。
【００８２】
さらに、１．５Ｔ-ＭＲＩでの予測の場合、７Ｔ-ＭＲトレーニングデータのＳＴＮ及びＲＮ（７Ｔ-ＭＲテストデータ上に位置合わせされた）並びに７Ｔ-ＭＲテストデータで手動分割されたＳＴＮ及びＲＮは、７Ｔでの軸方向Ｔ２Ｗと１．５Ｔ-ＭＲテストセットでのＴ１Ｗ画像との間での既存の変換行列を使用して、１．５ＴでのＴ１Ｗテスト画像に位置合わせされる。特に、１．５Ｔ-ＭＲテストデータの７Ｔの対応する対が利用可能ではない場合、一対の７Ｔ及び１．５Ｔ-ＭＲ参照データが利用される。図２（ｂ）に示されるように、７Ｔ-ＭＲトレーニングデータのＳＴＮ及びＲＮは、ＡＮＮＴＳ又は任意の他のパッケージ（市販又は組織内）を使用して、変換行列（Ｍ＿（ｉ，７Ｔ））を用いて７Ｔ-ＭＲ参照データに位置合わせされ、次に、２つのデータ間の変換行列（Ｍ＿（７Ｔ→１．５Ｔ））によって対応する１．５Ｔ-ＭＲ参照データに位置合わせされる。したがって、変換行列（Ｍ＿１．５Ｔ）を使用して１．５Ｔ-ＭＲ参照データに位置合わせされた、１．５Ｔ-ＭＲテストデータからのＳＴＮ及びＲＮ（７Ｔ-ＭＲ画像の対応する対を有さない）は、１．５Ｔ-ＭＲ参照データに位置合わせされた７Ｔ-ＭＲトレーニングセットのＳＴＮ及びＲＮと同じ座標空間に正規化され位置合わせされる。
【００８３】
［統計学的形状モデル］
統計学的形状モデルは、トレーニング形状セットから解剖学的構造の形状のばらつきについての事前情報を利用できるようにする。一般的な統計学的形状モデルは点分布モデル（ＰＤＭ：point distributioon model）に基づき、ＰＤＭは、表面に沿って分布する１組のランドマーク点（すなわち、メッシュ内の頂点）によって形状を表し、形状のばらつきをモデリングする。特に、トレーニング形状のランドマーク点の中の対応性が、各形状のばらつきを捕捉し、形状予測の回帰モデルを構築するために必要とされる。最小記述長（ＭＤＬ：minimum description length）ベースの手法が、有効な形状対応性方法として認識されている。トレーニング形状セットにわたるそのようなランドマークを対応して得るために、ＭＤＬベースの方法を採用し、この方法はまず、トレーニング形状への球形メッシュパラメータ化の等角写像によってランドマークを生成し、次に、勾配降下アルゴリズムを適用することによって対応性を最適化する。５つの７Ｔ-ＭＲトレーニングセットにわたるＳＮ、ＳＴＮ、ＲＮ、ＧＰｉ、及びＴｈａ（これらの重要な構造を本発明の例として使用することを想起し、本技法は他の領域にも同様に適用される）の３Ｄ形状及び表面に沿って分布したランドマーク点は対応して、図１０に示されている。
【００８４】
本発明の予測枠組みでは、Ｎ個のトレーニングセットでのｎ_“ｘ”ランドマーク点の集合を有するＳＴＮ、およびｎ_“ｙ”ランドマーク点の集合を有するＳＴＮの予測子の三次元形状は、
【００８５】
【数１】
2020028721.tif 000003
として表される。
【００８６】
【数２】
2020028721.tif 000004
また、本明細書ではＳＴＮ（ここでは、予測する形状の例として使用される）の予測子の形状及び実際のデータでの予測ＳＴＮを、それぞれ以下のベクトル、
【００８７】
【数３】
2020028721.tif 000005
として示す。
【００８８】
【数４】
2020028721.tif 000006
トレーニングデータセットを実際のデータに位置合わせすると、トレーニングセット及び実際のデータ上のＳＴＮの予測子は位置合わせされ、それにより、それらの間の姿勢のばらつきを低減することができる（図２（ａ）及び（ｂ））。さらに、トレーニングセット及び実際のデータ上の予測子の表面ランドマーク点を大域的に位置合わせして、ＳＴＮをより確実に予測することができる。
【００８９】
【数５】
2020028721.tif 000007
【００９０】
【数６】
2020028721.tif 000008
トレーニングセットにおけるＳＴＮ及び予測子の形状パラメータ及び姿勢、並びに実際のデータでの予測子は、位置合わせされた表面ランドマーク点から得られ、以下に詳述するように、より低次元の空間での空間関係を抽出する。
【００９１】
【数７】
2020028721.tif 000009
【００９２】
【数８】
2020028721.tif 000010
ここでは、例として、ＰＣＡの非線形拡張であるカーネルＰＣＡ（ＫＰＣＡ）を利用し、トレーニング形状間の非線形関係を考慮することによって予測潜在性を高める。
【００９３】
【数９】
2020028721.tif 000011
特徴空間Ｆ内の線形ＰＣＡは、固有値問題Ｃｖｘ＝λｖｘを解くことを含み、式中、ｖｘ∈Ｆ＼｛０｝は固有ベクトルであり、λは共分散行列Ｃの非ゼロ固有値である。共分散行列Ｃを固有値問題に代入すると、以下が得られる。
【００９４】
【数１０】
2020028721.tif 000012
【００９５】
【数１１】
2020028721.tif 000013
次に、等価固有値問題が、
【００９６】
【数１２】
2020028721.tif 000014
として考慮される。
【００９７】
【数１３】
2020028721.tif 000015
【００９８】
【数１４】
2020028721.tif 000016
【００９９】
【数１５】
2020028721.tif 000017
【０１００】
【数１６】
2020028721.tif 000018
【０１０１】
【数１７】
2020028721.tif 000019
【０１０２】
【数１８】
2020028721.tif 000020
特徴空間Ｆの次元が無限である場合、共分散行列Ｃでの固有値問題は解くことができない。カーネル法により、φ写像を明示的に計算せずに、特徴空間内の主成分を得ることができる（「カーネルトリック」）。
【０１０３】
【数１９】
2020028721.tif 000021
【０１０４】
【数２０】
2020028721.tif 000022
【０１０５】
【数２１】
2020028721.tif 000023
ここで、ＳＴＮ形状パラメータβy（p）が、本発明の予測の枠組みにおいてトレーニングセットにおける｛βix｝と｛βiy｝との関係を使用して実際のデータでのβx（p）から予測されると仮定し、かつ、ＳＴＮの姿勢パラメータγy（p）も、トレーニングセットにおける｛γix｝と｛γiy｝との関係を使用して実際のデータでのγx（p）から予測されると仮定する。
【０１０６】
【数２２】
2020028721.tif 000024
【０１０７】
【数２３】
2020028721.tif 000025
【０１０８】
【数２４】
2020028721.tif 000026
【０１０９】
【数２５】
2020028721.tif 000027
次に、開示される方法は、実際の臨床データでの予測子の形状及び姿勢からＳＴＮ（又は他の関心のある脳構造）の形状及び姿勢の予測問題に注目し、これは、次の項に記載される回帰モデルを使用して、トレーニングセットにわたる形状パラメータと姿勢との関係を調べることによって行われる。
【０１１０】
［部分最小二乗回帰ベースの形状予測］
本開示の方法は、ＳＴＮ（又は予測する形状／領域）の形状パラメータ及び姿勢と、トレーニングセットでの予測子の形状パラメータ及び姿勢との依存性を見つける回帰問題に更に対処する。形状予測は劣決定（すなわち、低サンプルサイズ及び高次元問題）であるため、バイアス線形回帰が好ましいが、他の予測ツールも同様に使用可能である。
【０１１１】
本開示の方法は、可能な一実施形態において、トレーニングデータ、特に劣決定系でのベクトル集合間に高い相関を有する多重共線性問題であっても広く適用されているＰＬＳＲ技法の利用を提案する。ＰＬＳＲは、多重線形回帰を拡張し、予測子の集合Ｘ内の主因子（すなわち、潜在ベクトル）を見つけて、共分散が最大であるような応答ベクトルの集合Ｙを予測する。
【０１１２】
【数２６】
2020028721.tif 000028
予測子が複数の構造からなる場合、各予測子の姿勢及び形状パラメータを抽出し、複数の予測子の姿勢及び形状パラメータベクトルをそれぞれ連結し、予測の枠組みにおいて各予測子の姿勢及び形状についての局所情報を利用する。
【０１１３】
【数２７】
2020028721.tif 000029
ＰＬＳＲは、潜在行列Ｔ及びＦを見つけることによってＸ及びＹを分解する。ここで、Ｆは、Ｘ＝ＴＰT、Ｙ＝ＦＧTとすると、Ｆ＝ＴＢ＋ε…（１６）であり、式中、ＰはＸの負荷行列であり、ＧはＹの重み行列であり、Ｂは回帰重みの対角行列であり、εは回帰誤差である。潜在行列Ｔ及びＦの潜在ベクトルｔ及びｆはそれぞれＸ及びＹの列の線形結合が最大共分散を有するような２つの重みベクトルｗ及びｇを見つけることによって次式、
ｔ＝Ｘｗ，ｆ＝Ｙｇ、但し、ｗTｗ＝１，ｔTｔ＝１…（１７）
によって得られ、Ｘ及びＹの列の線形結合が最大共分散を有することはｔTｆが最大であることを意味する。
【０１１４】
特に、非線形反復部分最小二乗法によって、ｗ及びｇは反復更新され、次に、得られた潜在ベクトルｔ及びｆはＸ及びＹからそれぞれ減算される。反復は、Ｘがヌル行列になるまで続けられる。
【０１１５】
【数２８】
2020028721.tif 000030
【０１１６】
【数２９】
2020028721.tif 000031
［ＳＴＮの予測子としての近傍構造］
形状及び姿勢を予測する枠組みが、開示される方法について導出されている。予測方法によって使用される実際の形状について以下に更に開示する。
【０１１７】
回帰ベースのモデルでの予測正確性は、ＳＴＮと予測子構造との相関の程度による影響を受ける。近年の研究では、隣接する脳構造対が高度に相関することが報告されている。皮質下構造内のＳＴＮの予測子は、ＳＴＮとの空間隣接性と、特定の被験者からの実際のＭＲデータ（１．５Ｔ／３Ｔ又は７Ｔ）での可視性（容易に検出可能）とを考慮することによって導入される。
【０１１８】
ＳＴＮは、後外側領域における運動領域（例えば、パーキンソン病でのＤＢＳ手術の標的部位として）と、前内側領域の辺縁領域とに分割される小さく概して凸形の構造であり、幾つかの皮質下領域−不確帯、内側毛帯、エディンガー・ウェストファル核、内包、ＳＮ、ＲＮ、ＧＰｉ、及びＴｈａによって囲まれている。
【０１１９】
高磁場（７Ｔ）ＭＲデータでの優れたコントラスト及び高分解能は、大脳基底核及び視床、特にＳＴＮ及びその隣接構造内の皮質下構造を直接視覚化し識別することを可能にしたが、ＳＴＮの自動位置特定はそれでもなお問題である。さらに、皮質下構造の体積モデルは、７Ｔ磁石で取得されるＳＷＩ及びＴ２Ｗ画像のマルチコントラストを利用し、それにより、構造の空間情報−位置、次元、及び向きを提供することによって構築されている。
【０１２０】
図９は、軸方向又は冠状Ｔ２Ｗ-７Ｔ-ＭＲＩ（本開示に使用されたこの緩急での被験者の一人からの）での手動分割されたＳＴＮ及び隣接構造、即ち、ＳＮ、Ｔｈａ、ＧＰｉ、及びＲＮの３Ｄ構造を示す。ＳＮ、Ｔｈａ、ＧＰｉ、及びＲＮが、ＳＴＮと解剖学的近接性を有し、高磁場（７Ｔ）ＭＲデータ上のＳＴＮの他の近傍構造よりも線引きがより実現可能であることに留意されたい。他の関心のある構造についても、同様の関係が利用可能である。したがって、本発明では、ＳＴＮの潜在的な予測子としてＳＮ、Ｔｈａ、ＧＰｉ、及びＲＮを考え、空間関係を誘導するためのトレーニングセットとして、７Ｔ-ＭＲデータで得られたＳＴＮ及び予測子の３Ｄモデル（この場合では手動分割であるが、自動又は半自動も同様に使用可能である）を利用する。実際の７Ｔ-ＭＲデータ（別の被験者から）での分割されたＳＮ、Ｔｈａ、ＧＰｉ、ＲＮ、又はそれら全てを所与として、前項で紹介した技法に従い、７Ｔ-ＭＲ上の予測子からＳＴＮを自動的に予測することができる。それにもかかわらず、７Ｔ-ＭＲデータは、臨床での使用に広く利用可能ではないことがあるため、従来のＭＲデータ（１．５Ｔ）でのＳＴＮ予測を更に考慮する必要がある。
【０１２１】
図１１は、１．５Ｔ及び７Ｔ-ＭＲ撮像をそれぞれ用いて得られた軸方向Ｔ２Ｗ画像でのＳＮ、ＳＴＮ、Ｔｈａ、ＧＰ、及びＲＮ領域を示す。軸方向Ｔ２Ｗ-７Ｔ-ＭＲデータは、比較のためにＴ２Ｗ-１．５Ｔ-ＭＲデータに位置合わせされた。
【０１２２】
１．５Ｔ（臨床）ＭＲ撮像でのＳＮ、ＧＰ領域内のＧＰｉ、及びＴｈａから、それらの領域内の強度が低くＳＮＲが低いことに起因して、ＳＴＮを識別することが難しく、一方、７Ｔ-ＭＲ撮像はそれらの領域に明確な視覚化を提供することに気付く。また、ＲＮ領域は、１．５Ｔ-ＭＲ撮像よりも７Ｔ-ＭＲ撮像でより明確に示されている。しかし、ＲＮが１．５Ｔ-ＭＲ撮像において他の構造よりも識別が容易であることに留意する。さらに、幾つかの研究は、ＤＢＳ中にＳＴＮ内に電極を移植するためのガイダンスとしてＲＮ（例えば、ＲＮへのＳＴＮの相対位置）を利用している。その結果、１．５Ｔ-ＭＲデータでは、ＲＮ構造が程々に可視であり、且つＳＴＮの近傍にあるため、ＳＴＮの自動予測子としてＲＮ構造を選ぶ。次に、実際の１．５Ｔ-ＭＲデータでのＳＴＮが、７Ｔ-ＭＲトレーニングデータで誘導された空間関係を利用することにより、同じデータで得られたＲＮから自動的に予測される。
【０１２３】
【数３０】
2020028721.tif 000032
［実験検証］及び［予測測定］
本開示で紹介される特定の実施形態の幾つかの価値を強調するために、これより、例示の形態で、上で説明され例示した実装形態を使用して得た実験結果の幾つかを説明する。
【０１２４】
予測正確性が、ダイス係数（ＤＣ：Dice Coefficient）形状類似性、表面ランドマーク点間の平均距離、並びに予測ＳＴＮとグラウンドトルース（すなわち、手動分割ＳＴＮ）との間でのセントロイド、寸法、及び向きの誤差によって測定される。ＤＣ値は、２（ＶＡ∩ＶＢ）／（ＶＡ＋ＶＢ）による形状の類似性測定として計算され、ここで、ＶＡ及びＶＢは比較されている構造Ａ及びＢの各体積である。
【０１２５】
【数３１】
2020028721.tif 000033
さらに、予測ＳＴＮ（すなわち、ＳＴＮの予測表面ランドマーク点から再構築された３Ｄ形状）及び手動分割それぞれの慣性楕円体の冠状、矢状、及び軸方向でのセントロイドｇ＝（ｇcor，ｇsag，ｇaxi）、３つの半短径（半径）の長さｌ＝（ｌx，ｌy，ｌz）、及び方位角ｏ＝（φ，θ，ψ）を計算して、幾何学的形状情報の類似性を分析する。より具体的には、ｇは、ランドマーク点の各座標の平均として計算される。ｌは、ｎy中心化ランドマーク点に対してＰＣＡを実行することによって得られる最初の３つの最大固有値から計算される。ｏは、対応する固有ベクトルからなる回転行列からオイラー角を計算して得られる。
【０１２６】
【数３２】
2020028721.tif 000034
［予測子の評価］
本開示の方法のＳＴＮの予測が、７Ｔ-ＭＲからの各構造−ＳＮ、ＲＮ、ＧＰｉ、又はＴｈａ−を使用して実行され、そのような予測子の予測性能を評価した。各データセットでの構造のトレーニングセットは、１個抜き（leave-one-out）法を使用して構築される。例えば、データセット１のトレーニングセットは、データセット１からの構造を除いたデータセット２、３、４、及び５からの各構造のランドマーク点ベクトルからなる。ここで５つのセットが例示のために使用されるが、本開示の構成要素は任意の数のトレーニングデータセットに適用される。
【０１２７】
【数３３】
2020028721.tif 000035
【０１２８】
【数３４】
2020028721.tif 000036
【０１２９】
【数３５】
2020028721.tif 000037
【０１３０】
【数３６】
2020028721.tif 000038
【０１３１】
【数３７】
2020028721.tif 000039
図１２では、５つのデータセットの１つでの予測子としてＳＮ及びＲＮ構造をそれぞれ使用することによって予測されたＳＴＮの対応するセントロイドを有する３Ｄ形状及び慣性楕円体が示され、対応する手動分割のものが重ねられている。データセットでの予測ＳＴＮのセントロイド、寸法、及び向きがなお手動分割と同等であり、形状がその他のデータセットよりもわずかに正確であることが認められる。
【０１３２】
［１．５Ｔ-ＭＲデータセットでの予測結果］
高磁場（７Ｔ）ＭＲ撮像は、優れたコントラスト及びより高い分解能に起因して、ＳＴＮの視覚化に成功し、隣接構造ＳＮを分離した。しかし、７Ｔ-ＭＲデータは、標準の臨床用途に常に利用可能であるとは限らない。したがって、従来の臨床１．５Ｔ（又は３Ｔ）データでのＳＴＮ予測は、ＤＢＳ標的及び手術後プログラミングについて極めて重要な情報を提供することができる。４つの（説明のための例として）１．５Ｔ-ＭＲデータセットでのＳＴＮの予測は、臨床１．５Ｔ-ＭＲＩで可視のＲＮ構造のみを使用して実行した。この比較のために、対応する（同じ被験者の）７Ｔ-ＭＲデータセットでのＲＮを使用して得られた予測結果と、１．５Ｔ及び７Ｔ-ＭＲテストセットの両方に位置合わせされたトレーニングデータセットにわたるＳＴＮの形状の平均とを提示する。
【０１３３】
【数３８】
2020028721.tif 000040
さらに、測定パラメータは、１．５Ｔ及び７ＴのＭＲデータセットにわたるＲＮへのＳＴＮの相対位置を分析するために、予測ＳＴＮ、ＳＴＮの平均形状、及び手動分割ＳＴＮのセントロイドのそれぞれとＲＮとの間の距離を含む。なお、本発明での予測枠組みにおいて、ＲＮとＳＴＮとの空間関係を利用した。複数の脳構造間のそのような相対姿勢は、神経変性疾患の診断についての有用情報を提供し得る。さらに、ＲＮは、近年の研究において、ＳＴＮのＤＢＳ標的の参照として利用されている。
【０１３４】
予測ＳＴＮが、１．５Ｔ及び７ＴのＭＲデータセットにわたり両側でＲＮのセントロイドの５ｍｍ〜８ｍｍ前方、３．９ｍｍ〜６．７ｍｍ側方、且つ０．７ｍｍ〜３．５ｍｍ上側（特に、手動分割ＳＴＮでは、４ｍｍ〜７．４ｍｍ前方、５．２ｍｍ〜７．５ｍｍ側方、及び０．２ｍｍ〜４．２ｍｍ上側）に位置したことが観測された。１．５Ｔ-ＭＲデータセットでのＲＮへの予測ＳＴＮ、ＳＴＮの平均形状、及び手動分割ＳＴＮの相対位置は、７Ｔ-ＭＲデータセットでのものと同様であった。また、１．５Ｔ及び７ＴのＭＲデータセットの両方で、ＲＮへの予測ＳＴＮの全体相対位置が、ＲＮへのＳＴＮの平均形状のものよりも、ＲＮへの手動分割ＳＴＮのものにより近かった。
【０１３５】
【数３９】
2020028721.tif 000041
ＳＴＮの予測結果、平均形状、及び手動分割の実際の寸法（ｍｍ３）での体積も測定した。１．５Ｔ及び７ＴのＭＲデータセットでの手動分割の体積は、２つのデータセット間の異なる体積はＤＣ類似値の不正確な比較に繋がるおそれがあるため、同様であるはずである。１．５Ｔでの手動分割ＳＴＮの平均体積は、７Ｔ-ＭＲデータセットでのものよりも１７ｍｍ３（左）大きく、１１ｍｍ３（右）大きい（１．５Ｔ-ＭＲデータへの７Ｔ-ＭＲデータの位置合わせ後の再サンプリングアーチファクトの閾値処理により、２つのデータセット間の体積差を生じさせることがある）。予測ＳＴＮ及び平均形状の平均体積は、１．５Ｔ及び７ＴのＭＲデータセットの両方で、手動分割ＳＴＮのものよりも大きい。また、７Ｔ-ＭＲデータセットでの予測ＳＴＮの左側の平均体積が、１．５Ｔ-ＭＲデータセットのものよりも１４ｍｍ３大きいことも認められた。これは、類似性測定にも拘わらず、１．５Ｔ-ＭＲデータセットでの予測ＳＴＮの左側のＤＣ値が、７Ｔ-ＭＲデータセットでのものよりも７％高かったかの理由を説明する。換言すれば、７Ｔ-ＭＲデータセットでの左ＳＴＮの手動分割と予測との体積差は、１．５Ｔ-ＭＲデータセットでのものよりも大きかった（すなわち、７Ｔでの６０ｍｍ３及び１．５Ｔ-ＭＲデータセットでの３０ｍｍ３）。両ＭＲデータセットでの予測ＳＴＮの体積は、ＳＴＮの平均形状のものと同様であった。データセットが、サイズ及び多様性（例えば、コントラスト）の両方の点で大きいほど、強力な予測が予期される。
【０１３６】
１．５Ｔ-ＭＲ及び７Ｔ-ＭＲのデータセットからの全体予測結果は、臨床で実現可能な測定値に完全に基づいて、トレーニングデータセットにわたりＳＴＮの平均形状よりも手動分割ＳＴＮにはるかに近く、開示される枠組み及び本発明の予測力を証明する。これは、臨床ＭＲＩから、開示される発明を使用して、神経外科のために重要構造をいかに位置特定することができるかを明らかに示す。
【０１３７】
図１３及び図１４は、特定の被験者からの１．５Ｔ-ＭＲデータセット及び７Ｔ-ＭＲデータセットのそれぞれでの３つの直交する平面でのＲＮ構造を使用することによって予測されたＳＴＮ、手動分割ＳＴＮ、及びトレーニングデータセットにわたるＳＴＮの平均の輪郭、３Ｄ形状、及びセントロイドを有する対応する慣性楕円体を示す。また、冠状、矢状、及び軸方向での手動分割ＳＴＮのものと比較した、予測ＳＴＮ及びＳＴＮの平均形状のセントロイド位置も提供する。次に、ＲＮへのＳＴＮのセントロイドの相対位置を測定する。
【０１３８】
【数４０】
2020028721.tif 000042
さらに、１．５Ｔ-ＭＲデータセットでの予測結果は、１．５Ｔ-ＭＲデータセットに位置合わせされたトレーニングセットにわたるＳＴＮの平均形状よりもはるかに良好であり、一方、７Ｔ-ＭＲデータセットでの予測結果は、この７Ｔ-ＭＲデータセットに位置合わせされたトレーニングセットにわたるＳＴＮの平均形状と同様であった。特に、１．５Ｔ-ＭＲデータセットでの右ＳＴＮの予測結果が、７Ｔ-ＭＲデータセットでのものよりも手動分割に近かったが、１．５Ｔ-ＭＲデータセットでのトレーニングセットにわたる平均形状がかなり不正確であった（図１３参照）ことに留意する。これは、この１．５Ｔ-ＭＲデータでの右ＳＴＮの予測が、対応する７Ｔ-ＭＲデータセットでのＳＴＮよりも効率的であったことを示す。
【０１３９】
本開示により、ＳＴＮ領域内コントラストが低く（例えば、１．５Ｔ-ＭＲ撮像）又はオクルージョン（例えば、ノイズ）を有するデータであっても、ＳＴＮの位置、寸法、及び向きを妥当に推定可能なことに留意されたい。これは、高磁場（７Ｔ）ＭＲトレーニングデータから学習されたＳＴＮと予測子との空間依存性を利用して行われる。特に、予測ＳＴＮの慣性楕円体のセントロイド、寸法、及び方位角は、領域内のより精密なＤＢＳ標的についての非常に重要な情報を提供し得る。
【０１４０】
７Ｔ-ＭＲトレーニングデータセットでのＳＴＮ予測テストから、開示される本発明は、ＳＮ及びＲＮ構造がＧＰｉ及びＴｈａよりもＳＴＮの良好な予測子であることを実証し、より高いＤＣ類似値並びにセントロイド、寸法、及び向きにおいてより低い平均誤差を有し、それらの構造がＳＴＮの形状及び姿勢についてより有意な情報を提供することを示す。さらに、ＳＮ、ＲＮ、ＧＰｉ、及びＴｈａの全てが７Ｔ-ＭＲテストデータで利用可能である場合、複数の構造は、同等の測定を有するＳＴＮの妥当な予測子として見なすことができる。
【０１４１】
【数４１】
2020028721.tif 000043
なお、Ｔｈａ構造は比較的より大きな体積を有し、したがって、ＳＴＮとＴｈａとのセントロイドの距離は他の予測子よりもはるかに大きい。また、Ｔｈａの手動分割は他の予測子の手動分割ほど精密ではないことがある。これらの要因は、予測ＳＴＮの不正確な姿勢に繋がり得る。本発明をこれら及び他の脳構造に適用する場合、このすべてを考慮に入れるべきである。
【０１４２】
さらに、方法は、ＲＮ構造が１．５Ｔ-ＭＲデータセットでＳＴＮの良好な予測子としての役割を果たすことも実証する。１．５Ｔ-ＭＲ撮像でのＲＮからのＳＴＮの全体予測結果は、７Ｔ-ＭＲ撮像でのＲＮからの得られた全体予測結果と同様であった。予測ＳＴＮの完全な形状は手動分割ＳＴＮほど正確ではないが、予測ＳＴＮの姿勢情報、即ち、セントロイド及び寸法は、手動分割ＳＴＮと同等であった。
【０１４３】
１．５Ｔ-ＭＲ撮像での予測ＳＴＮの平均ＤＣ値は、７Ｔ-ＭＲ撮像よりも良好であるが、１．５Ｔ-ＭＲデータセットでの方位角の距離誤差はより高く、その理由は、１．５Ｔ-ＭＲ撮像での予測ＳＴＮの平均体積がＳＴＮの手動分割により近いためである。また、本開示の方法によれば、１．５Ｔ-ＭＲ撮像で予測されたＳＴＮの表面ランドマーク点の平均数は、７Ｔ-ＭＲ撮像で予測されたＳＴＮのものよりもはるかに少ない（例えば、１．５Ｔ-ＭＲデータセットでの平均１４６点及び７Ｔ-ＭＲデータセットでの平均４３６点）。これは、１．５Ｔ-ＭＲデータセットでのＳＴＮのわずかにより正確な形状予測をもたらし得るが、より大きな方位角誤差をもたらし得る。
【０１４４】
しかし、７Ｔ-ＭＲデータセットで分割されたＲＮ構造が１．５Ｔ-ＭＲデータセットに位置合わせされ、次に、この例ではＳＴＮの予測子として利用されたことに留意されたい（７Ｔ-ＭＲＩが臨床用途に利用可能ではない場合、ＲＮは１．５Ｔ-ＭＲＩで完全に分割し、１．５Ｔ-ＭＲ参照データに位置合わせすることができる）。ここで使用されたＲＮは、１．５Ｔ-ＭＲデータよりも７Ｔ-ＭＲデータでより明確に識別することができるため、１．５Ｔ-ＭＲデータでのＳＴＮのより信頼できる予測子として見なされる。
【０１４５】
１．５Ｔ及び７ＴのＭＲデータセットの両方でのＲＮへの予測ＳＴＮの相対位置は、ＳＴＮとＲＮとの平均形状間よりも手動分割ＳＴＮとＲＮとの間により近かった。なお、本発明での予測の枠組みではＳＴＮとＲＮとの空間的関係（特に、神経変性疾患の診断に有用であり得る）を学習したが、ＲＮの相対姿勢はＳＴＮの平均形状において考慮されなかった。
【０１４６】
さらに、１．５Ｔ及び７ＴのＭＲデータセットの両方でのＳＴＮの予測は、トレーニングデータセットにわたるＳＴＮの平均形状がＳＴＮの手動分割により近い場合、より正確な結果を示した。特に、方位角の距離誤差は主に、トレーニングデータセットにわたる平均形状の方位角の距離誤差による影響を受ける。これらの発見は、テストデータセットでの予測子へのトレーニングデータセットでの予測及びＳＴＮの大域的な位置合わせにより（それにより、トレーニングセットにわたるＳＴＮの平均形状の姿勢が、手動分割により近づくことができる）、トレーニングセットでの予測子とテストセットでの予測子との間での姿勢のばらつきをそれぞれ低減して、ＳＴＮの予測が更に改善可能なことを暗に示す。
【０１４７】
方法は、高磁場（７Ｔ）ＭＲトレーニングデータセットから得られるＳＴＮとＲＮとの空間依存性を利用して、１．５Ｔ-ＭＲ撮像でのＲＮからＳＴＮを予測した。１．５Ｔ（臨床）Ｔ１Ｗ画像及び高磁場７Ｔ-Ｔ２Ｗ-ＭＲ撮像は、異なる分解能を有し（例えば、それぞれ０．９８×０．９８×１ｍｍ３及び０．４×０．４×２ｍｍ３）、これは、位置合わせがないことを意味し、予測での不正確な姿勢に繋がるおそれがある。したがって、本発明では、各７Ｔ-ＭＲトレーニングデータセットを７Ｔ-ＭＲテストデータに位置合わせし、次に、７Ｔ-ＭＲテストデータを１．５Ｔ-ＭＲテストデータに位置合わせすることによって得られる変換行列を使用して、各７Ｔ-ＭＲトレーニングデータでのＳＴＮ及びＲＮ構造を１．５Ｔ-ＭＲテストデータの共通座標空間に正規化した。さらに、１．５Ｔ-ＭＲ及び７Ｔ-ＭＲの参照データは、１．５Ｔ-ＭＲテストデータに対応する７Ｔ-ＭＲテストデータが利用可能ではない場合、図２（ｂ）でのように７Ｔ-ＭＲトレーニングデータセット及び１．５Ｔ-ＭＲテストセットの共通座標空間として見なされる。
【０１４８】
方法は、本開示の手法が、可視であることがあるが、それでもなお自動的に分割することが難しいか、又は視野の制限に起因して全く見えないこともあるＳＴＮの形状及び姿勢を直接予測可能なことを実証する。さらに、ＳＴＮを識別することが可能ではない臨床１．５Ｔ-ＭＲデータであっても、程々に可視であり、ＳＴＮに隣接するＲＮ構造を使用して予測可能なことが示された。
【０１４９】
本開示の方法によるＳＴＮの予測は、５つのトレーニングセットを使用して例示された。ここで報告された結果は、高磁場（７Ｔ）画像のより大きなデータセットを使用して更に改善することができる。さらに、データベース及びサンプルサイズが増大するにつれて、特定の事例のサブグループ（例えば、年齢、性別、又は疾患及び／又は症状固有）が、予測モデルをグループに適応させる場合、予測子の正確性を改善することを証明し得る。本開示の一実施形態は、そのような情報に基づくデータのグループ化を含む。最後に、トレーニング形状セットにわたる第３脳室幅及び両被殻距離等のＳＴＮのばらつきに関連する解剖学的情報は、追加の予測子として見なすことができ、提案される本開示に組み込むことができる。
【０１５０】
本開示のシステム及び方法の実施形態は、神経外科医が、はるかに洗練され正確な手術及びＤＢＳ移植処置に従事できるようにする。特に、血管及び重要な神経構造を回避する能力並びにＤＢＳ移植のより正確な初期標的及び経路を選択する能力は、ＤＢＳ移植に必要な時間及び合併症率を低減する潜在性を有する。開示されるシステム及び方法は、ロボット手術等の自動処置に従事する能力を強化する。そのような状況では、患者固有アトラスは、ロボットシステムに標的及び全侵入経路を指示し得る。遠隔手術もより高品質でより柔軟な撮像技法の使用から恩恵を受ける。白質等の撮像が難しい組織の描写を改善する能力により、神経外科医は、従来は届かなかった組織を標的とすることができる。
【０１５１】
システム及び方法が、その趣旨又は中心の特徴から逸脱せずに他の特定の形態で実施し得ることが理解される。したがって、本例及び実施形態は、あらゆる点で限定ではなく例示として見なされるべきであり、システム及び方法は、本明細書の詳細な説明に限定されるべきではない。
【図１】
2020028721.tif 000044
【図２】
2020028721.tif 000045
【図３】
2020028721.tif 000046
【図４】
2020028721.tif 000047
【図５】
2020028721.tif 000048
【図６】
2020028721.tif 000049
【図７】
2020028721.tif 000050
【図８】
2020028721.tif 000051
【図９】
2020028721.tif 000052
【図１０】
2020028721.tif 000053
【図１１】
2020028721.tif 000054
【図１２】
2020028721.tif 000055
【図１３】
2020028721.tif 000056
【図１４】
2020028721.tif 000057
【外国語明細書】
2020028721000001.pdf