(57)【要約】脳内電流シミュレーション方法は，患者の断層画像データのうち少なくとも脳の一部を含む頭部画像データを提供する第１の工程；第１の工程で提供された頭部画像データのうち脳を構成している少なくとも１つの領域を微小要素に分割したそれぞれの微小多面体単位からなる３次元脳モデルを形成する第２の工程；患者の頭部上にコイルを配置し，コイルに電流を印加することによって患者の脳に磁気刺激を与え，磁気刺激に対する患者の反応が観察されたときの条件のうち，少なくともコイルの位置及び向きの条件と，コイルに印加された電流の条件と，コイルの生成磁界に関わる構造の条件と，を含む第１の情報を提供する第３の工程；３次元脳モデルのそれぞれの微小多面体単位内に誘導される渦電流又は電界を，第３の工程で提供された第１の情報，及び微小多面体単位ごとに割り付けられた導電率を含む第２の情報に基づいて計算する第４の工程を含む。
2015122369.tif 000002
【特許請求の範囲】
【請求項１】
患者の断層画像データのうち少なくとも脳の一部を含む頭部画像データを提供する第１の工程と，
前記第１の工程で提供された前記頭部画像データのうち前記脳を構成している少なくとも１つの領域を微小要素に分割したそれぞれの微小多面体単位からなる，３次元脳モデルを形成する第２の工程と，
前記患者の頭部上にコイルを配置し，前記コイルに電流を印加することによって前記患者の脳に磁気刺激を与え，前記磁気刺激に対する患者の反応が観察されたときの条件のうち，少なくとも前記コイルの位置及び向きの条件と，前記コイルに印加された電流の条件と，前記コイルの生成磁界に関わる構造の条件と，を含む第１の情報を提供する第３の工程と，
前記３次元脳モデルの前記それぞれの微小多面体単位内に誘導される渦電流又は電界を，前記第３の工程で提供された前記第１の情報，及び前記微小多面体単位ごとに割り付けられた導電率を含む第２の情報に基づいて計算する第４の工程を含む，脳内電流シミュレーション方法。
【請求項２】
前記第４の工程で計算された渦電流又は電界の分布を視覚的に表示する第５の工程を含む，請求項１の脳内電流シミュレーション方法。
【請求項３】
前記それぞれの微小多面体単位が，灰白質，白質，及び脳脊髄液のいずれかの導電率を有するように割り付けられた，請求項１又は２の脳内電流シミュレーション方法。
【請求項４】
前記第１の情報が，前記電流の電流値または電圧値の少なくともいずれかを含む，請求項１〜３のいずれかの脳内電流シミュレーション方法。
【請求項５】
前記第４の工程において，前記渦電流又は電界の計算がスカラーポテンシャル有限差分法で行われることを特徴とする，請求項１〜4のいずれかの脳内電流シミュレーション方法。
【請求項６】
前記コイルに印加された電流が交流又は脈流である請求項１〜５のいずれかの脳内電流シミュレーション方法。
【請求項７】
患者の断層画像データのうち少なくとも脳の一部を含む頭部画像データを提供する第１の手段と，
前記第１の手段で提供された前記頭部画像データのうち前記脳を構成している少なくとも１つの領域を微小要素に分割したそれぞれの微小多面体単位からなる，３次元脳モデルを形成する第２の手段と，
前記患者の頭部上にコイルを配置し，前記コイルに電流を印加することによって前記患者の脳に磁気刺激を与え，前記磁気刺激に対する患者の反応が観察されたときの条件のうち，少なくとも前記コイルの位置及び向きの条件と，前記コイルに印加された電流の条件と，前記コイルの生成磁界に関わる構造の条件と，を含む第１の情報を提供する第３の手段と，
前記３次元脳モデルの前記それぞれの微小多面体単位内に誘導される渦電流又は電界を，前記第３の手段で提供された前記第１の情報，及び前記微小多面体単位ごとに割り付けられた導電率を含む第２の情報に基づいて計算する第４の手段を含む，脳内電流シミュレーション装置。
【請求項８】
前記第４の手段で計算された渦電流の分布を視覚的に表示する第５の手段を含む，請求項７の脳内電流シミュレーション装置。
【請求項９】
前記それぞれの微小多面体単位が，灰白質，白質，及び脳脊髄液のいずれかの導電率を有するように割り付けられた，請求項７又は８の脳内電流シミュレーション装置。
【請求項１０】
前記第１の情報が，前記電流の電流値または電圧値の少なくともいずれかを含む，請求項７〜９のいずれかの脳内電流シミュレーション装置。
【請求項１１】
前記第４の手段において，前記渦電流又は電界の計算がスカラーポテンシャル有限差分法で行われることを特徴とする，請求項７〜１０のいずれかの脳内電流シミュレーション装置。
【請求項１２】
前記コイルに印加された電流が交流又は脈流である請求項７〜１１のいずれかの脳内電流シミュレーション装置。
【請求項１３】
請求項７〜１２のいずれかの脳内電流シミュレーション装置を含む経頭蓋磁気刺激システム。
【請求項１４】
患者の断層画像データのうち少なくとも脳の一部を含む頭部画像データを提供する第１の手段と，
前記第１の手段で提供された前記頭部画像データのうち前記脳を構成している少なくとも１つの領域を微小要素に分割したそれぞれの微小多面体単位からなる，３次元脳モデルを形成する第２の手段と，
前記患者の頭部上にコイルを配置して前記コイルに電流を印加することによって前記患者の脳に磁気刺激を与える条件のうち，少なくとも前記コイルの位置及び向きの条件と，前記コイルに印加された電流の条件と，前記コイルの生成磁界に関わる構造の条件と，を含む第１の情報を提供する第３の手段と，
前記３次元脳モデルの前記それぞれの微小多面体単位内に誘導される渦電流又は電界を，前記第３の手段で提供された前記第１の情報，及び前記微小多面体単位ごとに割り付けられた導電率を含む第２の情報に基づいて計算する第４の手段を含む，脳内電流シミュレーション装置。
【請求項１５】
前記第４の手段で計算された渦電流又は電界の分布を視覚的に表示する第５の手段を含む，請求項１４の脳内電流シミュレーション装置。
【請求項１６】
前記それぞれの微小多面体単位が，灰白質，白質，及び脳脊髄液のいずれかの導電率を有するように割り付けられた，請求項１４又は１５の脳内電流シミュレーション装置。
【請求項１７】
前記第１の情報が，前記電流の電流値または電圧値の少なくともいずれかを含む，請求項１４〜１６のいずれかの脳内電流シミュレーション装置。
【請求項１８】
前記第４の手段において，前記渦電流又は電界の計算がスカラーポテンシャル有限差分法で行われることを特徴とする，請求項１４〜１７のいずれかの脳内電流シミュレーション装置。
【請求項１９】
前記コイルに印加された電流が交流又は脈流である請求項１４〜１８のいずれかの脳内電流シミュレーション装置。
【請求項２０】
請求項１４〜１９のいずれかの脳内電流シミュレーション装置を含む経頭蓋磁気刺激システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は，脳内に誘導される電流又は電界をシミュレーションする方法と装置に関する。本発明はまた，その装置を組み入れた経頭蓋磁気刺激システムに関する。
【背景技術】
【０００２】
経頭蓋磁気刺激は，電磁誘導によって脳内に電流又は電界を生じさせ，ニューロンを刺激する手法である。この手法によれば，図１，２に示すように，頭部の皮膚上に置いた刺激コイルに電流（例えば，交流）を印加することで変動磁場を形成するとともにその変動磁場の影響を受けて脳内にコイル電流とは逆向きの渦電流又は電界を誘導し，その渦電流又は電界でニューロンを刺激することによって活動電位を発生させる。【０００３】
図３に刺激コイル駆動回路の一例を示す。瞬間的な電流をコイルに発生させる原理としては，まず電源（交流電源，電源回路，昇圧回路を含む。）からコンデンサへ電荷を蓄積する。その後，サイリスタをターンオンすることにより刺激コイルに電流を流す。刺激コイルとコンデンサの共振回路にダイオードを通して電流が流れたのち，サイリスタがオフになる。これにより，図４に示す正弦波１周期分の電流が刺激コイルに流れる。【０００４】
経頭蓋磁気刺激は神経伝導速度の計測を始めとした臨床検査や脳機能研究に用いられている。【０００５】
近年においては，神経障害性疼痛やパーキンソン病，うつ病などの治療的応用として磁気刺激が注目されている。このような病気においては薬剤による治療では効果が見られないケースがあり，そのような場合電極を脳に埋め込むことで脳に電気刺激を与えるという治療法がある。しかしその場合，開頭手術を必要とすることから希望しない患者が多い。そこで，手術を要さない非侵襲的な磁気刺激を反復して行う反復経頭蓋磁気刺激が治療法として研究されている。例えば難治性神経障害性疼痛において，大脳の一次運動野に磁気刺激を行った後１日間ほど除痛効果をえられたことが報告されている。【０００６】
しかし，従来の磁気刺激装置は約７０Ｋｇの重量があり，また設置のために電気工事が必要となるため，設備の整った医療機関でのみ利用可能となっている。また，実際の治療時には患者のＭＲＩデータを参照しながら刺激位置を決定するため，熟練した医療従事者による治療が必要である。難治性神経障害性疼痛の治療では，標的となる一次運動野の上に，１ｍｍの単位でコイルの位置決めを行う。しかし，除痛効果を継続して得るためには医療機関に毎日通う必要がある。そこで，図５に示す，非医療従事者による操作のみで利用可能な，在宅治療に用いる新たな磁気刺激装置の開発が行われている。
【先行技術文献】
【特許文献】
【０００７】
本発明者らは，図５に示す磁気刺激装置を開発し，８字型を改良した磁場発生コイルや位置決めに関しては既に特許出願を行っている(WO2010/147064，特許公開公報2012-125546)。
【発明の概要】
【発明が解決しようとする課題】
【０００８】
経頭蓋磁気刺激療法では，脳内のターゲット部位に誘起される渦電流によりニューロンを刺激する。
この渦電流の強度や分布は，コイルの構造や特性，印加電流等の磁場生成能と，患者に対するコイルの位置により決定されると考えられる。
しかしながら，磁気治療中の対象患者の脳内に実際に生成された渦電流を電極などを用いて実測することは，現実的ではない。
そこで，コイルにより生成される磁場分布から，脳内に誘起される渦電流をシミュレーションにより解析する試みがなされている。
これまでの脳内渦電流分布シミュレーションは，標準モデルを用いた計算であり，脳を一様な導電体と見なし電流分布を求めるものであった。
しかし，実際の脳は患者個人毎に大きさや，形状が異なると共に，組織により導電率が異なるため，患者個別の電流強度分布の計算は困難であると共に，シミュレーション計算自体も複雑で時間を要する点が問題であった。
【課題を解決するための手段】
【０００９】
そこで，本発明は，個々の患者の断層像から脳形状を抽出した脳モデル，及びコイルをモデル化し，任意の脳モデルとコイルの位置関係を決定した上で，渦電流解析を行う手法を用いることにより，経頭蓋磁気刺激治療を行おうとしている患者の，脳内で生成される渦電流の大きさや分布を実測値に極めて近い正確さでシミュレートする方法，装置を確立することを目的とする。
この手法では，電磁界解析にスカラーポテンシャル有限差分法（Scalar Potential Finite Difference）（以下，「SPFD法」という。）を利用することが好ましい。
この手法によれば，刺激コイルの装着位置に伴う電流強度分布をシミュレーションにより観察することが可能となり，刺激強度と刺激部位における渦電流強度や，刺激を加えたい部位に対するコイルの最適位置を，実際に刺激を加えることなく検証が可能となる。
また，SPFD法を使用することにより従来の有限要素法に比して大幅に時間短縮が可能となるため医師の研究，治療において効率的な対応が可能となる。
さらに，疾患に応じて最適な誘導電流分布のコイルを選定することも容易となり，治療上の意義は大きい。【００１０】
脳内の渦電流分布の解析には，市販の渦電流解析ソフトにおいて用いられている有限要素法ではなくSPFD法を用い，脳を含む解析対象の立体を微小要素の微小多面体単位（例えば，立方体）に分割し，各々の微小多面体単位に発生する渦電流を求めた。SPFD法を用いた解析用のソフトウェア（プログラム）は独自に開発した。
本発明の実施に当たっては有限要素法など，他のシミュレーション算法を用いてもよいことは言うまでもないが，SPFD法を用いるオリジナルソフトウェアの利点は，以下の３つである。【００１１】
第１の利点は，有限要素法を用いた市販のソフトウェアでは分割モデルの要素数が１００万ほどに制限されているが，オリジナルソフトウェアに要素数の制限がない，ということである。分割数，ボクセルサイズを任意に指定することができるので，モデルをより細かい要素に分割することができ，細かい解析を行うことが可能である。第２の利点は，空気層を作成する必要がないことである。市販のソフトウェアは，脳モデルのメッシュに合わせた空気層を脳モデルの周囲に作成する必要があったが，本発明は脳モデルのみを作るため，解析用モデルの作成に要する時間が大幅に短くなる。第３の利点は，計算手法を変えたことにより，計算時間が市販のソフトウェアの約1/20に短縮されたことがあげられる。【００１２】
以下，解析に用いたSPFD法の詳細を説明する。
ファラデーの法則により，脳内に誘導される電場Eと，外部磁束密度Bは，数式（１）を満たす。
【数１】
2015122369.tif 000003
磁気ベクトルポテンシャルAおよび電気スカラーポテンシャルφを，数式（２）および（３）を満たすように選ぶと，これらのポテンシャルから得られる電場および磁場は，式（１）を満たす。
【数２】
2015122369.tif 000004
【数３】
2015122369.tif 000005
【００１３】
ここで，数式（３）の総電場Eは，式（４）から（６）に示すように，ベクトルポテンシャルに由来する電場E1と，スカラーポテンシャルに由来する電場E2に，分解できる。
2015122369.tif 000006
数式（５）において，１行目から２行目への変形は，ビオ・サバールの法則によった。式中のμ0は真空の透磁率，rは空間の位置ベクトル，r'は磁場発生コイル上の位置ベクトル，I(r')は位置r'におけるコイル電流である。ベクトルポテンシャルに由来する電場E1は，コイルの巻線形状と電流のみによって決まり，数式（５）を使って，比較的容易に計算可能である。スカラーポテンシャルに由来する磁場E2は，被験者の頭部形状に依存し，これを求めるためには，SPFD法や有限要素法などの高度な数値解析手法が必要である。総電場Eも同様に，高度な数値解析手法を使って求める必要がある。【００１４】
数式（３）とオームの法則の数式（７）を電流連続の数式（８）に代入すると，数式（９）が得られる。
2015122369.tif 000007
数式（７）において，Jは誘導電流密度，σは生体の導電率を表す。【００１５】
数式（９）について８個の微小多面体単位が共有する接点での離散化を行うと，数式（９）式は数式（１０）になる。
2015122369.tif 000008
数式（１０）において，φnは接点nにおける電気スカラーポテンシャル，Anは接点0と接点nを結ぶボクセルの辺nに平行な外部磁気ベクトルポテンシャルの成分で辺の中心における値，lnは辺nの長さ，Snは辺nのコンダクタンスである。【００１６】
anを辺nに垂直な直方体の面の面積，σnを辺nで接する４つの直方体の平均導電率として数式（１０）に数式（１１）を適用して方程式を解くことで電気スカラーポテンシャルを求めた。求めた電気スカラーポテンシャルをもとに，数式（３），（７）を用いて渦電流密度を求めた。
2015122369.tif 000009
上に述べたように，本シミュレーション方法では，SPFD法を利用することによって，被験者の頭部形状に依存する総電場Eと渦電流密度Jを，ともに解析することができる。コイルの巻線形状と電流によって決まる電場E1も，解析の過程で求めることができる。【００１７】
本願発明者は，SPFD法による計算の正確性を検証するため，図２１Ａに示す半径75mmの球モデルと，図２１Ｂに示す辺長120mmの立方体モデルそれぞれに対し，重畳８字コイルにより発生される誘導電流のシミュレーションを行った。「重畳８字コイル」は、２つの渦巻コイルがそれぞれの端部で一部重なるように８字の形に配置されたコイルである。
なお，球モデルと立方体モデルとを用いた理由は，有限要素法，SPFD法の計算アルゴリズム（モデル化の考え方）の違いから，シミュレーション対象の幾何学的性質が及ぼす影響（誤差など）が無いかどうかを調べるためである。
コイルはモデル表面より1cmの位置に，5.3kA，3.4kHzの電流を設定した。図２２Ａに示す，球モデルにおける結果，図２２Ｂに示す，立方体モデルにおける結果にみられるように，有限要素法による計算結果と比較して，モデル表面部から一定長以内の球内に含まれる要素の平均誘導電流量が，SPFD法による計算結果においても，ほぼ差異無く算出されていることが確認された。【００１８】
この解析においては，各個人の頭部MRIデータ（画像データ）のうち，灰白質，白質，脳脊髄液の画像データのみを取り出して解析を行った。図６に示すような脳のMRIデータから灰白質，白質，脳脊髄液の３つの脳組織を抽出したモデルが図７である。青色が灰白質，緑色が白質，赤色が脳脊髄液を表す。実際に数値人体モデルを用いて，皮膚や筋肉などの全ての生体組織を対象としたモデルと灰白質，白質，脳脊髄液の３要素だけを取り出したモデルに同じ条件で電気刺激を与えたときの灰白質表面の渦電流密度の解析結果を図８に示す。図示するように，２つのモデルの間には電流密度に大きな差がなく，渦電流密度の分布も同様であることが分かる。以上より，灰白質，白質，脳脊髄液のみを抽出した解析が有効であることが分かった。【００１９】
個々の患者の脳形状について，大脳の一次運動野上の最適刺激位置における運動閾値と同程度の電気刺激を最適刺激位置周辺に与えたときの，一次運動野での電流密度を解析した。ここで，運動閾値は，被験者の最適刺激位置にコイルを設置して脳に磁気刺激を与えたときに筋肉の反応が50 %以上の確率で見られる磁場強度をいう。具体的には，複数人の被験者に対して頭部のMRI画像を取得し，その被験者の最適刺激位置にコイルを設置して運動閾値を調べた。この実験では，内径25 mm，外径97 mm，高さ6 mm，１０回巻の８字コイルを使用し，パルス幅は280 μsとした。また，最適刺激位置を中心として左右（中心溝に沿う方向）に半径25 mm程度，前後（中心溝に垂直な方向）に半径15 mm程度の楕円形領域の内側の３０箇所にコイルを設置し，それぞれの場所で運動閾値と，MRIデータに対応したコイルの３次元位置情報を記録した。次に，得られた運動閾値と同じ大きさの刺激を脳に与えたときの，最適刺激位置の周囲5 mmの平均渦電流密度を計算した。併せて，最適刺激位置と先に述べた楕円形領域内の２，３ヶ所での運動閾値を調べた。【００２０】
図９に，６人の被験者に対して解析を行った結果を示す。各結果は６人の被験者に対してコイルを最適位置に置き，運動閾値に相当する刺激を与えた際の灰白質表面における渦電流Jと電場Eの分布を示す。数式（７）に示したとおり，渦電流と電場は比例関係にある。解析結果のカラースケールは全ての被験者に対して同じである。この結果から，１次運動野の周辺にも高い電流密度が現れる箇所があり，さらに１次運動野の刺激閾値に相当する電流密度は，被験者によってばらつきがあることがわかる。実際に６人の被験者に対する解析結果に対して標的部位を中心とした半径5 mmにおける平均渦電流密度を求めた結果を表１に示す。脳組織が興奮に達する渦電流密度（被験者の平均値）は17.19 A/m2であった。また，刺激閾値に相当する脳内の渦電流密度は，被験者によってばらつきがあることが分かった。【００２１】
【表１】
2015122369.tif 000010
【００２２】
図１０は，最適刺激位置からずれた場所にコイルを置いた際の刺激閾値（筋肉の反応が50 %以上の確率で見られた電圧振幅の実測値）を，最適刺激位置での運動閾値を100として表したグラフを示す。この結果から，コイルが最適刺激位置からずれることにしたがってtwitch（筋のけいれん，単収縮）を出すために必要な磁場刺激強度が大きくなった。【００２３】
図１１に，コイル位置を変化させたときの最適刺激位置における電流密度の違いを，本発明のシミュレーション方法に従って算出した結果を示す。原点を最適刺激位置とし，最適刺激位置からの“ずれ（移動量）”を，X軸方向を中心溝に平行方向，Y軸を中心溝に垂直方向としてグラフに示した。最適刺激位置と“ずれ”は実測値，電流密度は解析値である。この図から，中心溝に対して平行方向は垂直方向に対してコイルの位置がずれることによる最適刺激位置の渦電流密度の変化への影響が大きく，20 mm離れた場所においては渦電流密度の大きさが半分ほどになることが分かった。【００２４】
図１２に最適刺激位置を刺激したときと20 mm離れたときの脳の渦電流分布を示す。この図からも20 mmコイル位置がずれることにより，脳内の電流分布に大きな差が生じることが分かる。
以上のように，治療コイルの位置を変化させた際の，脳全体を俯瞰した渦電流密度の分布，あるいは脳内特定位置における渦電流密度を本発明にかかるシミュレーション方法を用いて算出する態様の他，それに加えてあるいは独立して，コイルの向きの情報，コイル印加電流の情報，あるいは印加電圧の情報の内の少なくともいずれかを変化させた際の，脳全体を俯瞰した渦電流密度の分布，あるいは脳内特定位置における渦電流密度を，先に説明を行った本発明にかかるシミュレーション方法を用いて算出するようにしてもよい。そのようにシミュレーション結果を算出し，または表示を行うための具体的な構成は別に説明を行う。【００２５】
最適刺激位置の周囲にコイルを置いたときに運動閾値に達する渦電流密度を標的部位に発生させるために必要な電源を解析結果から求めた値と実際に測定した結果との間には数％の差があることが分かった。両者の差の原因としては，渦電流密度を求めた領域が大きかったことが考えられる。そこで，平均渦電流密度を求める関心領域の半径を2 mmから10 mmに変化させた際の，最適刺激位置の周囲にコイルを置いたときに標的部位に運動閾値に達する渦電流を発生させるために必要な電流スルーレート（磁場を発生するためにコイルに正弦波の一周期分の電流を与えるとき，この電流が立ち上がる傾斜（又は時間当たりの電流の増加）をいう。）と実測値の差を各被験者において求めた。結果を図１３に示す。各曲線は各被験者を表す。この図から，関心領域（評価対象領域）の半径は結果に影響を与えないことが分かった。【００２６】
脳内の渦電流分布（解析結果）を図１４に示す。図内の点線は半径5 mmの関心領域を表し，赤い色の部分にはより大きい渦電流が発生している。この図から，関心領域内には渦電流密度が大きく発生している灰白質の部分だけではなく，渦電流があまり発生していない部分も含まれていることがわかる。したがって，コイルを最適位置に置いた際に渦電流が大きく発生した部分のみを関心領域とした結果を比較することで，より実測値と解析結果との差を小さくできる可能性があると考えられる。【００２７】
ところで，個々の患者に対して磁気刺激治療の条件を決定するためには，一次運動野を刺激し，それに反応して例えば指などが動くトゥイッチtwitch（筋のけいれん，単収縮）が出る閾値（運動閾値）を求めることが必要となる。しかし，刺激装置の出力を最大まで上げてもtwitchが現れない患者が存在し，そのような場合には磁気刺激治療を実施できなかった。また，疾患によっては一次運動野以外の部位が刺激の対象となることがあり，そのような部位への刺激ではtwitchを観察できないので，刺激強度は間接的な方法で決定せざるを得なかった。オリジナルソフトウェアの応用例として，このような患者への適用が考えられる。例えば，一次運動野に対する刺激でtwitchを観察できなかった場合，平均的な神経興奮閾値にあたる17.19 A/m2を一次運動野に誘導できるように，コンピュータシミュレーションによって刺激強度を提示することができる。また，うつ病の治療においては，左前頭前野背外側部への経頭蓋磁気刺激において効果が見られるとの報告がある。しかし，前頭前野への刺激ではtwitchを観察できないため，あらかじめ運動野への刺激を行って運動閾値に当たる磁場強度を求め，その磁場強度を基準として前頭前野における治療条件が定められている。本発明によるソフトウェアを用いれば，前頭前野に発生する渦電流密度を解析することにより推定することができるため，前頭前野への渦電流が適正値になるように治療条件を決定することもできるようになる。【００２８】
実際に表１の被験者Dの脳モデルにおいて左前頭前野背外側部に表１と同じ条件(電流の振幅は3431 A)で刺激したときの渦電流分布の解析結果を図１５に示す。図内の赤い枠は実際に刺激を与えた領域を示す。この結果では，先の解析と同様に，刺激位置の周囲半径5 mmにおける平均渦電流密度は21.49 A/m2であった。この結果と表１から，3233 A（数式（８）を参照）の振幅を与えることで運動野において運動閾値に相当する渦電流密度と同じ渦電流を左前頭前野背外側部に発生することができるものと理解できる。
【数４】
2015122369.tif 000011
【００２９】
同じ磁場強度の刺激を実施したとしても，刺激部位が変わると生起される誘導電流は異なる可能性がある。シミュレーションにより，予め対象部位に生起される誘導電流を評価することにより，より適正な刺激強度を設定出来る可能性がある。
以下に，本願発明にかかるシミュレーションの方法で上記の知見を確認した例を示す。
図２３Ａは当該被験者の頭部ＭＲＩ画像であり，矢印は前頭前野における刺激位置を示す。一次運動野を刺激位置とし，治療コイルの電流値を3430Ａとした条件の下で，先に示したシミュレーション方法（但し，球内平均を算出するための球の半径は１０ｍｍを用いた）により脳の誘導電流を算出したところ，図２３Ｂに示す分布となった。また刺激位置における半径１０ｍｍ内の平均渦電流密度は，１９．９Ａ／ｍ2  であった。
次に，同一の被験者について，治療コイルの電流値を同じく3430Ａとした条件の下で同様のシミュレーションを，前頭前野を刺激位置として行ったところ，脳の誘導電流は図２３Ｃに示す分布となった。また刺激位置における半径１０ｍｍ内の平均渦電流密度は，１７．６Ａ／ｍ2  であった。
このように，同一の被験者に対し，同一のコイル電流で刺激を行った場合，脳の刺激部位が異なると，生起される誘導電流分布や，刺激部位における渦電流値が異なることがわかる。
すなわち，一次運動野刺激により決定したRMTを基準とし，例えば刺激強度を110%や120%とする現在一般的に行われている方法では，意図した刺激強度とは異なる刺激が前頭前野に対して加えられている可能性が推測できる。
これらのことから，一次運動野以外の対象部位に関して，予め本願発明の脳内誘導電流シミュレーション手法を用いることにより，治療有効性に貢献可能な新しい刺激強度の検討を行うことが可能となる。【００３０】
＜導入＞
経頭蓋磁気刺激法（TMS）は神経学及び精神医学的な疾病の治療方法として，近年広く用いられてきた手法の一つである。特に，より強度の刺激を断続的に行う反復経頭蓋磁気刺激法(rTMS)は，より明確な治療効果があることで知られている。一方，経頭蓋磁気刺激に用いられるコイルの刺激強度は，意図せぬ副作用やコイルの不要な加熱を防ぐため，同様の治療効果を得つつもできるだけ小さく抑えたい，という課題があった。
この課題を解決する一つのアプローチとして，従来，実際に経頭蓋磁気刺激を行ったいくつかの研究において「被験者個々人に適切なコイル刺激角度」の存在が示されていることに本願発明者は着目した。
即ち，経頭蓋磁気刺激の施術前に，この最適なコイル刺激角度を電磁気的な計算によって明らかにすることで，より少ない電流で刺激を行うことができ，これを施術のサポートとすることができるのではないか，という仮説を本願発明者は設定した。
そして，以下に示す本願発明者による研究結果から，まず脳溝及び脳溝内に脳脊髄液を伴った簡易的な経頭蓋磁気刺激モデルにおいて，脳脊髄液による影響が電磁気学的にどのように作用しているのかということが明らかとなった。更に，これらの結果を元に，脳の各刺激部位における適切な刺激角度を簡易的な計算によって算出する方法を新たに確立したので，以下に説明を行う。【００３１】
＜方法＞
本願発明者は最初に，図２４Ａに示す，脳溝を含む簡易的な脳モデルを作成し，経頭蓋磁気刺激コイルの脳溝に対する回転角による影響をシミュレーションした。なお脳モデルは導電率0.11S/mの灰白質と導電率1.79S/mの2mm厚の脳脊髄液により成り，脳溝は2cmの深さ，2mmの幅として作成した。また刺激コイルは外径5.1cm内径1.1cmの8字形状とし，これに5.3kAの電流を印加した。コイルは脳モデル表面から1cm上に固定し，脳溝に対し0度，30度，45度，60度，90度の角度で刺激を行い，その結果について，コイルの中心点からの球内にふくまれるエレメントの誘導電流を平均するという形で比較を行った。【００３２】
次に，本願発明者は，被験者からのMRI画像として得られた脳形状データを灰白質，白質，脳脊髄液の3要素として形成し，これに対して経頭蓋磁気刺激のシミュレーションを行った。加えて，このＭＲＩ画像に基づいた脳形状データを用いて，脳溝など脳の形状は現実の脳と同じであるものの仮想的に脳脊髄液のみからなる中空モデルを作成し，同様のシミュレーションを行った。磁気刺激コイルは初期刺激方向より10度ずつ180度まで回転させ，3要素モデルにおける「灰白質＋白質 」と「脳脊髄液」それぞれにおいてコイル中心からの球内にふくまれる要素の誘導電流の平均，及び中空モデルにおける「脳脊髄液」の同平均を算出し，比較を行った。【００３３】
＜結果＞
図２４Ａに示す簡易的な脳モデルにおいては，灰白質と脳脊髄液における誘導電流に逆相関(r=-0.99)の関係が見られた。このような関係が見られる理由は，脳脊髄液は頭部において最も導電率の高い物質であり，無視できないコイル損失を生んだことに由来すると考えて良い。特に，図２４Ｂのようにコイルが脳溝に対し平行に当てられた際，コイルによる刺激強度の減少が引き起こされた。従って，最適なコイルの刺激方向は，図２４Ｃのように脳溝に対し垂直な方向であることが推察できる。【００３４】
尚，上に記した，コイルと脳溝との幾何学的な配置の説明において，「コイルが脳溝に対して平行」とは，コイルが生成する誘導電流の主な流れ方向が脳溝の溝方向にして平行であることを意味し，同様に「コイルが脳溝に対して垂直」とは，コイルが生成する誘導電流の主な流れ方向が脳溝の溝方向にして垂直であることを意味している。【００３５】
また，図２５Ａに示す，被験者のＭＲＩ画像から生成された現実の脳モデルでのシミュレーション（図２５Ｂに，誘導電流を平均する計算のための球を示す）においても，灰白質と白質における誘導電流と，脳脊髄液における誘導電流との間でも，逆相関の関係(r=-0.79)がみられ，刺激角度による誘導電流の変化が見られた（図２５Ｄ）。【００３６】
加えて，図２５Ｃに示すような、形状は患者の脳ＭＲＩ画像に由来する現実脳の形状でありながら，仮想的に脳脊髄液のみからなる中空モデルにおける脳脊髄液内の誘導電流について，形状が脳ＭＲＩ画像に由来し且つ灰白質，白質および脳脊髄液の3要素モデル（現実の脳モデル）における脳脊髄液内に発生する誘導電流とほぼ差異のない結果が得られた(ニ乗平均誤差=5.34V/m)。従って，最適なコイルの刺激角度は上記の中空モデルを用いた計算によっても十分な精度のもとで推測が可能であり，この結果，3要素モデルを用いた場合に比べ計算量を70％削減できることがわかった。【００３７】
以上の手法により，運動野以外の領域における神経興奮閾値を推定することができると考えられる。
さらに，コイル位置のずれによる脳内の渦電流密度分布を解析したことにより，コイルの位置の違いにより脳内の渦電流密度は変化し，それにしたがって，運動閾値も変化することが確かめられると同時に，脳内の渦電流密度の分布は異方性があることが分かった。脳内渦電流分布の解析を進めると共に，運動閾値の計測においては最適刺激位置により近い位置における計測も必要であると考えられる。【００３８】
本発明は，以上の知見に基づいてなされたもので，脳内電流シミュレーション方法は，
患者の断層画像データのうち少なくとも脳の一部を含む頭部画像データを提供する第１の工程と，
前記第１の工程で提供された前記頭部画像データのうち前記脳を構成している少なくとも１つの領域を微小要素に分割したそれぞれの微小多面体単位からなる，３次元脳モデルを形成する第２の工程と，
前記患者の頭部上にコイルを配置し，前記コイルに電流を印加することによって前記患者の脳に磁気刺激を与え，前記磁気刺激に対する患者の反応が観察されたときの条件のうち，少なくとも前記コイルの位置及び向きの条件と，前記コイルに印加された電流の条件と，前記コイルの生成磁界に関わる構造の条件と，を含む第１の情報を提供する第３の工程と，
前記３次元脳モデルの前記それぞれの微小多面体単位内に誘導される渦電流又は電界を，前記第３の工程で提供された前記第１の情報，及び前記微小多面体単位ごとに割り付けられた導電率を含む第２の情報に基づいて計算する第４の工程を含む。【００３９】
このような構成を備えた脳内電流シミュレーション方法によれば，個々の患者の脳形状に対応したモデルが形成され，それをもとに脳内電流が解析されるので，実際の脳内に発生する渦電流又は電界の分布を正確に把握することができる。また，コイルの位置が最適刺激位置からずれたときの脳内渦電流分布をより正確に把握できる。そして，それら渦電流分布又は電界分布の解析結果を用いて，より効率的に脳を刺激できる経頭蓋磁気刺激システムを設計できる。【００４０】
本発明の他の形態の脳内電流シミュレーション方法は，前記第４の工程で計算された渦電流又は電界の分布を視覚的に表示する第５の工程を含む。この方法によれば，脳内の渦電流又は電界の分布を視覚的に把握できる。【００４１】
前記それぞれの微小多面体単位は，灰白質，白質，及び脳脊髄液のいずれかの導電率を有するように割り付けられていることが好ましい。また，前記第１の情報は，前記電流の電流値または電圧値の少なくともいずれかを含むことが好ましい。さらに，渦電流又は電界の計算はスカラーポテンシャル有限差分法で行われることが好ましい。【００４２】
本発明に係る脳内電流シミュレーション装置及びその装置を組み入れた経頭蓋磁気刺激システムは，
患者の断層画像データのうち少なくとも脳の一部を含む頭部画像データを提供する第１の手段と，
前記第１の手段で提供された前記頭部画像データのうち前記脳を構成している少なくとも１つの領域を微小要素に分割したそれぞれの微小多面体単位からなる，３次元脳モデルを形成する第２の手段と，
前記患者の頭部上にコイルを配置して前記コイルに電流を印加することによって前記患者の脳に磁気刺激を与え，前記磁気刺激に対する患者の反応が観察されたときの条件のうち，少なくとも前記コイルの位置及び向きの条件と，前記コイルに印加された電流の条件と，前記コイルの生成磁界に関わる構造の条件と，を含む第１の情報を提供する第３の手段と，
前記３次元脳モデルの前記それぞれの微小多面体単位内に誘導される渦電流又は電界を，前記第３の手段で提供された前記第１の情報，及び前記微小多面体単位ごとに割り付けられた導電率を含む第２の情報に基づいて計算する第４の手段を含む。【００４３】
この脳内電流シミュレーション装置及び経頭蓋磁気刺激システムはまた，前記第４の手段で計算された渦電流又は電界の分布を視覚的に表示する第５の手段を含むことが好ましい。【００４４】
さらに，本発明に係る脳内電流シミュレーション装置及びその装置を組み入れた経頭蓋磁気刺激システムは，
患者の断層画像データのうち少なくとも脳の一部を含む頭部画像データを提供する第１の手段と，
前記第１の手段で提供された前記頭部画像データのうち前記脳を構成している少なくとも１つの領域を微小要素に分割したそれぞれの微小多面体単位からなる，３次元脳モデルを形成する第２の手段と，
前記患者の頭部上にコイルを配置して前記コイルに電流を印加することによって前記患者の脳に磁気刺激を与える条件のうち，少なくとも前記コイルの位置及び向きの条件と，前記コイルに印加された電流の条件と，前記コイルの生成磁界に関わる構造の条件と，を含む第１の情報を提供する第３の手段と，
前記３次元脳モデルの前記それぞれの微小多面体単位内に誘導される渦電流又は電界を，前記第３の手段で提供された前記第１の情報，及び前記微小多面体単位ごとに割り付けられた導電率を含む第２の情報に基づいて計算する第４の手段を含む。【００４５】
なお，本発明において，前記コイルに印加される電流は交流又は脈流のいずれであってもよい。【００４６】
ところで，図９に示すように，コイルの磁界によって脳内に誘導される渦電流と電界（電場）強度は比例関係にある。したがって，本願発明に関連する説明において「渦電流」は「電界」又は「電場」に読み換えて読むことができ，そのように読み換えた内容も本願発明の技術的範囲に属するものと理解すべきである。
【図面の簡単な説明】
【００４７】
【図１】患者に磁気刺激治療を施している状態を示す図である。
【図２】８字コイルによって誘導される渦電流を示す図である。
【図３】コイルに電流を印加する駆動回路を示す図である。
【図４】駆動回路から出力される電流波形を示す図である。
【図５】磁気刺激装置の使用状態を示す図である。
【図６】人の脳のＭＲＩデータを示す図である。
【図７】人の脳の灰白質，白質，脳脊髄液の脳組織を抽出したモデルを示す図である。
【図８】人の脳の灰白質，白質，脳脊髄液の脳組織を抽出したモデルに電気刺激を与えたときの灰白質表面の渦電流密度の解析結果を示す図である。
【図９】６人の被験者に対してコイルを最適位置に置き，運動閾値に相当する刺激を与えたときの灰白質表面の電流密度を示す図である。
【図１０】最適刺激位置から移動した場所にコイルを置いたときの，刺激閾値の実測値を，最適刺激位置に対応する運動閾値を１００として表したグラフである。
【図１１】コイルの位置を変化させたときの最適刺激位置における電流密度を示す図である。
【図１２】最適刺激位置とそこから２０ｍｍ離れた位置を刺激したときの脳内の渦電流の分布を示す図。
【図１３】平均渦電流密度を求める関心領域の半径を変化させたとき，標的部位に運動閾値に達する渦電流を発生させるために必要な電流スルーレートと実測値との差を各被験者について示した図である。
【図１４】脳内渦電流分布の解析結果を示す図である。
【図１５】被験者Ｄの脳モデルにおいて左前頭前野外側を刺激したときに発生する脳内渦電流密度の解析結果を示す図。
【図１６】本発明に係る経頭蓋磁気刺激システムの構成を模式的に示す図である。
【図１７】図１６のシステムに組み込まれたコイル駆動回路を示す図である。
【図１８】図１６のシステムに組み込まれた制御回路を示す図である。
【図１９】渦電流密度解析シミュレーションのプロセスを示す図である。
【図２０】３次元脳モデルと，コイルの巻線形状，コイルの位置を示す図である。
【図２１Ａ】半径７５ｍｍの球モデルについて，重畳８字コイルにより発生される誘導電流を示す図である。
【図２１Ｂ】辺長１２０ｍｍの立方体モデルについて，重畳８字コイルにより発生される誘導電流を示す図である。
【図２２Ａ】半径７５ｍｍの球モデルについて，モデル表面から１ｃｍの位置に配置した重畳８字コイルに5.3kA，3.4kHzの電流を加えたとき，モデル表面部から一定長以内の球内に含まれる要素における平均誘導電流量を有限要素法とSPFD法でそれぞれ計算した結果を示すグラフである。
【図２２Ｂ】辺長１２０ｍｍの立方体モデルについて，モデル表面から１ｃｍの位置に配置した重畳８字コイルに5.3kA，3.4kHzの電流を加えたとき，モデル表面部から一定長以内の球内に含まれる要素における平均誘導電流量を有限要素法とSPFD法でそれぞれ計算した結果を示すグラフである。
【図２３Ａ】被験者の頭部ＭＲＩ画像を示す図である。
【図２３Ｂ】図２３Ａに示す脳の一次運動野を刺激位置として治療コイルに所定電流を流したとき，図２３Ａに示す脳に誘導される電流の分布を示す図である。
【図２３Ｃ】図２３Ａに示す脳の前頭前野を刺激位置として治療コイルに所定電流を流したとき，図２３Ａに示す脳に誘導される電流の分布を示す図である。
【図２４Ａ】脳溝を含む簡易的な脳モデルを示す図。
【図２４Ｂ】図２４Ａに示す脳溝に平行にコイルを当てたときの誘導電流の流れ方向を示す図である。
【図２４Ｃ】図２４Ａに示す脳溝に垂直にコイルを当てたときの誘導電流の流れ方向を示す図である。
【図２５Ａ】被験者のＭＲＩ画像から生成された現実の脳モデルを示す図である。
【図２５Ｂ】図２５Ａに示す現実の脳モデルにおいて、誘導電流を平均化する計算のための球を示す図である。
【図２５Ｃ】形状は患者の脳ＭＲＩ画像に由来する現実脳の形状であって、仮想的に脳脊髄液のみからなる中空モデルにおいて、誘導電流を平均化する計算のための球を示す図である。
【図２５Ｄ】灰白質，白質，脳脊髄液における誘導電流と刺激角度との関係を示す図である。
【発明を実施するための形態】
【実施形態１】
【００４８】
図１６を参照すると，経頭蓋磁気刺激システム（以下，単に「システム」という。）１は，支持機構（例えば，椅子２又はベッド）に支持された患者３の脳に磁気刺激を与える磁気刺激装置４を有する。
【００４９】
磁気刺激装置４は，患者３の脳に磁気刺激を加える動磁場を形成するために，コイルユニット（コイル装置）５と制御ユニット６を有する。
【００５０】
図示するように，コイルユニット５は，患者３の頭部表面に沿って自由に移動できるとともに任意の位置に位置決めできるように，適当な位置決めユニット７で支持することが好ましい。コイルユニット５は，コイル８とこのコイル８を囲む電気絶縁材料からなるケーシング９を備えている。ケーシング９は，ケーシング９に一体的に形成されたホルダ１０を備えており，ホルダ１０を介して位置決めユニット７に保持されている。コイル８は，１つの環状コイル，導線を８字状に配した８字コイル（例えば，特開２０１２−１２５５４６号公報に開示のコイル）等の任意の公知のコイルが利用できる。ケーシング９は，３つ又はそれ以上の観察対象（例えば，マーク１１又は突起などのターゲット）を一体的に備えている。これらの観察対象は，患者頭部に対するコイル８の相対的な位置と方向を求めるために利用される。
【００５１】
制御ユニット６は，箱型のハウジング１２を備えている。ハウジング１２は入力部１３と出力部１４を備えている。
【００５２】
入力部１３は，システム１の駆動条件（例えば，コイル８に印加する電圧，電流，周波数）を設定する駆動条件設定部１５；断層画像撮影装置（例えば，ＭＲＩ，ＣＴ，ＰＥＴ）〕１６で生成された人体（特に頭部）断層画像データを受信するデータ受信部１７；コイルユニット５のケーシング９に設けたマーク１１と患者３が装着した眼鏡等の装着品（例えば，眼鏡）又は患者３の皮膚に設けた３つ又はそれ以上の観察対象（例えば，マーク１８又は突起）を同時に撮影するステレオ撮像式光学的３次元位置センシングカメラ（以下，単に「カメラ」という。）１９からの画像データを受信するデータ受信部２０を備えている。図示しないが，カメラ１９は，位置決めユニット７又はシステム１が収容される居室の固定部に取り付けられる。
【００５３】
なお，本件発明において，コイルに印加する電流は，時間とともに周期的に流れの方向が変化する電流（交流）だけでなく，流れる方向が一定で，大きさが周期的に変動する電流（いわゆる「脈流」）も含むものと理解すべきである。
【００５４】
出力部１４は，液晶表示装置等のディスプレイ２１，又はディスプレイを備えたコンピュータ（図示せず）に接続され，制御ユニット６から出力されるデータ（例えば，画像データ）をディスプレイ２１に出力してそこに対応する画像を表示することができるように構成されている。
【００５５】
ハウジング１２の内部には図１７に示すコイル駆動回路２５が収容されており，このコイル駆動回路２５がケーブル２６を介してコイル８と電気的に接続されている。
【００５６】
ハウジング１２の内部にはまた図１８に示す制御回路３０が収容されている。制御回路３０は，中央処理装置（以下，「ＣＰＵ」という。）３１と，ＣＰＵ３１に接続された第１記憶部３２，第２記憶部３３，及び演算部３４を備えている。
【００５７】
第１記憶部３２は種々のソフトウェアを格納している。例えば，断層画像データをもとに３次元脳モデル（３次元マップ）を作成する３次元マッピングソフトウェア４１（例えば，サイバネットシステム株式会社から販売されている医用画像処理ソフトウェア「Real INTAGE」）；入力部１３を介して入力された駆動条件に基づいてコイルに印加する電流を決定するコイル駆動条件決定ソフトウェア４２；カメラ１９で撮影された画像に含まれるマーク（コイルユニットに設けたマークと，患者装着品又は患者に設けたマーク）の情報をもとに患者頭部（患者の脳）に対するコイルの相対的な方向と位置を決定するコイル位置決定ソフトウェア４３；コイル位置決定ソフトウェア４３で決定されたコイルの方向と位置，及びコイル駆動条件決定ソフトウェア４２で決定されたコイルの駆動条件をもとに脳内に発生する渦電流密度を計算するとともに計算した渦電流密度の情報を３次元マッピングソフトウェアで作成された３次元脳モデルに重ねて渦電流密度マップ（図１４参照）を作成する渦電流密度マッピングソフトウェア４３が含まれる。
【００５８】
患者の脳に対するコイルの相対的な位置と方向は，例えば，ＷＯ２００７／１２３１４７Ａに開示された手法により決定することができる。この手法によれば，患者の頭部に固定されたマーク又は患者を固定する器具（例えば，椅子，ベッド）に少なくとも３つのターゲット（患者ターゲット）が取り付けられる。患者の頭部（脳）に対する患者ターゲットの相対的な位置は決められている。したがって，患者ターゲットの位置情報は，患者の頭部断層画像データ（３次元座標データ）に合成される。カメラは，患者ターゲットと，コイルユニット５に固定された少なくとも３つのターゲット（コイルターゲット）を撮影する。撮影された画像はコイル位置決定ソフトウェアによって処理され，患者ターゲットに対するコイルターゲットの相対的な位置と方向が求められるとともに，その情報をもとに患者頭部（脳）に対するコイルの相対的な位置と方向が求められる。なお，患者頭部（脳）に対するコイルの相対的な位置はリアルタイムに計算し，また，計算した結果をディスプレイ２１に表示できる。
【００５９】
演算部３４は，ＣＰＵ３１の指示に基づいて上述のソフトウェアを実行する機能を有する。
【００６０】
第２記憶部３３は種々のデータを格納する。例えば，入力部１３を介して入力された人体（頭部）断層画像データ５１及びコイル駆動条件５２，また，３次元マッピングソフトウェア４１を実行して得られた３次元脳モデルのデータ５３，コイル駆動条件決定ソフトウェア４２を実行して得られたコイル駆動電流のデータ５４，コイル位置決定ソフトウェア４３を実行して得られたコイルの方向と位置のデータ５５，及び渦電流密度マッピングソフトウェア４３を実行して得られた渦電流密度のデータ５６とそれをマッピングした渦電流マップ用のデータ５７が格納される。第２記憶部３３はまた，システム１を使って患者の脳に磁気刺激を与えたときの条件等〔コイルに印加した電流・電圧・周波数，頭部に対するコイルの相対的な位置と方向）とそのとき観察された運動閾値（最適刺激位置にコイルを設置して脳に磁気刺激を与えたときに筋肉の反応が50 %以上の確率で観察された電圧振幅）〕の磁気刺激情報５８が格納されている。
【００６１】
以上の構成を備えたシステム１を用いて患者を治療する場合，カメラ１９で撮影された画像をもとに，コイル位置決定ソフトウェアによって，患者頭部に対するコイル８の位置が求められる。患者頭部に対するコイル８の相対位置はディスプレイ２１に表示される。これにより，患者頭部の目的の場所（例えば，最適刺激位置）にコイル８を設置できる。その後，入力部１５を通じて入力されたコイル駆動条件をもとにコイル駆動回路２５が駆動し，患者３の脳に磁気刺激を与える。コイル駆動回路２５は，図１７に示すように，電源６１の出力電圧を所望電圧に変換する電源回路６２，電源回路６２の出力を昇圧する昇圧回路６３，昇圧回路６３から出力を利用して電荷を蓄積するコンデンサ６４，コンデンサ６４に流れる電流を調整する抵抗６５，コンデンサ６４からの出力を所定のタイミングで動作して交流を形成する半導体スイッチ６６を有し，ＣＰＵ３１の出力に基づいて半導体スイッチ１６を駆動して得られた電流がコイル８に印加される。
【００６２】
次に，システム１を使って，患者の脳内に発生する渦電流密度をシミュレーションする技術を，図１９に示すプロセスに従って説明する。
【００６３】
まず，シミュレーションにあたって，入力部１３を通じて，患者３又は被験者の頭部断層画像データ５１（例えば，ＭＲＩデータ）が入力される（ステップ＃１）。入力された頭部断層画像データ５１は，第２記憶部３３に記憶される。次に，ＣＰＵ３１の指示に基づいて第１記憶部３２に格納されている３次元マッピングソフトウェア４１と第２記憶部３３に格納されている頭部断層画像データ５１を利用して，患者３の脳の３次元脳モデル５３を作成する（ステップ＃２）。このとき作成する３次元脳モデル５３は，脳のすべての部位を対象とする必要はなく，少なくとも灰白質，白質，及び脳脊髄液のいずれか１つであってもよい。作成された３次元脳モデル５３は第２記憶部３３に格納される。３次元脳モデルは，必要に応じて，出力部１４を通じてディスプレイ２１に出力して表示することができる。次に，ＣＰＵ３１は，第２記憶部３３に格納されている，以前にシステム１を使って患者３の脳に磁気刺激を与えたときに得られた磁気刺激情報５８〔コイルに印加した電流・電圧・周波数，頭部に対するコイルの相対的な位置と方向）とコイルの巻線形状データ５９（図２０参照），さらに，そのとき観察された運動閾値（最適刺激位置にコイルを設置して脳に磁気刺激を与えたときに筋肉の反応が50 %以上の確率で観察された電圧振幅）〕を読み出す（ステップ＃３）。また，ＣＰＵ３１は，第１記憶部３２に記憶されている渦電流密度マッピングソフトウェア４３と第２記憶部３３に記憶されている脳の３次元脳モデル５３を読み出し，前記磁気刺激情報５８を参考にして，渦電流密度マッピングソフトウェア４３を使って，３次元脳モデルの各微小多面体単位に誘導される渦電流密度を計算する（ステップ＃４）。この計算は，上述したスカラーポテンシャル有限差分法によって行われる。計算された渦電流密度のデータ５６は第２記憶部３３に記憶される。最後に，ＣＰＵ３１は，計算した渦電流密度のデータ５６をもとに渦電流密度マップ５７（図１４参照）を作成する。図１５に示すように，計算された渦電流密度のデータ５６にはそのレベルに応じた色情報を付加し，この色情報をもって渦電流密度のレベルをディスプレイ２１に表示することが好ましい。
【００６４】
また，上記に説明を行った，第２記憶部３３に格納されている，以前にシステム１を使って患者３の脳に磁気刺激を与えたときに得られた磁気刺激情報５８〔コイルに印加した電流・電圧・周波数，頭部に対するコイルの相対的な位置と方向）を用いて渦電流密度を計算，表示する態様の他に，所定の設定手段を用いて設定入力を行う治療コイルの位置を変化させた際の，脳全体を俯瞰した渦電流密度の分布，あるいは脳内特定位置における渦電流密度を本発明にかかるシミュレーション方法を用いて算出する態様，それに加えてあるいは独立して，コイルの向きの情報，コイル印加電流の情報，あるいは印加電圧の情報の内の少なくともいずれかを変化させた際の，脳全体を俯瞰した渦電流密度の分布，あるいは脳内特定位置における渦電流密度を，上記のシミュレーション方法を用いて算出，表示するようにしてもよい。また，コイルの移動に併せて，脳内特定位置における渦電流密度をリアルタイムに表示することができる。なお，ディスプレイ２１には，複数のシミュレーションの結果と条件を同時に表示することで，それらの結果と条件を対比して観察することができる。
【００６５】
このように，本発明によれば，実際に患者の脳を刺激して得られた磁気刺激情報に基づいて，患者の脳を再現した３次元脳モデルについて渦電流密度を計算して渦電流密度マップを作成することで，刺激を与えた患者の脳のどの部分にどの程度の渦電流が誘導されたか，また，どれほどの渦電流が誘導されたことによって患者に反応（例えば，twitch）が表れたのかを確認できる。
【符号の説明】
【００６６】
１：経頭蓋磁気刺激システム
２：椅子
３：患者
４：磁気刺激装置
５：コイルユニット（コイル装置）
６：制御ユニット
７：位置決めユニット
８：コイル
９：ケーシング
１０：ホルダ
１１：マーク
１２：ハウジング
１３：入力部
１４：出力部
１９：カメラ
３０：制御回路
【図１６】
2015122369.tif 000012
【図１】
2015122369.tif 000013
【図２】
2015122369.tif 000014
【図３】
2015122369.tif 000015
【図４】
2015122369.tif 000016
【図５】
2015122369.tif 000017
【図６】
2015122369.tif 000018
【図７】
2015122369.tif 000019
【図８】
2015122369.tif 000020
【図９】
2015122369.tif 000021
【図１０】
2015122369.tif 000022
【図１１】
2015122369.tif 000023
【図１２】
2015122369.tif 000024
【図１３】
2015122369.tif 000025
【図１４】
2015122369.tif 000026
【図１５】
2015122369.tif 000027
【図１７】
2015122369.tif 000028
【図１８】
2015122369.tif 000029
【図１９】
2015122369.tif 000030
【図２０】
2015122369.tif 000031
【図２１Ａ】
2015122369.tif 000032
【図２１Ｂ】
2015122369.tif 000033
【図２２Ａ】
2015122369.tif 000034
【図２２Ｂ】
2015122369.tif 000035
【図２３Ａ】
2015122369.tif 000036
【図２３Ｂ】
2015122369.tif 000037
【図２３Ｃ】
2015122369.tif 000038
【図２４Ａ】
2015122369.tif 000039
【図２４Ｂ】
2015122369.tif 000040
【図２４Ｃ】
2015122369.tif 000041
【図２５Ａ】
2015122369.tif 000042
【図２５Ｂ】
2015122369.tif 000043
【図２５Ｃ】
2015122369.tif 000044
【図２５Ｄ】
2015122369.tif 000045
【国際調査報告】
2015122369.tif 000046
2015122369.tif 000047
2015122369.tif 000048
2015122369.tif 000049
2015122369.tif 000050
2015122369.tif 000051