The systolic timings of myocardia at various positions of a heart are analyzed and evaluated. An image processing apparatus 1 obtains myocardial thicknesses at various positions of a heart within a plurality of three dimensional images V1 through VK, obtained by imaging a heart at a plurality of temporal phases within a single cardiac cycle, at each of the temporal phases. The image processing apparatus 1 obtains representative values that represent the systolic phase of the myocardia at each position, based on the obtained myocardial thicknesses at each temporal phase, and outputs the obtained representative values.