A computer-implemented method for image mosaicing, wherein the method comprises: providing a mosaic ( M k-1) and an image ( I k), wherein the mosaic and the image each represent at least part of a scene, estimating a motion V ij k of blocks B ij k of the image ( I k), wherein the motion of a given block of a given image represents how the scene content of the given block has moved during a time interval before taking the given image, integrating (S31, S32) part of the image (I k) to the mosaic ( M k-1) according to the estimation, the method further comprising providing a plurality of images and iterating the providing, the estimating and the integrating over the plurality of images, wherein the integrating comprises: selecting an optimal set C J k of blocks B iJ k of the image ( I k), wherein the optimal set C J k is a set of blocks of the image ( I k) that minimizes a predetermined cost function, wherein the cost function takes as argument one set of blocks of an image and increasingly depends on a cumulated motion of the argument set, wherein the cumulated motion of a given set of blocks of a given image represents how the scene content of the given image at the level of the given set of blocks has moved during a time interval before taking the given image, and stitching the optimal set C j k to the mosaic ( M k), wherein the cumulated motion is a weighted sum of the motion V ij k of all the blocks B ij k of the image ( I k), wherein the weight of a given block decreasingly depends on a distance between the given block and the argument set. providing a mosaic ( M k-1) and an image ( I k), wherein the mosaic and the image each represent at least part of a scene, estimating a motion V ij k of blocks B ij k of the image ( I k), wherein the motion of a given block of a given image represents how the scene content of the given block has moved during a time interval before taking the given image, integrating (S31, S32) part of the image (I k) to the mosaic ( M k-1) according to the estimation, the method further comprising providing a plurality of images and iterating the providing, the estimating and the integrating over the plurality of images, wherein the integrating comprises: selecting an optimal set C J k of blocks B iJ k of the image ( I k), wherein the optimal set C J k is a set of blocks of the image ( I k) that minimizes a predetermined cost function, wherein the cost function takes as argument one set of blocks of an image and increasingly depends on a cumulated motion of the argument set, wherein the cumulated motion of a given set of blocks of a given image represents how the scene content of the given image at the level of the given set of blocks has moved during a time interval before taking the given image, and stitching the optimal set C j k to the mosaic ( M k), wherein the cumulated motion is a weighted sum of the motion V ij k of all the blocks B ij k of the image ( I k), wherein the weight of a given block decreasingly depends on a distance between the given block and the argument set. selecting an optimal set C J k of blocks B iJ k of the image ( I k), wherein the optimal set C J k is a set of blocks of the image ( I k) that minimizes a predetermined cost function, wherein the cost function takes as argument one set of blocks of an image and increasingly depends on a cumulated motion of the argument set, wherein the cumulated motion of a given set of blocks of a given image represents how the scene content of the given image at the level of the given set of blocks has moved during a time interval before taking the given image, and stitching the optimal set C j k to the mosaic ( M k), wherein the cumulated motion is a weighted sum of the motion V ij k of all the blocks B ij k of the image ( I k), wherein the weight of a given block decreasingly depends on a distance between the given block and the argument set.