PROBLEM TO BE SOLVED: To highly accurately detect the frequency component of a cyclic pattern even if a radiographic image includes a region in which the cyclic pattern is inconspicuous.SOLUTION: An image acquiring section 31 acquires a radiographic image P0 from a radiation detector 2, and a region setting section 32 sets a plurality of small regions for analyzing frequency to the radiation image P0. An abundance setting section 33 detects an irradiation out-of-field region, a direct radiation region, a high absorbing body region and a high noise region, and sets the abundance of the cyclic pattern caused by a grid to each of these regions and other regions. A frequency analysis section 34 applies frequency analysis to each of the small regions A10 to acquire frequency spectrum, averages the acquired frequency spectrum by weighing in accordance with the abundance set in the region to which the small region belongs, and calculates the frequency spectrum of the radiation image P0. A filtering processing section 35 removes the frequency component of the cyclic pattern due to the grid from the radiation image P0.COPYRIGHT: (C)2013,JPO&INPIT【課題】放射線画像に周期的パターンが目立たない領域が含まれていても、周期的パターンの周波数成分を精度よく検出する。【解決手段】画像取得部31が放射線検出器2から放射線画像P0を取得し、領域設定部32が、周波数解析を行うための複数の小領域を放射線画像P0に設定する。存在度合設定部33が、照射野外領域、直接放射線領域、高吸収体領域および高ノイズ領域を検出し、これらの領域およびこれらの領域以外の領域のそれぞれに、グリッドに起因する周期的パターンの存在度合を設定する。周波数解析部34が、小領域A10のそれぞれについて周波数解析を行って周波数スペクトルを取得し、取得した周波数スペクトルを小領域が属する領域に設定した存在度合に応じて重みづけ平均して、放射線画像P0の周波数スペクトルを算出する。フィルタリング処理部35が、放射線画像P0からグリッドに起因する周期的パターンの周波数成分を除去する。【選択図】図2