Method and computer program product for registering biomedical images with reduced imaging arefacts caused by object movement

A method of registering biomedical images to reduce imaging artifacts caused by object movement is disclosed, wherein a certain number of features or landmarks is automatically defined and tracked in a first and second image to be registered to determine the optical flow vector between the first and second image. Registration then takes place by applying the inverse optical flow to the second image. The automatic feature selection step is carried out by determining the mean signal intensity in a defined neighborhood for every pixel, which mean signal intensity is then compared to a predetermined threshold. If the mean signal intensity value of said neighborhood is higher than the predetermined threshold the pixel is defined as a feature and is added to a list of features to be tracked.