(57)【要約】    【課題】  経頭蓋的脳機能測定・刺激法で得られたデータを脳表上に正確に投影、表示できるようにする。３次元頭部画像がない場合は、標準脳の脳表上にデータを投影、表示できるようにする。    【解決手段】  頭表座標を極小距離探索法などを用いて、その直下の脳表座標に変換するとともに、複数の被験者のデータで正規化された標準脳において、頭表座標の脳表投影点の座標とその確率分布を求めた。【選択図】      図２
2005143764.tif 000002
【特許請求の範囲】
【請求項１】
頭表位置に設置した複数のマーカーと脳表面画像とを同時に撮像して得られた画像データを用い、前記データにおける３次元頭部画像上の前記マーカーの頭表上での位置をその直下の脳表位置に投影して投影点の３次元座標位置を求めるステップを備えたことを特徴とする頭表座標を脳表座標に変換する方法。
【請求項２】
複数の被験者について求めた脳表座標を標準脳に正規化するステップをさらに備えた請求項1に記載の方法。
【請求項３】
頭表位置をその直下の脳表位置に投影する前記ステップは、前記各マーカーの位置する頭表位置直下の脳表位置を極小距離探索法により求めるものである請求項1又は２に記載の方法。
【請求項４】
頭表位置をその直下の脳表位置に投影する前記ステップは、前記各マーカーの位置する頭表位置直下の脳表位置を垂線投射法により求めるものである請求項1又は２に記載の方法。
【請求項５】
頭表位置をその直下の脳表位置に投影する前記ステップは、前記各マーカーの位置する頭表位置直下の脳表位置を頭表・脳内部参照点線分連結法により求めるものである請求項1又は２に記載の方法。
【請求項６】
前記各マーカーは、頭表位置に実際に設置したものではなく、頭部形態情報から算出して仮想的に設置したものある請求項１から５のいずれかの方法。
【請求項７】
頭表上の基準点となる位置の直下の脳表位置を請求項1から６のいずれかの方法により予め求めておくとともに、頭表上の基準点との相対的位置から頭表上の任意の点または点の集合の座標を算出し、それらの脳表上の投影点または点の集合を求める方法。
【請求項８】
頭表上の基準点となる位置の直下の脳表位置を請求項1から６のいずれかの方法により予め求めておくとともに、複数の被験者についてのデータからそれらの基準点の脳表上投影点の確率分布を求めておき、それらの基準点を基にして頭表上の任意の点を脳表上に投影した脳表座標とその確率的誤差情報を求めることを特徴とする方法。
【請求項９】
頭表座標と請求項1から７のいずれかの方法により求めた脳表座標から頭表と脳表間の距離分布を求める距離分布測定方法。
【請求項１０】
請求項１から９に記載の方法を実現するソフトウェアプログラム。
【請求項１１】
被検体の頭表面から内部に電磁波又は放射線を照射する照射点と、その照射された電磁波又は放射線と脳との相互作用を前記頭表面上で検出する検出点とを備えたプローブと、前記プローブの検出点が検出した信号をもとに脳の状態を解析するデータ処理装置とを備えた経頭蓋的脳機能測定装置において、
前記データ処理装置は、請求項1から８のいずれかに記載の方法により頭表座標を脳表座標に変換したデータをもって頭表面上の位置を脳表座標に変換する座標変換部を備え、前記検出点が検出した信号による解析データを、前記照射点と検出点とにより定まる頭表面上の位置を前記座標変換部により脳表座標に変換するとともに、前記検出点が検出した信号による解析データを変換された脳表座標上に表示するものであることを特徴とする経頭蓋的脳機能測定装置。
【請求項１２】
前記経頭蓋的脳機能測定装置は光計測装置であり、前記プローブは光を照射する送光点と被検体中を透過及び／又は反射した後に外部に放出される光を受光する受光点とを備えている請求項１１に記載の経頭蓋的脳機能測定装置。
【請求項１３】
前記送光点と受光点とを頭表面上にそれぞれ複数個ずつ配置したマルチチャンネル光計測装置である請求項１２に記載の経頭蓋的脳機能測定装置。
【請求項１４】
前記送光点から被検体に照射される光は近赤外領域の光である請求項１２又は１３に記載の経頭蓋的脳機能測定装置。
【請求項１５】
前記照射点と検出点とにより定まる頭表面上の位置はそれらの２点を結ぶ直線の中央の位置である請求項１１から１４のいずれかに記載の経頭蓋的脳機能測定装置。
【請求項１６】
被検体の頭表面から内部に電磁波又は放射線を照射する照射点と、その照射された電磁波又は放射線と脳との相互作用を前記頭表面上で検出する検出点とを備えたプローブと、前記プローブの検出点が検出した信号をもとに脳の状態を解析するデータ処理装置とを備えた経頭蓋的脳機能測定装置において、
前記プローブにおける照射点と検出点との距離は、請求項９により求められた頭表と脳表間の距離分布に基づいて、照射点からの電磁波又は放射線が脳の表面で相互作用して検出点で検出される感度が最大になるように設定されていることを特徴とする経頭蓋的脳機能測定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、近赤外分光分析法（ＮＩＲＳ）や経頭蓋磁気刺激装置（ＴＭＳ）及びその他の経頭蓋的脳機能測定・刺激法において測定・刺激にかかる脳表座標を推定する方法とそのソフトウェア・プログラムと、この方法により推定された脳表座標を用いて測定・刺激結果を表示する経頭蓋的脳機能測定装置に関するものであり、より詳細には、頭表の任意の点又はその集合としての頭表面に展開された測定データを、その直下の脳表上に投影して表現する方法とそのソフトウェア・プログラムと装置に関するものである。
【背景技術】
【０００２】
経頭蓋的脳機能測定・刺激法において、測定・刺激の対象となる脳表位置を知るためには、その都度、核磁気共鳴装置（ＭＲＩ）等の脳画像撮像法を用いて脳の構造画像を得る必要があるが、このような脳画像撮像は煩雑であり、特殊な設備を要するため、より簡便に頭表と脳表の位置的対応を知る方法が望まれている。
【０００３】
これまでにも脳波研究分野において、頭表上の電極装着位置とその直下にある脳表構造を対応付けようとする研究はなされてきた。これらの方法で最も代表的なものは球形フィッティング法である（非特許文献１，２参照。）。球形フィッティング法は脳表及び頭表を同一中心を持つ球に当てはめる方法である。球形フィッティング法によれば、頭表を表す球の半径をｒ１、脳表を表す球の半径をｒ２とすると、任意の頭表上の点は（ｒ１，Φ，θ）という極座標点として表現でき、この頭表点に対応する脳表点はｒ１をｒ２に置き換えることによって求められる。
【０００４】
この球形フィッティング法は脳波研究に適したものである。しかし、球形フィッティング法は、非球形の頭表及び脳表を球に当てはめるため、空間的に歪みが生じるという問題があった。脳波研究では脳波の信号源推定は３次元空間で低空間分解能に行なわれるため、脳表と頭表の対応にずれがあっても特に実用上の問題はないが、他の経頭蓋的脳機能測定・刺激法においては、さらに高い空間分解能を持つ投影法が望まれている。
ＭＲＩ等の撮像により脳画像が存在する場合でも、任意の頭表上の点とその直下の脳表上の点を対応させる基本的方法は確立されておらず、正確な投影法の開発が望まれている。
【非特許文献１】Towle, V. L., Bolanos, J., Suarez, D., Tan, K., Grzeszczuk, R., Levin, D. N., Cakmur, R., Frank, S. A. and Spire, J. P. 1993. The spatial location of EEG electrodes: locating the best-fitting sphere relative to cortical anatomy. Electroencephalogr. Clin. Neurophysiol. 86: 1-6.
【非特許文献２】Lagerlund, T. D., Sharbrough, F. W., Jack, C. R. Jr., Erickson, B. J., Strelow, D. C., Cicora, K. M. and Busacker, N. E. 1993. Determination of 10-20 system electrode locations using magnetic resonance image scanning with markers. Electroencephalogr. Clin. Neurophysiol. 86: 7-14.
【発明の開示】
【発明が解決しようとする課題】
【０００５】
本発明は、このような技術の現状を鑑みてなされたものであって、第１の目的は、頭表と脳表を含む頭部３次元画像上における頭表上の任意の点に対応する脳表上の対応投影点を求めることにより、経頭蓋的脳機能測定・刺激法で得られたデータを脳表上に正確に投影できるようにすることである。
【０００６】
本発明の第２の目的は、被験者の頭部３次元画像が得られない状態においても、頭表上の基準点を媒介として、頭表上の位置を確率的分布として任意の３次元脳座標上に変換できるようにすることである。
本発明の第３の目的は、経頭蓋的脳機能測定装置において測定データを脳表上に正確に表示できるようにすることである。
【課題を解決するための手段】
【０００７】
第１の目的を達成するための本発明は、３次元頭部画像上の任意の頭表位置をその直下の脳表位置に投影し、投影点の３次元座標位置を求めることを特徴とする頭表座標を脳表座標に変換する方法とそのソフトウェア・プログラムである。
脳と頭表の構造には個人差があるため、本発明の方法で頭表上の点と脳表上の点に対応化がなされた場合でも、異なる被験者間において、経頭蓋的脳機能測定・刺激法で得られたデータを統合するためには、複数被験者間でのデータ統合を可能とできることが望ましい。
【０００８】
そのために、脳表上の位置を標準化可能なフォーマットで表現するために、複数の被験者について求めた脳表座標を標準脳に正規化することが好ましい。
そのような標準脳としては、当該分野で一般的に用いられているＭＮＩ（Montreal Neurological Institute）座標やTalairach座標を利用するのが好ましいが、他の座標系を用いることもできる
【０００９】
任意の頭表位置をその直下の脳表位置に投影する方法をいろいろな角度から検討する中で、極小距離探索法、垂線投射法、頭表・基準点線分連結法により、任意の頭表上の点をその直下の脳表位置に正確に対応付けることが可能であり、ソフトウェア・プログラムとして実行できるという新知見を得た。それぞれの方法は、用途に応じて、最適なものを用いることができる。
【００１０】
また、任意の頭表位置をその直下の脳表位置に投影する方法をいろいろな角度から検討する中で、任意の頭表上の点を、頭表上の基準点に対する相対的位置関係として記述することが可能であり、ソフトウェア・プログラムとして実行できるという新知見を得た。この方法は、単独で用いて任意の点の頭表上に表現することも可能であり、極小距離探索法、垂線投射法、頭表・脳内部参照点線分連結法と組み合わせて、脳表上に投影された位置として表現することも可能である。
【００１１】
本発明の好ましい一態様は、頭表位置をその直下の脳表位置に投影するために、頭表位置に設置した複数のマーカーと脳表面画像とを同時に撮像し、各マーカーの位置する頭表位置直下の脳表位置を極小距離探索法により求める方法に基づくソフトウェア・プログラムである。
【００１２】
極小距離探索法においては、３次元頭部画像上の頭表上の任意の点から半径の異なる等距離球を描き、その球と脳表上の接点を求める。頭表が脳表の凸部近傍に位置する場合は、１点の極小距離を持つ点を求めることができる。頭表が脳表の凹部近傍に位置する場合は２点以上の極小距離を持つ点が求まる。この場合は、これらの点の重心を仮想極小点と想定し、当該頭表点から仮想極小点を通る直線を引き、その直線と脳表の交点を、当該等表点の直下の点と定める。これらの工程は全てソフトウェア・プログラムとして実行可能である。極小距離探索法は特に未処理の頭部画像データを処理する場合に有効である。
【００１３】
極小距離探索法は、未処理の頭部画像データ処理に限定されるものではなく、スムージング処理された頭表、脳表画像処理にも有効である。スムージング処理はＭＲＩｃｒｏ等の一般的に普及している脳画像処理プログラムや、他の３次元画像処理プログラムで行なうことが可能である。
【００１４】
本発明は、極小距離探索法を活用するための、頭表及び脳画像の凸多面体・凸包フィッティング、並びに多面体・包体フィッティング法を含んでいる。頭表及び脳画像を凸多面体・凸包にフィッティングさせた場合、頭表にフィットした凸多面体・凸包の任意の点に対して、脳表にフィットした凸多面体・凸包に最短距離を持つ１点が定まる。このとき得られた最短距離を直線として表し、フィッティング前の頭表及び脳表元画像に重ね合わせて、直線との交点を求めることによって、より正確な脳表、頭表点の座標を求めることも可能である。
【００１５】
一般的に、多面体・包体の任意の点から、その内部にある凸多面体・凸包に対して最短距離を取る点が１点定まることを利用し、頭表を凸多面体・凸包よりも複雑な多面体・包体にフィッティングすることもできる。この場合も、極小距離検索法により、頭表上の任意の点に対応する脳表上の１点を求めることが可能である。
【００１６】
脳表構造は頭表構造に比較して複雑であるため、脳表を多面体・包体にフィッティングし、頭表を凸多面体・凸包にフィッティングすることも有効である。この場合は、脳表にフィットした多面体・包体を拡大し、頭表にフィットした凸多面体・凸包を内部に有するような位相変換を行なうことにより、極小距離探索法を用いて、頭表上の任意の点に対応する脳表上の１点を求めることが可能である。これらの工程は全てソフトウェア・プログラムとして実行可能である。
【００１７】
本発明の好ましい他の態様は、頭表位置をその直下の脳表位置に投影するために、頭表位置に設置した複数のマーカーと脳表面画像とを同時に撮像し、各マーカーの位置する頭表位置直下の脳表位置を垂線投射法により求める方法に基づくソフトウェア・プログラムである。
【００１８】
垂線投射法においては、３次元頭部画像上の頭表の任意の点に接する平面を描き、その平面上の接点から脳表上に垂線を下ろし、脳表との交点を投影点として求める。本方法は、未処理頭部画像及びスムース化された頭部画像にも利用できるが、凸多面体・凸包にフィットした頭表画像を利用する際により有効である。凸多面体・凸包にフィットした頭表画像上の任意の点から脳表上に垂線を下ろす際、脳表の元画像に重ね合わせて、垂線との交点を求めることによって、より正確な脳表、頭表点の座標を求めることも可能である。これらの工程は全てソフトウェア・プログラムとして実行可能である。
【００１９】
本発明の好ましいさらに他の態様は、頭表位置をその直下の脳表位置に投影するために、頭表位置に設置した複数のマーカーと脳表面画像とを同時に撮像し、各マーカーの位置する頭表位置直下の脳表位置を頭表・脳内部参照点線分連結法により求める方法に基づくソフトウェア・プログラムである。
【００２０】
頭表・基準点連結法においては、３次元頭部画像上の頭表の任意の点から、頭表で囲まれる脳内部の参照点に対して直線を描き、その直線と脳表の交点を当該頭表点の投影点と定める。脳内部の参照点は任意の点、又は点の集合である。例えば、頭表あるいは脳表の重心点を参照点に設定することも有効である。あるいは、参照点を前交連等の特定脳構造に設定することも可能である。あるいは、参照点を任意の頭表上の点を中心とする近傍頭表面の重心点に設定し、当該等表点をずらすことにより、基準点の集合を得ることも可能である。本方法は、未処理頭部画像、スムース化された頭部画像、多面体・包体にフィットした頭表画像のいずれにも利用可能である。これらの工程は全てソフトウェア・プログラムとして実行可能である。
【００２１】
本発明において、上記の３次元頭部画像上での任意の頭表点の脳表への投影法は、１点のみならず、点の集合として頭表表面の一定領域さらには全体を投影するように拡張でき、この工程をソフトウェア・プログラムとして実行することもできる。
本発明において、上記の３次元頭部画像上での任意の頭表点を脳表へ投影する方法は、任意の脳表点を頭表上に投影する際にも応用可能であり、この工程をソフトウェア・プログラムとして実行することもできる。
【００２２】
第２の目的を達成するために、これらの方法を用いて、予め頭表上の基準点を各種３次元脳表座標上に表現しておくとともに、複数の被験者についてそれらの基準点および脳表上投影点の確率分布を求めておくことによって、補助的脳画像撮像なしに、任意の頭表位置を任意の３次元座標上に確率的誤差情報を付加した形で表現し、ソフトウェア・プログラムに組み入れることを可能にした。
【００２３】
頭表上の基準点としては、好ましくはInternational Federation of Societies for Electroencephalography and Clinical Neurophysiologyが推奨する国際１０−２０システム基準点、及びその拡張法を用いるが、他の相対脳表座標系も用いることもできる。
本発明において、上記の３次元頭部画像上での任意の頭表点を脳表へ投影する方法を用いて、頭表と脳表の距離分布を得ることも可能である。
【００２４】
頭表上の基準点以外の任意の点についても、基準点との相対的位置関係から記述することができる。好ましくはInternational Federation of Societies for Electroencephalography and Clinical Neurophysiologyが推奨する国際１０−２０システム基準点、及びその拡張法を用いるが、他の相対脳表座標系における基準点も用いることもできる。また、上記では頭表位置に複数のマーカーを設置しているが、マーカーは設置しないで、３次元頭部画像上から耳根部、後頭結節などの特徴点を抽出した後、国際１０−２０システムの基準点、その他の相対座標系の基準点を３次元頭部画像上で求めて、これを仮想的なマーカーとすることもできる。
【００２５】
頭表上の基準点以外の任意の点に対する脳表上の投影点を、頭表上基準点の脳表上投影点との相対的位置関係から記述することができる。頭表上の基準点以外の任意の点を頭表座標系に表現し、これを脳表に投影することも可能であり、また、基準点の脳表上投影点を利用して、脳表上の座標系の位置として記述することも可能である。
【００２６】
第３の目的を達成するための本発明の経頭蓋的脳機能測定装置は、被検体の頭表面から内部に電磁波又は放射線を照射する照射点と、その照射された電磁波又は放射線と脳との相互作用を前記頭表面上で検出する検出点とを備えたプローブと、前記プローブの検出点が検出した信号をもとに脳の状態を解析するデータ処理装置とを備えており、前記データ処理装置は、本発明の方法により頭表座標を脳表座標に変換したデータをもって頭表面上の位置を脳表座標に変換する座標変換部を備え、前記検出点が検出した信号による解析データを、前記照射点と検出点とにより定まる頭表面上の位置を前記座標変換部により脳表座標に変換するとともに、前記検出点が検出した信号による解析データを変換された脳表座標上に表示するものである。
【００２７】
経頭蓋的脳機能測定装置の一例は光計測装置であり、そこではプローブは光を照射する送光点と被検体中を透過及び／又は反射した後に外部に放出される光を受光する受光点とを備えている。
好ましい光計測装置の一例は、送光点と受光点とを頭表面上にそれぞれ複数個ずつ配置したマルチチャンネル光計測装置である。
【００２８】
送光点から被検体に照射される光は、生体に対する透過性の高い近赤外領域の光であることが好ましい。
照射点と検出点とが離れている場合、その検出点が検出する脳表上の位置に対応する頭表面上の位置は照射点と検出点の２点を結ぶ直線の中央の位置とすることができる。
本発明の経頭蓋的脳機能測定装置の他の好ましい態様は、プローブにおける照射点と検出点との距離が、本発明の方法により求められた頭表と脳表間の距離分布に基づいて、照射点からの電磁波又は放射線が脳の表面で相互作用して検出点で検出される感度が最大になるように設定されているものである。
【発明の効果】
【００２９】
本発明により頭表座標を脳表座標に変換することにより、頭表上の任意の点を脳表上に投影することが可能となり、経頭蓋的脳機能測定・刺激法で得られたデータを脳表上に高精度で投影することがソフトウェア・プログラム上で可能となる。このように、本発明は、従来の球形フィッティング法とその応用法に伴う歪みの問題を解決することができる。
【００３０】
本発明は、経頭蓋的脳機能測定・刺激法で得られたデータ処理に極めて高い拡張性をもたらすものである。本発明による方法を用いて、経頭蓋的脳機能測定・刺激法で得られたデータを脳表上に投影することにより、脳機能研究で通常用いられる正規化可能なフォーマットにデータを変換できるようになった。この正規化を経て、複数被験者間でのデータ統合も簡便に行なえるようになり、さらにはこのデータを他の脳機能測定法で得られたデータと統合、比較することも可能となった。このように本発明は、ＮＩＲＳやＴＭＳ等の経頭蓋的脳機能測定・刺激法の普及に対して、基盤となる新技術を提供するものであり、脳機能研究とその臨床応用の発展に大きな貢献をもたらすものである。
【００３１】
また本発明において、頭表上の基準点とそれらの３次元脳表投影点の対応付けを予め行なっておくことによって、３次元頭部画像撮像なしに、経頭蓋的脳機能測定・刺激法で得られたデータを脳表上に投影することが可能となる。しかも、その投影に関する空間誤差情報を明確に把握することが可能となる。
【００３２】
頭表基準点の脳表３次元座標への変換は、複数被験者の測定データを基にして、確率分布として表現できる。また、確率分布の中心を用いれば、１つの基準脳表点に対して、１点の対応脳表３次元座標位置を求めることもできる。
【００３３】
頭表上の任意の点について、近傍３点の頭表基準点との相対位置関係を求めれば、当該基準点の脳表投影点の位置関係を参照して、当該頭表を直下の脳表上に投影させることが可能である。この工程もソフトウェア・プログラム上で実行できる。
このようにして本発明に基づくソフトウェア・プログラムによれば、３次元頭部画像が得られる場合でも得られない場合でも、頭表上の点を脳表上へ投影することが可能となる。３次元頭部画像撮像なしで、任意の頭表位置を脳表上に投影する方法はこれまでなく、本発明は全く新規な効果をもたらすものである。これまでにもＮＩＲＳ測定データの２次元マップを脳表上に表現するソフトウェア・プログラムは存在していたが、２次元マップにおける各ピクセルを３次元座標上のボクセルに正確に対応付ける方法には基づいておらず、本発明はより正確な脳機能画像解析を実現する上で、極めて有効である。
【００３４】
本発明の経頭蓋的脳機能測定装置によれば、解析データを脳表座標上に正しく表示できるようになる。例えば、ＮＩＲＳデータの場合、２次元平面上の点もしくはマップとして測定データが展開されているという現状において、２次元相対座標上に表現されたＮＩＲＳ測定データを正確に３次元空間の脳表上に投影し、３次元座標軸上のデータとして展開することができるので、極めて有効な表示となる。
【発明を実施するための最良の形態】
【００３５】
図１は経頭蓋的脳機能測定装置の実施例としてのマルチチャンネル光計測装置を概略的に表わしたものであり、被検体に光を照射する送光点４と被検体中を透過及び／又は反射した後に外部に放出される光を受光する受光点６とが配置面上にそれぞれ複数個ずつ配置されたプローブ１２と、プローブ１２の受光点６が受光した光をもとに解析するデータ処理装置１８とを備えている。
【００３６】
１４は各送光点に測定用の光を供給する送光点ごとの光源であり、ＬＥＤ（発光ダイオード）やＬＤ（レーザダイオード）などの発光素子を用いることができる。各光源１４からの光はそれぞれ光ファイバなどの導光路を介してそれぞれの送光点４に導かれる。この場合、送光点４は導光路の光送光端面となる。
【００３７】
１６は各受光点が受光した光を検出する検出器であり、フォトダイオードやフォトトランジスタなどの受光素子のほか、光電子増倍管などを使用することができる。各受光点からの光は光ファイバなどの導光路を介して各検出器１６に導かれる。この場合、受光点は導光路の光入射端面となる。
ただし、送光点４に発光素子を配置したり、受光点６に受光素子を配置してもよい。
【００３８】
プローブ１２における送光点４と受光点６の配置は特に限定されるものではなく、ここでは一例として配置面を互いに隣接する三角形からなる領域に分割するグリッド２２上の点に送光点４又は受光点６が位置し、各三角形において一つの辺については送光点４どおし又は受光点６どうしが位置し、他の二つの辺については送光点４と受光点６が１個ずつ位置するように設定されている。この実施例では、グリッド２２により分割される配置面上の領域は正三角形である。
【００３９】
１８はプローブ１２の受光点６が受光した光をもとに解析するデータ処理装置であり、２０はデータ処理装置１８で処理された結果を表示する表示装置である。
送光点４と受光点６は頭表面上に配置される。頭表面上の基準点の位置座標と、それぞれの基準点の位置の直下にある脳表面上の位置座標がデータとしてデータ処理装置１８に保持されている。一対の送光点４と受光点６において、送光点４から照射された光が脳表面で反射し受光点６で受光されたとき、その脳表面上の点はその一対の送光点４と受光点６点を結ぶ頭表面上の点の直下の点である。その脳表面上の点の座標は、データ処理装置１８において頭表面上の基準点の位置座標とその直下にある脳表面上の位置座標との関係を示すデータデータを基にして求められる。
【００４０】
この実施例のマルチチャンネル光計測装置において、送光点４と受光点６の各組のデータ収集は、光源１４と検出器１６の作動を送光点４と受光点６の各組ごとに順次切り替えることにより、互いに混信することなく実行することができる。
【００４１】
図２は、本方法による統合的経頭蓋脳機能測定・刺激データ解析ソフトウェア・プログラムのフローチャートである。このフローチャートで示されたソフトウェア・プログラムは本発明の好ましい実施例である。
【００４２】
いま、解析データとしてＮＩＲＳ二次元マップデータ、任意の測定点のデータ及びＴＭＳの刺激位置データが用意されているものとする。まず、どれかの解析データを選択し、その解析データ中の測定点と頭表基準点との相対的位置を計算する。次に頭表へのデータ投影方法を選択する。頭表へのデータ投影方法としては、極小距離探索法、垂線投影法及び頭表・参照点線分連結法が用意されている。これらのデータ投影方法の詳細は図３〜図５により後述する。投影先脳画像として標準テンプレート又はＭＲＩ画像を選択し、その選択した脳画像の脳表上に、先に選択したデータ投影方法によりデータを投影する。その投影したデータを表示するために、確率的データ表現法として最適推定位置表現法と任意の確率分布表現法の内から選択し、その表現法により脳表での投影データを表示する。
【００４３】
図３により極小距離探索法を説明する。
いま、ＭＲＩ画像データが３次元の画素値データとして用意されているものとする。まず、ＭＲＩ画像データから脳表及び頭表の輪郭となる部分を抽出する。次に、各々の画素値データを座標データに変換して３次元の数値データのリストを得る。さらに、脳表輪郭の座標データについて、脳表輪郭を包み込む、凸包フィッティングを行う。この凸包について、頭表上の任意の点から最小距離を持つような凸包上の点を検索して、線分で結ぶ。この結果得られた最小距離を持つ線分を脳表方向に延長して直線を得る。さらにまた、上記で得られた脳表輪郭の座標データを参照して、この直線と最も近接する脳表輪郭点を探索する。この探索の結果得られた脳表輪郭点を、任意の頭表点に対応する脳表投影点として決定する
【００４４】
図４により垂線投影法を説明する。
いま、ＭＲＩ画像データが３次元の画素値データとして用意されているものとする。まず、ＭＲＩ画像データから脳表及び頭表の輪郭となる部分を抽出する。次に、各々の画素値データを座標データに変換して３次元の数値データを得る。さらに、頭表輪郭の座標データについて、頭表上の任意の点に対する接平面を求める。この接平面に対する法線を生成する。さらにまた、上記で得られた脳表輪郭の座標データを参照して、この法線と最も近接する脳表輪郭点を探索する。この探索の結果得られた脳表輪郭点を、任意の頭表点に対応する脳表投影点として決定する。
【００４５】
図５により頭表・参照点線分連結法を説明する。
いま、ＭＲＩ画像データが３次元の画素値データとして用意されているものとする。また、脳の内部参照点の位置が座標データとして用意されているものとする。まず、ＭＲＩ画像データから脳表及び頭表の輪郭となる部分を抽出する。次に、各々の画素値データを座標データに変換して３次元の数値データを得る。さらに、頭表輪郭の座標データについて、頭表上の任意の点と脳の内部参照点とを結ぶ線分を求める。上記で得られた脳表輪郭の座標データを参照して、この線分と最も近接する脳表輪郭点を探索する。この探索の結果得られた脳表輪郭点を、任意の頭表点に対応する脳表投影点として決定する。
【００４６】
このソフトウェア・プログラムは単独で使用することも可能であるが、一部をモジュールとして他のプログラムに組み込むことも可能であり、あるいは、プログラム全体を他のプログラムに組み込んで、より統合的なソフトウェア・プログラムを構築することも可能である。
【実施例１】
【００４７】
頭表上の国際１０−２０システムの１９基準点について、脳表上の投影点を求めた。すなわち、被験者の頭表上の１０−２０システム基準点にマーカーを置き、ＭＲＩ画像撮像を行い、頭部３次元画像を得た。この頭部３次元画像について、極小距離探索法を用いて、頭表上基準点に対応する脳表上の投影点を求めた。
【００４８】
さらに、この測定を１７人の被験者に対して行い、全ての被験者間のデータを統合して、ＭＮＩ標準脳座標で確率分布として表した。すなわち、各被験者の頭部画像をＭＮＩ標準脳に正規化し、正規化変換後の１０−２０システム基準点脳表投影点の位置の確率分布を標準偏差として表現した。表１にその結果を示す。
【表１】
2005143764.tif 000003
【００４９】
上記のように、本発明によれば、各被験者において、任意の頭表位置を脳表上に高精度で投影させることが実現できた。また、複数被験者間のデータ統合も実現できた。このデータは、ＭＮＩ標準脳座標上に展開された理想的な脳において、各基準点の脳表上の投影点は表に示された座標を中心に分布し、標準偏差を半径として表現される球の中に６１％の点が分布する、という確率的推定を与えるものである。
【００５０】
このような誤差情報は、研究、臨床上の応用において、極めて有用な情報である。すなわち、徒に最適値を提示するのではなく、推定方法の有効範囲を示すことが重要である。本発明によって始めて、ＭＲＩ画像撮像なしに、頭表上の点を脳表に投影させる際の誤差情報が得られた。
【実施例２】
【００５１】
本発明による方法で、図１の測定装置を用いて、実際のＮＩＲＳ測定によって得られた脳活性データを、ＭＲＩ画像撮像なしに、脳表上に投影させた。
脳活性データとしては１０名の被験者のリンゴの皮むき時における前頭前野の酸化ヘモグロビン濃度変化を用いた。ＮＩＲＳの測定プローブの装着位置と国際１０−２０システムの１９基準点に対しての相対位置を計測しておき、ＮＩＲＳ測定プローブの頭表上における位置を得た。
【００５２】
次に、１０名の被験者について、その座標情報を統合し、頭表上におけるＮＩＲＳ測定プローブの平均位置座標を得た。これらの平均位置を、既に脳表への投影位置が確定している国際１０−２０システムの１９基準点との相対的位置情報を参照にして、代表的被験者の脳表画像上に投影させた。図６にその結果を示す。脳活性の度合いは右側の活性バーで表示された酸素化ヘモグロビン濃度の上昇で示されている。
【００５３】
本発明以外の方法により、ＮＩＲＳ測定データを脳表上に重ね合わせるソフトウェアも販売されているが、投影方法の設定が不明瞭であり、見栄えよく重ね合わせるだけの機能にとどまっている。一方、本発明による方法では、実測値に基づいた厳密な投影法を実現しており、研究、臨床目的に使用しうる品質が得られている。
また、表１で示したデータを参照すれば、当該領域においては、５ｍｍから１０ｍｍの空間精度で投影法が有効であるという情報が得られ、得られた結果について、空間誤差を検討することが可能となった。
【実施例３】
【００５４】
本発明による方法で、ＮＩＲＳと機能的核磁気共鳴装置（ｆＭＲＩ）の同時測定データの比較を行なった。ＮＩＲＳの測定プローブと装着位置を示すための位置マーカーを装着した状態で、被験者にリンゴの皮を剥く真似を行なった際の運動野の活性をＮＩＲＳとｆＭＲＩで同時測定した。本発明における極小距離検索法を用いて、ＮＩＲＳ測定領域を被験者の脳表画像上に投影させた。ｆＭＲＩによる測定では、還元型ヘモグロビン濃度に有意な低下がみられた領域を調べた。ＮＩＲＳ測定では、酸化ヘモグロビン濃度の変化を測定した。なお、ｆＭＲＩでは、還元型ヘモグロビン濃度低下を、またＮＩＲＳでは、酸化ヘモグロビン濃度の上昇によって、脳の活性化を測定する。両者の間には定性的な相関があることが知られている。図７にその結果を示す。左側がｆＭＲＩデータであり、濃灰色で示された領域に有意な活性が認められる。右側がＮＲＩＳの活性データで、脳活性の度合いは右側の活性バーで表示された酸素化ヘモグロビン濃度の上昇で示されている。
【００５５】
本発明以外の方法により、ＮＩＲＳ測定データとｆＭＲＩの測定データを同じ脳構造画像に高精度で重ねあわせることに成功した。複数の測定法を用いた脳機能測定は近年極めて盛んになりつつあるが、ＮＩＲＳデータを他の脳機能測定法と比較する標準的な方法はまだ確立されておらず、本発明は、その確立を実現するための基本技術を提供しうることが示された。
【実施例４】
【００５６】
本発明による方法で、頭表と脳表の距離分布を得た。まず、ＭＲＩ撮像により、被験者頭部の３次元画像を得た。次に、その頭表上の全ての点から極小距離探索法を用いて、頭表と脳表の距離を測定し、その距離情報を脳表上にプロットした。図８にその結果を示す。
【００５７】
頭表と脳表の距離情報は、ＮＩＲＳ計測に有用である。一般的にＮＩＲＳ測定において、プローブ間隔（送光点と受光点の間の間隔）を広く設定するほど、脳の深部の活性を測定することが可能である。本発明によって、頭表と脳表の距離を得たことによって、ＮＩＲＳプローブ間隔設定に有用な情報が得られた。また、ＴＭＳの磁気刺激強度は頭表と脳表の距離が１ｃｍ増すごとに３０％減衰することが知られている。ＴＭＳ磁気刺激強度の最適化において、頭表と脳表の距離情報を得ることは極めて有用であり、本発明はこの最適化に有用な情報を提供することが示された。
【産業上の利用可能性】
【００５８】
本発明は、近赤外分光分析法（ＮＩＲＳ）や経頭蓋磁気刺激装置（ＴＭＳ）及びその他の経頭蓋的脳機能測定・刺激法において測定された結果を脳表座標を用いて表示する方法や装置に利用することができる。
【図面の簡単な説明】
【００５９】
【図１】一実施例のマルチチャンネル光計測装置を概略的に示すブロック図である。
【図２】一実施例のソフトウェア・プログラムを概略的に示すフロー・チャート図である。
【図３】一極小距離探索法を概略的に示すフロー・チャート図である。
【図４】一垂線投射法を概略的に示すフロー・チャート図である。
【図５】一頭表・参照点線分連結法を概略的に示すフロー・チャート図である。
【図６】ＭＲＩ画像撮像なしで、頭表基準点との相対位置を参照することによって、脳表上に投影させたＮＩＲＳ脳活性データを示す図である。
【図７】ＮＩＲＳとｆＭＲＩの同時測定時における脳活性データを示す。左側がｆＭＲＩデータで、右側がＮＲＩＳの活性データである。
【図８】頭表と脳表の距離分布を示す図であり、黒点は国際１０−２０基準点を示す。
【符号の説明】
【００６０】
４      送光点
６      受光点
１２      プローブ
１８      データ処理装置
１４      光源
１６      検出器
２０      表示装置
【図１】
2005143764.tif 000004
【図２】
2005143764.tif 000005
【図３】
2005143764.tif 000006
【図４】
2005143764.tif 000007
【図５】
2005143764.tif 000008
【図６】
2005143764.tif 000009
【図７】
2005143764.tif 000010
【図８】
2005143764.tif 000011