(57)【要約】【課題】温度変化で共鳴周波数が変化しにくく、磁場変化で共鳴周波数が変化する特性が新たに発見された物質を用いて被検体内部の温度変化を計測できる温度計測方法を提供する。【解決手段】温度計測方法は、水の共鳴周波数を含む高周波磁場の印加で被検体から検出される磁気共鳴信号ｆＨ（ｔ）を周波数スペクトルＦＨ（ω）に変換し、リン化合物の共鳴周波数を含む高周波磁場の印加で被検体から検出される磁気共鳴信号ｆＰ（ｔ）を周波数スペクトルＦＰ（ω）に変換するステップと、周波数スペクトルＦＨ（ω）の中から水の共鳴周波数に相当する周波数ｆ１Ｈを取得し、周波数スペクトルＦＰ（ω）の中からリン化合物の共鳴周波数に相当する周波数ｆ１Ｐを取得するステップと、水の共鳴周波数の変化量と温度との間の相関関係を利用して、周波数ｆ１Ｈの変化量ΔｆＨと、周波数ｆ１Ｐの変化量ΔｆＰとの差に応じた温度変化を検出するステップとを備える。【選択図】図８
2016168087.tif 000002
【特許請求の範囲】
【請求項１】
静磁場中に存在する被検体に対して、水の共鳴周波数を含む第１高周波磁場を印加し、前記第１高周波磁場の印加後に、リン化合物の共鳴周波数を含む第２高周波磁場を印加するステップと、
前記第１高周波磁場の印加によって前記被検体から検出される磁気共鳴信号を、周波数ごとの信号強度を示す第１周波数スペクトルに変換するとともに、前記第２高周波磁場の印加によって前記被検体から検出される磁気共鳴信号を第２周波数スペクトルに変換するステップと、
前記第１周波数スペクトルの中から信号強度が相対的に高く、水の共鳴周波数に相当する第１周波数を取得するとともに、前記第２周波数スペクトルの中から信号強度が相対的に高く、リン化合物の共鳴周波数に相当する第２周波数を取得するステップと、
水の共鳴周波数の変化量と温度との間の相関関係を利用して、前記第１周波数の変化量と、前記第２周波数の変化量との差に応じた温度変化を検出するステップとを備える、温度計測方法。
【請求項２】
前記印加するステップは、前記第１高周波磁場の印加の停止後に連続して、前記第２高周波磁場を印加する、請求項１に記載の温度計測方法。
【請求項３】
前記検出するステップは、前記第１周波数の変化量から前記第２周波数の変化量を差分し、当該差分の結果に応じた温度変化を検出する、請求項１または２に記載の温度計測方法。
【請求項４】
前記リン化合物は、ホスホクレアチンを含む、請求項１〜３のいずれか１項に記載の温度計測方法。
【請求項５】
静磁場中に存在する被検体に対して、リン化合物の共鳴周波数を含む高周波磁場を印加するステップと、
前記高周波磁場の印加によって前記被検体から検出される磁気共鳴信号を、周波数ごとの信号強度を示す周波数スペクトルに変換するステップと、
前記周波数スペクトルの中から信号強度が相対的に高く、リン化合物の共鳴周波数に相当する周波数を取得するステップと、
リン化合物の共鳴周波数の変化量と磁場強度との間の相関関係を利用して、前記周波数の変化量に応じた磁場変化を検出するステップとを備える、磁場計測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、被検体内部の非侵襲的な温度計測方法および磁場計測方法に関し、特に、リン化合物の特性に基づいた温度計測方法および磁場計測方法に関する。
【背景技術】
【０００２】
近年、被検体内部の温度を非侵襲的に計測をすることが可能なＭＲＩ（Magnetic  Resonance  Imaging）が開発されている。たとえば、非特許文献１は、静磁場中にある水の共鳴周波数が温度に依存して変化することを利用して脳内の温度を非侵襲的に計測する方法を開示している。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】Rounald J.T.Corbett他,“Validation of a Noninvasive Method to Measure Brain Temperature In Vivo Using 1H NMR Spectroscopy”，Journal of Neurochemistry, 1995 Mar, 64(3), 1224-30
【発明の概要】
【発明が解決しようとする課題】
【０００４】
水の共鳴周波数は、温度の他にも様々な要因に依存して変化する。他の要因としては、たとえば、磁場変化が挙げられる。非特許文献１に開示される温度計測方法は、このような磁場変化による共鳴周波数の変化を排除するために、Ｎ−アセチルアスパラギン酸（ＮＡＡ）を用いる。ＮＡＡの共鳴周波数は、温度変化の影響を受けにくく、磁場変化の影響を受けやすい。一方で、水の共鳴周波数は、上述したように、温度変化と磁場変化との両方の影響を受けやすい。そのため、ＮＡＡの共鳴周波数を基準として、水の共鳴周波数に着目すれば、磁場変化等の影響を少なくすることができる。これにより、非特許文献１に開示される温度計測方法は、脳内の温度を高精度に計測する。
【０００５】
しかしながら、脳内のＮＡＡは、場所や個人差によって絶対量が増減するため、場合によっては、ＮＡＡの磁気共鳴信号の検出が困難であることがある。また、脳以外のＮＡＡが存在しない部位に対しては、この方法で温度を計測することができない。そのため、ＮＡＡと同様の特性を有する他の物質の発見が望まれている。すなわち、温度変化によって共鳴周波数が変化しにくく、磁場変化によって共鳴周波数が変化し、生体内に広く存在する他の物質（以下、「基準物質」ともいう。）の発見が望まれている。
【０００６】
本開示は上述のような問題点を解決するためになされたものであって、その目的は、温度変化によって共鳴周波数が変化しにくく、磁場変化によって共鳴周波数が変化する特性が新たに発見された物質を用いて被検体の内部の温度変化を計測することが可能な温度計測方法を提供することである。他の局面における目的は、温度変化によって共鳴周波数が変化しにくく、磁場変化によって共鳴周波数が変化する特性が新たに発見された物質を用いて被検体を含む領域の磁場変化を計測することが可能な磁場計測方法を提供することである。
【課題を解決するための手段】
【０００７】
ある局面に従うと、温度計測方法は、静磁場中に存在する被検体に対して、水の共鳴周波数を含む第１高周波磁場を印加し（選択的高周波磁場の印加によるボクセル信号収集を含む）、第１高周波磁場の印加後に、リン化合物の共鳴周波数を含む第２高周波磁場を印加するステップと、第１高周波磁場の印加によって被検体から検出される磁気共鳴信号を、周波数ごとの信号強度を示す第１周波数スペクトルに変換するとともに、第２高周波磁場の印加によって被検体から検出される磁気共鳴信号を第２周波数スペクトルに変換するステップと、第１周波数スペクトルの中から信号強度が相対的に高く、水の共鳴周波数に相当する第１周波数を取得するとともに、第２周波数スペクトルの中から信号強度が相対的に高く、リン化合物の共鳴周波数に相当する第２周波数を取得するステップと、水の共鳴周波数の変化量と温度との間の相関関係を利用して、第１周波数の変化量と、第２周波数の変化量との差に応じた温度変化を検出するステップとを備える。
【０００８】
好ましくは、印加するステップは、第１高周波磁場の印加の停止後に連続して、第２高周波磁場を印加する。
【０００９】
好ましくは、検出するステップは、第１周波数の変化量から第２周波数の変化量を差分し、当該差分の結果に応じた温度変化を検出する。
【００１０】
好ましくは、リン化合物は、ホスホクレアチンを含む。
他の局面に従うと、磁場計測方法は、静磁場中に存在する被検体に対して、リン化合物の共鳴周波数を含む高周波磁場を印加するステップと、高周波磁場の印加によって被検体から検出される磁気共鳴信号を、周波数ごとの信号強度を示す周波数スペクトルに変換するステップと、周波数スペクトルの中から信号強度が相対的に高く、リン化合物の共鳴周波数に相当する周波数を取得するステップと、リン化合物の共鳴周波数の変化量と磁場強度との間の相関関係を利用して、周波数の変化量に応じた磁場変化を検出するステップとを備える。
【発明の効果】
【００１１】
ある局面において、温度変化によって共鳴周波数が変化しにくく、磁場変化によって共鳴周波数が変化する特性が新たに発見された物質を用いて被検体の内部の温度変化を計測することができる。
【００１２】
本発明の上記および他の目的、特徴、局面および利点は、添付の図面と関連して理解される本発明に関する次の詳細な説明から明らかとなるであろう。
【図面の簡単な説明】
【００１３】
【図１】第１の実施の形態に従うＭＲＩの外観を示した図である。
【図２】温度を変化させ、かつ磁場を一定にした場合における、水、ＮＡＡ、および有機リン酸の化学シフトを示した図である。
【図３】温度を変化させ、かつ磁場を一定にした場合における、水、ＮＡＡ、および無機リン酸の化学シフトの時間的な変化を示した図である。
【図４】温度を一定にし、かつ磁場を変化させた場合における、有機リン酸およびＮＡＡの共鳴周波数の時間的な変化を示した図である。
【図５】温度を一定にし、かつ磁場を変化させた場合における、無機リン酸およびＮＡＡの共鳴周波数の時間的な変化を示した図である。
【図６】第１の実施の形態に従うＭＲＩの主要な構成を示した図である。
【図７】被検体に印加する高周波磁場のタイミングを示す図である。
【図８】第１の実施の形態に従う計算機システムによる温度計測方法を概略的に示した概念図である。
【図９】人の頭部におけるリン化合物の周波数スペクトルをグラフとして示す図である。
【図１０】第１の実施の形態に従う温度計測方法による温度計測結果を示す図である。
【図１１】第１の実施の形態に従う温度計測方法で検出された脳内の温度分布を示した図である。
【図１２】第１の実施の形態に従う計算機システムが実行する処理の一部を表わすフローチャートである。
【図１３】第１の実施の形態に従う計算機システムの主要なハードウェア構成を示すブロック図である。
【図１４】第２の実施の形態に従う計算機システムによる磁場計測方法を概略的に示した概念図である。
【図１５】第２の実施の形態に従う計算機システムが実行する処理の一部を表わすフローチャートである。
【発明を実施するための形態】
【００１４】
以下、図面を参照しつつ、本発明に従う各実施の形態について説明する。以下の説明では、同一の部品および構成要素には同一の符号を付してある。それらの名称および機能も同じである。したがって、これらについての詳細な説明は繰り返さない。また、以下で説明する各実施の形態または変形例は、適宜選択的に組み合わされてもよい。
【００１５】
＜新たな基準物質＞
上述したように、温度変化によって共鳴周波数が変化しにくく、磁場変化によって共鳴周波数が変化する、ＮＡＡに代わる物質（すなわち、基準物質）が望まれている。発明者らは、リン化合物が新たな基準物質となり得ることを発見した。
【００１６】
図１〜図５を参照して、リン化合物が新たな基準物質となり得る根拠について説明する。図１は、ＭＲＩ１００の外観を示した図である。発明者らは、ＭＲＩ１００を用いてリン化合物が基準物質となり得ることについての実験を行なった。当該実験では、ＭＲＩ１００として、シーメンス（登録商標）社製のＭＡＧＥＮＴＯＭ  ７Ｔが用いられた。
【００１７】
その結果、発明者らは、リン化合物の共鳴周波数が、温度変化によらずほぼ一定であることを新たに発見した。さらに、発明者らは、磁場変化によるリン化合物の共鳴周波数の変化が、ＮＡＡの共鳴周波数の変化と類似することを新たに発見した。
【００１８】
図２〜図５には、その実験結果が示されている。図２は、温度を変化させ、かつ磁場を一定にした場合における、水（典型的には、プロトン（１Ｈ））、ＮＡＡ、および有機リン酸の共鳴周波数の時間的な変化（以下、「化学シフト」ともいう。）を示した図である。図２の実験では、ＮＡＡと、有機リン酸の一種であるホスホクレアチン（ＰＣｒ）との混合水溶液がファントムとして使用された。
【００１９】
図２のグラフ５１の横軸は、時間変化を示す。グラフ５１の縦軸は、化学シフトを示す。グラフ５１には、時間「０」を基準とする各物質の化学シフトが実験結果として示されている。グラフ５１に示されるように、ＮＡＡおよび有機リン酸の化学シフトは、水の化学シフトよりも小さい。また、有機リン酸の化学シフトは、ＮＡＡの化学シフトと同様に、温度変化によらずほぼ一定である。
【００２０】
図３は、温度を変化させ、かつ磁場を一定にした場合における、水、ＮＡＡ、および無機リン酸の化学シフトの時間的な変化（すなわち、化学シフト）を示した図である。図３の実験では、ＮＡＡと、無機リン酸の一種であるリン酸二水素ナトリウム二水和物との混合水溶液がファントムとして使用された。
【００２１】
図３のグラフ５３の横軸は、時間変化を示す。グラフ５３の縦軸は、化学シフトを示す。グラフ５３には、時間「０」を基準とする各物質の化学シフトが実験結果として示されている。グラフ５３に示されるように、ＮＡＡおよび無機リン酸の化学シフトは、水の化学シフトよりも小さい。また、無機リン酸の化学シフトは、ＮＡＡの化学シフトと同様に、温度変化によらずほぼ一定である。
【００２２】
図２および図３の実験結果から、有機リン酸や無機リン酸等のリン化合物の共鳴周波数は、温度に依存せずほぼ一定であることが確認された。
【００２３】
図４は、温度を一定にし、かつ磁場の自然ドリフト時における、有機リン酸およびＮＡＡの共鳴周波数の時間的な変化を示した図である。図４のグラフ５５，５７の横軸は、時間を示す。グラフ５５の縦軸は、化学シフトを示す。グラフ５７の縦軸は、有機リン酸の化学シフトとＮＡＡの化学シフトとの差を示す。
【００２４】
グラフ５５に示されるように、有機リン酸およびＮＡＡは、磁場が自然ドリフトする環境下では、化学シフトの変化の傾向が類似している。また、グラフ５７に示されるように、有機リン酸の化学シフトとＮＡＡの化学シフトとの差は、±０．００２ｐｐｍで収まっている。このことからも、有機リン酸およびＮＡＡで化学シフトの変化の傾向がほぼ同じであることを確認できる。
【００２５】
図５は、温度を一定にし、かつ磁場の自然ドリフト時における、無機リン酸およびＮＡＡの共鳴周波数の時間的な変化を示した図である。図５のグラフ５９，６１の横軸は、時間を示す。グラフ５９の縦軸は、化学シフトを示す。グラフ６１の縦軸は、有機リン酸の化学シフトとＮＡＡの化学シフトとの差を示す。
【００２６】
グラフ５９に示されるように、無機リン酸およびＮＡＡは、磁場が自然ドリフトする環境下では、化学シフトの変化の傾向が類似している。また、グラフ６１に示されるように、無機リン酸の化学シフトとＮＡＡの化学シフトとの差は、±０．０００７ｐｐｍで収まっている。このことからも、無機リン酸およびＮＡＡで化学シフトの変化の傾向がほぼ同じであることを確認できる。
【００２７】
図４および図５の実験結果から、有機リン酸や無機リン酸等のリン化合物の共鳴周波数は、磁場が変化する環境下ではＮＡＡと同じように変化することとが確認された。
【００２８】
以上の実験結果から、リン化合物の共鳴周波数は、ＮＡＡと同様に、温度不変性と磁場可変性とを有することが示された。このことは、リン化合物がＮＡＡに代わる新たな基準物質となり得ることを示している。
【００２９】
発明者らの知る限り、リン化合物の上記特性を示す先行技術文献は存在しない。このことは、以下のような事情が背景にあったものと考えられる。
【００３０】
リン化合物の磁気共鳴信号は非常に微弱であり、この磁気共鳴信号を正確に捉えること自体が困難であった。しかしながら、今回の実験では、７Ｔ（テラス）の超高磁場の静磁場強度を発生することが可能なＭＲＩが用いられ、リン化合物の磁気共鳴信号を高精度に捉えることができた。このような世界最先端のＭＲＩは、装置自体が高価である等の理由により、日本国内でもほとんど普及していない。超高性能なＭＲＩを用いなければリン化合物の上記特性を確認することができなかったことが、当該特性が発見されていなかった理由の一つである考えられる。
【００３１】
リン化合物において新たに発見された特性は、様々なことに応用され得る。たとえば、当該特性は、被検体内部の温度計測方法や磁場計測方法に応用され得る。以下では、当該温度計測方法を利用する第１の実施の形態に従うＭＲＩ１００と、当該磁場計測方法を利用する第２の実施の形態に従うＭＲＩ１００とについて順に説明する。
【００３２】
＜第１の実施の形態＞
［ＭＲＩ１００］
図６および図７を参照して、第１の実施の形態に従うＭＲＩ１００の装置概要について説明する。図６は、ＭＲＩ１００の主要な構成を示した図である。図７は、被検体Ｓに印加する高周波磁場のタイミングを示す図である。
【００３３】
図６に示されるように、ＭＲＩ１００は、静磁場磁石１と、傾斜磁場コイル２と、送信部３と、送信コイル４と、受信コイル５と、受信部６と、電源１０と、シーケンス制御部１３と、寝台２０と、天板２１と、計算機システム２００とを備える。
【００３４】
静磁場磁石１は、中空の円筒形状に形成された磁石であり、内部の空間に一様な静磁場を発生する。静磁場磁石１としては、たとえば、永久磁石、超伝導磁石等が用いられる。好ましくは、静磁場磁石１としては、７Ｔ程度の静磁場強度を発生するものが用いられる。なお、静磁場磁石１の静磁場強度は、７Ｔ程度に限定されない。静磁場強度は、条件や環境が最適化されれば、７Ｔ以下（たとえば、１．５Ｔ程度、３Ｔ程度）であってもよい。あるいは、静磁場強度は、７Ｔ以上であってもよい。
【００３５】
傾斜磁場コイル２は、中空の円筒形状に形成されたコイルであり、静磁場磁石１の内側に配置される。この傾斜磁場コイル２は、互いに直交するＸ，Ｙ，Ｚ軸にそれぞれ対応する３つのコイルが組み合わされて形成されている。これら３つのコイルは、電源１０から個別に電流供給を受けて、Ｘ，Ｙ，Ｚ軸のそれぞれに沿って磁場強度が変化する傾斜磁場を発生させる。なお、Ｚ軸方向は、静磁場と同方向とする。
【００３６】
傾斜磁場コイル２によって発生するＸ，Ｙ，Ｚ軸の各々の傾斜磁場は、たとえば、スライス選択用傾斜磁場Ｇｓ、位相エンコード用傾斜磁場Ｇｅおよびリードアウト用傾斜磁場Ｇｒにそれぞれ対応する。スライス選択用傾斜磁場Ｇｓは、任意に撮像断面を決めるために利用される。位相エンコード用傾斜磁場Ｇｅは、空間的位置に応じて磁気共鳴信号の位相を変化させるために利用される。リードアウト用傾斜磁場Ｇｒは、空間的位置に応じて磁気共鳴信号の周波数を変化させるために利用される。
【００３７】
送信部３は、シーケンス制御部１３による制御のもと、ラーモア周波数に対応する高周波信号を送信コイル４に送信する。送信コイル４は、傾斜磁場コイル２の内側に配置され、高周波信号（ＲＦ（Radio  Frequency）信号）を送信部３から受けて、被検体Ｓに高周波磁場を印加する。
【００３８】
本実施の形態においては、送信部３は、水の共鳴周波数を含む高周波磁場（以下、「水用高周波磁場」ともいう。）と、リン化合物の共鳴周波数を含む高周波磁場（以下、「リン用高周波磁場」ともいう。）とを送信コイル４に交互に発生させる。
【００３９】
好ましくは、送信コイル４は、水用高周波磁場の印加時には、水の共鳴周波数を基準とする所定バンド幅（たとえば、±２０００Ｈｚ）で水用高周波磁場を印加する。これにより、水以外の物質から放射される磁気共鳴信号が排除され、後述の受信コイル５は、水の磁気共鳴信号を精度よく受信することができる。また、送信コイル４は、リン用高周波磁場の印加時には、リン化合物の共鳴周波数を中心とする所定バンド幅（たとえば、±３０００Ｈｚまたは±４０００Ｈｚ）で水用高周波磁場を印加する。これにより、リン化合物以外の物質から放射される磁気共鳴信号が排除され、後述の受信コイル５は、リン化合物の磁気共鳴信号を精度よく受信することができる。
【００４０】
好ましくは、図７に示されるように、送信部３は、水用高周波磁場の印加の停止後に連続して、リン用高周波磁場を印加する。これにより、後述の受信コイル５は、水とリン化合物との磁気共鳴信号を、ほぼ同じタイミングで、ほぼ同じ環境下で取得することができる。送信コイル４は、たとえば、水用高周波磁場およびリン用高周波磁場を５秒ずつ交互に印加する。なお、リン化合物の種類によっては磁気共鳴信号が微弱である可能性があるので、リン用高周波磁場は、複数回印加されてもよい。あるいは、リン用高周波磁場の印加時間は、水用高周波磁場の印加時間よりも長くてもよい。
【００４１】
受信コイル５は、複数のチャンネルの磁気共鳴信号を受信することが可能なものが用いられる。本実施の形態では、受信コイル５は、水用高周波磁場の印加によって被検体Ｓから放出される磁気共鳴信号と、リン用高周波磁場の印加によって被検体Ｓから放出される磁気共鳴信号との少なくとも２つの磁気共鳴信号を受信する。受信コイル５は、受信した磁気共鳴信号を受信部６に出力する。一例として、受信コイル５には、ＱＥＤ社製のMulti-Nuclear Loop Coilが用いられる。
【００４２】
受信部６は、シーケンス制御部１３による制御のもと、受信コイル５から出力される磁気共鳴信号をＡ／Ｄ（Analog-to-Digital）変換することで、時間領域の信号データを生成する。受信部６によって生成される信号データには、上述したスライス選択用傾斜磁場Ｇｓ、位相エンコード用傾斜磁場Ｇｅ、およびリードアウト用傾斜磁場Ｇｒによって、ｋ空間におけるＰＥ（Phase  Encode）方向、ＲＯ（Read  Out）方向、ＳＥ（Slice  Encode）方向の位置を示す成分がエンコードされる。
【００４３】
シーケンス制御部１３は、計算機システム２００から送信されるシーケンス実行データにしたがって送信部３、受信部６、および電源１０を駆動することによって、被検体Ｓの画像を生成するためのデータ収集を行なう。シーケンス実行データとは、被検体Ｓから生データを収集するためのパルスシーケンスを定義する情報である。より具体的には、シーケンス実行データは、電源１０によって傾斜磁場コイル２に供給される電源の強さや電力を供給するタイミング、送信部３が送信コイル４に送信する高周波信号のレベルや高周波信号を送信するタイミング、受信部６が磁気共鳴信号を検出するタイミング等、データ収集を行なうための手順を定義した情報である。
【００４４】
シーケンス制御部１３は、電源１０、送信部３および受信部６を駆動して被検体Ｓをスキャンした結果、受信部６から生データが送信されると、その生データを計算機システム２００へ転送する。
【００４５】
寝台２０には、被検体Ｓが載置される天板２１が取り付けられている。天板２１は、被検体Ｓが載置された状態で傾斜磁場コイル２の空洞（撮像口）内へ挿入される。通常、この寝台２０は、長手方向が静磁場磁石１の中心軸と平行になるように設置される。
【００４６】
計算機システム２００は、ＭＲＩ１００の全体制御を行なう。たとえば、計算機システム２００は、上述した各部を駆動することで、データ収集や画像再構成等を行なう。計算機システム２００は、収集したデータに基づいて、被検体の温度計測や磁場計測を行なう。
【００４７】
［温度計測方法］
図８を参照して、リン化合物の上記特性を利用した温度計測方法について説明する。図８は、計算機システム２００による温度計測方法を概略的に示した概念図である。計算機システム２００は、機能構成として、ＦＦＴ（Fast Fourier Transform）部２５１，２５２と、抽出部２５３，２５４と、算出部２５５，２５６と、差分部２５７と、検出部２５８と、特定部２６０とを含む。
【００４８】
ＦＦＴ部２５１は、水用高周波磁場の印加によって被検体から検出される磁気共鳴信号ｆＨ（ｔ）を、周波数ごとの信号強度を示す周波数スペクトルＦＨ（ω）に変換する。すなわち、ＦＦＴ部２５１は、時間領域の磁気共鳴信号ｆＨ（ｔ）を、周波数領域の周波数スペクトルＦＨ（ω）に変換する。一例として、当該変換処理には、高速フーリエ変換が用いられる。
【００４９】
ＦＦＴ部２５２は、リン用高周波磁場の印加によって被検体から検出される磁気共鳴信号ｆＰ（ｔ）を周波数スペクトルＦＰ（ω）に変換する。すなわち、ＦＦＴ部２５２は、時間領域の磁気共鳴信号ｆＰ（ｔ）を、周波数領域の周波数スペクトルＦＰ（ω）に変換する。一例として、当該変換処理には、高速フーリエ変換が用いられる。
【００５０】
抽出部２５３は、周波数スペクトルＦＨ（ω）の中から信号強度が相対的に高く、水の共鳴周波数に相当する周波数ｆ１Ｈを抽出する。より具体的には、抽出部２５３は、水の共鳴周波数を中心とする所定範囲のスペクトルの中から最大の信号強度を選択し、選択された信号強度に対応する周波数を、水の共鳴周波数に相当する周波数ｆ１Ｈとして抽出する。
【００５１】
抽出部２５４は、周波数スペクトルＦＰ（ω）の中から信号強度が相対的に高く、リン化合物の共鳴周波数に相当する周波数ｆ１Ｐを抽出する。より具体的には、抽出部２５４は、リン化合物の共鳴周波数を中心とする所定範囲のスペクトルの中から最大の信号強度を選択し、選択された信号強度に対応する周波数を、リン化合物の共鳴周波数に相当する周波数ｆ１Ｐとして抽出する。
【００５２】
算出部２５５は、水の共鳴周波数の時間的な変化量ΔｆＨ（すなわち、化学シフト）を算出する。たとえば、算出部２５５は、周波数ｆ１Ｈよりも前に取得された水の共鳴周波数である周波数ｆ０Ｈを周波数ｆ１Ｈから差分することで変化量ΔｆＨを算出する。
【００５３】
算出部２５６は、リン化合物の時間的な変化量ΔｆＰ（すなわち、化学シフト）を算出する。たとえば、算出部２５６は、周波数ｆ１Ｐよりも前に取得されたリン化合物の共鳴周波数である周波数ｆ０Ｐを周波数ｆ１Ｐから差分することで変化量ΔｆＰを算出する。
【００５４】
なお、周波数ｆ０Ｈと周波数ｆ０Ｐとは、同じタイミングまたはほぼ同じタイミングで検出される。あるいは、周波数ｆ０Ｈと周波数ｆ０Ｐとは、予め設定されていてもよい。周波数ｆ１Ｈと周波数ｆ１Ｐとは、同じタイミングまたはほぼ同じタイミングで検出される。
【００５５】
ここで、図９を参照して、算出部２５６による算出処理についてさらに詳細に説明する。図９は、周波数スペクトルＦＰ（ω）をグラフ６４として示した図である。グラフ６４の横軸は、化学シフトを示す。グラフ６４の縦軸は、磁気共鳴信号の強度を示す。
【００５６】
グラフ６４の横軸は、過去の所定タイミングにおけるリン化合物の共鳴周波数を基準とする化学シフトを示す。すなわち、０ｐｐｍからのズレがリン化合物の化学シフトに相当する。算出部２５６は、０ｐｐｍを中心とする所定範囲（たとえば、±５ｐｐｍ）のスペクトルの中から最大となる信号強度を選択する。図９の例では、最大スペクトル６４Ａが選択される。算出部２５６は、最大スペクトル６４Ａに対応する化学シフトの値を変化量ΔｆＰとする。
【００５７】
再び図８を参照して、差分部２５７は、算出部２５５によって算出されたΔｆＨと算出部２５６によって算出されたΔｆＰとの差を変化量ΔｆＨ−Ｐ算出する。上述したように、水の共鳴周波数は、温度と磁場との両方に依存する。リン化合物の共鳴周波数は、上述の実験結果で示したように、温度には依存せず、磁場には依存する。そのため、リン化合物の共鳴周波数の変化量を基準として、水の共鳴周波数の変化量に着目すれば、磁場変化の影響を除外することができる。すなわち、変化量ΔｆＨ−Ｐは、温度変化のみによって生じた共鳴周波数の変化量を示す。
【００５８】
検出部２５８は、水の共鳴周波数の変化量と温度との間の相関関係を利用して、変化量ΔｆＨ−Ｐに応じた温度変化を検出する。水の共鳴周波数は、１℃の変化で約０．０１ｐｐｍ（より具体的には、０．０１０４ｐｐｍ）変化することが知られている。検出部２５８は、変化量ΔｆＨ−Ｐを相関値「０．０１ｐｐｍ／℃」で除算し、その除算結果を被検体の温度変化として検出する。
【００５９】
特定部２６０は、共鳴周波数ｆＨ（ｔ）が放出された被検体の位置を特定する。より具体的には、傾斜磁場コイル２（図６参照）は、上述したように、スライス選択用傾斜磁場Ｇｓ、位相エンコード用傾斜磁場Ｇｅおよびリードアウト用傾斜磁場Ｇｒを発生する。特定部２６０は、スライス選択用傾斜磁場Ｇｓ、位相エンコード用傾斜磁場Ｇｅ、およびリードアウト用傾斜磁場Ｇｒによって、ｋ空間におけるＰＥ方向、ＲＯ方向、ＳＥ方向の位置を示す成分をエンコードする。特定部２６０は、エンコード結果に基づいて、共鳴周波数ｆＨ（ｔ）が放出された被検体の位置を特定する。当該位置情報は、たとえば、３次元の座標（ｘ，ｙ，ｚ）で示される。
【００６０】
特定部２６０により、被検体内（たとえば、脳内）の単位領域ごとに、水およびリン化合物の共鳴周波数を特定することができる。これにより、検出部２５８は、被検体内の単位領域ごとに温度変化を検出することができる。
【００６１】
［実験結果］
図１０を参照して、本実施の形態に従う温度計測方法による温度計測結果について説明する。図１０は、本温度計測方法による温度計測結果を示す図である。
【００６２】
図１０（Ａ）には、本温度計測方法による温度計測結果と、光ファイバー温度計による温度計測結果との比較結果が示されている。グラフ６３の横軸は、時間変化を示す。グラフ６３の縦軸は、時間「０」を基準とする無機リン酸水溶液の温度変化を示す。温度変化６５は、光ファイバー温度計で無機リン酸水溶液の温度を計測した結果である。温度変化６７は、本温度計測方法で無機リン酸水溶液の温度を計測した結果である。
【００６３】
グラフ６３に示されるように、光ファイバー温度計によって検出された温度変化６５と、本温度計測方法によって検出された温度変化６７とは、同じような傾向を示している。温度変化６５，６７と間で若干の差が生じているが、本温度計測方法を利用すると、少なくとも、相対的な温度変化が生じていることは検出され得る。温度算出時に用いた相関値（上述の「０．０１ｐｐｍ／℃」）が本温度計測処理に合わせて最適化されると、温度検出精度は、さらに改善され得る。グラフ６３の結果により、リン化合物の上記特性が温度計測処理に応用できることが示された。
【００６４】
図１０（Ｂ）には、有機リン酸を基準物質とする温度計測方法と、ＮＡＡを基準物質とする温度計測方法とが示されている。グラフ７１の横軸は、時間変化を示す。グラフ７１の縦軸は、時刻「０」を基準とする水溶液の温度変化を示す。当該水溶液には、有機リン酸とＮＡＡとが含まれる。温度変化７３は、有機リン酸を基準物質とした温度計測結果である。温度変化７５は、ＮＡＡを基準物質とする温度計測結果である。
【００６５】
グラフ７１に示されるように、温度変化７３と温度変化７５とは、同じような傾向を示している。グラフ７１の結果により、有機リン酸がＮＡＡの代わりの基準物質となり得ることが示された。
【００６６】
［温度計測結果の出力態様］
図１１を参照して、本温度計測処理による温度計測結果の出力態様の一例について説明する。図１１は、本温度計測方法で検出された脳内の温度分布を示した図である。
【００６７】
温度計測結果は、たとえば、計算機システム２００（図６参照）に備えられる表示部２０５に表示される。一例として、表示部２０５には、異なるタイミングで計測された脳内の温度分布２２０Ａ，２２０Ｂが表示される。
【００６８】
たとえば、温度分布２２０Ａ，２２０Ｂは、温度ごとに関連付けられている色で示される。たとえば、温度が高いほど暖色系の色（たとえば、赤色）が用いられ、温度が低いほど寒色系の色（たとえば、青色）が用いられる。温度分布２２０Ａ，２２０Ｂが並べて表示されることで、脳内の時間的な温度変化が一見で理解され得る。
【００６９】
また、表示部２０５には、指標画像２２２が表示される。指標画像２２２は、温度と色との対応関係を示す。診断医は、指標画像２２２に示される色と、温度分布２２０Ａ，２２０Ｂとを比較することで、脳内における温度を一見して理解することができる。
【００７０】
なお、図１１では、温度分布２２０Ａ，２２０Ｂが並べて表示されている例が示されているが、温度分布２２０Ａ，２２０Ｂのいずれか一方だけが表示されてもよい。この場合、診断医は、脳内における空間的な温度変化を一見して理解することができる。
【００７１】
また、図１１では、異なるタイミングで得られた温度分布２２０Ａ，２２０Ｂが並べて表示されているが、脳内の異なる領域の温度分布が並べて表示されてもよい。これにより、診断医は、脳内の場所ごとの温度分布を比較しやすくなる。
【００７２】
さらに、図１１では、脳内の絶対的な温度が検出結果として表示されているが、脳内の相対的な温度が検出結果として表示されてもよい。たとえば、診断医は、温度分布２２０Ａ，２２０Ｂに対してタッチ操作等で特定の位置を指示する。計算機システム２００は、指定位置を基準とした脳内の相対的な温度を計算し、その計算結果に応じて温度分布２２０Ａ，２２０Ｂの表示を更新する。一例として、指定位置よりも温度が高い部分については暖色系の色（たとえば、赤色）を表示する。指定位置よりも温度が低い部分については寒色系の色（たとえば、青色）を表示する。このとき、指標画像２２２に示される温度と色との対応関係も更新される。
【００７３】
［計算機システム２００の制御構造］
図１２を参照して、第１の実施の形態に従う計算機システム２００の制御構造について説明する。図１２は、第１の実施の形態に従う計算機システム２００が実行する処理の一部を表わすフローチャートである。図１２の処理は、計算機システム２００のＣＰＵ（Central  Processing  Unit）２０２が温度計測プログラムを実行することにより実現される。他の局面において、処理の一部または全部が、回路素子その他のハードウェアによって実行されてもよい。
【００７４】
ステップＳ２０において、ＣＰＵ２０２は、リン化合物の共鳴周波数を含む高周波磁場（すなわち、リン用高周波磁場）を被検体に印加するための制御信号を送信部３（図６参照）に送る。これにより、リン用高周波磁場が被検体に印加される。
【００７５】
ステップＳ２２において、ＣＰＵ２０２は、リン用高周波磁場の印加によって被検体から放出される磁気共鳴信号を受信する。ＣＰＵ２０２は、上述のＦＦＴ部２５２（図８参照）として、受信した時間領域の磁気共鳴信号を周波数領域の周波数スペクトルに変換する。
【００７６】
ステップＳ２４において、ＣＰＵ２０２は、上述の抽出部２５４（図８参照）として、周波数スペクトルの中から信号強度が相対的に高く、リン化合物の共鳴周波数に相当する周波数を抽出する。ＣＰＵ２０２は、上述の算出部２５６（図８参照）として、過去に得られたリン化合物の共鳴周波数と、現在のリン化合物の共鳴周波数との間の差を算出する。当該差がリン化合物の共鳴周波数の化学シフトに相当する。
【００７７】
ステップＳ３０において、ＣＰＵ２０２は、水の共鳴周波数を含む高周波磁場（すなわち、水用高周波磁場）を被検体に印加するための制御信号を送信部３（図６参照）に送る。これにより、水用高周波磁場が被検体に印加される。
【００７８】
ステップＳ３２において、ＣＰＵ２０２は、水用高周波磁場の印加によって被検体から放出される磁気共鳴信号を受信する。ＣＰＵ２０２は、上述のＦＦＴ部２５１（図８参照）として、受信した時間領域の磁気共鳴信号を周波数領域の周波数スペクトルに変換する。
【００７９】
ステップＳ３４において、ＣＰＵ２０２は、上述の抽出部２５３（図８参照）として、周波数スペクトルの中から信号強度が相対的に高く、水の共鳴周波数に相当する周波数を抽出する。ＣＰＵ２０２は、上述の算出部２５５（図８参照）として、過去に得られた水の共鳴周波数と、現在の水の共鳴周波数との間の差を算出する。当該差が水の共鳴周波数の化学シフトに相当する。
【００８０】
ステップＳ４０において、ＣＰＵ２０２は、上述の差分部２５７（図８参照）として、ステップＳ３４で算出された水の化学シフトから、ステップＳ２４で算出されたリン化合物の化学シフトを差分する。
【００８１】
ステップＳ５０において、ＣＰＵ２０２は、上述の検出部２５８（図８参照）として、水の共鳴周波数の変化量と温度との間の相関関係（たとえば、「０．０１ｐｐｍ／℃」）を利用して、ステップＳ４０で算出された差に応じた温度変化を検出する。典型的には、ＣＰＵ２０２は、当該差を相関係数「０．０１ｐｐｍ／℃」で除算する。
【００８２】
ステップＳ６０において、ＣＰＵ２０２は、本実施の形態に従う温度計測処理を終了するか否かを判断する。たとえば、ＣＰＵ２０２は、検査対象である全領域に対する温度計測処理を終えたときに本温度計測処理を終了する。ＣＰＵ２０２は、検査対象である全領域に対する温度計測処理を終えたと判断した場合（ステップＳ６０においてＹＥＳ）、本温度計測処理を終了すると判断する。そうでない場合には（ステップＳ６０においてＮＯ）、ＣＰＵ２０２は、制御をステップＳ２０に戻す。
【００８３】
なお、水用高周波磁場およびリン用高周波磁場を印可する順番は、任意である。すなわち、図１２に示されるように、リン用高周波磁場が印可された後に、水用高周波磁場が印可されてもよいし、図７に示されるように、水用高周波磁場が印可された後に、リン用高周波磁場が印可されてもよい。
【００８４】
［計算機システム２００のハードウェア構成］
図１３を参照して、第１の実施の形態に従う計算機システム２００のハードウェア構成の一例について説明する。図１３は、計算機システム２００の主要なハードウェア構成を示すブロック図である。図１３に示されるように、計算機システム２００は、ＲＯＭ（Read  Only  Memory）２０１と、ＣＰＵ２０２と、ＲＡＭ（Random  Access  Memory）２０３と、ネットワークＩ／Ｆ（インタフェース）２０４と、表示部２０５と、記憶装置２１０とを含む。
【００８５】
ＲＯＭ２０１は、オペレーティングシステム（ＯＳ：Operating  System）、計算機システム２００で実行される制御プログラム等を格納する。ＣＰＵ２０２は、オペレーティングシステムや計算機システム２００の制御プログラム等の各種プログラムを実行することで、計算機システム２００の動作を制御する。ＲＡＭ２０３は、ワーキングメモリとして機能し、プログラムの実行に必要な各種データを一時的に格納する。
【００８６】
ネットワークＩ／Ｆ２０４は、アンテナ（図示しない）を介して、他の通信機器との間でデータを送受信する。他の通信機器は、たとえば、サーバ、その他の通信機能を有する装置等を含む。計算機システム２００は、アンテナを介して、本実施の形態に従う各種の処理を実現するためのプログラムをダウンロードできるように構成されてもよい。
【００８７】
表示部２０５は、たとえば、液晶ディスプレイ、有機ＥＬ（Electro  Luminescence）ディスプレイ、またはその他の表示機器等を含む。表示部２０５は、タッチセンサ（図示しない）と組み合わされてタッチパネルとして構成されてもよい。これにより、診断医は、ＭＲＩ１００に対する操作をタッチパネル上で行なうことが可能になる。
【００８８】
記憶装置２１０は、たとえば、ハードディスクや外付けの記憶装置等の記憶媒体である。一例として、記憶装置２１０は、本実施の形態に従う各種の処理を実現するための温度計測プログラム２１２とを保持する。
【００８９】
温度計測プログラム２１２は、単体のプログラムとしてではなく、任意のプログラムの一部に組み込まれて提供されてもよい。この場合、任意のプログラムと協働して本実施の形態に従う温度計測処理が実現される。このような一部のモジュールを含まないプログラムであっても、本実施の形態に従う計算機システム２００の趣旨を逸脱するものではない。さらに、本実施の形態に従う温度計測プログラム２１２によって提供される機能の一部または全部は、専用のハードウェアによって実現されてもよい。さらに、計算機システム２００とサーバとが協働して、本実施の形態に従う温度計測処理を実現するようにしてもよい。さらに、少なくとも１つのサーバが本実施の形態に従う温度計測処理を実現する、所謂クラウドサービスの形態で計算機システム２００が構成されてもよい。
【００９０】
［小括］
上述の実験によって、発明者らは、温度が変化する環境下ではリン化合物の共鳴周波数がほぼ一定であることを発見した。また、発明者らは、磁場が変化する環境下ではリン化合物の共鳴周波数がＮＡＡと同じように変化することを発見した。一方で、水の共鳴周波数は、温度変化と磁場変化との両方の影響を受けやすい。そのため、リン化合物の共鳴周波数の変化量を基準として、水の共鳴周波数の変化量に着目すれば、磁場変化の影響を除外することができる。
【００９１】
本実施の形態に従う温度計測方法は、このような原理に基づいて、水の共鳴周波数の変化量と、リン化合物の共鳴周波数の変化量との差を算出する。そして、当該温度計測方法は、水の共鳴周波数の変化量と温度との間の相関関係を利用して、当該差に応じた温度変化を検出する。これにより、被検体内においてＮＡＡが存在しない場所であっても、リン化合物が存在する場所であれば、正確に温度が計測され得る。
【００９２】
＜第２の実施の形態＞
［概要］
第１の実施の形態では、リン化合物において新たに発見された特性を温度計測処理に応用していた。これに対して、第２の実施の形態では、リン化合物の当該特性を磁場計測処理に応用する。
【００９３】
上述したように、リン化合物の共鳴周波数は、磁場の変化に応じて変化する。一方で、リン化合物の共鳴周波数は、温度に依存せずほぼ一定である。そのため、リン化合物の共鳴周波数の変化量は、磁場の強度を示す指標となり得る。特に、温度が変化する環境下では、リン化合物の共鳴周波数の変化量が有効な指標になり得る。第２の実施の形態に従う磁場計測方法は、リン化合物の共鳴周波数における温度不変性と磁場可変性とを利用して、磁場の変化量を計測する。
【００９４】
本磁場計測方法は、図４および図５の実験結果を応用したものであり、磁場計測装置等に応用され得る。当該磁場計測装置は、たとえば、ＭＲＩ１００等の磁場発生装置の点検時に利用される。一例として、ＭＲＩ１００の静磁場磁石１（図６参照）が点検対象となる。上述したように、静磁場磁石１としては、超伝導磁石や永久磁石等が用いられる。超伝導磁石や永久磁石によって印加される静磁場は、時間の経過とともに劣化する。上記磁場計測装置は、このような劣化を点検することができ、その点検結果に基づいて、磁場発生装置が正常に動作しているか否かを判断することができる。当該点検においては、たとえば、リン化合物を含んだ水溶液や人体に存在するリン化合物が利用される。リン化合物を含んだ水溶液が利用される場合には、静磁場磁石１の安定性が点検され得る。人体に存在するリン化合物が利用される場合には、ＭＲＩ１００の動作時における性能が評価され得る。
【００９５】
また、磁場計測装置は、磁場発生装置の使用時に異常が発生していることを検出するために応用され得る。これにより、磁場発生装置に異常が発生したことを早期に発見することができ、磁場計測装置を停止する等の対処を早急に行なうことができる。
【００９６】
［磁場計測方法］
図１４を参照して、リン化合物の上記特性を利用した磁場計測方法について説明する。図１４は、第２の実施の形態に従う計算機システム２００による磁場計測方法を概略的に示した概念図である。計算機システム２００は、機能構成として、ＦＦＴ部２８２と、抽出部２８４と、算出部２８６と、検出部２８８と、特定部２９０とを含む。
【００９７】
ＦＦＴ部２８２は、リン用高周波磁場の印加によって被検体から検出される磁気共鳴信号ｆＰ（ｔ）を周波数スペクトルＦＰ（ω）に変換する。すなわち、ＦＦＴ部２８２は、時間領域の磁気共鳴信号ｆＰ（ｔ）を、周波数領域の周波数スペクトルＦＰ（ω）に変換する。一例として、当該変換処理には、高速フーリエ変換が用いられる。
【００９８】
抽出部２８４は、周波数スペクトルＦＰ（ω）の中から信号強度が相対的に高く、リン化合物の共鳴周波数に相当する周波数ｆ１Ｐを抽出する。より具体的には、抽出部２８４は、リン化合物の共鳴周波数を中心とする所定範囲のスペクトルの中から最大となる信号強度を選択し、選択された信号強度に対応する周波数を、リン化合物の共鳴周波数に相当する周波数ｆ１Ｐとして抽出する。
【００９９】
算出部２８６は、リン化合物の時間的な変化量ΔｆＰ（すなわち、化学シフト）を算出する。たとえば、算出部２８６は、周波数ｆ１Ｐよりも前に取得されたリン化合物の共鳴周波数である周波数ｆ０Ｐを周波数ｆ１Ｐから差分する。
【０１００】
検出部２８８は、リン化合物の共鳴周波数の変化量と磁場強度との間の相関関係を利用して、変化量ΔｆＰに応じた磁場変化を検出する。
【０１０１】
特定部２９０は、共鳴周波数ｆＰ（ｔ）が放出された被検体の位置を特定する。より具体的には、傾斜磁場コイル２（図６参照）は、上述したように、スライス選択用傾斜磁場Ｇｓ、位相エンコード用傾斜磁場Ｇｅおよびリードアウト用傾斜磁場Ｇｒを発生する。特定部２９０は、スライス選択用傾斜磁場Ｇｓ、位相エンコード用傾斜磁場Ｇｅ、およびリードアウト用傾斜磁場Ｇｒによって、ｋ空間におけるＰＥ方向、ＲＯ方向、ＳＥ方向の位置を示す成分をエンコードする。特定部２９０は、エンコード結果に基づいて、共鳴周波数ｆＰ（ｔ）が放出された被検体の位置を特定する。当該位置情報は、たとえば、３次元の座標（ｘ，ｙ，ｚ）で示される。
【０１０２】
特定部２９０により、被検体の単位領域ごとにリン化合物の共鳴周波数を特定することができる。これにより、検出部２８８は、被検体の単位領域ごとに磁場変化を検出することができる。
【０１０３】
［計算機システム２００の制御構造］
図１５を参照して、第２の実施の形態に従う計算機システム２００の制御構造について説明する。図１５は、第２の実施の形態に従う計算機システム２００が実行する処理の一部を表わすフローチャートである。図１５の処理は、計算機システム２００のＣＰＵ２０２が温度計測プログラムを実行することにより実現される。他の局面において、処理の一部または全部が、回路素子その他のハードウェアによって実行されてもよい。
【０１０４】
ステップＳ７０において、ＣＰＵ２０２は、リン化合物の共鳴周波数を含む高周波磁場（すなわち、リン用高周波磁場）を被検体に印加するための制御信号を送信部３（図６参照）に送る。これにより、リン用高周波磁場が被検体に印加される。
【０１０５】
ステップＳ７２において、ＣＰＵ２０２は、リン用高周波磁場の印加によって被検体から放出される磁気共鳴信号を受信する。ＣＰＵ２０２は、上述のＦＦＴ部２８２（図１４参照）として、受信した時間領域の磁気共鳴信号を周波数領域の周波数スペクトルに変換する。
【０１０６】
ステップＳ７４において、ＣＰＵ２０２は、上述の抽出部２８４（図１４参照）として、周波数スペクトルの中から信号強度が相対的に高く、リン化合物の共鳴周波数に相当する周波数を抽出する。ＣＰＵ２０２は、上述の算出部２８６（図１４参照）として、過去に得られたリン化合物の共鳴周波数と、現在のリン化合物の共鳴周波数との間の差を算出する。当該差がリン化合物の共鳴周波数の化学シフトに相当する。
【０１０７】
ステップＳ８０において、ＣＰＵ２０２は、上述の検出部２８８（図１４参照）として、リン化合物の共鳴周波数の変化量と磁場強度との間の相関関係を利用して、周波数の変化量に応じた磁場変化を検出する。
【０１０８】
ステップＳ９０において、ＣＰＵ２０２は、本実施の形態に従う磁場計測処理を終了するか否かを判断する。たとえば、ＣＰＵ２０２は、検査対象である全領域に対する磁場計測処理を終えたときに本磁場計測処理を終了する。ＣＰＵ２０２は、検査対象である全領域に対する磁場計測処理を終えたと判断した場合（ステップＳ８０においてＹＥＳ）、本磁場計測処理を終了すると判断する。そうでない場合には（ステップＳ８０においてＮＯ）、ＣＰＵ２０２は、制御をステップＳ７０に戻す。
【０１０９】
［小括］
以上のようにして、第２の実施の形態では、リン化合物の共鳴周波数における温度不変性と磁場可変性とを利用して、磁場の変化量を計測する。第２の実施の形態における磁場計測方法は、温度変化に影響されないので、正確に磁場を計測することが可能になる。
【０１１０】
今回開示された実施の形態は全ての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は上記した説明ではなくて特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内での全ての変更が含まれることが意図される。
【符号の説明】
【０１１１】
１  静磁場磁石、２  傾斜磁場コイル、３  送信部、４  送信コイル、５  受信コイル、６  受信部、１０  電源、１３  シーケンス制御部、２０  寝台、２１  天板、５１，５３，５５，５７，５９，６１，６３，６４，７１  グラフ、６４Ａ  最大スペクトル、６５，６７，７３，７５  温度変化、１００  ＭＲＩ、２００  計算機システム、２０１  ＲＯＭ、２０２  ＣＰＵ、２０３  ＲＡＭ、２０４  ネットワークＩ／Ｆ、２０５  表示部、２１０  記憶装置、２１２  温度計測プログラム、２２０Ａ，２２０Ｂ  温度分布、２２２  指標画像、２５１，２５２，２８２  ＦＦＴ部、２５３，２５４，２８４  抽出部、２５５，２５６，２８６  算出部、２５７  差分部、２５８，２８８  検出部、２６０，２９０  特定部。
【図１】
2016168087.tif 000003
【図２】
2016168087.tif 000004
【図３】
2016168087.tif 000005
【図４】
2016168087.tif 000006
【図５】
2016168087.tif 000007
【図６】
2016168087.tif 000008
【図７】
2016168087.tif 000009
【図８】
2016168087.tif 000010
【図１０】
2016168087.tif 000011
【図１３】
2016168087.tif 000012
【図１４】
2016168087.tif 000013
【図９】
2016168087.tif 000014
【図１１】
2016168087.tif 000015
【図１２】
2016168087.tif 000016
【図１５】
2016168087.tif 000017