"X線を照射するX線源と、前記X線の照射軸方向と直交する方向に複数のスリットが配列されて構成された第1格子及び第2格子と、被写体台と、照射されたX線に応じて電気信号を生成する変換素子が2次元状に配置され、当該変換素子により生成された電気信号を画像信号として読み取るX線検出器と、を有する撮影装置を備え、前記X線源により照射され、前記第1格子、前記被写体台に載置された被写体、及び前記第2格子を透過したX線を前記X線検出器で読み取ることにより得られたモアレ画像に基づいて、被写体の再構成画像を生成するX線撮影システムであって、前記被写体台に被写体を載置して撮影することにより取得された少なくとも一つの被写体有りのモアレ画像及び前記被写体台に前記被写体を載置せずに撮影することにより取得された少なくとも一つの被写体無しのモアレ画像に基づいて、前記被写体の再構成画像を作成する再構成画像作成手段と、前記作成された前記被写体の再構成画像に対して、前記被写体有りのモアレ画像を撮影したときと前記被写体無しのモアレ画像を撮影したときの撮影条件の差異に起因するアーチファクトの補正を行う再構成画像補正手段と、備えるX線撮影システム。前記撮影条件は、前記第1格子、前記第2格子の配置である請求項1に記載のX線撮影システム。前記撮影装置は、前記第1格子と前記第2格子を一定周期間隔で相対移動させ、前記第1格子と前記第2格子の相対移動毎に、前記X線源により照射されたX線に応じて前記X線検出器が画像信号の読み取る処理を繰り返すことにより複数のモアレ画像を取得する請求項1又は2に記載のX線撮影システム。前記撮影装置は、更に、前記X線の照射軸方向と直交する方向に複数のスリットが配列されて構成されたマルチスリットを備え、前記マルチスリットを前記X線の照射方向と直交する方向に移動させ、前記マルチスリットが一定周期間隔で移動する毎に、前記X線源により照射されたX線に応じて前記X線検出器が画像信号の読み取る処理を繰り返すことにより複数のモアレ画像を取得する請求項1又は2に記載のX線撮影システム。前記再構成画像作成手段は、前記被写体有りのモアレ画像、及び前記被写体無しのモアレ画像からフーリエ変換法を用いて前記被写体の診断用再構成画像を作成する請求項1又は2に記載のX線撮影システム。前記第1格子及び前記第2格子は一次元格子である請求項1~5の何れか一項に記載のX線撮影システム。前記第1格子及び前記第2格子は二次元格子である請求項5に記載のX線撮影システム。ところで、被写体有りのモアレ画像と被写体無しのモアレ画像を同時に撮影することはできないため、それぞれ別個に撮影を行うが、例えば、撮影室の空調が両者の撮影の間で稼働しはじめたりすると、これらの撮影間で撮影条件が変化し、再構成画像に撮影条件の変化に起因したアーチファクトが出てしまい、再構成画像のコントラストが悪く見えてしまうという問題がある。本発明の課題は、被写体有りの撮影と被写体無しの撮影との間で撮影条件が変化したことにより生じる画質の低下を抑制し、診断に良好な被写体再構成画像を提供できるようにすることである。上記課題を解決するため、請求項1に記載の発明は、X線を照射するX線源と、前記X線の照射軸方向と直交する方向に複数のスリットが配列されて構成された第1格子及び第2格子と、被写体台と、照射されたX線に応じて電気信号を生成する変換素子が2次元状に配置され、当該変換素子により生成された電気信号を画像信号として読み取るX線検出器と、を有する撮影装置を備え、前記X線源により照射され、前記第1格子、前記被写体台に載置された被写体、及び前記第2格子を透過したX線を前記X線検出器で読み取ることにより得られたモアレ画像に基づいて、被写体の再構成画像を生成するX線撮影システムであって、前記被写体台に被写体を載置して撮影することにより取得された少なくとも一つの被写体有りのモアレ画像及び前記被写体台に前記被写体を載置せずに撮影することにより取得された少なくとも一つの被写体無しのモアレ画像に基づいて、前記被写体の再構成画像を作成する再構成画像作成手段と、前記作成された前記被写体の再構成画像に対して、前記被写体有りのモアレ画像を撮影したときと前記被写体無しのモアレ画像を撮影したときの撮影条件の差異に起因するアーチファクトの補正を行う再構成画像補正手段と、備える。請求項2に記載の発明は、請求項1に記載の発明において、前記撮影条件は、前記第1格子、前記第2格子の配置である。請求項3に記載の発明は、請求項1又は2に記載の発明において、前記撮影装置は、前記第1格子と前記第2格子を一定周期間隔で相対移動させ、前記第1格子と前記第2格子の相対移動毎に、前記X線源により照射されたX線に応じて前記X線検出器が画像信号の読み取る処理を繰り返すことにより複数のモアレ画像を取得する。請求項4に記載の発明は、請求項1又は2に記載の発明において、前記撮影装置は、更に、前記X線の照射軸方向と直交する方向に複数のスリットが配列されて構成されたマルチスリットを備え、前記マルチスリットを前記X線の照射方向と直交する方向に移動させ、前記マルチスリットが一定周期間隔で移動する毎に、前記X線源により照射されたX線に応じて前記X線検出器が画像信号の読み取る処理を繰り返すことにより複数のモアレ画像を取得する。請求項5に記載の発明は、請求項1又は2に記載の発明において、前記再構成画像作成手段は、前記被写体有りのモアレ画像、及び前記被写体無しのモアレ画像からフーリエ変換法を用いて前記被写体の診断用再構成画像を作成する。請求項6に記載の発明は、請求項1~5の何れか一項に記載の発明において、前記第1格子及び前記第2格子は一次元格子である。請求項7に記載の発明は、請求項5に記載の発明において、前記第1格子及び前記第2格子は二次元格子である。本発明によれば、被写体有りの撮影と被写体無しの撮影との間で撮影条件が変化したことにより生じる画質の低下を抑制し、診断に良好な被写体再構成画像を提供することが可能となる。以下、図面を参照して本発明の実施形態について説明する。図1に、本実施形態に係るX線撮影システムを示す。X線撮影システムは、X線撮影装置1とコントローラ5を備える。X線撮影装置1はタルボロー干渉計によるX線撮影を行い、コントローラ5は当該X線撮影により得られたモアレ画像を用いて被写体の再構成画像を作成する。本実施形態においては、X線撮影装置1は、手指を被写体として撮影する装置として説明するが、これに限定されるものではない。X線撮影装置1は、図1に示すように、X線源11、マルチスリット12、被写体台13、第1格子14、第2格子15、X線検出器16、保持部17、本体部18等を備える。X線撮影装置1は縦型であり、X線源11、マルチスリット12、被写体台13、第1格子14、第2格子15、X線検出器16は、この順序に重力方向であるz方向に配置される。X線源11の焦点とマルチスリット12間の距離をd1mm、X線源11の焦点とX線検出器16間の距離をd2mm、マルチスリット12と第1格子14間の距離をd3mm、第1格子14と第2格子15間の距離をd4mmで表す。距離d1は好ましくは5~500mmであり、さらに好ましくは5~300mmである。距離d2は、一般的に放射線科の撮影室の高さは3m程度又はそれ以下であることから、少なくとも3000mm以下であることが好ましい。なかでも、距離d2は400~2500mmが好ましく、さらに好ましくは500~2000mmである。X線源11の焦点と第1格子14間の距離d1d3は、好ましくは300~5000mmであり、さらに好ましくは400~1800mmである。X線源11の焦点と第2格子15間の距離d1d3d4は、好ましくは400~5000mmであり、さらに好ましくは500~2000mmである。それぞれの距離は、X線源11から照射されるX線の波長から、第2格子15上に第1格子14による格子像自己像が重なる最適な距離を算出し、設定すればよい。X線源11、マルチスリット12、被写体台13、第1格子14、第2格子15、X線検出器16は、同一の保持部17に一体的に保持され、z方向における位置関係が固定されている。保持部17はC型のアーム状に形成され、本体部18に設けられた駆動部18aによりz方向に移動昇降可能に本体部18に取り付けられている。X線源11は、緩衝部材17aを介して保持されている。緩衝部材17aは、衝撃や振動を吸収できる材料であれば何れの材料を用いてもよいが、例えばエラストマー等が挙げられる。X線源11はX線の照射によって発熱するため、X線源11側の緩衝部材17aは加えて断熱素材であることが好ましい。X線源11はX線管を備え、当該X線管によりX線を発生させてz方向重力方向にX線を照射する。X線管としては、例えば医療現場で広く一般に用いられているクーリッジX線管や回転陽極X線管を用いることができる。陽極としては、タングステンやモリブデンを用いることができる。X線の焦点径は、0.03~3mmが好ましく、さらに好ましくは0.1~1mmである。マルチスリット12は回折格子であり、図2Aに示すように複数のスリットが所定間隔で配列されて設けられている。この複数のスリットは、X線照射軸方向図1のz方向と直交する方向図2Aに白矢印で示すに配列されている。マルチスリット12はシリコンやガラスといったX線の吸収率が低い材質の基板上に、タングステン、鉛、金といったX線の遮蔽力が大きい、つまりX線の吸収率が高い材質により形成される。例えば、フォトリソグラフィーによりレジスト層がスリット状にマスクされ、UVが照射されてスリットのパターンがレジスト層に転写される。露光によって当該パターンと同じ形状のスリット構造が得られ、電鋳法によりスリット構造間に金属が埋め込まれて、マルチスリット12が形成される。マルチスリット12のスリット周期は1~60μmである。スリット周期は、図2Aに示すように隣接するスリット間の距離を1周期とする。スリットの幅各スリットのスリット配列方向の長さはスリット周期の1~60の長さであり、さらに好ましくは10~40である。スリットの高さz方向の高さは1~500μmであり、好ましくは1~150μmである。マルチスリット12のスリット周期をw0μm、第1格子14のスリット周期をw1μmとすると、スリット周期w0は下記式により求めることができる。w0=w1d3d4/d4当該式を満たすように周期w0を決定することにより、マルチスリット12及び第1格子14の各スリットを通過したX線により形成される自己像が、それぞれ第2格子15上で重なり合い、いわばピントが合った状態とすることができる。マルチスリット12は、図2Bに示すように、ラック12aを有するホルダー12bに保持されている。ラック12aは、マルチスリット12のスリット配列方向に設けられている。ラック12aは、後述する駆動部122のピニオン122cと係合し、ピニオン122cの回転位相角に応じてホルダー12bに保持されたマルチスリット12をスリット配列方向に移動させるためのものである。本実施形態において、X線撮影装置1には、マルチスリット回転部121及び駆動部122が設けられている。マルチスリット回転部121は、第1格子14及び第2格子15のX線照射軸周りの回転に応じてマルチスリット12をX線照射軸周りに回転させるための機構である。駆動部122は、複数のモアレ画像の撮影のためにマルチスリット12をスリット配列方向に移動させるための機構である。図3に、マルチスリット回転部121及び駆動部122の平面図及びA-A ́断面図を示す。図3に示すように、マルチスリット回転部121は、モータ部121a、ギア部121b、ギア部121c、支持部121d等を備えて構成されている。モータ部121a、ギア部121b、ギア部121cは、支持部121dを介して保持部17に保持されている。モータ部121aは、マイクロステップ駆動に切り替え可能なパルスモータであり、制御部181図8参照からの制御に応じて駆動され、ギア部121bを介してギア部121cをX線照射軸図3に一点鎖線Rで示すを中心として回転させる。ギア部121cは、ホルダー12bに保持されたマルチスリット12を装着するための開口部121eを有している。ギア部121cを回転させることにより、開口部121eに装着されたマルチスリット12をX線照射軸周りに回転させ、マルチスリット12のスリット配列方向を可変することができる。なお、撮影において、マルチスリット12は0°~90°程度回転できればよいので、ギア部121cは全周にある必要はなく、図3に2点鎖線で示す範囲正逆回転方向にそれぞれ90°で回転できればよい。開口部121eは、ホルダー12bに保持されたマルチスリット12を上部から嵌め込むことが可能な形状及びサイズとなっている。ここでは、開口部121eにおけるスリット配列方向のサイズW4はホルダー12bにおけるスリット配列方向のサイズW2より若干大きくなっており、マルチスリット12をスリット配列方向にスライドさせることが可能となっている。なお、開口部121eにおけるスリット配列方向に直交する方向のサイズW3は、ホルダー12bにおけるスリット配列方向に直交する方向のサイズW1との精密嵌合可能な寸法としており、ホルダー12bを開口部121eに装着すると、ホルダー12bに設けられたラック12aは開口部121eの外に、後述するピニオン122cと係合可能に配置される。駆動部122は、マルチスリット周期に応じ0.1μm~数十μm単位でマルチスリット12をスリット配列方向に移動させる精密減速機等を備えて構成される。駆動部122は、例えば、図3に示すように、モータ部122a、ギア部122b、ピニオン122c等を備えて構成され、図示しないL字型板金等によりマルチスリット回転部121のギア部121cに固定されている。これにより、マルチスリット12と駆動部122は一体的に回転されるようになっている。モータ部122aは、例えば、制御部181からの制御に応じて駆動され、ギア部122bを介してピニオン122cを回転させる。ピニオン122cは、マルチスリット12のラック12aと係合して回転することで、マルチスリット12をスリット配列方向に移動させる。図1に戻り、被写体台13は、被写体である手指を載置するための台である。被写体台13は、患者の肘が載置できる高さに設けられていることが好ましい。このように、患者の肘まで載置できるように構成することで、患者は楽な姿勢となり、比較的長時間にわたる撮影の間に、指先の撮影部位の動きを低減させることができる。また、被写体台13には、被写体を固定するための被写体ホルダー130が設けられている。図4Aに示すように、被写体ホルダー130は、手のひらで掴みやすいマウスのような楕円形状131のついた板状の部材である。上記楕円形状131は、その断面を側面から観察すると、図4Bに示すように、手のひらサイズのなだらかな凸曲面となっており、患者が手のひらで楕円形状131を掴むことで、被写体が疲れにくい状態で被写体の下方への動きを抑制することができる。被写体ホルダー130が場所によってX線複素屈折率の不均一な形状又は厚みを有している場合、X線検出器16に到達するX線量は被写体ホルダー130のX線複素屈折率が不均一であることによってムラが生じる。被写体ホルダー130上には、更に被写体姿勢を安定させるため、指間スペーサ133を備えることが好ましい。また、患者毎に手や指間の大きさは異なるので、患者毎の手のひらの形状に合わせて被写体ホルダー130を作成しておき、撮影時には、その患者用の被写体ホルダー130を被写体台13にマグネット等で取り付けることが好ましい。腕から手首までの荷重は被写体台13が支えるので、被写体ホルダー130は指先部分の加重と患者が上方から押さえる力に耐えるものであればよく、安価で量産が可能なプラスチック成形とすることが可能である。図1に戻り、第1格子14は、マルチスリット12と同様にX線照射軸方向であるz方向と直交する方向に複数のスリットが配列されて設けられた回折格子である。第1格子14は、マルチスリット12と同様にUVを用いたフォトリソグラフィーによって形成することもできるし、いわゆるICP法によりシリコン基板に微細細線で深掘加工を行い、シリコンのみで格子構造を形成することとしてもよい。第1格子14のスリット周期は1~20μmである。スリットの幅はスリット周期の20~70であり、好ましくは35~60である。スリットの高さは1~100μmである。第1格子14として位相型を用いる場合、スリットの高さz方向の高さはスリット周期を形成する2種の素材、つまりX線透過部とX線遮蔽部の素材による位相差がπ/8~15×π/8となる高さとする。好ましくは、π/4~3×π/4となる高さである。第1格子14として吸収型を用いる場合、スリットの高さはX線遮蔽部によりX線が十分吸収される高さとする。第1格子14が位相型である場合、第1格子14と第2格子15間の距離d4は、次の条件をほぼ満たすことが必要である。d4=m1/2w12/λなお、mは整数であり、λはX線の波長である。上記条件は第1格子がπ/2型格子であること、つまり、第1格子のX線遮蔽部とX線透過部の素材による位相差がπ/2近傍にある場合の例を述べたが、第1格子はπ型を用いても良く、使用する格子の型式に応じた条件を演算すればよい。第2格子15は、マルチスリット12と同様に、X線照射軸方向であるz方向と直交する方向に複数のスリットが配列されて設けられた回折格子である。第2格子15もフォトリソグラフィーにより形成することができる。第2格子15のスリット周期は1~20μmである。スリットの幅はスリット周期の30~70であり、好ましくは35~60である。スリットの高さは1~100μmである。本実施形態では第1格子14及び第2格子15は、それぞれの格子面がz方向に対し垂直x-y平面内で平行であり、第1格子14のスリット配列方向と第2格子15のスリット配列方向とは、x-y平面内で所定角度だけ傾けて配置されているが、両者を平行な配置としても良い。また、本実施形態において、第1格子14及び第2格子15は、円盤状である。上記マルチスリット12、第1格子14、第2格子15は、例えば下記のように構成することができる。X線源11のX線管の焦点径;300μm、管電圧:40kVp、付加フィルタ:アルミ1.6mmX線源11の焦点からマルチスリット12までの距離d1:240mmマルチスリット12から第1格子14までの距離d3:1110mmマルチスリット12から第2格子15までの距離d3d4:1370mmマルチスリット12のサイズ:10mm四方、スリット周期:22.8μm第1格子14のサイズ:50mm四方、スリット周期:4.3μm第2格子15のサイズ:50mm四方、スリット周期:5.3μm本実施形態において、第1格子14及び第2格子15は格子回転部210に装着されている。図5に、格子回転部210の平面図を示す。図6に、第1格子14及び第2格子15が装着された状態の格子回転部210の平面図及びD-D ́断面図を示す。図5に示すように、格子回転部210は、回転トレイ212上に、ハンドル211、相対角調整部213、ストッパー214等を備えて構成されている。回転トレイ212は、第1格子14及び第2格子15を保持するための開口部212aを有している。ここで、本実施形態において、第1格子14は、複数のスリットが配列されてなる円形状の格子部140と、この格子部140を開口部212aに取り付けるための第1ホルダー部141及び第2ホルダー部142と、により構成されている図6参照。第1ホルダー部141は、格子部140の外周に取り付けられた、開口部212aと略同じ半径外周の半径の部材であり、第1格子14の装着時に開口部212aと嵌合する。第2ホルダー部142は、第1ホルダー部141より更に外側に取り付けられた、開口部212aよりやや半径外周の半径の大きい部材である。第2ホルダー部142は、その外周の一部がギア加工されている。また、第2ホルダー部142の外周の所定位置には、突起部142aが設けられている。第2格子15は、複数のスリットが配列されてなる円形状の格子部150と、この格子部150を開口部212aに取り付けるためのホルダー部151と、により構成されている。ホルダー部151は、開口部212aの半径と略同様の半径を有する円盤状の部材である。格子部150はホルダー部151の中央部上面に保持されている図6参照。回転トレイ212に第1格子14及び第2格子15を装着する際には、まず、第2格子15が開口部212aの底面に嵌め込まれる。次に、第1格子14が第2格子15の上から開口部212aに嵌め込まれる。これにより、図6に示す状態で回転トレイ212に第1格子14及び第2格子15が保持される。開口部212aに保持された第1格子14及び第2格子15は、撮影モードに応じて相対角調整部213により互いのスリット方向の相対角度が調整される。ここで、X線撮影装置1は、縞走査法による再構成画像用に複数ステップで撮影を行う第1の撮影モードと、フーリエ変換法による再構成画像用に1又は2方向で撮影を行う第2の撮影モードとを有している。縞走査法用の撮影において要求される第1格子14のスリット方向と第2格子15のスリット方向との相対角度は第2格子の周期と画像サイズと縞の本数に依存する。縞走査法においては、モアレ画像における干渉縞の本数が少ないほど、また、干渉縞が鮮明なほど、このモアレ画像に基づいて作成される再構成画像が鮮明となることが知られている山田朝治、横関俊介編著モアレ縞干渉縞応用計測法、コロナ社、1996年12月10日参照。そこで、第2格子の周期を5.3μm、60mm角の画像の中に干渉縞が0~3本程度にすると仮定すると、相対角度は0度~±0.015度とする必要がある。一方、フーリエ変換法用の撮影において要求される第1格子14のスリット方向と第2格子15のスリット方向との相対角度はX線検出器16の画素ピッチと空間分解能に依存する。一般的に使用される検出器空間分解能30μm~200μmとすると、その相対角度は0.4度~3度とする必要がある。よって、第1の撮影モードと第2の撮影モードを切り替えて撮影を行うためには、撮影モードに応じて第1格子と第2格子の相対角度を調整する必要がある。しかし、例えば、縞走査法では0.00005度縞1本分の誤差~0.00015度縞3本分の誤差の精度での調整が要求されるため、手動で第1格子14及び第2格子15のスリット方向の相対角度を調整するのは困難である。そこで、X線撮影装置1においては、相対角調整部213により、操作部182で設定された撮影モードに応じて第1格子14と第2格子15の相対角度を自動的に調整することが可能となっている。図5、6に示すように、相対角調整部213は、モータ部213a、第1ギア213b、第2ギア213c、レバー213dにより構成されている。モータ部213aは、第2ギア213cと係合し、制御部181からの制御に応じて第2ギア213cを回転させる。第2ギア213cは、レバー213dを介してその中心が第1ギア213bの中心と接続されており、その円周が第1ギア213bと係合している。第2ギア213cがモータ部213aの駆動に応じて回転すると、第1ギア213bが第2ギア213cの中心を回転軸として第2ギア213cの周囲に沿って回転移動して第1格子14の第2ホルダー142のギア部分と係合し、第2格子15を回転させることなく第1格子14をX線照射軸周りに回転させることが可能となっている。本実施形態においては、工場出荷時に、第2ホルダー部142の突起部142aが回転トレイ212上に設けられたストッパー凸状の突起214に突き当たったときに第1格子14と第2格子15のスリット方向の相対角度が第1の撮影モード時縞走査法用の撮影モードに最適な相対角となるように、ストッパー214の位置及び第1格子14と第2格子15の相対角度が予め調整されて開口部212aに装着されている。第2の撮影モードフーリエ変換用の撮影モードが設定されると、第1格子14と第2格子の相対角度が第2の撮影モードに最適となるように制御部181により相対角調整部213のパルスモータを採用したモータ部213aが駆動通電制御される。これにより、第2ギア213cを介して第1ギア213bが回転して第2ホルダー部142のギア部分に係合し、第1格子14と第2格子15のスリット方向の相対角度が第2の撮影モードに最適となるように第2ホルダー部142が回転される。その後、当該パルスモータへの通電状態を可変させ、後述するバネ力に打ち勝つモータ自己保持力励磁力を発揮させる程度の通電状態駆動時の定格電流の50未満等とすることで、第2ホルダー部142を当該位相に維持することができる。なお、このときの回転角度は1度程度と小さいため、まず、モータ部213aのパルスモータにより、突起部142aが基準位置215にくるまで反時計周りに第2ホルダー部142を回転させ、突起部142aが基準位置215に来たことが図示しないセンサにより検知されたら、第2ホルダー部142の回転方向を時計回りに切り替えてマイクロステップ駆動より第2ホルダー部142を回転させることが好ましい。第2ホルダー部142は図示しないバネに付勢されており、モータ部213aの駆動により第1ギア213bと第2ホルダー部142との係合が解除されると、バネの付勢力により突起部142aがストッパー214の位置に戻る。即ち、第1格子14と第2格子15が第1の撮影モードに最適な相対角度に戻る。以上のようにして、第1格子14と第2格子15の相対角度が撮影モードに応じた角度に調整される。格子回転部210は、また、相対角度が調整された第1格子14及び第2格子15を被写体に対してX線照射軸周り図6に点線Rで示すに一体的に回転させることができる。ここで、1次元格子スリットを用いるタルボ干渉計及びタルボロー干渉計では、第1格子14及び第2格子15のスリット方向と平行に線状に延びる構造物は鮮明に撮影することができないという特性がある。よって、被写体の注目すべき構造物の配置方向に応じて、第1格子14及び第2格子15のスリット方向の角度を調整する必要がある。格子回転部210は、以下の機構により第1格子14及び第2格子15をその相対角度を維持したまま一体的にX線照射軸周りに回転させ、被写体の注目すべき構造物の配置方向に対する第1格子14及び第2格子15のスリット方向の角度を調整することができる。回転トレイ212には、上述のように、ハンドル211が設けられている。ハンドル211は、撮影技師等のオペレータがX線照射軸図6に点線Rで示すを軸として回転トレイ212を手動で回転させるための突起である。また、回転トレイ212は、回転トレイ212の回転角度を固定するための凹部212b~212eを有している。凹部212b~212eは、予め0°と定められた位置ここでは凹部212bがトレイ固定部材171bのボールと対向する位置を0°の位置とするから所定の回転角度にある位置ここでは、0°、30°、60°、90°に設けられている。凹部212b~212eのそれぞれには、角度検知センサSE1~SE4が設けられており、トレイ固定部材171bと係合したことを検知して制御部181にその検知信号を出力する。このように、回転トレイ212を手動で回転させるので、患者が触れる範囲に第1格子14及び第2格子15を一体的に回転させるための電気コード等を設ける必要がなく、安全性を確保することができる。なお、本実施形態では、回転トレイ212が0°に設定されたときの第1格子14及び第2格子15の位置角度をホームポジションとする。また、第1格子14及び第2格子15がホームポジションであるときの第1格子14のスリット方向とマルチスリット12のスリット方向が平行である位置角度をマルチスリット12のホームポジションとする。図7Aは、保持部17における格子回転部210の保持部分171を拡大して示した平面図であり、図7Bは、図7AにおけるE-E ́断面図である。図7Cは、保持部17に格子回転部210を保持した状態を示す図である。図7A、図7Bに示すように、保持部分171には、格子回転部210の回転トレイ212と精密嵌合するサイズであり、回転トレイ212を回転可能に保持する開口部171aと、回転トレイ212の回転角度を固定するためのトレイ固定部材171bと、が設けられている。開口部171aの底部とX線検出器16の載置部の間は、X線の透過を妨げないように、中空とするか又はX線透過率の高いアルミやカーボン等とすることが好ましい。トレイ固定部材171bは、凹部212b~212eの何れかがトレイ固定部材171bと対向するように位置したときにその対向する凹部に係合するボールと、ボールを図7A、図7Bの矢印方向に誘導するための図示しないスライドガイド押圧バネのガイドにより構成されている。凹部212b~212eの何れかがトレイ固定部材171bと対向する位置で回転トレイ212の回転が停止すると、トレイ固定部材171bのスライドガイドにより、対向している凹部にボールが係合するとともに、凹部に設けられた角度検知センサSE1~SE4の何れかによりボールの係合が検知されて制御部181に検知信号が出力される。これにより、制御部181は、回転トレイ212の回転角度、即ち、第1格子14及び第2格子15の回転角度を検知できるようになっている。また、図7Dに示すように、回転トレイ212の開口部212aの下部に、X線検出器16の装着部212fを設け、第1格子14及び第2格子15とX線検出器16とを一体として回転させることができるようにしてもよい。このようにすれば、X線検出器16の縦横方向の鮮鋭性の異方性の影響を受けることがないので、再構成画像の縦横の鮮鋭性を第1格子14及び第2格子15の回転角度によらずに概ね一定とすることができる。図1に戻り、X線検出器16は、照射されたX線に応じて電気信号を生成する変換素子が2次元状に配置され、当該変換素子により生成された電気信号を画像信号として読み取る。X線検出器16の画素サイズは10~300μmであり、さらに好ましくは50~200μmである。X線検出器16は第2格子15に当接するように保持部17に位置を固定することが好ましい。第2格子15とX線検出器16間の距離が大きくなるほど、X線検出器16により得られるモアレ画像がボケるからである。X線検出器16としては、FPDFlatPanelDetectorを用いることができる。FPDにはX線をシンチレータを介して光電変換素子により電気信号に変換する間接変換型、X線を直接的に電気信号に変換する直接変換型があるが、何れを用いてもよい。間接変換型は、CsIやGd2O2等のシンチレータプレートの下に、光電変換素子がTFT薄膜トランジスタとともに2次元状に配置されて各画素を構成する。X線検出器16に入射したX線がシンチレータプレートに吸収されると、シンチレータプレートが発光する。この発光した光により、各光電変換素子に電荷が蓄積され、蓄積された電荷は画像信号として読み出される。直接変換型は、アモルファスセレンの熱蒸着により、100~1000μmの膜圧のアモルファスセレン膜がガラス上に形成され、2次元状に配置されたTFTのアレイ上にアモルファスセレン膜と電極が蒸着される。アモルファスセレン膜がX線を吸収するとき、電子正孔対の形で物質内に電圧が遊離され、電極間の電圧信号がTFTにより読み取られる。なお、CCDChargeCoupledDevice、X線カメラ等の撮影手段をX線検出器16として用いてもよい。X線撮影時のFPDによる一連の処理を説明する。まずFPDはリセットを行い、前回の撮影読取以降に残存する不要な電荷を取り除く。その後、X線の照射が開始するタイミングで電荷の蓄積が行われ、X線の照射が終了するタイミングで蓄積された電荷が画像信号として読み取られる。なお、リセットの直後や画像信号の読み取り後等に、オフセット補正用のダーク読み取りを少なくとも1回行う。本体部18は、図8に示すように、制御部181、操作部182、表示部183、通信部184、記憶部185を備えて構成されている。制御部181は、CPUCentralProcessingUnitやRAMRandomAccessMemory等から構成され、記憶部185に記憶されているプログラムとの協働により、X線撮影装置1の各部を制御するとともに各種処理を実行する。例えば、制御部181は、後述する撮影制御処理をはじめとする各種処理を実行する。操作部182は曝射スイッチや撮影条件等の入力操作に用いるキー群の他、表示部183のディスプレイと一体に構成されたタッチパネルを備え、これらの操作に応じた操作信号を生成して制御部181に出力する。表示部183は制御部181の表示制御に従って、ディスプレイに操作画面やX線撮影装置1の動作状況等を表示する。通信部184は通信インターフェイスを備え、ネットワーク上のコントローラ5と通信する。例えば、通信部184はX線検出器16によって読み取られ、記憶部185に記憶されたモアレ画像をコントローラ5に送信する。記憶部185は、制御部181により実行されるプログラム、プログラムの実行に必要なデータを記憶している。また、記憶部185はX線検出器16によって得られたモアレ画像を記憶する。コントローラ5は、オペレータによる操作に従ってX線撮影装置1の撮影動作を制御する。また、コントローラ5は、X線撮影装置1により得られたモアレ画像を用いて診断用の被写体再構成画像を作成する画像処理部として機能する。コントローラ5は、図9に示すように、制御部51、操作部52、表示部53、通信部54、記憶部55を備えて構成されている。制御部51は、CPUCentralProcessingUnitやRAMRandomAccessMemory等から構成され、記憶部55に記憶されているプログラムとの協働により、後述する縞走査法による再構成画像作成処理、フーリエ変換法による再構成画像作成処理をはじめとする各種処理を実行する。操作部52は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードで押下操作されたキーの押下信号とマウスによる操作信号とを、入力信号として制御部51に出力する。表示部53のディスプレイと一体に構成されたタッチパネルを備え、これらの操作に応じた操作信号を生成して制御部51に出力する構成としてもよい。表示部53は、例えば、CRTCathodeRayTubeやLCDLiquidCrystalDisplay等のモニタを備えて構成されており、制御部51の表示制御に従って、操作画面、X線撮影装置1の動作状況、作成された被写体再構成画像等を表示する。通信部54は、通信インターフェイスを備え、ネットワーク上のX線撮影装置1やX線検出器16と有線又は無線により通信する。例えば、通信部54は、X線撮影装置1に撮影条件や制御信号を送信したり、X線撮影装置1又はX線検出器16からモアレ画像を受信したりする。記憶部55は、制御部51により実行されるプログラム、プログラムの実行に必要なデータを記憶している。例えば、記憶部55は、RIS、HIS等や図示しない予約装置より予約されたオーダを示す撮影オーダ情報を記憶している。撮影オーダ情報は、患者名、撮影部位、撮影モード等の情報である。記憶部55は、X線検出器16によって得られたモアレ画像、モアレ画像に基づき作成された診断用の被写体再構成画像を撮影オーダ情報に対応付けて記憶する。また、記憶部55は、X線検出器16に対応するゲイン補正データ、欠陥画素マップ等を予め記憶する。コントローラ5においては、操作部52の操作により撮影オーダ情報の一覧表示が指示されると、制御部51により、記憶部55から撮影オーダ情報が読み出されて表示部53に表示される。操作部52により撮影オーダ情報が指定されると、指定された撮影オーダ情報に応じた撮影条件撮影モードを含むの設定情報やX線源11のウォームアップの指示等が通信部54によりX線撮影装置1に送信される。これにより、X線撮影装置1に撮影モードが設定される。即ち、コントローラ5は、撮影モードを設定する設定部として機能する。また、X線検出器16がケーブルレスのカセッテ型FPD装置である場合には、制御部51は、内部バッテリ消耗防止の為のスリープ状態から、撮影可能状態に起動せしめる。X線撮影装置1においては、通信部184によりコントローラ5から撮影条件の設定情報等が受信されると、X線撮影準備が実行される。上記X線撮影装置1のタルボロー干渉計によるX線撮影方法第1の撮影モードの撮影方法を説明する。図10に示すように、X線源11から照射されたX線が第1格子14を透過すると、透過したX線がz方向に一定の間隔で像を結ぶ。この像を自己像といい、自己像が形成される現象をタルボ効果という。自己像を結ぶ位置に第2格子15が平行に配置され、当該第2格子15はその格子方向が第1格子14の格子方向と平行な位置からわずかに傾けられているので、第2格子15を透過したX線によりモアレ画像Mが得られる。X線源11と第1格子14間に被写体Hが存在すると、被写体HによってX線の位相がずれるため、図10に示すようにモアレ画像M上の干渉縞は被写体Hの辺縁を境界に乱れる。この干渉縞の乱れを、モアレ画像Mを処理することによって検出し、被写体像を画像化することができる。これがタルボ干渉計及びタルボロー干渉計の原理である。X線撮影装置1では、X線源11と第1格子14との間のX線源11に近い位置に、マルチスリット12が配置され、タルボロー干渉計によるX線撮影が行われる。タルボ干渉計はX線源11が理想的な点線源であることを前提としているが、実際の撮影にはある程度焦点径が大きい焦点が用いられるため、マルチスリット12によってあたかも点線源が複数連なってX線が照射されているかのように多光源化する。これがタルボロー干渉計によるX線撮影法であり、焦点径がある程度大きい場合にも、タルボ干渉計と同様のタルボ効果を得ることができる。従来のタルボロー干渉計では、マルチスリット12は上述のように多光源化と照射線量の増大を目的に用いられ、複数のモアレ画像を得るためには第1格子14又は第2格子15を相対移動させていた。しかし、本実施形態では、第1格子14又は第2格子15を相対移動させるのではなく、第1格子14及び第2格子15の位置は固定したまま、第1格子14及び第2格子15に対してマルチスリット12を移動させることで一定周期間隔のモアレ画像を複数得る。なお、第2の撮影モードによりモアレ画像を得る場合には、マルチスリット12の移動は行わず、1回撮影又は被写体とスリット方向を90度回転させて2回撮影を行う。図11は、X線撮影装置1の制御部181により実行される撮影制御処理を示すフローチャートである。撮影制御処理は、制御部181と記憶部185に記憶されているプログラムの協働により実行される。まず、コントローラ5から受信された設定情報に基づいて、第1の撮影モード縞走査法用と第2の撮影モードフーリエ変換法用の何れの撮影モードが設定されているかが判断されるステップS1。第1の撮影モードが設定されていると判断されるとステップS1;第1の撮影モード、第1の撮影モード処理が実行されるステップS2。一方、第2の撮影モードが設定されていると判断されるとステップS1;第2の撮影モード、第2の撮影モード処理が実行されるステップS3。図12は、図11のステップS2においてX線撮影装置1の制御部181により実行される第1の撮影モード処理を示すフローチャートである。第1の撮影モード処理は、制御部181と記憶部185に記憶されているプログラムの協働により実行される。ここで、第1の撮影モードでのX線撮影には上述のタルボロー干渉計によるX線撮影方法が用いられ、被写体像の再構成には縞走査法が用いられる。X線撮影装置1では、制御部181の制御により駆動部122が駆動及び停止されることによりマルチスリット12が等間隔毎に複数ステップ移動され、ステップ毎に撮影が行われて、各ステップのモアレ画像が得られる。ステップ数は2~20、さらに好ましくは3~10である。視認性の高い再構成画像を短時間で得るという観点からすれば、5ステップが好ましい参照文献1KHibinoBFOrebandDIFarrantPhaseshiftingfornonsinusoidalwaveformswithphaseshifterrorsJOptSocAmAVol12761768、参照文献2AMomoseWYashiroYTakedaYSuzukiandTHattoriPhaseTomographybyXrayTalbotInterferometetryforbiologicalimagingJpnJApplPhysVol4552545262。図12に示すように、まず、制御部181により、X線源11がウォームアップ状態に切り替えられるステップS101。次いで、格子回転部210の相対角調整部213が制御され、第1格子14と第2格子15の相対角度が第1の撮影モードに最適となるように突起部142aがストッパー214に接触する位置にくるように第1格子14が回転される。これにより、第1格子14と第2格子15の相対角度が調整されるステップS102。次いで、オペレータの操作に応じて第1格子14及び第2格子15が一体的に回転され、被写体に対する第1格子14及び第2格子15のスリット方向が設定されるステップS103。即ち、撮影技師等のオペレータは、格子回転部210のハンドル211を回転させ、被写体台13に載置された被写体の注目すべき構造物の配置方向に応じて第1格子14及び第2格子15のスリット方向を設定する。ハンドル211の回転が停止し、トレイ固定部材171bのバネ附勢されたボールの係合により位置固定されると、角度検知センサSE1~SE4の何れかから制御部181に検知信号が出力され、制御部181において、設定されたスリット方向に対応する、格子回転部210の回転トレイ212即ち、第1格子14及び第2格子15のホームポジションからの回転角度が取得される。次いで、第1格子14及び第2格子15の回転角度に応じて、マルチスリット回転部121のモータ部121aがパルスにより制御され、第1格子14及び第2格子15の回転角度に応じてマルチスリット12が回転されるステップS104。例えば、モータ部121aのパルスモータが制御され、マルチスリット12のホームポジションからの回転角度が回転トレイ212回転角度近傍例えば、回転トレイ212が30°に設定された場合は29°ぐらいまで一気に回転される。次いで、モータ部121aがマイクロステップ精密制御に切り替えられ、マルチスリット12を少しずつ回転させながら複数の回転角度で撮影が行われて調整用の複数のモアレ画像が生成されるステップS105。例えば、回転トレイ212が30°に設定された場合は、マルチスリット12を29.5°、30°、30.5°の3つの回転角度に設定してX線源11により低線量X線が照射され、撮影が行われる。これにより、調整用の3つのモアレ画像が取得される。なお、ステップS105においては、被写体を被写体台13に載置しない状態で撮影が行われる。撮影された調整用の複数のモアレ画像は、マルチスリット12の回転角度に対応付けて、並べて表示部183に表示されるステップS106。ここで、上述のように、第1格子14と第2格子15の相対角度は干渉縞本数が最小となるようにステップS102で調整されているので、ステップS103においては、回転トレイ212の回転によりその相対角度を保ったまま第1格子14及び第2格子15が回転される。しかし、第1格子14及び第2格子15を載置した回転トレイ212が回転し、マルチスリット12と第1格子14及び第2格子15との相対角度が変化すると、干渉縞すなわちモアレの鮮明性が変化してしまう。そこで、マルチスリット12と第1格子14及び第2格子15、即ち、これらを載置した回転トレイ212との相対角度を調整する必要がある。一般的には、マルチスリット12と第1格子14との相対角度が少ないほど、干渉縞の鮮明性の高いモアレ画像が得られる。しかし、マルチスリット12は発熱部であるX線源11近傍に配置されるので熱影響を受けやすい。そのため、マルチスリット12の変形等を考慮して、マルチスリット12を回転トレイ212と同じ角度だけ回転させるだけでなく、モータ部121aをマイクロステップ駆動させてステップS105~S108における微調整を行うことが有効である。オペレータは、ステップS106で表示部183に表示されたモアレ画像を観察し、干渉縞が最も鮮明な回転角度を撮影に用いる回転角度として選択する。なお、ここでは、干渉縞の鮮明性はオペレータの目視により観察するが、干渉縞の鮮明性の度合いを示す鮮明度は、後述するサインカーブ図17参照における極大値をMAX、極小値をMINとした場合、下記の式で表すことができる。この干渉縞の鮮明度を用いてオペレータではなく、プログラムで自動的に最大値となる回転角度を設定しても良い。干渉縞の鮮明度=MAX-MIN/MAXMIN=振幅/平均値操作部182により、マルチスリット12の回転角度が入力されるとステップS107;YES、モータ部121aが再駆動され、マルチスリット12のホームポジションからの回転角度が入力された回転角度となるようにマルチスリット12の位置が微調整されるステップS108。マルチスリット12の回転角度の調整後、被写体台13に被写体が載置され、オペレータにより曝射スイッチがON操作されるとステップS109;YES、駆動部122によりマルチスリット12がそのスリット配列方向に移動され、複数ステップの撮影が実行され、被写体有りの複数のモアレ画像が生成されるステップS110。まず、マルチスリット12が停止した状態でX線源11によるX線の照射が開始される。X線検出器16ではリセット後、X線照射のタイミングに合わせて電荷が蓄積され、X線の照射停止のタイミングに合わせて蓄積された電荷が画像信号として読み取られる。これが1ステップ分の撮影である。1ステップ分の撮影が終了するタイミングで、制御部181の制御により駆動部122が起動され、マルチスリット12の移動が開始される。所定量移動すると駆動部122が停止されることによりマルチスリットの移動が停止され、次のステップの撮影が行われる。このようにして、マルチスリット12の移動と停止が所定のステップ数分だけ繰り返され、マルチスリット12が停止したときにX線の照射と画像信号の読み取りが行われる。読み取られた画像信号はモアレ画像として本体部18に出力される。例えば、マルチスリット12のスリット周期を22.8μmとし、5ステップの撮影を10秒で行うとする。マルチスリット12がそのスリット周期の1/5に該当する4.56μm移動し停止する毎に撮影が行われる。従来のように第2格子15又は第1格子14を移動させる場合、第2格子15のスリット周期は比較的小さく、各ステップの移動量も小さくなるが、マルチスリット12のスリット周期は第2格子15よりも比較的大きく、各ステップの移動量も大きい。例えば、スリット周期5.3μmの第2格子15のステップ毎の移動量は1.06μmであるのに対し、スリット周期22.8μmのマルチスリット12の移動量は4.56μmと約4倍の大きさである。同一の駆動伝達系駆動源、減速伝達系を含むを使用し、各ステップの撮影に際し、駆動部122の起動と停止を繰り返して撮影を行った場合、移動用のパルスモータ駆動源の制御量駆動パルス数に対応した実際の移動量に占める、起動時及び停止時の駆動部122のバックラッシュ等の影響による移動量誤差の割合は、本実施形態のようにマルチスリット12を移動させる方式の方が小さくなる。これは、後述するサインカーブに沿ったモアレ画像を得やすく、起動及び停止を繰り返しても高精細な再構成画像が得られることを示している。或いは、従来方式による画像でも充分診断に適合する場合には、モータ駆動源を含む駆動伝達系全体の精度特に、起動特性及び停止特性を緩和し、駆動伝達系を構成する部品のコストダウンが可能であることを示している。各ステップの撮影が終了すると、本体部18の通信部184からコントローラ5に、各ステップのモアレ画像が送信されるステップS111。本体部18からコントローラ5に対しては各ステップの撮影が終了する毎に1枚ずつ被写体有りのモアレ画像が送信される。次いで、X線検出器16においてダーク読み取りが行われ、被写体有り画像データ補正用のダーク画像オフセット補正データが取得されるステップS112。ダーク読み取りは、少なくとも1回行われる。又は、複数回のダーク読み取りを行ってその平均値をダーク画像として取得してもよい。ダーク画像は、通信部184からコントローラ5に送信されるステップS113。当該ダーク読取に基づくオフセット補正データは、各モアレ画像信号の補正に共通に用いられる。尚、ダーク画像の取得は、各ステップのモアレ画像取得後に、当該ステップのダーク読取を行って、各ステップ専用のオフセット補正データを生成することとしても良い。次いで、オペレータによる曝射スイッチのON操作待ち状態となるステップS114。ここで、オペレータは、被写体無しのモアレ画像を作成できるように、被写体台13から被写体を取り除いて患者を退避させる。被写体なしの撮影の準備が完了したら、曝射スイッチを押下する。曝射スイッチが押下されるとステップS114;YES、駆動部122によりマルチスリット12がそのスリット配列方向に移動され、被写体なしで複数ステップの撮影が実行され、被写体無しの複数のモアレ画像が生成されるステップS115。各ステップの撮影が終了すると、本体部18の通信部184からコントローラ5に、各ステップのモアレ画像が送信されるステップS116。本体部18からコントローラ5に対しては各ステップの撮影が終了する毎に通信部184により1枚ずつ被写体無しのモアレ画像が送信される。次いで、X線検出器16においてダーク読み取りが行われ、被写体無しのダーク画像が取得されるステップS117。ダーク読み取りは、少なくとも1回行われる。又は、複数回のダーク読み取りを行ってその平均値をダーク画像として取得してもよい。ダーク画像は、通信部184からコントローラ5に送信されステップS118、一つの撮影オーダに対する一連の撮影は終了する。尚、ダーク画像の取得は、各ステップのモアレ画像取得後に、当該ステップのダーク読取を行って、各ステップ専用のオフセット補正データを生成することとしても良い。なお、被写体無しの複数のモアレ画像及びダーク読取は被写体有りの撮影直後に実施されることが精度的に最も好ましいが、被写体画像の再構成までの時間短縮をはかるために、始業開始等に事前に取得しておいたデータを用いることとしても良い。コントローラ5においては、通信部54によりモアレ画像が受信されると、受信されたモアレ画像が撮影開始時に指定された撮影オーダ情報と対応付けて記憶部55に記憶される。図13は、図11のステップS3においてX線撮影装置1の制御部181により実行される第2の撮影モード処理を示すフローチャートである。第2の撮影モード処理は、制御部181と記憶部185に記憶されているプログラムの協働により実行される。図13に示すように、まず、制御部181により、X線源11がウォームアップ状態に切り替えられるステップS201。次いで、格子回転部210の相対角調整部213が制御され、第1格子14と第2格子15の相対角度が第2の撮影モードに最適となるように突起部142aがホームポジションから所定角度回転した位置にくるように調整されるステップS202。次いで、ステップS203~ステップS208の処理が行われる。ステップS203~S208の処理は、図12のステップS103~108で説明したものと同様であるので説明を援用する。被写体台13に被写体が載置され、オペレータにより曝射スイッチがON操作されるとステップS209;YES、撮影が実行され、被写体有りのモアレ画像が生成されるステップS210。即ち、X線源11から放射線が照射され、X線検出器16において読み取りが行われる。なお、第2の撮影モードでは、駆動部122は停止されたままの状態で、マルチスリット12を移動させることなく、1枚のみの撮影が行われる。撮影が終了すると、本体部18の通信部184からコントローラ5に撮影により得られたモアレ画像が送信されるステップS211。次いで、X線検出器16においてダーク読み取りが行われ、被写体有り画像データ補正用のダーク画像オフセット補正データが取得されるステップS212。ダーク読み取りは、少なくとも1回行われる。又は、複数回のダーク読み取りを行ってその平均値をダーク画像として取得してもよい。ダーク画像は、通信部184からコントローラ5に送信されるステップS213。当該ダーク読取に基づくオフセット補正データは、モアレ画像信号の補正に共通に用いられる。次いで、オペレータによる曝射スイッチのON操作待ち状態となるステップS214。ここで、オペレータは、被写体無しのモアレ画像を作成できるように、被写体台13から被写体を取り除いて患者を退避させる。被写体なしの撮影の準備が完了したら、曝射スイッチを押下する。曝射スイッチが押下されるとステップS214;YES、被写体無しで撮影が実行され、被写体無しのモアレ画像が生成されるステップS215。ステップS215においてもステップS210と同様に、駆動部122は停止されたままの状態で、マルチスリット12を移動させることなく、1枚のみの撮影が行われる。撮影が終了すると、本体部18の通信部184からコントローラ5に、モアレ画像が送信されるステップS216。次いで、X線検出器16においてダーク読み取りが行われ、被写体無しのダーク画像が取得されるステップS217。ダーク読み取りは、少なくとも1回行われる。又は、複数回のダーク読み取りを行ってその平均値をダーク画像として取得してもよい。ダーク画像は、通信部184からコントローラ5に送信されステップS218、一つの撮影オーダに対する一連の撮影は終了する。なお、被写体無しの複数のモアレ画像及びダーク読取は被写体有りの撮影直後に実施されることが精度的に最も好ましいが、被写体画像の再構成までの時間短縮をはかるために、始業開始等に事前に取得しておいたデータを用いることとしても良い。コントローラ5の制御部51においては、通信部54によりモアレ画像が受信されると、現在処理対象となっている撮影オーダ情報において設定されている撮影モードが第1の撮影モードである場合は縞走査法による再構成画像作成処理が実行され、第2の撮影モードである場合はフーリエ変換法による再構成画像作成処理が実行される。図14Aは、制御部51により実行される縞走査法による再構成画像作成処理を示すフローチャートである。縞走査法による再構成画像作成処理は、制御部51と記憶部55に記憶されているプログラムとの協働により実行される。まず、ステップS11~S13においては、被写体有りの複数のモアレ画像及び被写体なしの複数のモアレ画像のそれぞれについて、X線検出器16の各画素のバラツキを補正するための補正処理が実行される。具体的には、オフセット補正処理ステップS11、ゲイン補正処理ステップS12、欠陥画素補正処理ステップS13が実行される。ステップS11においては、被写体有り画像データ補正用のダーク画像に基づいて、被写体有りの各モアレ画像にオフセット補正が施される。被写体無し画像データ補正用のダーク画像に基づいて、被写体無しの各モアレ画像にオフセット処理が施される。ステップS12においては、撮影に用いられたX線検出器16に対応するゲイン補正データが記憶部55から読み出され、読み出されたゲイン補正データに基づいて、各モアレ画像にゲイン補正が施される。ステップS13においては、撮影に用いられたX線検出器16に対応する欠陥画素マップ欠陥画素位置を示すデータが記憶部55から読み出され、各モアレ画像における欠陥画素位置マップで示す位置の画素値信号値が周辺画素により補間算出される。次いで、複数のモアレ画像間でX線強度変動補正トレンド補正が行われるステップS14。縞走査法では、複数のモアレ画像に基づいて1枚の被写体再構成画像が作成される。そのため、各モアレ画像の撮影において照射されるX線強度にゆらぎ変動があると精巧な被写体再構成画像が得られず、微細な信号の変化が見落とされてしまう可能性がある。そこで、ステップS14においては、複数のモアレ画像における撮影時のX線強度変動による信号値差を補正する処理が行われる。具体的な処理としては、各モアレ画像の予め定められた1点の画素の信号値を用いて補正する方法、各モアレ画像間におけるX線検出器16の所定方向の信号値差を補正する一次元補正する方法、各モアレ画像間における2次元方向の信号値差を補正する二次元補正する方法、の何れであってもよい。1点の画素の信号値を用いて補正する方法では、まず、図15に示すように複数のモアレ画像のそれぞれについて、X線検出器16のモアレ画像領域被写体配置領域161外の直接X線領域に対応する予め定められた位置Pにある画素の信号値が取得される。次いで、1枚目のモアレ画像例えば、被写体有りの最初に撮影されたモアレ画像が2枚目以降の上記取得された位置Pの画素の平均信号値で規格化され、規格化後の位置Pの値に基づいて2枚目以降の各モアレ画像の補正係数が算出される。そして、2枚目以降の各モアレ画像に補正係数が乗算されることにより、X線強度変動が補正される。この補正方法では、各撮影間の全体的なX線強度の変動を容易に補正することができる。なお、X線検出器16の裏側に、X線照射量を検知するセンサ等の検知手段を設け、検知手段から出力される各モアレ画像撮影時のX線照射量に基づいて、各モアレ画像間における撮影時のX線強度変動に起因する信号値差を補正することとしてもよい。一次元補正では、まず、複数のモアレ画像のそれぞれについて、予め定められた行L1行は、X線検出器16における読み取りライン方向をさすの画素の平均信号値が算出される。次いで、1枚目のモアレ画像が2枚目以降の画素の平均信号値で規格化され、規格化後の行L1と2枚目以降の行L1の各画素の信号値に基づいて、2枚目以降の各モアレ画像の行方向の補正係数が算出される。そして、2枚目以降の各モアレ画像に行方向の位置に応じた補正係数が乗算されることにより、行方向のX線強度変動が補正される。この補正方法では、各撮影間の一次元方向のX線強度の変動を容易に補正することができる。例えば、或る撮影において、X線源11による照射タイミングとX線検出器16の読み取りタイミングのずれが生じた場合に、これにより生じるX線検出器16の読み取りライン方向のX線強度変動等を補正することができる。二次元補正では、まず、複数のモアレ画像のそれぞれについて、予め定められた行L1、列L2列は、X線検出器16における読み取りライン方向と直交する方向をさすのそれぞれにおける画素の平均信号値が算出される。次いで、1枚目のモアレ画像が2枚目以降の行L1の画素の平均信号値で規格化され、規格化後の行L1と2枚目以降の行L1の各画素の信号値に基づいて、2枚目以降の各モアレ画像の行方向の補正係数が算出される。同様に、1枚目のモアレ画像が2枚目以降の列L2の画素の平均信号値で規格化され、規格化後の列L2と2枚目以降の列L2の各画素の信号値に基づいて、2枚目以降の各モアレ画像の列方向の補正係数が算出される。そして行方向と列方向の補正係数が掛け合わされて2枚目以降の各モアレ画像の各画素の補正係数が算出される。そして、各画素に行方向及び列方向の補正係数が乗算されることにより、二次元方向のX線強度変動が補正される。この補正方法では、各撮影間の二次元方向のX線強度の変動を容易に補正することができる。次いで、モアレ画像の解析が行われステップS15、再構成画像の作成に使用できるか否かが判断されるステップS16。理想的な送り精度によりマルチスリット12を一定の送り量で移動できた場合、図16に示すように、5ステップの撮影でマルチスリット12のスリット周期1周期分のモアレ画像5枚が得られる。各ステップのモアレ画像は0.2周期という一定周期間隔毎に縞走査をした結果であるので、各モアレ画像の任意の1画素に注目すると、その信号値を正規化して得られるX線相対強度は、図17に示すようにサインカーブを描く。よって、コントローラ5は得られた各ステップのモアレ画像のある画素に注目してX線相対強度を求める。各モアレ画像から求められたX線相対強度が、図17に示すようなサインカーブを形成すれば、一定周期間隔のモアレ画像が得られているので、再構成画像の作成に使用できると判断することができる。なお、上記サインカーブ形状は、マルチスリット12の開口幅、第1格子14及び第2格子15の周期、及び第1格子及び第2格子の格子間距離に依存し、また、放射光のようなコヒーレント光の場合には三角波形状となるが、マルチスリット効果によりX線が準コヒーレント光として作用する為、サインカーブを描くものとなる。ステップS15の解析は、被写体有りのモアレ画像と被写体無しのモアレ画像についてそれぞれ行われる。各ステップのモアレ画像の中にサインカーブを形成できないモアレ画像がある場合、再構成画像の作成に使用できないと判断されステップS16;NO、撮影のタイミングを変更して再撮影するよう指示する制御情報がコントローラ5からX線撮影装置1に送信されるステップS17。例えば、図17に示すように、3ステップ目は本来0.4周期のところ、周期がずれて0.35周期のモアレ画像が得られた場合であれば、駆動部122の送り精度の低下が原因例えば、パルスモータの駆動パルスへのノイズ重畳等と考えられる。よって、0.05周期分だけ撮影のタイミングを早めて3ステップ目のみ再撮影を行うよう指示すればよい。或いは、5ステップ全てについて再撮影し、3ステップ目のみ0.05周期分の撮影時間を早めるように指示してもよい。5ステップ全てのモアレ画像が所定量ずつサインカーブからずれている場合、駆動部122の起動から停止までの駆動パルス数を増やすか、或いは減らすように指示してもよい。X線撮影装置1では、当該制御情報に従って撮影のタイミングが調整され、再撮影が実行される。一方、再構成画像の作成にモアレ画像を使用できると判断された場合ステップS16;YES、被写体有りと被写体なしのそれぞれの複数のモアレ画像を用いて、被写体有りの再構成画像と被写体無しの再構成画像の作成が行われるステップS18~ステップS20。具体的には、複数のモアレ画像を加算することにより吸収画像が作成されるステップS18。また、縞走査法の原理を用いて干渉縞の位相が計算され、微分位相画像が作成されるステップS19。また、縞走査法の原理を用いて干渉縞のVisibilityが計算されVisibility=2×振幅÷平均値、小角散乱画像が作成されるステップS20。次いで、被写体無しの再構成画像を用いて、被写体有りの再構成画像から、干渉縞の位相の除去と、画像ムラアーチファクトを除去するための補正処理が行われるステップS21。これにより、被写体再構成画像が作成される。ステップS21の処理には、撮影時のマルチスリット12や第1格子14及び第2格子15のスリット方向変更に起因するX線の線量分布のムラ、当該スリットの製造バラつき起因の線量分布のムラ、及び、主に被写体ホルダー130の画像への写り込みによるムラ、を含む画像ムラアーチファクトを除去するための処理が含まれる。例えば、被写体有りの再構成画像が微分位相画像である場合には、被写体有りの微分位相画像の各画素の信号値から被写体無しの微分位相画像の対応する同じ位置の画素の信号値を減算する処理が行われる公知文献;TimmWeitkamp,AnaDiazand,ChristianDavidfranzPfeifferandMarcoStampanoniPeterCloetensandEricZieglerXrayPhaseImagingwithagratinginterferometerOPTICSEXPRESS,Vol13No1662966004、公知文献B;AtsushiMomoseWataruYashiroYoshihiroTakedaYoshioSuzukiandTadashiHattoriPhaseTomographybyXrayTalbotInterferometryforBiologicalImagingJapaneseJournalofAppliedPhysicsVol45No6A2006pp525452622006参照。被写体有りの再構成画像が吸収画像、小角散乱画像である場合には、公知文献Cに記載されているように、被写体有りの再構成画像の各画素の信号値を被写体無しの再構成画像の対応する画素の信号値で除算する割り算処理が行われる公知文献CFPfeifferMBech,OBunkPKraftEFEikenberryCHBroennimann,CGrunzweigandC.David,HardXraydarkfieldimagingusingagratinginterferometernaturematerialsVol7,1341372008。ただし、微分位相画像については、被写体有りのモアレ画像から各画素の微分位相信号値を計算する直前の中間値の段階で、被写体無しのモアレ画像から計算された中間値を除算し、その値のarctanで被写体の微分位相信号を得ることも可能であり第2の手法と呼ぶ、被写体有りのモアレ画像の撮影時と被写体無しのモアレ画像の撮影時とで撮影条件が完全に一致している場合を除いてはこの第2の手法を用いることが望ましい。図14Bに、第2の手法を用いた縞走査法による再構成画像作成処理Bのフローチャートを示す。縞走査法による再構成画像作成処理Bにおいては、まず、上述のステップS11~S16の処理が行われる。ステップS16において、被写体有りのモアレ画像及び被写体無しのモアレ画像が使用可能であると判断されると、被写体有りのモアレ画像と被写体無しのモアレ画像に基づいて、診断用の被写体再構成画像吸収画像、微分位相画像、小角散乱画像が作成されるステップS18B~ステップS20B。具体的に、ステップS18Bにおいては、被写体有りの複数のモアレ画像、被写体無しの複数のモアレ画像のそれぞれについて干渉縞を加算することにより被写体有りの吸収画像、被写体無しの吸収画像が作成される。そして、被写体無しの吸収画像を用いて被写体有りの吸収画像の画像ムラの除去が行われ上記公知文献c参照、診断用の被写体再構成画像としての吸収画像が作成される。ステップS19Bにおいては、縞走査法の原理を用いて、被写体有りのモアレ画像及び被写体無しのモアレ画像から位相が計算されることにより、診断用の被写体再構成画像としての微分位相画像が作成される。具体的には、上述のように、被写体有りのモアレ画像の各画素の微分位相信号値を計算する直前の中間値の段階で、被写体無しのモアレ画像の対応する画素から計算された中間値を除算し、その値のarctanで被写体有りの微分位相信号値を得る。この第2の手法によれば、被写体有りの撮影時と被写体無しの撮影時とで撮影条件が変化して被写体有りのモアレ画像と被写体無しのモアレ画像の干渉縞が変化してしまったような場合であっても、再構成画像を±πの範囲で得ることができる。ステップS20Bにおいては、縞走査法の原理を用いて、被写体有りのモアレ画像、被写体無しのモアレ画像のそれぞれに基づいて干渉縞のVisibilityが計算されることにより被写体有りの小角散乱画像及び被写体無しの小角散乱画像が作成される。そして、被写体無しの小角散乱画像を用いて被写体有りの小角散乱画像の画像ムラの除去が行われ公知文献c参照、診断用の被写体再構成画像としての小角散乱画像が作成される。図14A及び図14Bに示す縞走査法による再構成画像作成処理では、マルチスリット12、第1格子14及び第2格子15の各格子のスリット方向変更や被写体台特性に起因するX線の線量分布のムラだけでなく、撮影に用いられるX線検出器16の個々の画素の特性にバラツキがあっても、この影響を除去することができるので好ましい。従い、スリット方向を被写体に応じて可変としても、被写体に対するX線検出器16の配置方向を固定位置変更せずとすることができ、コントローラ5に表示される被写体再構成画像における被写体の表示向きは、コントローラ表示画面上で常に同一方向となるので、経過観察等で過去画像との比較読影を行う場合に、コントローラ5において被写体再構成画像の向きを揃える操作を行う必要がなくなるので、より好ましい。図18は、制御部51により実行されるフーリエ変換法による再構成画像作成処理を示すフローチャートである。フーリエ変換法による再構成画像作成処理は、制御部51と記憶部55に記憶されているプログラムとの協働により実行される。まず、ステップS31~S33においては、被写体有りの複数のモアレ画像及び被写体なしの複数のモアレ画像のそれぞれについて、X線検出器16の各画素のバラツキを補正するための補正処理が実行される。具体的には、オフセット補正処理ステップS31、ゲイン補正処理ステップS32、欠陥画素補正処理ステップS33が実行される。各処理の内容は、図14AのステップS11~S13で説明したものと同様であるので説明を援用する。次いで、被写体有りのモアレ画像と被写体無しのモアレ画像間のX線強度変動補正トレンド補正が行われるステップS34。X線強度変動補正の具体的な処理内容は、図14AのステップS14で説明したものと同様であるので説明を援用する。次いで、ステップS35以降の処理において、フーリエ変換法による被写体再構成画像の作成が行われる。フーリエ変換法による再構成画像の作成については、公知の手法により行うことができるMTakedaH.InaandS.KobayashiFourierTransformMethodeofFringePatternAnalysisforComputerBasedTopographyandInterferometryJOptSocAm721561982参照。まず、補正後の被写体有りのモアレ画像と被写体無しのモアレ画像のそれぞれがフーリエ変換二次元フーリエ変換されるステップS35。図19aに、第2の撮影モードで撮影された被写体有りのモアレ画像の一例を示す。図19aにおけるH1はマジックペンであり、H2はUSBメモリである。図19bに、図19aのモアレ画像を二次元フーリエ変換した結果を示す。図20aに、第2の撮影モードで撮影された被写体無しのモアレ画像の一例を示す。図20bに、図20aのモアレ画像を二次元フーリエ変換した結果を示す。フーリエ変換後の計算結果は複素数であるため、図19b、図20bでは実部と虚部のnorm振幅を表示している。図19b、図20bに示すように、1枚のモアレ画像をフーリエ変換すると、低周波成分0次成分と呼ぶと干渉縞周波数付近の成分1次成分と呼ぶ、又は、0次成分と1次成分に加えさらに高周波成分X線撮影装置1の干渉性に依存が並んで得られる。0次成分と1次成分の並ぶ方向はモアレ画像の縞の方向に関係しており、モアレ画像の縞の方向に対して略直角になる。ここで、マルチスリット12、第1格子14、第2格子15の格子の向きスリット方向と干渉縞及び0次成分、1次成分の並ぶ向きとの関係について説明する。例えば、図21のA1に示すように、マルチスリット12、第1格子14、第2格子15の格子の向きが縦であった場合、縞走査法用のモアレ画像第1格子14に対して第2格子をわずかに傾けて得られた画像の縞は、図21のA2に示すように横となる。フーリエ変換用のモアレ画像第2格子を更に傾けて得られた画像の縞は、図21のA3に示すように、図21のA2に比べて細かな横縞となる。フーリエ変換用のモアレ画像をフーリエ変換した画像は、図21のA4に示すように、0次成分、1次成分が縦に並んだ画像となる。図22のB1に示すように、マルチスリット12、第1格子14、第2格子15の格子の向きが横であった場合、縞走査法用のモアレ画像第1格子14に対して第2格子15をわずかに傾けて得られた画像の縞は、図22のB2に示すように縦となる。フーリエ変換用のモアレ画像第2格子15を更に傾けて得られた画像の縞は、図22のB3に示すように、図22のB2に比べて細かな縦縞となる。フーリエ変換用のモアレ画像をフーリエ変換した画像は、図22のB4に示すように、0次成分、1次成分が横に並んだ画像となる。図23のC1に示すように、マルチスリット12、第1格子14、第2格子15の向きが斜め45°であった場合、縞走査法用のモアレ画像第1格子14に対してわずかに第2格子を少し傾けて得られた画像の縞は、図23のC2に示すように斜め45°スリット方向とは逆方向の斜めとなる。フーリエ変換用のモアレ画像第2格子を更に傾けて得られた画像の縞は、図23のC3に示すように、図23のC2と同方向の、より細かな斜め縞となる。フーリエ変換用のモアレ画像をフーリエ変換した画像は、図23のC4に示すように、0次成分、1次成分が縞方向とは逆の斜め45°に並んだ画像となる。次いで、フーリエ変換により得られた画像被写体有り、被写体無しのそれぞれにおいて、0次成分が図24に示すHanning窓Wにより切り出されるステップS36。Hanning窓Wで切り出すことによりHanning窓Wの周辺部が0に落とされ、Hanning窓Wの中心部はそのまま通される。次いで、フーリエ変換により得られた画像において、1次成分が図25に示すようにキャリア周波数=モアレ周波数分シフトされ、Hanning窓Wで切り出されるステップS37。切り出しの窓関数はHanning窓に限定されず、用途に応じてHamming窓、ガウス窓等を使用しても良い。次いで、切り出された0次成分、1次成分のそれぞれが逆フーリエ変換されるステップS38。逆フーリエ変換が終了すると、逆フーリエ変換された0次成分、1次成分を用いて被写体有りと被写体無しのそれぞれの再構成画像の作成が行われるステップS39~ステップS41。具体的には、0次成分の振幅から吸収画像が作成されるステップS39。また、1次成分の位相から微分位相画像が作成されるステップS40。また、0次成分と1次成分の振幅の比=Visibilityから小角散乱画像が作成されるステップS41。次いで、被写体無しの再構成画像を用いて被写体有りの再構成画像から干渉縞の位相の除去と、画像ムラアーチファクトを除去するための補正処理が行われるステップS42。これにより、診断用の被写体再構成画像が作成される。ステップS42の処理は、図14AのステップS21で説明したものと同様であるので説明を援用する。画像ムラの補正が終了すると、フーリエ変換法による再構成画像作成処理は終了する。なお、上述の従来のフーリエ変換法では、0次成分、1次成分を切り出す際に、縦、横ともに高周波成分を捨てるため、空間分解能が落ちて全体的にぼけた画像となる。ここで、本願発明者らは、X線格子に一次元格子を用いる場合は、タルボ干渉計及びタルボロー干渉計の微分位相画像及び小角散乱画像の情報が格子マルチスリット12、第1格子14、第2格子15のスリット方向と直交する一方向のみであることに着目し、ステップS36、S37で用いる窓wを従来型の正方形ではなく、図26に示すように格子のスリット方向と直交する方向に延びる長方形とすれば、画像情報を含む格子のスリット方向と直交する方向の信号の高周波成分を落とさないように取り出すことができ、格子のスリット方向と直交する方向のぼけを低減することができることを見出した改良型フーリエ変換法と呼ぶ。本発明は、原理的に空間分解能の低下が避けられない格子のスリット方向と平行な方向は元々画像情報を含まないことを利用していることが特徴であり、後述する二回撮影時に二次元格子を用いたフーリエ変換法の撮影に対して大きなアドバンテージが得られる。図26は、図21のA4に長方形の窓Wを設定した例を示したものである。図27aに、縞走査法により得られた被写体の再構成画像の一例を示す。図27bに、改良型フーリエ変換法により得られた再構成画像の一例を示す。図27cに、従来のフーリエ変換法により得られた再構成画像の一例を示す。図27a~図27cの再構成画像は、格子のスリット方向を縦にして撮影して得られた微分位相画像である。図27aに示すように、縞走査法により得られた画像は、縦方向、横方向ともぼけが小さい。図27bに示すように、改良型フーリエ変換法により得られた画像は、縦方向のみぼけており、横方向にはぼけていない。図27cに示すように、従来のフーリエ変換法により得られた画像では、縦方向、横方向ともぼけている。なお、図27a~図27cにおいては、微分位相画像を示しているが、吸収画像、小角散乱画像においても各手法によるぼけの方向は同様である。このように、改良型フーリエ変換法では格子のスリット方向と平行な方向の信号成分のみがぼけるので、被写体の長手方向が格子のスリット方向と直交する方向となるように配置して1回目の撮影を行った後、被写体と格子との相対角度を90°回転させて2回目の撮影を行い、1回目と2回目の撮影で得られたモアレ画像を合成すれば、合成画像の周辺部を除く領域では、被写体の縦方向、横方向ともぼけの少ない2次元的な画像を取得することができる微分位相画像、小角散乱画像の場合。2次元格子を用いてフーリエ変換方式で撮影することも可能だが、縦横ともに1次成分が存在するため、切り出す窓wは縦横ともに狭い範囲に限定される。そのため空間分解能が大きく低下することを避けられない。一方、本方式によれば、1次元格子を使用して、分解能を大きく低下させずに2次元画像を生成することが可能となる。医用現場の放射線撮影に於いては、一般的に被写体の関心領域を撮影領域の中央部に位置させることが多いので、上記の周辺部の被写体情報欠損も問題なることが無い。また、撮影自体も2回で済むため、被写体の体動影響も抑制可能である。なお、撮影方向の変更被写体に対するスリット方向の変更する際には、縞走査法方式と同様に、マルチスリット12、第1格子14及び第2格子15を同時に90°回転させる必要がある。第1格子14と第2格子15の相対角度を90°変化させて1回目と2回目の撮影を行った場合、コントローラ5の制御部51は、1回目の撮影画像と2回目の撮影画像のそれぞれについて図18に示すフーリエ変換法による再構成画像作成処理を実行した後、2つの画像の合成を行う。なお、1回目の画像と2回目の画像とで被写体の同じ部分が同じ画素に描画されていない場合被写体の変形や移動があった場合は、何れか一方の画像を並行又は回転移動させて両画像の誤差が最も小さくなる位置に位置あわせを行った後、合成を行う。合成方法としては、様々な手法を用いることができる。例えば、1回目の撮影画像の画素をf1x、y、2回目の撮影画像の画素をf2x、y、合成画像の画素をgx、yとし、各画素で以下の計算をして二乗和の平方根をとってパワーの平均値をとる。gx、y=√f1x、y^2f2x、y^2また、例えば、1回目の撮影画像を赤、2回目の撮影画像を青にする等、カラー表示することとしてもよい。ところで、図27a~図27cの画像からもわかるように、縞走査法で得られる再構成画像のほうがフーリエ変換法で得られる再構成画像よりも鮮明でぼけが少ない。しかし、縞走査法用の撮影では連続して複数枚の撮影を行うため、撮影時間が検出器の取込時間やX線曝射前後の処理時間、機構動作時間等に応じて長くなり約1分、体動が生じやすい。これに対してフーリエ変換法では1回の撮影で1枚の撮影となるので、撮影時間はX線の曝射時間だけに依存し、5秒程度に抑えることができるため、体動抑制効果を期待できる。更に、改良型フーリエ変換法では、空間分解能の劣化小さく抑えられる。よって、例えば、1被写体の固定が可能なものは縞走査法、体動を抑制したい場合はフーリエ変換法とする、2簡易検査ではフーリエ変換法、より精密な検査となった場合には縞走査法とする、等、両者を併用することにより、目的に応じた画像の取得、及び、より患者の負担や再撮影の少ない撮影を行うことが可能となる。本実施形態においては、第1格子14及び第2格子15の相対角度を簡単に調整して縞走査法用の撮影とフーリエ変換法用の撮影を切り替えることができるので、撮影の目的に応じた最適な撮影を行うことが可能となる。再構成画像の作成が終了すると、コントローラ5の制御部51においては、図29、図31、図33に示す再構成画像補正処理の何れかが実行される。上述のように、被写体再構成画像は、被写体有りのモアレ画像と被写体無しのモアレ画像に基づいて作成される。ここで、被写体有りのモアレ画像と被写体無しのモアレ画像を同時に撮影することはできないため、それぞれ別個に撮影を行うが、これらの撮影間で撮影条件が変化した場合、被写体再構成画像にアーチファクトが出てしまい、被写体再構成画像のコントラストが悪く見えてしまうことがある。そこで、再構成画像補正処理を行うことで、被写体有りのモアレ画像の撮影時と被写体無しのモアレ画像の撮影時の撮影条件の変化に起因する被写体再構成画像におけるアーチファクトを補正する。微分位相値は、モアレ画像の位相シフトから計算される公知文献D参照。この位相シフトは±πラジアンの範囲で表わされ、微分位相値は位相シフト値に装置固有の定数を乗算したものになる。つまり、微分位相値は、定数×±πラジアンの範囲で表わされるが、ここでは簡単のため、微分位相値が±πで表わされていると仮定して記述する公知文献D:WataruYashiroYoshihiroTakedaandAtsushiMomoseEfficiencyofcapturingaphaseimageusingconebeamxrayTalbotinterferometryJournalofOpticalSocietyofAmericaAVol25No8pp20252039。撮影条件の変化としては、例えば、熱等の環境変化や被写体の接触等による格子マルチスリット12、第1格子14、第2格子15の配置の変化等が挙げられる。本願発明者等の検討によれば、格子配置にずれが生じた場合、被写体再構成画像のうち微分位相画像について、X線検出器面上の画素配列x、yに対して、下記の式1で示す一次関数又は式2で示す二次関数で表される位相誤差が各画素の信号値に加わることがわかった。そして、微分位相画像の各画素の信号値から下記式1又は式2により算出された補正値Icompを減算することにより、格子配置の変化に起因するアーチファクトが補正できることを見出した。一次関数:Icomp=axbyc式1二次関数:Icomp=ax^2bxycy^2dxeyf式2例えば、格子面の方向図28ではxy平面を維持したまま格子配置がz方向にずれると、主に、格子方向スリット方向と垂直方向図28ではx方向のアーチファクトが発生する。一方、z方向の格子配置を維持したまま格子面図28ではxy平面の角度がz軸周り図28ではθ方向に回転すると、主に、格子方向スリット方向と平行方向図28ではy方向のアーチファクトが発生する。これらの格子配置方向の変化によるアーチファクトは、式1の一次関数、式2の二次関数のいずれにおいても補正をすることができる。二次関数による補正では、z方向のずれ、z軸周りの回転によるアーチファクトが補正できることは勿論であるが、更に、x軸周りy軸周りの回転による誤差やこれらの組合せについても補正することができる。一次関数による補正まず、一次関数による補正について説明する。図29は、一次関数による再構成画像補正処理を示すフローチャートである。この処理は、制御部51と記憶部55に記憶されているプログラムとの協働により実行される。まず、微分位相画像において、被写体が描画されていない領域から少なくとも3点の測定点又は領域でもよいが設定され、その位置座標x1、y1~x3、y3、及び信号値I1~I3が取得されるステップS51。この3点は、3画素であってもよいし、3領域からそれぞれの近傍で平均値をとるのでもよい。信号値として平均値の場合は、その中心画素の座標を位置座標として用いる。図30に示すD1は、補正前の微分位相画像における測定点の設定例を示す図である。測定点領域は、例えば、被写体が描画されていない領域を画像解析により判別してその領域内から自動的に設定してもよいし、表示部53に再構成画像を表示して操作部52によりユーザが被写体領域を矩形で指定する等、手動で設定することとしてもよい。画像解析により被写体が描画されていない領域を判別して測定点領域を設定する場合は、吸収画像及び小角散乱画像の平均が1になる領域を被写体が描画されていない領域として判別することができる。また、例えば、X線検出器16の撮影可能範囲内X線照射可能範囲内に被写体を載置しない領域を予め設定しておくこととしてもよい。測定点領域は、格子配置によるx方向及びy方向の微分位相変化が測定できるように、図30のD1に示すようにx、y方向のそれぞれに直線上に離れた点を含むことが好ましい。なお、微分位相画像内に不連続点図32のD11参照が生じている場合は、基本的には不連続点を含まない領域内で測定点領域を設定する。しかし、不連続点をまたいで存在する被写体が描画されていない領域でも、画像上に基準位置を設け、その点と上記被写体が描出されていない領域を結ぶ曲線上の不連続点で微分位相がπから-πへ変化する場合には設定領域の画素の信号値に2π、不連続点で微分位相が-πからπへ変化する場合は設定領域の画素の信号値に-2πを加算したものを補正関数導出時の信号値とすることにより測定点領域として用いることができる。ここで、微分位相画像には誤差である成分Icompが含まれているため、微分位相画像内の画素の信号値は本来-πIcomp≦信号値≦πIcompに分布するはずであるが、信号値のとり得る範囲は-π≦信号値≦πに限られる。そのため、絶対値がπを超える信号値は2nπ間隔で折り返されて表現され、位相が不連続となる。このように、±πで折り返されて表現され、位相が不連続となっている画素を不連続点と呼ぶ。次いで、測定点の位置座標及び信号値に基づいて、下記の式3により一次関数の係数a、b、cが算出されることにより、一次の補正関数Icomp=axbycが算出されるステップS52。次いで、微分位相画像の各画素の信号値Iについて、補正値Icomp=axbycが計算され、各画素の信号値IからIcompが減算される。ただし、2n1π<Icomp<2n1πn≠0である場合、式4により補正値Icomp2が算出され、各画素の信号値Iから補正値Icomp2が減算されるステップS53。Icomp2=Icomp-2nπ式4次いで、補正済みの微分位相画像の各画素の信号値Iについて、I>πの場合は2π、I≦-πの場合は-2πの減算が行われ、信号値が±πの範囲内に収まるように補正されステップS54、一次関数による再構成画像補正処理は終了する。図30のD2に、図30のD1に示す微分位相画像を上記処理により補正した補正結果を示す。図30のD1に示す画像ではy軸方向の上部と下部とで明るさが不均一な画像となっていたが、図30のD2に示す画像では明るさが均一となるように補正されている。その結果、全体的にコントラストがはっきりした画像となっている。なお、図29に示すフローチャートでは、3点の測定点領域をとる場合を例にとり説明したが、3点以上の測定点領域をとって最小二乗法により一次関数の係数a、b、cを算出することとしてもよい。二次関数による補正次に、二次関数による補正について説明する。図31は、二次関数による再構成画像補正処理を示すフローチャートである。この処理は、制御部51と記憶部55に記憶されているプログラムとの協働により実行される。まず、微分位相画像において、被写体が描画されていない領域から少なくとも6点領域の測定点領域が設定されるステップS61。測定点領域の設定については、一次関数による再構成画像補正処理において説明したものと同様であるので説明を援用する。なお、測定点領域は、図32のD11に示す矩形のように、被写体及び不連続点を含まない領域で設定する。次いで、測定点領域の位置座標及び信号値に基づいて、二次関数の係数a、b、c、d、e、fが算出されることにより、二次の補正関数Icomp=ax^2bxycy^2dxeyfが算出されるステップS62。二次関数の係数a、b、c、d、e、fは、6以上の測定点領域から、これらの測定点領域の値を最も再現する二次の補正関数の係数を最小二乗法を用いて算出することにより求めることができる。次いで、微分位相画像の各画素の信号値について、補正値Icomp=ax^2bxycy^2dxeyfが計算され、各画素の信号値IからIcompが減算される。ただし、2n1π<Icomp<2n1πn≠0である場合、上記式4により補正値Icomp2が算出され、各画素の信号値Iから補正値Icomp2が減算されるステップS63。図32のD12に、図32のD11に示す微分位相画像にステップS63の減算処理を施した結果を示す。図32のD12に示すように、Icompの減算処理を行った結果、画像全体の明るさが略均一に補正され、全体的にコントラストがはっきりした画像となっている。また、不連続点の大部分の信号値が-π≦信号値≦πの間に収まり、D1において認識できなかった被写体が認識できるようになっている。次いで、補正済みの微分位相画像の各画素の信号値Iについて、I>πの場合は2π、I≦-πの場合は-2πの減算が行われ、信号値Iが±πの範囲内に収まるように補正されステップS64、二次関数による再構成画像補正処理は終了する。図32のD13に、図32のD12に示す減算処理済みの微分位相画像に対して更にステップS64の信号値の範囲の補正処理を行った結果を示す。ステップS64の処理の結果、全ての画素の信号値が±πの範囲内に収まるため、格子配置のずれによるアーチファクトをほぼ取り除くことができる。定数による補正格子が格子配列方向x方向にずれる、即ち、格子の周期の原点がずれると、微分位相値信号値の本来0であるべきところが0でなくなるという事態が起こる。そうすると、図34のD21に示すように微分位相画像の信号値全体がシフトしてしまい、診断に全く使用できない画像となる。このx方向の格子配置のずれによる影響の補正は、図29、図31に示す一次関数及び二次関数を用いた補正処理で補正することができるが、x方向の格子配置のずれによる影響の補正であれば、測定領域の中心値又は平均値をゼロへ移動するよう信号値全体を定数で補正する処理が有効であることを本件発明者等は見出した。図33に、定数による再構成画像補正処理のフローチャートを示す。この処理は、制御部51と記憶部55に記憶されているプログラムとの協働により実行される。まず、図34のD21に示すように、位相のずれを測定するための測定領域が設定されるステップS71。測定領域は、被写体が描画されている領域であっても被写体が描画されていない領域であってもよい。また、予め定められた領域を自動的に測定領域として設定してもよいし、ユーザが操作部52により手動で設定することとしてもよい。次いで、測定領域の画素の信号値のヒストグラムが作成されるステップS72。図35に、ステップS72において作成されるヒストグラムの一例を示す。図35の波形W11が、図34のR1に示す領域のヒストグラムである。図35に示すように、図34のR3に代表される画像D21では、信号値である微分位相値が全体的に約0.9πほどずれている。次いで、微分位相画像の各画素x、yの信号値から補正値Icompが減算されるステップS73。ここでは、Icompは以下の式5により求めることができる。Icomp=ヒストグラム中心値式5図34のD22に、図34のD21に示す微分位相画像にステップS73の減算処理を施した結果を示す。ステップS73の処理により、図35においてW11で示すヒストグラムをW12に示すヒストグラムにずらすことができる。その結果、図34のD22に示すように、画像の明るさを全体的にシフトさせることができる。次いで、補正済みの微分位相画像の各画素の信号値Iについて、I>πの場合は2π、I≦πの場合は-2πの減算が行われ、信号値Iが±πの範囲内に収まるように補正されステップS74、定数による再構成画像補正処理は終了する。図34のD23に、図34のD22に示す減算処理済みの微分位相画像に対して更にステップS74の処理を行った結果を示す。ステップS74の処理の結果、全ての画素の信号値が±πの範囲内に収まるため、格子配置のx方向のずれによるアーチファクトをほぼ取り除くことができる。なお上記処理は、微分位相値が±πの範囲で表わされる場合について記述しているが、0から2πで表わされる場合についても同様に考えることが可能である。上述の一次関数による再構成画像補正処理、二次関数による再構成画像補正処理定数による再構成画像補正処理何れを行うかは、格子配置が変化する方向の頻度や処理時間等に応じて操作部52によりユーザが選択することができるようにすることが好ましい。以上説明したように、X線撮影システムによれば、X線撮影装置1において撮影された被写体有りのモアレ画像と被写体無しのモアレ画像に基づいてコントローラ5において作成された被写体再構成画像に対して、被写体有りのモアレ画像を撮影したときと被写体無しのモアレ画像を撮影したときの撮影条件の差異に起因するアーチファクトの補正を行う。従って、被写体有りの撮影と被写体無しの撮影との間で、撮影条件、例えば、格子配置が変化したことにより生じる画質の低下を抑制し、診断に良好な被写体再構成画像を提供することが可能となる。なお、上記実施形態は本発明の好適な一例であり、これに限定されない。例えば、上記第1及び第2の実施形態では、X線源11、マルチスリット12、被写体台13、第1格子14、第2格子15、X線検出器16をこの順に配置以下、第1の配置と呼ぶしたが、X線源11、マルチスリット12、第1格子14、被写体台13、第2格子15、X線検出器16の配置以下、第2の配置と呼ぶとしても、第1格子14及び第2格子15は固定のまま、マルチスリット12の移動により、再構成画像を得ることが可能である。第2の配置においては、被写体の厚み分だけ、被写体中心と第1格子14は離れることになり、上記の実施形態に比べ感度の点でやや劣ることになるが、一方で、被写体への被曝線量低減を考慮すると、当該配置の方が第1格子14でのX線吸収分だけX線を有効に活用していることになる。また、被写体位置での実効的な空間分解能は、X線の焦点径、検出器の空間分解能、被写体の拡大率、被写体の厚さ等に依存するが、上記実施形態に於ける検出器の空間分解能が120μmガウスの半値幅以下の場合には、第1の配置よりも第2の配置の方が実効的な空間分解能は小さくなる。感度、空間分解能、及び、第1格子14でのX線吸収量等を考慮して、第1格子14、被写体台13の配置順をきめることが好ましい。また、上記実施形態においては、X線撮影装置1を、マルチスリット12を有し、マルチスリット12を第1格子14及び第2格子15に対して相対的に移動させることで縞走査法用の複数のモアレ画像を生成するタルボロー干渉計の構成とした例を説明したが、第1格子14と第2格子15とを一定周期間隔で相対移動させ、一定周期間隔での移動毎にX線源により照射されたX線に応じて放射線検出器が画像信号を読み取る処理を繰り返すことで縞走査法用の複数のモアレ画像を生成するタルボ干渉計としてもよい。そして、タルボ干渉計で生成された複数のモアレ画像を再構成することにより得られる微分位相画像を補正することとしてもよい。また、上記実施の形態においては、X線撮影装置1を、縞走査方式とフーリエ変換方式改良型含むの両方式が可能な装置として説明したが、これに限定されるものでは無く、縞走査方式専用の撮影装置としてもよいし、フーリエ変換方式改良型含む専用の装置を用いても良い。また、X線撮影装置1は、例えば、国際公開WO2010/050483号に記載のように、第1格子14及び第2格子15を2次元格子としたフーリエ変換方式専用の撮影装置としてもよいし、更に焦点位置近傍にマルチ格子2次元格子を併用したフーリエ変換方式専用の撮影装置を用いる構成としてもよい。このように、撮影条件に起因する被写体再構成画像の補正は、撮影装置がタルボ干渉計、タルボロー干渉計、フーリエ変換法式専用の撮影装置の何れにおいても適用することができる。また、撮影装置は一次元格子、二次元格子の何れを用いるものにおいても適用することができる。また、被写体有りの撮影と、被写体無しの撮影の順序は、上記実施形態に限定されず、何れを先としてもよい。被写体有りの再構成画像の作成と、被写体無しの再構成画像の作成の順序についても同様である。また、X線検出器16として、バッテリを内蔵し、無線により画像信号を本体部18に出力するケーブルレスのカセッテタイプFPDを用いてもよい。カセッテタイプFPDによれば、本体部18に接続するケーブル類を排除することができ、X線検出器16周辺の更なる小スペース化を図ることができる。小スペース化によって被写体の足下を広く構成し、より患者が接触し難い構成とすることができる。また、被写体台13は患者との接触により振動を伝えやすい。よって、被写体台13を高精度な位置関係が求められるマルチスリット、第1格子14、第2格子15等が含まれる保持部17と切り離し、別の保持部に保持することとしてもよい。図36は被写体台13を別の保持部13bにより保持したときの側面図、図37は平面図である。このように被写体台13を第1格子14、第2格子15等から離間させて別体構成とすることにより、マルチスリット12、第1格子14、第2格子15の位置関係に及ぶ影響をできるだけ減らし、当該位置関係の維持を図ることができる。被写体台13を別体構成とした場合、図36及び図37に示すように、被写体台13をz方向に移動させる駆動部13aを保持部13bに設ける。これにより、被写体の高さに合わせて、被写体台13の位置を調整することができる。被写体台13には患者の体重等の負荷がかかるが、被写体台13を保持部17と別体とすることにより、昇降する保持部17にかかる負荷を除去することができる。負荷に耐えるために保持部17を強化する必要がなく、コストを低減することができる。また、上記実施形態では、縞走査法用の撮影において、各ステップの撮影毎にマルチスリット12の移動と停止を繰り返す例を説明した。しかし、駆動部122の構成によっては、移動と停止を繰り返すことにより制御量と実際の移動量との誤差が累積拡大し、一定間隔毎のモアレ画像が得難いことが想定される場合には、連続的にマルチスリット12を移動させながら複数回の撮影を行う連続撮影方式が好ましい。曝射スイッチがONされると、マルチスリット12の移動を開始し、起動時の不安定移動領域を越え、安定移動領域に達した後、更に、マルチスリットを連続的に移動させて、所定量例えば4.56μm移動する毎にX線のパルス照射と画像信号の読み取りを繰り返す。連続撮影方式におけるX線源11にはパルス照射可能なX線管を用いることが好ましい。また、X線検出器16としては、対応できるフレームレート単位時間あたり撮影可能な回数が大きく、動画撮影が可能なFPDが好ましい。数百m秒~数秒の間に5回以上の撮影を行うことを想定すると、少なくとも10フレーム/秒のフレームレートが必要であり、好ましくは20フレーム/秒以上のフレームレートである。連続撮影方式の場合、各ステップの前後でさらに予備撮影を行うこととしてもよい。駆動部122が理想的な送り精度によりマルチスリット12を一定の送り量、つまり一定の移動速度で移動できた場合、図17に示すように各ステップのモアレ画像によりサインカーブを形成することができる。しかし、経年変化や駆動部122の起動時の慣性影響、グリスの粘性影響等によって送り量にずれが生じると、一定周期間隔のモアレ画像が得られない。例えば、図17に示すように、3ステップのモアレ画像は本来0.4周期に該当するが、3ステップのときの駆動部122の送り量がずれると、0.4周期前後のモアレ画像が得られる。このように各ステップのモアレ画像の周期がばらつくと、正確な位相が計算できず、再構成画像において被写体像を正確に再現できない。そこで、例えば各撮影時間±0.1秒の撮影時間で撮影を行う予備撮影を加えて合計15回の撮影を行う。これにより、各ステップにつきそれぞれ3枚のモアレ画像が得られるので、そのうちX線相対強度のサインカーブに最も近いモアレ画像を選択して用いる。これにより、駆動部122の送り量に誤差が生じたとしても、再構成画像の再現性の向上を図ることができる。予備撮影する調整時間として上記に挙げた±0.1秒は例示であり、調整時間はテスト撮影によって適宜決定すればよい。例えば、X線撮影装置1の設置時に、各ステップの撮影の前後で、±0.1秒、±0.2秒等、予備撮影時の調整時間を変えてテスト撮影を行い、最もサインカーブに一致しやすい調整時間を求めることとしてもよい。これにより、駆動部122の機器特性によって必要な調整時間が異なる場合にも対応することができる。また、被写体無しの一連5ステップの撮影のみを定期的に実施し、これらの各画像が前述するサインカーブにうまく合致するか否かを判断し、サインカーブからずれていると判断された場合には、装置メンテナンス必要性をコントローラ上で告知し、精密減速系等のメンテナンスを行わしめるものとすれば、高精細な診断用再構成画像を維持することが可能となる。また、X線撮影装置1の製品化にあたっては、被写体有りの撮影時と被写体無しの撮影時における格子配置の変化を抑制するため、下記の構成とすることが好ましい。例えば、図38に示すように、マルチスリット12を含むユニットマルチ格子ユニット120覆う第1のカバーユニット21、第1格子14を含むユニット第1格子ユニットU1及び第2格子15を含むユニット第2格子ユニットU2を覆う第2のカバーユニットが取り付けられた構成とすることが好ましい。図38は、第1格子ユニットU1、第2格子ユニットU2の周囲に設けられている第2のカバーユニット22を分解した状態を示す斜視図である。図38に示すように、第2のカバーユニット22は、第1格子ユニットU1、第2格子ユニットU2を前面側から覆うように配置される前面カバー部材221と、第1格子ユニットU1の上方を覆う上部カバー部材222と、第1格子ユニットU1、第2格子ユニットU2の両側部に配置される側面カバー部材223とを備えている。第2のカバーユニットを構成する前面カバー部材221、上部カバー部材222、側面カバー部材223は、例えば金属板をプレス加工等することにより形成される。上部カバー部材222は、中央部が切り欠かれた上面視コ字状の部材であり、被写体台13に載置される被写体、第1格子14及び第2格子15の上を上部カバー部材222の上面が覆わないようになっている。第2のカバーユニット22を取り付ける際は、側面カバー部材223を第1格子ユニットU1及び第2格子ユニットU2の両側部を覆うように基台部19をそれぞれ固定し、前面カバー部材221を第1格子ユニットU1及び第2格子ユニットU2を前面側から覆うように位置を合わせて側面カバー部材223の前面側端部にねじ止め固定する。さらに、上方から上部カバー部材222を被せて、基台部19、側面カバー部材223及び前面カバー部材221とそれぞれねじ止め固定する。なお、第2のカバーユニット22は、第1格子ユニットU1、第2格子ユニットU2に対する外部からの影響を遮断できるものであればよく、その材料や形状、構成、固定手法等はここに例示したものに限定されない。また、当該第2のカバーユニット22の少なくとも内面側には断熱部材が設けられていることが好ましい。また、被写体台13については、大きさは特に限定されないが、図1に示すように、撮影時に被写体となる患者の手指等をX線照射範囲内撮影可能領域内に置いた際に、患者の肘部分まで載置できる程度の長さ寸法を有し、腕レスト部を構成することが好ましい。手指から肘部分までを被写体台13に載置することにより、撮影対象である手指の位置及び姿勢を安定させることができ、撮影時の手ぶれ等を防ぐことができる。また、図39に示すように、被写体台13は、キャスタ131を備える脚部132を有し、支柱17や基台部19から独立していることが好ましい。脚部132は、支柱17と基台部19との間に配置されるようになっており、基台部19側の脚部132には、キャスタ131をロックするロック機構133が設けられていることが好ましい。なお、被写体台13の構成はここに例示したものに限定されない。例えば全ての脚部132にロック機構133を設けてもよいし、ロック機構133を設けず、支柱17又は基台部19の一端に被写体台13が固定されるようにしてもよい。なお、被写体台13は、支柱17又は基台部19に接触した際に衝撃を吸収することのできる衝撃吸収部材図示せずを備えていることが好ましい。その他、X線撮影システムを構成する各装置の細部構成及び細部動作に関しても、発明の趣旨を逸脱することのない範囲で適宜変更可能である。"