"放射線を射出する放射線源と、該放射線源から射出された放射線を透過する部分と遮蔽する部分とからなる格子構造が周期的に配置され、前記放射線源から射出された放射線を通過させて周期パターン像を形成する第1の格子と、該第1の格子により形成された周期パターン像を透過する部分と遮蔽する部分とからなる格子構造が周期的に配置され、位相情報パターンを形成することを可能とする第2の格子と、該第2の格子により形成された位相情報パターンを検出する放射線画像検出器とを備えた放射線位相画像撮影装置であって、前記放射線画像検出器が、前記第1および第2の格子の少なくとも一方の予め設定された配置からの配置ずれによって形成されるモアレ像を検出するものであり、前記放射線画像検出器によって検出されたモアレ像を表すモアレ画像信号を取得し、該取得したモアレ画像信号に対してフーリエ変換処理を施して前記モアレ像の周波数成分を算出するモアレ周波数算出部を備えたことを特徴とする放射線位相画像撮影装置。前記モアレ周波数算出部によって算出された周波数成分が、前記第1の格子と前記第2の格子との相対的な位置関係に基づいて予め設定された周波数成分に近づくように前記少なくとも一方の格子の配置を調整する調整部を備えたことを特徴とする請求項1記載の放射線位相画像撮影装置。前記第1および第2の格子が、前記第1の格子の延伸方向と前記第2の格子の延伸方向とが平行になるように設置されたものであり、前記調整部が、前記モアレ周波数算出部によって算出された周波数成分が低周波成分になるように前記少なくとも一方の格子の配置を調整するものであることを特徴とする請求項2記載の放射線位相画像撮影装置。前記調整部が、前記モアレ周波数算出部によって算出された周波数成分に基づいて、前記格子のピッチに基づいて予め設定された配置情報との配置ずれ量を取得し、該配置ずれ量に基づいて前記少なくとも一方の格子の配置を調整するものであることを特徴とする請求項2または3記載の放射線位相画像撮影装置。前記調整部が、前記モアレ周波数算出部によって算出された周波数成分と前記少なくとも一方の格子の配置情報との相関関係を取得し、該相関関係に基づいて、前記モアレ周波数算出部によって算出された周波数成分が前記予め設定された周波数成分に近づくような配置情報を取得し、該配置情報に基づいて前記少なくとも一方の格子の配置を調整するものであることを特徴とする請求項2または3記載の放射線位相画像撮影装置。前記調整部が、前記モアレ周波数算出部によって算出された周波数成分に基づいて、前記少なくとも一方の格子の配置ずれ量として回転ずれ量を取得するものであり、該回転ずれ量に基づいて前記少なくとも一方の格子の配置を調整するものであることを特徴とする請求項1から5いずれか1項記載の放射線位相画像撮影装置。前記格子の中心を原点として、前記格子の面内において直交する軸をX軸およびY軸とするとともに前記格子の面に対して垂直方向に延びる軸をZ軸とした場合において、前記調整部が、前記回転ずれ量として、前記X軸を回転軸とした回転ずれ量、前記Y軸を回転軸として回転ずれ量および前記Z軸を回転軸として回転ずれ量の少なくとも1つを取得するものであることを特徴とする請求項6記載の放射線位相画像撮影装置。前記第1の格子によって形成される周期パターン像の周波数が前記第2の格子の周波数と一致するように前記第1の格子と前記第2の格子のピッチが設定されているとともに、前記予め設定された配置が決定されていることを特徴とする請求項1から7いずれか1項記載の放射線位相画像撮影装置。前記モアレ周波数成分算出部が、前記第1の格子によって形成される周期パターン像の周波数と前記第2の格子の周波数との差からなるモアレ像の周波数成分を算出するものであることを特徴とする請求項3項記載の放射線位相画像撮影装置。前記調整部が、前記モアレ周波数算出部によって算出された周波数成分がゼロに収束するような前記少なくとも一方の格子の配置ずれ量を取得するものであることを特徴とする請求項3項記載の放射線位相画像撮影装置。前記第1の格子と前記第2の格子とが、タルボ干渉計を構成するように配置されていることを特徴とする請求項1から10いずれか1項記載の放射線位相画像撮影装置。そこで、第1の格子や第2の格子などの位置決めを実現する手順の確立が重要となる。研究レベルの装置では高精度の測定器を使用すれば当然これらの位置決めを行うことが可能であるが、製品化して実際に設置などを行う場合を想定すると、高度な測定器の使用を前提にするのは現実的でないという問題がある。一方、特許文献2においては、2つの格子の相互の位置ずれによって発生するモアレ縞が消えるように位置調整することによって2つの格子の位置決めをすることが提案されている。しかしながら、たとえばモアレ縞を目視によって観察しながら位置調整したのでは、低周波側の視認限界により位置調整の精度に限界があり、高精度な位置調整を行うことは困難である。また、モアレ縞の周期が放射線画像検出器の画素サイズ程度の高周波である場合には、放射線位相コントラスト画像に高周波ノイズが表れるような致命的な影響を与えてしまうが、高周波のモアレ縞の目視についても視認限界があり、上述したような画素サイズ程度の高周波のモアレ縞は目視観測することができない。したがって、モアレ縞の目視による位置調整では、目視観測することができない高周波のモアレ縞が存在しているにも関わらず、視認することができないためにモアレ縞が見えなくなったと誤認してしまう可能性があり、適切な位置調整をすることができない。本発明は、上記の事情に鑑み、第1の格子および第2の格子などの位置調整を高精度に行うことができる放射線位相画像撮影装置を提供することを目的とする。本発明の放射線位相画像撮影装置は、放射線を射出する放射線源と、放射線源から射出された放射線を透過する部分と遮蔽する部分とからなる格子構造が周期的に配置され、放射線源から射出された放射線を通過させて周期パターン像を形成する第1の格子と、第1の格子により形成された周期パターン像を透過する部分と遮蔽する部分とからなる格子構造が周期的に配置され、位相情報パターンを形成することを可能とする第2の格子と、第2の格子により形成された位相情報パターンを検出する放射線画像検出器とを備えた放射線位相画像撮影装置であって、放射線画像検出器が、第1および第2の格子の少なくとも一方の予め設定された配置からの配置ずれによって形成されるモアレ像を検出するものであり、放射線画像検出器によって検出されたモアレ像を表すモアレ画像信号を取得し、その取得したモアレ画像信号に対してフーリエ変換処理を施してモアレ像の周波数成分を算出するモアレ周波数算出部を備えたことを特徴とする。また、上記本発明の放射線位相画像撮影装置においては、モアレ周波数算出部によって算出された周波数成分が、第1の格子と第2の格子との相対的な位置関係に基づいて予め設定された周波数成分に近づくように少なくとも一方の格子の配置を調整する調整部を設けることができる。また、第1および第2の格子を、第1の格子の延伸方向と第2の格子の延伸方向とが平行になるように設置するとともに、調整部を、モアレ周波数算出部によって算出された周波数成分周波数成分が低周波成分になるように少なくとも一方の格子の配置を調整するものとできる。また、調整部を、モアレ周波数算出部によって算出された周波数成分に基づいて、格子のピッチに基づいて予め設定された配置情報との配置ずれ量を取得し、その配置ずれ量に基づいて少なくとも一方の格子の配置を調整するものとできる。また、調整部を、モアレ周波数算出部によって算出された周波数成分と少なくとも一方の格子の配置情報との相関関係を取得し、その相関関係に基づいて、モアレ周波数算出部によって算出された周波数成分が上記予め設定された周波数成分に近づくような配置情報を取得し、その配置情報に基づいて少なくとも一方の格子の配置を調整するものとできる。また、調整部を、モアレ周波数算出部によって算出された周波数成分に基づいて、少なくとも一方の格子の配置ずれ量として回転ずれ量を取得するものであり、その回転ずれ量に基づいて少なくとも一方の格子の配置を調整するものとすることができる。また、格子の中心を原点として、格子の面内において直交する軸をX軸およびY軸とするとともに格子の面に対して垂直方向に延びる軸をZ軸とした場合において、調整部を、回転ずれ量として、X軸を回転軸とした回転ずれ量、Y軸を回転軸として回転ずれ量およびZ軸を回転軸として回転ずれ量の少なくとも1つを取得するものとできる。また、第1の格子によって形成される周期パターン像の周波数が第2の格子の周波数と一致するように第1の格子と第2の格子のピッチを設定するとともに、予め設定された配置を決定することができる。また、モアレ周波数成分算出部を、第1の格子によって形成される周期パターン像の周波数と第2の格子の周波数との差からなるモアレ像の周波数成分を算出するものとできる。また、調整部を、モアレ周波数算出部によって算出された周波数成分がゼロに収束するような少なくとも一方の格子の配置ずれ量を取得するものとできる。また、第1の格子と第2の格子とを、タルボ干渉計を構成するように配置することができる。本発明の放射線位相画像撮影装置によれば、第1および第2の格子の少なくとも一方の予め設定された配置からの配置ずれによって形成されるモアレ像を放射線画像検出器によって検出し、そのモアレ像を表すモアレ画像信号に対してフーリエ変換処理を施してモアレ像の周波数成分を算出するようにしたので、格子の配置ずれをモアレ像の周波数成分によって定量的に扱うことができ、高精度な調整を行うことができる。また、フーリエ変換処理によって算出されたモアレ像の周波数成分に基づいて、第1および第2の格子の配置を調整するようにしたので、人間の視認限界によって目視観測することができない低周波側や高周波側のモアレ像についても、フーリエ変換処理後の周波数空間上で認識することができ、たとえば、その周波数成分に基づいて配置ずれ量を算出することによって第1および第2の格子を適切な配置に調整することができる。以下、図面を参照して本発明の放射線位相画像撮影装置の一実施形態について説明する。図1に本実施形態の放射線位相画像撮影装置の概略構成を示す。図2に図1に示す放射線位相画像撮影装置の上面図X-Z断面図を示す。図2の紙面厚さ方向が図1のY方向である。放射線位相画像撮影装置は、図1に示すように、放射線源1、第1の格子2、第2の格子3、放射線画像検出器4、移動機構5、位相画像取得部6、モアレ周波数算出部7、調整部8を備えている。放射線源1は、被写体10に向けて放射線を射出するものであり、第1の格子2に被写体10を透過した放射線を照射したとき、タルボ効果を発生させうるだけの空間的干渉性を有するものである。たとえば、放射線の発光点のサイズつまり、放射線源の開口径を30ミクロン程度とし、その放射線源から約5メートル以上の位置での空間的干渉性がそれに相当する。放射線源1としては、マイクロフォーカスX線管やプラズマX線源を利用することができる。また、比較的サイズの大きな通常の放射線源を用いる場合は、放射線が通過するマルチスリットを放射線の射出側に設置して使用することができる。この場合の詳細な構成は、たとえば、FranzPfeifferTimmWeikampOliverBunkChristianDavidNaturePhysics2258261LettersPhaseretrievalanddifferentialphasecontrastimagingwithlowbrillianceXraysourcesに記されているが、そのスリットのピッチP0は以下の式を満たすような大きさとする必要がある。なお、P2は第2の格子3のピッチ、Z1は放射線源1マルチスリットを用いる場合はマルチスリットの位置から第1の格子2までの距離、Z2は第1の格子2から第2の格子3までの距離である図2参照。第1の格子2は、図3に示すように、放射線を透過する基板21と、基板21上に設けられた複数の部材22とを備えている。複数の部材22は、放射線を遮蔽するものであり、いずれもX方向図3中紙面の厚さ方向に延びる線状形状で形成されている。複数の部材22どうしの間隔つまり、格子の周期P1は、本実施形態では一定とされている。部材22の素材としては、たとえば、金やシリコンを用いることができる。また、部材22としては、照射される放射線に対して約90°または約180°の位相変調を与える、いわゆる位相変調型格子を構成するものであることが望ましい。たとえば、通常の医療診断用のX線エネルギー領域において必要な金の厚さは1μm~数μm程度になる。また、振幅変調型格子を用いることもできる。この場合、部材22は放射線を十分に吸収する厚さが必要である。たとえば、通常の医療診断用のX線エネルギー領域において必要な金の厚さは10μm~数10μm程度になる。第2の格子3は、図4に示すように、第1の格子2と同様に、放射線を透過する基板31と、基板31に設けられた複数の部材32とを備えている。複数の部材32は、放射線を遮蔽するものであり、いずれもX方向図4中紙面の厚さ方向に延びる線状形状で形成されている。複数の部材32どうしの間隔つまり、格子の周期P2は、この実施形態では一定とされている。複数の部材32の素材としては、たとえば、金を用いることができる。第2の格子3については、部材32をより厚くした振幅変調型格子であることが望ましい。このとき、部材32は放射線を十分に吸収する厚さが必要である。たとえば、通常の医療診断用のX線エネルギー領域において必要な金の厚さは10μm~数10μm程度になる。ここで、放射線源1から照射される放射線が、平行ビームではなく、コーンビームである場合には、第1の格子2を通過して形成される第1の格子2の自己像は、放射線源1からの距離に比例して拡大される。そして、本実施形態においては、第2の格子3の格子ピッチP2と間隔d2は、そのスリット部が、第2の格子3の位置における第1の格子2の自己像の明部の周期パターンとほぼ一致するように決定される。すなわち、放射線源1の焦点から第1の格子2までの距離をZ1、第1の格子2から第2の格子3までの距離をZ2とした場合、第2の格子ピッチP2および間隔d2は、次式2および次式3の関係を満たすように決定される。なお、放射線源1から照射される放射線が平行ビームである場合には、P2=P1,d2=d1を満たすように決定される。放射線画像検出器4は、第1の格子2に入射した放射線が形成する第1の格子2の自己像を第2の格子3によって強度変調された画像信号として検出するものである。このような放射線画像検出器4としては、直接変換型および間接変換型のフラットパネル検出器、イメージングプレート、増感スクリーンとフィルムの組合せなど従来の放射線位相画像撮影装置に使われているものを用いることができる。そして、たとえばフラットパネル検出器を用いた場合には、図1に示すX-Y方向に2次元状に配列された各検出素子画素からそれぞれ画素信号が読み出され、イメージングプレートやフィルムの場合には、所定の読取装置によってX-Y方向に沿って走査されて読み出されることによって、イメージングプレートやフィルム上の各画素に対応する画素信号が取得される。そして、放射線源1、第1の格子2、第2の格子3および放射線画像検出器4によって放射線位相コントラスト画像を取得可能な放射線位相画像撮影装置が構成されるが、本構成をタルボ干渉計として機能させるためには、さらにいくつかの条件をほぼ満たさねばならない。その条件について以下に説明する。まず、第1の格子2と第2の格子3とのグリッド面が、図1に示すX-Y平面に平行であることが必要である。そして、さらに、第1の格子2と第2の格子3との距離Z2は、第1の格子2が90°の位相変調を与える位相変調型格子である場合、次の条件をほぼ満たさなければならない。ただし、λは放射線の波長通常はピーク波長、mは0か正の整数、P1は上述した第1の格子2の格子ピッチ、P2は上述した第2の格子3の格子ピッチである。また、第1の格子2が180°の位相変調を与える位相変調型格子である場合には、次の条件をほぼ満たさなければならない。ただし、λは放射線の波長通常はピーク波長、mは0か正の整数、P1は上述した第1の格子2の格子ピッチ、P2は上述した第2の格子3の格子ピッチである。さらに、第1の格子2が振幅変調型格子である場合には、次の条件をほぼ満たさなければならない。ただし、λは放射線の波長通常はピーク波長、mは正の整数、P1は上述した第1の格子2の格子ピッチ、P2は上述した第2の格子3の格子ピッチである。なお、上式4,5,6は、放射線源1により照射される放射線がコーンビームである場合であり、放射線が平行ビームである場合には、上式4に代えて下式7、上式5に代えて下式8、上式6に代えて下式9となる。移動機構5は、第2の格子3をY方向に移動させるものであるが、たとえば、第2の格子3の格子ピッチP2の1/nnは2以上の整数ずつ動かして、それぞれのnに対する位置で放射線画像の撮影を行なうことにより第2の格子3によって強度変調されたn種類n枚の画像信号を取得することができる。移動機構5は、たとえば公知のアクチュエーターを用いて構成することができる。位相画像取得部6は、放射線画像検出器4により検出されたn種類n枚の画像信号に基づいて放射線位相コントラスト画像を取得するものである。具体的には、n枚の画像信号からは、たとえば、KAStetsonandWRBrohinskyApplOct243631に記載の縞走査法を用いることにより、各画素ごとの位相シフトの微分量、すなわち被写体10によって生じた放射線の屈折角度に応じた量を復元することができ、様々な画像処理とその表現を介して、いわゆる位相コントラスト画像として取得することができる。モアレ周波数算出部7は、第1の格子2および第2の格子3の少なくとも一方が、上述した条件を満たすような予め設定された配置からずれることによって放射線画像検出器4により検出されるモアレ像を表すモアレ画像信号を取得し、その取得したモアレ画像信号に対して高速フーリエ変換処理を施してモアレ像の周波数成分を算出するものである。放射線画像検出器4により検出されるモアレ像については、後で詳述する。調整部8は、モアレ周波数算出部7によって算出されたモアレ像の周波数成分に基づいて、第1の格子2と第2の格子3の位置が、上述した条件を満たすような予め設定された配置となるように第1の格子2と第2の格子3の配置を変更するものである。すなわち調整部8は、第1の格子2または第2の格子3の配置ずれに基づいて放射線画像検出器4上に形成されるモアレ像の周波数成分が、低周波成分となるように第1の格子2または第2の格子3の配置を調整し、望ましくはモアレ像の周波数成分がゼロとなるように調整するものである。調整部8は、第1の格子2と第2の格子3の配置を変更可能なように機械的に構成されたものであり、公知のアクチュエーターを用いて構成することができる。調整部8による調整方法については、後で詳述する。次に、本実施形態の放射線位相画像撮影装置の作用について説明する。本実施形態の放射線位相画像撮影装置においては、被写体10の放射線位相コントラスト画像の撮影の前に、放射線画像検出器4上に形成されるモアレ像の周波数成分に基づいて、第1の格子2と第2の格子3の配置ずれ量を検出し、その配置ずれ量に基づいて調整部8によって第1の格子2と第2の格子3との配置を予め設定された配置となるように調整する。以下、第1の格子2と第2の格子3の配置ずれを調整する方法について具体的に説明する。まず、被写体10を配置していない状態で放射線源1から放射線が射出され、第1の格子2と第2の格子3を透過した放射線が放射線画像検出器4に照射され、放射線画像検出器4によって放射線画像信号が検出される。ここで、放射線源1から第1の格子2までの距離をZ1、第1の格子2から第2の格子3までの距離をZ2、第1の格子2の周波数周期P1の逆数をf1、第2の格子3の周波数周期P2をf2とする。上述したような放射線の照射によって第2の格子3上に射影される第1の格子2の自己像の周波数はf1’=f1×Z1Z1Z2となる。そして、第1の格子2と第2の格子3の配置が予め設定された正しい配置である場合に、周波数f1’が周波数f2と一致するように第1の格子2の周期P1と第2の格子の周期P2とが設定されている。そして、一般的に、2つの周波数パターンfと周波数パターンgとで形成されるモアレ像の周波数は、|f±g|であり、第1の格子2と第2の格子3との配置ずれによって放射線画像検出器4上に形成されるモアレ像の周波数もこれと同様に、fm=|f1’±f2|となる。ここで、f1’とf2との和によるモアレ像の周波数は、格子と同程度の周波数になるが、放射線画像検出器の画素サイズは、これよりも十分に大きいため、拡大率を考慮したとしても画像として観測できない上、放射線位相コントラスト像にも影響を与えない。なお、格子の周期と画素サイズとの関係は、たとえば、格子の周期5μm程度に対して画素サイズは100μm程度である。したがって、以下、f1’とf2との差によるモアレ像のみを取り上げて説明する。モアレ像の周波数は上述したようにfm=|f1’-f2|となるが、第1の格子2と第2の格子3との配置が完全に正しい配置である場合にはこの周波数はゼロになり、モアレ像は形成されないことになる。したがって、第1の格子2と第2の格子3との配置が調整され、最終的にモアレ像の周波数が十分に低周波数となるか、望ましくはゼロとなれば第1の格子2と第2の格子3とが正しい配置に調整されたことを意味する。そこで、本実施形態においては、放射線画像検出器4によってモアレ像を検出し、そのモアレ像を表すモアレ画像信号を用いて第1の格子2と第2の格子3の配置を調整する。なお、第1の格子2と第2の格子3のX方向またはY方向の並進ずれについては、放射線位相コントラスト画像に影響を与えないため、おおよそ配置があっていれば問題ない。したがって、本実施形態においては、第1の格子2と第2の格子3のZ方向の位置ずれと、図5Aに示すようなZ軸を回転軸とする回転ずれθzと、図5Bに示すようなX軸を回転軸とする回転ずれθxと、図5Cに示すようなY軸を回転軸とする回転ずれθyとを調整することになる。次に、放射線画像検出器4上に形成されるモアレ像について説明する。ここでは、説明をわかりやすくするため、第2の格子3を除いた放射線源1、第1の格子2、放射線画像検出器4などの位置決めを要する部材は、正しく配置されているものとする。このとき第2の格子3の配置がずれている場合にはモアレ像が形成されるが、たとえば、Z方向の位置ずれのみが生じていた場合には、図6に示すようなモアレ像が形成されることになる。また、たとえば、回転ずれθzのみが生じていた場合には、図7に示すようなモアレ像が形成されることになる。そして、このモアレ像は放射線画像検出器4によって検出され、モアレ像を表すモアレ画像信号が放射線画像検出器4からモアレ周波数算出部7に出力される。モアレ周波数算出部7は、入力されたモアレ画像信号に対して高速フーリエ変換処理を施し、モアレ画像信号の周波数空間上における周波数成分を算出する。具体的には、たとえば、図6に示すようなZ方向の位置ずれのみによるモアレ像の場合には、横方向周波数がないため、図8の黒点P1のように縦軸上に存在するピーク輝点が算出される。また、図7に示すような回転ずれθzのみによるモアレ像の場合には、斜め方向の一定周波数の縞であるので、図8の白点P2のように横方向の周波数成分と縦方向の周波数成分の両方をもったピークが算出される。そして、モアレ周波数算出部7で算出されたピークの周波数成分は調整部8に出力され、調整部8において、入力されたモアレ像の周波数成分に基づいて配置ずれ量が算出される。具体的には、たとえば、Z方向の位置ずれのみの場合には、モアレ像の周波数成分fmは、fm=f1×Z1Z1Z2-f2により算出されるので、モアレ周波数算出部7において算出されたモアレ像の周波数成分を代入することによってZ2を算出する。そして、予め設定された正しいZ2の値との差を求めることによってZ方向の配置ずれ量を算出する。そして、調整部8は、算出した配置ずれ量に基づいて第2の格子3の配置を調整して正しい配置とする。この調整部8による調整によってモアレ像の周波数を十分な低周波、望ましくはゼロにすることができる。また、回転ずれθzのみの場合には、調整部8は、入力されたモアレ像の縦方向の周波数成分と横方向の周波数成分との比率から回転ずれ量θzを算出する。そして、算出した回転ずれ量θzに基づいて第2の格子3をZ軸を回転軸として回転させて調整し、正しい配置とする。ただし、ここで算出した回転ずれ量θzは、第1の格子2のθzが正しいという前提での推定値である。なお、調整部8における配置ずれ量の算出方法としては、上記の方法に限らず、要するに、モアレ周波数算出部7で算出されたピークの周波数成分が、周波数空間上でゼロに近づくような算出方法であれば如何なる算出方法を採用してもよい。たとえば、モアレ周波数算出部7において算出されたモアレ像の周波数成分を示すピークの原点からの位置rと、第2の格子3のZ方向の所定の位置Z1,Z2,Z3とに基づいて、図9に示すようなこれらの相関関係を示す関数rZを求め、この関数rZがゼロに収束するようなZ方向の位置Z4を求めて配置ずれ量を算出するようにしてもよい。なお、回転ずれθxと回転ずれθyによるモアレ像の周波数成分については図示していないが、面内で連続的にZ方向の位置ずれが生じていることと同じであるので、θxやθyのずれが生じていると、周波数空間上では図8に示すような鋭いピークとならず、広がりを有する輝点の分布となる。したがって、周波数空間上でこの輝点の分布が小さくなるように第2の格子3の配置を調整するようにすればよい。実際には、輝点の分布の広がりを表現する量を算出し、この量がゼロに収束するようにθxないしθyの調整量を決定するようにすればよい。輝点の分布の広がりを表現する量としては、たとえば輝点のピークの半値幅などを用いればよい。なお、θxまたはθyの回転ずれによる輝点の幅は、その回転ずれによって生じたZ方向の位置ずれの最大振れ幅を反映したものとなるので、輝点の幅に基づいて調整量を算出することが可能である。具体的には、たとえば、輝点の分布が、図10に示すような分布となっている場合には、まず、その輝点の分布のピークを検出し、そのピーク近辺のプロファイルを縦方向と横方向とについてそれぞれ取得する。図11Aに縦方向のプロファイルの一例を、図11Bに横方向のプロファイルの一例を示す。そして、この縦方向のプロファイルの半値幅Wと、横方向のプロファイルの半値幅W’とをそれぞれ取得する。ここで、たとえば、θyについてΔθの角度ずれがある場合、図12に示すように、格子の中央領域でΔZ=±L4tanΔθ程度のZ方向の配置ずれが生じていることになる。そして、上述したようにモアレ周波数はfm=f1×Z1Z1Z2-f2により算出されるので、Z2をZ2±ΔZとしたときのモアレ周波数fmの振れ幅が、上述した半値幅Wに相当することになる。すなわち、上式のZ2をZ2±L4tanΔθとしたときのモアレ周波数fmの振れ幅がWに相当することになる。したがって、下式を演算することによってΔθを求めることができる。||f1×Z1{Z1Z2L4tanΔθ}-f2|-|f1×Z1{Z1Z2-L4tanΔθ}-f2||=Wなお、上記説明では、θyの角度ずれΔθの算出方法について説明したが、θxの角度ずれΔθ’を算出する場合には、格子の見かけ上の周波数が変わるので、モアレ周波数fmの式のf2をf2tanθに置き換え、そのモアレ周波数fmの振れ幅が半値幅W’に相当するものとして算出するようにすればよい。また、上記説明では、θxとθyのずれ量を算出するようにしたが、必ずしもこのようにθxとθyのずれ量を算出しなくてもよく、たとえば、調整部8によってθyを変化させながら、そのθyの変化に対する半値幅Wの変化量を取得し、その半値幅Wがゼロになる、もしくは最小値になるようなθyを正しい位置として設定するようにしてもよい。また、θxについても、調整部8によってθxを変化させながら、そのθxの変化に対する半値幅W’の変化量を取得し、その半値幅W’がゼロになる、もしくは最小値になるようなθxを正しい位置として設定するようにしてもよい。また、θxやθyの回転ずれがある場合、図13に示すように、面内で一様でない歪んだモアレ像が発生する。θxやθyの回転ずれは、このように実画像上ではモアレ像の歪みとして観測され、周波数空間上では輝点の広がりとして現れるので、このままでは上述したZ方向の位置ずれやθzの回転ずれを調整することが困難である。したがって、まず、上述したように輝点の分布の広がりを表現する量を算出し、この量がゼロに収束するようにθxないしθyを調整することによって面内で一様なモアレ像となるようにし、その後、Z方向の位置ずれとθzの回転ずれを調整するようにすればよい。また、本実施形態の説明においては、モアレ像の周波数成分に基づいて配置ずれ量を算出し、その配置ずれ量に基づいて第1の格子2または第2の格子3の配置を調整部8によって自動的に調整するようにしたが、これに限らず、たとえば、モアレ像の周波数成分をプロットした周波数空間をグラフとして所定の表示部に表示したり、モアレ像の周波数成分の数値を表示部に表示したりし、操作者がその表示されたグラフや数値を観察しながら、モアレ像の周波数成分がゼロに近づくように第1の格子2または第2の格子3を手動で調整するようにしてもよい。そして、上述したような第1の格子2と第2の格子3との配置の調整の後、実際の放射線位相コントラスト画像の撮影が行われる。具体的には、まず、図1に示すように、放射線源1と第1の格子2との間に、被写体10が配置される。次に、放射線源1から放射線が射出され、その放射線は被写体10を透過した後、第1の格子2に照射される。照射された放射線は、第1の格子2で回折されることにより、第1の格子2から放射線の光軸方向において所定の距離において、タルボ干渉像を形成する。これをタルボ効果と呼び、光波が格子を通過したとき、格子から所定の距離において、格子の自己像を形成する。たとえば、第1の格子2が、90°の位相変調を与える位相変調型格子の場合、上式4または上式7180°の位相変調型格子の場合は上式5または上式8、強度変調型格子の場合は上式6または上式9で与えられる距離において第1の格子2の自己像を形成する一方、被写体10によって、第1の格子2に入射する放射線の波面は歪むため、第1の格子2の自己像はそれに従って変形している。続いて、放射線は、第2の格子3を通過する。その結果、上記の変形した第1の格子2の自己像は第2の格子3との重ね合わせにより、強度変調を受け、上記波面の歪みを反映した画像信号として放射線画像検出器4により検出される。続いて順次、移動機構5により第2の格子3がY方向について第2の部材のピッチP2の1/nnは2以上の整数ずつ動かされ、放射線画像検出器4により、それぞれのnに対する波面の歪みを反映した画像信号が検出される。n種の波面の歪みを反映した画像信号は位相画像取得部6に入力され、位相画像取得部6は、n種の波面の歪みを反映した画像を解析することにより、放射線位相コントラスト画像を生成する。波面の歪みは、被写体10による屈折効果によって放射線が曲げられた角度に比例している。したがって、放射線画像検出器4により検出されたn種類の波面の歪みを反映した画像から各画素での位相変調量を解析することにより、被写体10の内部の屈折率分布に依存した量を検出することができる。また、上記実施形態の放射線位相画像撮影装置は、第1の格子2から第2の格子3までの距離Z2がタルボ干渉距離となるようにしたが、これに限らず、第1の格子2が入射放射線を回折せずに投影させる構成とするようにしてもよい。このように構成すれば第1の格子2を通過して射影される投影像が、第1の格子2の後方の全ての位置で相似的に得られるため、第1の格子2から第2の格子3までの距離Z2を、タルボ干渉距離を無関係に設定することができる。具体的には、第1の格子2と第2の格子3とを、ともに吸収型振幅変調型格子として構成するとともに、タルボ干渉効果の有無に関わらず、スリット部を通過した放射線を幾何学的に投影するように構成する。より詳細には、第1の格子2の間隔d1と第2の格子3の間隔d2とを、放射線源1から照射される放射線のピーク波長より十分大きな値とすることで、照射放射線に含まれる大部分をスリット部で回折せずに、直進性を保ったまま通過するように構成することができる。たとえば、放射線源のターゲットとしてタングステンを用い、管電圧を50kVとした場合には、放射線のピーク波長は約0.4Åである。この場合には、第1の格子2の間隔d1と第2の格子3の間隔d2を、1μm~10μm程度とすればスリット部で大部分の放射線が回折されずに幾何学的に投影される。なお、第1の格子2の格子ピッチP1と第2の格子3の格子ピッチP2との関係と、第1の格子2の間隔d1と第2の格子3の間隔d2との関係とについては、上記第1の実施形態と同様である。そして、上記のような構成の放射線位相画像撮影装置においては、第1の格子2と第2の格子3との距離Z2を、上式6においてm=1とした場合の最小のタルボ干渉距離より短い値に設定することができる。すなわち、上記距離Z2が、次式10を満たす範囲の値に設定する。なお、第1の格子2の部材22と第2の格子3の部材32とは、コントラストの高い周期パターン像を生成するためには、放射線を完全に遮蔽吸収することが好ましいが、上述した放射線吸収に優れる材料金、白金等を用いたとしても、吸収されずに透過する放射線が少なからず存在する。このため、放射線の遮蔽性を高めるためには、部材22,32のそれぞれの厚みを、可能な限り厚くすることが好ましい。部材22,32による遮蔽は、照射放射線の90以上であることが好ましく、たとえば、放射線源1の管電圧が50kVの場合には、厚みは、金Au換算で30μm以上であることが好ましい。上記のような構成の放射線位相画像撮影装置によれば、第1の格子2と第2の格子3との距離Z2をタルボ干渉距離よりも短くすることができるので、一定のタルボ干渉距離を確保しなければならない上記実施形態の放射線位相画像撮影装置と比較すると、撮影装置をより薄型化することができる。また、上記実施形態の放射線位相画像撮影装置においては、移動機構5によって第2の格子3を並進移動させるとともに、複数回の撮影を行うことによって、放射線位相コントラスト画像を生成するための複数種類の画像信号を取得するようにしたが、このように第2の格子3を並進移動させることなく、一回の撮影によって複数種類の画像信号を取得する方法もある。具体的には、図14に示すように、第1の格子2と第2の格子3とが、第1の格子2の延伸方向と第2の格子3の延伸方向とが相対的に傾くように配置されるようにする。そして、このように配置された第1の格子2と第3の格子3に対して、放射線画像検出器4によって検出される画像信号の各画素の主走査方向図14のX方向の主画素サイズDxと副走査方向の副画素サイズDyとが、図14に示すような関係となるようにする。主画素サイズDxは、たとえば、放射線画像検出器として、多数の線状電極を有し、その線状電極の延伸方向に直交する方向に延設された線状読取光源によって走査されて画像信号が読み出される、いわゆる光読取方式の放射線画像検出器を用いた場合には、その放射線画像検出器の線状電極の配列ピッチによって決定されるものである。また、副画素サイズDyは、線状読取光源によって放射線画像検出器に照射される線状の読取光の幅によって決定されるものである。また、いわゆるTFT読取方式の放射線画像検出器やCMOSを用いた放射線画像検出器を用いた場合には、主画素サイズDxは、画像信号が読み出されるデータ電極の配列方向の画素回路の配列ピッチによって決定され、副画素サイズDyは、ゲート電圧が出力されるゲート電極の配列方向の画素回路の配列ピッチによって決定される。そして、放射線位相コントラスト画像を生成するための複数種類の画像の数をMとすると、M個の副画素サイズDyが位相コントラスト画像の副走査方向の1つの画像解像度Dとなるように第1の格子2が第2の格子3に対して傾けられる。具体的には、図15に示すように、第2の格子3のピッチおよび第1の格子2によって第2の格子3の位置に形成される第1の格子2の自己像G1のピッチをp、第2の格子3に対する第1の格子2の自己像のX-Y面内の相対的な回転角をθ、位相コントラスト画像の副走査方向の画像解像度をD=Dy×Mとすると、回転角θを下式11を満たすように設定することによって、副走査方向の画像解像度Dの長さに対して、第1の格子2の自己像G1と第2の格子3の位相がn周期分ずれることになる。なお、図15においては、M=5、n=1の場合を示している。したがって、放射線位相コントラスト画像の副走査方向の画像解像度DをM分割したDx×Dyの各画素によって、第1の格子2の自己像のn周期分の強度変調をM分割した画像信号が検出できることになる。図15に示す例では、n=1としているので、副走査方向の画像解像度Dの長さに対して、第1の格子2の自己像G1と第2の格子3の位相が1周期分ずれることになる。もっとわかり易く言えば、第1の格子2の自己像G1の1周期分の第2の格子3を通過する範囲が、副走査方向の画像解像度Dの長さにわたって変化する。そして、M=5としているので、Dx×Dyの各画素によって第1の格子2の自己像の1周期の強度変調を5分割した画像信号が検出できることになり、すなわち、Dx×Dyの各画素によって互いに異なる5つの画像の画像信号をそれぞれ検出することができることになる。なお、位相コントラスト画像の主走査方向の画像解像度Dxと副走査方向の画像解像度D=Dy×Mとは同じになるようしてもよいし、これに限らず任意の主副比としてもよい。さらに、本実施形態においては、M=5としているが、Mは3以上であればよく、5以外であってもよい。また、上記説明ではn=1としたが、nは0以外の整数であれば1以外の整数でもよい。すなわち、nが負の整数の場合には上述した例に対して反対周りの回転となり、また、nを±1以外の整数としてn周期分の強度変調としてもよい。ただし、nがMの倍数の場合は、1組のM個の副走査方向画素Dyの間で第1の格子2の自己像G1と第2の格子3の位相が等しくなり、異なるM個の縞画像とならないため除外するものとする。また、第2の格子3に対する第1の格子2の自己像の回転角θの調整については、たとえば、放射線画像検出器4と第2の格子3の相対回転角を固定した後、第1の格子2を回転させることによって行うことができる。たとえば、上式11でp=5μm、D=50μm、n=1とすると、理論上の回転角θは約5.7°である。そして、第2の格子3に対する第1の格子2の自己像の実際の回転角θ’は、たとえば、第1の格子の自己像と第2の格子3によるモアレのピッチによって検出することができる。具体的には、図16に示すように、実際の回転角をθ’、回転によって生じたX方向への見た目の自己像のピッチP’とすると、観測されるモアレのピッチPmは、1/Pm=|1/P’-1/P|であるので、P’=P/cosθ’を上式に代入することによって実際の回転角θ’を求めることができる。なお、モアレのピッチPmについては、放射線画像検出器4によって検出された画像信号に基づいて求めるようにすればよい。そして、理論上の回転角θと実際の回転角θ’とを比較し、その差の分だけで自動または手動で第1の格子2の回転角を調整するようにすればよい。そして、上記のように構成された放射線位相画像撮影装置においては、放射線画像検出器4から読み出された1フレーム全体の画像信号が位相コントラスト画像生成部61に記憶された後、その記憶された画像信号に基づいて、互いに異なる5つの画像の画像信号が取得される具体的には、図15に示すように、放射線位相コントラスト画像の副走査方向の画像解像度Dを5分割し、第1の格子2の自己像の1周期の強度変調を5分割した画像信号が検出できるように第1の格子2を第2の格子3に対して傾けるようにした場合には、図17に示すように、第1読取ラインから読み出された画像信号が第1の縞画像信号M1として取得され、第2読取ラインから読み出された画像信号が第2の縞画像信号M2として取得され、第3読取ラインから読み出された画像信号が第3の縞画像信号M3として取得され、第4読取ラインから読み出された画像信号が第4の縞画像信号M4として取得され、第5読取ラインから読み出された画像信号が第5の縞画像信号M5として取得される。なお、図17に示す第1~第5読取ラインは、図15に示す副画素サイズDyに相当する。また、図17においては、Dx×Dy×5の読取範囲しか示していないが、その他の読取範囲についても、上記と同様にして第1~第5の縞画像信号が取得される。すなわち、図18に示すように、副走査方向について4画素間隔毎の画素行読取ラインからなる画素行群の画像信号が取得されて1フレームの1つの縞画像信号が取得される。より具体的には、第1読取ラインの画素行群の画像信号が取得されて1フレームの第1の縞画像信号が取得され、第2読取ラインの画素行群の画像信号が取得されて1フレームの第2の縞画像信号が取得され、第3読取ラインの画素行群の画像信号が取得されて1フレームの第3の縞画像信号が取得され、第4読取ラインの画素行群の画像信号が取得されて1フレームの第4の縞画像信号が取得され、第5読取ラインの画素行群の画像信号が取得されて1フレームの第5の縞画像信号が取得される。なお、ここでは、画像解像度Dを5分割する場合について説明したが、M分割する場合には、図14に示すDy1,Dy1M,Dy12M,Dy13M,,,の系列、Dy2,Dy2M,Dy22M,Dy23M,,,の系列、Dy3,Dy3M,Dy32M,Dy33M,,,の系列などがそれぞれ1フレームの縞画像信号として取得されることになる。そして、上記第1~第5の縞画像信号に基づいて、位相画像取得部6において放射線位相コントラスト画像が生成される。ここで、上述したような1回の撮影によって複数の縞画像信号を取得する実施形態において、第1および第2の格子23と放射線画像検出器4との相対的な位置ずれを調整する方法について以下に説明する。まず、上記実施形態と同様に、放射線位相コントラスト画像の撮影前に、被写体10を配置していない状態で放射線源1から放射線が射出されてプレ曝射が行われ、第1の格子2と第2の格子3との相対的な位置ずれによって生じたモアレ像の撮影が行われる。なお、この配置ずれによるモアレ像の周波数成分は、上述したように第1の格子2と第2の格子3とが相対的に傾くように予め設置することによって生じたモアレ像の周波数成分と異なるものである。そして、放射線画像検出器4によって検出されたモアレ像の画像は、モアレ周波数算出部7に入力される。モアレ周波数算出部7は、位相画像取得部6において複数の縞画像信号を取得したときと同様に、モアレ像の画像について、副走査方向にM画素毎のDyの系列をまとめる。たとえば、図14に示すようなDy1,Dy1M,Dy12M,Dy13M,,,の系列や、Dy2,Dy2M,Dy22M,Dy23M,,,の系列である。つまり、画像解像度Dを単位として同じ位置のDyを1つの系列にまとめる。そして、これらの系列の少なくとも1つの系列に高速フーリエ変換処理を施すことによってモアレ像の周波数成分を算出し、その周波数成分の空間分布を取得する。なお、モアレ像の周波数成分を複数の系列毎に取得し、さらに正確な位置ずれを算出するようにしてもよい。そして、上記実施形態においては、第1の格子2の自己像の延伸方向と第2の格子3の延伸方向とが平行になるように設置するようにしたので、モアレ像の周波数成分が最終的に十分に低周波数となるか、望ましくはゼロとなるように第1の格子2または第2の格子3の位置を調整するようにしたが、上述したように第1の格子2と第2の格子3とを相対的に傾くように予め設置するようにした場合には、上述したプレ曝射によって取得されたモアレ像の周波数成分が、予め設定されたモアレ像の周波数成分に近づくか、望ましくはその周波数成分に等しくなるように第1の格子2または第2の格子3の位置が調整されることになる。この予め設定されたモアレ像の周波数成分とは、上述したように第1の格子2と第2の格子3とを相対的に傾くように予め設置することによって生じたモアレ像の周波数成分である。なお、第1の格子2または第2の格子3の配置ずれ量の算出方法については、基本的には上記実施形態と同様であり、最終的に調整されるモアレ像の周波数成分を変更することによって算出するようにすればよい。また、上記説明では、図14に示すように、第1の格子2の延伸方向と第2の格子3の延伸方向とが相対的に傾くように配置された状態で撮影された1枚の画像から、互いに異なる画素行群の画像信号を取得することによって複数の縞画像信号を取得し、その複数の縞画像信号を用いて放射線位相コントラスト画像を生成するようにしたが、上記のようにして撮影された1枚の画像に基づいて複数の縞画像信号を生成するのではなく、上記のようにして撮影した1枚の画像に対してフーリエ変換を施すことによっても放射線位相コントラスト画像を生成することができ、このような方法を採用してもよい。具体的には、まず、第1の格子2の延伸方向と第2の格子3の延伸方向とが相対的に傾くように配置された状態で撮影された1枚の画像に対してフーリエ変換処理を施すことによって、その画像に含まれる被検体による吸収情報と位相情報とを分離する。そして、周波数空間上において被検体による位相情報の部分のみを抽出して周波数空間上の中心原点位置に移動した後、その抽出した位相情報に対して逆フーリエ変換処理を施し、各画素に対して、その結果の虚部を実部で除算したものの逆正接関数arctan虚部/実部を演算することによって、被検体による放射線の屈折角を求めることができる。そして、この屈折角に基づいて被検体の位相シフト分布の微分量、すなわち、位相微分像を取得することができる。なお、上述のフーリエ変換を用いた位相コントラスト画像の生成方法では、第1の格子2の延伸方向と第2の格子3の延伸方向とが相対的に傾くように配置された状態で撮影された1枚の画像を用いることとしたが、このように第1の格子2と第2の格子3とを相対的に傾けた場合に限らず、第1の格子2の自己像と第2の格子3の重ね合わせによってモアレを生じさせ、該モアレが検出されている少なくとも1枚の画像縞画像を用いてもよい。"