X-ray images are recorded of a patient's heart and the heartbeat phase is registered as that is done. The heartbeat phases are coarsely divided into intervals and all X-ray images that have been assigned heartbeat phase from the interval are used for reconstructing a 3D image dataset. The movement fields of the other 3D image datasets are then calculated for one of said 3D image datasets. Movement fields are vector fields indicating the movements of similar structures from one local area to the other. A departure is then made from the coarse interval division, and for each heartbeat phase a movement field is interpolated individually or at least for fairly short intervals from the movement fields determined in advance, which field is used for generating a deformed 3D image dataset that has been imaged onto a reference heartbeat phase. The deformed 3D image datasets are then added together.