JP2002000582A - 磁気的核スピン共鳴により生物の機能的脳活動を時間および位置的に分解して表示する方法および装置 - Google Patents
磁気的核スピン共鳴により生物の機能的脳活動を時間および位置的に分解して表示する方法および装置Info
- Publication number
- JP2002000582A JP2002000582A JP2001146352A JP2001146352A JP2002000582A JP 2002000582 A JP2002000582 A JP 2002000582A JP 2001146352 A JP2001146352 A JP 2001146352A JP 2001146352 A JP2001146352 A JP 2001146352A JP 2002000582 A JP2002000582 A JP 2002000582A
- Authority
- JP
- Japan
- Prior art keywords
- stimulus
- function
- organism
- correlation
- functions
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/4806—Functional imaging of brain activation
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Neurosurgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
(57)【要約】
【課題】 多数の刺激関数に基づく生物内の活動変化を
検知し、それらを開始する刺激関数に正しく対応付ける
多くの並列の刺激関数によりfMRI実験を実行する方法お
よびこの方法を実行するための装置を提供する。 【解決手段】 少なくとも２つの相関関係になくかつ互
いに直交する刺激関数ｆ S1（ｔ）、ｆS2（ｔ）、ｆ
Sn（ｔ）により生物内の生理学的過程を刺激する過程、
核共鳴信号を励起しかつ読出すためパルスシーケンスに
より時間および位置的に分解された核共鳴信号Ｓ
MR（ｔ）を発生する過程、発生した核共鳴信号Ｓ
MR（ｔ）と刺激関数ｆS1（ｔ）、ｆS2（ｔ）、ｆ
Sn（ｔ）に対し各々直交する相関関数ｆorth（ｔ）との
相互相関を求める過程および少なくとも２つの刺激関数
ｆS1（ｔ）、ｆS2（ｔ）、ｆSn（ｔ）の非線形結合に基
づく生物内の活動変化P（τ）を検知すべく相関関数ｆ
orth（ｔ）を変更する過程を含む。
検知し、それらを開始する刺激関数に正しく対応付ける
多くの並列の刺激関数によりfMRI実験を実行する方法お
よびこの方法を実行するための装置を提供する。 【解決手段】 少なくとも２つの相関関係になくかつ互
いに直交する刺激関数ｆ S1（ｔ）、ｆS2（ｔ）、ｆ
Sn（ｔ）により生物内の生理学的過程を刺激する過程、
核共鳴信号を励起しかつ読出すためパルスシーケンスに
より時間および位置的に分解された核共鳴信号Ｓ
MR（ｔ）を発生する過程、発生した核共鳴信号Ｓ
MR（ｔ）と刺激関数ｆS1（ｔ）、ｆS2（ｔ）、ｆ
Sn（ｔ）に対し各々直交する相関関数ｆorth（ｔ）との
相互相関を求める過程および少なくとも２つの刺激関数
ｆS1（ｔ）、ｆS2（ｔ）、ｆSn（ｔ）の非線形結合に基
づく生物内の活動変化P（τ）を検知すべく相関関数ｆ
orth（ｔ）を変更する過程を含む。
Description
【０００１】
【発明の属する技術分野】本発明は、生物内の活動変化
を検知するため磁気的核共鳴を用いて得た核共鳴信号を
時間および位置的に分解して評価するための、生物内に
少なくとも２つの刺激関数の使用の下に生理学的な過程
を刺激する装置および方法に関する。
を検知するため磁気的核共鳴を用いて得た核共鳴信号を
時間および位置的に分解して評価するための、生物内に
少なくとも２つの刺激関数の使用の下に生理学的な過程
を刺激する装置および方法に関する。
【０００２】
【従来の技術】刺激により惹起される人間の大脳皮質内
の脳活動を核スピントモグラフにより検知できることが
確認された。この種の実験は、例えば視覚的な刺激によ
り、また指運動による一次モーターコーテックスの周り
の刺激により実行された。機能的な脳検査は他の技術、
例えばPET（陽子放出トモグラフィ）又はEEGによっても
実行される。しかし核スピントモグラフィによれば、は
るかに良い位置的な分解能が得られる。
の脳活動を核スピントモグラフにより検知できることが
確認された。この種の実験は、例えば視覚的な刺激によ
り、また指運動による一次モーターコーテックスの周り
の刺激により実行された。機能的な脳検査は他の技術、
例えばPET（陽子放出トモグラフィ）又はEEGによっても
実行される。しかし核スピントモグラフィによれば、は
るかに良い位置的な分解能が得られる。
【０００３】典型的な方法で核スピントモグラフィを用
いて実行される実験は所謂fRMI-BOLD実験である。その
際、BOLDはBlood Oxigenation Dependent（血液酸素
含有量に依存する）を意味する。刺激に基づく組織内の
活動は一時的な酸素不足を、組織を囲む血液内に発生す
る。これはオーガニズムスにより検出される。周囲の血
管を介して新しい酸素が供給される。こうして組織内の
急な活動の際に先ず酸素含有量の短時間の軽い低下が生
じ、それにオーガニズムスの調節反応の結果として血液
中の酸素含有量が長時間にわたり減衰する明らかな過振
動が続く。局所的な血液酸素含有量のこの時間的な変化
が、fRMI（機能的な磁気共鳴トモグラフィ）BOLD実験の
際に核スピントモグラフを用いて測定かつ局所化され
る。このような局所的な血液酸素濃度の時間的経過の例
を図４に示す。
いて実行される実験は所謂fRMI-BOLD実験である。その
際、BOLDはBlood Oxigenation Dependent（血液酸素
含有量に依存する）を意味する。刺激に基づく組織内の
活動は一時的な酸素不足を、組織を囲む血液内に発生す
る。これはオーガニズムスにより検出される。周囲の血
管を介して新しい酸素が供給される。こうして組織内の
急な活動の際に先ず酸素含有量の短時間の軽い低下が生
じ、それにオーガニズムスの調節反応の結果として血液
中の酸素含有量が長時間にわたり減衰する明らかな過振
動が続く。局所的な血液酸素含有量のこの時間的な変化
が、fRMI（機能的な磁気共鳴トモグラフィ）BOLD実験の
際に核スピントモグラフを用いて測定かつ局所化され
る。このような局所的な血液酸素濃度の時間的経過の例
を図４に示す。
【０００４】実時間でのデータ取得に関しては、核共鳴
技術における取得速度の制限に基づき時間分解能も制限
される。従って、機能的な撮像のためのデータ取得を、
刺激により開始することが提案されている。その際、刺
激に応じ完全な像データセットに対して必要な生データ
の一部分のみが取得される。従って、データ取得を、脳
活動を始動するタスクの周期的な反復と同期化すること
が提案された。比較可能な方法は、既に心臓運動の“フ
ィルム撮影”のために使用されている。
技術における取得速度の制限に基づき時間分解能も制限
される。従って、機能的な撮像のためのデータ取得を、
刺激により開始することが提案されている。その際、刺
激に応じ完全な像データセットに対して必要な生データ
の一部分のみが取得される。従って、データ取得を、脳
活動を始動するタスクの周期的な反復と同期化すること
が提案された。比較可能な方法は、既に心臓運動の“フ
ィルム撮影”のために使用されている。
【０００５】機能的撮像の問題は、刺激又は脳活動によ
り生ずる信号変化を他の信号変化、例えば運動により惹
起される信号変化から分離することにある。そのため特
に、各画素に対し、刺激関数と、得られた時間的信号経
過との相関係数を計算することが提案された。その際刺
激関数として、周期的に反復し、休止時間により隔てら
れる刺激が使用される。しかし周期的な刺激関数は多く
の欠点を有する。−周期的な擾乱プロセス（例えば心
拍、呼吸）を活動信号から分離できず、“生理学的なノ
イズ”として現れる。反復周期の整数倍の遅延を持つプ
ロセスは同様に正しく認識されない。この実験の延長
は、これらの場合のいずれにおいてもより良い擾乱抑制
に通じない。−さらに周期的な刺激関数は不均等な周波
数スペクトルを持つ。従って若干のスペクトル構成要素
は刺激により弱くしか励起されず、又は全く励起されな
い。このことは系統的な誤りをシステム同定に、即ち数
学的モデルのパラメータの決定に導入してしまう。
り生ずる信号変化を他の信号変化、例えば運動により惹
起される信号変化から分離することにある。そのため特
に、各画素に対し、刺激関数と、得られた時間的信号経
過との相関係数を計算することが提案された。その際刺
激関数として、周期的に反復し、休止時間により隔てら
れる刺激が使用される。しかし周期的な刺激関数は多く
の欠点を有する。−周期的な擾乱プロセス（例えば心
拍、呼吸）を活動信号から分離できず、“生理学的なノ
イズ”として現れる。反復周期の整数倍の遅延を持つプ
ロセスは同様に正しく認識されない。この実験の延長
は、これらの場合のいずれにおいてもより良い擾乱抑制
に通じない。−さらに周期的な刺激関数は不均等な周波
数スペクトルを持つ。従って若干のスペクトル構成要素
は刺激により弱くしか励起されず、又は全く励起されな
い。このことは系統的な誤りをシステム同定に、即ち数
学的モデルのパラメータの決定に導入してしまう。
【０００６】上記問題の解決のためドイツ特許第195 29
639号明細書は、患者の機能的な脳活動を時間および位
置的に分解して表示する方法を提案する。その際、刺激
関数により生理学的な過程が患者内に生起される。刺激
関数は非周期的であり、かつ自己相関関数に可能な限り
僅かな副ピークを有する。核共鳴信号を励起し読出すた
めのパルスシーケンスにより時間および位置的に分解さ
れた核共鳴信号が得られ、像情報に変換される。こうし
て得られた情報と刺激関数との時間的な相関により、時
間および位置的に分解された患者内の活動変化が検知さ
れる。
639号明細書は、患者の機能的な脳活動を時間および位
置的に分解して表示する方法を提案する。その際、刺激
関数により生理学的な過程が患者内に生起される。刺激
関数は非周期的であり、かつ自己相関関数に可能な限り
僅かな副ピークを有する。核共鳴信号を励起し読出すた
めのパルスシーケンスにより時間および位置的に分解さ
れた核共鳴信号が得られ、像情報に変換される。こうし
て得られた情報と刺激関数との時間的な相関により、時
間および位置的に分解された患者内の活動変化が検知さ
れる。
【０００７】必要なＭＲデータは、時間および位置的な
分解を考慮に入れて可能な限り速く取得せねばならな
い。従って先ず第１に速いパルスシーケンスが使用され
る。現在の技術による最も速いＭＲ撮像シーケンスは、
所謂EPI(Echo Planer Imaging)シーケンスである。し
かしこのシーケンスは実行案でしかなく、例えばターボ
スピンエコーシーケンス、FISP又はFLASHシーケンスの
ような他の速いパルスシーケンスも考察の対象となる。
分解を考慮に入れて可能な限り速く取得せねばならな
い。従って先ず第１に速いパルスシーケンスが使用され
る。現在の技術による最も速いＭＲ撮像シーケンスは、
所謂EPI(Echo Planer Imaging)シーケンスである。し
かしこのシーケンスは実行案でしかなく、例えばターボ
スピンエコーシーケンス、FISP又はFLASHシーケンスの
ような他の速いパルスシーケンスも考察の対象となる。
【０００８】EPIシーケンスでは先ず高周波パルスが入
射される。同時に、高周波パルスの周波数スペクトルお
よび層選択勾配の強さに関係して検査対象物の唯一の層
が励起されるように、層選択勾配が作用する。層選択勾
配の正の部分パルスに負の部分パルスが続き、それによ
って正の部分パルスにより惹起されるデフェージングが
再び取り消される。
射される。同時に、高周波パルスの周波数スペクトルお
よび層選択勾配の強さに関係して検査対象物の唯一の層
が励起されるように、層選択勾配が作用する。層選択勾
配の正の部分パルスに負の部分パルスが続き、それによ
って正の部分パルスにより惹起されるデフェージングが
再び取り消される。
【０００９】層選択勾配の負の部分パルスと同時に、２
つのプリフェージングパルスが位相コーディング方向又
は読出し方向に入射される。
つのプリフェージングパルスが位相コーディング方向又
は読出し方向に入射される。
【００１０】続いて交互の極性を有する読出し勾配が入
射される。読出し勾配の交互の符号により核共鳴信号は
常に再びリフェージングをされ、その際、読出し勾配の
各部分パルスの下に信号が生ずる。
射される。読出し勾配の交互の符号により核共鳴信号は
常に再びリフェージングをされ、その際、読出し勾配の
各部分パルスの下に信号が生ずる。
【００１１】信号は、位相を信号から信号へと、信号間
の小さい位相コード化パルスにより進めることによっ
て、それぞれ相異なって位相コード化される。
の小さい位相コード化パルスにより進めることによっ
て、それぞれ相異なって位相コード化される。
【００１２】信号は位相判別をして復調され、かつラス
ター内でディジタル化される。信号毎に、得たディジタ
ル値が生データマトリックスの行内に書込まれる。EPI
法の最も速い変形例、所謂“シングルショットEPI”で
は、単一励起後に、像に対する完全な生データセットを
作成するに十分な多くの信号が取得される。像は公知の
方法で生データマトリックスから二次元のフーリエ変換
で得られる。
ター内でディジタル化される。信号毎に、得たディジタ
ル値が生データマトリックスの行内に書込まれる。EPI
法の最も速い変形例、所謂“シングルショットEPI”で
は、単一励起後に、像に対する完全な生データセットを
作成するに十分な多くの信号が取得される。像は公知の
方法で生データマトリックスから二次元のフーリエ変換
で得られる。
【００１３】機能的な撮像のためには、信号の位置的な
分解だけでなく時間的な分解も行わなければならない。
そのため、次々と異なる時点に対応する像データが得ら
れるように、シーケンスが可能な限り速く反復される。
分解だけでなく時間的な分解も行わなければならない。
そのため、次々と異なる時点に対応する像データが得ら
れるように、シーケンスが可能な限り速く反復される。
【００１４】像データセットの最小の要素は画素と呼ば
れる。機能的な撮像では、一般に従来通常の核スピント
モグラフィ像の際よりも荒い分解能、例えば２５６×２
５６画素の典型的な分解能で済ますことができる。
れる。機能的な撮像では、一般に従来通常の核スピント
モグラフィ像の際よりも荒い分解能、例えば２５６×２
５６画素の典型的な分解能で済ますことができる。
【００１５】図５は、ドイツ特許第195 29 639号明細書
による従来の技術による方法の進行の概要を示す。その
際、核共鳴信号ＡMR（ｋ）を励起しかつ読出すためのパ
ルスシーケンスと刺激関数ｆS1（ｋ）とは互いに無関係
に進行する。これらの双方はシステム・人間／磁気共鳴
トモグラフ５１に作用し、また図示しない中央の制御計
算機によりクロックされるが、例えば刺激関数ｆ
S1（ｋ）によるパルスシーケンスＡMR（ｋ）の始動は行
われない。パルスシーケンスＡMR（ｋ）に基づき生デー
タセットＳMR（ｋ）が、またこれらから再び二次元のフ
ーリエ変換５２により像データセットＢ（ｔ）が取得さ
れる。生データマトリックスＳMR（ｋ）内の各要素又は
像データマトリックスＢ（ｔ）内の各画素に対し、時間
的な信号経過が得られる。続いてこの信号経過Ｓ
MR（ｋ）又はＢ（ｔ）と刺激関数ｆS1（ｋ）との間の相
互相関５４が行われる。刺激関数ｆS1（ｋ）は前もって
遅延τを持つ適当な遅延要素５３を通過する。本方法の
実行のため、生データセットＳMR（ｋ）又は像データセ
ットＢ（ｔ）との相互相関がフーリエ変換５２の前に行
われるか、後に行われるかはどちらでもよいことであ
る。
による従来の技術による方法の進行の概要を示す。その
際、核共鳴信号ＡMR（ｋ）を励起しかつ読出すためのパ
ルスシーケンスと刺激関数ｆS1（ｋ）とは互いに無関係
に進行する。これらの双方はシステム・人間／磁気共鳴
トモグラフ５１に作用し、また図示しない中央の制御計
算機によりクロックされるが、例えば刺激関数ｆ
S1（ｋ）によるパルスシーケンスＡMR（ｋ）の始動は行
われない。パルスシーケンスＡMR（ｋ）に基づき生デー
タセットＳMR（ｋ）が、またこれらから再び二次元のフ
ーリエ変換５２により像データセットＢ（ｔ）が取得さ
れる。生データマトリックスＳMR（ｋ）内の各要素又は
像データマトリックスＢ（ｔ）内の各画素に対し、時間
的な信号経過が得られる。続いてこの信号経過Ｓ
MR（ｋ）又はＢ（ｔ）と刺激関数ｆS1（ｋ）との間の相
互相関５４が行われる。刺激関数ｆS1（ｋ）は前もって
遅延τを持つ適当な遅延要素５３を通過する。本方法の
実行のため、生データセットＳMR（ｋ）又は像データセ
ットＢ（ｔ）との相互相関がフーリエ変換５２の前に行
われるか、後に行われるかはどちらでもよいことであ
る。
【００１６】上記の欠点を避けるため、刺激関数ｆ
S1（ｋ）は周期的であってはならず、かつその自己相関
関数の副ピークを考慮に入れて最適化されていなければ
ならない。例えば適当な２進コードが考慮の対象にな
る。
S1（ｋ）は周期的であってはならず、かつその自己相関
関数の副ピークを考慮に入れて最適化されていなければ
ならない。例えば適当な２進コードが考慮の対象にな
る。
【００１７】このようなコードの自己相関関数内での副
ピークの最小化は、平らなパワースペクトルおよび擾乱
源の関数の長さと共に増大する最適な抑制と等しい意義
がある。さらに自己相関関数内の副ピークの最小化に非
周期性が反映する。このような刺激関数と各画素に対し
MRデータから取得される時間に関係する関数との相互相
関により、MRデータから刺激の影響が抽出される。例え
ば運動（呼吸、心拍、脈動するCSF）による擾乱プロセ
スは相互相関の際に殆ど無くなる。
ピークの最小化は、平らなパワースペクトルおよび擾乱
源の関数の長さと共に増大する最適な抑制と等しい意義
がある。さらに自己相関関数内の副ピークの最小化に非
周期性が反映する。このような刺激関数と各画素に対し
MRデータから取得される時間に関係する関数との相互相
関により、MRデータから刺激の影響が抽出される。例え
ば運動（呼吸、心拍、脈動するCSF）による擾乱プロセ
スは相互相関の際に殆ど無くなる。
【００１８】相互相関の結果P1（τ）は各々の画素に対
してモニターの上に表示される。図６には各画素に対す
る信号の時間的経過の概要を示す。既に図５に関連して
説明したように、相関解析により刺激関数とそれにより
開まる脳活動との間の関連が確立される。この相関解析
の結果P1（τ）は、次いで画素毎にモニター上に表示さ
れる。その際特定の脳領域、即ち対応する画素を選び出
し、かつこれらの脳領域に対し相関解析の結果を考察す
ることができる。
してモニターの上に表示される。図６には各画素に対す
る信号の時間的経過の概要を示す。既に図５に関連して
説明したように、相関解析により刺激関数とそれにより
開まる脳活動との間の関連が確立される。この相関解析
の結果P1（τ）は、次いで画素毎にモニター上に表示さ
れる。その際特定の脳領域、即ち対応する画素を選び出
し、かつこれらの脳領域に対し相関解析の結果を考察す
ることができる。
【００１９】
【発明が解決しようとする課題】本発明の課題は、多数
の刺激関数に帰すべき、考察される生物内の活動変化を
検知し、かつそれらを開始する刺激関数に正しく対応付
ける多くの並列の刺激関数によりfMRI実験を実行するこ
とを可能とする方法を提供することである。さらに本発
明の課題は、この方法を実行するための装置を提供する
ことである。
の刺激関数に帰すべき、考察される生物内の活動変化を
検知し、かつそれらを開始する刺激関数に正しく対応付
ける多くの並列の刺激関数によりfMRI実験を実行するこ
とを可能とする方法を提供することである。さらに本発
明の課題は、この方法を実行するための装置を提供する
ことである。
【００２０】
【課題を解決するための手段】この課題は本発明の請求
項１、８、１３および２０に記載の方法および装置によ
り解決される。本発明の実施例は従属請求項に記載され
ている。
項１、８、１３および２０に記載の方法および装置によ
り解決される。本発明の実施例は従属請求項に記載され
ている。
【００２１】本発明の１つの形態によれば、磁気的核ス
ピン共鳴により生物の機能的脳活動を時間および位置的
に分解して表示するための方法において、 ａ）少なくとも２つの相関関係になくかつ互いに直交す
る刺激関数により生物内の生理学的過程を刺激する過程
と、 ｂ）核共鳴信号を励起しかつ読出すためパルスシーケン
スにより時間および位置的に分解された核共鳴信号を発
生する過程と、 ｃ）各刺激関数に帰すべき生物内の活動変化を検知する
ため発生された核共鳴信号と刺激関数との相互相関を求
める過程と を含むことを特徴とする方法が提案される。
ピン共鳴により生物の機能的脳活動を時間および位置的
に分解して表示するための方法において、 ａ）少なくとも２つの相関関係になくかつ互いに直交す
る刺激関数により生物内の生理学的過程を刺激する過程
と、 ｂ）核共鳴信号を励起しかつ読出すためパルスシーケン
スにより時間および位置的に分解された核共鳴信号を発
生する過程と、 ｃ）各刺激関数に帰すべき生物内の活動変化を検知する
ため発生された核共鳴信号と刺激関数との相互相関を求
める過程と を含むことを特徴とする方法が提案される。
【００２２】少なくとも２つの刺激関数が相関せず、か
つ互いに直交するので、核共鳴信号と刺激関数との後続
の時間的な相関を利用して、刺激関数の１つに帰すべき
生物内の線形の活動変化を一義的に特定の刺激関数に対
応付けることができる。信号の伝搬時間差をなくすた
め、刺激関数は適当な遅延要素を通過する。
つ互いに直交するので、核共鳴信号と刺激関数との後続
の時間的な相関を利用して、刺激関数の１つに帰すべき
生物内の線形の活動変化を一義的に特定の刺激関数に対
応付けることができる。信号の伝搬時間差をなくすた
め、刺激関数は適当な遅延要素を通過する。
【００２３】こうして、本発明の方法により、多くのfM
RI実験を同時に、即ちfMRI実験が検査すべき生物内の生
理学的な過程の際に少なくとも２つの刺激関数により刺
激され、実験の相互の擾乱なしに実行できる。この結
果、実験を実行するのに必要な時間と、その際にかかる
コストの大幅な節減とにつながる。
RI実験を同時に、即ちfMRI実験が検査すべき生物内の生
理学的な過程の際に少なくとも２つの刺激関数により刺
激され、実験の相互の擾乱なしに実行できる。この結
果、実験を実行するのに必要な時間と、その際にかかる
コストの大幅な節減とにつながる。
【００２４】本発明の第２の形態によれば、磁気的核ス
ピン共鳴により生物の機能的脳活動を時間および位置的
に分解して表示する方法において、 ａ）少なくとも２つの相関関係になくかつ互いに直交す
る刺激関数により生物内の生理学的過程を刺激する過程
と、 ｂ）核共鳴信号を励起しかつ読出すためパルスシーケン
スにより時間および位置的に分解された核共鳴信号を発
生する過程と、 ｃ）発生した核共鳴信号と、刺激関数に対してそれぞれ
直交する相関関数との相互相関を求める過程と、 ｄ）少なくとも２つの刺激関数の非線形結合に帰すべき
生物内の活動変化を検知するため相関関数を変更する過
程と を含むことを特徴とする方法が提案される。
ピン共鳴により生物の機能的脳活動を時間および位置的
に分解して表示する方法において、 ａ）少なくとも２つの相関関係になくかつ互いに直交す
る刺激関数により生物内の生理学的過程を刺激する過程
と、 ｂ）核共鳴信号を励起しかつ読出すためパルスシーケン
スにより時間および位置的に分解された核共鳴信号を発
生する過程と、 ｃ）発生した核共鳴信号と、刺激関数に対してそれぞれ
直交する相関関数との相互相関を求める過程と、 ｄ）少なくとも２つの刺激関数の非線形結合に帰すべき
生物内の活動変化を検知するため相関関数を変更する過
程と を含むことを特徴とする方法が提案される。
【００２５】少なくとも２つの刺激関数が相関せず、か
つ互いに直交しており、さらに相関関数も刺激関数に対
して直交するので、本発明によれば、核共鳴信号と刺激
関数に対し各々直交する共通の相関関数との時間的な相
関を利用して、刺激関数の組み合わせに基づく生物内の
より高い次数の活動変化を検知できる。その際相関関数
は、必ずしも単一の特定の関数ではなく、上記条件を満
足する関数の一群として生ずる。刺激関数に対し各々直
交する相関関数の相関の際、線形作用は抑制されるの
で、相関の後に得られる信号は２つ又はより多くの刺激
関数の非線形結合にのみ基づくものである。結果に結び
付く相関関数は後験的に判っているので、刺激関数のど
の組み合わせが生物内で測定された活動変化を生じたか
の帰納的推理が行える。さらに、例えば検査される生物
の意図せざる運動のような擾乱信号が除かれる。信号の
伝搬時間差を除くため、ここでも相関関数を、適当な遅
延要素を通過させる。こうして本発明の方法により、よ
り高い次数の作用、即ち多くの刺激（＝刺激関数）の同
時存在に応答し、かつこうして直接的な刺激には応答し
ない脳機能（所謂より高い次数の脳機能）に基づく生物
内の活動変化を、機能的な磁気共鳴トモグラフィを用い
て検知できる。脳領域が、例えば２つの刺激が同時に存
在する際にのみ応答する場合、刺激の非線形な結合とも
呼ばれる。
つ互いに直交しており、さらに相関関数も刺激関数に対
して直交するので、本発明によれば、核共鳴信号と刺激
関数に対し各々直交する共通の相関関数との時間的な相
関を利用して、刺激関数の組み合わせに基づく生物内の
より高い次数の活動変化を検知できる。その際相関関数
は、必ずしも単一の特定の関数ではなく、上記条件を満
足する関数の一群として生ずる。刺激関数に対し各々直
交する相関関数の相関の際、線形作用は抑制されるの
で、相関の後に得られる信号は２つ又はより多くの刺激
関数の非線形結合にのみ基づくものである。結果に結び
付く相関関数は後験的に判っているので、刺激関数のど
の組み合わせが生物内で測定された活動変化を生じたか
の帰納的推理が行える。さらに、例えば検査される生物
の意図せざる運動のような擾乱信号が除かれる。信号の
伝搬時間差を除くため、ここでも相関関数を、適当な遅
延要素を通過させる。こうして本発明の方法により、よ
り高い次数の作用、即ち多くの刺激（＝刺激関数）の同
時存在に応答し、かつこうして直接的な刺激には応答し
ない脳機能（所謂より高い次数の脳機能）に基づく生物
内の活動変化を、機能的な磁気共鳴トモグラフィを用い
て検知できる。脳領域が、例えば２つの刺激が同時に存
在する際にのみ応答する場合、刺激の非線形な結合とも
呼ばれる。
【００２６】両方法において、本発明で用いる刺激関数
が非周期的で、かつその自己相関に可能な限り僅かな副
ピークを持つのが望ましい。そうすれば信号の時間的な
相関により、特に心拍又はパルスシーケンスに基づく周
期的な擾乱信号を除ける。
が非周期的で、かつその自己相関に可能な限り僅かな副
ピークを持つのが望ましい。そうすれば信号の時間的な
相関により、特に心拍又はパルスシーケンスに基づく周
期的な擾乱信号を除ける。
【００２７】具体的には、本発明では刺激関数として２
進コードを使用するとよい。これは特に刺激関数が生物
に、活動に関する指令を前以て与えるときに有利であ
る。
進コードを使用するとよい。これは特に刺激関数が生物
に、活動に関する指令を前以て与えるときに有利であ
る。
【００２８】２進コードは生物の簡単な刺激を可能とす
る（明白な信号作用関係、例えば“光信号の際に小指を
動かす”）。さらに、公知の方法で容易に刺激関数と直
交する関数を導き出せる。
る（明白な信号作用関係、例えば“光信号の際に小指を
動かす”）。さらに、公知の方法で容易に刺激関数と直
交する関数を導き出せる。
【００２９】核共鳴信号を励起しかつ読出すパルスシー
ケンスと、刺激関数との反作用を除くため、本発明で
は、パルスシーケンスと刺激関数とを互いに無関係にす
る。
ケンスと、刺激関数との反作用を除くため、本発明で
は、パルスシーケンスと刺激関数とを互いに無関係にす
る。
【００３０】２つの刺激の非線形な結合を検知する場
合、本発明は、少なくとも２つの刺激関数の非線形結合
に基づく生物内の活動変化の検知を、統計的な方法を用
いて実行する。その際、相互相関の結果は、例えば“hi
gher order statistics”のような適当な統計的な方
法を用いてさらに評価される。これらの方法は、特に信
号内の非線形な項の検出および評価のために適してい
る。
合、本発明は、少なくとも２つの刺激関数の非線形結合
に基づく生物内の活動変化の検知を、統計的な方法を用
いて実行する。その際、相互相関の結果は、例えば“hi
gher order statistics”のような適当な統計的な方
法を用いてさらに評価される。これらの方法は、特に信
号内の非線形な項の検出および評価のために適してい
る。
【００３１】２つの刺激の非線形結合を検知するには、
相関関数が主要な役割を演ずる。本発明では、シフトレ
ジスタを用いて相関関数を変更する。シフトレジスタを
用いて、全てが個々の刺激関数に各々直交する適当な一
群の関数を通過させる。
相関関数が主要な役割を演ずる。本発明では、シフトレ
ジスタを用いて相関関数を変更する。シフトレジスタを
用いて、全てが個々の刺激関数に各々直交する適当な一
群の関数を通過させる。
【００３２】本発明の第３および第４の形態では、上記
の方法を実行する装置を開示する。
の方法を実行する装置を開示する。
【００３３】
【発明の実施の形態】以下、図１ないし６により本発明
を一層詳細に説明する。
を一層詳細に説明する。
【００３４】本発明を実現するfMRI実験の一般的な構成
を以下に簡潔に説明する。
を以下に簡潔に説明する。
【００３５】図３による概要図に相応し、患者３１が、
基本磁界を発生する磁石３２内に入れられる。この磁石
３２は、図を見易くするために図示しない基本磁界発生
のための内蔵物と、高周波信号の送信と受信のための高
周波アンテナとを有する。勾配コイルと高周波アンテナ
は、プロセッサユニット３４と接続されている。ＭＲ設
備の機能は、その構成が本質的に従来通常の設備と一致
しているので、周知であると仮定している。
基本磁界を発生する磁石３２内に入れられる。この磁石
３２は、図を見易くするために図示しない基本磁界発生
のための内蔵物と、高周波信号の送信と受信のための高
周波アンテナとを有する。勾配コイルと高周波アンテナ
は、プロセッサユニット３４と接続されている。ＭＲ設
備の機能は、その構成が本質的に従来通常の設備と一致
しているので、周知であると仮定している。
【００３６】パルス発生器３５から刺激関数が発生さ
れ、それにより例えば光発生器３３が駆動される。しか
し例えば電気的な刺激も使えるし、又は患者が、例えば
光学的な信号により、刺激関数に相応して運動を行うよ
う促される。プロセッサユニット３４により取得された
ＭＲデータおよびパルス発生器３５内で発生された刺激
関数は、相関計算ユニット３６内で相互の相関を求めら
れる。こうして得られたデータは、プロセス計算機３９
内でさらに編集され、モニター３８上に表示される。そ
の際、像信号への磁気共鳴トモグラフの出力信号の変換
は、公知の方法で二次元フーリエ変換を用いて行われ
る。
れ、それにより例えば光発生器３３が駆動される。しか
し例えば電気的な刺激も使えるし、又は患者が、例えば
光学的な信号により、刺激関数に相応して運動を行うよ
う促される。プロセッサユニット３４により取得された
ＭＲデータおよびパルス発生器３５内で発生された刺激
関数は、相関計算ユニット３６内で相互の相関を求めら
れる。こうして得られたデータは、プロセス計算機３９
内でさらに編集され、モニター３８上に表示される。そ
の際、像信号への磁気共鳴トモグラフの出力信号の変換
は、公知の方法で二次元フーリエ変換を用いて行われ
る。
【００３７】プロセッサユニット３４、相関計算ユニッ
ト３６およびパルス発生器３５は、中央の制御計算機３
７により制御される。
ト３６およびパルス発生器３５は、中央の制御計算機３
７により制御される。
【００３８】この構成により、患者／生物内で刺激関数
を用いて生理学的過程を刺激し、かつ患者／生物内の刺
激関数の結果としての活動変化を検知するfMRI実験が実
行可能である。
を用いて生理学的過程を刺激し、かつ患者／生物内の刺
激関数の結果としての活動変化を検知するfMRI実験が実
行可能である。
【００３９】図１は相互の擾乱なしに多くのfMRI実験を
並列的に、同時に実行するための本発明による方法の進
行の概要を示す。
並列的に、同時に実行するための本発明による方法の進
行の概要を示す。
【００４０】患者体内に生理学的過程を引起こすため、
人間／MRシステム１１（MRは磁気共鳴トモグラフを表
す）が、相関関係になくかつ互いに直交する２つの刺激
関数ｆ S1（ｔ）、ｆS2（ｔ）により励起される。システ
ム１１の出力端に、詳細は示さない内部処理ステップに
従い、核共鳴信号ＳMR（ｔ）が得られる。これら核共鳴
信号ＳMR（ｔ）は、２つの別々の相互相関１４、１５内
で各々刺激関数ｆS1（ｔ）、ｆS2（ｔ）の１つとの時間
的な相関を求められる。システムに付随する伝搬時間差
を等化するため、刺激関数ｆS1（ｔ）、ｆS2（ｔ）は先
ず各々遅延τ1、τ2を持つ適当な遅延回路１２、１３を
通過する。両方の相互相関１４、１５の出力端に生ずる
信号Ｐ1（τ1）、Ｐ2（τ2）は各刺激、即ち各刺激関数
ｆS1（ｔ）、ｆS2（ｔ）に基づく人間／患者内の活動変
化を示す。
人間／MRシステム１１（MRは磁気共鳴トモグラフを表
す）が、相関関係になくかつ互いに直交する２つの刺激
関数ｆ S1（ｔ）、ｆS2（ｔ）により励起される。システ
ム１１の出力端に、詳細は示さない内部処理ステップに
従い、核共鳴信号ＳMR（ｔ）が得られる。これら核共鳴
信号ＳMR（ｔ）は、２つの別々の相互相関１４、１５内
で各々刺激関数ｆS1（ｔ）、ｆS2（ｔ）の１つとの時間
的な相関を求められる。システムに付随する伝搬時間差
を等化するため、刺激関数ｆS1（ｔ）、ｆS2（ｔ）は先
ず各々遅延τ1、τ2を持つ適当な遅延回路１２、１３を
通過する。両方の相互相関１４、１５の出力端に生ずる
信号Ｐ1（τ1）、Ｐ2（τ2）は各刺激、即ち各刺激関数
ｆS1（ｔ）、ｆS2（ｔ）に基づく人間／患者内の活動変
化を示す。
【００４１】２つの刺激関数ｆS1（ｔ）、ｆS2（ｔ）は
相関関係になくかつ互いに直交するので、核共鳴信号Ｓ
MR（ｔ）と刺激関数ｆS1（ｔ）、ｆS2（ｔ）との時間的
な相関に基づき、刺激関数ｆS1（ｔ）、ｆS2（ｔ）の１
つに基づく人間／患者内の線形の活動変化を一義的に特
定の刺激関数ｆS1（ｔ）、ｆS2（ｔ）に対応付けること
ができる。こうして本発明の方法によって、２つのfMRI
実験を同時に、相互の擾乱なしに実行できる。それによ
り従来技術と比べて、実験を行うのに必要な時間と、そ
の際にかかるコストとの大幅な節減が可能である。
相関関係になくかつ互いに直交するので、核共鳴信号Ｓ
MR（ｔ）と刺激関数ｆS1（ｔ）、ｆS2（ｔ）との時間的
な相関に基づき、刺激関数ｆS1（ｔ）、ｆS2（ｔ）の１
つに基づく人間／患者内の線形の活動変化を一義的に特
定の刺激関数ｆS1（ｔ）、ｆS2（ｔ）に対応付けること
ができる。こうして本発明の方法によって、２つのfMRI
実験を同時に、相互の擾乱なしに実行できる。それによ
り従来技術と比べて、実験を行うのに必要な時間と、そ
の際にかかるコストとの大幅な節減が可能である。
【００４２】本発明による方法は、図３に示す装置によ
り実現される。その際、人間／ＭＲシステム１１は、患
者３１、磁石３２およびプロセッサユニット３４から成
る組み合わせに相当する。システム５１は、刺激関数ｆ
S1（ｔ）、ｆS2（ｔ）をパルス発生器３５と、１つ又は
それ以上の光発生器３３とを用いて発生する。核共鳴信
号ＳMR（ｔ）と刺激関数ｆS1（ｔ）、ｆS2（ｔ）との相
互相関１４、１５は相関計算ユニット３６内で行われ
る。遅延回路１２、１３はパルス発生器３５内で実現し
てもよい。例えばプロセス計算機３９内で行われるフー
リエ変換の後に、こうして得たデータがモニター３８の
上に表示される。プロセッサユニット３４、相関計算ユ
ニット３６およびパルス発生器３５の制御は、ここでも
中央の制御計算機３７により行われる。
り実現される。その際、人間／ＭＲシステム１１は、患
者３１、磁石３２およびプロセッサユニット３４から成
る組み合わせに相当する。システム５１は、刺激関数ｆ
S1（ｔ）、ｆS2（ｔ）をパルス発生器３５と、１つ又は
それ以上の光発生器３３とを用いて発生する。核共鳴信
号ＳMR（ｔ）と刺激関数ｆS1（ｔ）、ｆS2（ｔ）との相
互相関１４、１５は相関計算ユニット３６内で行われ
る。遅延回路１２、１３はパルス発生器３５内で実現し
てもよい。例えばプロセス計算機３９内で行われるフー
リエ変換の後に、こうして得たデータがモニター３８の
上に表示される。プロセッサユニット３４、相関計算ユ
ニット３６およびパルス発生器３５の制御は、ここでも
中央の制御計算機３７により行われる。
【００４３】図２は、磁気的核スピン共鳴により、所謂
より高い次数の活動変化、即ちより多くの刺激が同時に
存在する際にのみ応答する患者内の活動変化を検知する
ことを可能とする本発明による方法の進行の概要を示
す。
より高い次数の活動変化、即ちより多くの刺激が同時に
存在する際にのみ応答する患者内の活動変化を検知する
ことを可能とする本発明による方法の進行の概要を示
す。
【００４４】患者内の生理学的刺激過程を開始するた
め、人間／MRシステム２１が相関関係になくかつ互いに
直交する多くの刺激関数ｆS1（ｔ）、ｆS2（ｔ）、ｆSn
（ｔ）により励起される。本発明で用いる刺激関数ｆS1
（ｔ）、ｆS2（ｔ）、ｆSn（ｔ）は非周期的であり、か
つその自己相関に可能な限り僅かな副ピークを有する。
この実施例では刺激関数として２進コードを使用する。
患者と磁気共鳴トモグラフとを含むシステム２１の出力
端に核共鳴信号ＳMR（ｔ）が得られる。これら核共鳴信
号ＳMR（ｔ）は相互相関２３内で、刺激関数ｆ
S1（ｔ）、ｆS2（ｔ）、ｆ Sn（ｔ）に対し各々直交する
相関関数ｆorth（ｔ）との相互相関を求められる。信号
の伝搬時間差を除くため、相関関数ｆorth（ｔ）は遅延
τを持つ適当な遅延要素を通過する。相関関数ｆ
orth（ｔ）とは、個々の刺激関数ｆS1（ｔ）、ｆ
S2（ｔ）、ｆSn（ｔ）に各々直交する一群の関数から成
る任意の関数である。それは公知の数学的方法を用いて
求められる。類似の関数を決定する試みは“Advanced
Methods of Physiological System Modeling”（第
３巻、V.Z.Marmarelis編、ニューヨーク、Plenum、第８
７〜１１０頁；１９９４）の出版物、EthanA.Benardete
およびJonathan D.Victor “An extension of the
m-sequencetechnique for the analysis of Mul
ti-input Nonlinear systems”に記載されている。勿
論一群の適当な関数は他の方法でも導出できる。一群の
適当な相関関数ｆorth（ｔ）は、例えばシフトレジスタ
を通過し、こうして相関関数ｆorth（ｔ）が変更され
る。
め、人間／MRシステム２１が相関関係になくかつ互いに
直交する多くの刺激関数ｆS1（ｔ）、ｆS2（ｔ）、ｆSn
（ｔ）により励起される。本発明で用いる刺激関数ｆS1
（ｔ）、ｆS2（ｔ）、ｆSn（ｔ）は非周期的であり、か
つその自己相関に可能な限り僅かな副ピークを有する。
この実施例では刺激関数として２進コードを使用する。
患者と磁気共鳴トモグラフとを含むシステム２１の出力
端に核共鳴信号ＳMR（ｔ）が得られる。これら核共鳴信
号ＳMR（ｔ）は相互相関２３内で、刺激関数ｆ
S1（ｔ）、ｆS2（ｔ）、ｆ Sn（ｔ）に対し各々直交する
相関関数ｆorth（ｔ）との相互相関を求められる。信号
の伝搬時間差を除くため、相関関数ｆorth（ｔ）は遅延
τを持つ適当な遅延要素を通過する。相関関数ｆ
orth（ｔ）とは、個々の刺激関数ｆS1（ｔ）、ｆ
S2（ｔ）、ｆSn（ｔ）に各々直交する一群の関数から成
る任意の関数である。それは公知の数学的方法を用いて
求められる。類似の関数を決定する試みは“Advanced
Methods of Physiological System Modeling”（第
３巻、V.Z.Marmarelis編、ニューヨーク、Plenum、第８
７〜１１０頁；１９９４）の出版物、EthanA.Benardete
およびJonathan D.Victor “An extension of the
m-sequencetechnique for the analysis of Mul
ti-input Nonlinear systems”に記載されている。勿
論一群の適当な関数は他の方法でも導出できる。一群の
適当な相関関数ｆorth（ｔ）は、例えばシフトレジスタ
を通過し、こうして相関関数ｆorth（ｔ）が変更され
る。
【００４５】少なくとも２つの刺激関数ｆS1（ｔ）、ｆ
S2（ｔ）、ｆSn（ｔ）が相関関係になくかつ互いに直交
しており、かつ相関関数ｆorth（ｔ）も刺激関数ｆ
S1（ｔ）、ｆS2（ｔ）、ｆSn（ｔ）に対しそれぞれ直交
するので、本発明に従い、核共鳴信号ＳMR（ｔ）と相関
関数ｆorth（ｔ）の１つとの時間的な相関を利用して相
互相関２３の出力端に、刺激関数ｆS1（ｔ）、ｆ
S2（ｔ）、ｆSn（ｔ）の組み合わせに基づくより高い次
数の活動変化を人間／MRシステム２１に描出する信号Ｐ
（τ）が得られる。相関２３における線形の作用は、刺
激関数ｆS1（ｔ）、ｆS2（ｔ）、ｆSn（ｔ）に対して各
々直交する相関関数ｆorth（ｔ）により抑制されるの
で、得られた有意義な信号Ｐ（τ）は２つ又はそれ以上
の多くの刺激関数ｆS1（ｔ）、ｆS2（ｔ）、ｆSn（ｔ）
の非線形の結合のみに基づくものである。
S2（ｔ）、ｆSn（ｔ）が相関関係になくかつ互いに直交
しており、かつ相関関数ｆorth（ｔ）も刺激関数ｆ
S1（ｔ）、ｆS2（ｔ）、ｆSn（ｔ）に対しそれぞれ直交
するので、本発明に従い、核共鳴信号ＳMR（ｔ）と相関
関数ｆorth（ｔ）の１つとの時間的な相関を利用して相
互相関２３の出力端に、刺激関数ｆS1（ｔ）、ｆ
S2（ｔ）、ｆSn（ｔ）の組み合わせに基づくより高い次
数の活動変化を人間／MRシステム２１に描出する信号Ｐ
（τ）が得られる。相関２３における線形の作用は、刺
激関数ｆS1（ｔ）、ｆS2（ｔ）、ｆSn（ｔ）に対して各
々直交する相関関数ｆorth（ｔ）により抑制されるの
で、得られた有意義な信号Ｐ（τ）は２つ又はそれ以上
の多くの刺激関数ｆS1（ｔ）、ｆS2（ｔ）、ｆSn（ｔ）
の非線形の結合のみに基づくものである。
【００４６】刺激関数ｆS1（ｔ）、ｆS2（ｔ）、ｆ
Sn（ｔ）と、信号Ｐ（τ）を生み出す後験的に知り得る
相関関数ｆorth（ｔ）との時間的経過を考慮に入れ、刺
激関数のどの組み合わせが生物内の測定された活動変化
に基づくかの帰納的推理が行われる。こうして本発明の
提案する方法により、より高い次数の作用、即ち多くの
刺激が同時に存在する際に応答する作用を検知すること
が可能である。この場合、刺激の非線形な結合とも呼ば
れる。
Sn（ｔ）と、信号Ｐ（τ）を生み出す後験的に知り得る
相関関数ｆorth（ｔ）との時間的経過を考慮に入れ、刺
激関数のどの組み合わせが生物内の測定された活動変化
に基づくかの帰納的推理が行われる。こうして本発明の
提案する方法により、より高い次数の作用、即ち多くの
刺激が同時に存在する際に応答する作用を検知すること
が可能である。この場合、刺激の非線形な結合とも呼ば
れる。
【００４７】この実施例では、相互相関２３の出力信号
P（τ）内に含まれる非線形の項は統計的な方法２４（H
OS＝higher order statistics）を用いて抑制され
る。信号内の非線形の項を評価するこの方法は出版物
“Introducing higher orderstatistics（HOS） for
the detection of non-linearities”（UK Nonli
near News、１９９５年９月）、“Non-linearity det
ection for conditionmonitoring using higher o
rder statistics”および“The higher order stat
istics of speech signals゛(゛Techniques for S
peech SignalProcessing and their Application
s”に関するIEEEコロキューム、ロンドン、１９９４年
６月１日に発表；Digest No．１９９４／１３８ ｐ７
／１〜７／６）に開示されている。これら３つの論文の
主著者はJ.W.A.FackrellおよびS.McLaughlinである。
P（τ）内に含まれる非線形の項は統計的な方法２４（H
OS＝higher order statistics）を用いて抑制され
る。信号内の非線形の項を評価するこの方法は出版物
“Introducing higher orderstatistics（HOS） for
the detection of non-linearities”（UK Nonli
near News、１９９５年９月）、“Non-linearity det
ection for conditionmonitoring using higher o
rder statistics”および“The higher order stat
istics of speech signals゛(゛Techniques for S
peech SignalProcessing and their Application
s”に関するIEEEコロキューム、ロンドン、１９９４年
６月１日に発表；Digest No．１９９４／１３８ ｐ７
／１〜７／６）に開示されている。これら３つの論文の
主著者はJ.W.A.FackrellおよびS.McLaughlinである。
【００４８】実験の評価を簡単化するため、相互相関の
結果を相互相関図として表示すると有利である。これ
は、公知の方法で、相互相関の出力端に得られる信号の
二次元のフーリエ変換により行われる。この表示形態は
得られたデータの直観性を顕著に高め、かつ結果の迅速
かつ正確な評価を可能とする。
結果を相互相関図として表示すると有利である。これ
は、公知の方法で、相互相関の出力端に得られる信号の
二次元のフーリエ変換により行われる。この表示形態は
得られたデータの直観性を顕著に高め、かつ結果の迅速
かつ正確な評価を可能とする。
【００４９】本発明による方法は、図３の装置で実行で
きる。ここで人間／MRシステム２１は患者３１、磁石３
２およびプロセッサユニット３４から成る組み合わせに
相当する。刺激関数ｆS1（ｔ）、ｆS2（ｔ）を、システ
ム５１はパルス発生器３５と、１つ又はそれ以上の光発
生器３３とを用いて励起する。勿論音響的、熱的および
／又は電動的な発生器も使用できる。核共鳴信号Ｓ
MR（ｔ）と相関関数ｆorth（ｔ）との相互相関２３は、
相関計算ユニット３６内で行われる。相関関数ｆ
or th（ｔ）はプロセッサユニット３４（又は中央の制御
計算機３７）内で導き出され、かつ変更される。遅延回
路２２はパルス発生器３５（又は相関計算ユニット３
６）内で実現される。例えばプロセス計算機３９内で行
われるフーリエ変換の前又は後に、相互相関２３内で取
得された信号Ｐ（τ）は、データをモニター３８上に表
示する前に、プロセス計算機３９内で統計的な方法によ
りさらに処理される。プロセッサユニット３４、相関計
算ユニット３６およびパルス発生器３５の制御は、ここ
でも中央の制御計算機３７により行われる。勿論プロセ
ス計算機３９の制御も中央の制御計算機３７により行わ
れる。
きる。ここで人間／MRシステム２１は患者３１、磁石３
２およびプロセッサユニット３４から成る組み合わせに
相当する。刺激関数ｆS1（ｔ）、ｆS2（ｔ）を、システ
ム５１はパルス発生器３５と、１つ又はそれ以上の光発
生器３３とを用いて励起する。勿論音響的、熱的および
／又は電動的な発生器も使用できる。核共鳴信号Ｓ
MR（ｔ）と相関関数ｆorth（ｔ）との相互相関２３は、
相関計算ユニット３６内で行われる。相関関数ｆ
or th（ｔ）はプロセッサユニット３４（又は中央の制御
計算機３７）内で導き出され、かつ変更される。遅延回
路２２はパルス発生器３５（又は相関計算ユニット３
６）内で実現される。例えばプロセス計算機３９内で行
われるフーリエ変換の前又は後に、相互相関２３内で取
得された信号Ｐ（τ）は、データをモニター３８上に表
示する前に、プロセス計算機３９内で統計的な方法によ
りさらに処理される。プロセッサユニット３４、相関計
算ユニット３６およびパルス発生器３５の制御は、ここ
でも中央の制御計算機３７により行われる。勿論プロセ
ス計算機３９の制御も中央の制御計算機３７により行わ
れる。
【００５０】本発明で説明する方法および装置により、
多くの並列な刺激関数を有するfMRI実験を実行すること
ができる。考察する生物内の多数の刺激関数に基づく活
動変化を検知し、かつこれらの活動変化を、それらを開
始する刺激関数に正しく対応付けることが可能である。
多くの並列な刺激関数を有するfMRI実験を実行すること
ができる。考察する生物内の多数の刺激関数に基づく活
動変化を検知し、かつこれらの活動変化を、それらを開
始する刺激関数に正しく対応付けることが可能である。
【００５１】より正確に言えば、本発明が提案する方法
および装置により、多くの個々のfMRI実験を同時に並列
に実行することも、所謂より高い次数の脳機能、即ち直
接的な刺激には応答しない脳機能を検出することも可能
である。
および装置により、多くの個々のfMRI実験を同時に並列
に実行することも、所謂より高い次数の脳機能、即ち直
接的な刺激には応答しない脳機能を検出することも可能
である。
【図１】本発明方法の第１実施例の概要図。
【図２】本発明方法の第２実施例の概要図。
【図３】機能的な撮像のためのＭＲ設備の構成を示す概
要図。
要図。
【図４】組織内での急な活動時周囲の血液における血中
酸素含有量の経過を示す概要図。
酸素含有量の経過を示す概要図。
【図５】従来方法の概要図。
【図６】個々の画素内の信号経過を示す概要図。
【符号の説明】 １１、２１ 人間／MRシステム １２、１３、２２ 遅延回路 １４、１５、２３ 相互相関 ２４ 統計的な方法 ３１ 患者 ３２ 磁石 ３３ 光発生器 ３４ プロセッサユニット ３５ パルス発生器 ３６ 相関計算ユニット ３７ 制御計算機 ３８ モニター ３９ プロセス計算機
Claims (24)
- 【請求項１】 磁気的核スピン共鳴により生物の機能的
脳活動を時間および位置的に分解して表示する方法にお
いて、 ａ）少なくとも２つの相関関係になくかつ互いに直交す
る刺激関数（ｆs1（ｔ）、ｆs2（ｔ）、ｆsn（ｔ））に
より生物内の生理学的過程を刺激する過程と、 ｂ）核共鳴信号を励起しかつ読出すためパルスシーケン
スにより時間および位置的に分解した核共鳴信号（ＳMR
（ｔ））を発生する過程と、 ｃ）発生した核共鳴信号（ＳMR（ｔ））と、刺激関数
（ｆs1（ｔ）、ｆs2（ｔ）、ｆsn（ｔ））に対してそれ
ぞれ直交する相関関数（ｆorth（ｔ））との相互の相関
を求める過程と、 ｄ）少なくとも２つの刺激関数（ｆs1（ｔ）、ｆ
s2（ｔ）、ｆsn（ｔ））の非線形結合に帰すべき生物内
の活動変化（Ｐ（τ））を検知するため相関関数（ｆ
orth（ｔ））を変更する過程とを含むことを特徴とする
方法。 - 【請求項２】 刺激関数（ｆs1（ｔ）、ｆs2（ｔ）、ｆ
sn（ｔ））が非周期的であり、かつその自己相関に可能
な限り僅かな副ピークを有することを特徴とする請求項
１記載の方法。 - 【請求項３】 刺激関数（ｆs1（ｔ）、ｆs2（ｔ）、ｆ
sn（ｔ））として２進コードを使用することを特徴とす
る請求項１又は２記載の方法。 - 【請求項４】 刺激関数（ｆs1（ｔ）、ｆs2（ｔ）、ｆ
sn（ｔ））が生物に対し、活動に関する指令を前もって
与えることを特徴とする請求項１ないし３の１つに記載
の方法。 - 【請求項５】 少なくとも２つの刺激関数（ｆ
s1（ｔ）、ｆs2（ｔ）、ｆsn（ｔ））の非線形結合に帰
すべき生物内の活動変化（Ｐ（τ））の検知を、統計的
な方法（HOS）を用いて行うことを特徴とする請求項１
ないし４の１つに記載の方法。 - 【請求項６】 相関関数（ｆorth（ｔ））を、シフトレ
ジスタを用いて変更することを特徴とする請求項１ない
し５の１つに記載の方法。 - 【請求項７】 核共鳴信号を励起しかつ読出すためのパ
ルスシーケンスと、刺激関数（ｆs1（ｔ）、ｆ
s2（ｔ）、ｆsn（ｔ））とが互いに無関係であることを
特徴とする請求項１ないし６の１つに記載の方法。 - 【請求項８】 磁気的核スピン共鳴により生物の機能的
脳活動を時間および位置的に分解して表示する方法にお
いて、 ａ）少なくとも２つの相関関係になくかつ互いに直交す
る刺激関数（ｆs1（ｔ）、ｆs2（ｔ）、ｆsn（ｔ））に
より生物内の生理学的過程を刺激する過程と、 ｂ）核共鳴信号を励起しかつ読出すためパルスシーケン
スにより時間および位置的に分解された核共鳴信号（Ｓ
MR（ｔ））を発生する過程と、 ｃ）各刺激関数（ｆs1（ｔ）、ｆs2（ｔ））に基づく生
物内の活動変化（Ｐ1（τ1）、Ｐ2（τ2））を検知する
ため、発生した核共鳴信号（ＳMR（ｔ））と刺激関数
（ｆs1（ｔ）、ｆs2（ｔ））との相互相関を求める過程
とを含むことを特徴とする方法。 - 【請求項９】 刺激関数（ｆs1（ｔ）、ｆs2（ｔ））が
非周期的であり、かつその自己相関に可能な限り僅かな
副ピークを有することを特徴とする請求項８記載の方
法。 - 【請求項１０】 刺激関数（ｆs1（ｔ）、ｆs2（ｔ））
として２進コードを使用することを特徴とする請求項８
又は９記載の方法。 - 【請求項１１】 刺激関数（ｆs1（ｔ）、ｆs2（ｔ））
が、生物に対し活動に関する指令を前もって与えること
を特徴とする請求項８ないし１０の１つに記載の方法。 - 【請求項１２】 核共鳴信号を励起しかつ読出すためパ
ルスシーケンスおよび刺激関数（ｆs1（ｔ）、ｆ
s2（ｔ））が互いに無関係であることを特徴とする請求
項８ないし１１の１つに記載の方法。 - 【請求項１３】 磁気的核スピン共鳴により生物の機能
的脳活動を時間および位置的に分解して表示する装置に
おいて、 ａ）少なくとも２つの相関関係になくかつ互いに直交す
る刺激関数（ｆs1（ｔ）、ｆs2（ｔ）、ｆsn（ｔ））に
より生物内の生理学的過程を刺激するための関数発生器
（３３、３５）と、 ｂ）時間および位置的に分解された核共鳴信号（Ｓ
MR（ｔ））を発生するMR設備（３２、３４）と、 ｃ）MR設備（３２、３４）を用いて得た核共鳴信号（Ｓ
MRｔ））と、刺激関数（ｆs1（ｔ）、ｆs2（ｔ）、ｆsn
（ｔ））に対し各々直交する相関関数（ｆor th（ｔ））
との相互相関を求める相互相関計算ユニット（３６）
と、 ｄ）少なくとも２つの刺激関数（ｆs1（ｔ）、ｆ
s2（ｔ）、ｆsn（ｔ））の非線形結合に帰すべき生物
（３１）内の活動変化（Ｐ（τ））を検知するための評
価装置（３９）とを含むことを特徴とする装置。 - 【請求項１４】 刺激関数（ｆs1（ｔ）、ｆs2（ｔ）、
ｆsn（ｔ））が非周期的であり、かつその自己相関に可
能な限り僅かな副ピークを有することを特徴とする請求
項１３記載の装置。 - 【請求項１５】 刺激関数（ｆs1（ｔ）、ｆs2（ｔ）、
ｆsn（ｔ））として２進コードが使用されることを特徴
とする請求項１３又は１４記載の装置。 - 【請求項１６】 関数発生器（３３、３５）が生物（３
１）に対し活動に関する指令を前もって与えることを特
徴とする請求項１３ないし１５の１つに記載の装置。 - 【請求項１７】 少なくとも２つの刺激関数（ｆ
s1（ｔ）、ｆs2（ｔ）、ｆ sn（ｔ））の非線形結合に帰
すべき生物内の活動変化（Ｐ（τ））の検知が統計的な
方法（HOS）を用いて行われることを特徴とする請求項
１３ないし１６の１つに記載の装置。 - 【請求項１８】 評価装置（３９）および相関計算ユニ
ット（３６）がプロセス計算機又はパーソナルコンピュ
ータであることを特徴とする請求項１３ないし１７の１
つに記載の装置。 - 【請求項１９】 核共鳴信号を励起しかつ読出すためパ
ルスシーケンスおよび刺激関数（ｆs1（ｔ）、ｆ
s2（ｔ））が互いに無関係であることを特徴とする請求
項１３ないし１８の１つに記載の装置。 - 【請求項２０】 磁気的核スピン共鳴により生物の機能
的脳活動を時間および位置的に分解して表示する装置に
おいて、 ａ）少なくとも２つの相関関係になくかつ互いに直交す
る刺激関数（ｆs1（ｔ）、ｆs2（ｔ）、ｆsn（ｔ））に
より生物内の生理学的過程を刺激するための関数発生器
（３３、３５）と、 ｂ）時間および位置的に分解された核共鳴信号（Ｓ
MR（ｔ））を発生するＭＲ設備（３２、３４）と、 ｃ）ＭＲ設備（３２、３４）を用いて得た核共鳴信号
（ＳMRｔ））と、刺激関数（ｆs1（ｔ）、ｆs2（ｔ）、
ｆsn（ｔ））に対し各々直交する相関関数（ｆ
orth（ｔ））との相互相関を求める相関計算ユニット
（３６）と、 ｄ）相関関数（ｆorth（ｔ））を変更するための制御ユ
ニット（３７）と、 ｅ）少なくとも２つの刺激関数（ｆs1（ｔ）、ｆ
s2（ｔ）、ｆsn（ｔ））の非線形結合に帰すべき生物
（３１）内の活動変化（Ｐ（τ））を検知する評価装置
（３９）とを含むことを特徴とする装置。 - 【請求項２１】 刺激関数（ｆs1（ｔ）、ｆs2（ｔ）、
ｆsn（ｔ））が非周期的であり、かつその自己相関に可
能な限り僅かな副ピークを有することを特徴とする請求
項２０記載の装置。 - 【請求項２２】 刺激関数（ｆs1（ｔ）、ｆs2（ｔ）、
ｆsn（ｔ））として２進コードが使用されることを特徴
とする請求項２０又は２１記載の装置。 - 【請求項２３】 関数発生器（３３、３５）が生物（３
１）に対し、活動に関する指令を前もって与えることを
特徴とする請求項２０ないし２２の１つに記載の装置。 - 【請求項２４】 核共鳴信号を励起しかつ読出すためパ
ルスシーケンスおよび刺激関数（ｆs1（ｔ）、ｆ
s2（ｔ））が互いに無関係であることを特徴とする請求
項２０ないし２３の１つに記載の装置。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| DE10024488.2 | 2000-05-18 | ||
| DE10024488A DE10024488C2 (de) | 2000-05-18 | 2000-05-18 | fMRI-BOLD Experiment mit multiplen Stimulationsmustern |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2002000582A true JP2002000582A (ja) | 2002-01-08 |
| JP4804647B2 JP4804647B2 (ja) | 2011-11-02 |
Family
ID=7642608
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001146352A Expired - Fee Related JP4804647B2 (ja) | 2000-05-18 | 2001-05-16 | 磁気的核スピン共鳴により生物の機能的脳活動を時間および位置的に分解して表示する方法および装置 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US6539246B2 (ja) |
| JP (1) | JP4804647B2 (ja) |
| DE (1) | DE10024488C2 (ja) |
Families Citing this family (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE10040850C2 (de) * | 2000-08-21 | 2003-01-16 | Siemens Ag | Verfahren zum Betreiben eines Kernspintomographiegerätes mit verbessertem Offresonanzverhalten einer True-Fisp-Meßsequenz in Gegenwart zweier Spinkollektive und Kernspintomographiegerät zur Durchführung des Verfahrens |
| US20020103429A1 (en) * | 2001-01-30 | 2002-08-01 | Decharms R. Christopher | Methods for physiological monitoring, training, exercise and regulation |
| US20050283053A1 (en) * | 2002-01-30 | 2005-12-22 | Decharms Richard C | Methods for physiological monitoring, training, exercise and regulation |
| CA2434966A1 (en) * | 2001-01-30 | 2002-08-08 | Christopher R. Decharms | Methods for physiological monitoring, training, exercise and regulation |
| GB0122049D0 (en) * | 2001-09-12 | 2001-10-31 | Nycomed Imaging As | Method |
| US20040092809A1 (en) * | 2002-07-26 | 2004-05-13 | Neurion Inc. | Methods for measurement and analysis of brain activity |
| JP4241913B2 (ja) * | 2003-03-26 | 2009-03-18 | 株式会社日立製作所 | 訓練支援装置 |
| US20050033154A1 (en) * | 2003-06-03 | 2005-02-10 | Decharms Richard Christopher | Methods for measurement of magnetic resonance signal perturbations |
| US20050267357A1 (en) * | 2003-10-21 | 2005-12-01 | Rao Stephen M | fMRI system for detecting symptoms associated with Attention Deficit Hyperactivity Disorder |
| US7112965B2 (en) * | 2004-04-30 | 2006-09-26 | University Of Basel | Low-impact noise acquisition magnetic resonance imaging |
| AU2005307870A1 (en) * | 2004-11-15 | 2006-05-26 | Christopher Decharms | Stimulation of neural tissue with light |
| DE102005010093A1 (de) * | 2005-03-04 | 2006-09-14 | Siemens Ag | Verfahren zur Darstellung von Respirationsmustern |
| US20070167724A1 (en) * | 2005-12-09 | 2007-07-19 | Gadagkar Hrishikesh P | fMRI data acquisition system |
| WO2007124320A2 (en) * | 2006-04-19 | 2007-11-01 | University Of Florida Research Foundation, Inc. | Calibration of functional mri of living tissue involving a bold or boss contrast simulation phantom |
| US8600696B2 (en) | 2007-09-28 | 2013-12-03 | Oliver Zafiris | Method and system for determining a reaction signal for a selected location in an information processing system following the effect of at least one input signal |
| US8666475B2 (en) * | 2008-10-22 | 2014-03-04 | The Trustees Of Columbia University In The City Of New York | Images of language-sensitive neurocircuitry as a diagnostic for autism |
| WO2010048438A1 (en) * | 2008-10-22 | 2010-04-29 | The Trustees Of Columbia University In The City Of New York | Images of language-sensitive neurocircuitry as a diagnostic for autism |
| US8379947B2 (en) | 2010-05-28 | 2013-02-19 | International Business Machines Corporation | Spatio-temporal image reconstruction using sparse regression and secondary information |
| DE102012215718B4 (de) * | 2012-09-05 | 2022-05-12 | Siemens Healthcare Gmbh | Verfahren und Magnetresonanzanlage zur MR-Bildgebung eines vorbestimmten Volumenabschnitts eines lebenden Untersuchungsobjekts mittels Stimulieren des Untersuchungsobjekts |
| JP6773493B2 (ja) * | 2016-09-14 | 2020-10-21 | 株式会社東芝 | 検出装置、検出方法、および検出プログラム |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09117430A (ja) * | 1995-08-11 | 1997-05-06 | Siemens Ag | 患者の脳の機能活動の時間及び局所分解表示方法及び装置 |
| JPH09220213A (ja) * | 1996-02-19 | 1997-08-26 | Siemens Ag | 機能性画像生成方法及び核スピントモグラフィ装置 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5603322A (en) * | 1993-01-19 | 1997-02-18 | Mcw Research Foundation | Time course MRI imaging of brain functions |
| JP3512482B2 (ja) * | 1994-09-06 | 2004-03-29 | 株式会社東芝 | 磁気共鳴映像装置 |
| US6104943A (en) * | 1997-11-14 | 2000-08-15 | The Mclean Hospital Corporation | Phased array echoplanar imaging system for fMRI |
| DE19802337C2 (de) * | 1998-01-23 | 2000-05-31 | Forschungszentrum Juelich Gmbh | Computer zur Auswertung von Signalen aus der kernmagnetischen Resonanztomographie sowie mit dem Computer ausgestatteter Kernresonanztomograph |
| US6169403B1 (en) * | 1998-04-29 | 2001-01-02 | Siemens Aktiengesellschaft | Method and MR device for simulating electrical simulations in a subject by MR stimulation |
| DE19860037C2 (de) * | 1998-12-23 | 2001-01-25 | Siemens Ag | Verfahren zur ortsaufgelösten Messung der elektrischen Aktivität von Nervenzellen mittels magnetischer Resonanz und Vorrichtung zur Durchführung des Verfahrens |
| DE19962850B4 (de) * | 1999-12-24 | 2006-12-07 | Forschungszentrum Jülich GmbH | Spektroskopisches Bildgebungsverfahren |
| DE10100441C2 (de) * | 2001-01-08 | 2003-01-09 | Siemens Ag | Verfahren zum Betrieb eines Magnetresonanzgeräts bei Initiierung einer neuronalen Aktivität und Magnetresonanzgerät zum Durchführen des Verfahrens |
| DE10105387C2 (de) * | 2001-02-06 | 2003-04-10 | Siemens Ag | Verfahren zum Betrieb eines Magnetresonanzgeräts zur funktionellen Bildgebung sowie Magnetresonanzgerät zur Durchführung des Verfahrens |
-
2000
- 2000-05-18 DE DE10024488A patent/DE10024488C2/de not_active Expired - Fee Related
-
2001
- 2001-05-16 JP JP2001146352A patent/JP4804647B2/ja not_active Expired - Fee Related
- 2001-05-18 US US09/861,347 patent/US6539246B2/en not_active Expired - Lifetime
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09117430A (ja) * | 1995-08-11 | 1997-05-06 | Siemens Ag | 患者の脳の機能活動の時間及び局所分解表示方法及び装置 |
| JPH09220213A (ja) * | 1996-02-19 | 1997-08-26 | Siemens Ag | 機能性画像生成方法及び核スピントモグラフィ装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| US6539246B2 (en) | 2003-03-25 |
| DE10024488A1 (de) | 2001-11-29 |
| JP4804647B2 (ja) | 2011-11-02 |
| US20020045814A1 (en) | 2002-04-18 |
| DE10024488C2 (de) | 2003-04-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2002000582A (ja) | 磁気的核スピン共鳴により生物の機能的脳活動を時間および位置的に分解して表示する方法および装置 | |
| JP3638382B2 (ja) | 患者の脳の機能活動の時間及び局所分解表示装置 | |
| US5565777A (en) | Method/apparatus for NMR imaging using an imaging scheme sensitive to inhomogeneity and a scheme insensitive to inhomogeneity in a single imaging step | |
| CN110031786B (zh) | 磁共振图像重建方法、磁共振成像方法、设备及介质 | |
| US5377680A (en) | MRI cardiac image produced by temporal data sharing | |
| US6683454B2 (en) | Shifting of artifacts by reordering of k-space | |
| US7689263B1 (en) | Method and apparatus for acquiring free-breathing MR images using navigator echo with saturation RF pulse | |
| US5435303A (en) | MRA image produced by temporal flow data sharing | |
| US20090148021A1 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
| JPH0350545B2 (ja) | ||
| JPS61155740A (ja) | フ−リエ変換ｎｍｒ像作成法における運動ア−テイフアクトの低減方法 | |
| US10191132B2 (en) | Method and apparatus for multi-slice imaging of T2-relaxation time | |
| US6889071B2 (en) | Acquisition of high-temporal free-breathing MR images | |
| US6298258B1 (en) | Method and apparatus for spatially resolved measurement of the electrical activity of nerve cells using magnetic resonance | |
| US6845261B2 (en) | System for correlation of MR images with physiological data | |
| US20140303482A1 (en) | Magnetic resonance imaging method for imaging components with short transverse relaxation times (t2) in a human or an animal heart | |
| JP2523470B2 (ja) | 核磁気共鳴イメ−ジング方式 | |
| JP2005185732A (ja) | 磁気共鳴イメージング装置および磁気共鳴イメージング装置の信号処理方法 | |
| JPH08182661A (ja) | 磁気共鳴撮影方法 | |
| Lundervold et al. | Functional magnetic resonance imaging of primary visual processing using a 1.0 Tesla scanner | |
| JP3512115B2 (ja) | 磁気共鳴を用いた検査装置 | |
| Bresler et al. | Patient-adaptive spatio-temporal MRI: From PARADIGM to PARADISE and beyond | |
| JP2001340316A (ja) | 磁気共鳴イメージング装置 | |
| Sampath | Magnetic resonance pulse sequences for rapid imaging of myocardial motion and strain | |
| Al-Dayeh | Reduction of artifacts in functional magnetic resonance imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080512 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110322 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110621 |
|
| RD03 | Notification of appointment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7423 Effective date: 20110621 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20110712 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20110810 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140819 Year of fee payment: 3 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |