(57)【特許請求の範囲】
【請求項１】
一様な静磁場中に置かれた被検体に高周波磁場と勾配磁場を所定のパルスシーケンスに従って印加し、前記被検体からの磁気共鳴信号を検出して前記被検体の内部情報を第１の磁気共鳴画像と第２の磁気共鳴画像として映像化する磁気共鳴映像装置において、
前記第１の磁気共鳴画像の前記第２の磁気共鳴画像各々の空間周波数データの振幅情報と位相情報との少なくとも一方に基づいて、前記第１の磁気共鳴画像と前記第２の磁気共鳴画像との間で位置と大きさの少なくとも一方を整合することを特徴とする磁気共鳴映像装置。
【請求項２】
前記振幅情報と前記位相情報との少なくとも一方に関する前記第１の磁気共鳴画像と前記第２の磁気共鳴画像との間の相互相関係数に基づいて、前記整合のための補正値を求め、前記補正値に基づいて前記第１の磁気共鳴画像と前記第２の磁気共鳴画像との少なくとも一方を補正することを特徴とする請求項１記載の磁気共鳴映像装置。
【請求項３】
前記補正値は前記相互相関係数が最大となるように求められることを特徴とする請求項２記載の磁気共鳴映像装置。
【請求項４】
前記第１の磁気共鳴画像と前記第２の磁気共鳴画像各々の空間周波数データを極座標表現に変換し、前記極座標表現のもとで、スケール及び回転に関するずれを整合するための補正値を求め、前記補正値に基づいて前記第１の磁気共鳴画像と前記第２の磁気共鳴画像との少なくとも一方を補正することを特徴とする請求項１記載の磁気共鳴映像装置。
【請求項５】
前記極座標表現に変換された前記空間周波数データの径方向の座標スケールを対数的に変換して、前記補正値を求めることを特徴とする請求項４記載の磁気共鳴映像装置。
【請求項６】
前記空間周波数データから直交座標表現のもとで、前記整合のための補正値を求めることを特徴とする請求項１記載の磁気共鳴映像装置。
【請求項７】
前記整合を必要回数繰り返し実行することを特徴とする請求項１記載の磁気共鳴映像装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、一様な静磁場中に置かれた被検体に高周波磁場と勾配磁場を所定のパルスシーケンスに従って印加し、被検体からの磁気共鳴信号を検出して被検体の内部情報を磁気共鳴画像として映像化する磁気共鳴映像装置に関する。
【０００２】
【従来の技術】
磁気共鳴映像法はよく知られているように、固有の磁気モーメントを持つ核スピンの集団が一様な静磁場中に置かれたときに、特定の周波数で回転する高周波磁場のエネルギーを共鳴的に吸収する現象を利用して、物質の化学的および物理的な微視的情報を映像化する手法である。
【０００３】
この磁気共鳴映像法では、核スピンの縦緩和時間Ｔ1 を強調したコントラストの画像（以下、Ｔ1 画像という）、核スピンの横緩和時間Ｔ2 を強調したコントラストの画像（以下、Ｔ2 画像という）、核スピンの密度分布を強調したコントラストの画像（以下、密度画像という）、核スピンの横緩和時間Ｔ2 とボクセル内での微視的な磁場不均一性による核スピンの急激な位相変化を反映した実際的な横緩和時間Ｔ2 * を強調したコントラストの画像（以下、Ｔ2 * 画像という）といった種々のコントラストの画像を得る事ができる。
【０００４】
Magnetic Resonance in Medicine 23,37-45(1992) に記載されているように、生体組織内の局所的な血流量や血流速が変化すると磁気共鳴映像装置のある種の撮像法では生体組織の緩和時間（例えばＴ1 等）が見かけ上、変化したように観測され、画像コントラストが変化する。これらの撮像法や造影剤を利用することにより、生体内のperfusion やdiffusion に関する情報を得ることができる。perfusion やdiffusion の画像化法では、時間的な変化を調べるため複数の画像の差分を求めたり、あるいは微少な信号変化を検出するため加算平均を求めることが多いが、しばしば生体の動きの影響等で画像間での位置ずれが問題となる。
【０００５】
一方、Magnetic Resonance in Medicine 14,68-78(1990) に記載されているように、生体内の血中ヘモグロビンは動脈血中に多量に含まれる酸化ヘモグロビンは反磁性を示し、静脈血中に多く含まれる還元ヘモグロビンは常磁性を示す事が知られている。そしてMagnetic Resonance in Medicine 24,375-383(1992) に記載されているように、反磁性物質である酸化ヘモグロビンは局所的な磁場をあまり乱さないが（生体組織との帯磁率差０．０２ｐｐｍ）、常磁性物質である還元ヘモグロビンは周辺組織との帯磁率差が大きく（生体組織との帯磁率差０．１５ｐｐｍ）局所的に磁場を乱すためＴ2 * が短縮される。
【０００６】
上記性質を利用する事により、例えば光刺激に伴う脳の皮質部にある視覚領域の活性など、生体組織内の細胞活性等の生理機能に起因した酸素濃度の変化や血流の変化を画像化できる事が、Proc.Natl.Acad.Sci.USA 89,5675-5679(1992) 等に記載されている。これらの画像化に際して用いられる撮像法は、一般にグラジェントエコー法やエコープラナー法と呼ばれるパルスシーケンスである。
【０００７】
しかし、これらの画像化法によって得られる生体内の生理機能に伴って生ずる信号変化（画像コントラスト変化）は非常に微少である。そのため、この微少な信号変化を検出する方法として、生理機能現象が起こる前後の画像の差分をとる方法や統計的処理を行う方法が従来より用いられている。統計的データ処理法としては、Magnetic Resonance Imaging11,451-459(1993)に記載されているpaired-t検定法を用いた方法等がある。差分法を用いる場合には、ＳＮ比の高い画像を得る必要があり、また統計的処理を行う場合には複数の画像が必要となるため、撮影時間が長くなる。そのために、生体の動きの影響を受けやすい。
【０００８】
また、他の生体の動き補正法としては、画像自身の大きさや位置を補正する、いわゆるパターンマッチングを行う方法がある。パターンマッチングでよく知られている方法としては、画像間の相互相関を利用する方法がある。位相の相互相関を利用して磁気共鳴画像の補正を行った例はIEEE Tras.Pattern Anal. Machine Intel,16(12),1156-1168(1994)に記載されている。
【０００９】
また、静磁場分布が不均一な場合には画像歪を生じる事は良く知られているが、特に上記生体の細胞活性などの生理機能現象の検出に用いるＴ2 * コントラストの画像撮像法においては上記画像歪が顕著である。このような画像歪みを、アフィン変換などの方法を用いて補正する方法については、特願平５−２２７５９号に記載されている。
【００１０】
【発明が解決しようとする課題】
生体内のperfusion 、diffusion 、細胞活性等の情報を画像化するために複数の画像間でデータの処理を実施する場合に、生体の動きの影響等による画像の位置ずれや大きさの変化のために生じる画像の不一致があると、良好にperfusion 、diffusion 、細胞活性等の情報を画像化することが難しい。
【００１１】
実際に脳が心拍に同期して大きさや位置が変化する事が、Radiology,185,645-651(1992) にも記載されているように、よく知られている。
このように、従来の方法では呼吸や心拍等に伴う体動の影響等により、生体内の種々の生理機能情報を正確に検出する事ができないという問題がある。
【００１２】
本発明は上記問題点を解決し、ＭＲ画像間の位置や大きさを整合して、被検体内の生理機能情報を高精度で得ることができる磁気共鳴映像装置を提供することを目的とする。
【００１３】
【課題を解決するための手段】
本発明は、一様な静磁場中に置かれた被検体に高周波磁場と勾配磁場を所定のパルスシーケンスに従って印加し、上記被検体からの磁気共鳴信号を検出して上記被検体の内部情報を第１の磁気共鳴画像と第２の磁気共鳴画像として映像化する磁気共鳴映像装置において、上記第１の磁気共鳴画像の上記第２の磁気共鳴画像各々の空間周波数データの振幅情報と位相情報との少なくとも一方に基づいて、上記第１の磁気共鳴画像と上記第２の磁気共鳴画像との間で位置と大きさの少なくとも一方を整合することを特徴とする。
【００１４】
本発明により、第１と第２の磁気共鳴画像は位置と大きさの少なくとも一方が整合される。したがって、両画像間の差分等により与えられる生理機能情報が、被検体の動き等の影響を受けることなく高精度で得られる。
【００１５】
【発明の実施の形態】
以下、図面を参照して本発明の一実施の形態について説明する。
図１は、本発明の一実施の形態に係る磁気共鳴診断装置の構成を示すブロック図である。同図において、静磁場磁石１、磁場均一性調整コイル３および勾配磁場生成コイル５はそれぞれ励磁用電源２、磁場均一性調整コイル用電源４および勾配磁場生成コイル用電源６にて駆動される。これらにより被検体７には一様な静磁場とそれと同一方向で互いに直交する３方向に線形傾斜磁場分布を持つ勾配磁場が印加される。被検体７に高周波磁場を印加するために、送信部１０から高周波信号がプローブ９に送られる。ここでプローブ９は送受信両用でも、あるいは送受信別々に設けてもよい。プローブ９で受信された磁気共鳴信号は受信部１１で直交位相検波された後、データ収集部１３に転送されＡ／Ｄ変換後、電子計算機１４に送られる。以上、励磁用電源２、磁場均一性調整コイル用電源４、勾配磁場生成コイル用電源６、送信部１０、受信部１１、データ収集部１３はすべてシステムコントローラ１２によって制御されている。システムコントローラ１２は電子計算機１４を介してコンソール１５により制御される。
【００１６】
電子計算機１４ではデータ収集部１３から送られた磁気共鳴信号に基づいてＭＲ画像（磁気共鳴画像）を再構成し、ＭＲ画像データを得る。得られたＭＲ画像は画像表示装置１６に表示される。電子計算機１４および寝台８はコンソール１５により制御される。
【００１７】
本実施の形態における被検体７内のスライス面内のＭＲ画像データを収集するためのパルスシーケンス、及び被検体７に光や音等の刺激を与える刺激装置１７はシステムコントローラ１２によって制御される。
【００１８】
図２と図３はそれぞれ本実施の形態に係る被検体７内の生理機能画像化のためのフィールドエコー法によるパルスシーケンス、エコープラナー法によるパルスシーケンスである。図中のＲＦは高周波磁場であり、Ｇｓ，Ｇｒ，Ｇｅはそれぞれスライス用、読み出し用、位相エンコード用の勾配磁場であり、ＳＩＧ／ＡＤＣは磁気共鳴信号とデータ収集（サンプリング）のタイミングをそれぞれ示す。Ｇｓは被検体７内の所望の領域を励起するためのスライス用勾配磁場、Ｇｒは磁気共鳴信号の周波数情報に位置情報をエンコードするための読み出し用勾配磁場、Ｇｅは磁気共鳴信号の位相情報に位置情報をエンコードするための位相エンコード用勾配磁場である。
【００１９】
図２においては、はじめに高周波磁場パルスＲＦとスライス用勾配磁場Ｇｓとを印加して、所望の領域の磁化スピンを選択的に励起し、自由誘導減衰ＮＭＲ信号を発生させる。続いて、位相エンコード用勾配磁場Ｇｅ、読み出し用勾配磁場Ｇｒを順次印加し、読み出し用勾配磁場Ｇｒの印加中に発生するエコー信号ｅｃｈｏを収集する。そして、同様のパルスシーケンスを位相エンコード用勾配磁場Ｇｅの印加量を順次変えながら繰り返し時間ＴＲで繰り返し実行する。生理機能を画像化するための典型的な条件は、繰り返し時間ＴＲが５０〜１００ミリ秒、エコー時間ＴＥ（高周波磁場パルスＲＦの中心からデータを配列した際に中心となるデータまでの時間間隔）が３０〜７０ミリ秒である。また、高周波磁場パルスによるスピンの励起角（フリップ角）は１０〜４０°である。
【００２０】
図３においては、はじめに高周波磁場パルスＲＦとスライス用勾配磁場Ｇｓを印加して所望の領域の磁化スピンを選択的に励起し、自由誘導減衰ＮＭＲ信号を発生させる。続いて読み出し用勾配磁場Ｇｒを正負交互にスイッチング（極性反転）しながら複数のエコー信号を順次発生させ、その各々のエコー信号毎に位相エンコード用勾配磁場Ｇｅを印加する。そして、各読み出し用勾配磁場Ｇｒの中心付近で発生する複数のエコー信号ｅｃｈｏをそれぞれ収集する。この場合には、１回のスピンの励起で１画像分のデータを得ることができる。生理機能を画像化のための典型的な条件は、エコー時間ＴＥ（高周波磁場パルスＲＦの中心からデータを２次元配列した際に原点となるデータまでの時間間隔）が５０〜７０ミリ秒である。
【００２１】
図２あるいは図３のパルスシーケンスが繰り返され、適当な前処理をした後に、複素フーリエ変換することにより、時系列に連続した一連のＭＲ画像の群が収集される。このようにして得られるＭＲ画像はＴ2 * コントラストの画像であり、上述したように刺激や負荷に反応して脳細胞の特定部位が活性化され組織内酸素濃度や局所血流が変化することにより生じる活性化部位とその近傍での帯磁率変化に伴うＴ2 * コントラスト変化を捕えることができる。また、上記パルスシーケンスの条件によっては、刺激や負荷に反応した血流変化自体に伴うコントラスト変化を捕えることもできる。
【００２２】
実際に、前述したようなパルスシーケンスを用いて生理機能情報を得るためには、上記刺激や負荷を与えた場合と与えない場合との複数フレームのＭＲ画像、つまり刺激付与に対して時間的に前後して収集された少なくとも２フレームのＭＲ画像に対して、フレーム間で加算平均、差分、統計処理等の画像処理を行う必要がある。したがって、高精度で生理機能情報を得るためには、これらＭＲ画像間の位置やスケール（大きさ）の整合がとれている必要がある。なお、“位置”とは、勾配磁場方向に沿った空間的な位置のことをいう。本発明に係わる一連の複数のＭＲ画像間の位置や大きさを合わせる（整合する）ために必要とされる画像処理手段としては加算平均、差分、統計処理等の画像処理の前段階で電気計算機１３で実行される。
【００２３】
以下、複数のＭＲ画像間の位置や大きさを整合する画像処理手段について説明する。この画像処理手段は電子計算機１３に含まれる。
本実施の形態による画像処理手段の処理の流れを図４に示す。始めに、位置やスケールの基準とされるべき参照画像データ（第１のＭＲ画像データ）と、位置や大きさが参照画像に整合されるべき被補正画像データ（第２のＭＲ画像データ）とが当該画像処理手段に入力される。参照画像としては、オペレータによりコンソール１５を介して一連のＭＲ画像群の中から選択されるようにしてもよいし、一連のＭＲ画像群の中から最初に撮像されたＭＲ画像を参照画像としてもよい。ＭＲ画像としては、２次元の画像データでも、３次元の画像データでも良い。ここでは、２次元の画像データとして説明する。本実施の形態では、参照画像データはＲ（ｘ，ｙ）として表現され、被補正画像データはＦ（ｘ，ｙ）として表現されるものとする。
【００２４】
参照画像データと被補正画像データは、必要に応じて補間処理、フィルタ処理、マスク処理等の前処理に供せられる（Ｓ1 ，Ｓ11）。そして、ＦＦＴ（高速フーリエ変換）にかけられる（Ｓ2 ，Ｓ12）。これにより、参照画像データと被補正画像データはそれぞれ、空間周波数データｒ（ｋｘ，ｋｙ）、ｆ（ｋｘ，ｋｙ）に変換される。
【００２５】
次に、スケール・回転補正処理Ｓ100 が行われる。スケール・回転補正処理は、画像の回転の取扱いが容易になるように、まず、上記空間周波数データｒ（ｋｘ，ｋｙ）、ｆ（ｋｘ，ｋｙ）各々の絶対値が極座標形式のデータ｜Ｒ′（ｒ，θ）｜、｜Ｆ′（ｒ，θ）｜に変換される（Ｓ3 ，Ｓ13）。
【００２６】
更に、画像のスケールの取扱いが容易になるように、極座標形式のデータ｜Ｒ′（ｒ，θ）｜、｜Ｆ′（ｒ，θ）｜それぞれは、径方向の座標スケールが対数的に変換される（Ｓ4 ，Ｓ14）。これにより｜Ｒ″（ｌｎ（ｒ），θ）｜、｜Ｆ″（ｌｎ（ｒ），θ）｜が得られる。
【００２７】
そして、Ｈａｎｎｉｎｇフィルタ等の適当なフィルタ処理が行われる（Ｓ5 ，Ｓ15）。フィルタ処理された参照画像データと非補正画像データは（１）式のように与えられる。なお、原画像に対してｌｎ（ｒ）は画像のスケールを示し、θは回転を示す。
【００２８】
【数１】
0003672976.tif 000002
【００２９】
次に、フィルタ処理された両画像データ｜Ｒ″（ｌｎ（ｒ），θ）｜、｜Ｆ″（ｌｎ（ｒ），θ）｜との間で、位相の相互相関が算出される（Ｓ21）。まず、（１）式のフィルタ処理された参照画像データと非補正画像データをそれぞれＦＦＴにかけたものから、（２）式で定義される関係を用いて、係数ｐｃ１（ｌｎ（ｋｒ），θ）を求める。但し、＊は、複素共役を示す。
【００３０】
【数２】
0003672976.tif 000003
【００３１】
さらに、係数ｐｃ１（ｌｎ（ｋｒ），θ）をＩＦＦＴ（高速逆フーリエ変換）して、両画像間の位相の相互相関係数ＰＣ１（ｌｎ（ｒ），θ）が得られる。
この位相相関係数ＰＣ１（ｌｎ（ｒ），θ）が最大を示すような、ｌｎ（ｒ）とθをそれぞれスケールと回転の補正値として求める（Ｓ22）。この補正値をそれぞれλpc、θpcとする。但し、λpcは、周波数空間におけるスケールファクターであり、実空間上で画像の補正処理を行うために実空間におけるスケールファクターに変換する必要があり、具体的には１／λpcを実空間におけるスケールファクターとすればよい。
そして、被補正画像に対して（３）式の関係を用いてスケールと回転に関する補正を行う（Ｓ23）。
【００３２】
【数３】
0003672976.tif 000004
【００３３】
次に、上記スケール・回転補正を受けた被補正画像Ｆ （ｘ，ｙ）と参照画像Ｒ（ｘ，ｙ）に対して位置補正処理Ｓ101 を行う。まず、上記スケール・回転補正画像データＦ （ｘ，ｙ）をＦＦＴして空間周波数データｆ （ｋｘ，ｋｙ）に変換する。そして、ｆ （ｋｘ，ｋｙ）と、参照画像の空間周波数データｒ（ｋｘ，ｋｙ）に対して（４）式で定義される関係式にしたがって両画像間の位相相互相関ｐｃ２（ｋｘ，ｋｙ）を算出する（Ｓ31）。
【００３４】
【数４】
0003672976.tif 000005
【００３５】
次に、（４）式で得られた位相相関係数ｐｃ２（ｋｘ，ｋｙ）を、ＩＦＦＴ（高速逆フーリエ変換）してｐｃ２（ｘ，ｙ）を算出し、このｐｃ２（ｘ，ｙ）を最大とするような、ｘとｙをそれぞれ位置の補正値として求める（Ｓ32）。この補正値をそれぞれｘpc、ｙpcとする。
【００３６】
そして、被補正画像に対して（５）式の関係を用いて位置に関する補正を行い、場合によっては補間処理やフィルタ処理を行って、スケール、回転、位置が参照画像に対して整合された補正画像Ｆout （ｘ，ｙ）が得られる（Ｓ33）。
【００３７】
【数５】
0003672976.tif 000006
【００３８】
上記説明では、磁気共鳴画像の空間周波数データの位相のみの相互相関を利用したが、一般的な振幅と位相の両方の情報を利用した相互相関、あるいは振幅のみの相互相関であってもよい。両方の情報を利用した相互相関係数は、（２）式または（４）式の分母を１として分子だけを使って求め、また振幅のみの相互相関係数は、（２）式または（４）式の分母だけを使って求める。
【００３９】
また、図５に示すように、１回の整合処理Ｓ100 ，Ｓ101 で十分な補正画像が得られない場合には、上記１回目の補正画像Ｆout １（ｘ，ｙ）を被補正画像として、図４のスケール・回転補正処理Ｓ100 と位置補正処理Ｓ101 とを所望の補正画像が得られるまで繰り返し実行するようにしてもよい。
【００４０】
本発明は、上記以外にも主旨を逸脱しない範囲で種々変形して実施する事が可能である。例えば、上記説明では入力データを画像データとしているが、上記画像の原収集データ（ｋ空間に展開されるデータ）を利用する事も可能である。この場合には、一部ＦＦＴ等の処理を省略することができる。また、上記説明では、２次元の画像データに関する補正を示したが、３次元データに対しても同様の補正は可能である。３次元データの場合には、回転及び位置ずれ補正のパラメータがそれぞれ１つずつ増える。
【００４１】
また、本発明は基本的には上記一連の処理を実施するものであるが、上記処理の一部のみを実施するなど、種々変形して適用することも可能である。
更に、本発明は上記実施の形態以外にも、ｐｅｒｆｕｓｉｏｎやｄｉｆｆｕｓｉｏｎ画像をはじめ、種々の磁気共鳴画像に適用することができる。
【００４２】
【発明の効果】
本発明により、第１と第２の磁気共鳴画像間の位置と大きさの少なくとも一方が整合される。したがって、両画像間の差分等により与えられる生理機能情報が、被検体の動き等の影響を受けることなく高精度で得られる。
【図面の簡単な説明】
【図１】本発明の一実施の形態に係る磁気共鳴映像装置の構成を示すブロック図。
【図２】磁気共鳴画像を得るためのフィールドエコー法のパルスシーケンスを示す図。
【図３】磁気共鳴画像を得るためのエコープラナー法のパルスシーケンスを示す図。
【図４】スケール、回転、位置の整合処理の手順を示す図。
【図５】スケール、回転、位置の整合処理を繰り返す手順を示す図。
【符号の説明】
１…静磁場磁石、
２…励磁用電源、
３…静磁場均一性調整コイル、
４…静磁場均一性調整コイル用電源、
５…勾配磁場生成コイル、
６…勾配磁場生成コイル用電源、
７…被検体、
８…寝台、
９…プローブ、
１０…送信部、
１１…受信部、
１２…システムコントローラ、
１３…データ収集部、
１４…電子計算機、
１５…コンソール、
１６…画像ディスプレイ、
１７…刺激装置。
【図１】
0003672976.tif 000007
【図２】
0003672976.tif 000008
【図３】
0003672976.tif 000009
【図４】
0003672976.tif 000010
【図５】
0003672976.tif 000011