The digital image processing problem basing on velocity field determination is considered. An optimization method for constructing the velocity field basing on the study of the integral functional on the ensemble of trajectories is developed. The analytical form of functional variation allows using directed optimization methods to find the desired parameters. This method is considered for 3D tomographic image processing, in particular, for PET images for the movement correction in dynamic studies.