(57)【特許請求の範囲】
【請求項１】
タギングＭＲ画像における撮像対象の運動解析方法であって、
撮像対象の同じ運動に対して、異なるタグパターンよるＮL回（ＮLは２以上の正の整数）のｃｉｎｅ撮像をすることで、前記撮像対象の前記運動における複数の時相について撮像された時系列のタギングＭＲ画像をＮL個得る工程と、
各時相のＮL個のタギングＭＲ画像における同一画素についてのＮL個の画素値を、当該画素についての長さＮLの画素値系列とする工程と、
異なる時相において、前記画素値系列が同一の符号系列を構成する画素を検出することで、時系列のタギングＭＲ画像における撮像対象の運動を解析する工程と、
を含み、
前記タグパターンそれぞれは、前記画素値系列が、所定の符号系列を構成するようにパターンが形成されたものである、
ことを特徴とするタギングＭＲ画像における撮像対象の運動解析方法。
【請求項２】
撮像対象の運動を解析する前記工程では、前記画素値系列に基づいて、各画素において混在する複数の符号系列それぞれが占める割合を算出し、当該割合に基づいて、１画素の大きさ未満での撮像対象の運動量を算出する
請求項１記載の運動解析方法。
【請求項３】
前記符号系列は、直交符号系列である請求項１記載のタギングＭＲ画像における撮像対象の運動解析方法。
【請求項４】
前記符号系列は、拡散符号系列である請求項１記載のタギングＭＲ画像における撮像対象の運動解析方法。
【請求項５】
前記符号系列は、ＰＮ系列である請求項１記載のタギングＭＲ画像における撮像対象の運動解析方法。
【請求項６】
前記符号系列は、Ｍ系列である請求項１記載のタギングＭＲ画像における撮像対象の運動解析方法。
【請求項７】
タギングＭＲ画像における撮像対象の運動解析方法であって、
撮像対象に対し、所定のタグパターンでｃｉｎｅ撮像することで、前記撮像対象の運動における複数の時相について撮像された時系列のタギングＭＲ画像を得る工程と、
各時相のタギングＭＲ画像に含まれるＮL個（ＮLは２以上の正の整数）の画素それぞれの画素値を、当該ＮL個の画素からなる領域についての長さＮL個の画素値系列とする工程と、
異なる時相において、前記画素値系列が同一の符号系列を構成する領域を検出することで、時系列のタギングＭＲ画像における撮像対象の運動を解析する工程と、を含み、
前記所定のタグパターンは、前記画素値系列が、所定の符号系列を構成するようにパターンが形成されたものである、
ことを特徴とするタギングＭＲ画像における撮像対象の運動解析方法。
【請求項８】
撮像対象の運動を解析する前記工程では、各画素の画素値の大きさに基づいて、タグパターンの明領域又は暗領域が各画素の範囲内において本来占めるべき割合を算出し、当該割合に基づいて１画素の大きさ未満での撮像対象の運動量を算出する
請求項７記載の運動解析方法。
【請求項９】
ＭＲＩ装置であって、
撮像対象の同じ運動に対して、異なるタグパターンよるＮL回（ＮLは２以上の正数）のｃｉｎｅ撮像をすることで、前記撮像対象の前記運動における複数の時相について撮像された時系列のタギングＭＲ画像をＮL個得る手段と、
特定の時相のＮL個のタギングＭＲ画像における同一画素についてのＮL個の画素値からなる画素値系列を、前記特定の時相において当該画素を示す長さＮLの符号系列とみなし、前記符号系列と同じ符号系列となる画素値系列を持つ画素を、前記特定の時相以外の時相のＮL個のタギングＭＲ画像において検出することで、時系列のタギングＭＲ画像における撮像対象の運動を解析する手段と、
を備え、
前記タグパターンそれぞれは、前記画素値系列が、所定の符号系列を構成するようにパターンが形成されたものである、
ことを特徴とするＭＲＩ装置。
【請求項１０】
撮像対象の運動を解析する前記手段では、前記画素値系列に基づいて、各画素において混在する複数の符号系列それぞれが占める割合を算出し、当該割合に基づいて、１画素の大きさ未満での撮像対象の運動量を算出する
請求項９記載のＭＲＩ装置
【請求項１１】
前記符号系列は、直交符号系列である請求項９記載のＭＲＩ装置。
【請求項１２】
前記符号系列は、拡散符号系列である請求項９記載のＭＲＩ装置。
【請求項１３】
前記符号系列は、ＰＮ系列である請求項９記載のＭＲＩ装置。
【請求項１４】
前記符号系列は、Ｍ系列である請求項９記載のＭＲＩ装置。
【請求項１５】
ＭＲＩ装置であって、
撮像対象に対し、所定のタグパターンでｃｉｎｅ撮像することで、前記撮像対象の運動における複数の時相について撮像された時系列のタギングＭＲ画像を得る手段と、
特定の時相のタギングＭＲ画像に含まれるＮL個（ＮLは２以上の正の整数）の画素それぞれの画素値からなる画素値系列を、当該ＮL個の画素からなる領域を示す長さＮLの符号系列とみなし、前記符号系列と同じ符号系列となる画素値系列を持つ領域を、前記特定の時相以外の時相のタギングＭＲ画像において検出することで、時系列のタギングＭＲ画像における撮像対象の運動を解析する手段と、
を備え、
前記所定のタグパターンは、前記画素値系列が、所定の符号系列を構成するようにパターンが形成されたものである、
ことを特徴とするＭＲＩ装置。
【請求項１６】
撮像対象の運動を解析する前記手段では、各画素の画素値の大きさに基づいて、タグパターンの明領域又は暗領域が各画素の範囲内において本来占めるべき割合を算出し、当該割合に基づいて１画素の大きさ未満での撮像対象の運動量を算出する
請求項１５記載のＭＲＩ装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、タギングＭＲ画像における撮像対象の運動解析方法及びＭＲＩ装置に関するものである。
【背景技術】
【０００２】
心疾患の診断、治療には、心壁の運動を解析することが非常に重要である。心壁運動の詳細な解析には、人体の組織を磁気的に標識するタギングＭＲＩ法（非特許文献１参照）が用いられることが多い。タギングＭＲＩ法では、時系列タギングＭＲ画像（時系列撮影のＭＲ画像（ｃｉｎｅ  ＭＲＩ））の解析により、心壁の三次元的な運動を解析することが可能である。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】Elias A. Zerhouni, David M. parish, Walter J. Rogers, Andrew Yang,Edward P. Shapiro. Human Heart: Tagging with MR Imaging _ A Method forNoninvasive Assessment of Myocardial Motion. Radiology 169, 59-63, Oct, 1988
【非特許文献２】Albert Montillo, Dimitris Metaxas, Leon Axel. Extracting Tissuedeformation using Gabor filter banks. SPIE Vol.5369, Bellingham, WA, 2004
【発明の概要】
【発明が解決しようとする課題】
【０００４】
タギングＭＲＩ法において撮像対象の運動解析を行うには、タギングＭＲ画像の撮影後の後処理として、撮像対象である心壁領域の抽出、画像処理によるタグ位置の検出を行う必要がある。磁気的な格子状の標識部分はタグと呼ばれるが、時間とともにタグのコントラストが低下する。このため、タグ位置の検出処理には、低ＳＮ（Signal Noise ratio）画像に頑健な方法が必要である。【０００５】
非特許文献２には、ｇａｂｏｒフィルタバンクを利用する方法等が提案されている。この方法はタグの空間的な広がりを利用するため，組織の大域的な運動の解析には有効であるが、局所的なタグの運動解析では解析誤りが生じやすいという問題がある。【０００６】
本発明は、タギングＭＲ画像において、画像処理によってタグを抽出するのとは異なる発想によって、撮像対象の運動を解析することを可能とするための技術を提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明は、タギングＭＲ画像における撮像対象の運動解析方法であって、撮像対象の同じ運動に対して、異なるタグパターンよるＮL回（ＮLは２以上の正の整数）のｃｉｎｅ撮像をすることで、前記撮像対象の前記運動における複数の時相について撮像された時系列のタギングＭＲ画像をＮL個得る工程と、各時相のＮL個のタギングＭＲ画像における同一画素についてのＮL個の画素値を、当該画素についての長さＮLの画素値系列とする工程と、異なる時相において、前記画素値系列が同一の符号系列を構成する画素を検出することで、時系列のタギングＭＲ画像における撮像対象の運動を解析する工程と、を含み、前記タグパターンそれぞれは、前記画素値系列が、所定の符号系列を構成するようにパターンが形成されたものである、ことを特徴とするタギングＭＲ画像における撮像対象の運動解析方法である。上記本発明によれば、符号系列の追跡によって撮像対象の運動を解析することができる。【０００８】
ここで時相とは、撮像対象の運動におけるフェーズをいい、撮像の対象となる運動には、複数の時相（フェーズ）が含まれる。例えば、撮像対象に、移動・変形などの運動を複数回繰り返させることで、同じ運動を複数回ｃｉｎｅ撮像した場合、複数回のｃｉｎｅ撮像で得た複数のｃｉｎｅ画像において、同一の時相（運動フェーズ）にあるものは、撮像対象の同じ状態（位置・形状など）を示すことになる。また、複数のｃｉｎｅ画像において、異なる時相（運動フェーズ）にあるものは、撮像対象の運動によって異なる状態（位置・形状など）を示すことがある。【０００９】
なお、画素には、二次元のピクセルだけでなく３次元のボクセルも含む。また、画素値は、輝度値信号であってもよいし、位相値信号であってもよい。【００１０】
撮像対象の運動を解析する前記工程では、前記画素値系列に基づいて、各画素において混在する複数の符号系列それぞれが占める割合を算出し、当該割合に基づいて、１画素の大きさ未満での撮像対象の運動量を算出するのが好ましい。この場合、１画素の大きさ未満で運動を解析することができる。【００１１】
他の観点からみた本発明は、タギングＭＲ画像における撮像対象の運動解析方法であって、撮像対象に対し、所定のタグパターンでｃｉｎｅ撮像することで、前記撮像対象の運動における複数の時相について撮像された時系列のタギングＭＲ画像を得る工程と、各時相のタギングＭＲ画像に含まれるＮL個（ＮLは２以上の正の整数）の画素それぞれの画素値を、当該ＮL個の画素からなる領域についての長さＮL個の画素値系列とする工程と、異なる時相において、前記画素値系列が同一の符号系列を構成する領域を検出することで、時系列のタギングＭＲ画像における撮像対象の運動を解析する工程と、を含み、前記所定のタグパターンは、前記画素値系列が、所定の符号系列を構成するようにパターンが形成されたものである、ことを特徴とするタギングＭＲ画像における撮像対象の運動解析方法である。上記本発明においても、符号系列の追跡によって撮像対象の運動を解析することができる。【００１２】
前記符号系列は、耐ノイズ性に優れたものが好ましく、例えば，直交符号系列であるのが好ましく、また、拡散符号系列であるのが好ましく、さらには、ＰＮ系列であるのが好ましく、さらには，Ｍ系列であるのが好ましい。【００１３】
また、更に他の観点からみた本発明は、上記工程を実行することができるＭＲＩ装置である。
【発明の効果】
【００１４】
本発明によれば、符号系列の追跡によって撮像対象の運動を解析することができる。
【図面の簡単な説明】
【００１５】
【図１】ＭＲＩ装置の全体構成図である。
【図２】ＭＲＩ装置における処理手順を示すフローチャートである。
【図３】ｃｉｎｅ撮像のフローチャートである。
【図４】ＭＲＩ撮像のためのパルスシーケンス図である。
【図５】複数回のｃｉｎｅ撮像の仕方の説明図である。
【図６】撮像対象の運動解析のフローチャートである。
【図７】複数のｃｉｎｅ撮像から画素値系列を生成する仕方を説明する図である。
【図８】画素の移動追跡を示す図である。
【図９】タグパターンを示す図である。
【図１０】画素への符号系列の割り当て例を示す図である。
【図１１】実験結果を示す画像写真である。
【図１２】撮像対象（ファントム）の運動の仕方の説明図である。
【図１３】初期時相における画素境界と白黒境界を示す図である。
【図１４】第１時相における白黒境界の移動を示す図である。
【図１５】第２時相における白黒境界の移動を示す図である。
【図１６】符号系列相関値と移動量との関係を示す図である。
【図１７】符号系列相関値と移動量との関係を示す図である。
【図１８】符号系列相関値と移動量との関係を示す図である。
【図１９】サブピクセル解析による解析精度とピクセル解析による解析精度の違いを説明する概念図である。
【図２０】本発明の変形例（画像方向への符号化（１次元））を説明する図である。
【図２１】初期時相における画素境界と白黒境界を示す図である。
【図２２】第１時相における白黒境界の移動を示す図である。
【図２３】第２時相における白黒境界の移動を示す図である。
【図２４】白黒領域と移動量との関係を示す図である。
【図２５】変形例を２次元に拡張した例を説明する図である。
【図２６】図２５の例において、符号系列を構成する画像パターンの一覧を示す図である。
【図２７】変形例を２次元に拡張した他の例を説明する図である。
【図２８】図２７の例において、符号系列を構成する画像パターンの一覧を示す図である。
【図２９】図２７の例において、ＡとＡ’を区別するための画像パターンを示す図である。
【発明を実施するための形態】
【００１６】
以下、本発明の好ましい実施形態について添付図面を参照しながら説明する。
［１．ＭＲＩ装置の構成］
本実施形態に係るＭＲＩ装置１は、本発明の位置情報符号化ＭＲＩ法による撮像が行えるものである。本実施形態における位置情報符号化ＭＲＩ法は、撮像対象物（例えば、人体の組織）を磁気的に標識（タギング）するタギングＭＲＩ法を利用するため、前記ＭＲＩ装置１は、タギングＭＲＩ法による撮像が行えるように構成されている。なお、撮像対象を重畳されるタグとしては、任意の方向の縞状又は格子状のものを生成できる。
また、本実施形態では、解析対象の運動として、心壁の運動を例として説明するが、本発明における撮像対象物や解析対象の運動がこれらに限られるものではない。
【００１７】
図１に示すように、このＭＲＩ装置１は、静磁場や傾斜磁場を発生させるための磁場発生部２と、撮像対象物にＲＦパルスを発信するとともに撮像対象物で発生した信号を受信するためのＲＦコイル３と、を備えている。
また、ＭＲＩ装置１は、ＲＦコイル３から発信される所定のＲＦパルスを生成するためのＲＦ送信部４と、ＲＦコイル３にて受信した受信信号（ＭＲＩ信号）を処理するためのＲＦ受信部５と、を備え、ＲＦ送信部４で生成しＲＦコイル３から発信されるＲＦパルスと磁場発生部２で生成される傾斜磁場の組み合わせによりタグパルスシーケンスを生成できる。
【００１８】
磁場発生部２による磁場の発生とＲＦ送信部４によるＲＦパルスの送信は、制御部６によって制御される。
ＲＦ受信部５から出力されたＭＲＩ信号は、Ａ／Ｄ変換部７を介して、コンピュータ８に与えられる。コンピュータ８は、取得したＭＲＩ信号に基づいて、画像処理や後述の運動解析等の処理を行う。
また、コンピュータ８は、ＲＦ送信部４及び制御部６に対して、必要な指令を与えて、これらを制御する。
【００１９】
前記コンピュータ８は、処理プロセッサ９、内部記憶装置及び／又は外部記憶装置を有する記憶部１０、ディスプレイ等からなる表示部１１を備えている。記憶部１０には、ＭＲＩ装置１を制御するためのコンピュータプログラム、撮像対象の運動解析用のコンピュータプログラム、及びその他必要なコンピュータプログラムがインストールされており、これらのプログラムが処理プロセッサによって実行されることで、以下に説明する処理が実現される。
【００２０】
［２．運動解析処理］
図２は、ＭＲＩ装置１の撮像対象（心臓）の運動解析処理の手順を示している。この運動解析処理は、ＭＲＩ装置１が、複数回（ＮL回）のｃｉｎｅ撮像（動画撮像）を行う撮像工程（ステップＳ１）、及び、コンピュータ８が、撮像工程によって得たＮL個のｃｉｎｅＭＲ画像（時系列ＭＲ画像）に基づいて撮像対象の運動の解析処理を行う解析処理工程（ステップＳ２）に大別される。
【００２１】
［２．１．撮像工程］
ステップＳ１の撮像工程では、同一撮像対象の同じ動きを同一撮影範囲で、ＮL回のｃｉｎｅ撮像を繰り返し行う。これらのＮL回の撮像は、それぞれ、異なるＮL種類のタグパターンで行われる。つまり、ステップＳ１の撮像工程によって、撮像対象の同じ動きを示すものの、画像中のタグパターンが互いに異なるＮL個の時系列タギングＭＲ画像（ｃｉｎｅタギングＭＲ画像）が得られる。
【００２２】
ステップＳ１の撮像工程は、具体的には、図３に示すように、１回のｃｉｎｅ撮像（ステップＳ１２〜Ｓ１５）を、ｉが１からＮLになるまで（ＮL回）繰り返すためのループ処理（ステップＳ１１，Ｓ１６，Ｓ１７）を含んでおり、このループ内に、１回ごとのｃｉｎｅ撮像のための処理（ステップＳ１２〜Ｓ１５）が設けられている。
【００２３】
ｃｉｎｅ撮像処理（ステップＳ１２〜Ｓ１５）では、まず、コンピュータ８が、撮像に用いるタグパターンを選択する（ステップＳ１２）。コンピュータ８の記憶部１９には、第１タグパターンから第ＮLタグパターンまでＮL種類のタグパターンが記憶されており、ｉ（ｉ：１〜ＮL）回目のｃｉｎｅ撮像では、第ｉタグパターンが選択される。なお、ＮL種類のタグパターンの詳細については後述する。
【００２４】
そして、コンピュータ８は、選択された第ｉタグパターンのための第ｉタギングシーケンスを構成するＲＦパルスを送信部４によって生成させ、ＲＦコイル３から発信させるとともに、第ｉタギングシーケンスを構成する所定の傾斜磁場を磁場発生部２によって生成する処理を行う（ステップＳ１３）。続いて、コンピュータ８は、第ｉ回目のｃｉｎｅ撮像を行うための処理を実行する（ステップＳ１４）。このｃｉｎｅ撮像によって得られた時系列タギングＭＲ画像は、コンピュータ８の記憶部１０に保存される（ステップＳ１５）。
【００２５】
図４は、１回のｃｉｎｅ撮像処理（ステップＳ１２〜Ｓ１５）を行ったときのパルスシーケンスを示している。このパルスシーケンスは、ＦＬＡＳＨ法によるシーケンス（ＦＬＡＳＨシーケンス）であり、このＦＬＡＳＨシーケンスのプレパレーション部分に、選択されたタグパターンを生成するＤＡＮＴＥタギングシーケンスを付加したものである。
図４のＦＬＡＳＨシーケンスは、１回の静止タギングＭＲ画像を撮像するためのものであり、１回のｃｉｎｅ撮像あたり、図示のＦＬＡＳＨシーケンスがＭ回繰り返され、Ｍ個の静止タギングＭＲ画像からなる時系列タギングＭＲ画像が得られる。ここで、１個の静止タギングＭＲ画像を一つの時相とすると、Ｍ個の静止タギングＭＲ画像からなる時系列タギングＭＲ画像は、Ｍ時相分の時系列タギングＭＲ画像となる。
なお、ＦＬＡＳＨシーケンスは、例えば、エコー時間ＴＥ＝１２ｍｓとし、励起ＲＦパルスにフリップ角α＝３０度のｈｅｒｍｉｔｅパルスを利用することができる。
【００２６】
また、ＤＡＮＴＥタギングシーケンスにおける励起ＲＦパルスには矩形波を利用した。所望のタグパターンを得るには、ＤＡＮＴＥタギングシーケンスにおける矩形波のパルス間隔ｔintを適宜設定するとともに、オフセット周波数ｆ0または各ＲＦパルスに位相オフセットを適宜設定すればよい。ここで、パルス間隔ｔint及びＧxはタグ（画像中の低輝度の標識部分）の間隔ｄintとは反比例する。また、オフセット周波数ｆ0または各ＲＦパルスの位相オフセットは、タグの中心周波数位置を、任意の位置に移動させる。これらの性質を利用して、パルス間隔ｔint及びオフセット周波数ｆ0を適宜設定したタギングシーケンスを設計することで、所望のＮL種類のタグパターンが得られる。
なお、ＭＲＩの撮影シーケンスは、ＦＬＡＳＨシーケンスに限られるものではなく、タグパターンを生成するためのシーケンスもＤＡＮＴＥシーケンスに限られるものではない。
【００２７】
図５（ａ）は、図４に示すような１回分のｃｉｎｅ撮像処理（ステップＳ１２〜Ｓ１５）をＮL回行ったときのパルスシーケンスを示している。ここで、図５（ａ）における第１タギングシーケンスないし第ＮLタギングシーケンスは、それぞれ、図４のＤＡＮＴＥタギングシーケンスに対応する。ただし、第１タギングシーケンスないし第ＮLタギングシーケンスは、それぞれ異なるタグパターンを発生させるためのものである。
また、図５（ａ）における第１ｃｉｎｅ（＃１）から第ＮLｃｉｎｅ（＃ＮL）は、それぞれ、図４のＦＬＡＳＨシーケンスに対応し、第１回目のｃｉｎｅ撮像（＃１）から第ＮL回目のｃｉｎｅ撮像（＃ＮL）をするためのものである。
【００２８】
各回（１回目からＮL回目）のｃｉｎｅ撮像（第１ｃｉｎｅ（＃１）〜第ＮLｃｉｎｅ（＃ＮL））は、それぞれ、撮像対象の同じ運動を撮像するものである。
つまり、撮像対象が心臓である場合、図５（ｂ）に示すように繰り返される心拍において、例えば、１心拍分の時間の間において、１回のｃｉｎｅ撮像を行うものとする。そして、各回のｃｉｎｅ撮像は、心拍周期に関して同一のタイミング（時相）で行われる。
この結果、各回のｃｉｎｅ撮像では、心臓の同じ運動を撮像することになる。したがって、各回（１回目からＮL回目）の時系列タギングＭＲ画像において、画像中の心臓の位置・形状は、同一時相であれば同じ位置・形状となる。
【００２９】
［２．２．解析処理工程］
ステップＳ２の解析処理工程では、コンピュータ８が、ステップＳ１の撮像工程によって得たＮL個の時系列タギングＭＲ画像に基づいて撮像対象の運動の解析処理を行う。
具体的には、図６に示すように、まず、コンピュータ８は、タギング画像の各画素それぞれについて存在するＮL個の画素値を、当該画素についての長さＮLの画素値系列（長さＮLの信号）とする処理を行う（ステップＳ２１）。なお、ここでの画素値は、画素の輝度値であるが、画素の位相値であってもよい。
図７に示すように、ＮL個の時系列タギングＭＲ画像を、時相をそろえて並べた場合、各時相それぞれについて、タギング画像がＮL枚得られたことになる。したがって、タギング画像の個々の画素に着目すると、各時相の各画素について、ＮL個の画素値が存在することになり、これらＮL個の画素値を撮像順に並べて、長さＮLの画素値系列とする。
【００３０】
ここで、タグが存在する低輝度領域（黒い部分）を＋１とし、それ以外の領域（白い部分）を−１とみなすと、長さＮLの画素値系列（長さＮLの信号）は、長さＮLの符号系列で画素値を変調したものと考えることができる。つまり、第１時相における各画素に符号系列が割り当てられていることになる。したがって、撮像対象が運動しない場合、符号系列が特定できればその画素の位置を特定することができる。また、タグは、撮像対象自体に重畳されるものであり、撮像対象の運動に伴って、タグも移動する。したがって、コンピュータ８が、各画素について画素値系列に対応する符号系列を、後述の相関関数の計算によって決定し（ステップＳ２２）、各時相間において同一符号系列を構成する画素の位置を検出・追跡することで（ステップＳ２３）、タグの移動も追跡することができる。したがって、従来のようなタグ位置検出のための画像処理を行わなくても、撮像対象の運動解析が行える。
具体的には、各時相において、同一符号系列を持つ画素の位置を記憶部１０に記憶しておき、また、同一符号系列を持つ画素の移動軌跡を表示部１１に表示することで、装置１のユーザは、撮像対象の運動の様子を的確に把握することができる。
【００３１】
例えば、図８に示すように、第１時相における画素（ｘ1，ｙ1）が示す撮像対象の部位の移動を追跡したい場合を考える。第１時相における画素（ｘ1，ｙ1）の画素値系列に対応する符号系列がＬ1である場合、各時相において符号系列がＬ1となる画素値系列を持つ画素を特定することで、第１時相における画素（ｘ1，ｙ1）が示す撮像対象の部位が、第１時相から第ＮL時相までの各時相においてどのように移動したかを示す移動軌跡を求めることができる。
【００３２】
ここで、図８において、初期時相（第１時相）におけるタグパターンにおいて画素（ｘ1，ｙ1）に付与された位置情報（符号系列Ｌ1）は、撮像対象物の運動により、第ＭX時相においては画素（ｘ2，ｙ2）に存在するものとする。
この場合、初期時相における画素（ｘ1，ｙ1）が符号系列Ｌ1を送信した送信局であり、第ＭX時相における画素（ｘ2，ｙ2）が当該送信局から送信された符号系列Ｌ1を受信した受信局であるとみなすことができる。通信において、受信局は、受信信号（画素値系列）を復号することにより受信信号が示す符号系列を識別することができる。
したがって、本実施形態における上記の移動軌跡の追跡は、符号系列Ｌ1を受信した受信局（画素）を検出することと同様である。
また、タギングＭＲ画像では、時相の経過によって、タグコントラストが低下するため、第ＭX時相における画素の画素値系列は、低ＳＮ比となり、劣悪な通信環境に類似することになる。しかし、通信において、適切な符号系列を選択することで、劣悪な通信環境に対応できるのと同様に、適切な符号系列を選択することで、タギングＭＲ画像におけるタグコントラストの低下にも対応することができる。
【００３３】
さて、画像処理によってタグ領域を抽出する従来の手法では、タグ領域に対応する低輝度領域を検出していたため、正確にはタグ領域の境界のみが運動解析対象であり、それ以外の領域については、内挿による運動の推定を行う必要があった。従来手法で詳細な運動の推定のためには、タグ領域、非タグ領域共に小さくすることが必要となるが、タグ領域、非タグ領域の幅と画素サイズとの関係、あるいは、両領域間の境界が正確な矩形ではないことから、タグ領域、非タグ領域の幅を細くすると、両領域間のコントラストが低下して、検出が困難となる。したがって、心臓の壁運動解析をおこなうためのタギングＭＲＩでは、タグ領域と非タグ領域併せて６ｍｍ程度の間隔で繰り返すパターンを用いるのが一般的となっている。
しかし、心壁の厚さは１０ｍｍ程度であるため、このような間隔ではタグ領域が最大でも２回程度しか現れず、壁の厚み方向における運動の差違に関する詳細な検討は難しかった。
【００３４】
しかしながら、本実施形態の手法を用いると、全ての画素について位置情報を符号化するため、全ての画素に対し、画素サイズでの運動の解析が可能となる。したがって、詳細な運動解析が可能である。
【００３５】
［３．利用可能な符号系列］
前記符号系列として、耐雑音特性、耐混信特定の高いものを用いることで、符号系列が割り当てられた画素の位置を高い精度で特定することができる。例えば、符号系列としては拡散符号を用い、スペクトル拡散（Spread
Spectrum; ＳＳ）を行うことができる。スペクトル拡散とは、狭帯域変調信号を拡散符号によって数百〜数千倍のスペクトルを持つ広帯域信号に変調して伝送し、受信側で元の狭帯域信号に変換した後に復調する通信方式である。ＳＳには耐妨害性、劣悪な伝送路に強いという特長がある。
このように、スペクトラム拡散の劣悪な伝送路に強いという性質を利用することで、時相の経過によるタグコントラストの低下に対し、精度の高い運動解析が可能となる。
【００３６】
ＳＳに利用できる符号系列に必要な条件は二つ挙げられる。一つは，位相差τ＝０ において自己相関が鋭く、τ＝０以外において相関が十分に小さいことである．これはスペクトルを一様に拡散するために必要な条件である。もう一つは、符号系列間の相関が、全ての位相差において十分小さいことである。これは他局間干渉を抑圧するために必要な条件である．これらの条件を満たす系列の一つが，擬似雑音（Pseudo Noise;ＰＮ）系列である。狭義のＰＮ系列の定義は，自己相関値が２レベルの値だけを取る周期系列であることと、１周期中の＋１と−１の個数がたかだか１つだけしか異ならない周期系列であることである。主なＰＮ系列として、Ｍ系列が挙げられる。本実施形態では、符号系列として、このＭ系列を利用する。
なお、上記の拡散符号系列、ＰＮ系列、及びＭ系列は、いずれも符号がほぼ直交する直交性を有しており、直交符号系列でもある。
【００３７】
［４．画素値系列から符号系列を特定するための相関関数の計算］
［４．１．相関関数の計算］
長さＮL、系列数Ｎsの符号系列の集合をＬとする。符号系列集合Ｌのｉ番目の符号系列をＬi（１≦ｉ≦Ｎs）とし、各符号要素をｌij（１≦ｉ≦Ｎs，１≦ｊ≦ＮL）とする。ただし、ｌij∈｛１，−１｝とする。情報符号Ｓ＾∈｛１，−１｝を拡散符号系列Ｌiで拡散した結果をＳ(i)＝｛Ｓj(i)｝（１≦ｊ≦ＮL）とする。このとき、Ｓj(i)＝Ｓ＾×ｌij∈｛１，−１｝となる。
符号系列ＬkとＳj(i)との相関関数Φkは次のように計算される。ただし、本実施形態の手法では、τ＝０しか生じないため、相関関数もτ＝０でのみ計算するものとし、また、簡単のためＳ＾＝１とする。
【００３８】
【数１】
0005419024.tif 000002
【００３９】
このとき、ｉ＝ｋであれば、Φk＝ＮLとなり、ｉ≠ｋであれば、Φk＝−１となる。ゆえに、Φk（ｋ＝１，２，・・・，Ｎs）を計算して最大値を取るｋが、割り当てられた符号系列の番号となる。しかし、実際の通信環境ではＳj(i)のままで伝送を行うことができず、アナログ信号として変調を行う場合が多い。ここではアナログ変調の影響を加味した相関関数を考える。Ａを変調の深さ、Ｂを無変調時の信号強度とすると、変調された受信信号ｓj(i)は、Ｓj(i)＝−１のときＢ、Ｓj(i)＝１のときＢ−Ａとなるので、
【数２】
0005419024.tif 000003
と表すことができる。
【００４０】
このとき変調信号と符号系列Ｌkとの相関関数φkは、
【数３】
0005419024.tif 000004
となる。
【００４１】
ここで、ＬkがＭ系列であれば、
【数４】
0005419024.tif 000005
であるので、
【数５】
0005419024.tif 000006
となる。ただし、Ａ／２＝Ａ’とする。したがって、Ａ，Ｂが一定値を取れば、φk（ｋ＝１，２，・・・，Ｎs）を計算して最小値を取るｋが、割り当てられた符号系列の番号である。つまり、変調信号を利用して符号系列の特定を行うことが可能となる。
【００４２】
［４．２．実施形態における相関関数の計算］
あるＦＯＶ（Ｆｉｅｌｄ  Ｏｆ  Ｖｉｅｗ）について、解像度ＷMTX［ｐｉｘｅｌ］、時相間の間隔ｔph［ｓｅｃ］でＭ時相分のｃｉｎｅ撮像を行うものとする。時刻ｔ＝ｔph×ｐにおける画素（ｘ，ｙ）の画素値（画像信号）をａ（ｘ，ｙ，ｐ）と表す。ただし、初期時相（第１時相）の時刻はｔ＝ｔph×０とするため、０≦ｐ≦Ｍ−１である。
ここで、長さＮLの拡散符号系列Ｌi＝｛ｌij｝を考える。画素値ａ（ｘ，ｙ，ｐ）と拡散符号ｌij（ｊ＝１，２，・・・，ＮL）の積を取ると、スペクトル拡散変調された画素値ｇ(i)（ｘ，ｙ，ｐ，ｊ）が得られる。得られる画像の画素値ｇ(i)（ｘ，ｙ，ｐ，ｊ）を簡単のため、ｍj(i)とおき、ＭＲＩタギング法を用いた変調を考えると、
【数６】
0005419024.tif 000007
と表すことができる。
【００４３】
このとき、ｋ1はタグと非タグ領域の画素値コントラストであり、ｋ2は非タグ領域の画素値に対応する。これらの値は、理想的にはいずれも定数となる。ｌij＝１のとき、ｍj(i)＝ｋ2−ｋ1となり、画像では低輝度になる。また、ｌij＝−１のとき、ｍj(i)＝ｋ2となり、画像では高輝度になる。したがって、タグ領域を適切に設定することで、ＭＲＩの撮影における画素値をスペクトル拡散変調することができる。
【００４４】
このとき、画素値ｍj(i)（１≦ｊ≦ＮL）と符号系列Ｌkの相関関数φkは、
【数７】
0005419024.tif 000008
であり、ｋ1／２＝ｋ’1とすると、
【数８】
0005419024.tif 000009
【００４５】
よって、φkが最小となるｋがその画素の変調に利用された符号系列を示す。１〜ＮL時相までφkが最小となる画素を辿れば、初期時相で、（ｘ，ｙ）に位置していた画素の移動を追跡することが可能となる。
【００４６】
ただし、実際にＭＲＩ装置１で撮像を行うと、撮像の条件によって、符号系列の各要素に対応する画像ごとにｋ1，ｋ2は異なる。さらに、画素の位置によってもｋ1，ｋ2は変化する。そのため実際には、数６の式（４）は、
【数９】
0005419024.tif 000010
となる。
【００４７】
この場合、式５（数８）の式が成り立たず、割り当てられた符号系列Ｌiの特定が困難となる。そこで、画素位置や撮像条件による影響を除去するため、画素値ｍj(i)の除算を行うと、
【数１０】
0005419024.tif 000011
となる。
【００４８】
ｋ1（ｘ，ｙ，ｊ）／ｋ2（ｘ，ｙ，ｊ）＝ｋ3（ｘ，ｙ，ｊ）とすると、ｋ3は、タグのコントラストに相当する。ｋ3（ｘ，ｙ，ｊ）がｊに関わらずほぼ一定であると仮定すると、相関関数φkは、
【数１１】
0005419024.tif 000012
となり、やはり、φｋが最小となるｋを求めることにより、各画素に対応する符号系列を決定することができる。
【００４９】
［５．タグパターンの例］
ここでは、符号系列としてＭ系列を利用する。Ｍ系列は、ＮL＝２n−１の符号長、系列数を持ち、τ＝０における自己相関値と相互相関値の比が−ＮL：１となる。本実施形態において、画素値系列から符号系列を特定する際には、自己相関値がピークを持っていることが望ましいため、ＮLは大きいことが望まれる。一方、ＮLが大きいと、ＮL個の時系列タギングＭＲ画像を撮像する必要があり、撮像に時間を要する。そこで、ここでは、ＮL＝７とした。
ＮL＝７のＭ系列として、２次の原始多項式ｘ2＋ｘ＋１で生成される符号系列を用いた。前述のＤＡＮＴＥタギングシーケンス（又はＳＰＡＭＭタギングシーケンス）では、一定長さのタグ領域と一定長さの非タグ領域が交互に繰り返す周期的なタグパターンが実現できるため、生成した７つのＭ系列から符号系列を選択し、適切に空間配置することで、全ての撮影において周期的なタグパターンを成すことができる。
【００５０】
本実施形態では、７符号系列中６符号系列を選択して並べた符号配置を利用した。つまり、Ｎｓ＝６となる。下記表１は、６つのＭ系列Ｌ1〜Ｌ6を示している。
【表１】
0005419024.tif 000013
【００５１】
表１において、横方向の各列が、一つのタグパターンに対応する。表１の「−１」を高輝度領域、「１」をタグに対応する低輝度領域となるようにＤＡＮＴＥタギングシーケンスを設計した場合の７つのタグパターン（画像中の一部（１８×１８画素分）の拡大図）を図９に示す。
図９の各タグパターンは、対応する画像領域（１８×１８画素分）内の各画素に、図１０に示すように、符号系列Ｌ1〜Ｌ6を割り当てた場合のものである。各画素に割り当てられた符号系列Ｌ1〜Ｌ6は、画素の位置情報を示していることになる。
【００５２】
上記表１のような符号系列を用いると、例えば、撮像対象物の縦緩和時間Ｔ１の約２倍の時間まで位置情報を誤り無く符号することができる。
【００５３】
なお、図１０では、縦方向にはすべて同一の符号系列が割り当てられている。これは、単純化のため、複数のタグパターンすべてを縦縞状としたためである。しかし、図１０において縦方向も異なる符号系列が割り当てられるようにしてもよいのは当然である。横方向及び縦方向において隣接する画素の符号系列が異なるようにするには、複数のタグパターンの中に、横縞状のものや格子状のものを混在させることで実現できる。そのようなタグパターンもタギングシーケンスを適宜設定することで生成可能である。
【００５４】
［６．実験結果］
図１１（ａ）〜（ｃ）は、図１２に示すファントム（撮像対象）が、左右方向に正弦運動する様子をタギングＭＲＩ装置１にて撮影した結果を示している。
図１１（ａ）は、図９に示す７つのタグパターンにて７個のｃｉｎｅ画像をＭＲＩ装置１にて撮像し（ＮL＝７）、先に説明したＮL＝７のＭ系列にて、ファントムの位置の変位を判定した結果を示している。
【００５５】
図１１（ｂ）は、左右方向に正弦運動するファントムの時系列画像（第１時相から第３３時相まで）において、各時相の画像から、同一位置の領域（垂直方向に１画素×水平方向に約５０画素の領域）を切り出し、この短冊状の領域を、時刻（時相）に沿って、図１１の縦軸方向に並べたものである。つまり、図１１（ａ）の縦軸は、時間（第１時相から第３３時相）を示しており、横軸は、画素位置を示している。
図１１（ａ）においては、判定した符号系列Ｌ1〜Ｌ6を、画素の明暗で区別して示している。つまり、図１１（ａ）においては、明暗度が同じ画素は同一の符号系列であると判定された画素となっている。
【００５６】
図１１（ｂ）は、上記Ｍ系列を構成する７個のｃｉｎｅ画像のうち、タグ領域（低輝度領域）及び非タグ領域（高輝度領域）それぞれの左右幅が１ｍｍ（１画素分）のタグパターン（１画素の幅を１ｍｍとすると図９の第１タグパターンに相当）でファントムを撮影した画像である。
また、図１１（ｃ）は、上記Ｍ系列を構成する７個のｃｉｎｅ画像のうち、タグ領域（低輝度領域）及び非タグ領域（高輝度領域）それぞれの左右幅が３ｍｍ（３画素分）のタグパターン（１画素の幅を１ｍｍとすると図９の第３〜第５タグパターンのいずれかに相当）でファントムを撮影した画像である。
【００５７】
図１１（ｂ）及び図１１（ｃ）においても、各時相の画像から、同一位置の領域（垂直方向に１画素×水平方向に約５０画素の領域）が切り出され、この短冊状の領域が、時刻（時相）順に、図１１の縦軸方向に並べられている。
【００５８】
図１１（ａ）に示すように、Ｍ系列を用いたタギングＭＲ画像の判定結果では、第１３〜１４時相（上から１３〜１４行目）付近まで、１ｍｍ（１画素）間隔で、ファントムの移動（左右方向の正弦運動）をほぼ正しく判定できている。
これに対し、図１１（ｂ）に示す１ｍｍ幅のタギングＭＲＩ画像では、タグ領域及び非タグ領域が細いため、両領域のコントラストが低くなり、第１０時相までにタグパターンは消失し、それ以降の追跡ができない。
また、図１１（ｃ）に示す３ｍｍ幅のタギングＭＲＩ画像では、Ｍ系列による判定結果と同様に、第１３〜１４時相まで淡いコントラストでタグパターンを追跡できるが、太いタグパターンのため、位置判定の分解能は３ｍｍにすぎない。
【００５９】
このように、図１１（ｂ）及び（ｃ）のように元画像のタグパターンそのものを追跡する場合には、位置の検出可能な時間が短くなったり、あるいは検出精度が低下したりするのに対し、図１１（ａ）のＮL個の画像を用いた符号判定では、比較的長い時間、高い精度で撮像対象の位置を検出することができる。
【００６０】
［７．サブピクセル解析（Sub-Pixel Analysis）］
上記のＭ系列による判定では、１画素単位の解析（Pixel analysis）となり、位置検出の精度は１画素単位である。したがって、１画素未満の移動量を検出することはできない。
【００６１】
図１３は、タグパターンが付された撮像対象（ファントム）を撮像した時系列タギングＭＲ画像の初期時相（第１時相）における拡大図を示している。図１３では、横３画素（Ｘ１，Ｘ２，Ｘ３）×縦２画素（Ｙ１，Ｙ２）の計６画素を示している。
図１３に示すように、初期時相では、撮像対象に標識されたタグパターンにおける低輝度領域（黒領域）と高輝度領域（白領域）との境界は、画素の境界と一致している。図１３において、黒領域の横幅は、ちょうど２画素分である。
【００６２】
ところが、時間の経過に伴って、撮像対象が移動すると、各画素においては、撮像対象の低輝度領域（黒領域）と高輝度領域（白領域）とが混ざった状態で撮影されることになる。図１４は、撮像対象が図１３の位置から左に０．２５画素分移動したときの実際の位置と撮影した画像における画素値とを示している。また、図１５は、撮像対象が図１３の位置から左に０．５画素分移動したときの実際の位置と撮像した画層における画素値とを示している。
【００６３】
図１４においては、図１４（ａ）に示すように撮像対象は、左に０．２５画素分移動している。ただし、画素値（輝度値）は、画素単位でしか得られない。したがって、撮像画像における画素値は、図１４（ｂ）に示すように、Ｘ２の画素は黒領域となるものの、Ｘ１の画素は、白にわずかに黒みがかかった薄い灰色の領域となり、Ｘ３の画素は、黒にわずかに白みがかかった濃い灰色の領域となる。なお、図１３〜及び図１５の画素内に描かれたメッシュは、画素値（輝度）の大きさに対応しており、メッシュが荒いほど白（高輝度）に近く、メッシュが細かいほど黒（低輝度）に近いことを示している。
【００６４】
さらに、図１５においては、図１５（ａ）に示すように撮像対象は、左に０．５画素分移動している。この場合の画素値は、図１５（ｂ）に示すように、Ｘ２の画素は黒領域となるものの、Ｘ１及びＸ３の画素は、白と黒の中間の灰色の領域となる。
【００６５】
ここでの黒領域の横幅は、図１３に示すように、本来、２画素分であるが、撮像対象が１画素未満の移動を生じると、図１４（ｂ）及び図１５（ｂ）に示すように黒領域と白領域との境界がぼやけることになる。
【００６６】
このように、両領域の境界がはっきりしない場合においても、先に説明した符号系列の検出（符号判定）は、画素単位で行われる。したがって、図１３に示す初期時相において、画素Ｘ１が第１符号系列Ｌ1に、画素Ｘ２が第２符号系列Ｌ2に、画素Ｘ３が第３符号系列Ｌ3に対応しているものとすると、図１４の状態では、図１３と同様に、画素Ｘ１の画素値系列が第１符号系列Ｌ1、画素Ｘ２の画素値系列が第２符号系列Ｌ2、画素Ｘ３の画素値系列が第３符号系列Ｌ3であると検出される可能性が高い。この場合、撮像対象が左に０．２５画素分移動しているにもかかわらず、全く移動していないと判定されることになる。
【００６７】
また、図１５の状態では、図１４と同様に検出されるか、あるいは、画素Ｘ１の画素値系列が第２符号系列Ｌ2、画素Ｘ２の画素値系列が第３符号系列Ｌ3、画素Ｘ３の画素値系列が第４符号系列Ｌ4であると検出される可能性が高い。前者の場合、撮像対象は左に０．５画素分移動しているにもかかわらず、全く移動していないと判定されることになる。後者の場合、撮像対象は左に０．５画素分しか移動しているにもかかわらず、１画素分移動したと判定されることになる。
【００６８】
１画素未満の移動量を求めるサブピクセル解析（Sub-Pixel Analysis）を行う場合には、図６に示すステップＳ２２において、各画素の画素値系列から、唯一の符号系列を決定するのではなく、各画素において複数の符号系列（Ｌ1〜Ｌ6）が混在している可能性も含めて検出する。
【００６９】
各画素において複数の符号系列（Ｌ1〜Ｌ6）が混在している可能性も含めて検出するには、各画素の画素値系列と複数の符号系列（Ｌ1〜Ｌ6）それぞれとの相関値を用いる。
ここでは、画素値系列と特定の符号系列との相関値を０〜１で表し、「１」の場合、相関が最も高く、「０」の場合、相関が最も低いものとする。
【００７０】
図１３に示す初期時相のように、画素の境界と撮像対象の白黒領域の境界とが一致している場合、各画素には、複数の符号系列は混在せず、画素値系列は、唯一の符号系列に対応することになる。したがって、各画素の画素値系列と複数の符号系列（Ｌ1〜Ｌ6）それぞれとの相関値は、理論的には、図１３に示すようになる。具体的には、画素Ｘ１については、第１符号系列Ｌ1との相関が最大となり、他の符号系列Ｌ2〜Ｌ6との相関は最小となる。画素Ｘ２については、第２符号系列Ｌ2との相関が最大となり、他の符号系列Ｌ1，Ｌ3〜Ｌ6との相関は最小となる。画素Ｘ３については、第３符号系列Ｌ3との相関が最大となり、他の符号系列Ｌ1，Ｌ2，Ｌ4〜Ｌ6との相関は最小となる。
【００７１】
図１４（ａ）に示すように撮像対象が左に０．２５画素分移動すると、各画素には、複数の符号系列が混在することになる。したがって、各画素の画素値系列と複数の符号系列（Ｌ1〜Ｌ6）それぞれとの相関値は、理論的には、各画素において各符号系列が占める領域の面積の大小に応じた値となる。
図１４（ａ）に示すように、画素Ｘ１については、第１符号系列Ｌ1及び第２符号系列Ｌ2が混在している。このため、他の符号系列Ｌ3〜Ｌ6との相関は最小となるのに対し、第１符号系列Ｌ1との相関は０．７５、第２符号系列Ｌ2との相関は０．２５というように、各画素において符号系列Ｌ1，Ｌ2が占める領域の面積の大小に応じた値をとる。画素Ｘ２，Ｘ３についても同様に、混在する符号系列についてそれらの符号系列が占める領域の面積の大小に応じた値をとる。
【００７２】
さらに、図１５（ａ）に示すように撮像対象が左に０．５画素分移動すると、画素Ｘ１については、他の符号系列Ｌ3〜Ｌ6との相関は最小となるのに対し、第１符号系列Ｌ1との相関は０．５、第２符号系列Ｌ2との相関は０．５というように、各画素において符号系列Ｌ1，Ｌ2が占める領域の面積の大小に応じた値をとる。画素Ｘ２，Ｘ３についても同様に、混在する符号系列についてそれらの符号系列が占める領域の面積の大小に応じた値をとる。
【００７３】
処理プロセッサ９は、図１４及び図１５に記載のような各画素についての複数の符号系列（Ｌ1〜Ｌ6）それぞれとの相関値に基づいて、それぞれの画素に含まれる１又は複数の符号系列（Ｌ1〜Ｌ6）を特定するとともに、各画素において各符号系列が占める割合を求める。
【００７４】
具体的には、図１６〜図１８に示すように、複数の符号系列（Ｌ1〜Ｌ6）それぞれとの相関値と、移動量との関係とを規定した情報（移動量特定情報）が、記憶部１０に予め設定されており、処理プロセッサ９は、図１６〜図１８に示す情報を参照することで、符号系列の特定等を行う。
【００７５】
図１６は、ある画素に、符号系列Ｌ1だけが存在している状態（図１６の横軸の「０」の位置）から、符号系列Ｌ2が進入してきて、その画素が符号系列Ｌ2だけで占められた状態（図１６の横軸の「１画素」の位置）となるまでの間における、その画素の画素値系列と符号系列（Ｌ1〜Ｌ6）それぞれとの相関値の変化を示している。
【００７６】
図１７は、ある画素に、符号系列Ｌ2だけが存在している状態（図１７の横軸の「０」の位置）から、符号系列Ｌ3が進入してきて、その画素が符号系列Ｌ3だけで占められた状態（図１６の横軸の「１画素」の位置）となるまでの間における、その画素の画素値系列と符号系列（Ｌ1〜Ｌ6）それぞれとの相関値の変化を示している。
【００７７】
図１８は、ある画素に、符号系列Ｌ3だけが存在している状態（図１８の横軸の「０」の位置）から、符号系列Ｌ4が進入してきて、その画素が符号系列Ｌ4だけで占められた状態（図１８の横軸の「１画素」の位置）となるまでの間における、その画素の画素値系列と符号系列（Ｌ1〜Ｌ6）それぞれとの相関値の変化を示している。
【００７８】
なお、図１６〜図１８には、図１３〜図１５における撮像対象の動きに関連する３つの情報（移動量特定情報）だけを示したが、実際には、１つの画素に同時に存在する可能性がある符号系列の組み合わせの数だけの情報が記憶部１０に設定されている。
【００７９】
処理プロセッサ９は、例えば、画素Ｘ１について図１４（ａ）示すような各符号系列（Ｌ1〜Ｌ6）との相関値が得られた場合、記憶部１０に記憶されている情報（移動量特定情報）をサーチして、得られた相関値とフィットする情報を選択する。ここでは、１つ前の時相（初期時相；図１３）において、画素Ｘ１の符号系列は符号系列Ｌ1であるから、記憶部１０に記憶されている情報のうち、横軸＝０の位置において符号系列Ｌ1との相関値が最大となっている図１６の情報が選択される。
そして、符号系列Ｌ1との相関値＝０．７５及び符号系列Ｌ2との相関値＝０．２５という混在状況に基づいて、図１６に示す情報を参照して、符号系列Ｌ2が進入してきた移動量は「０．２５」であると判定する。
【００８０】
同様にして、図１４（ａ）に示す他の画素及び図１５（ａ）に示す画素についても、混在する符号系列及び進入してきた符号系列の移動量（その画素に占める割合）を求めることができる。
【００８１】
なお、図１３〜図１５のように、単純な左右方向の往復運動（正弦運動）の場合、各符号系列の占める割合だけからでも、撮像対象の移動方向を特定することが可能であるが、より複雑な運動となった場合には、各画素における符号系列の割合の変化の仕方を、連続する複数の時相間で検証することで、移動方向を推定することが可能である。
また、図１６〜図１８に示す情報（移動量特定情報）では、相関値と移動量との関係が線形となっているが、実際には、非線形となることが多い。相関値と移動量との関係が非線形性となる場合は、各符号系列の組み合わせ毎に、予め符号系列から計算により求めて置いた関係、あるいは予め実測によりもとめておいた関係を、移動量特定情報として用いるのが好ましい。
【００８２】
以上のようなサブピクセル解析を行うことで、１画素の大きさ未満での撮像対象の運動量を算出することができる。
図１９は、左右方向に正弦運動する撮像対象（ファントム）に対してサブピクセル解析とピクセル解析とを行った場合の解析結果のイメージを示している。撮像対象のある点に着目すると、図１９における実線（理論位置）で示すように、運動は正弦運動となる。この撮像対象に対して、ピクセル解析を行った場合、位置は、１ｍｍ単位でしか検出できないことになる。これに対し、サブピクセル解析では、１ｍｍ未満の移動も検出できるため、理論位置に近い位置を検出することが可能となる。
【００８３】
［８．変形例］
上記実施形態では、一つの画素それぞれに符号系列（位置情報）を割り当てたが、符号系列（位置情報）を割り当てる対象は、一つの画素に限らず、複数の画素からなる画像領域であってもよい。
【００８４】
しかも、符号系列（位置情報）を割り当てる対象を、１個の時系列タギングＭＲ画像におけるＮL個の画素からなる画像領域とした場合、符号系列は、ＮL方向ではなく、画像方向（２次元画像のＸ方向若しくはＹ方向、又は３次元画像のＸ方向，Ｙ方向若しくはＺ方向）に形成することができる。この場合、複数の時系列タギングＭＲ画像を撮像する必要がなくなり、撮像時間を短縮することができる。
【００８５】
図２０は、画像方向（Ｘ方向）に符号系列「１，−１，−１，１，−１，１，１」を割り当てた例を示している。この場合、時相（時刻）ｔ＝０において、画像の左端に位置していた前記符号系列を有する領域は、後の時相（時刻）ｔ＝ｔ１においては、撮像対象（人体組織）の運動により、画像の右方へ移動している。時相ｔ＝ｔ１において、前記符号系列が割り当てられている領域を検出するには、時相ｔ＝ｔ１の画像において、前記符号系列「１，−１，−１，１，−１，１，１」との相関をとり、相関値が大きくなる位置を検出することで行える。
例えば、図示のように、時相ｔ＝ｔ１の符号系列「１，−１，−１，１，−１，１，１」の位置と同じ位置で相関をとった場合、その位置から１画素分右へ移動した位置で相関をとった場合、２画素分右へ移動した位置で相関をとった場合、及び３画素分右へ移動した位置で相関をとった場合では、３画素分右へ移動した位置で相関をとった場合に相関が最も高くなり、その他の場合では、相関が−１（最小値）である。したがって、時相ｔ＝ｔ１においては、符号系列「１，−１，−１，１，−１，１，１」が割り当てられた領域は、３画素分右へ移動した位置にあることを検出することができる。
【００８６】
［９．変形例におけるサブピクセル解析］
この変形例においても、１画素単位の解析（ピクセル解析）だけでなく、１画素未満の解析（サブピクセル解析）を行うことができる。この変形例において、サブピクセル解析を行う場合には、各画素の画素値（輝度値）から移動量（白領域又は黒領域の占める割合）を求める。
つまり、タグパターンにおいて、白領域と黒領域の境界が１つの画素中に位置した場合、図２２（ｂ）及び図２３（ｂ）のように、その画素は、白領域と黒領域との間の中間色（中間輝度）となる。例えば、白領域部分だけが撮影された画素（例えば図２１の画素Ｘ１）に、隣の画素（例えば、図２１の画素Ｘ２）から黒領域部分が移動してきた場合、その画素Ｘ１の画素値は、その画素Ｘ１に占める黒領域の大きさに応じた強さの黒みがかかった値となる。
また、黒領域部分だけが撮影された画素（例えば図２１の画素Ｘ３）に、隣の画素から黒領域部分が移動してきた場合、その画素Ｘ３の画素値は、その画素Ｘ３に占める白領域の大きさに応じた強さの白みがかかった値となる。
【００８７】
タグパターンにおける黒領域の輝度を「０」、白領域の輝度の輝度を「１」とした場合、白領域画素に黒領域部分が移動してきた場合の画素値（輝度）と移動量とは、図２４（ａ）のように示される。また、黒領域画素に白領域部分が移動してきた場合の画素値（輝度）と移動量とは、図２４（ｂ）のように示される。
【００８８】
なお、図２４の画素値は、同じ撮像対象の同じ運動についてタグパターンが付されていない時系列ＭＲ画像における画素値で正規化したものである。撮像対象の位置によっては、同じ白領域又は黒領域であっても画素値が異なることがあり、また、時間が経つにつれてタグパターンが消失（例えば、黒領域も白領域となる）するが、上記の正規化を行っておくことで、条件により異なる画素値（輝度）が得られても、図２４に示す関係を適用することができる。また、正規化は、タグパターンが付されていない時系列ＭＲ画像における画素値ではなく、同一時相における白領域若しくは黒領域の画素値又はそれらのコントラストを用いて行っても良い。
【００８９】
図２４に示す関係（移動量特定情報）を用いて、移動量を求めるには、まず、画素値から、画素内において白領域が占める割合ｆwと、黒領域が占める割合ｆBを求める。これらの割合ｆw，ｆBは、ある画素の画素値を、正規化することによって得られる。
例えば、時間が経つにつれて、タグパターンの黒領域が消失する場合、時相ごとに黒領域がとる値は、白領域の値に近づいていくため変化する。これに対処するには、同じ撮像対象の同じ運動についてタグパターンが付されていない時系列ＭＲ画像において、割合ｆw，ｆBの算出対象となっている画素と、同一時相かつ同一位置の画素（対応画素）に着目する。正規化では、割合ｆw，ｆBの算出対象となっている画素の画素値を、対応画素の画素値で除算する。これにより、算出対象の画素値の画素値は、「１」（白領域のみ）から「０」（黒領域のみ）の値で正規化され、その値がｆWとなる。ｆBは、［１−ｆW］の演算で求めることができる。
【００９０】
割合ｆw，ｆBが求まると、図２４に示す関係を用いて、それらの  割合ｆw，ｆBから算出対象画素への白領域又は黒領域の進入量（移動量）を１画素未満の大きさで求めることができる。なお、移動方向は、過去の時相における白領域又は黒領域の配置や算出対象画素の周囲の画素の画素値に基づいて推定可能であるが、規則的な左右方向の正弦運動を前提とする場合、移動量がわかるだけで運動解析が可能である。
【００９１】
［１０．変形例の２次元（多次元）への拡張］
なお、図２０において、符号系列（位置情報）が割り当てられる画素の領域は、１次元（Ｘ方向のみ）の領域であったが、２次元又は３次元の領域であってもよい。また、図７に示す複数回（例えば、ＮL回）の撮像方向も加えて考えた場合、２次元画像であれば画像Ｘ方向、画像Ｙ方向、及び撮像方向の３次元における領域であってもよいし、３次元画像であれば、画像Ｘ方向、画像Ｙ方向、画像Ｚ方向、及び撮像方向の４次元における領域であってもよい。
【００９２】
図２５及び図２６は、図２０に示す変形例を２次元に拡張した例を示している。図２５及び図２６に示す例では、１つの符号系列（位置情報）が割り当てられる画素の領域は、２次元の領域（３画素×３画素）となっている。
図２５（ａ）に示すタグパターンは、変形例を２次元に拡張するためのパターンである。このパターンは、図２５（ｂ）にも示すように、横縞（水平の縞模様）と縦縞（垂直の縞模様）とを有するパターンである。図２５（ｂ）のパターンは、画像の縦方向にみたときに、横縞の白領域（白帯）幅：黒領域（黒帯）幅＝３：１であり、画像の横方向にみたときに、縦縞の白領域（白帯）幅：黒領域（黒帯）幅＝３：２である。なお、図２５（ｂ）の斜線部分が黒領域であり、白抜き領域が白領域である。
【００９３】
このようなタグパターンを撮像対象に印加すると、図２５（ａ）に示すように、２次元空間で、２０個（縦４画素×横５画素）の画素（図２５（ｂ）のＡ−Ｔ）を１単位として左右上下に繰り返す白黒パターンが形成される。
Ａ−Ｔの２０個の画素について、それぞれを中心とする例えば９個の画素（自身の画素及びその近傍の８個の画素；３×３画素の領域）の白黒パターンに注目すると、図２６に示すように、２０個の３×３画素の領域は、すべて異なる白黒パターンとなっている。これは、２０個の画素に特有の画像パターン（長さＮL＝９の符号系列；位置情報）とみなすことができ、タグパターンによって画像の各画素をコーディングできたことと等価である。
【００９４】
図２７〜図２９は、変形例を２次元に拡張した他の例を示している。
図２７（ｂ）に示すパターンも、横縞（水平の縞模様）と縦縞（垂直の縞模様）とを有するパターンである。図２７（ｂ）のパターンは、画像の縦方向にみたときに、横縞の白領域（白帯）幅：黒領域（黒帯）幅＝３：２であり、画像の横方向にみたときに、縦縞の第１白領域幅：第１黒領域幅：第２白領域幅：第２黒領域幅＝３：１：１：１である。なお、図２７（ｂ）の斜線部分が黒領域であり、白抜き領域が白領域である。
【００９５】
このようなタグパターンを撮像対象に印加すると、図２７（ａ）に示すように、２次元空間で、３０個（縦５画素×横６画素）の画素（図２７（ｂ）のＡ−Ｙ（Ａ’，Ｆ’，Ｋ’，Ｐ’，Ｕ’を含む））を１単位として左右上下に繰り返す白黒パターンが形成される。
Ａ−Ｙの３０個の画素について、それぞれを中心とする例えば９個の画素（自身の画素及びその近傍の８個の画素；３×３画素の領域）の白黒パターンに注目すると、Ａ’，Ｆ’，Ｋ’，Ｐ’，Ｕ’は、Ａ，Ｆ，Ｋ，Ｐ，Ｕと同じパターンになるため、図２８に示すように、２５個の白黒パターンが得られ、これら２５個の３×３画素の領域は、すべて異なる白黒パターンとなっている。これは、２５個の画素に特有の画像パターン（長さＮＬ＝９の符号系列；位置情報）とみなすことができ、タグパターンによって２５個の画素をコーディングできたことと等価である。
【００９６】
Ａ’，Ｆ’，Ｋ’，Ｐ’，Ｕ’の画素と、Ａ，Ｆ，Ｋ，Ｐ，Ｕの画素との区別については、図２９に示すように、９個の画素（自身の画素及びその近傍の８個の画素；３×３画素の領域）よりもさらに周辺の画素（例えば、自身の画素及びその近傍の２４個の画素；５×５画素の領域）に注目することにより判別することが可能である。なお、３０個の画素全部について、５×５画素の領域で符号化をしてもよい。
【００９７】
なお、以上説明した変形例において特に説明しない点については、図１〜図１０に関して説明した内容に即して実施することができる。
【００９８】
なお、上記において開示した事項は、例示であって、本発明を限定するものではなく、様々な変形が可能である。
例えば、画素信号が、輝度値に代えて、又は加えて、位相値であってもよい。
また、運動解析としては、単に、移動などを検出するのではなく、検出された移動量から撮像対象の弾性率を求めるものであってもよい。
さらに、本発明の撮像対象は、特に限定されない。ただし、解析精度の観点からは、血液などの流体のように基本的形態を持たない撮像対象よりも、心筋組織などのように移動・収縮・拡張などが生じるとしても基本的形態が保持される撮像対象の方が、良好な解析精度を得やすく好ましい。つまり、例えば血流を想定すると、血流方向の上流側に位置する部分が、下流側に位置する部分を追い越してしまい、ある時相における撮像対象（血流）の連続的な配置関係が、後の時相では崩れてしまう可能性があり、良好な解析精度を得るのを困難にする。一方、心筋組織は、運動（移動・収縮・拡張）しても、ある時相における心筋組織の連続的な配置関係が維持される連続体であるため、解析が容易となり、良好な解析精度を得るのが容易となる。
【符号の説明】
【００９９】
１  ＭＲＩ装置
２  磁場発生部
３  ＲＦコイル
４  ＲＦ送信部
５  ＲＦ受信部
６  制御部
７  Ａ／Ｄ変換部
８  コンピュータ
９  処理プロセッサ
１０  記憶部
１１  表示部
【図１】
0005419024.tif 000014
【図２】
0005419024.tif 000015
【図３】
0005419024.tif 000016
【図４】
0005419024.tif 000017
【図５】
0005419024.tif 000018
【図６】
0005419024.tif 000019
【図７】
0005419024.tif 000020
【図８】
0005419024.tif 000021
【図９】
0005419024.tif 000022
【図１０】
0005419024.tif 000023
【図１２】
0005419024.tif 000024
【図１４】
0005419024.tif 000025
【図１５】
0005419024.tif 000026
【図１６】
0005419024.tif 000027
【図１７】
0005419024.tif 000028
【図１８】
0005419024.tif 000029
【図２０】
0005419024.tif 000030
【図２２】
0005419024.tif 000031
【図２３】
0005419024.tif 000032
【図２４】
0005419024.tif 000033
【図２６】
0005419024.tif 000034
【図２８】
0005419024.tif 000035
【図２９】
0005419024.tif 000036
【図１１】
0005419024.tif 000037
【図１３】
0005419024.tif 000038
【図１９】
0005419024.tif 000039
【図２１】
0005419024.tif 000040
【図２５】
0005419024.tif 000041
【図２７】
0005419024.tif 000042