"磁気粒子を含む撮像対象を画像化する磁気粒子イメージング装置であって、前記撮像対象に磁場を付与する磁場付与手段と、前記磁場の付与方向を制御する付与方向制御手段と、前記撮像対象からの磁場変化を検知する検知手段と、前記検知結果に基づいて生成した投影データを処理する処理手段とを備え、前記付与方向制御手段は、前記撮像対象に対して付与する前記磁場の並進制御と前記磁場の付与角度の制御とを行う、磁気粒子イメージング装置。前記処理手段は、フィルタ補正逆投影法を適用することによって前記投影データを処理する、請求項1に記載の磁気粒子イメージング装置。前記処理手段は、最尤推定-期待値最大化法を適用することによって前記投影データを処理する、請求項1に記載の磁気粒子イメージング装置。前記処理手段は、順位付部分集合-期待値最大化法を適用することによって前記投影データを処理する、請求項1に記載の磁気粒子イメージング装置。前記撮像対象は、生体物質を含む、請求項1~請求項4のうちの一項に記載の磁気粒子イメージング装置。前記撮像対象は、工業製品を含む、請求項1~請求項4のうちの一項に記載の磁気粒子イメージング装置。磁気粒子を含む撮像対象を画像化する磁気粒子イメージング方法であって、前記撮像対象に磁場を付与する磁場付与ステップと、前記磁場の付与方向を制御する付与方向制御ステップと、前記撮像対象からの磁場変化を検知する検知ステップと、前記検知結果に基づいて生成した投影データを処理ステップとを包含し、前記付与方向制御ステップは、前記撮像対象に対して付与する前記磁場の並進制御と前記磁場の付与角度の制御とを実行する、磁気粒子イメージング方法。前記処理ステップは、フィルタ補正逆投影法を適用することによって前記投影データを処理する、請求項7に記載の磁気粒子イメージング方法。前記処理ステップは、最尤推定-期待値最大化法を適用することによって前記投影データを処理する、請求項7に記載の磁気粒子イメージング方法。前記処理ステップは、順位付部分集合-期待値最大化法を適用することによって前記投影データを処理する、請求項7に記載の磁気粒子イメージング方法。ところが、従来の方法では、FFPを用いるための空間分解能や感度がFFPの軌道に大きく依存し、その結果、空間分解能や感度の空間的均一性が低下する問題がある。本発明は、上記問題点に鑑みてなされたものであり、MPIにおいてFFP符号化スキームとは異なるスキームに画質改善手法を適応することで、空間分解能や感度の空間的均一性を向上し得る磁気粒子イメージング装置及び磁気粒子イメージング方法を提供することを目的とする。上記課題を解決するための本発明に係る磁気粒子イメージング装置は、磁気粒子を含む撮像対象を画像化する磁気粒子イメージング装置であって、前記撮像対象に磁場を付与する磁場付与手段と、前記磁場の付与方向を制御する付与方向制御手段と、前記撮像対象からの磁場変化を検知する検知手段と、前記検知結果に基づいて生成した投影データを処理する処理手段とを備え、前記付与方向制御手段は、前記撮像対象に対して付与する前記磁場の並進制御と前記磁場の付与角度の制御とを行う。本発明に係る磁気粒子イメージング装置によれば、前記処理手段は、フィルタ補正逆投影法を適用することによって前記投影データを処理し得る。本発明に係る磁気粒子イメージング装置によれば、前記処理手段は、最尤推定-期待値最大化法を適用することによって前記投影データを処理し得る。本発明に係る磁気粒子イメージング装置によれば、前記処理手段は、順位付部分集合-期待値最大化法を適用することによって前記投影データを処理し得る。本発明に係る磁気粒子イメージング装置によれば、前記撮像対象は、生体物質を含み得る。本発明に係る磁気粒子イメージング装置によれば、前記撮像対象は、工業製品を含み得る。上記課題を解決するための本発明に係る磁気粒子イメージング方法は、磁気粒子を含む撮像対象を画像化する磁気粒子イメージング方法であって、前記撮像対象に磁場を付与する磁場付与ステップと、前記磁場の付与方向を制御する付与方向制御ステップと、前記撮像対象からの磁場変化を検知する検知ステップと、前記検知結果に基づいて生成した投影データを処理ステップとを包含し、前記付与方向制御ステップは、前記撮像対象に対して付与する前記磁場の並進制御と前記磁場の付与角度の制御とを実行する。本発明に係る磁気粒子イメージング方法によれば、前記処理ステップは、フィルタ補正逆投影法を適用することによって前記投影データを処理し得る。本発明に係る磁気粒子イメージング方法によれば、前記処理ステップは、最尤推定-期待値最大化法を適用することによって前記投影データを処理し得る。本発明に係る磁気粒子イメージング方法によれば、前記処理ステップは、順位付部分集合-期待値最大化法を適用することによって前記投影データを処理し得る。本発明の磁気粒子イメージング装置は、3組以上のマックスウェルコイルを用いて磁場が0になる線フィールドフリーライン:FFLを生成し、FFLを並進移動及び回転移動させることによってFFL上に存在するMNPの画像化を行う。そして更にFFL上に存在するMNPの画像データに対して画質改善手法を適応することで、空間分解能や感度の空間的均一性を向上し得る。[実施の形態]図1は、本発明に係る磁気粒子イメージング装置100の模式図である。磁気粒子イメージング装置100は、磁気粒子を含む撮像対象Aを画像化する。磁気粒子イメージング装置100は、磁場制御手段110と、検知手段120と、処理手段130とを備える。磁場制御手段110は、撮像対象Aに磁場を付与する磁場付与手段112と、磁場の付与方向を制御する付与方向制御手段114とを備える。付与方向制御手段114は、画像再構成法により処理するための検知結果を得るように、磁場の付与方向を制御する。付与方向制御手段114は、撮像対象Aに対して付与する磁場の並進制御を行う並進制御手段114aと磁場の付与角度の制御を行う角度制御手段114bとを備える。検知手段120は、撮像対象Aからの磁場変化を検知する。検知手段120は、例えば、2組のコイルである。また、処理手段130は、検知手段120によって検知した検知結果に基づいて生成した投影データを処理する。処理手段130は、例えば、検知データを演算するために適した演算機PCである。さらに、投影データは、処理手段130がFFLに沿った方向に検知データを投影処理することによって得たデータである。磁気粒子イメージング装置100において、例えば、並進制御手段114aは2組のヘルムホルツコイルであり、角度制御手段114bは、4組のマックスウェルコイルである。ここで、並進制御手段114aは、撮像対象Aに磁場を付与する磁場付与手段112としても機能する。磁気粒子イメージング装置100において、磁場付与手段112によって撮像対象Aに磁場を付与しつつ、並進制御手段114aと角度制御手段114bとによって勾配磁場を制御し、FFLを生成する。例えば、磁場付与手段112並進制御手段114aは、撮像対象Aに交流磁場を付与する。並進制御手段114aは、交流磁場の付与状態を変更することで、FFLの並進制御を行う。例えば、角度制御手段114bは、撮像対象Aに直流磁場を付与する。角度制御手段114bは、直流磁場の付与状態を変更することで、FFLの角度の制御を行う。磁気粒子イメージング装置100によれば、並進制御手段114aと角度制御手段114bとを協働させることで勾配磁場を制御し、FFLを生成し、さらにFFLの並進回転を制御する。以下、図2を参照して、本発明に係る検知結果処理手法を説明する。図2は、本発明に係る検知結果処理手法を説明するための座標図である。Pr,θは、原点からFFLまでの距離がrであり、FFLの角度がθである場合の投影データを示す。Piは投影ビンiでの投影データを示す。処理手段130が、検知結果に基づいて生成した投影データを処理する。FFLは付与方向制御手段114並進制御手段114a2組のヘルムホルツコイルと角度制御手段114b4組のマックスウェルコイルにより生成し得、印加する電流を変えることにより、コイルを空間で静止させつつ、FFLを任意に回転、並進する。位置rにおける角度θ[Prθ]での投影データは、ただし、Unr,θは、図2に示すように、原点からFFLまでの距離がrであり、FFLの角度がθである場合に与えられるMPI信号のn番目の周波数成分であるにより生成される。本発明による画像再構成法として、フィルタ補正逆投影法、最尤推定-期待値最大化法又は順位付部分集合-期待値最大化法OS-EM法を適用する。処理手段130が、検知結果に基づいて生成した投影データをこれらの手法に適用する。処理手段130が、フィルタ補正逆投影法を適用することによって投影データを処理する場合には、x,y[cx,y]でのMNPの濃度が式1により計算される。フィルタ補正逆投影法によれば、投影データを補正した後に逆投影を行うため、画像は元データに近い状態で再構成される。処理手段130が、最尤推定-期待値最大化法を適用することによって投影データを処理する場合には、位置jでのk1反復ck1j後のMNPの濃度が式2ただし、iは投影ビン図2、Piはビンiにおける投影データ図2、Pijは画素jにおいてMPI信号が投影ビンiに寄与する確率により与えられる。最尤推定-期待値最大化法は、測定投影データから確率的に最も期待し得る最も可能性の高い再構成画像を逐次近似で推定する。最尤推定-期待値最大化法には、いくつかの魅力的な特長がある。復元画像が非負non-negativeである。全画像濃度は各反復において保存されている。収束は論理に仮定されている。復元画像は最大可能性推定に収束する。なお、本発明では、投影ビンiおよび画素jの両方が同一FFLに位置している場合は、式2におけるPijを合同と仮定している。同一FFLに位置していない場合は、Pijはゼロと仮定し得る。式2のPijに受信コイルの位置依存感度の効果を含めることで、復元画像の精度をさらに改良し得る。処理手段130は、順位付部分集合-期待値最大化法OS-EM法を適用することによっても投影データを処理することができる。OS-EM法orderedsubsets-expectationmaximizationは、最尤推定-期待値最大化法を投影データを部分集合に分割して計算に使う手法によって高速化した手法で、計算時間が早く臨床に使用できる。OS-EM法は、式3により計算される。ここで、Smはm番目の部分集合を表し、上記の計算を全ての部分集合について計算する。すると、計算時間はML-EM法と同じであるが、更新回数は部分集合の数だけ多くなり、それだけ実質的に収束の速度が上昇することになる。以上、図1と図2とを参照して、本発明の磁気粒子イメージング装置100を用いた磁気粒子イメージングの手法を説明した。なお、本発明の磁気粒子イメージング装置100によれば、付与方向制御手段114が磁場の付与方向を制御し得る限りは、ヘルムホルツコイルとマックスウェルコイルとの組み合わせに限定されない。さらに、例えば、並進制御手段114aは2組のヘルムホルツコイルで角度制御手段114bは4組のマックスウェルコイルであったが、並進制御手段114aが撮像対象Aに対して付与する磁場の並進制御を行い得る限りはヘルムホルツコイルの組の数は任意であり、角度制御手段114bが磁場の付与角度の制御を行い得る限りはマックスウェルコイルの組の数は任意である。以下、図3~図8を参照して、本発明の磁気粒子イメージング手法の妥当性と有効性とを説明する。本発明に係る磁気粒子イメージング装置100によれば、磁気粒子を含み得る限りは所望の撮像対象に適応し得る。例えば、生体物質である。磁気粒子イメージング装置100によれば、体外からの電磁波照射により生じる非線形応答を基に、生体物質に投与した磁性ナノ粒子からの微弱な磁気信号を検出する。磁気粒子イメージング装置100は、高感度、高空間分解能および速い描画速度で超常磁性酸化鉄SPIOのような磁気ナノ粒子MNPの空間分布を検知する。しかし、磁気粒子を含み得る限りは、撮像対象は、工業製品や、その他新規な材料であり得る。磁気粒子イメージング装置100によれば、電磁波照射により生じる非線形応答を基に、例えばガス管工業製品の一例の内部に投与した磁性ナノ粒子からの微弱な磁気信号を検出する。磁気粒子イメージング装置100は、ガス管内部の亀裂に侵入した磁気ナノ粒子MNPの空間分布を検知し、ガス管内部の亀裂を非破壊で認識し得る。以下、撮像対象Aとしてモデルを使用して本発明の磁気粒子イメージング手法の妥当性と有効性とを説明する。図3は、シミュレーション研究の中で使用されるモデル撮像対象Aを示す。モデルは、128×128画素を有する128×128/mm2の大きさである。図中の数字は、穴の直径mm単位であり、磁気ナノ粒子MNPが充填されていると仮定している。図4は、FFLを回転、並進した例を示す。図4aは、0、30、60、90および120度で回転したFFLの例であり、図4bは、60度回転したFFLを0、20、40、-20及び-40ピクセルで各々並進したFFLの例を示す図5は、磁界の勾配変化に対応する再構成画像の比較を示す図である。図5aはフィルタ補正逆投影法Shepp-Loganフィルタを使用を適用した再構成画像を示し、図5bは最尤推定-期待値最大化法を適用した再構成画像を示す。図中の数字は、磁界の勾配単位:T/mを示す。SN比を10、投影データの数を32、駆動磁界の強度を5mT、MNPの直径を30nmとした。図6は、MNPの直径変化に対応する再構成画像の比較を示す図である。図6aはフィルタ補正逆投影法を適用した再構成画像を示し、図6bは最尤推定-期待値最大化法を適用した再構成画像を示す。図中の数字は、MNPの直径単位:nmを示す。SN比を10、投影データの数を32、駆動磁界の強度を5mT、磁界の勾配を5T/mとした。図7は、投影データ数の変化に対応する再構成画像の比較を示す図である。図7aはフィルタ補正逆投影法を適用した再構成画像を示し、図7bは最尤推定-期待値最大化法を適用した再構成画像を示す。図中の数字は、投影データの数を示す。SN比を10、駆動磁界の強度を5mT、選択磁界の勾配を5T/m、MNPの直径を30nmとした。図8は、駆動磁界の強度の変化に対応する再構成画像の比較を示す図である。図8aはフィルタ補正逆投影法を適用した再構成画像を示し、図8bは最尤推定-期待値最大化法を適用した再構成画像を示す。図中の数字は、駆動磁界の強度単位:mTを示す。SN比を10、投影データの数を32、磁界の勾配を5mT、MNPの直径を30nmとした。以上、図3~図8を参照して本発明の磁気粒子イメージング手法の妥当性と有効性とを説明した。本発明に係る磁気粒子イメージング装置100によれば、画質の向上やデータ収集時間の短縮を実現し得る。例えば、現在臨床に広く普及している核磁気共鳴撮像装置MRIに比較して、感度が103~106大きく、また、コストも大幅に低減できるため、本発明の実用化によって臨床にも広く普及し得、経済効果も大きいものと期待し得る。"