In the mask image (i) the first sheet, is set the first recording time according to the X-ray imaging condition, in the mask image (ii) the second and subsequent mask image capturing time calculating section 36, taken in the past average luminance of the mask image to be taken and the next mask image that was to be the target luminance, and calculating the second shooting time of the mask image to be photographed next based on the luminance and photographing time of the mask image taken in the past The live image capturing time calculating section 38 calculates the recording time of the live image on the basis of the actual recording time of the mask image X-ray is actually irradiated from the X-ray irradiator based on the second shooting time or the first The image processing unit 6 calculates a subtraction image from a difference between the live image and the reference mask image obtained by averaging the plurality of mask images taken on the basis of the second shooting time or the first. ([] Selection Figure Figure 4)マスク画像撮影時間算出部36は、(i)1枚目のマスク画像において、X線撮影条件に応じて第1撮影時間を設定し、(ii)2枚目以降のマスク画像において、過去に撮影したマスク画像と次に撮影するマスク画像との平均輝度が目標輝度となるように、過去に撮影したマスク画像の撮影時間と輝度とに基づいて次に撮影するマスク画像の第2撮影時間を算出し、ライブ画像撮影時間算出部38は、第1または第2撮影時間に基づいてX線照射器から実際にX線が照射されたマスク画像の実撮影時間を基にライブ画像の撮影時間を算出し、画像処理部6は、第1または第2撮影時間に基づいて撮影された複数のマスク画像を平均化した基準マスク画像とライブ画像との差分によりサブトラクション像を算出する。(【選択図】 図4)
