To perform a high-precision and high-speed blood flow dynamic analysis by applying a tissue blood flow change model even if the number of actually captured images is as small as a three-phase image. A pixel value of a region of interest set in a position of a blood vessel and an organ that becomes a flow path of contrast blood flow is measured based on medical images of a plurality of patients, a reference normalized signal intensity curve in a predetermined blood vessel, and the curve A database 41 that stores parameters related to the relationship between the blood vessel and a predetermined blood vessel, a database 42 that stores coefficients of a pharmacokinetic analysis model in association with a plurality of simulation points determined based on the curve, and a target patient that has undergone three-phase imaging Search unit 13 for searching for a simulation point corresponding to the combination of the contrast intensity of the early phase and the late phase of the region of interest at the position of the blood vessel and organ that becomes the flow path of the contrast blood flow based on the medical image of An acquisition unit 14 that acquires the attached coefficient and an estimation unit 15 that estimates a temporal contrast intensity curve of a predetermined blood vessel based on the coefficient are provided. [Selection] Figure 2【課題】実際に撮影した画像が３相画像のように少ない画像であっても、組織の血流変化モデルを適用して高精度かつ高速な血流動態解析を行うこと。【解決手段】複数患者の医用画像に基づき造影血流の流路となる血管及び臓器の位置に設定した関心領域の画素値を測定し、所定の血管における基準正規化信号強度曲線と、この曲線と所定の血管との関係に関するパラメータを保存するデータベース４１と、この曲線に基づき定めた複数のシミュレーション点に対応付けて薬物動態解析モデルの係数を保存したデータベース４２と、３相撮影された対象患者の医用画像に基づき、造影血流の流路となる血管及び臓器の位置における関心領域の早期相と後期相の造影強度の組み合わせに対応するシミュレーション点を検索する検索部１３と、シミュレーション点に対応付けられた係数を取得する取得部１４と、係数に基づき所定の血管の時間造影強度曲線を推定する推定部１５とを備える。【選択図】図２
