Abstract: Method for processing a sequence of at least two image datasets each having at least two pixels, each pixel having an intensity value by calculating a structure tensor for each pixel in an image dataset included in the sequence of image datasets; calculating values in a summation kernel based on the structure tensor for each pixel in the image dataset; calculating a weighted intensity value for each pixel in the first image dataset, using as weights the values in the summation kernel; storing the weighted intensity value for each pixel in the image dataset as a processed intensity value for each corresponding pixel in a processed output image dataset; rotating a local coordinate system in which the summation kernel is described so that the coordinate axes of said local coordinate system coincide with the directions of the eigenvectors of the structure tensor.