A method for motion compensation between first and second images in a temporal sequence includes processing the first and second images in a reduction process for providing respective reduced resolution first and second images; deriving respective first and second feature maps from the respective reduced resolution first and second images, the feature maps including deriving the respective Laplacian of image data in the respective reduced resolution first and second images; deriving a displacement field by processing the first and second feature maps in accordance with a registration algorithm, the registration algorithm comprising solving, for each picture element or voxel, a local Gaussian weighted least mean square problem so as to derive respective vectors forming the displacement field; and warping the second image with the displacement field.