PROBLEM TO BE SOLVED: To simply select an optimum combination of CT images to be subjected to difference processing.SOLUTION: A storage part 13 extracts an X-ray high absorption area about X-ray high absorption substance from data on a non-contrast volume image of at least one time phase included in a non-contrast volume image group and data on a contrast volume image of a plurality of time phases included in a contrast volume image group. A parameter calculation part 17 calculates a feature amount that fluctuates in accordance with the movement of the X-ray high absorption area about each non-contrast volume image of at least one time phase and each contrast volume image of a plurality of time phases. An image selection part 19 selects a first non-contrast volume image and the contrast volume image in which image features of the X-ray high absorption areas are similar from among non-contrast volume images of at least one time phase and the contrast volume image of a plurality of time phases on the basis of the calculated feature amount.【課題】差分処理に供されるＣＴ画像の最適な組み合わせを簡便に選択すること。【解決手段】記憶部１３は、非造影ボリューム画像群に含まれる少なくとも一つの時相の非造影ボリューム画像のデータと造影ボリューム画像群に含まれる複数の時相の造影ボリューム画像のデータとから、Ｘ線高吸収物質に関するＸ線高吸収領域を抽出する。パラメータ算出部１７は、少なくとも一つの時相の非造影ボリューム画像の各々と複数の時相の造影ボリューム画像の各々とについて、Ｘ線高吸収領域の動きに応じて変動する特徴量を算出する。画像選択部１９は、算出された特徴量に基づいて、少なくとも一つの時相の非造影ボリューム画像と複数の時相の造影ボリューム画像との中から、Ｘ線高吸収領域の画像特徴が類似する第１の非造影ボリューム画像と造影ボリューム画像とを選択する。【選択図】 図１
