A method for referencing at least one first optical image of the surface of the Earth taken by an optical sensor on board a satellite or on board an aircraft comprising the steps of: obtaining a stereoscopic optical image pair (19, 23) including the first optical image; obtaining at least one reference radar image taken by a synthetic aperture radar sensor, the surface area covered on the ground by the at least one reference radar image including an overlapping area with the surface area covered on the ground by the images of the stereoscopic optical image pair (19, 23); selecting at least one area of interest on the overlapping area; the method comprising for each of the areas of interest: obtaining a 3D model on the area of interest from the stereoscopic optical image pair (19, 23); calculating at least one simulated radar image on the area of interest from the obtained 3D model and the acquisition parameters of the at least one reference radar image; estimating a geometric offset (di, dj) between the at least one simulated radar image and the at least one reference radar image; selecting at least one reference point on the 3D model of the area of interest; projecting the at least one reference point in the at least one reference radar image by a radar projection function (P rad) so as to obtain at least one radar connection point; correcting the at least one radar connection point by applying the estimated geometric offset (di, dj) so as to obtain at least one corrected radar connection point determining at least one pair of connection points (46', 46") of the stereoscopic optical image pair (19, 23) by projection of said at least one reference point in each image of said stereoscopic optical image pair (19, 23); the method comprising a step of georeferencing of said at least one first optical image from at least the pair of optical connection points (46', 46") and from the at least one corrected radar connection point. obtaining a stereoscopic optical image pair (19, 23) including the first optical image; obtaining at least one reference radar image taken by a synthetic aperture radar sensor, the surface area covered on the ground by the at least one reference radar image including an overlapping area with the surface area covered on the ground by the images of the stereoscopic optical image pair (19, 23); selecting at least one area of interest on the overlapping area; the method comprising for each of the areas of interest: obtaining a 3D model on the area of interest from the stereoscopic optical image pair (19, 23); calculating at least one simulated radar image on the area of interest from the obtained 3D model and the acquisition parameters of the at least one reference radar image; estimating a geometric offset (di, dj) between the at least one simulated radar image and the at least one reference radar image; selecting at least one reference point on the 3D model of the area of interest; projecting the at least one reference point in the at least one reference radar image by a radar projection function (P rad) so as to obtain at least one radar connection point; correcting the at least one radar connection point by applying the estimated geometric offset (di, dj) so as to obtain at least one corrected radar connection point determining at least one pair of connection points (46', 46") of the stereoscopic optical image pair (19, 23) by projection of said at least one reference point in each image of said stereoscopic optical image pair (19, 23); obtaining a 3D model on the area of interest from the stereoscopic optical image pair (19, 23); calculating at least one simulated radar image on the area of interest from the obtained 3D model and the acquisition parameters of the at least one reference radar image; estimating a geometric offset (di, dj) between the at least one simulated radar image and the at least one reference radar image; selecting at least one reference point on the 3D model of the area of interest; projecting the at least one reference point in the at least one reference radar image by a radar projection function (P rad) so as to obtain at least one radar connection point; correcting the at least one radar connection point by applying the estimated geometric offset (di, dj) so as to obtain at least one corrected radar connection point determining at least one pair of connection points (46', 46") of the stereoscopic optical image pair (19, 23) by projection of said at least one reference point in each image of said stereoscopic optical image pair (19, 23); the method comprising a step of georeferencing of said at least one first optical image from at least the pair of optical connection points (46', 46") and from the at least one corrected radar connection point.