(57)【要約】【課題】磁性微粒子の３次元分布を画像化する場合の撮影時間を短縮する。【解決手段】高周波コイルが、ゼロ磁場領域内の位置に応じて高周波磁場の方向を変化させる。また、複数の検出コイルが、磁性微粒子が磁気飽和することにより生じる磁場の変化を信号値として検出する。そして、制御装置１５の磁化検出部１５ｂが、高周波コイルにより印加される高周波磁場の方向と、ゼロ磁場領域の位置と、ゼロ磁場領域周辺の磁場分布とから磁気モーメントの方向を特定し、特定した磁気モーメントの方向および検出コイルにより検出された信号値に基づいて磁気モーメントの大きさを求めることによって、磁化の分布を計測する。【選択図】        図１
2010088683.tif 000002
【特許請求の範囲】
【請求項１】
被検体が置かれる撮影空間に無磁場領域を形成する無磁場領域形成手段と、
前記撮影空間に高周波磁場を印加する高周波磁場印加手段と、
前記高周波磁場印加手段により印加される高周波磁場によって前記無磁場領域内の磁性微粒子が磁気飽和することにより生じる磁化の分布を計測する磁化分布計測手段と、
前記磁化分布計測手段によって計測された磁化の分布を画像化する画像生成手段と、
を備え、
前記磁化分布計測手段は、前記無磁場領域内または前記無磁場領域周辺に位置する複数の点にそれぞれ生じる磁気モーメントの大きさを求めることによって前記磁化分布を計測することを特徴とする磁性微粒子イメージング装置。
【請求項２】
前記高周波磁場印加手段は、前記無磁場領域内の位置に応じて、前記高周波磁場の方向を変化させ、
前記磁化分布計測手段は、
前記磁性微粒子が磁気飽和することにより生じる磁場の変化を信号値として検出する複数の磁場検出手段と、
前記高周波磁場印加手段により印加される高周波磁場の方向と前記無磁場領域の位置と前記無磁場領域周辺の磁場分布とから前記磁気モーメントの方向を特定し、特定した磁気モーメントの方向および前記磁場検出手段により検出された信号値に基づいて前記磁気モーメントの大きさを求めることによって前記磁化の分布を計測する磁化検出手段と、
を有することを特徴とする請求項１に記載の磁性微粒子イメージング装置。
【請求項３】
前記高周波磁場印加手段および前記磁場検出手段はそれぞれコイルであり、
前記磁場検出手段は、前記高周波磁場印加手段によって印加される磁場に対してコイルの軸が直交するように配置されることを特徴とする請求項２に記載の磁性微粒子イメージング装置。
【請求項４】
前記無磁場領域形成手段は、線状の無磁場領域を形成し、
前記磁化分布計測手段は、前記無磁場領域の軸方向に沿った磁化の分布を計測することを特徴とする請求項１、２または３に記載の磁性微粒子イメージング装置。
【請求項５】
前記磁化分布計測手段は、前記無磁場領域の周方向に沿った磁化の分布を計測することを特徴とする請求項１〜４のいずれか一つに記載の磁性微粒子イメージング装置。
【請求項６】
前記高周波磁場印加手段は、前記撮影空間内で磁場が回転するように前記高周波磁場を印加することを特徴とする請求項１〜５のいずれか一つに記載の磁性微粒子イメージング装置。
【請求項７】
前記無磁場領域形成手段によって形成された無磁場領域で前記被検体内を走査する無磁場領域走査手段をさらに備えたことを特徴とする請求項１〜６のいずれか一つに記載の磁性微粒子イメージング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、被検体に注入された磁性微粒子の分布を画像化する磁性微粒子イメージング装置に関し、特に、磁性微粒子の３次元分布を画像化する場合の撮影時間を短縮することができる磁性微粒子イメージング装置に関する。
【背景技術】
【０００２】
近年、被検体内に造影剤として注入された磁性微粒子（例えば、超常磁性酸化鉄など）の分布を画像化する方法が考案されている（例えば、非特許文献１、特許文献１および２を参照）。この方法は、「磁性微粒子イメージング（ＭＰＩ：Magnetic Particle Imaging）」と呼ばれている。
【０００３】
図１７は、従来の磁性微粒子イメージングの原理を説明するための図である。同図に示すように、磁性微粒子イメージングでは、ゼロ磁場領域形成コイル（Selection coils）と、高周波コイル（drive field coils）と、検出コイル（Recording coils）とが用いられる。
【０００４】
ゼロ磁場領域形成コイルは、互いに反対向きの磁場を生成することによって、撮影空間の中央部にゼロ磁場領域（無磁場領域：Field-free point）を発生させる。高周波コイルは、例えば、２５kHz程度の高周波磁場を撮影空間に印加する。ここで、ゼロ磁場領域の外では、高周波磁場が印加されても磁性微粒子の磁束密度が常に飽和しており、磁束密度は変化しない。他方、ゼロ磁場領域内では、静磁場による飽和が起きておらず、高周波磁場が印加されると磁気飽和が起きる。したがって、ゼロ磁場領域内では、磁気飽和をともなう磁束密度が生じる。
【０００５】
ゼロ磁場領域内での磁束密度の変化は、検出コイルを鎖交する磁束の変化を生じる。この磁束の変化の大きさはゼロ磁場領域内の磁性微粒子の量に依存しており、検出コイルを鎖交する磁束の変化は検出コイルに誘起される起電力（電圧）の変化として現れる。すなわち、検出コイルに誘起される電圧はゼロ磁場領域内の磁性微粒子の量に応じて変化することになる。
【０００６】
このような原理を利用し、撮影空間内でゼロ磁場領域を少しずつ移動しながら、検出コイルに誘起される電圧の計測を行えば、磁性微粒子の分布が反映された画像を作成することができる。なお、かかる磁性微粒子イメージングを実現するための装置は「磁性微粒子イメージング装置」と呼ばれる。
【０００７】
【非特許文献１】B. Gleich, J. Borgert, J. Weizenecker, "Magnetic Particle Imaging (MPI), Philips Medic Mundi vol.50 no.1, 2006/5 [online], May 23 2007, retrieved from the Internet: <URL:http://www.medical.philips.com/main/news/assets/docs/medicamundi/mm_vol50_no1/12_Gleich.pdf>
【特許文献１】特開２００３−１９９７６７号公報
【特許文献２】特表２００６−５２３４９２号公報
【発明の開示】
【発明が解決しようとする課題】
【０００８】
ところで、上述した従来の磁性微粒子イメージングでは、磁性粒子の３次元分布を画像化する場合、被検体内でゼロ磁場領域を３軸方向に走査する必要があるが、ゼロ磁場領域が点状の領域であるため撮影時間が長時間となる。例えば、一点を１０kHzの高周波磁場の１００周期分で計測する場合、一点の計測にかかる所要時間は１０msとなる。この場合、３mmの分解能で一辺が２００mmの立方体の範囲を撮影するためには２９６０００点の計測が必要となるので、合計で約５０分の撮影時間が必要となる。
【０００９】
そのため、臨床への応用を考えた場合には、いかにして磁性微粒子の３次元分布を画像化する場合の撮影時間を短縮するかが大きな技術的課題となっている。
【００１０】
この発明は、上述した従来技術による課題を解決するためになされたものであり、磁性微粒子の３次元分布を画像化する場合の撮影時間を短縮することができる磁性微粒子イメージング装置を提供することを目的とする。
【課題を解決するための手段】
【００１１】
上述した課題を解決し、目的を達成するため、請求項１記載の本発明は、被検体が置かれる撮影空間に無磁場領域を形成する無磁場領域形成手段と、前記撮影空間に高周波磁場を印加する高周波磁場印加手段と、前記高周波磁場印加手段により印加される高周波磁場によって前記無磁場領域内の磁性微粒子が磁気飽和することにより生じる磁化の分布を計測する磁化分布計測手段と、前記磁化分布計測手段によって計測された磁化の分布を画像化する画像生成手段と、を備え、前記磁化分布計測手段は、前記無磁場領域内または前記無磁場領域周辺に位置する複数の点にそれぞれ生じる磁気モーメントの大きさを求めることによって前記磁化の分布を計測することを特徴とする。
【発明の効果】
【００１２】
請求項１記載の本発明によれば、ゼロ磁場領域を固定した状態でも同時に複数の点における磁化を計測することができるようになり、磁性微粒子の３次元分布を画像化する場合の撮影時間を短縮することができるという効果を奏する。
【発明を実施するための最良の形態】
【００１３】
以下に添付図面を参照して、この発明に係る磁性微粒子イメージング装置の好適な実施例を詳細に説明する。
【実施例】
【００１４】
本実施例に係る磁性微粒子イメージング装置は、被検体が置かれた撮影空間にゼロ磁場領域を形成し、形成したゼロ磁場領域に高周波磁場を印加する。そして、ゼロ磁場領域内の磁性微粒子が磁気飽和することによって生じる磁化の分布を計測し、計測した磁化の分布すなわち磁性微粒子の分布を画像化する。ここで、従来の磁性微粒子イメージング装置は、点状に形成したゼロ磁場領域のただ一点における磁化の大きさを計測するものであり、ゼロ磁場領域で被検体内を走査することによって磁化の分布を計測していた。
【００１５】
これに対し、本実施例に係る磁性微粒子イメージング装置は、ゼロ磁場領域内またはゼロ磁場領域周辺に位置する複数の点にそれぞれ生じる磁気モーメントの大きさを求めることによって、磁化の分布を計測する。これにより、ゼロ磁場領域を固定した状態でも同時に複数の点における磁化を計測することができるようになり、磁性微粒子の３次元分布を画像化する場合の撮影時間を短縮することが可能になる。
【００１６】
以下、かかる磁性微粒子イメージングについて具体的に説明する。まず、本実施例に係る磁性微粒子イメージング装置の構成について説明する。図１および２は、本実施例に係る磁性微粒子イメージング装置１０の構成を示す図である。図１は、磁性微粒子イメージング装置の斜視図であり、図２は、磁性微粒子イメージング装置の正面図である。
【００１７】
図１に示すように、本実施例に係る磁性微粒子イメージング装置１０は、ゼロ磁場領域形成磁石１１ａおよび１１ｂと、ゼロ磁場スキャンコイル１２ａおよび１２ｂと、寝台１３と、架台１４と、制御装置１５とを有する。さらに、図２に示すように、磁性微粒子イメージング装置１０は、高周波コイル１６ａ〜１６ｆを有する。
【００１８】
なお、図１および２では図示を省略しているが、磁性微粒子イメージング装置１０は、さらに複数の検出コイルを備える。各検出コイルは、磁性微粒子が磁気飽和することにより生じる磁場の変化を信号値として検出する。かかる検出コイルについては、後に詳細に説明する。
【００１９】
ゼロ磁場領域形成磁石１１ａおよび１１ｂは、被検体Ｐが置かれる撮影空間にゼロ磁場領域（無磁場領域）を形成する。具体的には、ゼロ磁場領域形成磁石１１ａおよび１１ｂは、それぞれ永久磁石であり、図１に示すように、磁化の方向を逆にしてそれぞれ配置されている。なお、ここでは、ゼロ磁場領域形成磁石１１ａおよび１１ｂとして永久磁石を用いた場合について説明するが、永久磁石の代わりに電磁石を用いてもよい。
【００２０】
図３および４は、ゼロ磁場領域形成磁石１１ａおよび１１ｂによって形成されるゼロ磁場領域を説明するための図である。図３に示すように、磁化の方向を逆にしてゼロ磁場領域形成磁石１１ａおよび１１ｂを配置することによって、２つの磁石の中間点で静磁場が相殺され、ゼロ磁場領域が作られる。このゼロ磁場領域は、図４に示すように、頭−足方向に伸びる直線状の領域となる。
【００２１】
図１にもどって、ゼロ磁場スキャンコイル１２ａおよび１２ｂは、撮影空間内でゼロ磁場領域を左−右方向に移動させる。具体的には、ゼロ磁場スキャンコイル１２ａがＮ極、ゼロ磁場スキャンコイル１２ｂがＳ極となるように、それぞれのコイルに電流を流せば、ゼロ磁場領域は左方向（図１における右側）へ移動し、逆向きに電流を流せば右（図１における左側）へ移動する。
【００２２】
寝台１３は、被検体Ｐが載置され、撮影空間内で被検体Ｐを上−下方向、左−右方向、頭−足方向に移動する。ゼロ磁場スキャンコイル１２ａおよび１２ｂでは、ゼロ磁場領域を上下方向へ移動させることはできないが、寝台１３が上下方向へ移動することによって、等価的に、被検体内でゼロ磁場領域を上下方向へ移動させることができる。
【００２３】
架台１４は、磁性微粒子イメージング装置１０を構成する各部を支持する。
【００２４】
制御装置１５は、磁性微粒子イメージングの全体制御を行い、特に、コイル制御部１５ａと、磁化検出部１５ｂと、画像生成部１５ｃとを有する。
【００２５】
コイル制御部１５ａは、磁性微粒子イメージング装置１０が有する各種コイルを制御する。例えば、コイル制御部１５ａは、ゼロ磁場領域を移動するためにゼロ磁場スキャンコイル１２ａおよび１２ｂを制御したり、ゼロ磁場領域に高周波磁場を印加するために高周波コイル１６ａ〜１６ｆを制御したりする。
【００２６】
磁化検出部１５ｂは、高周波コイル１６ａ〜１６ｆにより印加される高周波磁場の方向と、ゼロ磁場領域の位置と、ゼロ磁場領域周辺の磁場分布とから磁気モーメントの方向を特定し、特定した磁気モーメントの方向および検出コイルにより検出された信号値に基づいて磁気モーメントの大きさを求めることによって磁化の分布を計測する。この磁化検出部１５ｂによって行われる磁化分布の計測については後に詳細に説明する。
【００２７】
画像生成部１５ｃは、磁化分布計測手段によって計測された磁化の２次元分布あるいは３次元分布を画像化する。
【００２８】
高周波コイル１６ａ〜１６ｆは、撮影空間に高周波磁場を印加する。具体的には、高周波コイル１６ａ〜１６ｆは、それぞれ電磁石であり、１０kHzから１００kHz程度の高周波磁場を発生させる。被検体Ｐに磁性微粒子を分布させたうえで、高周波コイル１６ａ〜１６ｆによって高周波磁場を印加することによって、ゼロ磁場領域に存在する磁性微粒子の量に応じて、歪を含んだ高周波磁場がゼロ磁場領域から発生する。
【００２９】
図５−１および５−２は、高周波コイル１６ａ〜１６ｆによる高周波磁場の印加を説明するための図である。具体的には、図５−１に示すように、高周波コイル１６ｅおよび１６ｆは上下方向の高周波磁場を印加し、図５−２に示すように、高周波コイル１６ａ〜１６ｄは左−右方向の高周波磁場を印加する。これらを組み合わせることによって、斜め方向を含む任意の向きの高周波磁場を印加することができる。
【００３０】
ここでは、高周波コイル１６ａ〜１６ｆは、上下方向の高周波磁場の位相と左−右方向の高周波磁場の位相とを９０度ずらすことによって、撮影空間中の磁場の向きを１０kHzから１００kHzの周波数（以下、「変調周波数」と呼ぶ）で回転させるものとする。
【００３１】
図６は、高周波コイルの配置を説明するための図である。同図に示すように、磁性微粒子イメージング装置１０には、高周波コイル１６ａ〜１６ｆと同様のコイルの組である高周波コイル１６ｇ〜１６ｌが配置されている。ここで、高周波コイル１６ａ〜１６ｆと高周波コイル１６ｇ〜１６ｌとは、それぞれ異なる向きの磁場を発生させるものとする。
【００３２】
次に、上記で説明した構成によってゼロ磁場領域の近傍に形成される磁場について説明する。図７は、ゼロ磁場領域の近傍で形成される磁場を説明するための図である。同図は、ゼロ磁場領域形成磁石１１ａおよび１１ｂによって発生する静磁場の磁場分布を示す模式図である。なお、ここでは説明の便宜上、永久磁石によって発生する静磁場に近似するものとして、静磁場コイル（電磁石）によって発生した静磁場を示している。
【００３３】
同図において、グラフの目盛りが交差する点から伸びる直線の向きは、各点における静磁場の向きを示しており、直線の長さは、各点における静磁場の強さを示している。また、同図の上下にある静磁場コイルの矢印は、各静磁場コイルに流れる電流の向きを示しており、同図の中央にある円は、ゼロ磁場領域を示している。
【００３４】
図８は、図７に示したゼロ磁場領域の近傍を拡大した拡大図である。同図の（ａ）〜（ｄ）において、大きい円はゼロ磁場領域を示している。また、小さい円（網掛けされた円）は、それぞれ、高周波コイル１６ａ〜１６ｆによって実線の矢印の向きに磁場が印加された瞬間に磁気飽和する領域を示している。また、白抜きの矢印は、外部から近似的に観測される磁化Ｍの磁気モーメントを示している。
【００３５】
ここで、同図に示すように、高周波コイル１６ａ〜１６ｆによって印加される磁場の向きに応じて、磁気飽和が生じる位置はわずかに変化し、その変化に同期して磁気モーメントの向きも変化する。これは、磁化Ｍによって形成される磁場の分布から磁気モーメントの向きを知ることができれば、撮影空間内で磁気飽和が生じている位置を知ることができることを示しており、特に、ゼロ磁場領域内ではより詳細に磁気飽和が生じている位置を知ることができることを示している。
【００３６】
なお、前述した図５−１および５−２を用いた説明で、高周波コイル１６ａ〜１６ｆの組と高周波コイル１６ｇ〜１６ｌの組が、それぞれ異なる向きの磁場を発生し、変調周波数で磁場を回転させることを述べた。
【００３７】
図９は、ある瞬間において高周波コイル１６ａ〜１６ｌが作る磁場の様子を示す模式図である。同図は、図６に示した高周波コイル１６ａ〜１６ｌと位置関係が対応しており、手前側に高周波コイル１６ａ〜１６ｆの組が、奥側に高周波コイル１６ｇ〜１６ｌの組がそれぞれ配置されているものとする。
【００３８】
同図に示す瞬間では、手前側の高周波コイル１６ａ〜１６ｆは横向きの磁場を発生し、奥側の高周波コイル１６ｇ〜１６ｌは下向きの磁場を発生している。実際の磁場はこれらの合成なので、直線状のゼロ磁場領域に沿って磁場の向きが変化している。
【００３９】
図１０は、図９に示した瞬間において計測される磁気モーメントを示す模式図である。前述したように、高周波磁場の向きによって磁気モーメントの向きも変化するので、図９に示した瞬間では、図１０に示すように、磁気モーメントは直線状のゼロ磁場領域に沿って配置される。なお、同図において、左右に配列された小円は、それぞれ、各向きの磁気モーメントが作る磁場を検出し、分解するための検出コイルを示している。
【００４０】
従来の方法では、ゼロ磁場領域をある位置に置いたとき、計測できる磁化の量（＝磁性粒子の分布に比例）は一点のみであった。本実施例に係る磁性微粒子イメージング装置１０では、以下に説明するように、ゼロ磁場領域を固定した状態でも、ある範囲の磁化の分布を測定することができる。
【００４１】
次に、本実施例における磁化分布の計測方法（分解方法）について説明する。ここでは、ゼロ磁場領域の周方向に沿った磁化の分布を計測する方法と、ゼロ磁場領域の軸方向に沿った磁化の分布を計測する方法とをそれぞれ説明する。
【００４２】
まず、ゼロ磁場領域の周方向に沿った磁化の分布を計測する方法について説明する。図１１−１および１１−２は、ゼロ磁場領域の周方向に沿った磁化の分布を計測する方法を説明するための図であり、それぞれ、検出コイル１７の配置の一例を示している。
【００４３】
本実施例に係る磁性微粒子イメージング装置１０では、全ての検出コイル１７は、高周波コイル１６ａ〜１６ｆが作る磁場に対してコイルの軸が直交するように配置される。例えば、図１１−１および１１−２に示すように、検出コイル１７は、高周波コイル１６ｅおよび１６ｆが作る磁場および高周波コイル１６ａ〜１６ｄが作る磁場の両方に対してコイルの軸が直交するように配置されている。
【００４４】
このように検出コイル１７を配置することによって、どのような向きの磁場を高周波コイル１６ａ〜１６ｄが発生させていても、検出コイル１７が、高周波コイル１６ａ〜１６ｄによって発生する磁場を直接計測することはない。つまり、前述したように高周波コイル１６ａ〜１６ｄは回転磁場を発生させるが、検出コイル１７は、この回転磁場に対して感度を持たない。
【００４５】
他方、検出コイル１７は、図１１−１および１１−２に示すように、高周波コイル１６ａ〜１６ｄによってゼロ磁場領域の近傍で発生する磁化Ｍ（磁気モーメント）に対しては感度を持つよう配置されている。磁気モーメントの向きは、高周波コイル１６ａ〜１６ｄが発生させる回転磁場に伴って回転するが、この回転中、複数の検出コイル１７のうち、少なくとも一つは、磁気モーメントが作る磁場を計測することになる。
【００４６】
すなわち、磁気モーメントの向きに応じて、複数の検出コイル１７によって計測される信号値のパターンが異なる。この信号値のパターンは、ｂｉ（θ）と表すことができる。ここで、θは、高周波コイル１６ａ〜１６ｄによって発生する磁場の向きであり、ｉは、検出コイル１７のインデックスである。実際の検出信号ｓｉ（θ）とｂｉ（θ）の内積は、磁化Ｍの大きさに比例する。
【００４７】
また、磁化Ｍが存在する領域の位置は、ゼロ磁場領域の中心からθに関連する方向に少しだけずれた位置となる（図８を参照）。この位置は、高周波コイル１６ａ〜１６ｄの配置などによって異なる。例えば、図１１−１および１１−２の例では、磁化Ｍが存在する領域の位置は、ゼロ磁場領域の中心から約（９０度−θ）の方向にずれた位置となる。
【００４８】
そこで、磁化Ｍの方向（磁気モーメントの方向）をφ＝π／２−θとすると、ゼロ磁場領域の中心からφの方向にずれた位置における磁気モーメントの大きさａ（φ）は、以下に示す式（１）で求めることができる。これにより、ゼロ磁場領域の周方向に沿った磁化の分布が得られる。ここで、αは適当な比例係数である。
【００４９】
ａ（φ）＝αΣｉ｛ｓｉ（π／２−φ）ｂｉ（π／２−φ）｝      ・・・（１）
【００５０】
なお、ここで説明したゼロ磁場領域の周方向に沿った磁化分布の計測は、前述した制御装置１５の磁化検出部１５ｂによって行われるものとする。
【００５１】
次に、ゼロ磁場領域の軸方向に沿った磁化の分布を計測する方法について説明する。図１２は、ゼロ磁場領域の軸方向に沿った磁化の分布を計測する方法を説明するための図である。同図は、図１０に示した瞬間で検出される検出信号および磁気モーメントを示している。同図に示すように、例えば、Ｎ個の磁気モーメントの大きさＭ１〜ＭＮの総和が、Ｍ個の検出コイルによって検出信号ｓ１〜ｓＭとして計測される。
【００５２】
ここで、磁気モーメントが存在する領域をＮ個に離散化すれば、各検出コイルによって計測される検出信号ｓｉは、以下に示す式（２）で表すことができる。
【００５３】
【数１】
2010088683.tif 000003
【００５４】
さらに、線形逆問題理論を用いることによって、以下に示す式（３）のように、Ｎ個の磁気モーメントの大きさＭ１〜ＭＮ（磁化の分布）を求めることができる。これにより、ゼロ磁場領域の軸方向に沿った磁化の分布が得られる。ここで、Ａ＋は一般逆行列を表している。
【００５５】
【数２】
2010088683.tif 000004
【００５６】
上記の方程式の特異性はそれほど大きくないため、安定的に解くことは比較的容易であり、線形逆問題の解法としては公知の方法を用いればよい。例えば、ムーア・ペンローズ一般逆行列などを用いることができるほか、公知の正則化法を適用することもできる。
【００５７】
また、式（１）における係数行列Ａのランクは、分離できる磁気モーメントの数に関係する。ゼロ磁場領域内の全点に同じ方向の磁場を印加する場合には（磁気モーメントの向きが同じ場合）、各点における磁気モーメントが作る磁場分布の違いが小さく、結果として、行列Ａのランクはあまり大きくならない。これは、分解できる磁気モーメントの数が少ないことを意味する。
【００５８】
本実施例では、各点の磁気モーメントの方向が異なるように磁場を印加するので、行列Ａのランクが増加する。これにより、分解できる磁気モーメントの数を大きくすることができ、結果として、撮影する磁化分布の分解能を高めることが可能になる。
【００５９】
なお、ここで説明したゼロ磁場領域の軸方向に沿った磁化分布の計測は、前述した制御装置１５の磁化検出部１５ｂによって行われるものとする。
【００６０】
上述してきたように、本実施例では、高周波コイル１６ａ〜１６ｆが、ゼロ磁場領域内の位置に応じて高周波磁場の方向を変化させる。また、複数の検出コイル１７が、磁性微粒子が磁気飽和することにより生じる磁場の変化を信号値として検出する。そして、制御装置１５の磁化検出部１５ｂが、高周波コイル１６ａ〜１６ｆにより印加される高周波磁場の方向と、ゼロ磁場領域の位置と、ゼロ磁場領域周辺の磁場分布とから磁気モーメントの方向を特定し、特定した磁気モーメントの方向および検出コイル１７により検出された信号値に基づいて磁気モーメントの大きさを求めることによって、磁化の分布を計測する。これにより、ゼロ磁場領域を固定した状態でも同時に複数の点における磁化を計測することができるようになり、磁性微粒子の３次元分布を画像化する場合の撮影時間を短縮することが可能になる。
【００６１】
ところで、上記実施例では、直線状のゼロ磁場領域を形成する場合について説明したが、本発明はこれに限られるわけではなく、例えば、円環状のゼロ磁場領域を形成するようにしてもよい。図１３は、円環状のゼロ磁場領域を形成する場合のコイルの配置を示す図である。また、図１４は、図１３に示したコイルによって発生する静磁場の磁場分布を示す模式図であり、図１５は、円環状のゼロ磁場領域を示す図である。
【００６２】
例えば、図１３に示すように、３つのゼロ磁場領域形成コイル１８ａ、１８ｂおよび１８ｃを、それぞれの中心軸が一致するように配置する。このようにゼロ磁場領域形成コイル１８ａ、１８ｂおよび１８ｃを配置することによって、図１４に示すように静磁場が発生する。これにより、図１５に示す円環状のゼロ磁場領域を形成することができる。
【００６３】
ここで、ゼロ磁場領域の円環の直径は、ゼロ磁場領域形成コイル１８ａおよび１８ｂに流す電流とゼロ磁場領域形成コイル１８ｃに流す電流の比を変えることによって変化させることができる。なお、ゼロ磁場領域形成コイル１８ａおよび１８ｂは、コイル（電磁石）の代わりに永久磁石を用いて構成することもできる。その場合には、ゼロ磁場領域形成コイル１８ｃに流す電流によって円環の直径を変化させる。また、頭−足方向へのゼロ磁場領域の移動は、寝台１３を移動することによって行う。
【００６４】
また、上記実施例では、線状（直線状または円環状など）のゼロ磁場領域を形成する場合について説明したが、本発明はこれに限られるわけではなく、従来のように点状のゼロ磁場領域を形成する場合にも同様に適用することができる。
【００６５】
図１６は、点状のゼロ磁場領域を形成する場合の磁性微粒子イメージング装置２０の構成を示す図である。なお、ここでは、図１に示した各部と同様の役割を果たすものについては、同一の符号を付すこととしてその説明を省略する。また、ここでは、制御装置１５については図示を省略している。
【００６６】
点状のゼロ磁場領域を形成する場合、磁性微粒子イメージング装置２０には、同図に示すように、２つのゼロ磁場領域形成磁石２１ａおよび２１ｂ、８つのゼロ磁場スキャンコイル２２ａ〜２２ｈ、６つの高周波コイル２６ａ〜２６ｆがそれぞれ配設される。
【００６７】
ゼロ磁場領域形成磁石２１ａおよび２１ｂは、対向して設置され、互いに反対向きの静磁場を発生させる。これにより、ゼロ磁場領域形成磁石２１ａと２１ｂの中間点で静磁場が相殺され、点状のゼロ磁場領域が形成される。
【００６８】
ゼロ磁場スキャンコイル２２ａ〜２２ｈは、撮影空間内でゼロ磁場領域を左−右および足−頭の方向に移動させる。具体的には、ゼロ磁場スキャンコイル２２ａおよび２２ｅからゼロ磁場スキャンコイル２２ｂおよび２２ｆへ向かう磁場と、ゼロ磁場スキャンコイル２２ｃおよび２２ｇからゼロ磁場スキャンコイル２２ｄおよび２２ｈへ向かう磁場が発生するように各コイルに電流を流せば、ゼロ磁場領域は左方向へ移動する。また、これとは逆向きの磁場が発生するように各コイルに電流を流せば、ゼロ磁場領域は右方向へ移動する。
【００６９】
一方、ゼロ磁場スキャンコイル２２ａおよび２２ｂからゼロ磁場スキャンコイル２２ｅおよび２２ｆへ向かう磁場と、ゼロ磁場スキャンコイル２２ｃおよび２２ｄからゼロ磁場スキャンコイル２２ｇおよび２２ｈへ向かう磁場が発生するように各コイルに電流を流せば、ゼロ磁場領域は頭方向へ移動する。また、これとは逆向きの磁場が発生するように各コイルに電流を流せば、ゼロ磁場領域は足方向へ移動する。
【００７０】
高周波コイル２６ａ〜２６ｆは、図５−１、５−２および６に示した高周波コイル１６ａ〜１６ｆと同様に配設される。
【００７１】
また、図１６では図示を省略しているが、磁性微粒子イメージング装置２０には、検出コイルも配設される。点状のゼロ磁場領域を形成する場合も、検出コイルは、図１１−１および１１−２に示した検出コイル１７と同様に、高周波コイル２６ｅおよび２６ｆが作る磁場および高周波コイル２６ａ〜２６ｄが作る磁場の両方に対してコイルの軸が直交するように配置される。
【００７２】
上記の構成によれば、高周波コイル２６ａ〜２６ｆによって形成される磁場、および、検出コイルによって計測される検出信号は、それぞれ直線状のゼロ磁場領域の場合と同様になる。したがって、ゼロ磁場領域の周方向に沿った磁化の分布は、直線状のゼロ磁場領域の場合と同様の計測方法で計測することができる。
【００７３】
以上のように、本実施例に係る磁性微粒子イメージング装置は、高周波コイルによってゼロ磁場領域に印加される磁場の方向を変化させ、特に回転磁場を生成する。また、高周波コイルが作る磁場、信号発生領域（ゼロ磁場近傍でかつ、磁場が加算される領域から信号が発生する）、および、磁気モーメント（磁性微粒子の濃度と大きさが関連する）の方向を、高周波コイルの信号と連動して変化させるとともに、高周波コイルが作っている磁場の方向から磁気モーメントの方向、および、信号発生領域を決定する手段と、複数の検出コイルの信号値と、磁気モーメントの方向から、磁気モーメントの大きさを決定する。
【００７４】
これにより、本実施例に係る磁性微粒子イメージング装置では、ゼロ磁場領域周辺の磁化分布を計測することが可能になった。従来は、ゼロ磁場領域一点につき、一点の磁化の大きさを計測できるだけであったが、本実施例に係る磁性微粒子イメージング装置では、ゼロ磁場領域一点につき、ゼロ磁場領域近傍の磁化の分布を計測できるようになった。したがって、従来と同一の点数の磁化分布を撮影すると仮定した場合、測定点数を従来と比べて少なくすることができ、その結果、高速撮影が可能になる。
【００７５】
また、本実施例に係る磁性微粒子イメージング装置は、線状のゼロ磁場領域を形成し、複数の検出器の信号から、線状ゼロ磁場に沿った方向の磁化分布を求める。これにより、一度に複数点の磁化分布を計測でき、撮影時間を短縮することが可能になる。特に、線状ゼロ磁場領域内の位置に応じて高周波コイルが作る磁場の向きが異なるように高周波コイルを配置するので、同一方向の磁場を印加する場合と比べて分解能が向上する。
【産業上の利用可能性】
【００７６】
以上のように、本発明に係る磁性微粒子イメージング装置は、磁性微粒子の３次元分布を画像化する場合に有用であり、特に、撮影時間を短縮することが要求される場合に適している。
【図面の簡単な説明】
【００７７】
【図１】本実施例に係る磁性微粒子イメージング装置の構成を示す図（１）である。
【図２】本実施例に係る磁性微粒子イメージング装置の構成を示す図（２）である。
【図３】ゼロ磁場領域形成磁石によって形成されるゼロ磁場領域を説明するための図（１）である。
【図４】ゼロ磁場領域形成磁石によって形成されるゼロ磁場領域を説明するための図（２）である。
【図５−１】高周波コイルによる高周波磁場の印加を説明するための図（１）である。
【図５−２】高周波コイルによる高周波磁場の印加を説明するための図（２）である。
【図６】高周波コイルの配置を説明するための図である。
【図７】ゼロ磁場領域の近傍で形成される磁場を説明するための図である。
【図８】図７に示したゼロ磁場領域の近傍を拡大した拡大図である。
【図９】ある瞬間において高周波コイルが作る磁場の様子を示す模式図である。
【図１０】図９に示した瞬間において計測される磁気モーメントを示す模式図である。
【図１１−１】ゼロ磁場領域の周方向に沿った磁化の分布を計測する方法を説明するための図（１）である。
【図１１−２】ゼロ磁場領域の周方向に沿った磁化の分布を計測する方法を説明するための図（２）である。
【図１２】ゼロ磁場領域の軸方向に沿った磁化の分布を計測する方法を説明するための図である。
【図１３】円環状のゼロ磁場領域を形成する場合のコイルの配置を示す図である。
【図１４】図１３に示したコイルによって発生する静磁場の磁場分布を示す模式図である。
【図１５】円環状のゼロ磁場領域を示す図である。
【図１６】点状のゼロ磁場領域を形成する場合の磁性微粒子イメージング装置の構成を示す図である。
【図１７】従来の磁性微粒子イメージングの原理を説明するための図である。
【符号の説明】
【００７８】
１０，２０    磁性微粒子イメージング装置
１１ａ，１１ｂ，２１ａ，２１ｂ    ゼロ磁場領域形成磁石
１２ａ，１２ｂ，２２ａ〜２２ｈ    ゼロ磁場スキャンコイル
１３    寝台
１４    架台
１５    制御装置
１５ａ    コイル制御部
１５ｂ    磁化検出部
１５ｃ    画像生成部
１６ａ〜１６ｌ，２６ａ〜２６ｆ    高周波コイル
１７    検出コイル
【図１】
2010088683.tif 000005
【図２】
2010088683.tif 000006
【図３】
2010088683.tif 000007
【図４】
2010088683.tif 000008
【図５−１】
2010088683.tif 000009
【図５−２】
2010088683.tif 000010
【図６】
2010088683.tif 000011
【図７】
2010088683.tif 000012
【図８】
2010088683.tif 000013
【図９】
2010088683.tif 000014
【図１０】
2010088683.tif 000015
【図１１−１】
2010088683.tif 000016
【図１１−２】
2010088683.tif 000017
【図１２】
2010088683.tif 000018
【図１３】
2010088683.tif 000019
【図１４】
2010088683.tif 000020
【図１５】
2010088683.tif 000021
【図１６】
2010088683.tif 000022
【図１７】
2010088683.tif 000023