A method for geometric referencing of multi-spectral data in a sequence of images, the images containing distinct regions representing imaged radiation in respective distinct wavelengths. The method includes selecting a subset from the sequence of images, the subset containing a plurality of images, each image of which represents a field of view that overlaps with a field of view of at least one other image in the same subset; generating a 3D model by performing photogrammetric 3D reconstruction on the images in the subset; estimating a displacement and a yaw between pairs of successive images of the sequence of images; and using the 3D model and the estimated displacement and yaw to calculate photogrammetric parameters for images of the sequence of images that are not in the selected subset. Additionally, a corresponding system and a corresponding computer program product.