"放射線源から射出され、被検体を透過した放射線を検出して画像データを生成する放射線検出器と、前記放射線源と前記放射線検出器との間に配置された格子部と、前記画像データに基づき、アンラップエラーが生じやすいNG領域を検出するNG領域検出部と、前記画像データの前記NG領域にノイズ低減処理を施すノイズ低減処理部と、前記ノイズ低減処理後の前記画像データに基づき、所定の範囲にラップされた値で表現された位相微分画像を生成する位相微分画像生成部と、前記位相微分画像にアンラップ処理を施すアンラップ処理部と、を備えることを特徴とする放射線撮影装置。前記ノイズ低減処理は、前記画像データの前記NG領域内を平滑化する平滑化処理であることを特徴とする請求項1に記載の放射線撮影装置。前記平滑化処理は、前記NG領域内の各画素値を注目画素とし、注目画素の画素値を、注目画素を含む所定領域の画素値の平均値に置き換える移動平均処理であることを特徴とする請求項2記載の放射線撮影装置。オフセット画像として記憶するオフセット画像記憶部と、前記アンラップ処理後の位相微分画像から前記オフセット画像を減算するオフセット処理部とを備え、前記オフセット画像は、前記被検体を配置しない状態で行われるプレ撮影時において、前記位相微分画像生成部が、前記放射線検出器により生成される画像データに基づいて位相微分画像を生成し、前記アンラップ処理部がアンラップ処理を施したものであることを特徴とする請求項1から3いずれか1項に記載の放射線撮影装置。前記格子部は、放射線源からの放射線を通過させて第1の周期パターン像を生成する第1の格子と、前記第1の周期パターン像を部分的に遮蔽して第2の周期パターン像を生成する第2の格子と有し、前記放射線画像検出器は、前記第2の周期パターン像を検出して画像データを生成することを特徴とする請求項1から4いずれか1項に記載の放射線撮影装置。前記格子部は、前記第1の格子または第2の格子を所定の走査ピッチで移動させ、複数の走査位置に順に設定する走査機構を備え、前記放射線画像検出器は、前記各走査位置で前記第2の周期パターン像を検出して画像データを生成し、前記位相微分画像生成部は、前記放射線画像検出器により生成される複数の画像データに基づいて位相微分画像を生成することを特徴とする請求項5に記載の放射線撮影装置。前記走査機構は、前記第1の格子または第2の格子を、格子線に直交する方向に移動させることを特徴とする請求項6に記載の放射線撮影装置。前記走査機構は、前記第1の格子または第2の格子を、格子線に対して傾斜する方向に移動させることを特徴とする請求項6に記載の放射線撮影装置。前記位相微分画像生成部は、前記放射線検出器により得られる単一の画像データに基づいて前記位相微分画像を生成することを特徴とする請求項5に記載の放射線撮影装置。前記NG領域検出部は、画素値の強度変化を表す強度変調信号の平均強度、振幅、ビジビリティのうち1つまたは複数の組み合わせに基づいてNG領域を検出することを特徴とする請求項6から9いずれか1項に記載の放射線撮影装置。前記第1の格子は、吸収型格子であり、入射した放射線を幾何光学的に投影することにより前記第1の周期パターン像を生成することを特徴とする請求項5から10いずれか1項に記載の放射線撮影装置。前記第1の格子は、吸収型格子または位相型格子であり、入射した放射線にタルボ効果を生じさせて前記第1の周期パターン像を生成することを特徴とする請求項5から10いずれか1項に記載の放射線撮影装置。前記放射線源から放射された放射線を部分的に遮蔽して焦点を分散化するマルチスリットを備えることを特徴とする請求項1から12いずれか1項に記載の放射線撮影装置。放射線源と放射線検出器との間に格子部を配置して被検体の撮影を行うことにより得られた画像データに基づき、アンラップエラーが生じやすいNG領域を検出するNG領域検出ステップと、前記画像データの前記NG領域にノイズ低減処理を施すノイズ低減処理ステップと、前記ノイズ低減処理後の前記画像データに基づき、所定の範囲にラップされた値で表現された位相微分画像を生成する位相微分画像生成ステップと、前記位相微分画像にアンラップ処理を施すアンラップ処理ステップと、を備えることを特徴とする画像処理方法。しかしながら、被検体に骨部等のX線吸収能が高い高吸収体が含まれる場合には、該高吸収体がX線を大きく減衰させ、上記強度変調信号の強度や振幅を低下させるため、位相ズレ量の算出精度を低下させてしまう。これにより、位相微分画像の高吸収体領域では、アンラップエラーが生じやすくなる。このアンラップエラーには、本来不連続点でない箇所に不連続性が生じて不連続点と見なされることによりアンラップ処理が行われるケースと、本来不連続点である箇所の不連続性が低下して不連続点と見なされないことによりアンラップ処理が行われないケースとがある。例えば、アンラップ処理を行う経路上に骨部領域がある場合には、骨部領域ではアンラップエラーが生じやすく、アンラップエラーが生じた箇所以降の経路にエラー値上記関数の値域に相当する値が積算される。この結果、アンラップ処理後の位相微分画像にはアンラップ処理の経路方向に沿った筋状のノイズが生じる。このノイズが骨部領域の周囲の軟骨部領域軟部組織領域の一部に重なり、X線位相イメージングでの関心領域である肝心の軟部組織の画像化を阻害してしまうという問題がある。そこで、位相微分画像からアンラップエラーが生じやすい領域NG領域を判別し、NG領域を避けてそれ以外の領域OK領域のみにアンラップ処理経路を設定してアンラップ処理を行う方法が考えられる。しかしながら、この方法では、OK領域の形状に応じて、最適なアンラップ処理経路を逐一設定することが不可欠であり、複雑な経路設定処理が要されるという問題がある。本発明は、上記課題を鑑みてなされたものであり、アンラップ処理のための複雑な経路設定処理が不要であり、簡便にアンラップエラーを防止することを可能とする放射線撮影装置及び画像処理方法を提供することを目的とする。上記目的を達成するために、本発明の放射線撮影装置は、放射線源から射出され、被検体を透過した放射線を検出して画像データを生成する放射線検出器と、前記放射線源と前記放射線検出器との間に配置された格子部と、前記画像データに基づき、アンラップエラーが生じやすいNG領域を検出するNG領域検出部と、前記画像データの前記NG領域にノイズ低減処理を施すノイズ低減処理部と、前記ノイズ低減処理後の前記画像データに基づき、所定の範囲にラップされた値で表現された位相微分画像を生成する位相微分画像生成部と、前記位相微分画像にアンラップ処理を施すアンラップ処理部と、を備えることを特徴とする。前記ノイズ低減処理は、前記画像データの前記NG領域内を平滑化する平滑化処理であることが好ましい。また、前記平滑化処理は、前記NG領域内の各画素値を注目画素とし、注目画素の画素値を、注目画素を含む所定領域の画素値の平均値に置き換える移動平均処理であることが好ましい。オフセット画像として記憶するオフセット画像記憶部と、前記アンラップ処理後の位相微分画像から前記オフセット画像を減算するオフセット処理部とを備え、前記オフセット画像は、前記被検体を配置しない状態で行われるプレ撮影時において、前記位相微分画像生成部が、前記放射線検出器により生成される画像データに基づいて位相微分画像を生成し、前記アンラップ処理部がアンラップ処理を施したものであることが好ましい。前記格子部は、放射線源からの放射線を通過させて第1の周期パターン像を生成する第1の格子と、前記第1の周期パターン像を部分的に遮蔽して第2の周期パターン像を生成する第2の格子と有し、前記放射線画像検出器は、前記第2の周期パターン像を検出して画像データを生成することが好ましい。前記格子部は、前記第1の格子または第2の格子を所定の走査ピッチで移動させ、複数の走査位置に順に設定する走査機構を備え、前記放射線画像検出器は、前記各走査位置で前記第2の周期パターン像を検出して画像データを生成し、前記位相微分画像生成部は、前記放射線画像検出器により生成される複数の画像データに基づいて位相微分画像を生成することが好ましい。この場合、前記走査機構は、前記第1の格子または第2の格子を、格子線に直交する方向に移動させることが好ましい。また、前記走査機構は、前記第1の格子または第2の格子を、格子線に対して傾斜する方向に移動させることも好ましい。前記位相微分画像生成部は、前記放射線検出器により得られる単一の画像データに基づいて前記位相微分画像を生成することも好ましい。前記NG領域検出部は、画素値の強度変化を表す強度変調信号の平均強度、振幅、ビジビリティのうち1つまたは複数の組み合わせに基づいてNG領域を検出することが好ましい。前記第1の格子は、吸収型格子であり、入射した放射線を幾何光学的に投影することにより前記第1の周期パターン像を生成することが好ましい。また、前記第1の格子は、吸収型格子または位相型格子であり、入射した放射線にタルボ効果を生じさせて前記第1の周期パターン像を生成することも好ましい。前記放射線源から放射された放射線を部分的に遮蔽して焦点を分散化するマルチスリットを備えることが好ましい。本発明の画像処理方法は、放射線源と放射線検出器との間に格子部を配置して被検体の撮影を行うことにより得られた画像データに基づき、アンラップエラーが生じやすいNG領域を検出するNG領域検出ステップと、前記画像データの前記NG領域にノイズ低減処理を施すノイズ低減処理ステップと、前記ノイズ低減処理後の前記画像データに基づき、所定の範囲にラップされた値で表現された位相微分画像を生成する位相微分画像生成ステップと、前記位相微分画像にアンラップ処理を施すアンラップ処理ステップと、を備えることを特徴とする。本発明によれば、アンラップエラーが生じやすいNG領域を検出し、画像データのNG領域にノイズ低減処理を施したうえで、位相微分画像を生成してアンラップ処理を行うので、アンラップ処理のための複雑な経路設定処理が不要であり、どのような経路に沿ってアンラップ処理を行ったとしてもアンラップエラーは生じない。このように、本発明によれば、簡便にアンラップエラーを防止することができる。図1において、X線撮影装置10は、X線源11、格子部12、X線画像検出器13、メモリ14、画像処理部15、画像記録部16、撮影制御部17、コンソール18、及びシステム制御部19を備える。X線源11は、例えば、回転陽極型のX線管と、X線の照射野を制限するコリメータとを有し、撮影制御部17の制御に基づき、被検体Hに向けてX線を放射する。格子部12は、第1の格子21、第2の格子22、及び走査機構23を備える。第1及び第2の格子21,22は、X線照射方向であるz方向に関してX線源11に対向配置されている。X線源11と第1の格子21との間には、被検体Hが配置可能な間隔が設けられている。X線画像検出器13は、例えば、半導体回路を用いたフラットパネル検出器であり、第2の格子22の背後に、検出面13aがz方向に直交するように配置されている。第1の格子21は、z方向に直交する格子面内の一方向であるy方向に延伸された複数のX線吸収部21a及びX線透過部21bを備えた吸収型格子である。X線吸収部21a及びX線透過部21bは、z方向及びy方向に直交するx方向に交互に配列されており、縞状のパターンを形成している。第2の格子22は、第1の格子21と同様にy方向に延伸され、かつx方向に交互に配列された複数のX線吸収部22a及びX線透過部22bを備えた吸収型格子である。X線吸収部21a,22aは、金Au、白金Pt等のX線吸収性を有する材料により形成されている。X線透過部21b,22bは、シリコンSiや樹脂等のX線透過性を有する材料や空隙により形成されている。第1の格子21は、X線源11から放射されたX線を部分的に通過させて第1の周期パターン像以下、G1像というを生成する。第2の格子22は、第1の格子21により生成されたG1像を部分的に透過させて第2の周期パターン像以下、G2像というを生成する。被検体Hが配置されていない場合において、G1像は、第2の格子22の格子パターンとほぼ一致する。X線画像検出器13は、G2像を検出して画像データを生成する。メモリ14は、X線画像検出器13から読み出された画像データを一時的に記憶する。画像処理部15は、メモリ14に記憶された画像データに基づいて位相微分画像を生成し、この位相微分画像に基づいて位相コントラスト画像を生成する。画像記録部16は、位相微分画像と位相コントラスト画像とを記録する。走査機構23は、第2の格子22をx方向に並進移動させ、第1の格子21に対する第2の格子22の相対位置を順次に変更する。走査機構23は、圧電アクチュエータや静電アクチュエータにより構成され、後述する縞走査を実行するために、撮影制御部17の制御に基づいて駆動される。メモリ14には、縞走査の各走査位置でX線画像検出器13により得られる画像データが一括して記憶される。コンソール18は、操作部18a及びモニタ18bを備えている。操作部18aは、キーボードやマウス等により構成され、X線源11の管電圧、管電流、照射時間等の撮影条件の設定や、本撮影またはプレ撮影のモード選択、撮影実行指示等の操作入力を可能とする。本撮影とは、X線源11と第1の格子21との間に被検体Hを配置した状態で行う撮影モードである。プレ撮影とは、X線源11と第1の格子21との間に被検体Hを配置せずに行う撮影モードである。詳しくは後述するが、プレ撮影は、第1及び第2の格子21,22の製造誤差や配置誤差等により生じるバックグランド成分をオフセット画像として取得するために用いられる。モニタ18bは、撮影条件等の撮影情報や、画像記録部16に記録された位相微分画像及び位相コントラスト画像の表示を行う。システム制御部19は、操作部18aから入力される信号に応じて各部を統括的に制御する。図2において、X線画像検出器13は、入射X線により半導体膜図示せずに生じた電荷を収集する画素電極31と、画素電極31によって収集された電荷を読み出すためのTFTThinFilmTransistor32とを備えた画素部30が2次元状に多数配列されたものである。半導体膜は、例えば、アモルファスセレンにより形成されている。また、X線画像検出器13は、ゲート走査線33、走査回路34、信号線35、及び読み出し回路36を備える。ゲート走査線33は、画素部30の行ごとに設けられている。走査回路34は、TFT32をオン/オフするための走査信号を各ゲート走査線33に付与する。信号線35は、画素部30の列ごとに設けられている。読み出し回路36は、各信号線35を介して画素部30から電荷を読み出し、画像データに変換して出力する。各画素部30の詳細な層構成については、例えば、特開2002-26300号公報に記載されている層構成と同様である。読み出し回路36は、積分アンプ、A/D変換器、補正回路いずれも図示せず等を備える。積分アンプは、各画素部30から信号線35を介して出力された電荷を積分して画像信号を生成する。A/D変換器は、積分アンプにより生成された画像信号を、デジタル形式の画像データに変換する。補正回路は、画像データに対して、暗電流補正、ゲイン補正、リニアリティ補正等を行う。この補正後の画像データがメモリ14に記憶される。X線画像検出器13は、入射X線を半導体膜で直接電荷に変換する直接変換型に限られず、ヨウ化セシウムCsIやガドリウムオキシサルファイドGOS等のシンチレータで入射X線を可視光に変換し、可視光をフォトダイオードで電荷に変換する間接変換型であってもよい。さらに、X線画像検出器13を、シンチレータとCMOSセンサを組み合わせて構成してもよい。図3において、X線源11から照射されるX線は、X線焦点11aを発光点としたコーンビームである。第1の格子21は、タルボ効果が生じず、X線透過部21bを通過したX線を幾何光学的に投影するように構成される。具体的には、x方向へのX線透過部21bの幅を、X線源11から照射されるX線のピーク波長より十分大きな値とし、X線の大部分がX線透過部21bで回折しないようにすることで実現される。X線源11の回転陽極としてタングステンを用い、管電圧を50kVとした場合には、X線のピーク波長は約0.4Åである。この場合には、X線透過部21bの幅を1~10μm程度とすればよい。これにより、G1像は、第1の格子21からz方向下流への距離に依らず、常に第1の格子21の自己像となる。G1像は、X線焦点11aからz方向下流への距離に比例して拡大される。第2の格子22の格子ピッチp2は、前述のように、第2の格子22の格子パターンが第2の格子22の位置におけるG1像に一致するように設定されている。具体的には、第2の格子22の格子ピッチp2は、第1の格子21の格子ピッチp1、X線焦点11aと第1の格子21との間の距離L1、第1の格子21と第2の格子22との間の距離L2と、下式1をほぼ満たすように設定されている。G1像は、被検体HでX線に位相変化が生じて屈折することにより変調される。この変調量には、被検体HでのX線の屈折角φxが反映される。同図には、被検体HでのX線の位相変化を表す位相シフト分布Φxに応じて屈折するX線の経路が例示されている。符号X1は、被検体Hが存在しない場合にX線が直進する経路を示し、符号X2は、被検体Hにより屈折したX線の経路を示している。位相シフト分布Φxは、X線の波長をλ、被検体Hのxz面内での屈折率分布をnx,zとして、下式2で表される。屈折角φxは、位相シフト分布Φxと、下式3の関係にある。第2の格子22の位置において、X線は、屈折角φxに応じた量だけx方向に変位する。この変位量Δxは、X線の屈折角φxが微小であることに基づいて、近似的に下式4で表される。このように、変位量Δxは、位相シフト分布Φxの微分値に比例する。したがって、変位量Δxを後述する縞走査により検出することにより、位相シフト分布Φxの微分値が得られ、位相微分画像が生成される。縞走査は、格子ピッチp2をM個に分割した値p2/Mを走査ピッチとし、走査機構23により、この走査ピッチで第2の格子22を並進移動させ、第2の格子22を並進移動させるたびに、X線源11からX線を放射してG2像をX線画像検出器13により撮影することにより行われる。Mは3以上の整数であり、例えば、M=5であることが好ましい。上式1を僅かに満たさない場合や、第1の格子21と第2の格子22との間にz方向周りの回転や、xy平面に対する傾斜が僅かに生じている場合には、G2像にはモアレ縞が生じる。このモアレ縞は、第2の格子22の並進移動に伴って移動し、x方向への移動距離が格子ピッチp2に達すると元のモアレ縞に一致する。このモアレ縞の移動を確認することで、第2の格子22の並進移動量を検証することができる。上記縞走査により、X線画像検出器13の各画素部30について、M個の画素値が得られる。図4に示すように、M個の画素値Ikは、第2の格子22の走査位置kに対して周期的に変化する。走査位置kは、第2の格子22を一周期分並進移動させた場合の走査ピッチp2/Mごとの各位置である。走査位置kに対する画素値Ikの変化を表す信号を強度変調信号と呼ぶ。同図中の破線は、被検体Hを配置しない状態で得られる強度変調信号を示している。これに対して、実線は、被検体Hを配置した状態で、被検体Hにより位相ズレ量ψxが生じた強度変調信号を示している。この位相ズレ量ψxは、上記変位量Δxと下式5の関係にある。したがって、各画素部30について、縞走査で得られるM個の画素値Ikに基づき、強度変調信号の位相ズレ量ψxを求めることにより、位相微分画像が得られる。次に、位相ズレ量ψxの算出方法について説明する。強度変調信号は、一般に下式6で表される。ここで、A0は入射X線の平均強度を表し、Anは強度変調信号の振幅を表す。nは正の整数、iは虚数単位である。なお、図4に示すように、強度変調信号が正弦波を描く場合には、n=1である。本実施形態では、走査ピッチp2/Mが一定であるため、下式7が成立する。上式7を上式6に適用すると、位相ズレ量ψxは、下式8で表される。ここで、arg[...]は、複素数の偏角を抽出する関数である。また、位相ズレ量ψxは、逆正接関数を用いて下式9のように表すことも可能である。複素数の偏角は、値域が-πからπの範囲であるため、上式8に基づいて位相ズレ量ψxを算出した場合には、位相ズレ量ψxは、-πからπの範囲にラップされた値を取る。これに対して、逆正接関数は、通常、値域が-π/2からπ/2の範囲であるため、上式9に基づいて位相ズレ量ψxを算出した場合には、位相ズレ量ψxは、-π/2からπ/2の範囲にラップされた値を取る。なお、上式9において、逆正接関数内の分母及び分子の正負を判別することにより、値域を-πからπとすることができるため、-πからπの範囲で位相ズレ量ψxを算出することも可能である。本実施形態では、各画素部30について位相ズレ量ψxを算出することにより得られるデータを位相微分画像という。なお、位相ズレ量ψxに定数を乗じたり加算したりしたデータを位相微分画像としてもよい。以下、位相微分画像の画素値は、幅αを有する所定の値域例えば0からαの範囲にラップされているとする。図5において、画像処理部15は、NG領域検出部40、ノイズ低減処理部41、位相微分画像生成部42、アンラップ処理部43、オフセット画像記憶部44、オフセット処理部45、及び位相コントラスト画像生成部46を備える。画像処理部15には、本撮影によりメモリ14に格納されたM枚分の画像データ本画像データ50と、プレ撮影によりメモリ14に格納されたM枚分の画像データプレ画像データ51とがそれぞれ入力される。NG領域検出部40は、本画像データ50に基づき、位相微分画像を生成した場合にアンラップエラーが生じやすい領域以下、NG領域というを検出する。具体的には、NG領域検出部40は、本画像データ50に基づき、各画素部30について前述の強度変調信号を生成し、強度変調信号の平均強度A0が所定の閾値より低い領域、振幅A1が所定の閾値より低い領域、またはビジビリティA1/A0が所定の閾値より低い領域をNG領域として検出する。このNG領域は、被検体Hに含まれる高吸収体領域被検体Hが人体である場合には、X線吸収能が高い骨部等に相当する。これは、X線が高吸収体で吸収されることにより、平均強度A0、振幅A1、またはビジビリティA1/A0が低下することに基づいている。なお、平均強度A0、振幅A1、ビジビリティA1/A0のうち2以上を組み合わせてNG領域を検出してもよい。また、NG領域が散在して、ある程度の大きさを有する集合領域として得られない場合には、上記閾値を変化させてNG領域の大きさを調整すればよい。ノイズ低減処理部41は、M枚分の各本画像データ50に対して、NG領域検出部40により検出されたNG領域に移動平均処理を施すことによりノイズを低減する。具体的には、図6に示すように、ノイズ低減処理部41は、各本画像データ50のNG領域内の各画素を注目画素Qx,yとし、各注目画素Qx,yの画素値を、注目画素Qx,yの画素値と、注目画素Qx,yの周囲に隣接する8個の隣接画素Qx-1,y-1,Qx,y-1,Qx1,y-1,Qx-1,y,Qx1,y,Qx-1,y1,Qx,y1,Qx1,y1の画素値との平均値で置き換える移動平均処理を施す。これにより、各本画像データ50のNG領域内の画素値が平滑化され、ノイズが低減される。本実施形態では、注目画素Qx,yを含む9個の画素を用いて平均値を算出しているが、平均値を算出する画素数は任意である。但し、この画素数が大きい場合には、本画像データ50に基づいて生成される位相微分画像のNG領域の解像度が必要以上に低下してしまう。一方、この画素数が少ない場合には、ノイズが十分に低減されず、後述するアンラップ処理でアンラップエラーが発生する可能性がある。このため、移動平均処理において平均値を算出するための画素数は、解像度とアンラップエラーの発生度合いとの兼ね合いから適正に設定されることが好ましい。位相微分画像生成部42は、ノイズ低減処理部41によりNG領域にノイズ低減処理が施されたM枚分の各本画像データ50を用い、上式8または上式9に基づいて演算を行うことにより第1の位相微分画像を生成する。一方、プレ撮影時には、M枚分のプレ画像データ51は位相微分画像生成部42に直接入力され、位相微分画像生成部42は、上式8または上式9に基づいて演算を行うことにより第2の位相微分画像を生成する。アンラップ処理部43は、位相微分画像生成部42により生成された第1及び第2の位相微分画像にそれぞれアンラップ処理を施す。アンラップ処理部43は、プレ撮影時には、第2の位相微分画像にアンラップ処理を施したものを、オフセット画像としてオフセット画像記憶部44に記憶させる。なお、オフセット画像記憶部44は、新たにプレ撮影を行って新たなオフセット画像が入力された場合には、既に記憶しているオフセット画像を消去して、新たに入力されたオフセット画像を記憶する。また、アンラップ処理部43は、本撮影時には、第1の位相微分画像にアンラップ処理を施したものを、アンラップ済み位相微分画像としてオフセット処理部45に入力する。オフセット処理部45は、アンラップ処理部43から入力されたアンラップ済み位相微分画像から、オフセット画像記憶部44に記憶されたオフセット画像を減算するオフセット補正を行う。位相コントラスト画像生成部46は、オフセット処理部45によりオフセット補正が施されたアンラップ済み位相微分画像を、x方向に沿って積分処理することにより、位相シフト分布を表す位相コントラスト画像を生成する。オフセット補正後の位相微分画像と位相コントラスト画像とが画像記録部16に記録される。次に、図7を用いてアンラップ処理部43によるアンラップ処理方法を具体的に説明する。アンラップ処理部43は、第1の位相微分画像60について、端部の一辺に沿って起点群SP1~SPnを設定するとともに、起点群SP1~SPnの配列方向と直交する方向に直線経路UR1~URnを設定する。アンラップ処理部43は、起点群SP1~SPnをアンラップ処理を行った後、各直線経路UR1~URnに沿ってアンラップ処理を行う。このアンラップ処理は、アンラップ経路上において隣接する画素の画素値の変化量が所定量例えば、α/2以上である不連続点を検出し、検出した不連続点以降のアンラップ経路上に存在する画素の画素値に、前述の画素値がラップされた幅αを加算または減算する処理である。これにより、アンラップ経路上の画素値の変化がほぼ連続化される。第1の位相微分画像60のNG領域は、本画像データ50の段階でノイズ低減処理部41が低減されているため、画素値のノイズが少なく、アンラップエラーの発生が防止される。アンラップ処理部43は、第2の位相微分画像についても上記と同一のアンラップ処理を行う。なお、本実施形態では、起点群SP1~SPnの配列方向をx方向としているが、y方向としてもよい。以下、図8及び図9に示すフローチャートを参照しながらX線撮影装置10の作用を説明する。操作部18aを用いて撮影モードの選択がなされるとステップS10、選択された撮影モードがプレ撮影であるか否かの判定が行われるステップS11。プレ撮影である場合にはステップS11でYES、撮影指示の待受状態となるステップS12。操作部18aを用いて撮影指示がなされるとステップS12でYES、走査機構23により第2の格子22が所定の走査ピッチずつ並進移動されながら、各走査位置kにおいて、X線源11によるX線照射及びX線画像検出器13によるG2像の検出が行われるステップS13。この縞走査の結果、M枚分のプレ画像データ51が生成され、メモリ14に格納される。M枚分のプレ画像データ51は、画像処理部15に読み出される。画像処理部15内では、位相微分画像生成部42により、M枚分のプレ画像データ51に基づき第2の位相微分画像が生成されるステップS14。この第2の位相微分画像は、アンラップ処理部43によりアンラップ処理が施されステップS15、アンラップ処理後の第2の位相微分画像がオフセット画像としてオフセット画像記憶部44に記憶されるステップS16。プレ撮影動作は、以上で終了する。なお、このプレ撮影は、X線撮影装置10の立ち上げ時等に被検体Hを配置しない状態で少なくとも一度行われればよく、本撮影の前に毎回行われる必要はない。次に、被検体Hが配置され、ステップS10の撮影モードの選択により本撮影が選択された場合にはステップS11でNO、撮影指示の待受状態となるステップS20。操作部18aを用いて撮影指示がなされるとステップS20でYES、縞走査が行われステップS21、メモリ14にM枚分の本画像データ50が格納される。M枚分の本画像データ50は、画像処理部15に読み出される。画像処理部15内では、NG領域検出部40によりNG領域の検出が行われステップS22、ノイズ低減処理部41により各本画像データ50のNG領域に移動平均処理が施され、ノイズが低減されるステップS23。この後、位相微分画像生成部42により、ノイズ低減処理が施されたM枚分の本画像データ50に基づき第1の位相微分画像が生成されるステップS24。この第1の位相微分画像は、アンラップ処理部43によりアンラップ処理が施されステップS25、アンラップ処理後の第1の位相微分画像がアンラップ済み位相微分画像としてオフセット処理部45に入力される。オフセット処理部45では、アンラップ済み位相微分画像から、オフセット画像記憶部44に記憶されたオフセット画像を減算するオフセット補正が行われるステップS26。そして、位相コントラスト画像生成部46がオフセット補正後の位相微分画像を積分処理することにより、位相コントラスト画像が生成されステップS27、オフセット補正後の位相微分画像及び位相コントラスト画像が画像記録部16に記録された後、モニタ18bに画像表示されるステップS28。以上のように、本撮影時に本画像データ50からアンラップエラーが生じやすいNG領域を検出し、このNG領域にノイズ低減処理を施したうえで第1の位相微分画像を生成しているため、第1の位相微分画像のアンラップ処理時にはアンラップエラーの発生が防止される。このため、NG領域骨部領域等の周囲に位置する軟部組織領域図7参照は、従来のようにアンラップエラーによる筋状のノイズで画像化が阻害されることはなく、視認性が良好である。なお、上記実施形態では、図7に示すように、アンラップ処理部43は、NG領域を含む第1及び第2の位相微分画像の全体を対象としてアンラップ処理を行っているが、NG領域以外のOK領域に起点を設定して、OK領域のみを対象としてアンラップ処理を行ってもよい。また、上記実施形態では、ノイズ低減処理部41は、線形フィルタ処理である移動平均処理によってNG領域内のノイズの低減を行っているが、NG領域内の画素値を平滑化しノイズを低減することができる処理であれば、その他の平滑化処理方法を用いることが可能である。例えば、ノイズ低減処理として、周知のメディアンフィルタ処理や最頻値フィルタ処理等の非線形フィルタ処理を用いてもよい。メディアンフィルタ処理は、注目画素を含む所定領域内の画素値の中央値で注目画素の画素値を置き換える処理である。最頻値フィルタ処理は、注目画素を含む所定領域内の画素値の最頻値で注目画素の画素値を置き換える処理である。さらに、線形処理と非線形処理とを組み合わせてノイズ低減処理を行うことも可能である。例えば、フィルタ行列の重み付け線形処理と閾値判定非線形処理とによりノイズ低減処理を行う。また、上記実施形態では、ノイズ低減処理部41は、NG領域内で注目画素を1画素ずつ変更し、各注目画素の画素値を、それを含む領域の画素値の平均値、中央値、最頻値等で置き換える方法を用いているが、NG領域を複数の領域に区分し、区分された各領域内の画素値を、該領域内の画素値の平均値、中央値、最頻値等ですべて置き換える方法によりノイズの低減を行ってもよい。また、上記実施形態では、ノイズ低減処理部41は、強度変調信号の平均強度、振幅、またはビジビリティに基づいてNG領域を検出しているが、NG領域の検出方法はこれに限られず、強度変調信号の平均強度の画素間のばらつきすなわち、吸収画像の画素間のばらつきや、位相微分画像の画素間のばらつきが所定値より大きい領域をNG領域として検出してもよい。なお、この位相微分画像の画素間のばらつきは、第1及び第2の格子21,22の格子線に直交する方向x方向へのばらつきとすることが好ましい。また、位相微分画像の各画素について絶対値を取り、この絶対値が所定値を超える箇所を検出することにより、高吸収体領域のエッジ部分を検出することができるため、このエッジ部で囲われる領域をNG領域として検出してもよい。また、強度変調信号の平均強度や最大強度が所定値より大きく、強度変調信号に飽和が生じている領域をNG領域として検出してもよい。この強度変調信号の飽和は、被検体Hを透過せずに第1及び第2の格子21,22を介してX線画像検出器13に直接入射した画素領域素抜け領域で生じやすい。強度変調信号が飽和すると位相ズレ量ψxが正確に得られなくなるため、この素抜け領域もアンラップエラーが生じやすい領域である。以上の検出基準を適宜組み合わせてもよい。さらに、X線画像検出器13、第1の格子21、第2の格子22に欠陥が生じたり、ゴミなどが付着したりした場合には、特定の画素部30の画素値が常に高く、または低くなることがある。このような画素欠陥が生じた領域は、強度変調信号の平均強度、振幅、またはビジビリティが異常値を示すため、アンラップエラーが生じやすい領域となる。このような画素欠陥領域についても、上記の検出基準を適宜組み合わせることにより、NG領域として検出可能である。上記のような素抜けや画素欠陥に起因するNG領域が検出された場合においても、ノイズ低減処理部41によるノイズ低減処理によりアンラップエラーが防止される。このような素抜けや画素欠陥に起因するNG領域は、プレ画像データ51においても検出可能であるため、プレ画像データ51に対してもNG領域検出部40によるNG領域の検出、及びノイズ低減処理部41によるノイズ低減処理を施した後、位相微分画像生成部42により、第2の位相微分画像を生成することも好ましい。また、上記実施形態では、被検体HをX線源11と第1の格子21との間に配置しているが、被検体Hを第1の格子21と第2の格子22との間に配置してもよい。また、上記実施形態では、縞走査時に第2の格子22を格子線に直交する方向x方向に移動させているが、本出願人により特願2011-097090号として出願されているように、第2の格子22を格子線に対して傾斜する方向xy平面内でx方向及びy方向に直交しない方向に移動させてもよい。この移動方向は、xy平面内で、かつy方向以外であれば、いずれの方向であってもよい。この場合には、第2の格子22の移動のx方向成分に基づいて、走査位置kを設定すればよい。第2の格子22を格子線に対して傾斜する方向に移動させることにより、縞走査の一周期分の走査に要するストローク移動距離が長くなるため、移動精度が向上するといった利点がある。また、上記実施形態では、縞走査時に第2の格子22を移動させているが、第2の格子22に代えて、第1の格子21を格子線に直交する方向または傾斜する方向に移動させてもよい。また、上記第実施形態では、X線源11から射出されるコーンビーム状のX線を射出するX線源11を用いているが、平行ビーム状のX線を射出するX線源を用いることも可能である。この場合には、上式1に代えて、p2=p1をほぼ満たすように第1及び第2の格子21,22を構成すればよい。また、上記実施形態では、X線源11から射出されたX線を第1の格子21に入射させており、X線源11は単一焦点であるが、X線源11の射出側直後に、WO2006/131235号公報等に記されたマルチスリット線源格子を設けることにより、X焦点を分散化してもよい。これより、高出力のX線源を用いることが可能となり、X線量が向上するため、位相微分画像の画質が向上する。マルチスリットの格子線は、y方向に沿っている。この場合、マルチスリットのピッチp0は、下式10を満たす必要がある。ここで、距離L1は、マルチスリットから第1の格子21までの距離を表す。また、上記実施形態では、第1の格子21が入射X線を幾何光学的に投影するように構成しているが、WO2004/058070号公報等で知られているように、第1の格子21をタルボ効果が生じる構成としてもよい。第1の格子21でタルボ効果を生じさせるためには、X線の空間干渉性を高めるように、小焦点のX線光源を用いるか、上記マルチスリットを用いればよい。第1の格子21でタルボ効果を生じさせる場合には、第1の格子21を、吸収型格子に代えて、位相型格子とすることも好ましい。位相型格子は、吸収型格子のX線吸収部をX線位相形成部に置換することにより構成される。X線位相形成部は、隣接するX線透過部に対して所定の屈折率差を有する材料により形成される。この位相型格子としては、入射X線にπ/2の位相変調を与えて透過させるものと、入射X線にπの位相変調を与えて透過させるものがある。第1の格子21でタルボ効果が生じる場合には、第1の格子21の自己像G1像が、第1の格子21からz方向下流にタルボ距離Zmだけ離れた位置に生じるため、第1の格子21から第2の格子22までの距離L2をタルボ距離Zmとする必要がある。タルボ距離Zmは、第1の格子21の構成とX線のビーム形状とに依存する。第1の格子21が吸収型格子であり、X線源11から射出されるX線がコーンビーム状である場合には、タルボ距離Zmは、下式11で表される。ここで、mは正の整数である。また、第1の格子21がX線にπ/2の位相変調を与える位相型格子であり、X線源11から射出されるX線がコーンビーム状である場合には、タルボ距離Zmは、下式12で表される。ここで、mは0または正の整数である。また、第1の格子21がX線にπの位相変調を与える位相型格子であり、X線源11から射出されるX線がコーンビーム状である場合には、タルボ距離Zmは、下式13で表される。ここで、mは0または正の整数である。また、第1の格子21が吸収型格子であり、X線源11から射出されるX線が平行ビーム状である場合には、タルボ距離Zmは、下式14で表される。ここで、mは正の整数である。また、第1の格子21がX線にπ/2の位相変調を与える位相型格子であり、X線源11から射出されるX線が平行ビーム状である場合には、タルボ距離Zmは、下式15で表される。ここで、mは0または正の整数である。そして、第1の格子21がX線にπの位相変調を与える位相型格子であり、X線源11から射出されるX線が平行ビーム状である場合には、タルボ距離Zmは、下式16で表される。ここで、mは0または正の整数である。また、上記実施形態では、格子部12に第1及び第2の格子21,22の2つの格子を設けているが、第2の格子22を省略し、第1の格子21のみとすることも可能である。例えば、特開平2009-133823号公報に記されたX線画像検出器を用いることにより、第2の格子22を省略し、第1の格子21のみとすることが可能である。このX線画像検出器は、X線を電荷に変換する変換層と、変換層において変換された電荷を収集する電荷収集電極とを備えた直接変換型のX線画像検出器であり、各画素の電荷収集電極が複数の線状電極群を備える。1つの線状電極群は、一定の周期で配列された線状電極を互いに電気的に接続したものであり、他の線状電極群と互いに位相が異なるように配置されている。この線状電極群が第2の格子22として機能し、線状電極群が複数存在することにより、一度の撮影で位相の異なる複数のG2像の検出が行われる。したがって、この構成では、走査機構23を省略することが可能である。また、走査機構23を省略し、第1及び第2の格子21,22を介してX線画像検出器13により得られる単一の画像データに基づいて位相微分画像を生成する方法がある。この方法として、本出願人により特願2010-256241号として出願されている画素分割法がある。この画素分割法では、第1の格子21と第2の格子22とを、z方向の回りに僅かに回転させて、y方向に周期を有するモアレ縞をG2像に発生させる。X線画像検出器13により得られる単一の画像データを、該モアレ縞に対して互いに位相が異なる画素行x方向に並ぶ画素の群に分割し、分割された複数の画像データを、縞走査により互いに異なる複数のG2像に基づくものと見なして、上記縞走査法と同様な手順で位相微分画像を生成する。この画素分割法において、前述の強度変調信号は、単一の画像データに生じるモアレ縞の1周期分の画素値の強度変化として表される。さらに、画素分割法と同様に、走査機構23を省略し、第1及び第2の格子21,22を介してX線画像検出器13により得られる単一の画像データに基づいて位相微分画像を生成する方法として、WO2010/050483号公報に記載されたフーリエ変換法が知られている。このフーリエ変換法は、上記単一の画像データに対してフーリエ変換を行うことによりフーリエスペクトルを取得し、このフーリエスペクトルからキャリア周波数に対応したスペクトル位相情報を担うスペクトルを分離した後、逆フーリエ変換を行うことにより位相微分画像を生成する方法である。なお、このフーリエ変換法において、前述の強度変調信号は、画素分割法の場合と同様に、単一の画像データに生じるモアレ縞の1周期分の画素値の強度変化として表される。本発明は、医療診断用の放射線撮影装置の他に、工業用の放射線撮影装置等に適用することが可能である。また、放射線は、X線以外に、ガンマ線等を用いることも可能である。"