PROBLEM TO BE SOLVED: To check local abnormality of the lung as well as three-dimensional movement of the lung in aspiration.SOLUTION: A lung region extracting section 14 extracts a lung region from three-dimensional images Vi consisting of a plurality of time phases, and a positioning section 16 positions pixel positions corresponding to the lung region extracted from each of the three-dimensional images Vi, among the three-dimensional images Vi. Thereby, a displacement vector field in each time phase of the three-dimensional images Bi is calculated. A function calculating section 18 calculates a local ventilation amount function that represents temporal changes in ventilation amounts at each point of the displacement vector field, and a quantification section 20 calculates a difference function that is a function of a difference value of the local ventilation amount function and a reference ventilation amount function, as a quantitative value representing the difference between the local ventilation amount function and the reference ventilation amount function. A display control section 22 displays VR image of the lung in which the difference function is mapped, on a display 4 in chronological order.【課題】呼吸時における肺の３次元的な動きとともに、肺の局所的な異常を確認できるようにする。【解決手段】肺領域抽出部１４が、複数時相からなる３次元画像Ｖｉから肺領域を抽出し、位置合わせ部１６が、３次元画像Ｖｉのそれぞれから抽出した肺領域に対応する画素位置を、３次元画像Ｖｉ間において位置合わせする。これにより、３次元画像Ｖｉの各時相における変位ベクトル場が算出される。関数算出部１８が、変位ベクトル場の各点における換気量の時間的な変化を表す局所換気量関数を算出し、定量化部２０が、局所換気量関数および基準換気量関数の差分値の関数である差分関数を、局所換気量関数と基準換気量関数との相違を表す定量値として算出する。表示制御部２２が、差分関数がマッピングされた肺のＶＲ画像をディスプレイ４に時系列順に表示する。【選択図】図１
