PROBLEM TO BE SOLVED: To increase a spatial resolution and to improve phase recovery accuracy in radiation phase imaging of acquiring phase information of a subject using Fourier transformation and inverse Fourier transformation.SOLUTION: A radiographic system 10 performs Fourier transformation on a radiographic image including a periodic pattern modulated by a subject arranged in a radiation field to acquire a spatial frequency spectrum of the radiographic image, separates a first frequency domain B on a straight line L passing through peak coordinates of a fundamental frequency component 62 of the periodic pattern and a second frequency domain C on either side of the straight line, in a frequency domain A including the fundamental frequency component from the spatial frequency spectrum, combines the first frequency domain, the second frequency domain and a third frequency domain D formed by rotating the second frequency domain by 180° with respect to the peak coordinates and reversing the positive/negative of the phase of the spatial frequency spectrum included there into one frequency domain E, performs inverse Fourier transformation on the combined frequency domain, and generates a phase contrast image.COPYRIGHT: (C)2013,JPO&INPIT【課題】フーリエ変換及び逆フーリエ変換を用いて被写体の位相情報を取得する放射線位相イメージングにおいて、空間分解能を高めると共に位相復元精度を高める。【解決手段】放射線撮影システム10は、放射線照射野に配置される被写体によって変調を受けた周期パターンを含む放射線画像に対してフーリエ変換を行って該放射線画像の空間周波数スペクトルを取得し、前記周期パターンの基本周波数成分62を含む周波数領域Aのうち、前記基本周波数成分のピーク座標を通る直線L上の第1の周波数領域B、及び該直線に関していずれか片側にある第2の周波数領域Cを、前記空間周波数スペクトルから分離し、前記第1の周波数領域及び前記第2の周波数領域、並びに前記第2の周波数領域を前記ピーク座標に関して180°回転すると共にそこに含まれる空間周波数スペクトルの位相の正負を逆転してなる第3の周波数領域D、を一つの周波数領域Eに合成して、合成された前記周波数領域に対して逆フーリエ変換を行って位相コントラスト画像を生成する。【選択図】図12