A method of processing input digital image data representing an image into output digital image data, comprising: providing (100) input image data; providing and initializing (200) a level set function representation for approximating the output digital image data; providing (300) a speed function; determining (400) a propagation direction trend map relating to image points of the level set function, said map comprising one or more trend directions associated with said image points so that each one of said image points has an associated trend direction; updating (500, 520) the level set function using the speed function and the propagation direction trend map, so that the speed function is excluded from contributing in the updating of the level set function in image points for which the speed of the speed function is not in the associated trend direction; and providing (700) the output digital image data based on the level set function. By the exclusion, propagation in a coherent manner is enabled and speed penalty owing to local wiggling as observed in conventional methods can be avoided.